Molecular line radiative transfer in protoplanetary disks: Monte Carlo simulations versus approximate methods

We analyze the line radiative transfer in protoplanetary disks using several approximate methods and a well-tested Accelerated Monte Carlo code. A low-mass flaring disk model with uniform as well as stratified molecular abundances is adopted. Radiative transfer in low and high rotational lines of CO, C18O, HCO+, DCO+, HCN, CS, and H2CO is simulated. The corresponding excitation temperatures, synthetic spectra, and channel maps are derived and compared to the results of the Monte Carlo calculations. A simple scheme that describes the conditions of the line excitation for a chosen molecular transition is elaborated. We find that the simple LTE approach can safely be applied for the low molecular transitions only, while it significantly overestimates the intensities of the upper lines. In contrast, the Full Escape Probability (FEP) approximation can safely be used for the upper transitions (J_{\rm up} \ga 3) but it is not appropriate for the lowest transitions because of the maser effect. In general, the molecular lines in protoplanetary disks are partly subthermally excited and require more sophisticated approximate line radiative transfer methods. We analyze a number of approximate methods, namely, LVG, VEP (Vertical Escape Probability) and VOR (Vertical One Ray) and discuss their algorithms in detail. In addition, two modifications to the canonical Monte Carlo algorithm that allow a significant speed up of the line radiative transfer modeling in rotating configurations by a factor of 10--50 are described.Comment: 47 pages, 12 figures, accepted for publication in Ap

Investigating planet formation in circumstellar disks: CARMA observations of RY Tau and DG Tau

(Abridged) We present CARMA observations of the thermal dust emission from the circumstellar disks around the young stars RYTau and DGTau at wavelengths of 1.3mm and 2.8mm. The angular resolution of the maps is as high as 0.15arcsec, or 20AU at the distance of the Taurus cloud, which is a factor of 2 higher than has been achieved to date at these wavelengths. The unprecedented detail of the resulting disk images enables us to address three important questions related to the formation of planets. (1) What is the radial distribution of the circumstellar dust? (2) Does the dust emission show any indication of gaps that might signify the presence of (proto-)planets? (3) Do the dust properties depend on the orbital radius? We find that modeling the disk surface density in terms of either a classical power law or the similarity solution for viscous disk evolution, reproduces the observations well. The 1.3mm image from RYTau shows two peaks separated by 0.2arcsec with a decline in the dust emission toward the stellar position, which is significant at about 2-4sigma. For both RYTau and DGTau, the dust emission at radii larger than 15 AU displays no significant deviation from an unperturbed viscous disk model. In particular, no radial gaps in the dust distribution are detected. Under reasonable assumptions, we exclude the presence of planets more massive than 5 Jupiter masses orbiting either star at distances between about 10 and 60 AU. The radial variation of the dust opacity slope, beta, was investigated by comparing the 1.3mm and 2.8mm observations. We find mean values of beta of 0.5 and 0.7 for DGTau and RYTau respectively. Variations in beta are smaller than 0.7 between 20 and 70 AU. These results confirm that the circumstellar dust throughout these disks differs significantly from dust in the interstellar medium.Comment: ApJ in press

Structure and evolution of pre-main sequence circumstellar disks

We present new sub-arcsecond (0.7'') Combined Array for Research in Millimeter-wave Astronomy (CARMA) observations of the 1.3 mm continuum emission from circumstellar disks around 11 low and intermediate mass pre-main sequence stars. High resolution observations for 3 additional sources were obtained from literature. In all cases the disk emission is spatially resolved. We adopt a self consistent accretion disk model based on the similarity solution for the disk surface density and constrain the dust radial density distribution on spatial scales of about 40 AU. Disk surface densities appear to be correlated with the stellar ages where the characteristic disk radius increases from ~ 20 AU to 100 AU over about 5 Myr. This disk expansion is accompanied by a decrease in the mass accretion rate, suggesting that our sample disks form an evolutionary sequence. Interpreting our results in terms of the temporal evolution of a viscous $\alpha$-disk, we estimate (i) that at the beginning of the disk evolution about 60% of the circumstellar material was located inside radii of 25--40 AU, (ii) that disks formed with masses from 0.05 to 0.4 M_{\sun} and (iii) that the viscous timescale at the disk initial radius is about 0.1-0.3 Myr. Viscous disk models tightly link the surface density $\Sigma(R)$ with the radial profile of the disk viscosity $\nu(R) \propto R^{\gamma}$. We find values of $\gamma$ ranging from -0.8 to 0.8, suggesting that the viscosity dependence on the orbital radius can be very different in the observed disks. Adopting the $\alpha$ parameterization for the viscosity, we argue that $\alpha$ must decrease with the orbital radius and that it may vary between 0.5 and $10^{-4}$. (abridged)Comment: Accepted for publication in The Astrophysical Journal, 43 pages, 18 figures, Typo in the author name correcte

Chemistry in Disks. IV. Benchmarking gas-grain chemical models with surface reactions

Abridged: We detail and benchmark two sophisticated chemical models developed by the Heidelberg and Bordeaux astrochemistry groups. The main goal of this study is to elaborate on a few well-described tests for state-of-the-art astrochemical codes covering a range of physical conditions and chemical processes, in particular those aimed at constraining current and future interferometric observations of protoplanetary disks. We consider three physical models: a cold molecular cloud core, a hot core, and an outer region of a T Tauri disk. Our chemical network (for both models) is based on the original gas-phase osu_03_2008 ratefile and includes gas-grain interactions and a set of surface reactions for the H-, O-, C-, S-, and N-bearing molecules. The benchmarking is performed with the increasing complexity of the considered processes: (1) the pure gas-phase chemistry, (2) the gas-phase chemistry with accretion and desorption, and (3) the full gas-grain model with surface reactions. Using atomic initial abundances with heavily depleted metals and hydrogen in its molecular form, the chemical evolution is modeled within 10^9 years. The time-dependent abundances calculated with the two chemical models are essentially the same for all considered physical cases and for all species, including the most complex polyatomic ions and organic molecules. This result however required a lot of efforts to make all necessary details consistent through the model runs, e.g. definition of the gas particle density, density of grain surface sites, the strength and shape of the UV radiation field, etc. The reference models and the benchmark setup, along with the two chemical codes and resulting time-dependent abundances are made publicly available in the Internet: http://www.mpia.de/homes/semenov/Chemistry_benchmark/home.htmlComment: 12 pages, 5 figures, 3 tables, accepted for publication in A&A. The relevant web page is: http://www.mpia.de/homes/semenov/Chemistry_benchmark/home.htm

Chemistry in Disks. III. -- Photochemistry and X-ray driven chemistry probed by the ethynyl radical (CCH) in DM Tau, LkCa 15, and MWC 480

We studied several representative circumstellar disks surrounding the Herbig Ae star MWC 480 and the T Tauri stars LkCa 15 and DM Tau at (sub-)millimeter wavelengths in lines of CCH. Our aim is to characterize photochemistry in the heavily UV-irradiated MWC 480 disk and compare the results to the disks around cooler T Tauri stars. We detected and mapped CCH in these disks with the IRAM Plateau de Bure Interferome- ter in the C- and D-configurations in the (1-0) and (2-1) transitions. Using an iterative minimization technique, the CCH column densities and excitation conditions are con- strained. Very low excitation temperatures are derived for the T Tauri stars. These values are compared with the results of advanced chemical modeling, which is based on a steady-state flared disk structure with a vertical temperature gradient, and a gas- grain chemical network with surface reactions. Both model and observations suggest that CCH is a sensitive tracer of the X-ray and UV irradiation. The predicted radial dependency and source to source variations of CCH column densities qualitatively agree with the observed values, but the predicted column densities are too low by a factor of several. The chemical model fails to reproduce high concentrations of CCH in very cold disk midplane as derived from the observed low excitation condition for both the (1-0) and (2-1) transitions

HERSCHEL/PACS SURVEY OF PROTOPLANETARY DISKS IN TAURUS/AURIGA-OBSERVATIONS OF [O I] AND [C II], AND FAR-INFRARED CONTINUUM

<p>The Herschel Space Observatory was used to observe similar to 120 pre-main-sequence stars in Taurus as part of the GASPS Open Time Key project. Photodetector Array Camera and Spectrometer was used to measure the continuum as well as several gas tracers such as [O I] 63 mu m, [O I] 145 mu m, [C II] 158 mu m, OH, H2O, and CO. The strongest line seen is [O I] at 63 mu m. We find a clear correlation between the strength of the [O I] 63 mu m line and the 63 mu m continuum for disk sources. In outflow sources, the line emission can be up to 20 times stronger than in disk sources, suggesting that the line emission is dominated by the outflow. The tight correlation seen for disk sources suggests that the emission arises from the inner disk (</p>