Showing papers in "Astronomy and Astrophysics in 2017"

The GMRT 150 MHz all-sky radio survey - First alternative data release TGSS ADR1

Abstract: We present the first full release of a survey of the 150 MHz radio sky, observed with the Giant Metrewave Radio Telescope (GMRT) between April 2010 and March 2012 as part of the TIFR GMRT Sky Survey (TGSS) project. Aimed at producing a reliable compact source survey, our automated data reduction pipeline efficiently processed more than 2000 h of observations with minimal human interaction. Through application of innovative techniques such as image-based flagging, direction-dependent calibration of ionospheric phase errors, correcting for systematic offsets in antenna pointing, and improving the primary beam model, we created good quality images for over 95 percent of the 5336 pointings. Our data release covers 36 900 deg2 (or 3.6 π steradians) of the sky between −53° and +90° declination (Dec), which is 90 percent of the total sky. The majority of pointing images have a noise level below 5 mJy beam-1 with an approximate resolution of 25′′×25′′ (or 25′′×25′′/ cos(Dec−19°) for pointings south of 19° declination). We have produced a catalog of 0.62 Million radio sources derived from an initial, high reliability source extraction at the seven sigma level. For the bulk of the survey, the measured overall astrometric accuracy is better than two arcseconds in right ascension and declination, while the flux density accuracy is estimated at approximately ten percent. Within the scope of the TGSS alternative data release (TGSS ADR) project, the source catalog, as well as 5336 mosaic images (5°×5°) and an image cutout service, are made publicly available at the CDS as a service to the astronomical community. Next to enabling a wide range of different scientific investigations, we anticipate that these survey products will provide a solid reference for various new low-frequency radio aperture array telescopes (LOFAR, LWA, MWA, SKA-low), and can play an important role in characterizing the epoch-of-reionisation (EoR) foreground. The TGSS ADR project aims at continuously improving the quality of the survey data products. Near-future improvements include replacement of bright source snapshot images with archival targeted observations, using new observations to fill the holes in sky coverage and replace very poor quality observational data, and an improved flux calibration strategy for less severely affected observational data.

The Sloan Digital Sky Survey Quasar Catalog: Twelfth data release

Abstract: PRIN INAF; "Investissements d'Avenir" French Government [ANR-11-IDEX-0001-02]; Agence Nationale de la Recherche [ANR-08-BLAN-0222, ANR-12-BS05-0015]; NASA ADAP [NNX12AE38G]; NSF [1211112, 1515404, AST-1516784]; Alfred P. Sloan Foundation; National Science Foundation; US Department of Energy Office of Science; University of Arizona; Brazilian Participation Group; Brookhaven National Laboratory; Carnegie Mellon University; University of Florida; French Participation Group; German Participation Group; Harvard University; Instituto de Astrofisica de Canarias; Michigan State/Notre Dame/JINA Participation Group; Johns Hopkins University; Lawrence Berkeley National Laboratory; Max Planck Institute for Astrophysics; Max Planck Institute for Extraterrestrial Physics; New Mexico State University; New York University; Ohio State University; Pennsylvania State University; University of Portsmouth; Princeton University; Spanish Participation Group; University of Tokyo; University of Utah; Vanderbilt University; University of Virginia; University of Washington; Yale University; National Aeronautics and Space Administration; Office of Science of the US Department of Energy [DE-AC02-05CH11231]

The LOFAR Two-metre Sky Survey. I. Survey description and preliminary data release

Abstract: The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire northern sky. Each of the 3170 pointings will be observed for 8 h, which, at most declinations, is sufficient to produce ~5? resolution images with a sensitivity of ~100 ?Jy/beam and accomplish the main scientific aims of the survey, which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Owing to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-Area low-frequency survey. In excess of 44 000 sources are detected in the images that have a resolution of 25?, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45°00?00? to 57°00?00?).

AGN wind scaling relations and the co-evolution of black holes and galaxies

Abstract: This work was supported by ASI/INAF contract I/009/10/0 and INAF PRIN 2011, 2012 and 2014. MB acknowledges support from the FP7 Career Integration Grant “eEASy” (CIG 321913). LZ acknowledges support from ASI/INAF grant I/037/12/0. CF acknowledges funding from the European Union Horizon 2020 research and innovation programme under the Marie SklodowskaCurie grant agreement No 664931. CC acknowledges funding from the European Union Horizon 2020 research and innovation programme under the Marie SklodowskaCurie grant agreement No 664931 and support from Swiss National Science Foundation Grants PP00P2 138979 and PP00P2 166159. RM acknowledges the ERC Advanced Grant 695671 QUENCH and support from the Science and Technology Facilities Council (STFC).

First light for GRAVITY: Phase referencing optical interferometry for the Very Large Telescope Interferometer

Abstract: GRAVITY is a new instrument to coherently combine the light of the European Southern Observatory Very Large Telescope Interferometer to form a telescope with an equivalent 130 m diameter angular resolution and a collecting area of 200 m2. The instrument comprises fiber fed integrated optics beam combination, high resolution spectroscopy, built-in beam analysis and control, near-infrared wavefront sensing, phase-tracking, dual-beam operation, and laser metrology. GRAVITY opens up to optical/infrared interferometry the techniques of phase referenced imaging and narrow angle astrometry, in many aspects following the concepts of radio interferometry. This article gives an overview of GRAVITY and reports on the performance and the first astronomical observations during commissioning in 2015/16. We demonstrate phase-tracking on stars as faint as mK ≈ 10 mag, phase-referenced interferometry of objects fainter than mK ≈ 15 mag with a limiting magnitude of mK ≈ 17 mag, minute long coherent integrations, a visibility accuracy of better than 0.25%, and spectro-differential phase and closure phase accuracy better than 0.5°, corresponding to a differential astrometric precision of better than ten microarcseconds (μas). The dual-beam astrometry, measuring the phase difference of two objects with laser metrology, is still under commissioning. First observations show residuals as low as 50 μas when following objects over several months. We illustrate the instrument performance with the observations of archetypical objects for the different instrument modes. Examples include the Galactic center supermassive black hole and its fast orbiting star S2 for phase referenced dual-beam observations and infrared wavefront sensing, the high mass X-ray binary BP Cru and the active galactic nucleus of PDS 456 for a few μas spectro-differential astrometry, the T Tauri star S CrA for a spectro-differential visibility analysis, ξ Tel and 24 Cap for high accuracy visibility observations, and η Car for interferometric imaging with GRAVITY.

Photodissociation and photoionisation of atoms and molecules of astrophysical interest

Abstract: A new collection of photodissociation and photoionisation cross sections for 102 atoms and molecules of astrochemical interest has been assembled, along with a brief review of the basic physical processes involved. These have been used to calculate dissociation and ionisation rates, with uncertainties, in a standard ultraviolet interstellar radiation field (ISRF) and for other wavelength-dependent radiation fields, including cool stellar and solar radiation, Lyman-α dominated radiation, and a cosmic-ray induced ultraviolet flux. The new ISRF rates generally agree within 30% with our previous compilations, with a few notable exceptions. Comparison with other databases such as PHIDRATES is made. The reduction of rates in shielded regions was calculated as a function of dust, molecular and atomic hydrogen, atomic C, and self-shielding column densities. The relative importance of these shielding types depends on the atom or molecule in question and the assumed dust optical properties. All of the new data are publicly available from the Leiden photodissociation and ionisation database. Sensitivity of the calculated rates to variation of temperature and isotope, and uncertainties in measured or calculated cross sections, are tested and discussed. Tests were conducted on the new rates with an interstellar-cloud chemical model, and find general agreement (within a factor of two) in abundances obtained with the previous iteration of the Leiden database assuming an ISRF, and order-of-magnitude variations assuming various kinds of stellar radiation. The newly parameterised dust-shielding factors makes a factor-of-two difference to many atomic and molecular abundances relative to parameters currently in the UDfA and KIDA astrochemical reaction databases. The newly-calculated cosmic-ray induced photodissociation and ionisation rates differ from current standard values up to a factor of 5. Under high temperature and cosmic-ray-flux conditions the new rates alter the equilibrium abundances of abundant dark cloud abundances by up to a factor of two. The partial cross sections for H2 O and NH3 photodissociation forming OH, O, NH2 and NH are also evaluated and lead to radiation-field-dependent branching ratios.

The COSMOS2015 galaxy stellar mass function . Thirteen billion years of stellar mass assembly in ten snapshots

Abstract: We measure the stellar mass function (SMF) and stellar mass density of galaxies in the COSMOS field up to z ~ 6. We select them in the near-IR bands of the COSMOS2015 catalogue, which includes ultra-deep photometry from UltraVISTA-DR2, SPLASH, and Subaru/Hyper Suprime-Cam. At z > 2.5 we use new precise photometric redshifts with error σ_z = 0.03(1 + z) and an outlier fraction of 12%, estimated by means of the unique spectroscopic sample of COSMOS (~100 000 spectroscopic measurements in total, more than one thousand having robust z_(spec)> 2.5). The increased exposure time in the DR2, along with our panchromatic detection strategy, allow us to improve the completeness at high z with respect to previous UltraVISTA catalogues (e.g. our sample is >75% complete at 10^(10)M⊙ and z = 5). We also identify passive galaxies through a robust colour–colour selection, extending their SMF estimate up to z = 4. Our work provides a comprehensive view of galaxy-stellar-mass assembly between z = 0.1 and 6, for the first time using consistent estimates across the entire redshift range. We fit these measurements with a Schechter function, correcting for Eddington bias. We compare the SMF fit with the halo mass function predicted from ΛCDM simulations, finding that at z > 3 both functions decline with a similar slope in thehigh-mass end. This feature could be explained assuming that mechanisms quenching star formation in massive haloes become less effective at high redshifts; however further work needs to be done to confirm this scenario. Concerning the SMF low-mass end, it shows a progressive steepening as it moves towards higher redshifts, with α decreasing from -1.47^(+0.02)_(-0.02) at z ≃ 0.1 to -2.11^(+0.30)_(-0.13) at z ≃ 5. This slope depends on the characterisation of the observational uncertainties, which is crucial to properly remove the Eddington bias. We show that there is currently no consensus on the method to quantify such errors: different error models result in different best-fit Schechter parameters.

Measurement of baryon acoustic oscillation correlations at z = 2.3 with SDSS DR12 Ly α-Forests

Abstract: Author(s): Bautista, JE; Busca, NG; Guy, J; Rich, J; Blomqvist, M; Du Mas Des Bourboux, H; Pieri, MM; Font-Ribera, A; Bailey, S; Delubac, T; Kirkby, D; Le Goff, JM; Margala, D; Slosar, A; Vazquez, JA; Brownstein, JR; Dawson, KS; Eisenstein, DJ; Miralda-Escude, J; Noterdaeme, P; Palanque-Delabrouille, N; Pâris, I; Petitjean, P; Ross, NP; Schneider, DP; Weinberg, DH; Yeche, C | Abstract: We have used flux-transmission correlations in Lyα forests to measure the imprint of baryon acoustic oscillations (BAO). The study uses spectra of 157 783 quasars in the redshift range 2:1 ≤ z ≤ 3:5 from the Sloan Digital Sky Survey (SDSS) data release 12 (DR12). Besides the statistical improvements on our previous studies using SDSS DR9 and DR11, we have implemented numerous improvements in the analysis procedure, allowing us to construct a physical model of the correlation function and to investigate potential systematic errors in the determination of the BAO peak position. The Hubble distance, DH = c=H(z), relative to the sound horizon is DH(z = 2:33)=rd = 9:07 ± 0:31. The best-determined combination of comoving angular-diameter distance, DM, and the Hubble distance is found to be D0:7 H D0:3 M =rd = 13:94±0:35. This value is 1:028±0:026 times the prediction of the flat-ΛCDM model consistent with the cosmic microwave background (CMB) anisotropy spectrum. The errors include marginalization over the effects of unidentified high-density absorption systems and fluctuations in ultraviolet ionizing radiation. Independently of the CMB measurements, the combination of our results and other BAO observations determine the open-ΛCDM density parameters to be ΩM = 0:296 ± 0:029, ΩΛ = 0:699 ± 0:100 and Ωk = -0:002 ± 0:119.

Spectroscopy Made Easy: Evolution

Abstract: Context. The Spectroscopy Made Easy (SME) package has become a popular tool for analyzing stellar spectra, often in connection with large surveys or exoplanet research. SME has evolved significantl ...

The VLA-COSMOS 3 GHz Large Project: Continuum data and source catalog release

Abstract: We present the VLA-COSMOS 3 GHz Large Project based on 384 h of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 μ Jy beam-1 over the two square degrees at an angular resolution of 0.75″. To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10 830 radio sources down to 5σ , out of which 67 are combined from multiple components.Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01″ at the bright end (signal-to-noise ratio, S /N 3 GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4 GHz Joint Project catalog yields an expected median radio spectral index of α = −0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75″) angular resolution to date, bridging the gap between last-generation and next-generation surveys.

X-shooter spectroscopy of young stellar objects in Lupus. Accretion properties of class II and transitional objects

Abstract: The mass accretion rate, Ṁ acc , is a key quantity for the understanding of the physical processes governing the evolution of accretion discs around young low-mass ( M ⋆ ≲ 2.0 M ⊙ ) stars and substellar objects (YSOs). We present here the results of a study of the stellar and accretion properties of the (almost) complete sample of class II and transitional YSOs in the Lupus I, II, III and IV clouds, based on spectroscopic data acquired with the VLT/X-shooter spectrograph. Our study combines the dataset from our previous work with new observations of 55 additional objects. We have investigated 92 YSO candidates in total, 11 of which have been definitely identified with giant stars unrelated to Lupus. The stellar and accretion properties of the 81 bona fide YSOs, which represent more than 90% of the whole class II and transition disc YSO population in the aforementioned Lupus clouds, have been homogeneously and self-consistently derived, allowing for an unbiased study of accretion and its relationship with stellar parameters. The accretion luminosity, L acc , increases with the stellar luminosity, L ⋆ , with an overall slope of ~1.6, similar but with a smaller scatter than in previous studies. There is a significant lack of strong accretors below L ⋆ ≈ 0.1 L ⊙ , where L acc is always lower than 0.01 L ⋆ . We argue that the L acc − L ⋆ slope is not due to observational biases, but is a true property of the Lupus YSOs. The log Ṁ acc – log M ⋆ correlation shows a statistically significant evidence of a break, with a steeper relation for M ⋆ ≲ 0.2 M ⊙ and a flatter slope for higher masses. The bimodality of the Ṁ acc – M ⋆ relation is confirmed with four different evolutionary models used to derive the stellar mass. The bimodal behaviour of the observed relationship supports the importance of modelling self-gravity in the early evolution of the more massive discs, but other processes, such as photo-evaporation and planet formation during the YSO’s lifetime, may also lead to disc dispersal on different timescales depending on the stellar mass. The sample studied here more than doubles the number of YSOs with homogeneously and simultaneously determined L acc and luminosity, L line , of many permitted emission lines. Hence, we also refined the empirical relationships between L acc and L line on a more solid statistical basis.

Baryon acoustic oscillations from the complete SDSS-III Ly α -quasar cross-correlation function at z = 2.4

Abstract: We present a measurement of baryon acoustic oscillations (BAO) in the cross-correlation of quasars with the Lyα-forest flux transmission at a mean redshift of z = 2.40. The measurement uses the complete Sloan Digital Sky Survey (SDSS-III) data sample: 168 889 forests and 234 367 quasars from the SDSS data release DR12. In addition to the statistical improvement on our previous study using DR11, we have implemented numerous improvements at the analysis level enabling a more accurate measurement of this cross-correlation. We have also developed the first simulations of the cross-correlation that allow us to test different aspects of our data analysis and to search for potential systematic errors in the determination of the BAO peak position. We measure the two ratios DH(z = 2.40) /rd = 9.01 ± 0.36 and DM(z = 2.40) /rd = 35.7 ± 1.7, where the errors include marginalization over the non-linear velocity of quasars and the cross-correlation of metals and quasars, among other effects. These results are within 1.8σ of the prediction of the flat-ΛCDM model describing the observed cosmic microwave background anisotropies. We combine this study with the Lyα-forest auto-correlation function, yielding DH(z = 2.40) /rd = 8.94 ± 0.22 and DM(z = 2.40) /rd = 36.6 ± 1.2, within 2.3σ of the same flat-ΛCDM model. Key words: cosmological parameters / dark energy / large-scale structure of Universe

The MUSE Hubble Ultra Deep Field Survey. I. Survey description, data reduction, and source detection

Abstract: We present the MUSE Hubble Ultra Deep Survey, a mosaic of nine MUSE fields covering 90% of the entire HUDF region with a 10-h deep exposure time, plus a deeper 31-h exposure in a single 1.15 arcmin 2 field. The improved observing strategy and advanced data reduction results in datacubes with sub-arcsecond spatial resolution (0.65 at 7000 A) and accurate astrometry (0.07 rms). We compare the broadband photometric properties of the datacubes to HST photometry, finding a good agreement in zeropoint up to m AB = 28 but with an increasing scatter for faint objects. We have investigated the noise properties and developed an empirical way to account for the impact of the correlation introduced by the 3D drizzle interpolation. The achieved 3 σ emission line detection limit for a point source is 1.5 and 3.1 × 10 -19 erg s -1 cm -2 for the single ultra-deep datacube and the mosaic, respectively. We extracted 6288 sources using an optimal extraction scheme that takes the published HST source locations as prior. In parallel, we performed a blind search of emission line galaxies using an original method based on advanced test statistics and filter matching. The blind search results in 1251 emission line galaxy candidates in the mosaic and 306 in the ultradeep datacube, including 72 sources without HST counterparts ( m AB > 31). In addition 88 sources missed in the HST catalog but with clear HST counterparts were identified. This data set is the deepest spectroscopic survey ever performed. In just over 100 h of integration time, it provides nearly an order of magnitude more spectroscopic redshifts compared to the data that has been accumulated on the UDF over the past decade. The depth and high quality of these datacubes enables new and detailed studies of the physical properties of the galaxy population and their environments over a large redshift range.

ALMA unveils rings and gaps in the protoplanetary system HD 169142: signatures of two giant protoplanets

Abstract: The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has recently been detected by direct imaging in the near-infrared. To study the interaction between the protoplanet and the disk itself, observations of the gas and dust surface density structure are needed. This paper reports new ALMA observations of the dust continuum at 1.3 mm, 12 CO, 13 CO, and C18 O J = 2−1 emission from the system HD 169142 (which is observed almost face-on) at an angular resolution of ~ (~35 × 20 au). The dust continuum emission reveals a double-ring structure with an inner ring between (~20−35 au) and an outer ring between (~56−83 au). The size and position of the inner ring is in good agreement with previous polarimetric observations in the near-infrared and is consistent with dust trapping by a massive planet. No dust emission is detected inside the inner dust cavity (R ≲ 20 au) or within the dust gap (~35−56 au) down to the noise level. In contrast, the channel maps of the J = 2−1 line of the three CO isotopologs reveal gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission; it extends to ~1 (~180 au) in radius. This difference and the sharp drop of the continuum emission at large radii point to radial drift of large dust grains (>μ m size). Using the thermo-chemical disk code dali, we modeled the continuum and the CO isotopolog emission to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of ~30−40 inward of the dust gap. The gas and dust distribution indicate that two giant planets shape the disk structure through dynamical clearing (dust cavity and gap) and dust trapping (double-ring dust distribution).

Planetesimal formation starts at the snow line

Abstract: Context. The formation stage of planetesimals represents a major gap in our understanding of the planet formation process. Late-stage planet accretion models typically make arbitrary assumptions about planetesimal and pebble distribution, while dust evolution models predict that planetesimal formation is only possible at some orbital distances. Aims. We wish to test the importance of the water snow line in triggering the formation of the first planetesimals during the gas-rich phase of a protoplanetary disk, when cores of giant planets have to form. Methods. We connected prescriptions for gas disk evolution, dust growth and fragmentation, water ice evaporation and recondensation, the transport of both solids and water vapor, and planetesimal formation via streaming instability into a single one-dimensional model for protoplanetary disk evolution. Results. We find that processes taking place around the snow line facilitate planetesimal formation in two ways. First, because the sticking properties between wet and dry aggregates change, a “traffic jam” inside of the snow line slows the fall of solids onto the star. Second, ice evaporation and outward diffusion of water followed by its recondensation increases the abundance of icy pebbles that trigger planetesimal formation via streaming instability just outside of the snow line. Conclusions. Planetesimal formation is hindered by growth barriers and radial drift and thus requires particular conditions to take place. The snow line is a favorable location where planetesimal formation is possible for a wide range of conditions, but not in every protoplanetary disk model, however. This process is particularly promoted in large cool disks with low intrinsic turbulence and an increased initial dust-to-gas ratio.

Discovery of a warm, dusty giant planet around HIP 65426

Abstract: Aims. The SHINE program is a high-contrast near-infrared survey of 600 young, nearby stars aimed at searching for and characterizing new planetary systems using VLT/SPHERE’s unprecedented high-contrast and high-angular-resolution imaging capabilities. It is also intended to place statistical constraints on the rate, mass and orbital distributions of the giant planet population at large orbits as a function of the stellar host mass and age to test planet-formation theories. Methods. We used the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-contrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP 65426. It is a member of the ~17 Myr old Lower Centaurus-Crux association. Results. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP 65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2 μm indicate a warm, dusty atmosphere characteristic of young low-surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6–12 MJup, Teff = 1300–1600 K and R = 1.5 ± 0.1 RJup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log (g) = 4.0–5.0 with smaller radii (1.0–1.3 RJup). Conclusions. Given its physical and spectral properties, HIP 65426 b occupies a rather unique placement in terms of age, mass, and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution.

Concentrating small particles in protoplanetary disks through the streaming instability

Abstract: Laboratory experiments indicate that direct growth of silicate grains via mutual collisions can only produce particles up to roughly millimeters in size. On the other hand, recent simulations of the streaming instability have shown that mm/cm-sized particles require an excessively high metallicity for dense filaments to emerge. Using a numerical algorithm for stiff mutual drag force, we perform simulations of small particles with significantly higher resolutions and longer simulation times than in previous investigations. We find that particles of dimensionless stopping time τs = 10-2 and 10-3 - representing cm- and mm-sized particles interior of the water ice line - concentrate themselves via the streaming instability at a solid abundance of a few percent. We thus revise a previously published critical solid abundance curve for the regime of τs ≪ 1. The solid density in the concentrated regions reaches values higher than the Roche density, indicating that direct collapse of particles down to mm sizes into planetesimals is possible. Our results hence bridge the gap in particle size between direct dust growth limited by bouncing and the streaming instability. (Less)

Three-dimensional mapping of the local interstellar medium with composite data

Abstract: Context. Three-dimensional maps of the Galactic interstellar medium are general astrophysical tools. Reddening maps may be based on the inversion of color excess measurements for individual target stars or on statistical methods using stellar surveys. Three-dimensional maps based on diffuse interstellar bands (DIBs) have also been produced. All methods benefit from the advent of massive surveys and may benefit from Gaia data.Aims. All of the various methods and databases have their own advantages and limitations. Here we present a first attempt to combine different datasets and methods to improve the local maps.Methods. We first updated our previous local dust maps based on a regularized Bayesian inversion of individual color excess data by replacing Hipparcos or photometric distances with Gaia Data Release 1 values when available. Secondly, we complemented this database with a series of ≃5000 color excess values estimated from the strength of the λ 15273 DIB toward stars possessing a Gaia parallax. The DIB strengths were extracted from SDSS/APOGEE spectra. Third, we computed a low-resolution map based on a grid of Pan-STARRS reddening measurements by means of a new hierarchical technique and used this map as the prior distribution during the inversion of the two other datasets.Results. The use of Gaia parallaxes introduces significant changes in some areas and globally increases the compactness of the structures. Additional DIB-based data make it possible to assign distances to clouds located behind closer opaque structures and do not introduce contradictory information for the close structures. A more realistic prior distribution instead of a plane-parallel homogeneous distribution helps better define the structures. We validated the results through comparisons with other maps and with soft X-ray data.Conclusions. Our study demonstrates that the combination of various tracers is a potential tool for more accurate maps. An online tool makes it possible to retrieve maps and reddening estimations.

Lyman-α spectral properties of five newly discovered Lyman continuum emitters

Abstract: Aims. We have recently reported the discovery of five low redshift Lyman continuum (LyC) emitters (LCEs) with absolute escape fractions f esc LyC that range from 6 to 13%, higher than previously found and that more than double the number of low redshift LCEs. We use these observations to test theoretical predictions about a link between the characteristics of the Lyman-alpha (Lyα ) line from galaxies and the escape of ionizing photons. Methods. We analyse the Lyα spectra of eight LCEs of the local Universe observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope (our five leakers and three galaxies from the litterature), and compare their strengths and shapes to the theoretical criteria and comparison samples of local galaxies: the Lyman Alpha Reference Survey, Lyman-break analogs, Green Peas, and the high-redshift strong LyC leaker Ion2 .Results. Our LCEs are found to be strong Lyα emitters, with high equivalent widths, EW (Lyα ) > 70 A, and large Lyα escape fractions, f esc Lyα > 20%. The Lyα profiles are all double-peaked with a small peak separation, in agreement with our theoretical expectations. They also have no underlying absorption at the Lyα position. All these characteristics are very different from the Lyα properties of typical star-forming galaxies of the local Universe. A subset of the comparison samples (2−3 Green Pea galaxies) share these extreme values, indicating that they could also be leaking. We also find a strong correlation between the star formation rate surface density and the escape fraction of ionizing photons, indicating that the compactness of star-forming regions plays a role in shaping low column density paths in the interstellar medium of LCEs. Conclusions. The Lyα properties of LCEs are peculiar: Lyα can be used as a reliable tracer of LyC escape from galaxies, complementing other indirect diagnostics proposed in the literature.

Planetesimal formation near the snowline: in or out?

Abstract: Context. The formation of planetesimals in protoplanetary disks is not well-understood. Streaming instability is a promising mechanism to directly form planetesimals from pebble-sized particles, provided a high enough solids-to-gas ratio. However, local enhancements of the solids-to-gas ratio are difficult to realize in a smooth disk, which motivates the consideration of special disk locations such as the snowline – the radial distance from the star beyond which water can condense into solid ice.Aims. In this article we investigate the viability of planetesimal formation by streaming instability near the snowline due to water diffusion and condensation. We aim to identify under what disk conditions streaming instability can be triggered near the snowline.Methods. To this end, we adopt a viscous disk model, and numerically solve the transport equations for vapor and solids on a cylindrical, 1D grid. We take into account radial drift of solids, gas accretion on to the central star, and turbulent diffusion. We study the importance of the back-reaction of solids on the gas and of the radial variation of the mean molecular weight of the gas. Different designs for the structure of pebbles are investigated, varying in the number and size of silicate grains. We also introduce a semi-analytical model that we employ to obtain results for different disk model parameters.Results. We find that water diffusion and condensation can locally enhance the ice surface density by a factor 3–5 outside the snowline. Assuming that icy pebbles contain many micron-sized silicate grains that are released during evaporation, the enhancement is increased by another factor ~2. In this “many-seeds” model, the solids-to-gas ratio interior to the snowline is enhanced as well, but not as much as just outside the snowline. In the context of a viscous disk, the diffusion-condensation mechanism is most effective for high values of the turbulence parameter α (10-3 –10-2 ). Therefore, assuming young disks are more vigorously turbulent than older disks, planetesimals near the snowline can form in an early stage of the disk. In highly turbulent disks, tens of Earth masses can be stored in an annulus outside the snowline, which can be identified with recent ALMA observations.

Global simulations of protoplanetary disks with net magnetic flux - I. Non-ideal MHD case

Abstract: Context. The planet-forming region of protoplanetary disks is cold, dense, and therefore weakly ionized. For this reason, magnetohydrodynamic (MHD) turbulence is thought to be mostly absent, and another mechanism has to be found to explain gas accretion. It has been proposed that magnetized winds, launched from the ionized disk surface, could drive accretion in the presence of a large-scale magnetic field. Aims. The efficiency and the impact of these surface winds on the disk structure is still highly uncertain. We present the first global simulations of a weakly ionized disk that exhibits large-scale magnetized winds. We also study the impact of self-organization, which was previously demonstrated only in non-stratified models. Methods. We perform numerical simulations of stratified disks with the PLUTO code. We compute the ionization fraction dynamically, and account for all three non-ideal MHD effects: ohmic and ambipolar diffusions, and the Hall drift. Simplified heating and cooling due to non-thermal radiation is also taken into account in the disk atmosphere.Results. We find that disks can be accreting or not, depending on the configuration of the large-scale magnetic field. Magnetothermal winds, driven both by magnetic acceleration and heating of the atmosphere, are obtained in the accreting case. In some cases, these winds are asymmetric, ejecting predominantly on one side of the disk. The wind mass loss rate depends primarily on the average ratio of magnetic to thermal pressure in the disk midplane. The non-accreting case is characterized by a meridional circulation, with accretion layers at the disk surface and decretion in the midplane. Finally, we observe self-organization, resulting in axisymmetric rings of density and associated pressure “bumps”. The underlying mechanism and its impact on observable structures are discussed.

The global dust modelling framework THEMIS

Abstract: Here we introduce the interstellar dust modelling framework THEMIS (The Heterogeneous dust Evolution Model for Interstellar Solids), which takes a global view of dust and its evolution in response to the local conditions in interstellar media. This approach is built upon a core model that was developed to explain the dust extinction and emission in the diffuse interstellar medium. The model was then further developed to self-consistently include the effects of dust evolution in the transition to denser regions. The THEMIS approach is under continuous development and we are currently extending the framework to explore the implications of dust evolution in HII regions and the photon-dominated regions associated with star formation. We provide links to the THEMIS, DustEM and DustPedia websites where more information about the model, its input data and applications can be found.

Limb and gravity-darkening coefficients for the TESS satellite at several metallicities, surface gravities, and microturbulent velocities

Abstract: The Spanish MEC (AYA2015-71718-R and ESP2015-65712-C5-5-R) is gratefully acknowledged for its support during the development of this work.

The third data release of the Kilo-Degree Survey and associated data products

Abstract: Context. The Kilo-Degree Survey (KiDS) is an ongoing optical wide-field imaging survey with the OmegaCAM camera at the VLT Survey Telescope. It aims to image 1500 square degrees in four filters (ugri ). The core science driver is mapping the large-scale matter distribution in the Universe, using weak lensing shear and photometric redshift measurements. Further science cases include galaxy evolution, Milky Way structure, detection of high-redshift clusters, and finding rare sources such as strong lenses and quasars.Aims. Here we present the third public data release and several associated data products, adding further area, homogenized photometric calibration, photometric redshifts and weak lensing shear measurements to the first two releases.Methods. A dedicated pipeline embedded in the Astro-WISE information system is used for the production of the main release. Modifications with respect to earlier releases are described in detail. Photometric redshifts have been derived using both Bayesian template fitting, and machine-learning techniques. For the weak lensing measurements, optimized procedures based on the THELI data reduction and lens fit shear measurement packages are used. Results. In this third data release an additional 292 new survey tiles (≈300 deg2 ) stacked ugri images are made available, accompanied by weight maps, masks, and source lists. The multi-band catalogue, including homogenized photometry and photometric redshifts, covers the combined DR1, DR2 and DR3 footprint of 440 survey tiles (44 deg2 ). Limiting magnitudes are typically 24.3, 25.1, 24.9, 23.8 (5σ in a 2′′ aperture) in ugri , respectively, and the typical r -band PSF size is less than 0.7′′. The photometric homogenization scheme ensures accurate colours and an absolute calibration stable to ≈2% for gri and ≈3% in u . Separately released for the combined area of all KiDS releases to date are a weak lensing shear catalogue and photometric redshifts based on two different machine-learning techniques.

Lupus disks with faint CO isotopologues: low gas/dust or high carbon depletion?

Abstract: Context. An era has started in which gas and dust can be observed independently in protoplanetary disks, thanks to the recent surveys with the Atacama Large Millimeter/sub-millimeter Array (ALMA). The first near-complete high-resolution disk survey in both dust and gas in a single star-forming region has been carried out in Lupus, finding surprisingly low gas-to-dust ratios.Aims. The goal of this work is to fully exploit CO isotopologue observations in Lupus, comparing them with physical-chemical model results, in order to obtain gas masses for a large number of disks and compare gas and dust properties.Methods. We have employed the grid of physical-chemical models presented previously to analyze continuum and CO isotopologue (13 CO J = 3−2 and C18 O J = 3−2) observations of Lupus disks, including isotope-selective processes and freeze-out. We also employed the ALMA 13 CO-only detections to calculate disk gas masses for a total of 34 sources, which expands the sample of 10 disks reported earlier, where C18 O was also detected.Results. We confirm that overall gas-masses are very low, often lower than 1M J , when volatile carbon is not depleted. Accordingly, global gas-to-dust ratios are much lower than the expected interstellar-medium value of 100, which is predominantly between 1 and 10. Low CO-based gas masses and gas-to-dust ratios may indicate rapid loss of gas, or alternatively chemical evolution, for example, through sequestering of carbon from CO to more complex molecules, or carbon locked up in larger bodies.Conclusions. Current ALMA observations of 13 CO and continuum emission cannot distinguish between these two hypotheses. We have simulated both scenarios, but chemical model results do not allow us to rule out one of the two, pointing to the need to calibrate CO-based masses with other tracers. Assuming that all Lupus disks have evolved mainly as a result of viscous processes over the past few Myr, the previously observed correlation between the current mass accretion rate and dust mass implies a constant gas-to-dust ratio, which is close to 100 based on the observed M disk /Ṁ acc ratio. This in turn points to a scenario in which carbon depletion is responsible for the low luminosities of the CO isotopologue line.

Shadows and spirals in the protoplanetary disk HD 100453

Abstract: ESO; CNRS (France); MPIA (Germany); INAF (Italy); FINES (Switzerland); NOVA (Netherlands); European Commission, Optical Infrared Coordination Network for Astronomy (OPTICON) [RII3-Ct-2004-001566, 226604, 312430]; Programme National de Planetologie (PNP); Programme National de Physique Stellaire (PNPS) of CNRS-INSU; French Labex OSUG@2020 (Investissements d'avenir) [ANR10 LABX56]; NASA's Science Mission Directorate; DISCSIM project - European Research Council [341137]; Italian Ministry of Education, University, and Research

The GAPS Programme with HARPS-N at TNG - XIV. Investigating giant planet migration history via improved eccentricity and mass determination for 231 transiting planets

Abstract: We carried out a Bayesian homogeneous determination of the orbital parameters of 231 transiting giant planets (TGPs) that are alone or have distant companions; we employed differential evolution Markov chain Monte Carlo methods to analyse radial-velocity (RV) data from the literature and 782 new high-accuracy RVs obtained with the HARPS-N spectrograph for 45 systems over ~3 years. Our work yields the largest sample of systems with a transiting giant exoplanet and coherently determined orbital, planetary, and stellar parameters. We found that the orbital parameters of TGPs in non-compact planetary systems are clearly shaped by tides raised by their host stars. Indeed, the most eccentric planets have relatively large orbital separations and/or high mass ratios, as expected from the equilibrium tide theory. This feature would be the outcome of planetary migration from highly eccentric orbits excited by planet-planet scattering, Kozai-Lidov perturbations, or secular chaos. The distribution of α = a / a R , where a and a R are the semi-major axis and the Roche limit, for well-determined circular orbits peaks at 2.5; this agrees with expectations from the high-eccentricity migration (HEM), although it might not be limited to this migration scenario. The few planets of our sample with circular orbits and α > 5 values may have migrated through disc-planet interactions instead of HEM. By comparing circularisation times with stellar ages, we found that hot Jupiters with a 5 ≲ Q' p ≲ 10 9 , and that stellar Q' s ≳ 10 6 - 10 7 are required to explain the presence of eccentric planets at the same orbital distance. As aby-product of our analysis, we detected a non-zero eccentricity e = 0.104 -0.018 +0.021 for HAT-P-29; we determined that five planets that were previously regarded to be eccentric or to have hints of non-zero eccentricity, namely CoRoT-2b, CoRoT-23b, TrES-3b, HAT-P-23b, and WASP-54b, have circular orbits or undetermined eccentricities; we unveiled curvatures caused by distant companions in the RV time series of HAT-P-2, HAT-P-22, and HAT-P-29; we significantly improved the orbital parameters of the long-period planet HAT-P-17c; and we revised the planetary parameters of CoRoT-1b, which turned out to be considerably more inflated than previously found.

Planck intermediate results: LI. features in the cosmic microwave background temperature power spectrum and shifts in cosmological parameters

Abstract: The six parameters of the standard ΛCDM model have best-fit values derived from the Planck temperature power spectrum that are shifted somewhat from the best-fit values derived from WMAP data. These shifts are driven by features in the Planck temperature power spectrum at angular scales that had never before been measured to cosmic-variance level precision. We have investigated these shifts to determine whether they are within the range of expectation and to understand their origin in the data. Taking our parameter set to be the optical depth of the reionized intergalactic medium τ, the baryon density ωb, the matter density ωm, the angular size of the sound horizon θ∗, the spectral index of the primordial power spectrum, ns, and Ase− 2τ (where As is the amplitude of the primordial power spectrum), we have examined the change in best-fit values between a WMAP-like large angular-scale data set (with multipole moment l 800, or splitting at a different multipole, yields similar results. We examined the l 800 power spectrum data and find that the features there that drive these shifts are a set of oscillations across a broad range of angular scales. Although they partly appear similar to the effects of enhanced gravitational lensing, the shifts in ΛCDM parameters that arise in response to these features correspond to model spectrum changes that are predominantly due to non-lensing effects; the only exception is τ, which, at fixed Ase− 2τ, affects the l> 800 temperature power spectrum solely through the associated change in As and the impact of that on the lensing potential power spectrum. We also ask, “what is it about the power spectrum at l < 800 that leads to somewhat different best-fit parameters than come from the full l range?” We find that if we discard the data at l < 30, where there is a roughly 2σ downward fluctuation in power relative to the model that best fits the full l range, the l < 800 best-fit parameters shift significantly towards the l < 2500 best-fit parameters. In contrast, including l < 30, this previously noted “low-l deficit” drives ns up and impacts parameters correlated with ns, such as ωm and H0. As expected, the l < 30 data have a much greater impact on the l < 800 best fit than on the l < 2500 best fit. So although the shifts are not very significant, we find that they can be understood through the combined effects of an oscillatory-like set of high-l residuals and the deficit in low-l power, excursions consistent with sample variance that happen to map onto changes in cosmological parameters. Finally, we examine agreement between PlanckTT data and two other CMB data sets, namely the Planck lensing reconstruction and the TT power spectrum measured by the South Pole Telescope, again finding a lack of convincing evidence of any significant deviations in parameters, suggesting that current CMB data sets give an internally consistent picture of the ΛCDM model.Key words: cosmology: observations / cosmic background radiation / cosmological parameters / cosmology: theory

Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars. VI. Age and abundance structure of the stellar populations in the central sub-kpc of the Milky Way

Abstract: We present a detailed elemental abundance study of 90 F and G dwarf, turn-off, and subgiant stars in the Galactic bulge. Based on high-resolution spectra acquired during gravitational microlensing events, stellar ages and abundances for 11 elements (Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Zn, Y and Ba) have been determined. Four main findings are presented: (1) a wide metallicity distribution with distinct peaks at [Fe/H] = -1.09, -0.63, -0.20, + 0.12, + 0.41; (2) a highfraction of intermediate-age to young stars where at [Fe/H] > 0 more than 35% are younger than 8 Gyr, and for [Fe/H] ≲-0.5 most stars are 10 Gyr or older; (3) several episodes of significant star formation in the bulge has been identified: 3, 6, 8, and 11 Gyr ago; (4) tentatively the "knee" in the α-element abundance trends of the sub-solar metallicity bulge is located at a slightly higher [Fe/H] than in the local thick disk. These findings show that the Galactic bulge has complex age and abundance properties that appear to be tightly connected to the main Galactic stellar populations. In particular, the peaks in the metallicity distribution, the star formation episodes, and the abundance trends, show similarities with the properties of the Galactic thin and thick disks. At the same time, the star formation rate appears to have been slightly faster in the bulge than in the local thick disk, which most likely is an indication of the denser stellar environment closer to the Galactic centre. There are also additional components not seen outside the bulge region, and that most likely can be associated with the Galactic bar. Our results strengthen the observational evidence that support the idea of a secular origin for the Galactic bulge, formed out of the other main Galactic stellar populations present in the central regions of our Galaxy. Additionally, our analysis of this enlarged sample suggests that the (V-I)0 colour of the bulge red clump should be revised to 1.09. (Less)

Delay-time distribution of core-collapse supernovae with late events resulting from binary interaction

Abstract: Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, 15+9-8%, of core-collapse supernovae are “late”, that is, they occur 50–200 Myr after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass (4–8 M⊙). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by 14+15-14% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that φ Persei is a likely progenitor and that eccentric neutron star – white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution.