Search

We are performing a line survey of 8 planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from >20 C-, O,- N-, and S-bearing species. The observations in four spectral settings at 210-280 GHz with $1\sigma$ rms sensitivity of ~ 8-12 mJy/beam at 0.9" and 0.3 km/s resolution will be completed in 2024. The uv-visibilities are fitted with the DiskFit model to obtain key stellar and disk properties. In this paper, the combined $^{12}$CO, $^{13}$CO and C$^{18}$O $J = 2-1$ data are presented. We find that the CO fluxes and disk masses inferred from dust continuum tentatively correlate with the CO emission sizes. We constrain dynamical stellar masses, geometries, temperatures, the CO column densities and gas masses for each disk. The best-fit temperatures at 100 au are ~ 17-37 K, and decrease radially with the power-law exponent q ~ 0.05-0.76. The inferred CO column densities decrease radially with the power-law exponent p ~ 0.2-3.1. The gas masses estimated from $^{13}$CO (2-1) are ~ $0.001-0.2 M_\textrm{Sun}$. The best-fit CO column densities point to severe CO freeze-out in the disks. The DL Tau disk is an outlier, and has either stronger CO depletion or lower gas mass than the rest of the sample. The CO isotopologue ratios are roughly consistent with the observed values in disks and the low-mass star-forming regions., Comment: 47 pages, 10 figures, 6 tables, accepted for publication by A&A on 22/02/2024

Characterizing the dust thermal structure in protoplanetary disks is a fundamental task as the dust surface temperature can affect both the planetary formation and the chemical evolution. Since the temperature is dependent on many parameters, including the grain size, properly modeling the grain temperature structure can be challenging. Many chemistry disk models usually employ a sophisticated single dust structure designed to reproduce the effect of a realistic population presumably composed of a large diversity of sizes. This generally represents a good approximation in most cases. Nonetheless, this dilutes the effects of the complex radiative interactions between the different grain populations on the resulting dust temperature, and thus the chemistry. We seek to show that the radiative interactions between dust grains of different sizes can induce a non-trivial dust temperature structure that cannot be reproduced by a single dust population and that can significantly affect the chemical outcome. The disk thermal structures are computed using the Monte-Carlo radiative transfer code RADMC-3D. The thermal structures are post-processed using the gas-grain code NAUTILUS to calculate the evolution of the chemical abundance. We find that simultaneously using at least two independent dust grain populations in disk models produces a complex temperature structure due to the starlight intercepted by the upper layers of the disk. In particular, we find that micron-sized dust grains are warmer than larger grains and can even show a radial temperature bump in some conditions. This dust temperature spread between the grains populations results in the segregation of the CO snowline and the presence of an unexpected CO gas hole along the midplane. We compare the results with observed close to edge-on class I/II disks., Comment: 21 pages, 19 figures

The disk-outflow connection plays a key role in extracting excess angular momentum from a forming protostar. We have previously reported the discovery of a small molecular outflow from the edge-on T Tauri star in the Bok globule CB26 that shows a peculiar velocity pattern, reminiscent of an outflow that corotates with the disk. We report new, high-resolution mm-interferometric observations of CB26 with the aim of revealing the morphology and kinematics of the outflow at the disk-outflow interface. The IRAM PdBI was used to observe CO(2-1) at 1.3mm with a resolution of 0.5". Using a physical model of the disk, which was derived from the dust emission, we employed chemo-dynamical modeling combined with line radiative transfer to constrain kinematic parameters and to construct a model of the CO emission from the disk that allowed us to separate the emission of the disk from that of the outflow. Our observations confirm the disk-wind nature of the rotating molecular outflow from CB26. The new high-resolution data reveal an X-shaped morphology of the CO emission close to the disk, and vertical streaks extending from the disk surface with a small half-opening angle of ~7deg, which can be traced out to vertical heights of ~500au. We interpret this emission as the combination of the disk atmosphere and a well-collimated disk wind, which we trace down to vertical heights of 40au, where it is launched from the surface of the flared disk at radii of 20-45au. The observed CO outflow has a total momentum flux of 1e-5 Msun km/s/yr, which is nearly three orders of magnitude larger than the maximum thrust that can be provided by the luminosity of the central star. We conclude that photoevaporation cannot be the main driving mechanism for this outflow, but it must be predominantly an MHD disk wind. It is thus far the best-resolved rotating disk wind observed to be launched from a circumstellar disk., Comment: Accepted by A&A, 25 pages, 19 figures

Context. In the past few years, there has been a rise in the detection of streamers, asymmetric flows of material directed toward the protostellar disk with material from outside the star's natal core. It is unclear how they affect the process of mass accretion, in particular beyond the Class 0 phase. Aims. We investigate the gas kinematics around Per-emb-50, a Class I source in the crowded star-forming region NGC 1333. Our goal is to study how the mass infall proceeds from envelope to disk scales in this source. Results. We discover a streamer delivering material toward Per-emb-50 in H$_2$CO and C$^{18}$O emission. The streamer's emission can be well described by the analytic solutions for an infalling parcel of gas along a streamline with conserved angular momentum, both in the image plane and along the line of sight velocities. The streamer has a mean infall rate of $1.3 \times 10^{ -6}$ M$_{ \odot}$ yr$^{ -1}$, $5 -10$ times higher than the current accretion rate of the protostar. SO and SO$_2$ emission reveal asymmetric infall motions in the inner envelope, additional to the streamer around Per-emb-50. Furthermore, the presence of SO$_2$ could mark the impact zone of the infalling material. Conclusions. The streamer delivers sufficient mass to sustain the protostellar accretion rate and might produce an accretion burst, which would explain the protostar's high luminosity with respect to other Class I sources. Our results highlight the importance of late infall for protostellar evolution: streamers might provide a significant amount of mass for stellar accretion after the Class 0 phase., Comment: 20 pages, 14 figures, accepted for publication in A&A

Molecular line surveys are among the main tools to probe the structure and physical conditions in protoplanetary disks (PPDs), the birthplace of planets. The large radial and vertical temperature as well as density gradients in these PPDs lead to a complex chemical composition, making chemistry an important step to understand the variety of planetary systems. We aimed to study the chemical content of the protoplanetary disk surrounding GG Tau A, a well-known triple T Tauri system. We used NOEMA with the new correlator PolyFix to observe rotational lines at 2.6 to 4 mm from a few dozen molecules. We analysed the data with a radiative transfer code to derive molecular densities and the abundance relative to $^{13}$CO, which we compare to those of the TMC1 cloud and LkCa15 disk. We report the first detection of CCS in PPDs. We also marginally detect OCS and find 16 other molecules in the GG Tauri outer disk. Ten of them had been found previously, while seven others ($^{13}$CN, N$_2$H$^+$, HNC, DNC, HC$_3$N, CCS, and C$^{34}$S) are new detections in this disk. The analysis confirms that sulphur chemistry is not yet properly understood. The D/H ratio, derived from DCO$^{+}$/HCO$^{+}$, DCN/HCN, and DNC/HNC ratios, points towards a low temperature chemistry. The detection of the rare species CCS confirms that GG Tau is a good laboratory to study the protoplanetary disk chemistry, thanks to its large disk size and mass., Comment: Accepted by Astronomy and Astrophysics

Grain surface chemistry is key to the composition of protoplanetary disks around young stars. The temperature of grains depends on their size. We evaluate the impact of this temperature dependence on the disk chemistry. We model a moderately massive disk with 16 different grain sizes. We use POLARIS to calculate the dust grain temperatures and the local UV flux. We model the chemistry using the 3-phase astrochemical code NAUTILUS. Photoprocesses are handled using frequency-dependent cross-sections, and a new method to account for self and mutual shielding. The multi-grain model outputs are compared to those of single-grain size models (0.1 $\mu$m), with two different assumptions for their equivalent temperature. We find that the Langmuir-Hinshelwood (LH) mechanism at equilibrium temperature is not efficient to form H$_2$ at 3-4 scale heights ($H$), and adopt a parametric fit to a stochastic method to model H$_2$ formation instead. We find the molecular layer composition (1-3 $H$) to depend on the amount of remaining H atoms. Differences in molecular surface densities between single and multi-grain models are mostly due to what occurs above 1.5 $H$. At 100 au, models with colder grains produce H$_2$O and CH$_4$ ices in the midplane, and warmer ones produce more CO$_2$ ices, both allowing efficient depletion of C and O as soon as CO sticks on grain surfaces. Complex organic molecules (COMs) production is enhanced by the presence of warmer grains in the multi-grain models. Using a single grain model mimicking grain growth and dust settling fails to reproduce the complexity of gas-grain chemistry. Chemical models with a single grain size are sensitive to the adopted grain temperature, and cannot account for all expected effects. A spatial spread of the snowlines is expected to result from the ranges in grain temperature. The amplitude of the effects will depend on the dust disk mass., Comment: 39 pages, 33 figures

Molecular line surveys are among the main tools to probe the structure and physical conditions in protoplanetary disks (PPDs), the birthplace of planets. The large radial and vertical temperature as well as density gradients in these PPDs lead to a complex chemical composition, making chemistry an important step to understand the variety of planetary systems. We aimed to study the chemical content of the protoplanetary disk surrounding GG Tau A, a well-known triple T Tauri system. We used NOEMA with the new correlator PolyFix to observe rotational lines at 2.6 to 4 mm from a few dozen molecules. We analysed the data with a radiative transfer code to derive molecular densities and the abundance relative to $^{13}$CO, which we compare to those of the TMC1 cloud and LkCa15 disk. We report the first detection of CCS in PPDs. We also marginally detect OCS and find 16 other molecules in the GG Tauri outer disk. Ten of them had been found previously, while seven others ($^{13}$CN, N$_2$H$^+$, HNC, DNC, HC$_3$N, CCS, and C$^{34}$S) are new detections in this disk. The analysis confirms that sulphur chemistry is not yet properly understood. The D/H ratio, derived from DCO$^{+}$/HCO$^{+}$, DCN/HCN, and DNC/HNC ratios, points towards a low temperature chemistry. The detection of the rare species CCS confirms that GG Tau is a good laboratory to study the protoplanetary disk chemistry, thanks to its large disk size and mass., Comment: Accepted by Astronomy and Astrophysics

Grain surface chemistry is key to the composition of protoplanetary disks around young stars. The temperature of grains depends on their size. We evaluate the impact of this temperature dependence on the disk chemistry. We model a moderately massive disk with 16 different grain sizes. We use POLARIS to calculate the dust grain temperatures and the local UV flux. We model the chemistry using the 3-phase astrochemical code NAUTILUS. Photoprocesses are handled using frequency-dependent cross-sections, and a new method to account for self and mutual shielding. The multi-grain model outputs are compared to those of single-grain size models (0.1 $\mu$m), with two different assumptions for their equivalent temperature. We find that the Langmuir-Hinshelwood (LH) mechanism at equilibrium temperature is not efficient to form H$_2$ at 3-4 scale heights ($H$), and adopt a parametric fit to a stochastic method to model H$_2$ formation instead. We find the molecular layer composition (1-3 $H$) to depend on the amount of remaining H atoms. Differences in molecular surface densities between single and multi-grain models are mostly due to what occurs above 1.5 $H$. At 100 au, models with colder grains produce H$_2$O and CH$_4$ ices in the midplane, and warmer ones produce more CO$_2$ ices, both allowing efficient depletion of C and O as soon as CO sticks on grain surfaces. Complex organic molecules (COMs) production is enhanced by the presence of warmer grains in the multi-grain models. Using a single grain model mimicking grain growth and dust settling fails to reproduce the complexity of gas-grain chemistry. Chemical models with a single grain size are sensitive to the adopted grain temperature, and cannot account for all expected effects. A spatial spread of the snowlines is expected to result from the ranges in grain temperature. The amplitude of the effects will depend on the dust disk mass., Comment: 39 pages, 33 figures

Context. Planet formation is expected to take place in the first million years of a planetary system through various processes, which remain to be tested through observations. Aims. With the recent discovery, using ALMA, of two gaseous spiral arms inside the 120 au cavity and connected to dusty spirals, the famous protoplanetary disk around AB Aurigae presents a strong incentive for investigating the mechanisms that lead to giant planet formation. A candidate protoplanet located inside a spiral arm has already been claimed in an earlier study based on the same ALMA data. Methods. We used SPHERE at the Very Large Telescope (VLT) to perform near-infrared (IR) high-contrast imaging of AB Aur in polarized and unpolarized light in order to study the morphology of the disk and search for signs of planet formation. Results. SPHERE has delivered the deepest images ever obtained for AB Aur in scattered light. Among the many structures that are yet to be understood, we identified not only the inner spiral arms, but we also resolved a feature in the form of a twist in the eastern spiral at a separation of about 30 au. The twist of the spiral is perfectly reproduced with a planet-driven density wave model when projection effects are accounted for. We measured an azimuthal displacement with respect to the counterpart of this feature in the ALMA data, which is consistent with Keplerian motion on a 4-yr baseline. Another point sxce is detected near the edge of the inner ring, which is likely the result of scattering as opposed to the direct emission from a planet photosphere. We tentatively derived mass constraints for these two features. Conclusions. The twist and its apparent orbital motion could well be the first direct evidence of a connection between a protoplanet candidate and its manifestation as a spiral imprinted in the gas and dust distributions., Comment: accepted in Astronomy and Astrophysics, 8 pages, 5 figures

We used new ALMA $^{13}$CO and C$^{18}$O(3-2) observations obtained at high angular resolution ($\sim$0.2") together with previous CO(3-2) and (6-5) ALMA data and continuum maps at 1.3 and 0.8 mm in order to determine the gas properties (temperature, density, and kinematics) in the cavity and to a lesser extent in the outer disk of GG Tau A, the prototype of a young triple T Tauri star that is surrounded by a massive and extended Keplerian outer disk. By deprojecting, we studied the radial and azimuthal gas distribution and its kinematics. We also applied a new method to improve the deconvolution of the CO data and in particular better quantify the emission from gas inside the cavity. We perform local and nonlocal thermodynamic equilibrium studies in order to determine the excitation conditions and relevant physical parameters inside the ring and in the central cavity. Residual emission after removing a smooth-disk model indicates unresolved structures at our angular resolution, probably in the form of irregular rings or spirals. The outer disk is cold, with a temperature $<20$ K beyond 250 au that drops quickly (r$^{-1}$). The kinematics of the gas inside the cavity reveals infall motions at about 10% of the Keplerian speed. We derive the amount of gas in the cavity, and find that the brightest clumps, which contain about 10% of this mass, have kinetic temperatures 40$-$80 K, CO column densities of a few 10$^{17}$ cm$^{-2}$, and H$_2$ densities around 10$^7$ cm$^{-3}$. Although the gas in the cavity is only a small fraction of the disk mass, the mass accretion rate throughout the cavity is comparable to or higher than the stellar accretion rate. It is accordingly sufficient to sustain the circumstellar disks on a long timescale., Comment: 18 pages, 14 figures

Extreme outbursts in young stars may be a common stage of pre-main-sequence stellar evolution. These outbursts, caused by enhanced accretion and accompanied by increased luminosity, can also strongly impact the evolution of the circumstellar environment. We present ALMA and VLA observations of EX Lupi, a prototypical outburst system, at 100 GHz, 45 GHz, and 15 GHz. We use these data, along with archival ALMA 232 GHz data, to fit radiative transfer models to EX Lupi's circumstellar disk in its quiescent state following the extreme outburst in 2008. The best fit models show a compact disk with a characteristic dust radius of 45 au and a total mass of 0.01 M$_{\odot}$. Our modeling suggests grain growth to sizes of at least 3 mm in the disk, possibly spurred by the recent outburst, and an ice line that has migrated inward to $0.2-0.3$ au post-outburst. At 15 GHz, we detected significant emission over the expected thermal disk emission which we attribute primarily to stellar (gyro)synchrotron and free-free disk emission. Altogether, these results highlight what may be a common impact of outbursts on the circumstellar dust., Comment: Accepted to ApJ, 15 pages, 8 figures

We present a multi-wavelength analysis to reveal the nature of the enigmatic T Tauri triple star system. New optical and infrared measurements are coupled with archival X-ray, UV and mm datasets to show morphologies of disk material and outflow kinematics. A dark lane of obscuring material is seen in silhouette in several emission lines and in model-subtracted ALMA mm continuum dust residuals near the position of T Tau Sa+Sb, revealing the attenuating circumbinary ring around T Tau S. The flux variability of T Tau S is linked in part to the binary orbit; T Tau Sb brightens near orbital apastron as it emerges from behind circumbinary material. Outflow diagnostics confirm that T Tau N powers the blue-shifted western outflow, and the T Tau S binary drives the northwest-southeastern flow. Analysis of the southern outflow shows periodic arcs ejected from the T Tau system. Correlation of these arc locations and tangential kinematics with the orbit timing suggests that launch of the last four southern outflow ejections is contemporaneous with, and perhaps triggered by, the T Tau Sa+Sb binary periastron passage. We present a geometry of the T Tau triple that has the southern components foreground to T Tau N, obscured by a circumbinary ring, with mis-aligned disks and interacting outflows. Particularly, a wind from T Tauri Sa that is perpendicular to its circumstellar disk might interact with the circumbinary material, which may explain conflicting high contrast measurements of the system outflows in the literature. T Tauri is an important laboratory to understand early dynamical processes in young multiple systems. We discuss the historical and future characteristics of the system in this context., Comment: Accepted for publication in the Astrophysical Journal

Context. Planet formation is expected to take place in the first million years of a planetary system through various processes, which remain to be tested through observations. Aims. With the recent discovery, using ALMA, of two gaseous spiral arms inside the 120 au cavity and connected to dusty spirals, the famous protoplanetary disk around AB Aurigae presents a strong incentive for investigating the mechanisms that lead to giant planet formation. A candidate protoplanet located inside a spiral arm has already been claimed in an earlier study based on the same ALMA data. Methods. We used SPHERE at the Very Large Telescope (VLT) to perform near-infrared (IR) high-contrast imaging of AB Aur in polarized and unpolarized light in order to study the morphology of the disk and search for signs of planet formation. Results. SPHERE has delivered the deepest images ever obtained for AB Aur in scattered light. Among the many structures that are yet to be understood, we identified not only the inner spiral arms, but we also resolved a feature in the form of a twist in the eastern spiral at a separation of about 30 au. The twist of the spiral is perfectly reproduced with a planet-driven density wave model when projection effects are accounted for. We measured an azimuthal displacement with respect to the counterpart of this feature in the ALMA data, which is consistent with Keplerian motion on a 4-yr baseline. Another point sxce is detected near the edge of the inner ring, which is likely the result of scattering as opposed to the direct emission from a planet photosphere. We tentatively derived mass constraints for these two features. Conclusions. The twist and its apparent orbital motion could well be the first direct evidence of a connection between a protoplanet candidate and its manifestation as a spiral imprinted in the gas and dust distributions., Comment: accepted in Astronomy and Astrophysics, 8 pages, 5 figures

We used new ALMA $^{13}$CO and C$^{18}$O(3-2) observations obtained at high angular resolution ($\sim$0.2") together with previous CO(3-2) and (6-5) ALMA data and continuum maps at 1.3 and 0.8 mm in order to determine the gas properties (temperature, density, and kinematics) in the cavity and to a lesser extent in the outer disk of GG Tau A, the prototype of a young triple T Tauri star that is surrounded by a massive and extended Keplerian outer disk. By deprojecting, we studied the radial and azimuthal gas distribution and its kinematics. We also applied a new method to improve the deconvolution of the CO data and in particular better quantify the emission from gas inside the cavity. We perform local and nonlocal thermodynamic equilibrium studies in order to determine the excitation conditions and relevant physical parameters inside the ring and in the central cavity. Residual emission after removing a smooth-disk model indicates unresolved structures at our angular resolution, probably in the form of irregular rings or spirals. The outer disk is cold, with a temperature $<20$ K beyond 250 au that drops quickly (r$^{-1}$). The kinematics of the gas inside the cavity reveals infall motions at about 10% of the Keplerian speed. We derive the amount of gas in the cavity, and find that the brightest clumps, which contain about 10% of this mass, have kinetic temperatures 40$-$80 K, CO column densities of a few 10$^{17}$ cm$^{-2}$, and H$_2$ densities around 10$^7$ cm$^{-3}$. Although the gas in the cavity is only a small fraction of the disk mass, the mass accretion rate throughout the cavity is comparable to or higher than the stellar accretion rate. It is accordingly sufficient to sustain the circumstellar disks on a long timescale., Comment: 18 pages, 14 figures

The accuracy of masses of pre-main sequence (PMS) stars derived from their locations on the Hertzsprung-Russell Diagram (HRD) can be tested by comparison with accurate and precise masses determined independently. We present 29 single stars in the Taurus star-forming region (SFR) and 3 in the Ophiuchus SFR with masses measured dynamically to a precision of at least $10 \%$. Our results include 9 updated mass determinations and 3 that have not had their dynamical masses published before. This list of stars with fundamental, dynamical masses, M$_{dyn}$, is drawn from a larger list of 39 targets in the Taurus SFR and 6 in the Ophiuchus SFR. Placing the stars with accurate and precise dynamical masses on HRDs that do not include internal magnetic fields underestimates the mass compared to M$_{dyn}$ by about $30 \%$. Placing them on an HRD that does include magnetic fields yields mass estimates in much better agreement with M$_{dyn}$, with an average difference between M$_{dyn}$ and the estimated track mass of $0.01\pm0.02$~\msun. The ages of the stars, 3--10 MY on tracks that include magnetic fields, is older than the 1--3 MY indicated by the non-magnetic models. The older ages of T Tauri stars predicted by the magnetic models increase the time available for evolution of their disks and formation of the giant gas exoplanets. The agreement between our M$_{dyn}$ values and the masses on the magnetic field tracks provides indirect support for these older ages., Comment: Ap.J. to be published

The accuracy of masses of pre-main sequence (PMS) stars derived from their locations on the Hertzsprung-Russell Diagram (HRD) can be tested by comparison with accurate and precise masses determined independently. We present 29 single stars in the Taurus star-forming region (SFR) and 3 in the Ophiuchus SFR with masses measured dynamically to a precision of at least $10 \%$. Our results include 9 updated mass determinations and 3 that have not had their dynamical masses published before. This list of stars with fundamental, dynamical masses, M$_{dyn}$, is drawn from a larger list of 39 targets in the Taurus SFR and 6 in the Ophiuchus SFR. Placing the stars with accurate and precise dynamical masses on HRDs that do not include internal magnetic fields underestimates the mass compared to M$_{dyn}$ by about $30 \%$. Placing them on an HRD that does include magnetic fields yields mass estimates in much better agreement with M$_{dyn}$, with an average difference between M$_{dyn}$ and the estimated track mass of $0.01\pm0.02$~\msun. The ages of the stars, 3--10 MY on tracks that include magnetic fields, is older than the 1--3 MY indicated by the non-magnetic models. The older ages of T Tauri stars predicted by the magnetic models increase the time available for evolution of their disks and formation of the giant gas exoplanets. The agreement between our M$_{dyn}$ values and the masses on the magnetic field tracks provides indirect support for these older ages., Comment: Ap.J. to be published

Studying molecular species in protoplanetary disks is very useful to characterize the properties of these objects, which are the site of planet formation. We attempt to constrain the chemistry of S-bearing molecules in the cold parts of circumstellar disk of GG Tau A. We searched for H$_2$S, CS, SO, and SO$_2$ in the dense disk around GG Tau A with the NOrthem Extended Millimeter Array (NOEMA) interferometer. We detected H$_2$S emission from the dense and cold ring orbiting around GG Tau A. This is the first detection of H$_2$S in a protoplanetary disk. We also detected HCO$^+$, H$^{13}$CO$^+$, and DCO$^+$ in the disk. Upper limits for other molecules, CCS, SO$_2$, SO, HC$_3$N, and $c$-C$_3$H$_2$ are also obtained. The observed DCO$^+$/HCO$^+$ ratio is similar to those in other disks. The observed column densities, derived using our radiative transfer code DiskFit, are then compared with those from our chemical code Nautilus. The column densities are in reasonable agreement for DCO$^{+}$, CS, CCS, and SO$_2$. For H$_2$S and SO, our predicted vertical integrated column densities are more than a factor of 10 higher than the measured values. Our results reinforce the hypothesis that only a strong sulfur depletion may explain the low observed H$_2$S column density in the disk. The H$_2$S detection in GG Tau A is most likely linked to the much larger mass of this disk compared to that in other T Tauri systems., Comment: 6 pages, 7 figures, 3 tables

Context. Several sulfur-bearing molecules are observed in the interstellar medium and in comets, in strong contrast to protoplanetary disks where only CS, H$_2$CS and SO have been detected so far. Aims. We combine observations and chemical models to constrain the sulfur abundances and their sensitivity to physical and chemical conditions in the DM Tau protoplanetary disk. Methods. We obtained $0.5^{"}$ ALMA observations of DM Tau in Bands 4 and 6 in lines of CS, SO, SO$_2$, OCS, CCS, H$_2$CS and H$_2$S, achieving a $\sim 5$ mJy sensitivity. Using the non-LTE radiative transfer code RADEX and the forward-modeling tool DiskFit, disk-averaged CS column densities and upper limits for the other species were derived. Results. Only CS was detected with a derived column density of $\sim 2-6 \times 10^{12}$ cm$^{-2}$. We report a first tentative detection of SO$_2$ in DM Tau. The upper limits range between $\sim 10^{11}$ and $10^{14}$ cm$^{-2}$ for the other S-bearing species. The best-fit chemical model matching these values requires a gas-phase C/O ratio of > 1 at $r \sim 50-100$ au. With chemical modeling we demonstrate that sulfur-bearing species could be robust tracers of the gas-phase C/O ratio, surface reaction rates, grain size and UV intensities. Conclusions. The lack of detections of a variety of sulfur-bearing molecules in DM Tau other than CS implies a dearth of reactive sulfur in the gas phase, either through efficient freeze-out or because most of the elemental sulfur is in other large species, as found in comets. The inferred high CS/SO and CS/SO$_2$ ratios require a non-solar C/O gas-phase ratio of > 1, consistent with the recent observations of hydrocarbon rings in DM Tau. The stronger depletion of oxygen-bearing S-species compared to CS is likely linked to the low observed abundances of gaseous water in DM Tau and points to a removal mechanism of oxygen from the gas., Comment: Accepted by Astronomy & Astrophysics (15 pages, 15 figures, 8 tables)

Observations by the Atacama Large Millimeter/sub-millimeter Array of the dust continuum and $^{13}$CO(3-2) millimeter emissions of the triple stellar system GG Tau A are analysed, giving evidence for a rotating gas disc and a concentric and coplanar dust ring. The present work complements an earlier analysis (Tang et al. 2016) by exploring detailed properties of the gas disc. A 95% confidence level upper limit of 0.24 arcsec (34 au) is placed on the disc scale height at a distance of 1 arcsec (140 au) from the central stars. Evidence for Keplerian rotation of the gas disc is presented, the rotation velocity reaching ~3.1 kms$^{-1}$ at 1 arcsec from the central stars, and a 99% confidence level upper limit of 9% is placed on a possible in-fall velocity relative contribution. Variations of the intensity across the disc area are studied in detail and confirm the presence of a hot spot in the south-eastern quadrant. However several other significant intensity variations, in particular a depression in the northern direction, are also revealed. Variations of the intensity are found to be positively correlated to variations of the line width. Possible contributions to the measured line width are reviewed, suggesting an increase of the disc temperature and opacity with decreasing distance from the stars., Comment: 15 pages, 15 figures, 2 tables and 7 references, accepted for publication in RAA

Studying molecular species in protoplanetary disks is very useful to characterize the properties of these objects, which are the site of planet formation. We attempt to constrain the chemistry of S-bearing molecules in the cold parts of circumstellar disk of GG Tau A. We searched for H$_2$S, CS, SO, and SO$_2$ in the dense disk around GG Tau A with the NOrthem Extended Millimeter Array (NOEMA) interferometer. We detected H$_2$S emission from the dense and cold ring orbiting around GG Tau A. This is the first detection of H$_2$S in a protoplanetary disk. We also detected HCO$^+$, H$^{13}$CO$^+$, and DCO$^+$ in the disk. Upper limits for other molecules, CCS, SO$_2$, SO, HC$_3$N, and $c$-C$_3$H$_2$ are also obtained. The observed DCO$^+$/HCO$^+$ ratio is similar to those in other disks. The observed column densities, derived using our radiative transfer code DiskFit, are then compared with those from our chemical code Nautilus. The column densities are in reasonable agreement for DCO$^{+}$, CS, CCS, and SO$_2$. For H$_2$S and SO, our predicted vertical integrated column densities are more than a factor of 10 higher than the measured values. Our results reinforce the hypothesis that only a strong sulfur depletion may explain the low observed H$_2$S column density in the disk. The H$_2$S detection in GG Tau A is most likely linked to the much larger mass of this disk compared to that in other T Tauri systems., Comment: 6 pages, 7 figures, 3 tables

Context. Several sulfur-bearing molecules are observed in the interstellar medium and in comets, in strong contrast to protoplanetary disks where only CS, H$_2$CS and SO have been detected so far. Aims. We combine observations and chemical models to constrain the sulfur abundances and their sensitivity to physical and chemical conditions in the DM Tau protoplanetary disk. Methods. We obtained $0.5^{"}$ ALMA observations of DM Tau in Bands 4 and 6 in lines of CS, SO, SO$_2$, OCS, CCS, H$_2$CS and H$_2$S, achieving a $\sim 5$ mJy sensitivity. Using the non-LTE radiative transfer code RADEX and the forward-modeling tool DiskFit, disk-averaged CS column densities and upper limits for the other species were derived. Results. Only CS was detected with a derived column density of $\sim 2-6 \times 10^{12}$ cm$^{-2}$. We report a first tentative detection of SO$_2$ in DM Tau. The upper limits range between $\sim 10^{11}$ and $10^{14}$ cm$^{-2}$ for the other S-bearing species. The best-fit chemical model matching these values requires a gas-phase C/O ratio of > 1 at $r \sim 50-100$ au. With chemical modeling we demonstrate that sulfur-bearing species could be robust tracers of the gas-phase C/O ratio, surface reaction rates, grain size and UV intensities. Conclusions. The lack of detections of a variety of sulfur-bearing molecules in DM Tau other than CS implies a dearth of reactive sulfur in the gas phase, either through efficient freeze-out or because most of the elemental sulfur is in other large species, as found in comets. The inferred high CS/SO and CS/SO$_2$ ratios require a non-solar C/O gas-phase ratio of > 1, consistent with the recent observations of hydrocarbon rings in DM Tau. The stronger depletion of oxygen-bearing S-species compared to CS is likely linked to the low observed abundances of gaseous water in DM Tau and points to a removal mechanism of oxygen from the gas., Comment: Accepted by Astronomy & Astrophysics (15 pages, 15 figures, 8 tables)

Observations by the Atacama Large Millimeter/sub-millimeter Array of the dust continuum and $^{13}$CO(3-2) millimeter emissions of the triple stellar system GG Tau A are analysed, giving evidence for a rotating gas disc and a concentric and coplanar dust ring. The present work complements an earlier analysis (Tang et al. 2016) by exploring detailed properties of the gas disc. A 95% confidence level upper limit of 0.24 arcsec (34 au) is placed on the disc scale height at a distance of 1 arcsec (140 au) from the central stars. Evidence for Keplerian rotation of the gas disc is presented, the rotation velocity reaching ~3.1 kms$^{-1}$ at 1 arcsec from the central stars, and a 99% confidence level upper limit of 9% is placed on a possible in-fall velocity relative contribution. Variations of the intensity across the disc area are studied in detail and confirm the presence of a hot spot in the south-eastern quadrant. However several other significant intensity variations, in particular a depression in the northern direction, are also revealed. Variations of the intensity are found to be positively correlated to variations of the line width. Possible contributions to the measured line width are reviewed, suggesting an increase of the disc temperature and opacity with decreasing distance from the stars., Comment: 15 pages, 15 figures, 2 tables and 7 references, accepted for publication in RAA

We report new dynamical masses for 5 pre-main sequence (PMS) stars in the L1495 region of the Taurus star-forming region (SFR) and 6 in the L1688 region of the Ophiuchus SFR. Since these regions have VLBA parallaxes these are absolute measurements of the stars' masses and are independent of their effective temperatures and luminosities. Seven of the stars have masses $<0.6$ solar masses, thus providing data in a mass range with little data, and of these, 6 are measured to precision $< 5 \%$. We find 8 stars with masses in the range 0.09 to 1.1 solar mass that agree well with the current generation of PMS evolutionary models. The ages of the stars we measured in the Taurus SFR are in the range 1-3 MY, and $<1$ MY for those in L1688. We also measured the dynamical masses of 14 stars in the ALMA archival data for Akeson~\&~Jensen's Cycle 0 project on binaries in the Taurus SFR. We find that the masses of 7 of the targets are so large that they cannot be reconciled with reported values of their luminosity and effective temperature. We suggest that these targets are themselves binaries or triples., Comment: 20 pages

Determining the gas density and temperature structures of protoplanetary disks is a fundamental task to constrain planet formation theories. This is a challenging procedure and most determinations are based on model-dependent assumptions. We attempt a direct determination of the radial and vertical temperature structure of the Flying Saucer disk, thanks to its favorable inclination of 90 degrees. We present a method based on the tomographic study of an edge-on disk. Using ALMA, we observe at 0.5$"$ resolution the Flying Saucer in CO J=2-1 and CS J=5-4. This edge-on disk appears in silhouette against the CO J=2-1 emission from background molecular clouds in $\rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity provide a direct measure of the gas temperature as a function of radius and height above the disk mid-plane. The overall thermal structure is consistent with model predictions, with a cold ($< 15-12 $~K), CO-depleted mid-plane, and a warmer disk atmosphere. However, we find evidence for CO gas along the mid-plane beyond a radius of about 200\,au, coincident with a change of grain properties. Such a behavior is expected in case of efficient rise of UV penetration re-heating the disk and thus allowing CO thermal desorption or favoring direct CO photo-desorption. CO is also detected up to 3-4 scale heights while CS is confined around 1 scale height above the mid-plane. The limits of the method due to finite spatial and spectral resolutions are also discussed. This method appears to be very promising to determine the gas structure of planet-forming disks, provided that the molecular data have an angular resolution which is high enough, of the order of $0.3 - 0.1"$ at the distance of the nearest star forming regions., Comment: 14 pages + 11 figures

We present ALMA Cycle~2 observations at $0.5^{\prime\prime}$ resolution of TW Hya of CS $J=5-4$ emission. The radial profile of the integrated line emission displays oscillatory features outwards of $1.5^{\prime\prime}$ ($\approx 90$ au). A dip-like feature at $1.6^{\prime\prime}$ is coincident in location, depth and width with features observed in dust scattered light at near-infrared wavelengths. Using a thermochemical model indicative of TW Hya, gas-grain chemical modelling and non-LTE radiative transfer, we demonstrate that such a feature can be reproduced with a surface density depression, consistent with the modelling performed for scattered light observations of TW Hya. We further demonstrate that a gap in the dust distribution and dust opacity only cannot reproduce the observed CS feature. The outer enhancement at $3.1^{\prime\prime}$ is identified as a region of intensified desorption due to enhanced penetration of the interstellar FUV radiation at the exponential edge of the disk surface density, which intensifies the photochemical processing of gas and ices., Comment: 12 pages, 16 figures, published in ApJ

Accurate temperature calculations for circumstellar disks are particularly important for their chemical evolution. Their temperature distribution is determined by the optical properties of the dust grains, which, among other parameters, depend on their radius. However, in most disk studies, only average optical properties and thus an average temperature is assumed to account for an ensemble of grains with different radii. We investigate the impact of subdividing the grain radius distribution into multiple sub-intervals on the resulting dust temperature distribution and spectral energy distribution (SED). These quantities were computed for two different scenarios: (1) Radius distribution represented by 16 logarithmically distributed radius intervals, and (2) radius distribution represented by a single grain species with averaged optical properties (reference). Within the considered parameter range, i.e., of grain radii between 5 nm and 1 mm and an optically thin and thick disk with a parameterized density distribution, we obtain the following results: In optically thin disk regions, the temperature spread can be as large as ~63% and the relative grain surface below a certain temperature is lower than in the reference disk. With increasing optical depth, the difference in the midplane temperature and the relative grain surface below a certain temperature decreases. Furthermore, below ~20K, this fraction is higher for the reference disk than for the case of multiple grain radii, while it shows the opposite behavior for temperatures above this threshold. The thermal emission in the case of multiple grain radii at short wavelengths is stronger than for the reference disk. The freeze-out radius is a function of grain radius, spanning a radial range between the coldest and warmest grain species of ~30AU., Comment: 14 pages, 22 figures, 2 tables

Determining the gas density and temperature structures of protoplanetary disks is a fundamental task to constrain planet formation theories. This is a challenging procedure and most determinations are based on model-dependent assumptions. We attempt a direct determination of the radial and vertical temperature structure of the Flying Saucer disk, thanks to its favorable inclination of 90 degrees. We present a method based on the tomographic study of an edge-on disk. Using ALMA, we observe at 0.5$"$ resolution the Flying Saucer in CO J=2-1 and CS J=5-4. This edge-on disk appears in silhouette against the CO J=2-1 emission from background molecular clouds in $\rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity provide a direct measure of the gas temperature as a function of radius and height above the disk mid-plane. The overall thermal structure is consistent with model predictions, with a cold ($< 15-12 $~K), CO-depleted mid-plane, and a warmer disk atmosphere. However, we find evidence for CO gas along the mid-plane beyond a radius of about 200\,au, coincident with a change of grain properties. Such a behavior is expected in case of efficient rise of UV penetration re-heating the disk and thus allowing CO thermal desorption or favoring direct CO photo-desorption. CO is also detected up to 3-4 scale heights while CS is confined around 1 scale height above the mid-plane. The limits of the method due to finite spatial and spectral resolutions are also discussed. This method appears to be very promising to determine the gas structure of planet-forming disks, provided that the molecular data have an angular resolution which is high enough, of the order of $0.3 - 0.1"$ at the distance of the nearest star forming regions., Comment: 14 pages + 11 figures

We report new dynamical masses for 5 pre-main sequence (PMS) stars in the L1495 region of the Taurus star-forming region (SFR) and 6 in the L1688 region of the Ophiuchus SFR. Since these regions have VLBA parallaxes these are absolute measurements of the stars' masses and are independent of their effective temperatures and luminosities. Seven of the stars have masses $<0.6$ solar masses, thus providing data in a mass range with little data, and of these, 6 are measured to precision $< 5 \%$. We find 8 stars with masses in the range 0.09 to 1.1 solar mass that agree well with the current generation of PMS evolutionary models. The ages of the stars we measured in the Taurus SFR are in the range 1-3 MY, and $<1$ MY for those in L1688. We also measured the dynamical masses of 14 stars in the ALMA archival data for Akeson~\&~Jensen's Cycle 0 project on binaries in the Taurus SFR. We find that the masses of 7 of the targets are so large that they cannot be reconciled with reported values of their luminosity and effective temperature. We suggest that these targets are themselves binaries or triples., Comment: 20 pages

Accurate temperature calculations for circumstellar disks are particularly important for their chemical evolution. Their temperature distribution is determined by the optical properties of the dust grains, which, among other parameters, depend on their radius. However, in most disk studies, only average optical properties and thus an average temperature is assumed to account for an ensemble of grains with different radii. We investigate the impact of subdividing the grain radius distribution into multiple sub-intervals on the resulting dust temperature distribution and spectral energy distribution (SED). These quantities were computed for two different scenarios: (1) Radius distribution represented by 16 logarithmically distributed radius intervals, and (2) radius distribution represented by a single grain species with averaged optical properties (reference). Within the considered parameter range, i.e., of grain radii between 5 nm and 1 mm and an optically thin and thick disk with a parameterized density distribution, we obtain the following results: In optically thin disk regions, the temperature spread can be as large as ~63% and the relative grain surface below a certain temperature is lower than in the reference disk. With increasing optical depth, the difference in the midplane temperature and the relative grain surface below a certain temperature decreases. Furthermore, below ~20K, this fraction is higher for the reference disk than for the case of multiple grain radii, while it shows the opposite behavior for temperatures above this threshold. The thermal emission in the case of multiple grain radii at short wavelengths is stronger than for the reference disk. The freeze-out radius is a function of grain radius, spanning a radial range between the coldest and warmest grain species of ~30AU., Comment: 14 pages, 22 figures, 2 tables

We present ALMA Cycle~2 observations at $0.5^{\prime\prime}$ resolution of TW Hya of CS $J=5-4$ emission. The radial profile of the integrated line emission displays oscillatory features outwards of $1.5^{\prime\prime}$ ($\approx 90$ au). A dip-like feature at $1.6^{\prime\prime}$ is coincident in location, depth and width with features observed in dust scattered light at near-infrared wavelengths. Using a thermochemical model indicative of TW Hya, gas-grain chemical modelling and non-LTE radiative transfer, we demonstrate that such a feature can be reproduced with a surface density depression, consistent with the modelling performed for scattered light observations of TW Hya. We further demonstrate that a gap in the dust distribution and dust opacity only cannot reproduce the observed CS feature. The outer enhancement at $3.1^{\prime\prime}$ is identified as a region of intensified desorption due to enhanced penetration of the interstellar FUV radiation at the exponential edge of the disk surface density, which intensifies the photochemical processing of gas and ices., Comment: 12 pages, 16 figures, published in ApJ

The increased sensitivity of millimeter-wave facilities now makes possible the detection of low amounts of gas in debris disks. Some of the gas-rich debris disks harbor peculiar properties, with possible pristine gas and secondary generated dust. The origin of the gas in these hybrid disks is strongly debated and the current sample is too sparse to understand this phenomenon. More detections are necessary to increase the statistics on this population. Lying at the final stages of evolution of proto-planetary disks and at the beginning of the debris disk phase, these objects could provide new insight into the processes involved in the making of planetary systems. We carried out a deep survey of the 12CO(2-1) and 12CO(3-2) lines with the APEX and IRAM radiotelescopes in young debris disks selected according to hybrid disk properties. The survey is complemented with a bibliographic study of the ratio between the emission of the gas and the continuum (S_CO/F_cont) in CTTS, Herbig Ae, WTTS, hybrid, and debris disks. Our sub-mm survey comprises 25 stars, including 17 new targets, and we increase the sensitivity limit by a factor 2 on eight sources compared to similar published studies. We report a 4sigma tentative detection of a double-peaked 12CO(2-1) line around HD23642; an eclipsing binary located in the Pleiades. We also reveal a correlation between the emission of the CO gas and the dust continuum from CTTS, Herbig Ae and few debris disks. The observed trend of the gas to dust flux ratio suggests a concurrent dissipation of the dust and gas components. Hybrid disks systematically lie above this trend, suggesting that these systems may witness a transient phase, when the dust has evolved more rapidly than the gas, with a flux ratio S_CO/F_cont enhanced by a factor of between 10 and 100 compared to standard (proto-)planetary disks., Comment: Accepted for publication in Astronomy&Astrophysics

Protoplanetary disks are the target of many chemical studies (both observational and theoretical) as they contain the building material for planets. Their large vertical and radial gradients in density and temperature make them challenging objects for chemical models. In the outer part of these disks, the large densities and low temperatures provide a particular environment where the binding of species onto the dust grains can be very efficient and can affect the gas-phase chemical composition. We attempt to quantify to what extent the vertical abundance profiles and the integrated column densities of molecules predicted by a detailed gas-grain code are affected by the treatment of the molecular hydrogen physisorption at the surface of the grains. We performed three different models using the Nautilus gas-grain code. One model uses a H2 binding energy on the surface of water (440 K) and produces strong sticking of H2. Another model uses a small binding energy of 23 K (as if there were already a monolayer of H2), and the sticking of H$_2$ is almost negligible. Finally, the remaining model is an intermediate solution known as the encounter desorption mechanism. We show that the efficiency of molecular hydrogen binding (and thus its abundance at the surface of the grains) can have a quantitative effect on the predicted column densities in the gas phase of major species such as CO, CS, CN, and HCN., Comment: Accepted for publication as a Research Note in A&A

Deuterated species are unique and powerful tools in astronomy since they can probe the physical conditions, chemistry, and ionization level of various astrophysical media. Recent observations of several deuterated species along with some of their spin isomeric forms have rekindled the interest for more accurate studies on deuterium fractionation. This paper presents the first publicly available chemical network of multiply deuterated species along with spin chemistry implemented on the latest state-of-the-art gas-grain chemical code `NAUTILUS'. D/H ratios for all deuterated species observed at different positions of TMC-1 are compared with the results of our model, which considers multiply deuterated species along with the spin chemistry of light hydrogen bearing species H2, H2+, H3+ and their isotopologues. We also show the differences in the modeled abundances of non-deuterated species after the inclusion of deuteration and spin chemistry in the model. Finally, we present a list of potentially observable deuterated species in TMC-1 awaiting detection., Comment: Accepted for publication in MNRAS on Dec. 22nd

(abridged) We used the IRAM 30-m to perform a sensitive wideband survey of 30 protoplanetary disks in the Taurus Auriga region. We simultaneously observed HCO$^+$(3-2), HCN(3-2), C$_2$H(3-2), CS(5-4), and two transitions of SO. We combine the results with a previous survey which observed $^{13}$CO (2-1), CN(2-1), two o-H$_2$CO lines and one of SO. We use available interferometric data to derive excitation temperatures of CN and C$_2$H in several sources. We determine characteristic sizes of the gas disks and column densities of all molecules using a parametric power-law disk model. Our study is mostly sensitive to molecules at 200-400 au from the stars. We compare the derived column densities to the predictions of an extensive gas-grain chemical disk model, under conditions representative of T Tauri disks. This survey provides 20 new detections of HCO$^+$ in disks, 18 in HCN, 11 in C$_2$H, 8 in CS and 4 in SO. HCO$^+$ is detected in almost all sources, and its J=3-2 line is essentially optically thick, providing good estimates of the disk radii. The other transitions are (at least partially) optically thin. Variations of the column density ratios do not correlate with any specific property of the star or disk. Disks around Herbig Ae stars appear less rich in molecules than those around T Tauri stars, although the sample remains small. SO is only found in the (presumably younger) embedded objects, perhaps reflecting an evolution of the S chemistry due to increasing depletion with time. Overall, the molecular column densities, and in particular the CN/HCN and CN/C$_2$H ratios, are well reproduced by gas-grain chemistry in cold disks. This study provides a census of simple molecules in disks of radii $> 200-300$ au. Extending that to smaller disks, or searching for less abundant or more complex molecules requires a much more sensitive facility, i.e. NOEMA and ALMA., Comment: To appear in Astronomy & Astrophysics

The first hydrostatic core, also called the first Larson core, is one of the first steps in low-mass star formation, as predicted by theory. With recent and future high performance telescopes, details of these first phases become accessible, and observations may confirm theory and even bring new challenges for theoreticians. In this context, we study from a theoretical point of view the chemical and physical evolution of the collapse of prestellar cores until the formation of the first Larson core, in order to better characterize this early phase in the star formation process. We couple a state-of-the-art hydrodynamical model with full gas-grain chemistry, using different assumptions on the magnetic field strength and orientation. We extract the different components of each collapsing core (i.e., the central core, the outflow, the disk, the pseudodisk, and the envelope) to highlight their specific physical and chemical characteristics. Each component often presents a specific physical history, as well as a specific chemical evolution. From some species, the components can clearly be differentiated. The different core models can also be chemically differentiated. Our simulation suggests some chemical species as tracers of the different components of a collapsing prestellar dense core, and as tracers of the magnetic field characteristics of the core. From this result, we pinpoint promising key chemical species to be observed., Comment: Accepted for publication in ApJ

Dust determines the temperature structure of protoplanetary disks. However, dust temperature determinations almost invariably rely on a complex modeling of the Spectral Energy Distribution. We attempt a direct determination of the temperature of large grains emitting at mm wavelengths.} We observe the edge-on dust disk of the Flying Saucer, which appears in silhouette against the CO J=2-1 emission from a background molecular cloud in $\rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity allows us to directly measure the %absorbing dust temperature. The dust opacity can then be derived from the emitted continuum radiation. The dust disk absorbs the radiation from the CO clouds at several velocities. We derive very low dust temperatures, 5 to 7 K at radii around 100 au, which is much lower than most model predictions. The dust optical depth is $> 0.2$ at 230 GHz, and the scale height at 100 au is at least 8 au (best fit 13 au). However, the dust disk is very flat (flaring index -0.35), which is indicative of dust settling in the outer parts., Comment: 5 pages

Deuterated species are unique and powerful tools in astronomy since they can probe the physical conditions, chemistry, and ionization level of various astrophysical media. Recent observations of several deuterated species along with some of their spin isomeric forms have rekindled the interest for more accurate studies on deuterium fractionation. This paper presents the first publicly available chemical network of multiply deuterated species along with spin chemistry implemented on the latest state-of-the-art gas-grain chemical code `NAUTILUS'. D/H ratios for all deuterated species observed at different positions of TMC-1 are compared with the results of our model, which considers multiply deuterated species along with the spin chemistry of light hydrogen bearing species H2, H2+, H3+ and their isotopologues. We also show the differences in the modeled abundances of non-deuterated species after the inclusion of deuteration and spin chemistry in the model. Finally, we present a list of potentially observable deuterated species in TMC-1 awaiting detection., Comment: Accepted for publication in MNRAS on Dec. 22nd

The increased sensitivity of millimeter-wave facilities now makes possible the detection of low amounts of gas in debris disks. Some of the gas-rich debris disks harbor peculiar properties, with possible pristine gas and secondary generated dust. The origin of the gas in these hybrid disks is strongly debated and the current sample is too sparse to understand this phenomenon. More detections are necessary to increase the statistics on this population. Lying at the final stages of evolution of proto-planetary disks and at the beginning of the debris disk phase, these objects could provide new insight into the processes involved in the making of planetary systems. We carried out a deep survey of the 12CO(2-1) and 12CO(3-2) lines with the APEX and IRAM radiotelescopes in young debris disks selected according to hybrid disk properties. The survey is complemented with a bibliographic study of the ratio between the emission of the gas and the continuum (S_CO/F_cont) in CTTS, Herbig Ae, WTTS, hybrid, and debris disks. Our sub-mm survey comprises 25 stars, including 17 new targets, and we increase the sensitivity limit by a factor 2 on eight sources compared to similar published studies. We report a 4sigma tentative detection of a double-peaked 12CO(2-1) line around HD23642; an eclipsing binary located in the Pleiades. We also reveal a correlation between the emission of the CO gas and the dust continuum from CTTS, Herbig Ae and few debris disks. The observed trend of the gas to dust flux ratio suggests a concurrent dissipation of the dust and gas components. Hybrid disks systematically lie above this trend, suggesting that these systems may witness a transient phase, when the dust has evolved more rapidly than the gas, with a flux ratio S_CO/F_cont enhanced by a factor of between 10 and 100 compared to standard (proto-)planetary disks., Comment: Accepted for publication in Astronomy&Astrophysics

Protoplanetary disks are the target of many chemical studies (both observational and theoretical) as they contain the building material for planets. Their large vertical and radial gradients in density and temperature make them challenging objects for chemical models. In the outer part of these disks, the large densities and low temperatures provide a particular environment where the binding of species onto the dust grains can be very efficient and can affect the gas-phase chemical composition. We attempt to quantify to what extent the vertical abundance profiles and the integrated column densities of molecules predicted by a detailed gas-grain code are affected by the treatment of the molecular hydrogen physisorption at the surface of the grains. We performed three different models using the Nautilus gas-grain code. One model uses a H2 binding energy on the surface of water (440 K) and produces strong sticking of H2. Another model uses a small binding energy of 23 K (as if there were already a monolayer of H2), and the sticking of H$_2$ is almost negligible. Finally, the remaining model is an intermediate solution known as the encounter desorption mechanism. We show that the efficiency of molecular hydrogen binding (and thus its abundance at the surface of the grains) can have a quantitative effect on the predicted column densities in the gas phase of major species such as CO, CS, CN, and HCN., Comment: Accepted for publication as a Research Note in A&A

(abridged) We used the IRAM 30-m to perform a sensitive wideband survey of 30 protoplanetary disks in the Taurus Auriga region. We simultaneously observed HCO$^+$(3-2), HCN(3-2), C$_2$H(3-2), CS(5-4), and two transitions of SO. We combine the results with a previous survey which observed $^{13}$CO (2-1), CN(2-1), two o-H$_2$CO lines and one of SO. We use available interferometric data to derive excitation temperatures of CN and C$_2$H in several sources. We determine characteristic sizes of the gas disks and column densities of all molecules using a parametric power-law disk model. Our study is mostly sensitive to molecules at 200-400 au from the stars. We compare the derived column densities to the predictions of an extensive gas-grain chemical disk model, under conditions representative of T Tauri disks. This survey provides 20 new detections of HCO$^+$ in disks, 18 in HCN, 11 in C$_2$H, 8 in CS and 4 in SO. HCO$^+$ is detected in almost all sources, and its J=3-2 line is essentially optically thick, providing good estimates of the disk radii. The other transitions are (at least partially) optically thin. Variations of the column density ratios do not correlate with any specific property of the star or disk. Disks around Herbig Ae stars appear less rich in molecules than those around T Tauri stars, although the sample remains small. SO is only found in the (presumably younger) embedded objects, perhaps reflecting an evolution of the S chemistry due to increasing depletion with time. Overall, the molecular column densities, and in particular the CN/HCN and CN/C$_2$H ratios, are well reproduced by gas-grain chemistry in cold disks. This study provides a census of simple molecules in disks of radii $> 200-300$ au. Extending that to smaller disks, or searching for less abundant or more complex molecules requires a much more sensitive facility, i.e. NOEMA and ALMA., Comment: To appear in Astronomy & Astrophysics

The first hydrostatic core, also called the first Larson core, is one of the first steps in low-mass star formation, as predicted by theory. With recent and future high performance telescopes, details of these first phases become accessible, and observations may confirm theory and even bring new challenges for theoreticians. In this context, we study from a theoretical point of view the chemical and physical evolution of the collapse of prestellar cores until the formation of the first Larson core, in order to better characterize this early phase in the star formation process. We couple a state-of-the-art hydrodynamical model with full gas-grain chemistry, using different assumptions on the magnetic field strength and orientation. We extract the different components of each collapsing core (i.e., the central core, the outflow, the disk, the pseudodisk, and the envelope) to highlight their specific physical and chemical characteristics. Each component often presents a specific physical history, as well as a specific chemical evolution. From some species, the components can clearly be differentiated. The different core models can also be chemically differentiated. Our simulation suggests some chemical species as tracers of the different components of a collapsing prestellar dense core, and as tracers of the magnetic field characteristics of the core. From this result, we pinpoint promising key chemical species to be observed., Comment: Accepted for publication in ApJ

Dust determines the temperature structure of protoplanetary disks. However, dust temperature determinations almost invariably rely on a complex modeling of the Spectral Energy Distribution. We attempt a direct determination of the temperature of large grains emitting at mm wavelengths.} We observe the edge-on dust disk of the Flying Saucer, which appears in silhouette against the CO J=2-1 emission from a background molecular cloud in $\rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity allows us to directly measure the %absorbing dust temperature. The dust opacity can then be derived from the emitted continuum radiation. The dust disk absorbs the radiation from the CO clouds at several velocities. We derive very low dust temperatures, 5 to 7 K at radii around 100 au, which is much lower than most model predictions. The dust optical depth is $> 0.2$ at 230 GHz, and the scale height at 100 au is at least 8 au (best fit 13 au). However, the dust disk is very flat (flaring index -0.35), which is indicative of dust settling in the outer parts., Comment: 5 pages

Debris disks are usually thought to be gas-poor, the gas being dissipated by accretion or evaporation during the protoplanetary phase. HD141569A is a 5 Myr old star harboring a famous debris disk, with multiple rings and spiral features. We present here the first PdBI maps of the 12CO(2-1), 13CO(2-1) gas and dust emission at 1.3 mm in this disk. The analysis reveals there is still a large amount of (primordial) gas extending out to 250 au, i. e. inside the rings observed in scattered light. HD141569A is thus a hybrid disk with a huge debris component, where dust has evolved and is produced by collisions, with a large remnant reservoir of gas., Comment: To appear in "Young Stars and Planets Near the Sun", Proceedings of IAU Symposium No. 314 (Cambridge University Press), J.H. Kastner, B. Stelzer, S.A. Metchev, eds

The gas mass of protoplanetary disks, and the gas-to-dust ratio, are two key elements driving the evolution of these disks and the formation of planetary system. We explore here to what extent CO (or its isotopologues) can be used as a tracer of gas mass. We use a detailed gas-grain chemical model and study the evolution of the disk composition, starting from a dense pre-stellar core composition. We explore a range of disk temperature profiles, cosmic rays ionization rates, and disk ages for a disk model representative of T Tauri stars. At the high densities that prevail in disks, we find that, due to fast reactions on grain surfaces, CO can be converted to less volatile forms (principally s-CO$_2$, and to a lesser extent s-CH$_4$) instead of being evaporated over a wide range of temperature. The canonical gas-phase abundance of 10$^{-4}$ is only reached above about 30-35 K. The dominant Carbon bearing entity depends on the temperature structure and age of the disk. The chemical evolution of CO is also sensitive to the cosmic rays ionization rate. Larger gas phase CO abundances are found in younger disks. Initial conditions, such as parent cloud age and density, have a limited impact. This study reveals that CO gas-phase abundance is heavily dependent on grain surface processes, which remain very incompletely understood so far. The strong dependence on dust temperature profile makes CO a poor tracer of the gas-phase content of disks., Comment: 11 pages, 9 figures, Accepted for publication in A&A

We attempt to determine the molecular composition of disks around young low-mass stars in the $\rho$ Oph region and to compare our results with a similar study performed in the Taurus-Auriga region. We used the IRAM 30 m telescope to perform a sensitive search for CN N=2-1 in 29 T Tauri stars located in the $\rho$ Oph and upper Scorpius regions. $^{13}$CO J=2-1 is observed simultaneously to provide an indication of the level of confusion with the surrounding molecular cloud. The bandpass also contains two transitions of ortho-H$_2$CO, one of SO, and the C$^{17}$O J=2-1 line, which provides complementary information on the nature of the emission. Contamination by molecular cloud in $^{13}$CO and even C$^{17}$O is ubiquitous. The CN detection rate appears to be lower than for the Taurus region, with only four sources being detected (three are attributable to disks). H$_2$CO emission is found more frequently, but appears in general to be due to the surrounding cloud. The weaker emission than in Taurus may suggest that the average disk size in the $\rho$ Oph region is smaller than in the Taurus cloud. Chemical modeling shows that the somewhat higher expected disk temperatures in $\rho$ Oph play a direct role in decreasing the CN abundance. Warmer dust temperatures contribute to convert CN into less volatile forms. In such a young region, CN is no longer a simple, sensitive tracer of disks, and observations with other tracers and at high enough resolution with ALMA are required to probe the gas disk population., Comment: 18 pages, 5 figures, accepted for publication in A&A

Debris disks are usually thought to be gas-poor, the gas being dissipated by accretion or evaporation during the protoplanetary phase. HD141569A is a 5 Myr old star harboring a famous debris disk, with multiple rings and spiral features. We present here the first PdBI maps of the 12CO(2-1), 13CO(2-1) gas and dust emission at 1.3 mm in this disk. The analysis reveals there is still a large amount of (primordial) gas extending out to 250 au, i. e. inside the rings observed in scattered light. HD141569A is thus a hybrid disk with a huge debris component, where dust has evolved and is produced by collisions, with a large remnant reservoir of gas., Comment: To appear in "Young Stars and Planets Near the Sun", Proceedings of IAU Symposium No. 314 (Cambridge University Press), J.H. Kastner, B. Stelzer, S.A. Metchev, eds

The gas mass of protoplanetary disks, and the gas-to-dust ratio, are two key elements driving the evolution of these disks and the formation of planetary system. We explore here to what extent CO (or its isotopologues) can be used as a tracer of gas mass. We use a detailed gas-grain chemical model and study the evolution of the disk composition, starting from a dense pre-stellar core composition. We explore a range of disk temperature profiles, cosmic rays ionization rates, and disk ages for a disk model representative of T Tauri stars. At the high densities that prevail in disks, we find that, due to fast reactions on grain surfaces, CO can be converted to less volatile forms (principally s-CO$_2$, and to a lesser extent s-CH$_4$) instead of being evaporated over a wide range of temperature. The canonical gas-phase abundance of 10$^{-4}$ is only reached above about 30-35 K. The dominant Carbon bearing entity depends on the temperature structure and age of the disk. The chemical evolution of CO is also sensitive to the cosmic rays ionization rate. Larger gas phase CO abundances are found in younger disks. Initial conditions, such as parent cloud age and density, have a limited impact. This study reveals that CO gas-phase abundance is heavily dependent on grain surface processes, which remain very incompletely understood so far. The strong dependence on dust temperature profile makes CO a poor tracer of the gas-phase content of disks., Comment: 11 pages, 9 figures, Accepted for publication in A&A

We attempt to determine the molecular composition of disks around young low-mass stars in the $\rho$ Oph region and to compare our results with a similar study performed in the Taurus-Auriga region. We used the IRAM 30 m telescope to perform a sensitive search for CN N=2-1 in 29 T Tauri stars located in the $\rho$ Oph and upper Scorpius regions. $^{13}$CO J=2-1 is observed simultaneously to provide an indication of the level of confusion with the surrounding molecular cloud. The bandpass also contains two transitions of ortho-H$_2$CO, one of SO, and the C$^{17}$O J=2-1 line, which provides complementary information on the nature of the emission. Contamination by molecular cloud in $^{13}$CO and even C$^{17}$O is ubiquitous. The CN detection rate appears to be lower than for the Taurus region, with only four sources being detected (three are attributable to disks). H$_2$CO emission is found more frequently, but appears in general to be due to the surrounding cloud. The weaker emission than in Taurus may suggest that the average disk size in the $\rho$ Oph region is smaller than in the Taurus cloud. Chemical modeling shows that the somewhat higher expected disk temperatures in $\rho$ Oph play a direct role in decreasing the CN abundance. Warmer dust temperatures contribute to convert CN into less volatile forms. In such a young region, CN is no longer a simple, sensitive tracer of disks, and observations with other tracers and at high enough resolution with ALMA are required to probe the gas disk population., Comment: 18 pages, 5 figures, accepted for publication in A&A