Search

Circumstellar disks do not evolve in isolation, as about half of solar-type stars were born in binary or multiple systems. Resolving disks in binary systems provides the opportunity to examine the influence of stellar companions on the outcomes of planet formation. We aim to investigate and compare disks in stellar multiple systems with near-infrared scattered-light imaging as part of the Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS) program. We used polarimetric differential imaging with SPHERE/IRDIS at the VLT to search for scattered light from the circumstellar disks in three multiple systems, CHX 22, S CrA, and HP Cha. We performed astrometric and orbit analyses for the stellar companions using archival HST, VLT/NACO, and SPHERE data. Combined with the age and orbital constraints, the observed disk structures provide insights into the evolutionary history and the impact of the stellar companions. The small grains in CHX 22 form a tail-like structure surrounding the close binary, which likely results from a close encounter and capture of a cloudlet. S CrA shows intricate structures (tentative ringed and spiral features) in the circumprimary disk as a possible consequence of perturbations by companions. The circumsecondary disk is truncated and connected to the primary disk via a streamer, suggesting tidal interactions. In HP Cha, the primary disk is less disturbed and features a tenuous streamer, through which the material flows towards the companions. The comparison of the three systems spans a wide range of binary separation (50 - 500 au) and illustrates the decreasing influence on disk structures with the distance of companions. This agrees with the statistical analysis of exoplanet population in binaries, that planet formation is likely obstructed around close binary systems, while it is not suppressed in wide binaries., 19 pages, 6 figures, accpeted for publication in A&A

This is a correction for 2021 AJ 161 38

The physical and chemical conditions in Class 0/I protostars are fundamental in unlocking the protostellar accretion process and its impact on planet formation. The aim is to determine which physical components are traced by different molecules at sub-arcsecond scales (100 - 400 au). We use a suite of Atacama Large Millimeter/submillimeter Array (ALMA) datasets in Band 6 (1 mm), Band 5 (1.8 mm) and Band 3 (3 mm) at spatial resolutions 0.5 - 3" for 16 protostellar sources. The protostellar envelope is well traced by C$^{18}$O, DCO$^+$ and N$_2$D$^+$, with the freeze-out of CO governing the chemistry at envelope scales. Molecular outflows are seen in classical shock tracers like SiO and SO, but ice-mantle products such as CH$_3$OH and HNCO released with the shock are also observed. The molecular jet is prominent not only in SiO and SO but also occasionally in H$_2$CO. The cavity walls show tracers of UV-irradiation such as C$_2$H c-C$_3$H$_2$ and CN. The hot inner envelope, apart from showing emission from complex organic molecules (COMs), also presents compact emission from small molecules like H$_2$S, SO, OCS and H$^{13}$CN, most likely related to ice sublimation and high-temperature chemistry. Sub-arcsecond millimeter-wave observations allow to identify those (simple) molecules that best trace each of the physical components of a protostellar system. COMs are found both in the hot inner envelope (high excitation lines) and in the outflows (lower-excitation lines) with comparable abundances. COMs can coexist with hydrocarbons in the same protostellar sources, but they trace different components. In the near future, mid-IR observations with JWST-MIRI will provide complementary information about the hottest gas and the ice mantle content, at unprecedented sensitivity and at resolutions comparable to ALMA for the same sources., Comment: 38 pages, 33 figures, accepted for publication to A\&A

FU Orionis-type objects (FUors) are low-mass pre-main sequence stars undergoing a temporary, but significant increase of mass accretion rate from the circumstellar disk onto the protostar. It is not yet clear what triggers the accretion bursts and whether the disks of FUors are in any way different from disks of non-bursting young stellar objects. Motivated by this, we conducted a 1.3 mm continuum survey of ten FUors and FUor-like objects with ALMA, using both the 7 m array and the 12 m array in two different configurations to recover emission at the widest possible range of spatial scales. We detected all targeted sources and several nearby objects as well. To constrain the disk structure, we fit the data with models of increasing complexity from 2D Gaussian to radiative transfer, enabling comparison with other samples modeled in a similar way. The radiative transfer modeling gives disk masses that are significantly larger than what is obtained from the measured millimeter fluxes assuming optically thin emission, suggesting that the FUor disks are optically thick at this wavelength. In comparison with samples of regular Class II and Class I objects, the disks of FUors are typically a factor of 2.9-4.4 more massive and a factor of 1.5-4.7 smaller in size. A significant fraction of them (65-70%) may be gravitationally unstable., Comment: 34 pages, 8 tables, 29 figures, accepted for publication in the Astrophysical Journal Supplement Series

Context. Eta Carinae is a highly eccentric, massive binary system (semimajor axis ~15.5 au) with powerful stellar winds and a phase-dependent wind-wind collision (WWC) zone. The primary star, η Car A, is a luminous blue variable (LBV); the secondary, η Car B, is a Wolf-Rayet or O star with a faster but less dense wind. Aperture-synthesis imaging allows us to study the mass loss from the enigmatic LBV η Car. Understanding LBVs is a crucial step toward improving our knowledge about massive stars and their evolution. Aims. Our aim is to study the intensity distribution and kinematics of η Car’s WWC zone. Methods. Using the VLTI-MATISSE mid-infrared interferometry instrument, we perform Brα imaging of η Car’s distorted wind. Results. We present the first VLTI-MATISSE aperture-synthesis images of η Car A’s stellar windin several spectral channels distributed across the Brα 4.052 μm line (spectral resolving power R ~ 960). Our observations were performed close to periastron passage in February 2020 (orbital phase ~ 14.0022). The reconstructed iso-velocity images show the dependence of the primary stellar wind on wavelength or line-of-sight (LOS) velocity with a spatial resolution of 6 mas (~14 au). The radius of the faintest outer wind regions is ~26 mas (~60 au). At several negative LOS velocities, the primary stellar wind is less extended to the northwest than in other directions. This asymmetry is most likely caused by the WWC. Therefore, we see both the velocity field of the undisturbed primary wind and the WWC cavity. In continuum spectral channels, the primary star wind is more compact than in line channels. A fit of the observed continuum visibilities with the visibilities of a stellar wind CMFGEN model (CMFGEN is an atmosphere code developed to model the spectra of a variety of objects) provides a full width at half maximum fit diameter of the primary stellar wind of 2.84 ± 0.06 mas (6.54 ± 0.14 au). We comparethe derived intensity distributions with the CMFGEN stellar wind model and hydrodynamic WWC models.

Context. Carbon is one of the most abundant components in the Universe. While silicates have been the main focus of solid phase studies in protoplanetary discs (PPDs), little is known about the solid carbon content especially in the planet-forming regions (~0.1–10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as well as polycyclic aromatic hydrocarbons), which generate infrared (IR) features that can be used to trace the solid carbon reservoirs. The new mid-IR instrument MATISSE, installed at the Very Large Telescope Interferometer (VLTI), can spatially resolve the inner regions (~1–10 au) of PPDs and locate, down to the au-scale, the emission coming from carbon grains. Aims. Our aim is to provide a consistent view on the radial structure, down to the au-scale, as well as basic physical properties and the nature of the material responsible for the IR continuum emission in the inner disk region around HD 179218. Methods. We implemented a temperature-gradient model to interpret the disk IR continuum emission, based on a multiwavelength dataset comprising a broadband spectral energy distribution and VLTI H-, L-, and N-bands interferometric data obtained in low spectral resolution. Then, we added a ring-like component, representing the carbonaceous L-band features-emitting region, to assess its detectability in future higher spectral resolution observations employing mid-IR interferometry. Results. Our temperature-gradient model can consistently reproduce our dataset. We confirmed a spatially extended inner 10 au emission in H- and L-bands, with a homogeneously high temperature (~1700 K), which we associate with the presence of stochastically heated nano-grains. On the other hand, the N-band emitting region presents a ring-like geometry that starts at about 10 au with a temperature of 400 K. Moreover, the existing low resolution MATISSE data exclude the presence of aromatic carbon grains (i.e., producing the 3.3 μm feature) in close proximity tothe star (≲1 au). Future medium spectral resolution MATISSE data will confirm their presence at larger distances. Conclusions. Our best-fit model demonstrates the presence of two separated dust populations: nano-grains that dominate the near- to mid-IR emission in the inner 10 au region and larger grains that dominate the emission outward. The presence of such nano-grains in the highly irradiated inner 10 au region of HD 179218 requires a replenishment process. Considering the expected lifetime of carbon nano-grains from The Heterogeneous dust Evolution Model for Interstellar Solids (THEMIS model), the estimated disk accretion inflow of HD 179218 could significantly contribute to feed the inner 10 au region in nano-grains.Moreover, we also expect a local regeneration of those nano-grains by the photo-fragmentation of larger aggregates.

Context. Infrared spectroscopy of star and planet forming regions is at the dawn of a new age with the upcoming James Webb Space Telescope. In support of these observations, laboratory spectra are required to identify complex organic molecules in the ices that cover the dust grains in these regions. Aims. This study aims to provide reference spectra to firmly detect icy methyl formate in the different stages of star and planet forming regions. Methyl formate is mixed in astronomically relevant matrices, and the peak positions, FWHMs, and relative band intensities are characterized for different temperatures to provide an analytical tool for astronomers. Methods. Methyl formate is deposited at 15 K under high-vacuum conditions. Specifically, methyl formate is deposited pure and mixed with CO, H$_2$CO, CH$_3$OH, H$_2$O, and CO:H$_2$CO:CH$_3$OH combined. Throughout the experiment infrared spectra are acquired with a FTIR spectrometer in the range from 4000-500 cm$^{-1}$ (2.5-20 $\mu$m) at a spectral resolution of 0.5 cm$^{-1}$. Results. We present the characterization of five solid-state methyl formate vibrational modes in pure and astronomically relevant ice matrices. The five selected vibrational modes, namely the C=O stretch, C$-$O stretch, CH$_3$ rocking, O$-$CH$_3$ stretching, and OCO deformation, are best suited for a JWST identification of methyl formate. For each of these vibrational modes, and each of the mixtures the TvS heatmaps, peak position versus FWHM, and relative band intensities are given. Additionally, the acquired reference spectra of methyl formate are compared with Spitzer observations of HH 46. A tentative detection of methyl formate provides an upper limit to the column density of $1.7\times10^{17}$ cm$^{-2}$, corresponding to an upper limit relative to water of $\leq 2.2\%$ and $\leq 40\%$ with respect to methanol., Comment: Accepted for publication in A&A

We sincerely thank the referee (Erik Rosolowsky) for the careful reading and useful comments to improve our manuscript. We would also like to show our gratitude to him for the kind assistance with the usage of CPROPS in the early stages of this work. REM was supported by the ALMA Japan Research Grant of NAOJ ALMA Project, NAOJ-ALMA-222. DE was supported by JSPS KAKENHI grant number JP 17K14254. DE was supported by the ALMA Japan Research Grant of NAOJ ALMA Project, NAOJ-ALMA-0093. MINK was supported by JSPS KAKENHI grant number JP 15J04974. KK was supported by JSPS KAKENHI grant number JP17H06130 and the NAOJ ALMA Scientific Research grant number 2017-06B. SV acknowledges support by the research projects AYA2014-53506-P and AYA2017-84897-P from the Spanish Ministerio de Economia y Competitividad, and by the Consejer ' ia de Conocimiento, Investigaci ' on y Universidad, Junta de Andaluc ' ia (FQM108) and European Regional Development Fund (ERDF)". This study has been partially financed by the Consejer ' ia de Conocimiento, Investigaci ' on y Universidad, Junta de Andaluc ' ia and European Regional Development Fund (ERDF), ref. SOMM17/6105/UGR. Part of this work was achieved using the grant of Visiting Scholars Program supported by the Research Coordination Committee, National Astronomical Observatory of Japan (NAOJ), National Institutes ofNatural Sciences (NINS). SM would like to thank the Ministry of Science and Technology (MOST) of Taiwan, MOST 107-2119-M-001-020. This research has made use of NASA's Astrophysics Data System. This research has made use of Astropy, a community-developed core PYTHON (http://www.python.org) package for Astronomy (Astropy Collaboration 2013, 2018); IPYTHON (Perez & Granger 2007); MATPLOTLIB (Hunter 2007); APLPY, an open-source plotting package for PYTHON (Robitaille & Bressert 2012), and NUMPY (van derWalt, Colbert & Varoquaux 2011). Data analysis was in part carried out on the open use data analysis computer system at the Astronomy Data Center, ADC, of the National Astronomical Observatory of Japan. This research has made use of the NASA/IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00803.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic ofKorea), in cooperationwith theRepublic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The NationalRadio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc., We present the first census of giant molecular clouds (GMCs) complete down to 106M and within the inner 4 kpc of the nearest giant elliptical and powerful radio galaxy, Centaurus A. We identified 689 GMCs using CO(1–0) data with 1 arcsec spatial resolution (∼20 pc) and 2 kms−1 velocity resolution obtained with the Atacama Large Millimeter/submillimeter Array. The I(CO)-N(H2) conversion factor based on the virial method is XCO = (2 ± 1) × 1020 cm−2(K km s−1)−1 for the entire molecular disc, consistent with that of the discs of spiral galaxies including the Milky Way, and XCO = (5 ± 2) × 1020 cm−2(K km s−1)−1 for the circumnuclear disc (CND; within a galactocentric radius of 200 pc). We obtained the GMC mass spectrum distribution and find that the best truncated power-law fit for the whole molecular disc, with index γ −2.41 ± 0.02 and upper cut-off mass ∼1.3 × 107M , is also in agreement with that of nearby disc galaxies. A trend is found in the mass spectrum index from steep to shallow as we move to inner radii. Although the GMCs are in an elliptical galaxy, the general GMC properties in the molecular disc are as in spiral galaxies. However, in the CND, large offsets in the line-width-size scaling relations (∼0.3 dex higher than those in the GMCs in the molecular disc), a different XCO factor, and the shallowest GMC mass distribution shape (γ = −1.1 ± 0.2) all suggest that there the GMCs are most strongly affected by the presence of the active galactic nucleus and/or shear motions., ALMA Japan Research Grant of NAOJ ALMA Project NAOJ-ALMA-222 NAOJ-ALMA-0093, Ministry of Education, Culture, Sports, Science and Technology, Japan (MEXT) Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research (KAKENHI) JP 17K14254 JP 15J04974 JP17H06130, NAOJ ALMA Scientific Research grant 2017-06B, Spanish Government AYA2014-53506-P AYA2017-84897-P, Junta de Andalucia FQM108, European Commission SOMM17/6105/UGR, Research Coordination Committee, National Astronomical Observatory of Japan (NAOJ), National Institutes ofNatural Sciences (NINS), Ministry of Science and Technology, Taiwan MOST 107-2119-M-001-020 2013.1.00803.S

Quantifying the composition of the material in protoplanetary disks is paramount to determining the potential for exoplanetary systems to produce and support habitable environments. A key complex organic molecule (COM) to detect is methanol (CH3OH). CH3OH primarily forms at low temperatures via the hydrogenation of CO ice on the surface of icy dust grains and is a necessary basis for the formation of more complex species like amino acids and proteins. We report the detection of CH3OH in a disk around a young, luminous A-type star HD100546. This disk is warm and therefore does not host a significant CO ice reservoir. We argue that the CH3OH cannot form in situ, and hence, this disk has likely inherited COMs rich ice from an earlier cold dark cloud phase. This is strong evidence that at least some of the organic material survives the disk formation process and can then be incorporated into forming planets, moons and comets. Therefore, crucial pre-biotic chemical evolution already takes place in dark star-forming clouds., subm. to Nature Astronomy

Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. We use geometric models to fit the data. Our models include a smoothed ring, a flat disk with inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetries as well. We also perform numerical hydro-dynamical simulations of the inner edge of the disk. Results. Our modeling reveals a significant brightness asymmetry in the L-band disk emission. The brightness maximum of the asymmetry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio in the azimuthal variation is $3.5 \pm 0.2$. Comparing our result on the location of the asymmetry with other interferometric measurements, we confirm that the morphology of the $r, accepted for publication in A&A

Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resolve the disk down to scales of ~7 au. In addition to the new SPHERE observations, we utilize VLT/NACO, HST/STIS and ALMA archival data. The new SPHERE data show the disk around SU Aur and extended dust structures in unprecedented detail. We resolve several dust tails connected to the Keplerian disk. By comparison with ALMA data, we show that these dust tails represent material falling onto the disk. The disk itself shows an intricate spiral structure and a shadow lane, cast by an inner, misaligned disk component. Our observations suggest that SU Aur is undergoing late infall of material, which can explain the observed disk structures. SU Aur is the clearest observational example of this mechanism at work and demonstrates that late accretion events can still occur in the class II phase, thereby significantly affecting the evolution of circumstellar disks. Constraining the frequency of such events with additional observations will help determine whether this process is responsible for the spin-orbit misalignment in evolved exoplanet systems., 18 pages, 12 figures, published in ApJL on 18-02-2021

The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. We first use photometric observations in various bands and Spitzer Space Telescope spectroscopy to constrain the IR excess of $\ell$ Car. Then, we analyze the VLTI/MATISSE measurements at a specific phase of observation, in order to determine the flux contribution, the size and shape of the environment of the star in the L band. We finally test the hypothesis of a shell of ionized gas in order to model the IR excess. We report the first detection in the L band of a centro-symmetric extended emission around l Car, of about 1.7$R_\star$ in FWHM, producing an excess of about 7.0\% in this band. In the N band, there is no clear evidence for dust emission from VLTI/MATISSE correlated flux and Spitzer data. On the other side, the modeled shell of ionized gas implies a more compact CSE ($1.13\pm0.02\,R_\star$) and fainter (IR excess of 1\% in the L band). We provide new evidences for a compact CSE of $\ell$ Car and we demonstrate the capabilities of VLTI/MATISSE for determining common properties of CSEs. While the compact CSE of $\ell$ Car is probably of gaseous nature, the tested model of a shell of ionized gas is not able to simultaneously reproduce the IR excess and the interferometric observations. Further Galactic Cepheids observations with VLTI/MATISSE are necessary for determining the properties of CSEs, which may also depend on both the pulsation period and the evolutionary state of the stars., Comment: 13 pages, 8 figures, accepted in Astronomy and Astrophysics

Connecting the composition of planet-forming disks with that of gas giant exoplanet atmospheres, in particular through C/O ratios, is one of the key goals of disk chemistry. Small hydrocarbons like $\rm C_2H$ and $\rm C_3H_2$ have been identified as tracers of C/O, as they form abundantly under high C/O conditions. We present resolved $\rm C_3H_2$ observations from the TW Hya Rosetta Stone Project, a program designed to map the chemistry of common molecules at $15-20$ au resolution in the TW Hya disk. Augmented by archival data, these observations comprise the most extensive multi-line set for disks of both ortho and para spin isomers spanning a wide range of energies, $E_u=29-97$ K. We find the ortho-to-para ratio of $\rm C_3H_2$ is consistent with 3 throughout extent of the emission, and the total abundance of both $\rm C_3H_2$ isomers is $(7.5-10)\times10^{-11}$ per H atom, or $1-10$% of the previously published $\rm C_2H$ abundance in the same source. We find $\rm C_3H_2$ comes from a layer near the surface that extends no deeper than $z/r=0.25$. Our observations are consistent with substantial radial variation in gas-phase C/O in TW Hya, with a sharp increase outside $\sim30$ au. Even if we are not directly tracing the midplane, if planets accrete from the surface via, e.g., meridonial flows, then such a change should be imprinted on forming planets. Perhaps interestingly, the HR 8799 planetary system also shows an increasing gradient in its giant planets' atmospheric C/O ratios. While these stars are quite different, hydrocarbon rings in disks are common, and therefore our results are consistent with the young planets of HR 8799 still bearing the imprint of their parent disk's volatile chemistry., 15 pages, 8 figures, Accepted in ApJ

The Origins Space Telescope (Origins) is one of four science and technology definition studies selected by the National Aeronautics and Space Administration (NASA) in preparation of the 2020 Astronomy and Astrophysics Decadal survey in the US. Origins will trace the history of our origins from the time dust and heavy elements permanently altered the cosmic landscape to present-day life. It is designed to answer three major science questions: How do galaxies form stars, make metals, and grow their central supermassive black holes from reionization? How do the conditions for habitability develop during the process of planet formation? Do planets orbiting M-dwarf stars support life? Origins operates at mid- to far-infrared wavelengths from ~ 2.8 μm to 588 μm, and is more than 1000 times more sensitive than prior far-IR missions due to its cold (~ 4.5 K) aperture and state-of-the-art instruments.

Molecular D/H ratios are frequently used to probe the chemical past of Solar System volatiles. Yet it is unclear which parts of the Solar Nebula hosted an active deuterium fractionation chemistry. To address this question, we present 0".2-0".4 ALMA observations of DCO+ and DCN 2-1, 3-2 and 4-3 towards the nearby protoplanetary disk around TW Hya, taken as part of the TW Hya Rosetta Stone project, augmented with archival data. DCO+ is characterized by an excitation temperature of ~40 K across the 70 au radius pebble disk, indicative of emission from a warm, elevated molecular layer. Tentatively, DCN is present at even higher temperatures. Both DCO+ and DCN present substantial emission cavities in the inner disk, while in the outer disk the DCO+ and DCN morphologies diverge: most DCN emission originates from a narrow ring peaking around 30~au, with some additional diffuse DCN emission present at larger radii, while DCO+ is present in a broad structured ring that extends past the pebble disk. Based on parametric disk abundance models, these emission patterns can be explained by a near-constant DCN abundance exterior to the cavity, and an increasing DCO+ abundance with radius. There appears to be an active deuterium fractionation chemistry in multiple disk regions around TW Hya, but not in the cold planetesimal-forming midplane and in the inner disk. More observations are needed to explore whether deuterium fractionation is actually absent in these latter regions, and if its absence is a common feature, or something peculiar to the old TW Hya disk., Comment: Accepted for publication in The Astronomical Journal

The thermal structure of protoplanetary disks is a fundamental characteristic of the system that has wide reaching effects on disk evolution and planet formation. In this study, we constrain the 2D thermal structure of the protoplanetary disk TW Hya structure utilizing images of seven CO lines. This includes new ALMA observations of 12CO J=2-1 and C18O J=2-1 as well as archival ALMA observations of 12CO J=3-2, 13CO J=3-2, 6-5, C18O J= 3-2, 6-5. Additionally, we reproduce a Herschel observation of the HD J=1-0 line flux, the spectral energy distribution, and utilize a recent quantification of CO radial depletion in TW Hya. These observations were modeled using the thermochemical code RAC2D, and our best fit model reproduces all spatially resolved CO surface brightness profiles. The resulting thermal profile finds a disk mass of 0.025 Msun and a thin upper layer of gas depleted of small dust with a thickness of approx 1.2% of the corresponding radius. Using our final thermal structure, we find that CO alone is not a viable mass tracer as its abundance is degenerate with the total H2 surface density. Different mass models can readily match the spatially resolved CO line profiles with disparate abundance assumptions. Mass determination requires additional knowledge and, in this work, HD provides the additional constraint to derive the gas mass and supports the inference of CO depletion in the TW Hya disk. Our final thermal structure confirms the use of HD as a powerful probe of protoplanetary disk mass. Additionally, the method laid out in this paper is an employable strategy for extraction of disk temperatures and masses in the future., Comment: 28 pages, 12 main figures, 14 appendix figures. Accepted to ApJ

Context. The star formation rates and stellar densities found in young massive clusters suggest that these stellar systems originate from gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Until today, however, the physical characterization of this ultra high density material remains largely unconstrained in observations. Aims. We investigated the density properties of the star-forming gas in the OMC-1 region located in the vicinity of the Orion Nebula Cluster (ONC). Methods. We mapped the molecular emission at 652 GHz in OMC-1 as part of the APEX-SEPIA660 Early Science. Results. We detect bright and extended N$_2$H$^+$ (J=7-6) line emission along the entire OMC-1 region. Comparisons with previous ALMA data of the (J=1-0) transition and radiative transfer models indicate that the line intensities observed in this N$_2$H$^+$ (7-6) line are produced by large mass reservoirs of gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Conclusions. The first detection of this N$_2$H$^+$ (7-6) line at parsec-scales demonstrates the extreme density conditions of the star-forming gas in young massive clusters such as the ONC. Our results highlight the unique combination of sensitivity and mapping capabilities of the new SEPIA660 receiver for the study of the ISM properties at high frequencies., 8 pages, 5 Figures. Accepted for publication by A&A

Context: ALMA observations of Venus at 267 GHz have been presented in the literature that show the apparent presence of phosphine (PH3) in its atmosphere. Phosphine has currently no evident production routes on the planet's surface or in its atmosphere. Aims: The aim of this work is to assess the statistical reliability of the line detection by independent re-analysis of the ALMA data. Methods: The ALMA data were reduced as in the published study, following the provided scripts. First the spectral analysis presented in the study was reproduced and assessed. Subsequently, the spectrum was statistically evaluated, including its dependence on selected ALMA baselines. Results: We find that the 12th-order polynomial fit to the spectral passband utilised in the published study leads to spurious results. Following their recipe, five other >10 sigma lines can be produced in absorption or emission within 60 km/s from the PH3 1-0 transition frequency by suppressing the surrounding noise. Our independent analysis shows a feature near the PH3 frequency at a ~2 sigma level, below the common threshold for statistical significance. Since the spectral data have a non-Gaussian distribution, we consider a feature at such level as statistically unreliable that cannot be linked to a false positive probability. Conclusions: We find that the published 267-GHz ALMA data provide no statistical evidence for phosphine in the atmosphere of Venus., Accepted for publication in A&A after minor revisions

Context. The evolution of protoplanetary disks is dominated by the conservation of angular momentum, where the accretion of material onto the central star is fed by the viscous expansion of the outer disk or by disk winds extracting angular momentum without changing the disk size. Studying the time evolution of disk sizes therefore allows us to distinguish between viscous stresses or disk winds as the main mechanism of disk evolution. Observationally, estimates of the size of the gaseous disk are based on the extent of CO submillimeter rotational emission, which is also affected by the changing physical and chemical conditions in the disk during the evolution. Aims. We study how the gas outer radius measured from the extent of the CO emission changes with time in a viscously expanding disk. We also investigate to what degree this observable gas outer radius is a suitable tracer of viscous spreading and whether current observations are consistent with viscous evolution. Methods. For a set of observationally informed initial conditions we calculated the viscously evolved density structure at several disk ages and used the thermochemical code DALI to compute synthetic emission maps, from which we measured gas outer radii in a similar fashion as observations. Results. The gas outer radii (RCO, 90%) measured from our models match the expectations of a viscously spreading disk: RCO, 90% increases with time and, for a given time, RCO, 90% is larger for a disk with a higher viscosity αvisc. However, in the extreme case in which the disk mass is low (Mdisk ≤ 10−4 M⊙) and αvisc is high (≥10−2), RCO, 90% instead decreases with time as a result of CO photodissociation in the outer disk. For most disk ages, RCO, 90% is up to ~12× larger than the characteristic size Rc of the disk, and RCO, 90%/Rc is largest for the most massive disk. As a result of this difference, a simple conversion of RCO, 90% to αvisc overestimates the true αvisc of the disk by up to an order of magnitude. Based on our models, we find that most observed gas outer radii in Lupus can be explained using viscously evolving disks that start out small (Rc(t = 0) ≃ 10 AU) and have a low viscosity (αvisc = 10−4−10−3). Conclusions. Current observations are consistent with viscous evolution, but expanding the sample of observed gas disk sizes to star-forming regions, both younger and older, would better constrain the importance of viscous spreading during disk evolution.

Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures., 13 pages + 6 pages appendix, 7 figures, accepted for publication in aanda

Formaldehyde (H2CO) is an important precursor to organics like methanol (CH3OH). It is important to understand the conditions that produce H2CO and prebiotic molecules during star and planet formation. H2CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold (≲20 K) dust grains. To understand which pathway dominates, gaseous H2CO's ortho-to-para ratio (OPR) has been used as a probe, with a value of 3 indicating "warm" conditions and 12 cm−2 and 25 ± 2 K, (4.4 ± 0.3) × 1011 cm−2 for ortho- and para-H2CO, respectively, and an OPR of 2.49 ± 0.23. A radially resolved analysis shows that the observed H2CO emits mostly at rotational temperatures of 30–40 K, corresponding to a layer with z/R ≥ 0.25. The OPR is consistent with 3 within 60 au, the extent of the pebble disk, and decreases beyond 60 au to 2.0 ± 0.5. The latter corresponds to a spin temperature of 12 K, well below the rotational temperature. The combination of relatively uniform emitting conditions, a radial gradient in the OPR, and recent laboratory experiments and theory on OPR ratios after sublimation, led us to speculate that gas-phase formation is responsible for the observed H2CO across the TW Hya disk.

Recent ALMA surveys of protoplanetary disks have shown that for most disks the extent of the gas emission is greater than the extent of the thermal emission of the millimeter-sized dust. Both line optical depth and the combined effect of radially dependent grain growth and radial drift may contribute to this observed effect. For a sample of 10 disks from the Lupus survey we investigate how well dust-based models without radial dust evolution reproduce the observed 12CO outer radius, and determine whether radial dust evolution is required to match the observed gas-dust size difference. We used the thermochemical code DALI to obtain 12CO synthetic emission maps and measure gas and dust outer radii (Rco, Rmm) using the same methods as applied to the observations, which were compared to observations on a source-by-source basis. For 5 disks we find that the observed gas-dust size difference is larger than the gas-dust size difference due to optical depth, indicating that we need both dust evolution and optical depth effects to explain the observed gas-dust size difference. For the other 5 disks the observed gas-dust size difference can be explained using only line optical depth effects. We also identify 6 disks not included in our initial sample but part of a survey of the same star-forming region that show significant 12CO emission beyond 4 x Rmm. These disks, for which no Rco is available, likely have gas-dust size differences greater than 4 and are difficult to explain without substantial dust evolution. Our results suggest that radial drift and grain growth are common features among both bright and fain disks. The effects of radial drift and grain growth can be observed in disks where the dust and gas radii are significantly different, while more detailed models and deeper observations are needed to see this effect in disks with smaller differences., Comment: 17 pages, 11 figures, accepted in A&A

We present ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature., 21 pages, 11 figures, accepted by ApJ

Erratum to: A&A, 627, L6 (2019), https://doi.org/10.1051/0004-6361/201935633

The chemical composition of high-mass protostars reflects the physical evolution associated with different stages of star formation. In this study, the molecular inventory of the forming high-mass star AFGL 4176 is studied in detail at high angular resolution (~0.35 arcsec) using ALMA. This high resolution makes it possible to separate the emission associated with the inner hot envelope and disk around the forming star from that of its cool outer envelope. In addition, the high sensitivity of ALMA makes it possible to identify weak and optically thin lines and allows for many isotopologues to be detected, providing a more complete and accurate inventory of the source. A total of 23 different molecular species and their isotopologues are detected in the spectrum towards AFGL 4176. The most abundant species is methanol (CH3OH), remaining species are present at levels between 0.003 % and 15 % with respect to CH3OH. Hints that N-bearing species peak slightly closer to the location of the peak continuum emission than the O-bearing species are seen. AFGL 4176 comprises a rich chemical inventory including many complex species present on disk-scales. On average, the derived column density ratios with respect to methanol of O-bearing species are higher than those derived for N-bearing species by a factor of three. This may indicate that AFGL 4176 is a relatively young source since nitrogen chemistry generally takes longer to evolve in the gas-phase. Taking methanol as a reference, the composition of AFGL 4176 more closely resembles that of the low-mass protostar IRAS 16293-2422B than that of high-mass star-forming regions located near the Galactic centre. This similarity hints that the chemical composition of complex species is already set in the cold cloud stage and implies that AFGL 4176 is a young source whose chemical composition has not yet been strongly processed by the central protostar., Accepted for publication in Astronomy & Astrophysics

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales 150K, this does not preclude that the substructures formed earlier through disc fragmentation., 9 Total, inc references, appendix, long author list, 4 figures (inc appendix)

The cyanide radical CN is abundant in protoplanetary disks, with line fluxes often comparable to those of $^{13}$CO. It is known to be sensitive to UV irradiation of the upper disk atmosphere, with models predicting ring-shaped emission. We seek to characterize the CN emission from 94 Class II disks in the Lupus star forming region, compare it to observations in other regions, and interpret our observations with a grid of models. The CN emission morphology is discussed for two primordial disks, Sz 71 and Sz 98, and modeled in more detail. ALMA observed CN N=3-2 in Lupus disks down to sensitivities better than previous surveys. Models constructed with the physico-chemical code DALI are used to study the integrated fluxes of the disks and resolved emission of CN in disks without (dust) substructures. We find that CN N=3-2 is bright, and detected in 38% of sources, but its disk-integrated flux is not strongly correlated to either $^{13}$CO or continuum flux. Compared to pre-ALMA single-dish surveys, no significant difference in the CN flux distributions in Lupus and Taurus-Auriga is found, although $\rho$ Ophiuchus disks may be fainter on average. We find ring-shaped CN emission with peak radii of ~50 AU in two resolved disks. We conclude that a large fraction of sources is faint in CN; only exponential gas surface density cutoffs at $R_{c} \leq 15$ AU can reconcile observations with models. This is the first observational evidence of such a compact gas disk population in Lupus. Absolute intensities and the emission morphology of CN are reproduced by DALI models without the need for any continuum substructure; they are unrelated to the CO snowline location. These observations and the successful modeling of these rings provide a new probe of the structure and conditions in disks, and particularly their incident UV radiation field, if disk size is determined by the data., Comment: 22 pages, 14 figures, A&A accepted

Methanol (CH$_3$OH) is at the root of organic ice chemistry in protoplanetary disks. However, its weak emission has made detections difficult. To date, gas-phase CH$_3$OH has been detected in only one Class II disk, TW Hya. We use the Atacama Large Millimeter/submillimeter Array (ALMA) to search for a total of four CH$_3$OH emission lines in bands 6 and 7 toward the disk around the young Herbig Ae star HD 163296. The disk-averaged column density of methanol and its related species formaldehyde (H$_2$CO) are estimated assuming optically thin emission in local thermodynamic equilibrium. We compare these results to the gas-phase column densities of the TW Hya disk. No targeted methanol lines were detected individually nor after line stacking. The 3$��$ disk-integrated intensity upper limits are $< 51$ mJy km s$^{-1}$ for the band 6 lines and $< 26$ mJy km s$^{-1}$ for the band 7 lines. The band 7 lines provide the strictest 3$��$ upper limit on disk-averaged column density with $N_{\mathrm{avg}} < 5.0 \times 10^{11}$ cm$^{-2}$. The methanol-to-formaldehyde ratio is CH$_3$OH/H$_2$CO $< 0.24$ in the HD 163296 disk compared to a ratio of $1.27$ in the TW Hya disk. Differences in the stellar irradiation of Herbig disks compared to T Tauri disks likely influence the gaseous methanol and formaldehyde content. Possible reasons for the lower HD 163296 methanol-to-formaldehyde ratio include: a higher than expected gas-phase formation of H$_2$CO in the HD 163296 disk, uncertainties in the grain surface formation efficiency of CH$_3$OH and H$_2$CO, and differences in the disk structure and/or CH$_3$OH and H$_2$CO desorption processes that release the molecules from ice mantles back into the gas phase. These results provide observational evidence that the gas-phase chemical complexity found in disks may be strongly influenced by the spectral type of the host star., Accepted to Astronomy & Astrophysics, 11 pages, 3 figures

We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previously seen, up to a radius of 4 kpc. We find that the global star formation rate is $\sim$1 \Msol yr$^{-1}$, which yields a star formation efficiency (SFE) of 0.6 Gyr$^{-1}$ (depletion time $\tau =$1.5 Gyr), similar to those in disk galaxies. We show the most detailed view to date (40\,pc resolution) of the relation between molecular gas and star formation within the stellar component of an elliptical galaxy, from several kpc scale to the circumnuclear region close to the powerful radio jet. Although on average the SFEs are similar to those of spiral galaxies, the circumnuclear disk (CND) presents SFEs of 0.3 Gyr$^{-1}$, lower by a factor of 4 than the outer disk. The low SFE in the CND is in contrast to the high SFEs found in the literature for the circumnuclear regions of some nearby disk galaxies with nuclear activity, probably as a result of larger shear motions and longer AGN feedback. The higher SFEs in the outer disk suggests that only central molecular gas or filaments with sufficient density and strong shear motions will remain in $\sim$1 Gyr, which will later result in the compact molecular distributions and low SFEs usually seen in other giant ellipticals with cold gas., Comment: Minor update. 20 pages, 11 figures, accepted for publication in ApJ

Hydrogen deuteride (HD) rotational line emission can provide reliable protoplanetary disk gas mass measurements, but it is difficult to observe and detections have been limited to three T-Tauri disks. No new data have been available since the \emph{Herschel} Space Observatory mission ended in 2013. We set out to obtain new disk gas mass constraints by analysing upper limits on HD 1-0 emission in \emph{Herschel}/PACS archival data from the DIGIT key programme. With a focus on the Herbig Ae/Be disks, whose stars are more luminous than T Tauris, we determine upper limits for HD in data previosly analysed for its line detections. Their significance is studied with a grid of models run with the DALI physical-chemical code, customised to include deuterium chemistry. Nearly all the disks are constrained to $M_{\rm gas}\leq0.1\,$M$_{\odot}$, ruling out global gravitational instability. A strong constraint is obtained for the HD 163296 disk mass, $M_{\rm gas} \leq 0.067\,$M$_{\odot}$, implying $\Delta_{\rm g/d}\leq100$. This HD-based mass limit is towards the low end of CO-based mass estimates for the disk, highlighting the large uncertainty in using only CO and suggesting that gas-phase CO depletion in HD 163296 is at most a factor of a few. The $M_{\rm gas}$ limits for HD 163296 and HD 100546, both bright disks with massive candidate protoplanetary systems, suggest disk-to-planet mass conversion efficiencies of $M_{\rm p}/(M_{\rm gas} + M_{\rm p})\approx10$ to $40\,$% for present-day values. Near-future observations with SOFIA/HIRMES will be able to detect HD in the brightest Herbig~Ae/Be disks within $150\,$pc with $\approx10\,$h integration time., Comment: Accepted for publication in A&A

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation., Comment: Accepted in A&A

In the search for the building blocks of life, nitrogen-bearing molecules are of particular interest since nitrogen-containing bonds are essential for the linking of amino acids and ultimately the formation of larger biological structures. The elusive molecule methylamine (CH$_3$NH$_2$) is thought to be a key pre-biotic species but has so far only been securely detected in the giant molecular cloud Sgr B2. We identify CH$_3$NH$_2$ and other simple nitrogen-bearing species towards three hot cores in NGC 6334I. Column density ratios are derived in order to investigate the relevance of the individual species as precursors of biotic molecules. Observations obtained with ALMA were used to study transitions of CH$_3$NH$_2$, CH$_2$NH, NH$_2$CHO, and the $^{13}$C- and $^{15}$N-methyl cyanide (CH$_3$CN) isotopologues. Column densities are derived for each species assuming LTE and excitation temperatures in the range 220-340 K for CH$_3$NH$_2$, 70-110 K for the CH$_3$CN isotopologues, and 120-215 K for NH$_2$CHO and CH$_2$NH. We report the first detections of CH$_3$NH$_2$ towards NGC 6334I with column density ratios with respect to CH$_3$OH of 5.9$\times$10$^{-3}$, 1.5$\times$10$^{-3}$, and 5.4$\times$10$^{-4}$ for the three hot cores MM1, MM2, and MM3, respectively. These values are slightly lower than the values derived for Sgr B2 but higher by more than order of magnitude as compared with the values derived for the low-mass protostar IRAS 16293-2422B. The detections of CH$_3$NH$_2$ in the hot cores of NGC 6334I hint that CH$_3$NH$_2$ is generally common in the interstellar medium, albeit high-sensitivity observations are essential for its detection. The good agreement between model predictions of CH$_3$NH$_2$ ratios and the observations towards NGC 6334I indicate a main formation pathway via radical recombination on grain surfaces., Comment: Accepted for publication in A&A

Protoplanetary disks are thought to evolve viscously, where the disk mass - the reservoir available for planet formation - decreases over time as material is accreted onto the central star. Observations show a correlation between dust mass and the stellar accretion rate, as expected from viscous theory. However, the gas mass inferred from 13CO and C18O line fluxes, which should be a more direct measure, shows no such correlation. Using thermochemical DALI models, we investigate how 13CO and C18O J=3-2 line fluxes change over time in a viscously evolving disk. We also investigate if the chemical conversion of CO through grain-surface chemistry combined with viscous evolution can explain the observations of disks in Lupus. The 13CO and C18O 3-2 line fluxes increase over time due to their optically thick emitting regions growing in size as the disk expands viscously. The C18O 3-2 emission is optically thin throughout the disk for only a subset of our models (Mdisk (t = 1 Myr) < 1e-3 Msun). For these disks the integrated C18O flux decreases with time, similar to the disk mass. The C18O 3-2 fluxes for the bulk of the disks in Lupus (with Mdust < 5e-5 Msun) can be reproduced to within a factor of ~2 with viscously evolving disks in which CO is converted into other species through grain-surface chemistry driven by a cosmic-ray ionization rate zeta_cr ~ 5e-17 - 1e-16 s^-1. However, explaining the stacked C18O upper limits requires a lower average abundance than our models can produce and they cannot explain the observed 13CO fluxes, which, for most disks, are more than an order of magnitude fainter than what our models predict. Reconciling the 13CO fluxes of viscously evolving disks with the observations requires either a combination of efficient vertical mixing and a high zeta_cr or low mass disks (Mdust < 3e-5 Msun) being much thinner and/or smaller than their more massive counterparts., 21 pages, 14 figures, accepted in A&A

MATISSE is the second-generation mid-infrared spectrograph and imager for the Very Large Telescope Interferometer (VLTI) at Paranal. This new interferometric instrument will allow significant advances in various fundamental research fields: studying the planet-forming region of disks around young stellar objects, understanding the surface structures and mass loss phenomena affecting evolved stars, and probing the environments of black holes in active galactic nuclei. As a first breakthrough, MATISSE will enlarge the spectral domain of current optical interferometers by offering the L and M bands in addition to the N band. This will open a wide wavelength domain, ranging from 2.8 to 13 μm, exploring angular scales as small as 3 mas (L band) / 10 mas (N band). As a second breakthrough, MATISSE will allow mid-infrared imaging - closure-phase aperture-synthesis imaging - with the four Unit Telescopes (UT) or Auxiliary Telescopes (AT) of the VLTI. Moreover, MATISSE will offer a spectral resolution range from R ~ 30 to R ~ 5000. Here, we remind the concept, the instrumental design, and the main features of MATISSE. We also describe the last months of preparation, the status of the instrument, which was shipped to Cerro Paranal on the site of the ESO Very Large Telescope in October 2017, and the expected schedule for the opening to the community. The instrument is currently in its Commissioning phase. A complementary dedicated article details the Commissioning results, which include the first performance estimates on sky.

Transition disks with large dust cavities around young stars are promising targets for studying planet formation. Previous studies have revealed the presence of gas cavities inside the dust cavities hinting at recently formed, giant planets. However, many of these studies are biased towards the brightest disks in the nearby star forming regions, and it is not possible to derive reliable statistics that can be compared with exoplanet populations. We present the analysis of 11 transition disks with large cavities (>20 AU radius) from a complete disk survey of the Lupus star forming region, using ALMA Band 7 observations at 0.3" (22-30 AU radius) resolution of the 345 GHz continuum, 13CO and C18O 3-2 observations and the Spectral Energy Distribution of each source. Gas and dust surface density profiles are derived using the physical-chemical modeling code DALI. This is the first study of transition disks of large cavities within a complete disk survey within a star forming region. The dust cavity sizes range from 20-90 AU radius and in three cases, a gas cavity is resolved as well. The deep drops in gas density and large dust cavity sizes are consistent with clearing by giant planets. The fraction of transition disks with large cavities in Lupus is ~11%, which is inconsistent with exoplanet population studies of giant planets at wide orbits. Furthermore, we present a hypothesis of an evolutionary path for large massive disks evolving into transition disks with large cavities., 29 pages, 15 figures, Accepted by ApJ

The abundance of deuterated molecules in a star-forming region is sensitive to the environment in which they are formed. Deuteration fractions therefore provide a powerful tool for studying the physical and chemical evolution of a star-forming system. While local low-mass star-forming regions show very high deuteration ratios, much lower fractions are observed towards Orion and the Galactic Centre. We derive methanol deuteration fractions at a number of locations towards the high-mass star-forming region NGC 6334I, located at a mean distance of 1.3 kpc, and discuss how these can shed light on the conditions prevailing during its formation. We use high sensitivity, high spatial and spectral resolution observations obtained with ALMA to study transitions of the less abundant, optically thin, methanol-isotopologues: (13)CH3OH, CH3(18)OH, CH2DOH and CH3OD, detected towards NGC 6334I. Assuming LTE and excitation temperatures of 120-330 K, we derive column densities for each of the species and use these to infer CH2DOH/CH3OH and CH3OD/CH3OH fractions. Interestingly, the column densities of CH3OD are consistently higher than those of CH2DOH throughout the region. All regions studied in this work show CH2DOH/CH3OH as well as CH2DOH/CH3OD ratios that are considerably lower than those derived towards low-mass star-forming regions and slightly lower than those derived for the high-mass star-forming regions in Orion and the Galactic Centre. The low ratios indicate a grain surface temperature during formation ~30 K, for which the efficiency of the formation of deuterated species is significantly reduced., Comment: Accepted for publication in Astronomy & Astrophysics

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency., Comment: 39 pages (54 including the Appendix), 9 figures, 3 tables, in press in ApJ

VLA 1623$-$2417 is a triple protostellar system deeply embedded in Ophiuchus A. Sources A and B have a separation of 1.1", making their study difficult beyond the submillimeter regime. Lack of circumstellar gas emission suggested that VLA 1623$-$2417 B has a very cold envelope and is much younger than source A, generally considered the prototypical Class 0 source. We explore the consequences of new ALMA Band 9 data on the spectral energy distribution (SED) of VLA 1623$-$2417 and their inferred nature. Using dust continuum observations spanning from centimeter to near-infrared wavelengths, the SED of each component in VLA 1623$-$2417 is constructed and analysed. The ALMA Band 9 data presented here show that the SED of VLA 1623$-$2417 B does not peak at 850 $\mu$m as previously expected, but instead presents the same shape as VLA 1623$-$2417 A at wavelengths shorter than 450 $\mu$m. The results presented here indicate that the previous assumption that the flux in $Herschel$ and Spitzer observations is solely dominated by VLA 1623$-$2417 A is not valid, and instead, VLA 1623$-$2417 B most likely contributes a significant fraction of the flux at $\lambda~, Comment: 6 pages, 2 figures, accepted to A&A letters to the editor

We present ALMA Band 6 observations of a complete sample of protoplanetary disks in the young (1-3 Myr) Lupus star-forming region, covering the 1.33 mm continuum and the 12CO, 13CO, and C18O J=2-1 lines. The spatial resolution is 0.25 arcsec with a medium 3-sigma continuum sensitivity of 0.30 mJy, corresponding to M_dust ~ 0.2 M_earth. We apply "Keplerian masking" to enhance the signal-to-noise ratios of our 12CO zero-moment maps, enabling measurements of gas disk radii for 22 Lupus disks; we find that gas disks are universally larger than mm dust disks by a factor of two on average, likely due to a combination of the optically thick gas emission as well as the growth and inward drift of the dust. Using the gas disk radii, we calculate the dimensionless viscosity parameter, alpha_visc, finding a broad distribution and no correlations with other disk or stellar parameters, suggesting that viscous processes have not yet established quasi-steady states in Lupus disks. By combining our 1.33 mm continuum fluxes with our previous 890 micron continuum observations, we also calculate the mm spectral index, alpha_mm, for 70 Lupus disks; we find an anti-correlation between alpha_mm and mm flux for low-mass disks (M_dust < 5), followed by a flattening as disks approach alpha_mm = 2, which could indicate faster grain growth in higher-mass disks, but may also reflect their larger optically thick components. In sum, this work demonstrates the continuous stream of new insights into disk evolution and planet formation that can be gleaned from unbiased ALMA disk surveys., 22 pages, 10 figures; accepted to ApJ

This work aims to understand which midplane conditions are probed by the DCO$^+$ emission in the disk around the Herbig Ae star HD 169142. We explore the sensitivity of the DCO$^+$ formation pathways to the gas temperature and the CO abundance. The DCO$^+$ $J$=3-2 transition was observed with ALMA at a spatial resolution of 0.3". The HD 169142 DCO$^+$ radial intensity profile reveals a warm, inner component at radii 120 AU. The warm deuterium fractionation pathway is implemented as a constant DCO$^+$ abundance of 2.0$\times$10$^{-12}$ between 30-70 K. The DCO$^+$ emission probes a reservoir of cold material in the HD 169142 outer disk that is not revealed by the millimeter continuum, the SED, nor the emission from the 12CO, 13CO, or C18O $J$=2-1 lines., Accepted to Astronomy & Astrophysics

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources., Comment: 23 pages, 13 fig (inc 3 Appendix figures)

We aim to reproduce the DCO$^+$ emission in the disk around HD163296 using a simple 2D chemical model for the formation of DCO$^+$ through the cold deuteration channel and a parametric treatment of the warm deuteration channel. We use data from ALMA in band 6 to obtain a resolved spectral imaging data cube of the DCO$^+$ $J$=3--2 line in HD163296 with a synthesized beam of 0."53$\times$ 0."42. We adopt a physical structure of the disk from the literature that reproduces the spectral energy distribution. We then apply a simplified chemical network for the formation of DCO$^+$ that uses the physical structure of the disk as parameters along with a CO abundance profile, a constant HD abundance and a constant ionization rate. Finally, from the resulting DCO$^+$ abundances, we calculate the non-LTE emission using the 3D radiative transfer code LIME. The observed DCO$^+$ emission is reproduced by a model with cold deuteration producing abundances up to $1.6\times 10^{-11}$. Warm deuteration, at a constant abundance of $3.2\times 10^{-12}$, becomes fully effective below 32 K and tapers off at higher temperatures, reproducing the lack of DCO$^+$ inside 90 AU. Throughout the DCO$^+$ emitting zone a CO abundance of $2\times 10^{-7}$ is found, with $\sim$99\% of it frozen out below 19 K. At radii where both cold and warm deuteration are active, warm deuteration contributes up to 20\% of DCO$^+$, consistent with detailed chemical models. The decrease of DCO$^+$ at large radii is attributed to a temperature inversion at 250 AU, which raises temperatures above values where cold deuteration operates. Increased photodesorption may also limit the radial extent of DCO$^+$. The corresponding return of the DCO$^+$ layer to the midplane, together with a radially increasing ionization fraction, reproduces the local DCO$^+$ emission maximum at $\sim$260 AU., Comment: 9 pages, 5 figures, accepted 7th July 2018

Context. Deuterium fractionation has been used to study the thermalhistory of prestellar environments. Their formation pathways tracedifferent regions of the disk and may shed light into the physicalstructure of the disk, including locations of important features forplanetary formation. Aims: We aim to constrain the radial extentof the main deuterated species; we are particularly interested inspatially characterizing the high and low temperature pathways forenhancing deuteration of these species. Methods: We observed thedisk surrounding the Herbig Ae star HD 163296 using ALMA in Band 6 andobtained resolved spectral imaging data of DCO+ (J = 3 - 2),DCN (J = 3 - 2) and N2D+ (J = 3 - 2) withsynthesized beam sizes of 0.̋53 × 0.̋42, 0.̋53× 0.̋42, and 0.̋50 × 0.̋39, respectively. Weadopted a physical model of the disk from the literature and use the 3Dradiative transfer code LIME to estimate an excitation temperatureprofile for our detected lines. We modeled the radial emission profilesof DCO+, DCN, and N2D+, assuming theiremission is optically thin, using a parametric model of their abundancesand our excitation temperature estimates. Results:DCO+ can be described by a three-region model withconstant-abundance rings centered at 70 AU, 150 AU, and 260 AU. The DCNradial profile peaks at about 60 AU and N2D+ isseen in a ring at 160 AU. Simple models of both molecules using constantabundances reproduce the data. Assuming reasonable average excitationtemperatures for the whole disk, their disk-averaged column densities(and deuterium fractionation ratios) are 1.6-2.6×1012cm-2 (0.04-0.07), 2.9-5.2×1012cm-2 ( 0.02), and 1.6-2.5×1011cm-2 (0.34-0.45) for DCO+, DCN, andN2D+, respectively. Conclusions: Our simplebest-fit models show a correlation between the radial location of thefirst two rings in DCO+ and the DCN andN2D+ abundance distributions that can beinterpreted as the high and low temperature deuteration pathwaysregimes. The origin of the third DCO+ ring at 260 AU isunknown but may be due to a local decrease of ultraviolet opacityallowing the photodesorption of CO or due to thermal desorption of CO asa consequence of radial drift and settlement of dust grains. The derivedDf values agree with previous estimates of 0.05 forDCO+/HCO+ and 0.02 for DCN/HCN in HD 163296, and0.3-0.5 for N2D+/N2H+ in AS209, a T Tauri disk. The highN2D+/N2H+ confirmsN2D+ as a good candidate for tracing ionization inthe cold outer disk.The reduced images (FITS files) are only available at the CDS viaanonymous ftp to http://cdsarc.u-strasbg.fr (http://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/606/A125

Three formaldehyde lines were observed (H$_2$CO 3$_{03}$--2$_{02}$, H$_2$CO 3$_{22}$--2$_{21}$, and H$_2$CO 3$_{21}$--2$_{20}$) in the protoplanetary disk around the Herbig Ae star HD 163296 with ALMA at 0.5 arcsecond (60 AU) spatial resolution. H$_2$CO 3$_{03}$--2$_{02}$ was readily detected via imaging, while the weaker H$_2$CO 3$_{22}$--2$_{21}$ and H$_2$CO 3$_{21}$--2$_{20}$ lines required matched filter analysis to detect. H$_2$CO is present throughout most of the gaseous disk, extending out to 550 AU. An apparent 50 AU inner radius of the H$_2$CO emission is likely caused by an optically thick dust continuum. The H$_2$CO radial intensity profile shows a peak at 100 AU and a secondary bump at around 300 AU, suggesting increased production in the outer disk. Different parameterizations of the H$_2$CO abundance were compared to the observed visibilities with $��^2$ minimization, using either a characteristic temperature, a characteristic radius or a radial power law index to describe the H$_2$CO chemistry. Similar models were applied to ALMA Science Verification data of C$^{18}$O. In all modeling scenarios, fits to the H$_2$CO data show an increased abundance in the outer disk. The overall best-fit H$_2$CO model shows a factor of two enhancement beyond a radius of 270$\pm$20 AU, with an inner abundance of $2\!-\!5 \times 10^{-12}$. The H$_2$CO emitting region has a lower limit on the kinetic temperature of $T > 20$ K. The C$^{18}$O modeling suggests an order of magnitude depletion in the outer disk and an abundance of $4\!-\!12 \times 10^{-8}$ in the inner disk. The increase in H$_2$CO outer disk emission could be a result of hydrogenation of CO ices on dust grains that are then sublimated via thermal desorption or UV photodesorption, or more efficient gas-phase production beyond about 300 AU if CO is photodisocciated in this region.

This paper presents observations of the Lupus I molecular cloud at 450 and 850 μm with\ud Submillimetre Common User Bolometer Array (SCUBA-2) as part of the James Clerk Maxwell\ud Telescope Gould Belt Survey (JCMT GBS). Nine compact sources, assumed to be the discs of\ud young stellar objects (YSOs), 12 extended protostellar, pre-stellar and starless cores, and one\ud isolated, low-luminosity protostar, are detected in the region. Spectral energy distributions,\ud including submillimetre fluxes, are produced for 15 YSOs, and each is fitted with the models\ud of Robitaille et al. The proportion of Class 0/I protostars is higher than that seen in other\ud Gould Belt regions such as Ophiuchus and Serpens. Circumstellar disc masses are calculated\ud for more evolved sources, while protostellar envelope masses are calculated for protostars.\ud Up to four very low luminosity objects are found; a large fraction when compared to other\ud Spitzer c2d regions. One YSO has a disc mass greater than the minimum mass solar nebula.\ud 12 starless/protostellar cores are detected by SCUBA-2 and their masses are calculated. The\ud stability of these cores is examined using both the thermal Jeans mass and a turbulent virial\ud mass when possible. Two cores in Lupus I are super-Jeans and contain no known YSOs. One\ud of these cores has a virial parameter of 1.1 ± 0.4, and could therefore be pre-stellar. The high\ud ratio of Class 0/I to Class III YSOs (1:1), and the presence of a pre-stellar core candidate,\ud provides support for the hypothesis that a shock recently triggered star formation in Lupus I.

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that if modelled as thermally supported Bonnor–Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival 13CO velocity dispersion measurements and find that our cores are typically pressure confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields., Monthly Notices of the Royal Astronomical Society, 464 (4), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

The protoplanetary system HD 169142 is one of the few cases where a potential candidate protoplanet has been recently detected via direct imaging. 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 C$^{18}$O $J=2-1$ emission from the system HD 169142 at angular resolution of $\sim 0".18 - 0".28$ ($\sim 20\,$au$ - 33\,$au). The dust continuum emission reveals a double-ring structure with an inner ring between $0".17-0".28$ ($\sim 20 - 35\,$au) and an outer ring between $0".48-0".64$ ($\sim 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 \lesssim 20\,$au) or within the dust gap ($\sim 35 - 56\,$au). In contrast, the channel maps of the $J=2-1$ line of the three CO isotopologues reveal the presence of gas inside the dust cavity and dust gap. The gaseous disk is also much larger than the compact dust emission extending to $\sim 1'.5$ ($\sim 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 ($>$ \micron-size). Using the thermo-chemical disk code \textsc{dali}, the continuum and the CO isotopologues emission are modelled to quantitatively measure the gas and dust surface densities. The resulting gas surface density is reduced by a factor of $\sim 30-40$ inward of the dust gap. The gas and dust distribution hint at the presence of multiple planets shaping the disk structure via dynamical clearing (dust cavity and gap) and dust trapping (double ring dust distribution)., Accepted for publication in A&A

FU Orionis-type objects are young, low-mass stars with large outbursts in visible light that last for several years or decades. They are thought to represent an evolutionary phase during the life of every young star when accretion from the circumstellar disk is enhanced during recurring time periods. These outbursts are able to rapidly build up the star while affecting the circumstellar disk and thus the ongoing or future planet formation. In many models infall from a circumstellar envelope seems to be necessary to trigger the outbursts. We observed the J=1$-$0 rotational transition of $^{13}$CO and C$^{18}$O towards eight northern FU Orionis-type stars (V1057 Cyg, V1515 Cyg, V2492 Cyg, V2493 Cyg, V1735 Cyg, V733 Cep, RNO 1B and RNO 1C) and derive temperatures and envelope masses and discuss the morphology and kinematics of the circumstellar material. We detected extended CO emission associated with all our targets. Smaller scale CO clumps were found to be associated with five objects with radii of 2000$-$5000 AU and masses of 0.02$-$0.5 $M_{\odot}$; these are clearly heated by the central stars. Three of these envelopes are also strongly detected in the 2.7 mm continuum. No central CO clumps were detected around V733 Cep and V710 Cas but there are many other clumps in their environments. Traces of outflow activity were observed towards V1735 Cyg, V733 Cep and V710 Cas. The diversity of the observed envelopes enables us to set up an evolutionary sequence between the objects. We find their evolutionary state to range from early, embedded Class I stage to late, Class II-type objects with very low-mass circumstellar material. The results reinforce the idea of FU Orionis-type stars as representatives of a transitory stage between embedded Class I young stellar objects and classical T-Tauri stars., 17 pages, 11 figures; accepted for publication in A&A

An era has started in which gas and dust can be observed independently in protoplanetary disks, thanks to the recent surveys with 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/dust ratios. The goal of this work is to fully exploit CO isotopologues observations in Lupus, comparing them with physical-chemical model results, in order to obtain gas masses for a large number of disks. We have employed physical-chemical models to analyze continuum and CO isotopologues observations of Lupus disks, including isotope-selective processes and freeze-out. Employing also the ALMA 13CO-only detections, disk gas masses have been calculated for a total of 34 sources, expanding the sample of 10 disks studied by Ansdell et al. (2016), where also C18O was detected. We confirm that overall gas-masses are very low, often smaller than 1 $M_{\rm J}$, if volatile carbon is not depleted. Accordingly, global gas/dust ratios predominantly between 1 and 10. Low CO-based gas masses and gas/dust ratios may indicate rapid loss of gas, or alternatively chemical evolution, e.g. via sequestering of carbon from CO to more complex molecules, or carbon locked up in larger bodies. Current ALMA observations 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. Assuming that all Lupus disks have evolved mainly due to viscous processes over the past few Myr, the observed correlation between the current mass accretion rate and dust mass found by Manara et al. (2016) implies a constant gas-to-dust ratio, which is close to 100 based on the observed $M_{\rm disk}/\dot{M}_{\rm acc}$ ratio. This in turn points to a scenario in which carbon depletion is responsible for the low CO isotopologue line luminosities., 11 pages, 9 figures, A&A accepted