Search

{Sulfur Monoxide is known to be a good shock tracer in molecular clouds and protostar environments, but its abundance is difficult to reproduce even using state-of-the-art astrochemical models. } {We investigate the properties of the observed SO emission in the protoplanetary disk of AB Auriga, a Herbig Ae star of 2.4 $\Msun$ located at 156 pc. The AB Aur system is unique because it exhibits a dust trap and at least one young putative planet orbiting at about 30 au from the central star.} {We reduced ALMA archival data (projects 2019.1.00579.S, 2021.1.00690.S and 2021.1.01216.S) and analyzed the three detected SO lines (SO $6_5-5_4$, $6_7-5_6$ and $5_6-4_5$). We also present C$^{17}$O and C$^{18}$O 2-1 data to complement the interpretation of the SO data.} {For the three SO lines, the maximum SO emission in the ring is not located in the dust trap. Moreover, the inner radius of the SO ring is significantly larger than the CO emission inner radius, $\sim 160$ au versus $\sim 90$ au. The SO emission traces gas located in part beyond the dust ring. This emission likely originates from shocks at the interface of the outer spirals, observed in CO and scattered light emission, and the molecular and dust ring. SO is also detected inside the cavity, at a radius $\sim 20-30$ au and with a rotation velocity compatible with the proto-planet P1. We speculate that this SO emission originates from accretion shocks onto the circumplanetary disk of the putative proto-planet P1.} {These observations confirm that SO is a good tracer of shocks in protoplanetary disks and could be a powerful, new tool to detect embedded (proto-)planets.} Conclusions. These observations confirm that SO is a good tracer of shocks in protoplanetary disks and could be a powerful, new tool to detect embedded (proto-)planets.

We aim to constrain the gas density and temperature distributions as well as gas masses in several T Tauri protoplanetary disks located in Taurus. We use the 12CO, 13CO, and C18O (2-1) isotopologue emission observed at 0.9 with the IRAM NOrthern Extended Millimeter Array (NOEMA) as part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution PIs: P. Caselli & Th. Henning). Our sample consists of Class II disks with no evidence of strong radial substructures. We use thesedata to constrain the thermal and chemical structure of these disks through theoretical models for gas emission. To fit the combined optically thick and thin CO line data in Fourier space, we developed the DiskCheF code, which includes the parameterized disk physical structure, machine-learning (ML) accelerated chemistry, and the RADMC-3D line radiative transfer module. A key novelty of DiskCheF is the fast and feasible ML-based chemistry trained on the extended grid of the disk physical-chemical models precomputed with the ANDES2 code. This ML approach allows complex chemical kinetics models to be included in a time-consuming disk fitting without the need to run a chemical code. We present a novel approach to incorporate chemistry into disk modeling without the need to explicitly calculate a chemical network every time. Using this new disk modeling tool, we successfully fit the 12CO, 13CO, and C18O (2-1) data from the CI, CY, DL, DM, DN, and IQ Tau disks. The combination of optically thin and optically thick CO lines allows us to simultaneously constrain the disk temperature and mass distribution, and derive the CO-based gas masses. These values are in reasonable agreement with the disk dust masses rescaled by a factor of 100 as well as with other indirect gas measurements., Comment: 25 pages, 19 figures

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

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

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

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

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

Aims: We attempt to determine the masses of single or multiple young T Tauri and HAeBe stars from the rotation of their Keplerian disks. Methods:We used the IRAM PdBI interferometer to perform arcsecond resolution images of the CN N=2-1 transition with good spectral resolution. Integrated spectra from the 30-m radiotelescope show that CN is relatively unaffected by contamination from the molecular clouds. Our sample includes 12 sources, among which isolated stars like DM Tau and MWC 480 are used to demonstrate the method and its accuracy. We derive the dynamical mass by fitting a disk model to the emission, a process giving M/D the mass to distance ratio. We also compare the CN results with higher resolution CO data, that are however affected by contamination. Results: All disks are found in nearly perfect Keplerian rotation. We determine accurate masses for 11 stars, in the mass range 0.5 to 1.9 solar masses. The remaining one, DG Tau B, is a deeply embedded object for which CN emission partially arises from the outflow. With previous determination, this leads to 14 (single) stars with dynamical masses. Comparison with evolutionary tracks, in a distance independent modified HR diagram, show good overall agreement (with one exception, CW Tau), and indicate that measurement of effective temperatures are the limiting factor. The lack of low mass stars in the sample does not allow to distinguish between alternate tracks., Comment: Astronomy & Astrophysics in press, 20 pages

We aim at unveiling the observational imprint of physical mechanisms that govern planetary formation in young, multiple systems. In particular, we investigate the impact of tidal truncation on the inner circumstellar disks. We observed the emblematic system GG Tau at high-angular resolution: a hierarchical quadruple system composed of low-mass T Tauri binary stars surrounded by a well-studied, massive circumbinary disk in Keplerian rotation. We used the near-IR 4-telescope combiner PIONIER on the VLTI and sparse-aperture-masking techniques on VLT/NaCo to probe this proto-planetary system at sub-AU scales. We report the discovery of a significant closure-phase signal in H and Ks bands that can be reproduced with an additional low-mass companion orbiting GG Tau Ab, at a (projected) separation rho = 31.7 +/- 0.2mas (4.4 au) and PA = 219.6 +/- 0.3deg. This finding offers a simple explanation for several key questions in this system, including the missing-stellar-mass problem and the asymmetry of continuum emission from the inner dust disks observed at millimeter wavelengths. Composed of now five co-eval stars with 0.02 <= Mstar <= 0.7 Msun, the quintuple system GG Tau has become an ideal test case to constrain stellar evolution models at young ages (few 10^6yr)., Comment: 5pages, 3 figures, 1 appendix (online material)

(abridged) Most Class II sources (of nearby star forming regions) are surrounded by disks with weak millimeter continuum emission. These "faint" disks may hold clues to the disk dissipation mechanism. We attempt to determine the characteristics of such faint disks around classical T Tauri stars, and to explore the link between disk faintness and the proposed disk dispersal mechanisms (accretion, viscous spreading, photo-evaporation, planetary system formation). We performed high-angular resolution (0.3") imaging of a small sample of disks (9 sources) with low 1.3mm continuum flux (mostly <30 mJy) with the IRAM Plateau de Bure interferometer and simultaneously searched for 13CO (or CO) J=2-1 line emission. Using a simple parametric disk model, we determine characteristic sizes of the disks, in dust and gas, and we constrain surface densities in the central 50 AU. All disks are much smaller than the bright disks imaged so far, both in continuum and 13CO lines (5 detections). In continuum, half of the disks are very small, with characteristic radii less than 10AU, but still have high surface density values. Small sizes appear to be the main cause for the low disk luminosity. Direct evidence for grain growth is found for the three disks that are sufficiently resolved. Low continuum opacity is attested in two systems only, but we cannot firmly distinguish between a low gas surface density and a lower dust emissivity resulting from grain growth. We report a tentative discovery of a 20 AU radius cavity in DS Tau, bringing the proportion of transitional disks to a similar value to that of brighter sources, but cavities cannot explain the low mm flux. This study highlights a category of very compact dust disks, still exhibiting high surface densities, which may represent up to 25 % of the whole disk population, but its origin is unclear with the current data alone., Comment: 15 pages Accepted for publication in Astronomy & Astrophysics

Investigating the dynamical evolution of dust grains in proto-planetary disks is a key issue to understand how planets should form. We identify under which conditions dust settling can be constrained by high angular resolution observations at mm wavelengths, and which observational strategies are suited for such studies. Exploring a large range of models, we generate synthetic images of disks with different degrees of dust settling, and simulate high angular resolution (~ 0.05-0.3") ALMA observations of these synthetic disks. The resulting data sets are then analyzed blindly with homogeneous disk models (where dust and gas are totally mixed) and the derived disk parameters are used as tracers of the settling factor. Our dust disks are partially resolved by ALMA and present some specific behaviors on radial and mainly vertical directions, which can be used to quantify the level of settling. We find out that an angular resolution better than or equal to ~ 0.1" (using 2.3 km baselines at 0.8mm) allows us to constrain the dust scale height and flaring index with sufficient precision to unambiguously distinguish between settled and non-settled disks, provided the inclination is close enough to edge-on (i >= 75{\deg}). Ignoring dust settling and assuming hydrostatic equilibrium when analyzing such disks affects the derived dust temperature and the radial dependency of the dust emissivity index. The surface density distribution can also be severely biased at the highest inclinations. However, the derived dust properties remain largely unaffected if the disk scale height is fitted separately. ALMA has the potential to test some of the dust settling mechanisms, but for real disks, deviations from ideal geometry (warps, spiral waves) may provide an ultimate limit on the dust settling detection., Comment: 14 pages, 11 figures

Turbulence is thought to be a key driver of the evolution of protoplanetary disks, regulating the mass accretion process, the transport of angular momentum, and the growth of dust particles. We intend to determine the magnitude of the turbulent motions in the outer parts (> 100 AU) of the disk surrounding DM Tau. Turbulent motions can be constrained by measuring the nonthermal broadening of line emission from heavy molecules. We used the IRAM Plateau de Bure interferometer to study emission from the CS molecule in the disk of DM Tau. High spatial (1.4 x 1 ") and spectral resolution (0.126 km/s) CS J=3-2 images provide constraints on the molecule distribution and velocity structure of the disk. A low sensitivity CS J=5-4 image was used in conjunction to evaluate the excitation conditions. We analyzed the data in terms of two parametric disk models, and compared the results with detailed time-dependent chemical simulations. The CS data confirm the relatively low temperature suggested by observations of other simple molecules. The intrinsic linewidth derived from the CS J=3-2 data is much larger than expected from pure thermal broadening. The magnitude of the derived nonthermal component depends only weakly on assumptions about the location of the CS molecules with respect to the disk plane. Our results indicate turbulence with a Mach number around 0.4 - 0.5 in the molecular layer. Geometrical constraints suggest that this layer is located near one scale height, in reasonable agreement with chemical model predictions., Comment: 9 pages Astronomy & Astrophysics, in press

We use the IRAM 30-m telescope to perform a sensitive search for CN N=2-1 in 42 T Tauri or Herbig Ae systems located mostly in the Taurus-Auriga region. $^{13}$CO J=2-1 is observed simultaneously to indicate 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 provide complementary information on the nature of the emission. While $^{13}$CO is in general dominated by residual emission from the cloud, CN exhibits a high disk detection rate $> 50$% in our sample. We even report CN detection in stars for which interferometric searches failed to detect $^{12}$CO, presumably because of obscuration by a foreground, optically thick, cloud. Comparison between CN and o-H$_2$CO or SO line profiles and intensities divide the sample in two main categories. Sources with SO emission are bright and have strong H$_2$CO emission, leading in general to [H$_2$CO/CN]$ > 0.5$. Furthermore, their line profiles, combined with a priori information on the objects, suggest that the emission is coming from outflows or envelopes rather than from a circumstellar disk. On the other hand, most sources have [H$_2$CO/CN]$ < 0.3$, no SO emission, and some of them exhibit clear double-peaked profiles characteristics of rotating disks. In this second category, CN is likely tracing the proto-planetary disks. From the line flux and opacity derived from the hyperfine ratios, we constrain the outer radii of the disks, which range from 300 to 600 AU. The overall gas disk detection rate (including all molecular tracers) is $\sim 68%$, and decreases for fainter continuum sources. This study shows that gas disks, like dust disks, are ubiquitous around young PMS stars in regions of isolated star formation, and that a large fraction of them have $R > 300$ AU., Comment: 31 pages (including 59 figures)

The radio-quiet quasar BR1202-0725 (z=4.695) is a remarkable source with a bright Northwest (NW) companion detected at submm and radio wavelengths but invisible in the optical. In the absence of amplification by gravitational lensing, BR1202-0725 would be the most luminous binary CO and FIR source in the Universe. In this paper, we report observations with the IRAM Plateau de Bure interferometer of BR1202-0725 in the redshifted emission of the CO(5-4) and (7-6) lines, the [C I](3P2-3P1) line, a high angular resolution (0.3 x 0.8 arcsec) 1.3 mm map of the rest-frame, far-IR dust continuum, and a search for the CO(11-10) line. We compare these results with recent ALMA data in the [C II] line. Both the quasar host galaxy and its NW companion are spatially resolved in the molecular line emission and the dust continuum. The CO profile of the NW companion is very broad with a full width at half maximum of 1000 +/- 130 km/s, compared to 360 +/- 40 km/s for the quasar host galaxy to the Southeast (SE). The difference in linewidths and center velocities, and the absence of any lens candidate or arc-like structure in the field, at any wavelength, show that the obscured NW galaxy and the SE quasar host galaxy cannot be lensed images of the same object. Instead, we find morphological and kinematic evidence for sub-structures in both the NW and SE sources. We interpret these results as strong indications that the BR1202-0725 complex is a group of young, interacting, and highly active starburst galaxies., Comment: Minor changes to agree with published version

Molecular line emission from protoplanetary disks is a powerful tool to constrain their physical and chemical structure. Nevertheless, only a few molecules have been detected in disks so far. We take advantage of the enhanced capabilities of the IRAM 30m telescope by using the new broad band correlator (FTS) to search for so far undetected molecules in the protoplanetary disks surrounding the TTauri stars DM Tau, GO Tau, LkCa 15 and the Herbig Ae star MWC 480. We report the first detection of HC3N at 5 sigma in the GO Tau and MWC 480 disks with the IRAM 30-m, and in the LkCa 15 disk (5 sigma), using the IRAM array, with derived column densities of the order of 10^{12}cm^{-2}. We also obtain stringent upper limits on CCS (N < 1.5 x 10^{12} cm^{-3}). We discuss the observational results by comparing them to column densities derived from existing chemical disk models (computed using the chemical code Nautilus) and based on previous nitrogen and sulfur-bearing molecule observations. The observed column densities of HC3N are typically two orders of magnitude lower than the existing predictions and appear to be lower in the presence of strong UV flux, suggesting that the molecular chemistry is sensitive to the UV penetration through the disk. The CCS upper limits reinforce our model with low elemental abundance of sulfur derived from other sulfur-bearing molecules (CS, H2S and SO)., Comment: 8 pages, 4 figures, 3 tables, Accepted for publication in ApJ

Circumstellar disks are characteristic for star formation and vanish during the first few Myr of stellar evolution. During this time planets are believed to form in the dense midplane by growth, sedimentation and aggregation of dust. Indicators of disk evolution, such as holes and gaps, can be traced in the spectral energy distribution (SED) and spatially resolved images. We aim to construct a self-consistent model of HH 30 by fitting all available continuum observations simultaneously. New data sets not available in previous studies, such as high-resolution interferometric imaging with the Plateau de Bure Interferometer (PdBI) at lambda = 1.3 mm and SED measured with IRS on the Spitzer Space Telescope in the mid-infrared, put strong constraints on predictions and are likely to provide new insights into the evolutionary state of this object. A parameter study based on simulated annealing was performed to find unbiased best-fit models for independent observations made in the wavelength domain lambda ~ 1 micron ... 4 mm. The method essentially creates a Markov chain through parameter space by comparing predictions generated by our self-consistent continuum radiation transfer code MC3D with observations. We present models of the edge-on circumstellar disk of HH 30 based on observations from the near-infrared to mm-wavelengths that suggest the presence of an inner depletion zone with about 45 AU radius and a steep decline of mm opacity beyond 140 AU. Our modeling indicates that several modes of dust evolution such as growth, settling, and radial migration are taking place in this object. High-resolution observations of HH 30 at different wavelengths with next-generation observatories such as ALMA and JWST will enable the modeling of inhomogeneous dust properties and significantly expand our understanding of circumstellar disk evolution., Comment: Accepted for publication in Astronomy and Astrophysics

The chemistry of proto-planetary disks is thought to be dominated by two major processes: photodissociation near the disk surface, and depletion on dust grains in the disk mid-plane, resulting in a layered structure with molecules located in a warm layer above the disk mid-plane. We attempt here to confront this warm molecular layer model prediction with the distribution of two key molecules for dissociation processes: CN and HCN. Using the IRAM Plateau de Bure interferometer, we obtained high spatial and spectral resolution images of the CN J=2-1 and HCN J=1-0 lines in the disks surrounding the two T-Tauri DM Tau and LkCa 15 and the Herbig Ae MWC 480. Disk properties are derived assuming power law distributions. The hyperfine structure of the observed transitions allows us to constrain the line opacities and excitation temperatures. We compare the observational results with predictions from existing chemical models, and use a simple PDR model (without freeze-out of molecules on grains and surface chemistry) to illustrate dependencies on UV field strength, grain size and gas-to-dust ratio. We also evaluate the impact of Lyman alpha radiation. The temperature ordering follows the trend found from CO lines, with DM Tau being the coldest object and MWC 480 the warmest. Although CN indicates somewhat higher excitation temperatures than HCN, the derived values in the T-Tauri disks are very low (8-10 K). They agree with results obtained from CCH, and are in contradiction with thermal and chemical model predictions. These very low temperatures, as well as geometrical constraints, suggest that substantial amounts of CN and HCN remain in the gas phase close to the disk mid-plane, and that this mid-plane is quite cold. The observed CN/HCN ratio (5-10) is in better agreement with the existence of large grains, and possibly also a substantial contribution of Lyman alpha radiation., Comment: 22 pages, 20 figures, 5 tables. To be published in Astronomy and Astrophysics

The structure in density and temperature of protoplanetary disks surrounding low-mass stars is not yet well known. The protoplanetary disks mid-planes are expected to be very cold and thus depleted in molecules in gas phase, especially CO. Recent observations of molecules at very low apparent temperature (~ 6 K) challenge this current picture of the protoplanetary disk structures. We aim at constraining the physical conditions, and in particular the gas-phase CO abundance in the mid-plane of protoplanetary disks. The light molecule H2D+,is a tracer of cold and CO-depleted environment. It is therefore a good candidate to explore the disks mid-planes. We performed a deep search for H2D+ in the two well-known disks surrounding TW Hya and DM Tau using the APEX and JCMT telescopes. The analysis of the observations are done with DISKFIT, a radiative transfer code dedicated to disks. In addition, we used a chemical model describing deuterium chemistry to infer the implications of our observations on the level of CO depletion and on the ionization rate in the disks mid-plane. The ortho-H2D+(1(1,0))-1(1,1)) line at 372 GHz was not detected. Although our limit is three times better than previous observations, comparison with the chemical modeling indicates that it remains insufficient to put valuable constraints on the CO abundance in the disk mid-plane. Even with ALMA, the detection of H2D+ may not be straightforward, and H2D+ may not be a sufficiently sensitive tracer of the protoplanetary disks mid-plane, Comment: Submitted to A&A

Context: Dark cloud chemical models usually predict large amounts of O2, often above observational limits. Aims: We investigate the reason for this discrepancy from a theoretical point of view, inspired by the studies of Jenkins and Whittet on oxygen depletion. Methods: We use the gas-grain code Nautilus with an up-to-date gas-phase network to study the sensitivity of the molecular oxygen abundance to the oxygen elemental abundance. We use the rate coefficient for the reaction O + OH at 10 K recommended by the KIDA (KInetic Database for Astrochemistry) experts. Results: The updates of rate coefficients and branching ratios of the reactions of our gas-phase chemical network, especially N + CN and H3+ + O, have changed the model sensitivity to the oxygen elemental abundance. In addition, the gas-phase abundances calculated with our gas-grain model are less sensitive to the elemental C/O ratio than those computed with a pure gas-phase model. The grain surface chemistry plays the role of a buffer absorbing most of the extra carbon. Finally, to reproduce the low abundance of molecular oxygen observed in dark clouds at all times, we need an oxygen elemental abundance smaller than 1.6E-4. Conclusion: The chemistry of molecular oxygen in dense clouds is quite sensitive to model parameters that are not necessarily well constrained. That O2 abundance may be sensitive to nitrogen chemistry is an indication of the complexity of interstellar chemistry., Comment: 6 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

(Abridged) We attempt to characterize the radial distribution of dust in disks around a sample of young stars from an observational point of view, and, when possible, in a model-independent way, by using parametric laws. We used the IRAM PdBI interferometer to provide very high angular resolution (down to 0.4" in some sources) observations of the continuum at 1.3 mm and 3 mm around a sample of T Tauri stars in the Taurus-Auriga region. The sample includes single and multiple systems, with a total of 23 individual disks. We used track-sharing observing mode to minimize the biases. We fitted these data with two kinds of models: a "truncated power law" model and a model presenting an exponential decay at the disk edge ("viscous" model). ect evidence for tidal truncation is found in the multiple systems. The temperature of the mm-emitting dust is constrained in a few systems. Unambiguous evidence for large grains is obtained by resolving out disks with very low values of the dust emissivity index Beta. In most disks that are sufficiently resolved at two different wavelengths, we find a radial dependence of Beta, which appears to increase from low values (as low as 0) at the center to about 1.7 -- 2 at the disk edge. The same behavior could apply to all studied disks. It introduces further ambiguities in interpreting the brightness profile, because the regions with apparent Beta = 0 can also be interpreted as being optically thick when their brightness temperature is high enough. Despite the added uncertainty on the dust absorption coefficient, the characteristic size of the disk appears to increase with a higher estimated star age. These results provide the first direct evidence of the radial dependence of the grain size in proto-planetary disks., Comment: 27 pages plus 14 pages of figures

Gas and dust dissipation processes of proto-planetary disks are hardly known. Transition disks between Class II (proto-planetary disks) and Class III (debris disks) remain difficult to detect. We investigate the carbon chemistry of the peculiar CQ Tau gas disk. It is likely a transition disk because it exhibits weak CO emission with a relatively strong millimeter continuum, indicating that the disk might be currently dissipating its gas content. We used APEX to observe the two CI lines at 492GHz and 809 GHz in the disk orbiting CQ Tau. We compare the observations to several chemical model predictions. We focus our study on the influence of the stellar UV radiation shape and gas-to-dust ratio. We did not detect the CI lines. However, our upper limits are deep enough to exclude high-CI models. The only available models compatible with our limits imply very low gas-to-dust ratio, of the order of a few, only. These observations strengthen the hypothesis that CQ Tau is likely a transition disk and suggest that gas disappears before dust., Comment: 5 pages, 5 figures, accepted for publication in A&A

Recent sub-millimetric observations at the Plateau de Bure interferometer evidenced a cavity at ~ 46 AU in radius into the proto-planetary disk around the T Tauri star LkCa15 (V1079 Tau), located in the Taurus molecular cloud. Additional Spitzer observations have corroborated this result possibly explained by the presence of a massive (>= 5 MJup) planetary mass, a brown dwarf or a low mass star companion at about 30 AU from the star. We used the most recent developments of high angular resolution and high contrast imaging to search directly for the existence of this putative companion, and to bring new constraints on its physical and orbital properties. The NACO adaptive optics instrument at VLT was used to observe LkCa15 using a four quadrant phase mask coronagraph to access small angular separations at relatively high contrast. A reference star at the same parallactic angle was carefully observed to optimize the quasi-static speckles subtraction (limiting our sensitivity at less than 1.0). Although we do not report any positive detection of a faint companion that would be responsible for the observed gap in LkCa15's disk (25-30 AU), our detection limits start constraining its probable mass, semi-major axis and eccentricity. Using evolutionary model predictions, Monte Carlo simulations exclude the presence of low eccentric companions with masses M >= 6 M Jup and orbiting at a >= 100 AU with significant level of confidence. For closer orbits, brown dwarf companions can be rejected with a detection probability of 90% down to 80 AU (at 80% down to 60 AU). Our detection limits do not access the star environment close enough to fully exclude the presence of a brown dwarf or a massive planet within the disk inner activity (i.e at less than 30 AU). Only, further and higher contrast observations should unveil the existence of this putative companion inside the LkCa15 disk., Comment: 6 pages, 4 figures, accepted for publication in A&A

Because of new telescopes that will dramatically improve our knowledge of the interstellar medium, chemical models will have to be used to simulate the chemistry of many regions with diverse properties. To make these models more robust, it is important to understand their sensitivity to a variety of parameters. In this article, we report a study of the sensitivity of a chemical model of a cold dense core, with homogeneous and time-independent physical conditions, to variations in the following parameters: initial chemical inventory, gas temperature and density, cosmic-ray ionization rate, chemical reaction rate coefficients, and elemental abundances. From the results of the parameter variations, we can quantify the sensitivity of the model to each parameter as a function of time. Our results can be used in principle with observations to constrain some parameters for different cold clouds. We also attempted to use the Monte Carlo approach with all parameters varied collectively. Within the parameter ranges studied, the most critical parameters turn out to be the reaction rate coefficients at times up to 4e5 yr and elemental abundances at later times. At typical times of best agreement with observation, models are sensitive to both of these parameters. The models are less sensitive to other parameters such as the gas density and temperature. The improvement of models will require that the uncertainties in rate coefficients of important reactions be reduced. As the chemistry becomes better understood and more robust, it should be possible to use model sensitivities concerning other parameters, such as the elemental abundances and the cosmic ray ionization rate, to yield detailed information on cloud properties and history. Nevertheless, at the current stage, we cannot determine the best values of all the parameters simultaneously based on purely observational constraints., Comment: Accepted for publication in Astronomy & Astrophysics