Detection of H_2 Pure Rotational Line Emission from the GG Tauri Binary System

We present the first detection of the low-lying pure rotational emission lines of H_2 from circumstellar disks around T Tauri stars, using the Short Wavelength Spectrometer on the Infrared Space Observatory. These lines provide a direct measure of the total amount of warm molecular gas in disks. The J = 2 → 0 S(0) line at 28.218 μm and the J = 3 → 1 S(1) line at 17.035 μm have been observed toward the double binary system GG Tau. Together with limits on the J = 5 → 3 S(3) and J = 7 → 5 S(5) lines, the data suggest the presence of gas at T_(kin) ≈ 110 ± 10 K with a mass of (3.6 ± 2.0) × 10^(-3) M_☉ (±3 σ). This amounts to ~3% of the total gas + dust mass of the circumbinary disk as imaged by millimeter interferometry, but it is larger than the estimated mass of the circumstellar disk(s). Possible origins for the warm gas seen in H_2 are discussed in terms of photon and wind-shock heating mechanisms of the circumbinary material, and comparisons with model calculations are made

Interferometric Observations of Formaldehyde in the Protoplanetary Disk around LkCa15

Emission from the $2_{12}-1_{11}$ line of H$_2$CO has been detected and marginally resolved toward LkCa15 by the Nobeyama Millimeter Array. The column density of H$_2$CO is higher than that observed in DM Tau and than predicted by theoretical models of disk chemistry; also the line-intensity profile is less centrally peaked than that for CO. A similar behavior is observed in other organic gaseous molecules in the LkCa 15 disk.Comment: 5 pages, 4 figures. accepted to PASJ (Publication of Astronomical Society of Japan

The gas temperature in flaring disks around pre-main sequence stars

A model is presented which calculates the gas temperature and chemistry in the surface layers of flaring circumstellar disks using a code developed for photon-dominated regions. Special attention is given to the influence of dust settling. It is found that the gas temperature exceeds the dust temperature by up to several hundreds of Kelvins in the part of the disk that is optically thin to ultraviolet radiation, indicating that the common assumption that Tgas=Tdust is not valid throughout the disk. In the optically thick part, gas and dust are strongly coupled and the gas temperature equals the dust temperature. Dust settling has little effect on the chemistry in the disk, but increases the amount of hot gas deeper in the disk. The effects of the higher gas temperature on several emission lines arising in the surface layer are examined. The higher gas temperatures increase the intensities of molecular and fine-structure lines by up to an order of magnitude, and can also have an important effect on the line shapes.Comment: 14 pages, 10 figures, accepted for publication in A&

Detection of H2 pure rotational line emission from the GG~Tau binary system

We present the first detection of the low-lying pure rotational emission lines of H2 from circumstellar disks around T~Tauri stars, using the Short Wavelength Spectrometer on the Infrared Space Observatory. These lines provide a direct measure of the total amount of warm molecular gas in disks. The J=2->0 S(0) line at 28.218 mum and the J=3->1 S(1) line at 17.035 mum have been observed toward the double binary system GG Tau. Together with limits on the J=5->3 S(3) and J=7->5 S(5) lines, the data suggest the presence of gas at T_kin=110+-10 K with a mass of (3.6+-2.0)x10^-3 M_sol (3sigma). This amounts to ~3% of the total gas + dust mass of the circumbinary disk as imaged by millimeter interferometry, but is larger than the estimated mass of the circumstellar disk(s). Possible origins for the warm gas seen in H2 are discussed in terms of photon and wind-shock heating mechanisms of the circumbinary material, and comparisons with model calculations are made.Comment: 14 pages including 1 figure. To appear in Astrophysical Journal Letter

H_2 and CO Emission from Disks around T Tauri and Herbig Ae Pre-Main-Sequence Stars and from Debris Disks around Young Stars: Warm and Cold Circumstellar Gas

We present ISO Short-Wavelength Spectrometer observations of H_2 pure-rotational line emission from the disks around low- and intermediate-mass pre-main-sequence stars as well as from young stars thought to be surrounded by debris disks. The pre-main-sequence sources have been selected to be isolated from molecular clouds and to have circumstellar disks revealed by millimeter interferometry. We detect "warm" (T ≈100-200 K) H_2 gas around many sources, including tentatively the debris-disk objects. The mass of this warm gas ranges from ~ 10^(-4) M_☉ up to 8 x 10^(-3) and can constitute a nonnegligible fraction of the total disk mass. Complementary single-dish ^(12)CO 3-2/^(13)CO 3-2, and ^(12)CO 6-5 observations have been obtained as well. These transitions probe cooler gas at T ≈ 20-80 K. Most objects show a double-peaked CO emission profile characteristic of a disk in Keplerian rotation, consistent with interferometer data on the lower J lines. The ratios of the ^(12)CO 3-2/^(13)CO 3-2 integrated fluxes indicate that ^(12)CO 3-2 is optically thick but that ^(13)CO 3-2 is optically thin or at most moderately thick. The ^(13)CO 3-2 lines have been used to estimate the cold gas mass. If a H_2/CO conversion factor of 1 x 10^(-4) is adopted, the derived cold gas masses are factors of 10-200 lower than those deduced from 1.3 millimeter dust emission assuming a gas/dust ratio of 100, in accordance with previous studies. These findings confirm that CO is not a good tracer of the total gas content in disks since it can be photodissociated in the outer layers and frozen onto grains in the cold dense part of disks, but that it is a robust tracer of the disk velocity field. In contrast, H_2 can shield itself from photodissociation even in low-mass "optically thin" debris disks and can therefore survive longer. The warm gas is typically 1%-10% of the total mass deduced from millimeter continuum emission, but it can increase up to 100% or more for the debris-disk objects. Thus, residual molecular gas may persist into the debris-disk phase. No significant evolution in the H_2 CO, or dust masses is found for stars with ages in the range of 10^6-10^7 yr, although a decrease is found for the older debris-disk star β Pictoris. The large amount of warm gas derived from H_2 raises the question of the heating mechanism(s). Radiation from the central star as well as the general interstellar radiation field heat an extended surface layer of the disk, but existing models fail to explain the amount of warm gas quantitatively. The existence of a gap in the disk can increase the area of material influenced by radiation. Prospects for future observations with ground- and space-borne observations are discussed

Detection of DCO+ in a circumstellar disk

We report the first detection of DCO+ in a circumstellar disk. The DCO+ J=5-4 line at 360.169 GHz is observed with the 15m James Clerk Maxwell Telescope in the disk around the pre-main sequence star TW Hya. Together with measurements of the HCO+ and H13CO+ J=4-3 lines, this allows an accurate determination of the DCO+/HCO+ ratio in this disk. The inferred value of 0.035+-0.015 is close to that found in cold pre-stellar cores and is somewhat higher than that measured in the envelope around the low-mass protostar IRAS 16293 -2422. It is also close to the DCN/HCN ratio obtained for pristine cometary material in the jet of comet Hale-Bopp. The observed DCO+/HCO+ ratio for TW Hya is consistent with theoretical models of disks which consider gas-phase fractionation processes within a realistic 2-D temperature distribution and which include the effects of freeze-out onto grains

Warm Molecular Layers in Protoplanetary Disks

We have investigated molecular distributions in protoplanetary disks, adopting a disk model with a temperature gradient in the vertical direction. The model produces sufficiently high abundances of gaseous CO and HCO+ to account for line observations of T Tauri stars using a sticking probability of unity and without assuming any non-thermal desorption. In regions of radius R > 10 AU, with which we are concerned, the temperature increases with increasing height from the midplane. In a warm intermediate layer, there are significant amounts of gaseous molecules owing to thermal desorption and efficient shielding of ultraviolet radiation by the flared disk. The column densities of HCN, CN, CS, H2CO, HNC and HCO+ obtained from our model are in good agreement with the observations of DM Tau, but are smaller than those of LkCa15. Molecular line profiles from our disk models are calculated using a 2-dimensional non-local-thermal-equilibrium (NLTE) molecular-line radiative transfer code for a direct comparison with observations. Deuterated species are included in our chemical model. The molecular D/H ratios in the model are in reasonable agreement with those observed in protoplanetary disks.Comment: 11 pages, Latex (aa.cls), to be published in Astronomy and Astrophysic

Radiation thermo-chemical models of protoplanetary disks I. Hydrostatic disk structure and inner rim

This paper introduces a new disk code, called ProDiMo, to calculate the thermo-chemical structure of protoplanetary disks and to interpret gas emission lines from UV to sub-mm. We combine frequency-dependent 2D dust continuum radiative transfer, kinetic gas-phase and UV photo-chemistry, ice formation, and detailed non-LTE heating & cooling balance with the consistent calculation of the hydrostatic disk structure. We include FeII and CO ro-vibrational line heating/cooling relevant for the high-density gas close to the star, and apply a modified escape probability treatment. The models are characterized by a high degree of consistency between the various physical, chemical and radiative processes, where the mutual feedbacks are solved iteratively. In application to a T Tauri disk extending from 0.5AU to 500AU, the models are featured by a puffed-up inner rim and show that the dense, shielded and cold midplane (z/r<0.1, Tg~Td) is surrounded by a layer of hot (5000K) and thin (10^7 to 10^8 cm^-3) atomic gas which extends radially to about 10AU, and vertically up to z/r~0.5. This layer is predominantly heated by the stellar UV (e.g. PAH-heating) and cools via FeII semi-forbidden and OI 630nm optical line emission. The dust grains in this "halo" scatter the star light back onto the disk which impacts the photo-chemistry. The more distant regions are characterized by a cooler flaring structure. Beyond 100AU, Tgas decouples from Tdust even in the midplane and reaches values of about Tg~2Td. Our models show that the gas energy balance is the key to understand the vertical disk structure. Models calculated with the assumption Tg=Td show a much flatter disk structure.Comment: 24 pages, 14 figures, 120 equations, accepted by A&A, download a high-resolution version from http://www.roe.ac.uk/~ptw/prodimo1_article.pd

Vertical structure models of T Tauri and Herbig Ae/Be disks

In this paper we present detailed models of the vertical structure (temperature and density) of passive irradiated circumstellar disks around T Tauri and Herbig Ae/Be stars. In contrast to earlier work, we use full frequency- and angle-dependent radiative transfer instead of the usual moment equations. We find that this improvement of the radiative transfer has strong influence on the resulting vertical structure of the disk, with differences in temperature as large as 70 %. However, the spectral energy distribution (SED) is only mildly affected by this change. In fact, the SED compares reasonably well with that of improved versions of the Chiang & Goldreich (CG) model. This shows that the latter is a reasonable model for the SED, in spite of its simplicity. It also shows that from the SED alone, little can be learned about the vertical structure of a passive circumstellar disk. The molecular line emission from these disks is more sensitive to the vertical temperature and density structure, and we show as an example how the intensity and profiles of various CO lines depend on the adopted disk model. The models presented in this paper can also serve as the basis of theoretical studies of e.g. dust coagulation and settling in disks.Comment: 12 pages, 15 figures, accepted for publication in A&

LIME - a flexible, non-LTE line excitation and radiation transfer method for millimeter and far-infrared wavelengths

We present a new code for solving the molecular and atomic excitation and radiation transfer problem in a molecular gas and predicting emergent spectra. This code works in arbitrary three dimensional geometry using unstructured Delaunay latices for the transport of photons. Various physical models can be used as input, ranging from analytical descriptions over tabulated models to SPH simulations. To generate the Delaunay grid we sample the input model randomly, but weigh the sample probability with the molecular density and other parameters, and thereby we obtain an average grid point separation that scales with the local opacity. Our code does photon very efficiently so that the slow convergence of opaque models becomes traceable. When convergence between the level populations, the radiation field, and the point separation has been obtained, the grid is ray-traced to produced images that can readily be compared to observations. Because of the high dynamic range in scales that can be resolved using this type of grid, our code is particularly well suited for modeling of ALMA data. Our code can furthermore deal with overlapping lines of multiple molecular and atomic species.Comment: 13 pages, 12 figures, Accepted by A&A on 06/08/201