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.

The results of single-dish observations of low- and high-J transitions of selected molecules from protoplanetary disks around two TTauri stars (LkCa15 and TWHya) and two HerbigAe stars (HD163296 and MWC480) are reported. Simple molecules such as CO, 13CO, HCO+, CN and HCN are detected. Several lines of H2CO are found toward the TTauri star LkCa15 but not in other objects. No CH3OH has been detected down to abundances of 10E-9 - 10E-8 with respect to H2. SO and CS lines have been searched for without success. Line ratios indicate that the molecular emission arises from dense 10E6 - 10E8 cm-3 and moderately warm (T ~ 20-40K) intermediate height regions of the disk atmosphere, in accordance with predictions from models of the chemistry in disks. The abundances of most species are lower than in the envelope around the solar-mass protostar IRAS 16293-2422. Freeze-out and photodissociation are likely causes of the depletion. DCO+ is detected toward TWHya, but not in other objects. The high inferred DCO+/HCO+ ratio of ~0.035 is consistent with models of the deuterium fractionation in disks which include strong depletion of CO. The inferred ionization fraction in the intermediate height regions as deduced from HCO+ is at least 10E-11 - 10E-10, comparable to that derived for the midplane from recent H2D+ observations. (abridged abstract)

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.

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.

The nature of the tenuous disks around A stars has raised considerable controversy in the literature during the past decade. The debate whether or not the disk around beta Pictoris contains gaseous molecular hydrogen is only the most recent example. Since CO is in general a poor tracer for the gas content of these low mass disks, we discuss here detailed emission line calculations for alternative tracers like C and C+, based on recent optically thin disk models by Kamp & van Zadelhoff (2001). The [CII] 157.7 mum line was searched toward Vega and beta Pictoris -- the most prominent A stars with disks -- using ISO LWS data, and a tentative detection is reported toward the latter object. From a comparison with emission line observations as well as absorption line studies of both stars, the gas-to-dust ratio is constrained to lie between 0.5 and 9 for beta Pictoris. For Vega the [CII] observations indicate an upper limit of 0.2 M_Earth for the disk gas mass. Predicted line intensities of C+ and C are presented for a range of models and appear promising species to trace the gas content in the disks around A stars with future instrumental capabilities (SOFIA, Herschel, APEX and ALMA). Searches for CO emission should focus on the J=3-2 line.

A new two-dimensional axi-symmetric ultraviolet radiative transfer code is presented, which is used to calculate photodissociation and ionization rates for use in chemistry models of flaring circumstellar disks. Scattering and absorption of photons from the central star and from the interstellar radiation field are taken into account. The molecules are effectively photodissociated in the surface layer of the disk, but can exist in the intermediate, moderately warm layers. A comparison has been made with an approximate 2D ray-tracing method and it was found that the latter underestimates the ultraviolet field and thus the molecular photodissociation rates below the disk surface. The full 2D results show significantly higher abundances of radicals such as CN and C2H than previous work, partly due to the fact that CO is dissociated to greater depths. Results for different stellar radiation fields are also presented. The CN/HCN ratio shows a strong dependence on the stellar spectrum, whereas other ratios such as HCO+/CO show only little variation.

Comparison is made between a number of independent computer programs for radiative transfer in molecular rotational lines. The test models are spherically symmetric circumstellar envelopes with a given density and temperature profile. The first two test models have a simple power law density distribution, constant temperature and a fictive 2-level molecule, while the other two test models consist of an inside-out collapsing envelope observed in rotational transitions of HCO+. For the 2-level molecule test problems all codes agree well to within 0.2%, comparable to the accuracy of the individual codes, for low optical depth and up to 2% for high optical depths (tau=4800). The problem of the collapsing cloud in HCO+ has a larger spread in results, ranging up to 12% for the J=4 population. The spread is largest at the radius where the transition from collisional to radiative excitation occurs. The resulting line profiles for the HCO+ J=4-3 transition agree to within 10%, i.e., within the calibration accuracy of most current telescopes. The comparison project and the results described in this paper provide a benchmark for future code development, and give an indication of the typical accuracy of present day calculations of molecular line transfer.

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.

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.

Observations of submillimeter lines of CO (up to J=6-5), HCO+, HCN (up to J=4-3) and their isotopes from circumstellar disks around low mass pre-main sequence stars are presented. The different line ratios and optical depths indicate that most of the observed line emission arises from an intermediate disk layer with high densities of 1E6-1E8 cm-3 and moderately warm temperatures in the outer regions. The data are compared with three disk models from the literature using a full 2D Monte Carlo radiative transfer code. The abundances of the molecules are constrained from the more optically thin 13C species and indicate depletions of 1-30 for LkCa 15 and very high depletions of > 100 for TW Hya with respect to dark cloud abundances. Evidence for significant freeze-out (factors 10 or larger) of CO and HCO+ onto grain surfaces at temperatures below 22 K is found, but the abundances of these molecules must also be low in the warmer upper layer, most likely as a result of photodissociation. A warm upper layer near the surface of a flaring disk heated by stellar and interstellar radiation is an appropriate description of the observations of TW Hya. LkCa 15 seems to be cooler at the surface, perhaps due to dust settling.

We present ISO-SWS observations of H2 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. We detect `warm' (T ~ 100-200 K) H2 gas around many sources, including tentatively the debris-disk objects. The mass of this warm gas ranges from ~1E-4 Solar mass up to 8E-3 Solar mass, and can constitute a non-negligible fraction of the total disk mass. Complementary single-dish 12CO 3-2, 13CO 3-2 and 12CO 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 12CO 3-2/ 13CO 3-2 integrated fluxes indicate that 12CO 3-2 is optically thick but that 13CO 3-2 is optically thin or at most moderately thick. The 13CO 3-2 lines have been used to estimate the cold gas mass. If a H2/CO conversion factor of 1E4 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. The warm gas is typically 1-10 % of the total mass deduced from millimeter continuum emission, but 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 H2, CO or dust masses is found for stars with ages in the range of 1E6-1E7 years, although a decrease is found for the older debris-disk star beta Pictoris. Existing models fail to explain the amount of warm gas quantitatively.

In circumstellar disks or shells it is often assumed that gas and dust temperatures are equal where the latter is determined by radiative equilibrium. This paper deals with the question whether this assumption is applicable for tenous circumstellar disks around young A stars. In this paper the thin hydrostatic equilibrium models described by Kamp & Bertoldi (2000) are combined with a detailed heating/cooling balance for the gas. The most important heating and cooling processes are heating through infrared pumping, heating due to the drift velocity of dust grains, and fine structure and molecular line cooling. Throughout the whole disk gas and dust are not efficiently coupled by collisions and hence their temperatures are quite different. Most of the gas in the disk models considered here stays well below 300 K. In the temperature range below 300 K the gas chemistry is not much affected by T_gas and therefore the simplifying approximation T_gas = T_dust can be used for calculating the chemical structure of the disk. Nevertheless the gas temperature is important for the quantitative interpretation of observations, like fine structure and molecular lines.

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.