The innermost astronomical unit in protoplanetary disks is a key region for stellar and planet formation, as exoplanet searches have shown a large occurrence of close-in planets that are located within the first au around their host star. We aim to reveal the morphology of the disk inner rim using near-infrared interferometric observations with milli-arcsecond resolution provided by infrared interferometry. We provide reconstructed images of 15 objects selected from the Herbig AeBe survey carried out with PIONIER at the VLTI, using SPARCO. We find that 40% of the systems are centrosymmetric at the angular resolution of the observations. For the rest of the objects, we find evidence for asymmetric emission due to moderate-to-strong inclination of a disk-like structure for 30% of the objects and noncentrosymmetric morphology due to a nonaxisymmetric and possibly variable environment (30%). Among the systems with a disk-like structure, 20% show a resolved dust-free cavity. The image reconstruction process is a powerful tool to reveal complex disk inner rim morphologies. At the angular resolution reached by near-infrared interferometric observations, most of the images are compatible with a centrally peaked emission (no cavity). For the most resolved targets, image reconstruction reveals morphologies that cannot be reproduced by generic parametric models. Moreover, the nonaxisymmetric disks show that the spatial resolution probed by optical interferometers makes the observations of the near-infrared emission sensitive to temporal evolution with a time-scale down to a few weeks. The evidence of nonaxisymmetric emission that cannot be explained by simple inclination and radiative transfer effects requires alternative explanations, such as a warping of the inner disks. Interferometric observations can, therefore, be used to follow the evolution of the asymmetry of those disks at a sub-au scale.

Near-IR observations of protoplanetary disks provide information about the properties of the inner disk. High resolution spectra of abundant molecules such as CO can be used to determine the disk structure in the warm inner parts. The $v2/v1$ ro-vibrational ratio of $v_{1-0}$ and $v_{2-1}$ transitions has been recently observed to follow distinct trends with the CO emitting radius, in a sample of TTauri and Herbig disks; these trends have been empirically interpreted as due to inner disk depletion from gas and dust. In this work we use existing thermo-chemical disk models to explore the interpretation of these observed trends in ro-vibrational CO emission. We use the radiation thermo-chemical code ProDiMo, exploring a set of previously published models with different disk properties and varying one parameter at a time: the inner radius, the dust-to-gas mass ratio, the gas mass. In addition, we use models where we change the surface density power law index, and employ a larger set of CO ro-vibrational levels, including also fluorescence from the first electronic state. We investigate these models for both TTauri and Herbig star disks. Finally, we include a set of DIANA models for individual TTauri and Herbig disks which were constructed to reproduce a large set of multi-wavelength observations.

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

Context. Several new multi-object spectrographs are currently planned or under construction that are capable of observing thousands of Galactic and extragalactic objects simultaneously. Aims. In this paper we present a probabilistic fibre-to-target assignment algorithm that takes spectrograph targeting constraints into account and is capable of dealing with multiple concurrent surveys. We present this algorithm using the 4-metre Multi-Object Spectroscopic Telescope (4MOST) as an example. Methods. The key idea of the proposed algorithm is to assign probabilities to fibre-target pairs. The assignment of probabilities takes the fibre positioner's capabilities and constraints into account. Additionally, these probabilities include requirements from surveys and take the required exposure time, number density variation, and angular clustering of targets across each survey into account. The main advantage of a probabilistic approach is that it allows for accurate and easy computation of the target selection function for the different surveys, which involves determining the probability of observing a target, given an input catalogue. Results. The probabilistic fibre-to-target assignment allows us to achieve maximally uniform completeness within a single field of view. The proposed algorithm maximises the fraction of successfully observed targets whilst minimising the selection bias as a function of exposure time. In the case of several concurrent surveys, the algorithm maximally satisfies the scientific requirements of each survey and no specific survey is penalised or prioritised. Conclusions. The algorithm presented is a proposed solution for the 4MOST project that allows for an unbiased targeting of many simultaneous surveys. With some modifications, the algorithm may also be applied to other multi-object spectroscopic surveys.

Mid-infrared molecular line emission detected with the Spitzer Space Telescope is often interpreted using slab models. However, we need to understand the mid-infrared line emission in 2D disk models, such that we gain information about from where the lines are being emitted and under which conditions, such that we gain information about number densities, temperatures, and optical depths in both the radial and vertical directions. In this paper, we introduce a series of 2D thermochemical models of a prototypical T Tauri protoplanetary disk, in order to examine how sensitive the line-emitting regions are to changes in the UV and X-ray fluxes, the disk flaring angle, dust settling, and the dust-to-gas ratio. These all affect the heating of the inner disk, and thus can affect the mid-infrared spectral lines. Using the ProDiMo and FLiTs codes, we produce a series of 2D thermochemical disk models. We find that there is often a significant difference between the gas and dust temperatures in the line emitting regions, and we illustrate that the size of the line emitting regions is relatively robust against changes in the stellar and disk parameters (namely, the UV and X-ray fluxes, the flaring angle, and dust settling). These results demonstrate the potential for localized variations in the line-emitting region to greatly affect the resulting spectra and line fluxes, and the necessity of allowing for such variations in our models.

This paper investigates the impact of the increased reaction rate constant due to tunneling effects on planet-forming disks. Our aim is to quantify the astrophysical implications of atom tunneling for simple molecules that are frequently used to infer disk structure information or to define the initial conditions for planet (atmosphere) formation. We explain the tunneling effect on reaction rates by using a scholarly example in comparison to previous UMIST2012 rate constants. In a chemical network with 1299 reactions, we identify all reactions that could show atom tunneling. We devise a simple formulation of reaction rate constants that overestimates tunneling and screen a standard T Tauri disk model for changes in species abundances. For those reactions found to be relevant, we find values of the most recent literature for the rate constants including tunneling and compare the resulting disk chemistry to the standard disk models. The rate constants in the UMIST2012 database in many cases already capture tunneling effects implicitly. A rigorous screening procedure identified three neutral-neutral reactions where atom tunneling could change simple molecule abundances. However, by adopting recent values of the rate constants of these reactions and due to the layered structure of planet-forming disks, the effects are limited to a small region where cold neutral-neutral chemistry dominates. Abundances of water close to the midplane snowline can increase by a factor of two at most compared to previous results with UMIST2012 rates. Observables from the disk surface, such as high excitation (> 500 K) water line fluxes, decrease by 60% at most when tunneling effects are explicitly excluded. On the other hand, disk midplane quantities relevant for planet formation such as the C-to-O ratio and also the ice-to-rock ratio are clearly affected by these gas-phase tunneling effects.

The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium.

In this paper, we couple together the dust evolution code two-pop-py with the thermochemical disk modelling code ProDiMo. We create a series of thermochemical disk models that simulate the evolution of dust over time from 0.018 Myr to 10 Myr, including the radial drift, growth, and settling of dust grains. We examine the effects of this dust evolution on the mid-infrared gas emission, focussing on the mid-infrared spectral lines of C2H2, CO2, HCN, NH3, OH, and H2O that are readily observable with Spitzer and the upcoming E-ELT and JWST. The addition of dust evolution acts to increase line fluxes by reducing the population of small dust grains. We find that the spectral lines of all species except C2H2 respond strongly to dust evolution, with line fluxes increasing by more than an order of magnitude across the model series as the density of small dust grains decreases over time. The C2H2 line fluxes are extremely low due to a lack of abundance in the infrared line-emitting regions, despite C2H2 being commonly detected with Spitzer, suggesting that warm chemistry in the inner disk may need further investigation. Finally, we find that the CO2 flux densities increase more rapidly than the other species as the dust disk evolves. This suggests that the flux ratios of CO2 to other species may be lower in disks with less-evolved dust populations.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

Several detections of wide-orbit planet-mass/sub-stellar companions around young solar-like stars were reported in the last decade. The origin of those possible planets is still unclear but accretion tracers and VLT/SPHERE observations indicate that they are surrounded by circumplanetary material or even a circumplanetary disk. We want to investigate if the gas component of disks around wide-orbit companions is detectable with current and future (sub)mm telescopes and what constraints such gas observations can provide on the nature of the circumplanetary material and on the mass of the companion. We applied the radiation thermo-chemical disk code ProDiMo to model the dust and gas component of passive circumplanetary disks and produced realistic synthetic observables. We considered different companion properties, disk parameters and radiative environments and compared the resulting synthetic observables to telescope sensitivities and to existing dust observations. The main criterion for a successful detection is the size of the circumplanetary disk. At a distance of about 150 pc, a circumplanetary disk with an outer radius of about 10 au is detectable with ALMA in about 6 hours in optically thick CO lines. Other aspects such as the companion's luminosity, disk inclination and background radiation fields are also relevant and should be considered to optimize the observing strategy for detection experiments. For most of the known wide-orbit planet-mass companions, their maximum theoretical disk size of one third of the Hill radius would be sufficient to allow detection of CO lines. It is therefore feasible to detect their gas disks and constrain the mass of the companion through the kinematic signature. Even in the case of non-detections such observations will provide stringent constraints on disk size and gas mass, information crucial for formation theories.

Consistent modeling of protoplanetary disks requires the simultaneous solution of both continuum and line radiative transfer, heating/cooling balance between dust and gas and, of course, chemistry. Such models depend on panchromatic observations that can provide a complete description of the physical and chemical properties and energy balance of protoplanetary systems. Along these lines we present a homogeneous, panchromatic collection of data on a sample of 85 T Tauri and Herbig Ae objects for which data cover a range from X-rays to centimeter wavelengths. Datasets consist of photometric measurements, spectra, along with results from the data analysis such as line fluxes from atomic and molecular transitions. Additional properties resulting from modeling of the sources such as disc mass and shape parameters. dust size and PAH properties are also provided for completeness. Targets were selected based on their properties data availability. Data from more than 50 different telescopes and facilities were retrieved and combined in homogeneous datasets directly from public data archives or after being extracted from more than 100 published articles. X-ray data for a subset of 56 sources represent an exception as they were reduced from scratch and are presented here for the first time. Compiled datasets along with a subset of continuum and emission-line models are stored in a dedicated database and distributed through a publicly accessible online system. All datasets contain metadata descriptors that allow to backtrack them to their original resources. The graphical user interface of the online system allows the user to visually inspect individual objects but also compare between datasets and models. It also offers to the user the possibility to download any of the stored data and metadata for further processing.

Anomalies in the abundance measurements of short lived radionuclides in meteorites indicate that the protosolar nebulae was irradiated by a large number of energetic particles ($E\gtrsim10\,$MeV), often called solar cosmic rays. The particle flux of the contemporary Sun cannot explain these anomalies, but, similar to \mbox{T Tauri} stars, the young Sun was more active and probably produced enough high energy particles. However, the stellar particle (SP) flux of young stars is essentially unknown. We model the impact of high-energy ionization sources on the chemistry of the circumstellar environment (disks and envelopes). The model includes X-ray radiative transfer and makes use of particle transport models to calculate the individual molecular hydrogen ionization rates. We study the impact on the chemistry via the ionization tracers HCO$^+$ and N$_2$H$^+$. We argue that spatially resolved observations of those molecules combined with detailed models allow for disentangling the contribution of the individual high-energy ionization sources and to put constraints on the SP flux in young stars.

Molecular hydrogen (H2) is the main constituent of the gas in the planet-forming disks that surround many PMS stars. H2 can be incorporated in the atmosphere of the giant planets. HD has been detected in a few disks and can be considered the most reliable tracer of H2. We wish to form H2 and HD efficiently for the varied conditions encountered in protoplanetary disks: the densities vary from 1E4 to 1E16 cm^-3; the dust temperatures range from 5 to 1500 K, the gas temperatures go from 5 to a few 1000 Kelvin, and the ultraviolet field can be 1E7 stronger than the standard interstellar field. We implemented a comprehensive model of H2 and HD formation on cold and warm grain surfaces and via hydrogenated PAHs in the physico-chemical code ProDiMo. The H2 and HD formation can proceed via the Langmuir-Hinshelwood and Eley-Ridel mechanisms for physisorbed or chemisorbed H (D) atoms. H2 and HD also form by H (D) abstraction from hydrogenated neutral and ionised PAHs and via gas phase reactions. H2 and HD are formed efficiently on dust grain surfaces from 10 to 700 K. All the deuterium is converted into HD in UV shielded regions as soon as H2 is formed by gas-phase D abstraction reactions. The detailed model compares well with standard analytical prescriptions for H2 (HD) formation. At low temperatures, H2 is formed from the encounter of two physisorbed atoms. HD molecules form on the grain surfaces and in the gas-phase. At temperatures greater than 20 K, the meeting between a weakly bound H- (or D-) atom or a gas-phase H (D) atom and a chemisorbed atom is the most efficient H2 formation route. H2 formation through hydrogenated PAHs alone is efficient above 80 K. The contribution of hydrogenated PAHs to the overall H2 and HD formation is relatively low if chemisorption on silicate is taken into account and if a small hydrogen abstraction cross-section is used.

The origin of the reservoirs of water on Earth is debated. The Earth's crust may contain at least three times more water than the oceans. This crust water is found in the form of phyllosilicates, whose origin probably differs from that of the oceans. We test the possibility to form phyllosilicates in protoplanetary disks, which can be the building blocks of terrestrial planets. We developed an exploratory rate-based warm surface chemistry model where water from the gas-phase can chemisorb on dust grain surfaces and subsequently diffuse into the silicate cores. We apply the phyllosilicate formation model to a zero-dimensional chemical model and to a 2D protoplanetary disk model (ProDiMo). The disk model includes in addition to the cold and warm surface chemistry continuum and line radiative transfer, photoprocesses (photodissociation, photoionization, and photodesorption), gas-phase cold and warm chemistry including three-body reactions, and detailed thermal balance. Despite the high energy barrier for water chemisorption on silicate grain surfaces and for diffusion into the core, the chemisorption sites at the surfaces can be occupied by a hydroxyl bond (-OH) at all gas and dust temperatures from 80 to 700 K for a gas density of 2E4 cm^-3. The chemisorption sites in the silicate cores are occupied at temperatures between 250 and 700 K. At higher temperatures thermal desorption of chemisorbed water occurs. The occupation efficiency is only limited by the maximum water uptake of the silicate. The timescales for complete hydration are at most 1E5 years for 1 mm radius grains at a gas density of 1E8 cm^-3. Phyllosilicates can be formed on dust grains at the dust coagulation stage in protoplanetary disks within 1 Myr. It is however not clear whether the amount of phyllosilicate formed by warm surface chemistry is sufficient compared to that found in Solar System objects.

The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps. We define and apply a standardized modelling procedure to simultaneously fit all these data by state-of-the-art modelling codes (ProDiMo, MCFOST, MCMax) which solve the continuum and line radiative transfer, disk chemistry, and the heating & cooling balance for both the gas and the dust. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh-Jeans limit. Some line observations cannot be reproduced by the models, probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available at an online database (http://www.univie.ac.at/diana).

Context. Disks around pre-main-sequence stars evolve over time by turbulent viscous spreading. The main contender to explain the strength of the turbulence is the Magneto-Rotational-Instability (MRI) model, whose efficiency depends on the disk ionization fraction. Aims. We aim at computing self-consistently the chemistry including PAH charge chemistry, the grain charging and an estimate of an effective value of the turbulence alpha parameter in order to find observational signatures of disk turbulence. Methods. We introduced PAH and grain charging physics and their interplay with other gas-phase reactions in the physico-chemical code ProDiMo. Non-ideal magnetohydrodynamics effects such as Ohmic and ambipolar diffusion are parametrized to derive an effective value for the turbulent parameter alpha_eff . We explored the effects of turbulence heating and line broadening on CO isotopologue sub-millimeter lines. Results. The spatial distribution of alpha_eff depends on various unconstrained disk parameters such as the magnetic parameter beta_mag or the cosmic ray density distribution inside the protoplanetary disks. The inner disk midplane shows the presence of the so-called dead-zone where the turbulence is quasi-inexistent. The disk is heated mostly by thermal accommodation on dust grains in the dead-zone, by viscous heating outside the dead-zone up to a few hundred astronomical units, and by chemical heating in the outer disk. The CO rotational lines probe the warm molecular disk layers where the turbulence is at its maximum. However, the effect of turbulence on the CO line profiles is minimal and difficult to distinguish from the thermal broadening. Conclusions. Viscous heating of the gas in the disk midplane outside the dead-zone is efficient. The determination of alpha from CO rotational line observations alone is challenging.

We present a study of the CO isotopologues and the high-density tracers H$_{2}$CO, HCO$^{+}$, and CS in Class 0/I proto-brown dwarfs (proto-BDs). We have used the IRAM 30m telescope to observe the $^{12}$CO (2-1), $^{13}$CO (2-1), C$^{18}$O (2-1), C$^{17}$O (2-1), H$_{2}$CO (3-2), HCO$^{+}$ (3-2), and CS (5-4) lines in 7 proto-BDs. The hydrogen column density for the proto-BDs derived from the CO gas emission is $\sim$2-15 times lower than that derived from the dust continuum emission, indicating CO depletion from the gas-phase. The mean H$_{2}$CO ortho-to-para ratio is $\sim$3 for the proto-BDs and indicates gas-phase formation for H$_{2}$CO. We have investigated the correlations in the molecular abundances between the proto-BDs and protostars. Proto-BDs on average show a factor of $\sim$2 higher ortho-to-para H$_{2}$CO ratio than the protostars. Possible explanations include a difference in the H$_{2}$CO formation mechanism, spin-selective photo-dissociation, self-shielding effects, or different emitting regions for the ortho and para species. There is a tentative trend of a decline in the HCO$^{+}$ and H$_{2}$CO abundances with decreasing bolometric luminosity, while the CS and CO abundances show no particular difference between the proto-BDs and protostars. These trends reflect the scaled-down physical structures for the proto-BDs compared to protostars and differences in the peak emitting regions for these species. The C$^{17}$O isotopologue is detected in all of the proto-BDs as well as the more evolved Class Flat/Class II BDs in our sample, and can probe the quiescent gas at both early and late evolutionary stages.

We present results from a study of nitrogen chemistry in Class 0/I proto-brown dwarfs (proto-BDs). We have used the IRAM 30 m telescope to observe the CN (2-1), HCN (3-2), and HNC (3-2) lines in 7 proto-BDs. All proto-BDs show a large CN/HCN abundance ratio of >20, and a HNC/HCN abundance ratio close to or larger than unity. The enhanced CN/HCN ratios can be explained by high UV flux originating from an active accretion zone in the proto-BDs. The larger than unity HNC/HCN ratio for the proto-BDs is likely caused by a combination of low temperature and high density. Both CN and HNC show a flat distribution with CO, indicating that these species can survive in regions where CO is depleted. We have investigated the correlations in the molecular abundances of these species for the proto-BDs with Class 0/I protostars. We find tentative trends of CN (HCN) abundances being about an order of magnitude higher (lower) in the proto-BDs compared to protostars. HNC for the proto-BDs shows a nearly constant abundance unlike the large spread of ~2 orders of magnitude seen for the protostars. Also notable is a rise in the HNC/HCN abundance ratio for the lowest luminosity objects, suggesting that this ratio is higher under low-temperature environments. None of the relatively evolved Class Flat/Class II brown dwarfs in our sample show emission in HNC. The HNC molecule can be considered as an efficient tracer to search and identify early stage sub-stellar mass objects.

Endogenous or exogenous, dry or wet, various scenarios have been so far depicted for the origin of water on our Solar Systems rocky bodies. Hydrated silicates found in meteorites and in interplanetary dust particles together with observations of abundant water reservoirs in the habitable zone of protoplanetary disks are evidences that support aqueous alteration of silicate dust grains by water vapor condensation in a nebular setting. We investigate the thermodynamics (temperature and pressure dependencies) and kinetics (adsorption rates and energies, surface diffusion and cluster formation) of water adsorption on surfaces of forsterite grains. (abbreviated) The protoplanetary disk model (ProDiMo) code is tuned to simulate the thermochemical disk structure of the early solar nebula at three evolutionary stages. Pressure, temperature and water vapor abundance within 1 au from the protosun were extracted and used as input for a Monte Carlo code to model water associative adsorption using adsorption energies that resemble the forsterite [100] crystal lattice. Hydration of forsterite surfaces by water vapor adsorption could have occurred within the nebula lifetime already at a density of 1e8 cm-3, with increasing surface coverage for higher water vapor densities. Full surface coverage is attained for temperatures lower than 500 K, while for hotter grain surfaces water cluster formation plays a crucial role. (abbreviated) This work shows that water cluster formation enhances the water surface coverage and enables a stable water layer to form at high temperature and low water vapor density conditions. Finally, surface diffusion of physisorbed water molecules shortens the timescale for reaching steady state, enabling phyllosilicate formation within the solar nebula timescale.

We introduce a new modelling framework called FLiTs to simulate infrared line emission spectra from protoplanetary discs. This paper focuses on the mid-IR spectral region between 9.7 um to 40 um for T Tauri stars. The generated spectra contain several tens of thousands of molecular emission lines of H2O, OH, CO, CO2, HCN, C2H2, H2 and a few other molecules, as well as the forbidden atomic emission lines of SI, SII, SIII, SiII, FeII, NeII, NeIII, ArII and ArIII. In contrast to previously published works, we do not treat the abundances of the molecules nor the temperature in the disc as free parameters, but use the complex results of detailed 2D ProDiMo disc models concerning gas and dust temperature structure, and molecular concentrations. FLiTs computes the line emission spectra by ray tracing in an efficient, fast and reliable way. The results are broadly consistent with R=600 Spitzer/IRS observational data of T Tauri stars concerning line strengths, colour, and line ratios. In order to achieve that agreement, however, we need to assume either a high gas/dust mass ratio of order 1000, or the presence of illuminated disc walls at distances of a few au. The molecules in these walls cannot be photo-dissociated easily by UV because of the large densities in the walls which favour their re-formation. Most observable molecular emission lines are found to be optically thick, rendering a standard analysis with column densities difficult. We find that the difference between gas and dust temperatures in the disc surface is important for the line formation. We briefly discuss the effects of C/O ratio and choice of chemical rate network on these results. Our analysis offers new ways to infer the chemical and temperature structure of T Tauri discs from future JWST/MIRI observations, and to possibly detect secondary illuminated disc walls based on their specific mid-IR molecular signature.

The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the chemical structure of disks. Further, we want to investigate the importance of the dust for X-ray radiative transfer in disks. We present a new X-ray radiative transfer module for the radiation thermo-chemical disk code ProDiMo, which includes X-ray scattering and absorption by both the gas and dust component. For the X-ray radiative transfer, we consider irradiation by the star and by X-ray background fields. To study the impact of X-rays on the chemical structure of disks we use the well-established disk ionization tracers N2H+ and HCO+. For evolved dust populations, X-ray opacities are mostly dominated by the gas; only for photon energies $E\gtrsim5-10\,$keV, dust opacities become relevant. Consequently, the local disk X-ray radiation field is only affected in dense regions close to the disk midplane. X-ray background fields can dominate the local X-ray disk ionization rate for disk radii $r\gtrsim20\,$au. However, the N2H+ and HCO+ column densities are only significantly affected in case of low cosmic-ray ionization rates, or if the background flux is at least a factor of ten higher than the flux level expected for clusters typical for the solar vicinity. Observable signatures of X-ray background fields in low-mass star-formation regions, like Taurus, are only expected for cluster members experiencing a strong X-ray background field. For the majority of the cluster members, the X-ray background field has only little impact on the disk chemical structure.

Accurate measurements of the physical structure of protoplanetary discs are critical inputs for planet formation models. These constraints are traditionally established via complex modelling of continuum and line observations. Instead, we present an empirical framework to locate the CO isotopologue emitting surfaces from high spectral and spatial resolution ALMA observations. We apply this framework to the disc surrounding IM Lupi, where we report the first direct, i.e. model independent, measurements of the radial and vertical gradients of temperature and velocity in a protoplanetary disc. The measured disc structure is consistent with an irradiated self-similar disc structure, where the temperature increases and the velocity decreases towards the disc surface. We also directly map the vertical CO snow line, which is located at about one gas scale height at radii between 150 and 300 au, with a CO freeze-out temperature of $21\pm2$ K. In the outer disc ($> 300$ au), where the gas surface density transitions from a power law to an exponential taper, the velocity rotation field becomes significantly sub-Keplerian, in agreement with the expected steeper pressure gradient. The sub-Keplerian velocities should result in a very efficient inward migration of large dust grains, explaining the lack of millimetre continuum emission outside of 300 au. The sub-Keplerian motions may also be the signature of the base of an externally irradiated photo-evaporative wind. In the same outer region, the measured CO temperature above the snow line decreases to $\approx$ 15 K because of the reduced gas density, which can result in a lower CO freeze-out temperature, photo-desorption, or deviations from local thermodynamic equilibrium.

To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Our observations in the near-infrared K-band continuum reveal a disk around each binary component, with similar half-flux radii of about 0.1 au at d$\sim$130 pc, inclinations ($i=$28$\pm$3$^o$\ and $i=$22$\pm$6$^o$), and position angles (PA=0$^o\pm$6$^o$ and PA=-2$^o\pm$12$^o$), suggesting that they formed from the fragmentation of a common disk. The S CrA N spectrum shows bright HeI and Br$\gamma$ line emission exhibiting inverse P-Cygni profiles, typically associated with infalling gas. The continuum-compensated Br$\gamma$ line visibilities of S CrA N show the presence of a compact Br$\gamma$ emitting region the radius of which is about $\sim$0.06 au, which is twice as big as the truncation radius. This component is mostly tracing a wind. Moreover, a slight radius change between the blue- and red-shifted Br$\gamma$ line components is marginally detected. The presence of an inverse P-Cygni profile in the HeI and Br$\gamma$ lines, along with the tentative detection of a slightly larger size of the blue-shifted Br$\gamma$ line component, hint at the simultaneous presence of a wind and magnetospheric accretion in S CrA N.

Aims. We define a small and large chemical network which can be used for the quantitative simultaneous analysis of molecular emission from the near-IR to the submm. We revise reactions of excited molecular hydrogen, which are not included in UMIST, to provide a homogeneous database for future applications. Methods. We use the thermo-chemical disk modeling code ProDiMo and a standard T Tauri disk model to evaluate the impact of various chemical networks, reaction rate databases and sets of adsorption energies on a large sample of chemical species and emerging line fluxes from the near-IR to the submm wavelength range. Results. We find large differences in the masses and radial distribution of ice reservoirs when considering freeze-out on bare or polar ice coated grains. Most strongly the ammonia ice mass and the location of the snow line (water) change. As a consequence molecules associated to the ice lines such as N2H+ change their emitting region; none of the line fluxes in the sample considered here changes by more than 25% except CO isotopologues, CN and N2H+ lines. The three-body reaction N+H2+M plays a key role in the formation of water in the outer disk. Besides that, differences between the UMIST 2006 and 2012 database change line fluxes in the sample considered here by less than a factor 2 (a subset of low excitation CO and fine structure lines stays even within 25%); exceptions are OH, CN, HCN, HCO+ and N2H+ lines. However, different networks such as OSU and KIDA 2011 lead to pronounced differences in the chemistry inside 100 au and thus affect emission lines from high excitation CO, OH and CN lines. H2 is easily excited at the disk surface and state-to-state reactions enhance the abundance of CH+ and to a lesser extent HCO+. For sub-mm lines of HCN, N2H+ and HCO+, a more complex larger network is recommended. ABBREVIATED

Episodic accretion is an important process in the evolution of young stars and their environment. The observed strong luminosity bursts of young stellar objects likely have a long lasting impact on the chemical evolution of the disk and envelope structure. We want to investigate observational signatures of the chemical evolution in the post-burst phase for embedded sources. With such signatures it is possible to identify targets that experienced a recent luminosity burst. We present a new model for episodic accretion chemistry based on the 2D, radiation thermo-chemical disk code ProDiMo. We have extended ProDiMo with a proper treatment for envelope structures. For a representative Class I model, we calculated the chemical abundances in the post-burst phase and produced synthetic observables like intensity maps and radial profiles. During a burst many chemical species, like CO, sublimate from the dust surfaces. As the burst ends they freeze out again (post-burst phase). This freeze-out happens from inside-out due to the radial density gradient in the disk and envelope structure. This inside-out freeze-out produces clear observational signatures in spectral line emission, like rings and distinct features in the slope of radial intensity profiles. We fitted synthetic C18O J=2-1 observations with single and two component fits and find that post-burst images are much better matched by the latter. Comparing the quality of such fits allows identification of post-burst targets in a model-independent way. Our models confirm that it is possible to identify post-burst objects from spatially resolved CO observations. However, to derive proper statistics, like frequencies of bursts, from observations it is important to consider aspects like the inclination and structure of the target and also dust properties as those have a significant impact on the freeze-out timescale.

Anomalies in the abundance measurements of short lived radionuclides in meteorites indicate that the protosolar nebulae was irradiated by a high amount of energetic particles (E$\gtrsim$10 MeV). The particle flux of the contemporary Sun cannot explain these anomalies. However, similar to T Tauri stars the young Sun was more active and probably produced enough high energy particles to explain those anomalies. We want to study the interaction of stellar energetic particles with the gas component of the disk and identify possible observational tracers of this interaction. We use a 2D radiation thermo-chemical protoplanetary disk code to model a disk representative for T Tauri stars. We use a particle energy distribution derived from solar flare observations and an enhanced stellar particle flux proposed for T Tauri stars. For this particle spectrum we calculate the stellar particle ionization rate throughout the disk with an accurate particle transport model. We study the impact of stellar particles for models with varying X-ray and cosmic-ray ionization rates. We find that stellar particle ionization has a significant impact on the abundances of the common disk ionization tracers HCO$^+$ and N$_2$H$^+$, especially in models with low cosmic-ray ionization rates. In contrast to cosmic rays and X-rays, stellar particles cannot reach the midplane of the disk. Therefore molecular ions residing in the disk surface layers are more affected by stellar particle ionization than molecular ions tracing the cold layers/midplane of the disk. Spatially resolved observations of molecular ions tracing different vertical layers of the disk allow to disentangle the contribution of stellar particle ionization from other competing ionization sources. Modeling such observations with a model like the one presented here allows to constrain the stellar particle flux in disks around T Tauri stars.

The physical properties of brown dwarf discs, in terms of their shapes and sizes, are still largely unexplored by observations. To what extent brown dwarf discs are similar to scaled-down T Tauri discs is currently unknown, and this work is a step towards establishing a relationship through the eventual modelling of future observations. We use observations of the brown dwarf disc $\rho$ Oph 102 to infer a fiducial model around which we build a small grid of brown dwarf disc models, in order to model the CO, HCN, and HCO+ line fluxes and the chemistry which drives their abundances. These are the first brown dwarf models to be published which relate detailed, 2D radiation thermochemical disc models to observational data. We predict that moderately extended ALMA antenna configurations will spatially resolve CO line emission around brown dwarf discs, and that HCN and HCO+ will be detectable in integrated flux, following our conclusion that the flux ratios of these molecules to CO emission are comparable to that of T Tauri discs. These molecules have not yet been observed in sub-mm wavelengths in a brown dwarf disc, yet they are crucial tracers of the warm surface-layer gas and of ionization in the outer parts of the disc. We present the prediction that if the physical and chemical processes in brown dwarf discs are similar to those that occur in T Tauri discs -- as our models suggest -- then the same diagnostics that are used for T Tauri discs can be used for brown dwarf discs (such as HCN and HCO+ lines that have not yet been observed in the sub-mm), and that these lines should be observable with ALMA. Through future observations, either confirmation (or refutation) of these ideas about brown dwarf disc chemistry will have strong implications for our understanding of disc chemistry, structure, and subsequent planet formation in brown dwarf discs.

Context: Quantifying the gas content inside the dust gaps of transition disks is important to establish their origin. Aims: We seek to constrain the surface density of warm gas in the disk of HD 139614, a Herbig Ae star with a transition disk exhibiting a dust gap from 2.3 to 6 AU. Methods: We have obtained ESO/VLT CRIRES high-resolution spectra of CO ro-vibrational emission. We derived constraints on the disk's structure by modeling the line-profiles, the spectroastrometric signal, and the rotational diagrams using flat Keplerian disk models. Results: We detected v=1-0 12CO, 2-1 12CO, 1-0 13CO, 1-0 C18O, and 1-0 C17O ro-vibrational lines. 12CO v=1-0 lines have an average width of 14 km/s, Tgas of 450 K and an emitting region from 1 to 15 AU. 13CO and C18O lines are on average 70 and 100 K colder, 1 and 4 km/s narrower, and are dominated by emission at R>6 AU. The 12CO v=1-0 line-profile indicates that if there is a gap in the gas it must be narrower than 2 AU. We find that a drop in the gas surface density (delta_gas) at R<5-6 AU is required to simultaneously reproduce the line-profiles and rotational diagrams of the three CO isotopologs. Delta_gas can range from 10^-2 to 10^-4 depending on the gas-to-dust ratio of the outer disk. We find that at 1<R<6 AU the gas surface density profile is flat or increases with radius. We derive a gas column density at 1<R<6 AU of NH=3x10^19 - 10^21 cm^-2. We find a 5sigma upper limit on NCO at R<1 AU of 5x10^15 cm^-2 (NH<5x10^19 cm^-2). Conclusions: The dust gap in the disk of HD 139614 has gas. The gas surface density in the disk at R<6 AU is significantly lower than the surface density expected from HD 139614's accretion rate assuming a viscous alpha-disk model. The gas density drop, the non-negative density gradient of the gas inside 6 AU, and the absence of a wide (>2 AU) gas gap suggest the presence of an embedded <2 MJ planet at around 4 AU.

Context. It is now generally accepted that the near-infrared excess of Herbig AeBe stars originates in the dust of a circumstellar disk. Aims. The aims of this article are to infer the radial and vertical structure of these disks at scales of order one au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6micron) PIONIER/VLTI interferometer. The largest baselines allowed us to resolve (at least partially) structures of a few tenths of an au at typical distances of a few hundred parsecs. Dedicated UBVRIJHK photometric measurements were also obtained. Spectral and 2D geometrical parameters are extracted via fits of a few simple models: ellipsoids and broadened rings with azimuthal modulation. Model bias is mitigated by parallel fits of physical disk models. Sample statistics were evaluated against similar statistics for the physical disk models to infer properties of the sample objects as a group. Results. We find that dust at the inner rim of the disk has a sublimation temperature Tsub~1800K. A ring morphology is confirmed for approximately half the resolved objects; these rings are wide delta_r>=0.5. A wide ring favors a rim that, on the star-facing side, looks more like a knife edge than a doughnut. The data are also compatible with a the combination of a narrow ring and an inner disk of unspecified nature inside the dust sublimation radius. The disk inner part has a thickness z/r~0.2, flaring to z/r~0.5 in the outer part. We confirm the known luminosity-radius relation; a simple physical model is consistent with both the mean luminosity-radius relation and the ring relative width; however, a significant spread around the mean relation is present. In some of the objects we find a halo component, fully resolved at the shortest interferometer spacing, that is related to the HAeBe class.

Mid-IR water lines from protoplanetary disks around T Tauri stars have a detection rate of 50\%. Models have identified multiple physical properties of disks such as dust-to-gas mass ratio, dust size power law distribution, disk gas mass, disk inner radius, and disk scale height as potential explanation for the current detection rate. We search for a connection between mid-IR water line fluxes and the strength of the 10~$\mu$m silicate feature. We analyse observed water line fluxes from three blends and compute the 10~$\mu$m silicate feature strength from Spitzer spectra. We use a series of published models, exploring disk dust and gas properties, and the effects of different stars. The models also show that the increasing stellar luminosity enhance simultaneously the strength of this dust feature and the water lines fluxes. No correlation is found between the observed mid-IR water lines and the 10~$\mu$m silicate. Our sample shows the same difference in the peak strength between amorphous and crystalline silicates that was noted in earlier studies, but our models do not support this intrinsic difference in silicate peak strength. Individual properties of our model series are not able to reproduce the most extreme observations, suggesting that more complex dust properties are required. A parametrized settling prescription is able to boost the peak strength by a factor 2 for the standard model. Water line fluxes are unrelated to the composition of the dust. The pronounced regular trends seen in the model results are washed out in the data due to the larger diversity in stellar and disk properties compared to our model series. The disks with with weaker mid-IR water line fluxes are depleted in gas or enhanced in dust in the inner 10~au. In the case of gas depleted disks, settling produces very strong 10~$\mu$m silicate features, with strong peak strength.

Because of the very peculiar conditions of chemistry in many astrophysical gases (low densities, mostly low temperatures, kinetics-dominated chemical evolution), great efforts have been devoted to study molecular signatures and chemical evolution. While experiments are being performed in many laboratories, it appears that the efforts directed towards theoretical works are not as strong. This report deals with the present status of chemical physics/physical chemistry theory, for the qualitative and quantitative understanding of kinetics of molecular scattering, being it reactive or inelastic. By gathering several types of expertise, from applied mathematics to physical chemistry, dialog is made possible, as a step towards new and more adapted theoretical frameworks, capable of meeting the theoretical, methodological and numerical challenges of kinetics-dominated gas phase chemistry in astrophysical environments. A state of the art panorama is presented, alongside present-day strengths and shortcomings. However, coverage is not complete, being limited in this report to actual attendance of the workshop. Some paths towards relevant progress are proposed.

LkCa 15 is an extensively studied star in the Taurus region known for its pre-transitional disk with a large inner cavity in dust continuum and normal gas accretion rate. The most popular hypothesis to explain the LkCa 15 data invokes one or more planets to carve out the inner cavity, while gas continues to flow across the gap from the outer disk onto the central star. We present spatially unresolved HCO+ J=4-3 observations of the LkCa 15 disk from the JCMT and model the data with the ProDiMo code. We find that: (1) HCO+ line-wings are clearly detected, certifying the presence of gas in the cavity within <50 AU of the star. (2) Reproducing the observed line-wing flux requires both a significant suppression of cavity dust (by a factor >10^4 compared to the ISM) and a substantial increase in the gas scale-height within the cavity (H_0/R_0 ~ 0.6). An ISM dust-to-gas ratio (d:g=10^-2) yields too little line-wing flux regardless of the scale-height or cavity gas geometry, while a smaller scale-height also under predicts the flux even with a reduced d:g. (3) The cavity gas mass is consistent with the surface density profile of the outer disk extended inwards to the sublimation radius (corresponding to mass M_d ~ 0.03 M_sun), and masses lower by a factor >10 appear to be ruled out.

Aims. CO ro-vibrational lines observed from Herbig group II discs are seen to be often broad while the same lines observed from group I discs are often narrow. This difference is not well understood. In this paper we explore the underlying cause for this difference and provide a pathway for a better understanding of the geometry and structure of the inner discs around Herbig Ae/Be stars. Methods. High-spectral-resolution infrared spectra of CO ro-vibrational emission from six Herbig Ae/Be candidate stars were taken with CRIRES (CRyogenic high-resolution InfraRed Echelle Spectrograph) at the VLT (Very Large Telescope). From these spectra, we produce individual and co-added CO ro-vibrational line profiles. We investigate line profile shape differences, and we explore the FWHM (Full Width Half Maximum) variations with J quantum number in the context of disc geometry. Furthermore, we put our new sources into the context of earlier observed sources to study a large sample. For comparison, we also investigate the FWHM variations with J of modelled CO ro-vibrational lines from two typical disc geometries produced with the thermo-chemical disc modelling code ProDiMo. Results. - abbreviated - Comparing dust and gas inner disc geometries (inferred by SED and CO ro-vibrational emission) for the expanded sample of observed Herbig discs, we find no clear correspondence between the SED groups of the sources and their inner CO radius. Conclusions. The FWHM versus J is a potential new gas diagnostic for the inner disc with e.g. a constant FWHM versus J indicating the presence of a large gas hole or gap. Both models and observations indicate the potential of this new diagnostic. - abbreviated -

We present for the first time a direct comparison of multi-epoch (2001-2002 and 2012) CO ro-vibrational emission lines from HD163296. We find that both the line shapes and the FWHM (Full Width Half Maximum) differ between these two epochs. The FWHM of the median observed line profiles are 10-25 km/s larger in the earlier epoch, and confirmed double peaks are only present in high J lines from 2001-2002. The line wings of individual transitions are similar in the two epochs making an additional central component in the later epoch a likely explanation for the single peaks and the lower FWHM. Variations in NIR brightness have been reported and could be linked to the observed variations. Additionally, we use the thermo chemical disc code ProDiMo to compare for the first time the line shapes, peak separations, FWHM, and line fluxes, to those observed. The ProDiMo model reproduces the peak separations, and low and mid J line fluxes well. The FWHM however, are over predicted and high J line fluxes are under predicted. We propose that a variable non-Keplerian component of the CO ro-vibrational emission, such as a disc wind or an episodic accretion funnel, is causing the difference between the two data sets collected at different epochs, and between model and observations. Additional CO ro-vibrational line detections (with CRIRES/VLT or NIRSPEC/Keck) or [Ne II] line observations with VISIR/VLT could help to clarify the cause of the variability.

We present calculated rate coefficients for ro-vibrational transitions of CO in collisions with H atoms for a gas temperature range of 10 K $\leq T \leq$ 3000 K, based on the recent three-dimensional ab initio H-CO interaction potential of Song et al(2013). Rate coefficients for ro-vibrational $v=1,j=0-30 \rightarrow v'=0, j'$ transitions were obtained from scattering cross sections previously computed with the close-coupling method by Song et al(2015). Combining these with the rate coefficients for vibrational $v=1-5 \rightarrow v' < v$ quenching obtained with the infinite-order sudden approximation, we propose a new extrapolation scheme that yields the rate coefficients for ro-vibrational $v=2-5,j=0-30 \rightarrow v',j'$ de-excitation. Cross sections and rate coefficients for ro-vibrational $v=2, j=0-30 \rightarrow v'=1,j'$ transitions calculated with the close-coupling method confirm the effectiveness of this extrapolation scheme. Our calculated and extrapolated rates are very different from those that have been adopted in the modeling of many astrophysical environments. The current work provides the most comprehensive and accurate set of ro-vibrational de-excitation rate coefficients for the astrophysical modeling of the H-CO collision system. Application of the previously available and new data sets in astrophysical slab models shows that the line fluxes typically change by 20-70% in high temperature environments (800 K) with an H/H$_2$ ratio of 1; larger changes occur for lower temperatures.

Spatially resolving the inner dust cavity of the transitional disks is a key to understanding the connection between planetary formation and disk dispersal. The disk around the Herbig star HD 139614 is of particular interest since it presents a pretransitional nature with an au-sized gap, in the dust, that was spatially resolved by mid-IR interferometry. Using new NIR interferometric observations, we aim to characterize the 0.1-10~au region of the HD~139614 disk further and identify viable mechanisms for the inner disk clearing. We report the first multiwavelength radiative transfer modeling of the interferometric data acquired on HD~139614 with PIONIER, AMBER, and MIDI, complemented by Herschel/PACS photometries. We confirm a gap structure in the um-sized dust, extending from about 2.5 au to 6 au, and constrained the properties of the inner dust component: e.g., a radially increasing surface density profile, and a depletion of 10^3 relative to the outer disk. Since self-shadowing and photoevaporation appears unlikely to be responsible for the au-sized gap of HD~139614, we thus tested if dynamical clearing could be a viable mechanism using hydrodynamical simulations to predict the gaseous disk structure. Indeed, a narrow au-sized gap is expected when a single giant planet interacts with the disk. Assuming that small dust grains are well coupled to the gas, we found that a ~ 3~Mjup planet located at 4.5 au from the star could, in less than 1 Myr, reproduce most of the aspects of the dust surface density profile, while no significant depletion in gas occurred in the inner disk, in contrast to the dust. However, the dust-depleted inner disk could be explained by the expected dust filtration by the gap and the efficient dust growth/fragmentation in the inner disk regions. Our results support the hypothesis of a giant planet opening a gap and shaping the inner region of the HD~139614 disk.

We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

The mid-IR detection rate of water lines in disks around Herbig stars disks is about 5\%, while it is around 50\% for disks around TTauri stars. The reason for this is still unclear. In this study, we want to find an explanation for the different detection rates between low mass and high mass pre-main-sequence stars (PMSs) in the mid-IR regime. We run disk models with stellar parameters adjusted to spectral types B9 through M2, using the radiation thermo-chemical disk modeling code ProDiMo. We produce convolved spectra at the resolution of Spitzer IRS, JWST MIRI and VLT VISIR spectrographs. We apply random noise derived from typical Spitzer spectra for a direct comparison with observations. The strength of the mid-IR water lines correlates directly with the luminosity of the central star. We explored a small parameter space around a standard disk model, considering dust-to-gas mass ratio, disk gas mass, mixing coefficient for dust settling, flaring index, dust maximum size and size power law distribution index. The models show that it is possible to suppress the water emission, however, current observations are not sensitive enough to detect mid-IR lines in disks for most of the explored parameters. The presence of noise in the spectra, combined with the high continuum flux (noise level is proportional to the continuum flux), is the most likely explanation for the non detections towards Herbig stars. Mid-IR spectra with resolution higher than 20000 are needed to investigate water in protoplanetary disks. Intrinsic differences in disk structure, e.g. inner gaps, gas-to-dust ratio, dust size and distribution, and inner disk scale height, between Herbig and TTauri star disks are able to explain a lower water detection rate in disks around Herbig stars.

We investigate which properties of protoplanetary disks around T Tauri stars affect the physics and chemistry in the regions where mid- and far-IR water lines originate and their respective line fluxes. We search for diagnostics for future observations. With the code ProDiMo, we build a series of models exploring a large parameter space, computing rotational and rovibrational transitions of water in nonlocal thermodynamic equilibrium (non-LTE). We select a sample of transitions in the mid- IR regime and the fundamental ortho and para water transitions in the far-IR. We investigate the chemistry and the local physical conditions in the line emitting regions. We calculate Spitzer spectra for each model and compare far-IR and mid-IR lines. In addition, we use mid-IR colors to tie the water line predictions to the dust continuum. Parameters affecting the water line fluxes in disks by more than a factor of three are : the disk gas mass, the dust-to-gas mass ratio, the dust maximum grain size, ISM(InterStellarMedium) UV radiation field, the mixing parameter of Dubrulle settling, the disk flaring parameter, and the dust size distribution. The first four parameters affect the mid-IR lines much more than the far-IR lines. A key driver behind water spectroscopy is the dust opacity, which sets the location of the water line emitting region. We identify three types of parameters. Parameters, such as dust-to-gas ratio, ISM radiation field, and dust size distribution, affect the mid-IR lines more, while the far-IR transitions are more affected by the flaring index. The gas mass greatly affects lines in both regimes. Higher spectral resolution and line sensitivities, like from the James Webb Space Telescope, are needed to detect a statistically relevant sample of individual water lines to distinguish further between these types of parameters.

Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 {\mu}m (called "coreshine") seems to indicate the presence of grains with sizes above 0.5 {\mu}m. We aim to analyze 3.6 and 4.5 {\mu}m coreshine from molecular cloud cores to probe the largest grains in the size distribution. We analyzed dedicated deep Cycle 9 Spitzer IRAC observations in the 3.6 and 4.5 {\mu}m bands for a sample of 10 low-mass cores. We used a new modeling approach based on a combination of ratios of the two background- and foreground-subtracted surface brightnesses and observed limits of the optical depth. The dust grains were modeled as ice-coated silicate and carbonaceous spheres. We discuss the impact of local radiation fields with a spectral slope differing from what is seen in the DIRBE allsky maps. For the cores L260, ecc806, L1262, L1517A, L1512, and L1544, the model reproduces the data with maximum grain sizes around 0.9, 0.5, 0.65, 1.5, 0.6, and > 1.5 {\mu}m, respectively. The maximum coreshine intensities of L1506C, L1439, and L1498 in the individual bands require smaller maximum grain sizes than derived from the observed distribution of band ratios. Additional isotropic local radiation fields with a spectral shape differing from the DIRBE map shape do not remove this discrepancy. In the case of Rho Oph 9, we were unable to reliably disentangle the coreshine emission from background variations and the strong local PAH emission. Considering surface brightness ratios in the 3.6 and 4.5 {\mu}m bands across a molecular cloud core is an effective method of disentangling the complex interplay of structure and opacities when used in combination with observed limits of the optical depth.

Context : The properties of the inner disks of bright Herbig AeBe stars have been studied with near infrared (NIR) interferometry and high resolution spectroscopy. The continuum and a few molecular gas species have been studied close to the central star; however, sensitivity problems limit direct information about the inner disks of the fainter T Tauri stars. Aims : Our aim is to measure some of the properties of the inner regions of disks surrounding southern T Tauri stars. Methods : We performed a survey with the PIONIER recombiner instrument at H-band of 21 T Tauri stars. The baselines used ranged from 11 m to 129 m, corresponding to a maximum resolution of 3mas (0.45 au at 150 pc). Results : Thirteen disks are resolved well and the visibility curves are fully sampled as a function of baseline in the range 45-130 m for these 13 objects. A simple qualitative examination of visibility profiles allows us to identify a rapid drop-off in the visibilities at short baselines in 8 resolved disks. This is indicative of a significant contribution from an extended contribution of light from the disk. We demonstrate that this component is compatible with scattered light, providing strong support to a prediction made by Pinte et al. (2008). The amplitude of the drop-off and the amount of dust thermal emission changes from source to source suggesting that each disk is different. A by-product of the survey is the identification of a new milli-arcsec separation binary: WW Cha. Spectroscopic and interferometric data of AK Sco have also been fitted with a binary and disk model. Conclusions : Visibility data are reproduced well when thermal emission and scattering form dust are fully considered. The inner radii measured are consistent with the expected dust sublimation radii. Modelling of AK Sco suggests a likely coplanarity between the disk and the binary's orbital plane

Molecular gas constitutes the dominant mass component of protoplanetary discs. To date, these sources have not been studied comprehensively at the longest far-infrared and shortest submillimetre wavelengths. This paper presents Herschel SPIRE FTS spectroscopic observations toward 18 protoplanetary discs, covering the entire 450-1540 GHz (666-195 $\mu$m) range at R~400-1300. The spectra reveal clear detections of the dust continuum and, in six targets, a significant amount of spectral line emission primarily attributable to $^{12}$CO rotational lines. Other targets exhibit little to no detectable spectral lines. Low signal-to-noise detections also include signatures from $^{13}$CO, [CI] and HCN. For completeness, we present upper limits of non-detected lines in all targets, including low-energy transitions of H2O and CH$^+$ molecules. The ten $^{12}$CO lines that fall within the SPIRE FTS bands trace energy levels of ~50-500 K. Combined with lower and higher energy lines from the literature, we compare the CO rotational line energy distribution with detailed physical-chemical models, for sources where these are available and published. Our 13CO line detections in the disc around Herbig Be star HD 100546 exceed, by factors of ~10-30, the values predicted by a model that matches a wealth of other observational constraints, including the SPIRE $^{12}$CO ladder. To explain the observed $^{12}$CO/$^{13}$CO ratio, it may be necessary to consider the combined effects of optical depth and isotope selective (photo)chemical processes. Considering the full sample of 18 objects, we find that the strongest line emission is observed in discs around Herbig Ae/Be stars, although not all show line emission. In addition, two of the six T Tauri objects exhibit detectable $^{12}$CO lines in the SPIRE range.

Context: Observations of scattered light in the 1-5 $\mu$m range have revealed dust grains in molecular cores with sizes larger than commonly inferred for the diffuse interstellar medium. It is currently unclear whether these grains are grown within the molecular cores or are an ubiquitous component of the interstellar medium. Aims: We investigate whether the large grains necessary for efficient scattering at 1-5 mum are associated with the abundance of water-ice within molecular clouds and cores. Methods: We combined water-ice abundance measurements for sight lines through the Lupus IV molecular cloud complex with measurements of the scattered light at 3.6 mum for the same sight lines. Results: We find that there is a similar threshold for the cores in emission in scattered light at 3.6 mum (tau_9.7=0.15pm0.05, A_K=0.4pm0.2 as water-ice (tau_9.7=0.11pm0.01, A_K=0.19pm0.04) and that the scattering efficiency increases as the relative water-ice abundance increases. The ice layer increases the average grain size, which again strongly increases the albedo. Conclusions: The higher scattering efficiency is partly due to layering of ice on the dust grains. Although the layer can be relatively thin it can enhance the scattering substantially.

The disk around HD 142527 attracts a lot of attention, amongst others because of its resolved (sub) mm dust continuum that is concentrated into a horseshoe-shape towards the north of the star. In this manuscript we present spatially resolved ALMA detections of the HCN J=4-3 and CS J=7-6 emission lines. These lines give us a view deeper into the disk compared to the (optically thicker) CO isotopes. This is the first detection of CS J=7-6 coming from a protoplanetary disk. Both emission lines are azimuthally asymmetric and are suppressed under the horseshoe-shaped continuum emission peak. A possible mechanism to explain the decrease under the horseshoe-shaped continuum is the increased opacity coming from the higher dust concentration at the continuum peak. Lower {\gr dust and/or gas} temperatures and an optically thick radio-continuum reduce line emission by freeze-out and shielding of emission from the far side of the disk.

Investigating the evolution of protoplanetary disks is crucial for our understanding of star and planet formation. Several theoretical and observational studies have been performed in the last decades to advance this knowledge. FT Tauri is a young star in the Taurus star forming region that was included in a number of spectroscopic and photometric surveys. We investigate the properties of the star, the circumstellar disk, and the accretion and ejection processes and propose a consistent gas and dust model also as a reference for future observational studies. We performed a multi-wavelength data analysis to derive the basic stellar and disk properties, as well as mass accretion/outflow rate from TNG-Dolores, WHT-Liris, NOT-Notcam, Keck-Nirspec, and Herschel-Pacs spectra. From the literature, we compiled a complete Spectral Energy Distribution. We then performed detailed disk modeling using the MCFOST and ProDiMo codes. Multi-wavelengths spectroscopic and photometric measurements were compared with the reddened predictions of the codes in order to constrain the disk properties. This object can serve as a benchmark for primordial disks with significant mass accretion rate, high gas content and typical size.

HD 135344B is an accreting (pre-) transition disk that displays the emission of warm CO extending tens of AU inside its 30 AU dust cavity. We used the dust radiative transfer code MCFOST and the thermochemical code ProDiMo to derive the disk structure from the simultaneous modeling of the spectral energy distribution (SED), VLT/CRIRES CO P(10) 4.75 micron, Herschel/PACS [O I] 63 micron, Spitzer-IRS, and JCMT 12CO J=3-2 spectra, VLTI/PIONIER H-band visibilities, and constraints from (sub-)mm continuum interferometry and near-IR imaging. We found a disk model able to describe the current observations simultaneously. This disk has the following structure. (1) To reproduce the SED, the near-IR interferometry data, and the CO ro-vibrational emission, refractory grains (we suggest carbon) are present inside the silicate sublimation radius (0.08<R<0.2 AU). (2) The dust cavity (R<30 AU) is filled with gas, the surface density of this gas must increase with radius to fit the CO P(10) line profile, a small gap of a few AU in the gas is compatible with current data, and a large gap in the gas is not likely. (4) The gas/dust ratio inside the cavity is > 100 to account for the 870 micron continuum upper limit and the CO P(10) line flux. (5) The gas/dust ratio at 30<R<200 AU is < 10 to simultaneously describe the [O I] 63 micron line flux and the CO P(10) line profile. (6) In the outer disk, most of the mass should be located in the midplane, and a significant fraction of the dust is in large grains. Conclusions: Simultaneous modeling of the gas and dust is required to break the model degeneracies and constrain the disk structure. An increasing gas surface density with radius in the inner dust cavity echoes the effect of a migrating Jovian planet. The low gas mass (a few MJupiter) in the HD 135344B's disk suggests that it is an evolved disk that has already lost a large portion of its mass.

Transitional disks are protoplanetary disks characterized by reduced near- and mid-infrared emission with respect to full disks. This characteristic spectral energy distribution indicates the presence of an optically thin inner cavity within the dust disk believed to mark the disappearance of the primordial massive disk. We present new Herschel Space Observatory PACS spectra of [OI] 63 micron for 21 transitional disks. Our survey complements the larger Herschel GASPS program "Gas in Protoplanetary Systems" (Dent et al. 2013) by quadrupling the number of transitional disks observed with PACS at this wavelength. [OI] 63 micron traces material in the outer regions of the disk, beyond the inner cavity of most transitional disks. We find that transitional disks have [OI] 63 micron line luminosities two times fainter than their full disk counterparts. We self consistently determine various stellar properties (e.g. bolometric luminosity, FUV excess, etc.) and disk properties (e.g. disk dust mass, etc.) that could influence the [OI] 63 micron line luminosity and we find no correlations that can explain the lower [OI] 63 micron line luminosities in transitional disks. Using a grid of thermo-chemical protoplanetary disk models, we conclude that either transitional disks are less flared than full disks or they possess lower gas-to-dust ratios due to a depletion of gas mass. This result suggests that transitional disks are more evolved than their full disk counterparts, possibly even at large radii.

We discuss the chemical pre-conditions for planet formation, in terms of gas and ice abundances in a protoplanetary disk, as function of time and position, and the resulting chemical composition and cloud properties in the atmosphere when young gas giant planets form. Large deviations between the abundances of the host star and its gas giants seem likely to occur if the planet formation follows the core-accretion scenario. These deviations stem from the separate evolution of gas and dust in the disk, where the dust forms the planet cores, followed by the final run-away accretion of the left-over gas. ProDiMo protoplanetary disk models are used to predict the chemical evolution of gas and ice in the midplane. We find that cosmic rays play a crucial role in slowly un-blocking the CO, where the liberated oxygen forms water, which then freezes out quickly. Therefore, the C/O ratio in the gas phase is found to gradually increase with time, in a region bracketed by the water and CO ice-lines. In this regions, C/O is found to approach unity after about 5 Myrs, scaling with the cosmic ray ionisation rate. We then explore how the atmospheric chemistry and cloud properties in young gas giants are affected when the non-solar C/O ratios predicted by the disk models are assumed. The DRIFT cloud formation model is applied to study the formation of atmospheric clouds under the influence of varying primordial element abundances and its feedback onto the local gas. We demonstrate that element depletion by cloud formation plays a crucial role in converting an oxygen-rich atmosphere gas into carbon-rich gas when non-solar, primordial element abundances are considered as suggested by disk model.

The structure of protoplanetary disks is thought to be linked to the temperature and chemistry of their dust and gas. Whether the disk is flat or flaring depends on the amount of radiation that it absorbs at a given radius, and on the efficiency with which this is converted into thermal energy. The understanding of these heating and cooling processes is crucial to provide a reliable disk structure for the interpretation of dust continuum emission and gas line fluxes. Especially in the upper layers of the disk, where gas and dust are thermally decoupled, the infrared line emission is strictly related to the gas heating/cooling processes. We aim to study the thermal properties of the disk in the oxygen line emission region, and to investigate the relative importance of X-ray (1-120 Angstrom) and far-UV radiation (FUV, 912-2070 Angstrom) for the heating balance there. We use [OI] 63 micron line fluxes observed in a sample of protoplanetary disks of the Taurus/Auriga star forming region and compare it to the model predictions presented in our previous work. The data were obtained with the PACS instrument on board the Herschel Space Observatory as part of the Herschel Open Time Key Program GASPS (GAS in Protoplanetary diskS), published in Howard et al. (2013). Our theoretical grid of disk models can reproduce the [OI] absolute fluxes and predict a correlation between [OI] and the sum Lx+Lfuv. The data show no correlation between the [OI] line flux and the X-ray luminosity, the FUV luminosity or their sum. The data show that the FUV or X-ray radiation has no notable impact on the region where the [OI] line is formed. This is in contrast with what is predicted from our models. Possible explanations are that the disks in Taurus are less flaring than the hydrostatic models predict, and/or that other disk structure aspects that were left unchanged in our models are important. ..abridged..

Recently discovered scattered light at 3-5 $\mu$m from low-mass cores (so-called "coreshine") reveals the presence of grains around 1 $\mu$m, which is larger than the grains found in the low-density interstellar medium. But only about half of the 100+ cores investigated so far show the effect. This prompts further studies on the origin of this detection rate. From the 3D continuum radiative transfer equation, we derive the expected scattered light intensity from a core placed in an arbitrary direction seen from Earth. We use the approximation of single scattering, consider extinction up to 2nd-order Taylor approximation, and neglect spatial gradients in the dust size distribution. The impact of the directional characteristics of the scattering on the detection of scattered light from cores is calculated for a given grain size distribution, and local effects like additional radiation field components are discussed. The surface brightness profiles of a core with a 1D density profile are calculated for various Galactic locations, and the results are compared to the approximate detection limits. We find that for optically thin radiation and a constant size distribution, a simple limit for detecting scattered light from a low-mass core can be derived that holds for grains with sizes smaller than 0.5 $\mu$m. The extinction by the core prohibits detection in bright parts of the Galactic plane, especially near the Galactic center. For scattered light received from low-mass cores with grain sizes beyond 0.5 $\mu$m, the directional characteristics of the scattering favors the detection of scattered light above and below the Galactic center, and to some extent near the Galactic anti-center. We identify the local incident radiation field as the major unknown causing deviations from this simple scheme.

We present spectrally resolved Herschel/HIFI observations of the young multiple system T Tau in atomic and molecular lines. While CO, H2O, [C II], and SO lines trace the envelope and the outflowing gas up to velocities of 33 km/s with respect to systemic, the CN 5-4 hyperfine structure lines at 566.7, 566.9 GHz show a narrow double-peaked profile centered at systemic velocity, consistent with an origin in the outer region of the compact disk of T Tau N. Disk modeling of the T Tau N disk with the thermo-chemical code ProDiMo produces CN line fluxes and profiles consistent with the observed ones and constrain the size of the gaseous disk (R_out = 110 (+10, -20) AU) and its inclination (i = 25 \pm 5 degree). The model indicates that the CN lines originate in a disk upper layer at 40-110 AU from the star, which is irradiated by the stellar UV field and heated up to temperatures of 50-700 K. With respect to previously observed CN 2-1 millimeter lines, the CN 5-4 lines appear to be less affected by envelope emission, due to their larger critical density and excitation temperature. Hence, high-J CN lines are a unique confusion-free tracer of embedded disks, such as the disk of T Tau N.

We have studied the emission of CO ro-vibrational lines in the disc around the Herbig Be star HD100546 with the final goal of using these lines as a diagnostic to understand inner disc structure in the context of planet formation. High-resolution IR spectra of CO ro-vibrational emission at eight different position angles were taken with CRIRES at the VLT. From these spectra flux tables, CO ro-vibrational line profiles, and population diagrams were produced. We have investigated variations in the line profile shapes and line strengths as a function of slit position angle. We used the thermochemical disc modelling code ProDiMo based on the chemistry, radiation field, and temperature structure of a previously published model for HD100546. Comparing observations and the model, we investigated the possibility of disc asymmetries, the excitation mechanism (UV fluorescence), the geometry, and physical conditions of the inner disc. The observed CO ro-vibrational lines are largely emitted from the inner rim of the outer disc at 10-13 AU. The line shapes are similar for all v levels and line fluxes from all vibrational levels vary only within one order of magnitude. All line profile asymmetries and variations can be explained with a symmetric disc model to which a slit correction and pointing offset is applied. Because the angular size of the CO emitting region (10-13 AU) and the slit width are comparable the line profiles are very sensitive to the placing of the slit. The model reproduces the line shapes and the fluxes of the v=1-0 lines as well as the spatial extent of the CO ro-vibrational emission. It does not reproduce the observed band ratios of 0.5-0.2 with higher vibrational bands. We find that lower gas volume densities at the surface of the inner rim of the outer disc can make the fluorescence pumping more effcient and reproduce the observed band ratios.

Context. The largest grains (0.5-1 micron) in the interstellar size distribution are efficient in scattering near- and mid-infrared radiation. These wavelengths are therefore particularly well suited to probe the still uncertain high-end of the size distribution. Aims. We investigate the change in appearance of a typical low-mass molecular core from the Ks (2.2 micron) band to the Spitzer IRAC 3.6 and 8 micron bands, and compare with model calculations, which include variations of the grain size distribution. Methods. We combine Spitzer IRAC and ground-based near-infrared observations to characterize the scattered light observed at the near- and mid-infrared wavelengths from the core L260. Using a spherical symmetric model core, we perform radiative transfer calculations to study the impact of various dust size distributions on the intensity profiles across the core. Results. The observed scattered light patterns in the Ks and 3.6 micron bands are found to be similar. By comparison with radiative transfer models the two profiles place constraints on the relative abundance of small and large (more than 0.25 micron) dust grains. The scattered light profiles are found to be inconsistent with an interstellar silicate grain distribution extending only to 0.25 micron and large grains are needed to reach the observed fluxes and the flux ratios. The shape of the Ks band surface brightness profile limits the largest grains to 1-1.5 micron. Conclusions. In addition to observing coreshine in the Spitzer IRAC channels, the combination with ground-based near-infrared observations are suited to constrain the properties of large grains in cores.

At the distance of 99-116 pc, HD141569A is one of the nearest HerbigAe stars that is surrounded by a tenuous disk, probably in transition between a massive primordial disk and a debris disk. We observed the fine-structure lines of OI at 63 and 145 micron and the CII line at 157 micron with the PACS instrument onboard the Herschel Space Telescope as part of the open-time large programme GASPS. We complemented the atomic line observations with archival Spitzer spectroscopic and photometric continuum data, a ground-based VLT-VISIR image at 8.6 micron, and 12CO fundamental ro-vibrational and pure rotational J=3-2 observations. We simultaneously modeled the continuum emission and the line fluxes with the Monte Carlo radiative transfer code MCFOST and the thermo-chemical code ProDiMo to derive the disk gas- and dust properties assuming no dust settling. The models suggest that the oxygen lines are emitted from the inner disk around HD141569A, whereas the [CII] line emission is more extended. The CO submillimeter flux is emitted mostly by the outer disk. Simultaneous modeling of the photometric and line data using a realistic disk structure suggests a dust mass derived from grains with a radius smaller than 1 mm of 2.1E-7 MSun and from grains with a radius of up to 1 cm of 4.9E-6 MSun. We constrained the polycyclic aromatic hydrocarbons (PAH) mass to be between 2E-11 and 1..4E-10 MSun assuming circumcircumcoronene (C150H30) as the representative PAH. The associated PAH abundance relative to hydrogen is lower than those found in the interstellar medium (3E-7) by two to three orders of magnitude. The disk around HD141569A is less massive in gas (2.5 to 4.9E-4 MSun or 67 to 164 MEarth) and has a flat opening angle (<10%). [abridged]

The Herschel Space Observatory was used to observe ~ 120 pre-main-sequence stars in Taurus as part of the GASPS Open Time Key project. PACS was used to measure the continuum as well as several gas tracers such as [OI] 63 \mu m, [OI] 145 \mu m, [CII] 158 \mu m, OH, H2O and CO. The strongest line seen is [OI] at 63 \mu m. We find a clear correlation between the strength of the [OI] 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 ($<$ 50 AU) and lower surface layers of the disk where the gas and dust are coupled. The [OI] 63 \mu m is fainter in transitional stars than in normal Class II disks. Simple SED models indicate that the dust responsible for the continuum emission is colder in these disks, leading to weaker line emission. [CII] 158 \mu m emission is only detected in strong outflow sources. The observed line ratios of [OI] 63 \mu m to [OI] 145 \mu m are in the regime where we are insensitive to the gas-to-dust ratio, neither can we discriminate between shock or PDR emission. We detect no Class III object in [OI] 63 \mu m and only three in continuum, at least one of which is a candidate debris disk.

Context. This paper discusses the sensitivity of water lines to chemical processes and radiative transfer for the protoplanetary disk around TW Hya. The study focuses on the Herschel spectral range in the context of new line detections with the PACS instrument from the Gas in Protoplanetary Systems project (GASPS). Aims. The paper presents an overview of the chemistry in the main water reservoirs in the disk around TW Hya. It discusses the limitations in the interpretation of observed water line fluxes. Methods. ... (abbreviated) Results. We report new line detections of p-H2O (3_22-2_11) at 89.99 micron and CO J=18-17 at 144.78 micron for the disk around TW Hya. Disk modeling shows that the far-IR fine structure lines ([OI], [CII]) and molecular submm lines are very robust to uncertainties in the chemistry, while the water line fluxes can change by factors of a few. The water lines are optically thick, sub-thermally excited and can couple to the background continuum radiation field. The low-excitation water lines are also sensitive to uncertainties in the collision rates, e.g. with neutral hydrogen. The gas temperature plays an important role for the [OI] fine structure line fluxes, the water line fluxes originating from the inner disk as well as the high excitation CO, CH+ and OH lines. Conclusions. Due to their sensitivity on chemical input data and radiative transfer, water lines have to be used cautiously for understanding details of the disk structure. Water lines covering a wide range of excitation energies provide access to the various gas phase water reservoirs (inside and outside the snow line) in protoplanetary disks and thus provide important information on where gas-phase water is potentially located. Experimental and/or theoretical collision rates for H2O with atomic hydrogen are needed to diminish uncertainties from water line radiative transfer.

We present an extension of the code ProDiMo that allows for a modeling of processes pertinent to active galactic nuclei and to an ambient chemistry that is time dependent. We present a proof-of-concept and focus on a few astrophysically relevant species, e.g., H+, H2+ and H3+; C+ and N+; C and O; CO and H2O; OH+, H2O+ and H3O+; HCN and HCO+. We find that the freeze-out of water is strongly suppressed and that this affects the bulk of the oxygen and carbon chemistry occurring in AGN. The commonly used AGN tracer HCN/HCO+ is strongly time-dependent, with ratios that vary over orders of magnitude for times longer than 10^4 years. Through ALMA observations this ratio can be used to probe how the narrow-line region evolves under large fluctuations in the SMBH accretion rate. Strong evolutionary trends, on time scales of 10^4-10^8 years, are also found in species such as H3O+, CO, and H2O. These reflect, respectively, time dependent effects in the ionization balance, the transient nature of the production of molecular gas, and the freeze-out/sublimation of water.

Context. Circumstellar discs are the places where planets form, therefore knowledge of their evolution is crucial for our understanding of planet formation. The Herschel Space Observatory is providing valuable data for studying disc systems, thanks to its sensitivity and wavelength coverage. This paper is one of several devoted to analysing and modelling Herschel-PACS observations of various young stellar associations from the GASPS Open Time Key Programme. Aims. The aim of this paper is to elucidate the gas and dust properties of circumstellar discs in the 10 Myr TW Hya Association (TWA) using new far-infrared (IR) imaging and spectroscopy from Herschel-PACS. Methods. We obtained far-IR photometric data at 70, 100, and 160 microns of 14 TWA members; spectroscopic observations centred on the [OI] line at 63.18 microns were also obtained for 9 of the 14. The new photometry for each star was incorporated into its full spectral energy distribution (SED). Results. We detected excess IR emission that is characteristic of circumstellar discs from five TWA members, and computed upper limits for another nine. Two TWA members (TWA 01 and TWA 04B) also show [OI] emission at 63.18 microns. Discs in the TWA association display a variety of properties, with a wide range of dust masses and inner radii, based on modified blackbody modelling. Both transitional and debris discs are found in the sample. Models for sources with a detected IR excess give dust masses in the range from 0.15 Msun to 63 Msun.

GASPS is a far-infrared line and continuum survey of protoplanetary and young debris disks using PACS on the Herschel Space Observatory. The survey includes [OI] at 63 microns, as well as 70, 100 and 160um continuum, with the brightest objects also studied in [OI]145um, [CII]157um, H2O and CO. Targets included T Tauri stars and debris disks in 7 nearby young associations, and a sample of isolated Herbig AeBe stars. The aim was to study the global gas and dust content in a wide disk sample, systemically comparing the results with models. In this paper we review the main aims, target selection and observing strategy. We show initial results, including line identifications, sources detected, and a first statistical study. [OI]63um was the brightest line in most objects, by a factor of ~10. Detection rates were 49%, including 100% of HAeBe stars and 43% of T Tauri stars. Comparison with published dust masses show a dust threshold for [OI]63um detection of ~1e-5 M_solar. Normalising to 140pc distance, 32% with mass 1e-6 - 1e-5 M_solar, and a small number with lower mass were also detected. This is consistent with moderate UV excess and disk flaring. In most cases, continuum and line emission is spatially and spectrally unresolved, suggesting disk emission. ~10 objects were resolved, likely from outflows. Detection rates in [OI]145um, [CII]157um and CO J=18-17 were 20-40%, but [CII] was not correlated with disk mass, suggesting it arises instead from a compact envelope. [OI] detection rates in T Tauri associations of ages 0.3-4Myr were ~50%. ~2 stars were detectable in associations of 5-20Myr, with no detections in associations of age >20Myr. Comparing with the total number of young stars, and assuming a ISM-like gas/dust ratio, this indicates that ~18% of stars retain a gas-rich disk of total mass >1M_Jupiter for 1-4Myr, 1-7% keep such disks for 5-10Myr, and none remain beyond 10-20Myr.

We present far-IR/sub-mm imaging and spectroscopy of 49 Ceti, an unusual circumstellar disk around a nearby young A1V star. The system is famous for showing the dust properties of a debris disk, but the gas properties of a low-mass protoplanetary disk. The data were acquired with the Herschel Space Observatory PACS and SPIRE instruments, largely as part of the "Gas in Protoplanetary Systems" (GASPS) Open Time Key Programme. Disk dust emission is detected in images at 70, 160, 250, 350, and 500 \mu m; 49 Cet is significantly extended in the 70 \mu m image, spatially resolving the outer dust disk for the first time. Spectra covering small wavelength ranges centered on eight atomic and molecular emission lines were obtained, including [OI] 63 \mu m and [CII] 158 \mu m. The CII line was detected at the 5\sigma\ level - the first detection of atomic emission from the disk. No other emission lines were seen, despite the fact that the OI line is the brightest one observed in Herschel protoplanetary disk spectra (Meeus et al. 2012; Dent et al. 2013). We present an estimate of the amount of circumstellar atomic gas implied by the CII emission. The new far-IR/sub-mm data fills in a large gap in the previous spectral energy distribution (SED) of 49 Cet. A simple model of the new SED confirms the two-component structure of the disk: warm inner dust and cold outer dust that produces most of the observed excess. Finally, we discuss preliminary thermochemical modeling of the 49 Cet gas/dust disk and our attempts to match several observational results simultaneously. Although we are not yet successful in doing so, our investigations shed light on the evolutionary status of the 49 Cet gas, which might not be primordial gas but rather secondary gas coming from comets.

Circumstellar disc evolution is paramount for the understanding of planet formation. The GASPS program aims at determining the circumstellar gas and solid mass around ~250 pre-main-sequence Herbig Ae and TTauri stars. We aim to understand the origin and nature of the circumstellar matter orbiting 51 Oph, a young (<1 Myr) luminous B9.5 star. We obtained continuum and line observations with the PACS instrument on board the Herschel Space Observatory and continuum data at 1.2 mm with the IRAM 30m telescope. The SED and line fluxes were modelled using the physico-chemo radiative transfer code ProDiMo. We detected a strong emission by OI at 63 microns using the Herschel Space Observatory. The [OI] emission at 145 microns, the [CII] emission at 158 microns, the high-J CO emissions, and the warm water emissions were not detected. Continuum emission was detected at 1.2 mm. The continuum from the near- to the far-infrared and the [OI] emission are well explained by the emission from a compact hydrostatic disc model with a gas mass of 5E-6 MSun, 100 times that of the solid mass. However, this model fails to match the continuum millimeter flux, which hints at a cold outer disc with a mass in solids of 1E-6 MSun or free-free emission from a photoevaporative disc wind. This outer disc can either be devoid of gas and/or is to cold to emit in the [OI] line. A very flat extended disc model (Rout=400 AU) with a fixed vertical structure and dust settling matches all photometric points and most of the [O I] flux. The observations can be explained by an extended flat disc where dust grains have settled. However, a flat gas disc cannot be reproduced by hydrostatic disc models. The low mass of the 51 Oph inner disc in gas and dust may be explained either by the fast dissipation of an initial massive disc or by a very small initial disc mass.

Water is key in the evolution of protoplanetary disks and the formation of comets and icy/water planets. While high excitation water lines originating in the hot inner disk have been detected in several T Tauri stars (TTSs), water vapor from the outer disk, where most of water ice reservoir is stored, was only reported in the closeby TTS TW Hya. We present spectrally resolved Herschel/HIFI observations of the young TTS DG Tau in the ortho- and para- water ground-state transitions at 557, 1113 GHz. The lines show a narrow double-peaked profile, consistent with an origin in the outer disk, and are ~19-26 times brighter than in TW Hya. In contrast, CO and [C II] lines are dominated by emission from the envelope/outflow, which makes H2O lines a unique tracer of the disk of DG Tau. Disk modeling with the thermo-chemical code ProDiMo indicates that the strong UV field, due to the young age and strong accretion of DG Tau, irradiates a disk upper layer at 10-90 AU from the star, heating it up to temperatures of 600 K and producing the observed bright water lines. The models suggest a disk mass of 0.015-0.1 Msun, consistent with the estimated minimum mass of the solar nebula before planet formation, and a water reservoir of ~1e2-1e3 Earth oceans in vapour, and ~100 times larger in the form of ice. Hence, this detection supports the scenario of ocean delivery on terrestrial planets by impact of icy bodies forming in the outer disk.

The carbon monoxide rovibrational emission from discs around Herbig Ae stars and T Tauri stars with strong ultraviolet emissions suggests that fluorescence pumping from the ground X1 Sigma+ to the electronic A1 Pi state of CO should be taken into account in disc models. We implemented a CO model molecule that includes up to 50 rotational levels within nine vibrational levels for the ground and A excited states in the radiative photochemical code ProDiMo. We took CO collisions with hydrogen molecules, hydrogen atoms, helium, and electrons into account. We estimated the missing collision rates using standard scaling laws and discussed their limitations. UV fluorescence and IR pumping impact on the population of ro-vibrational v > 1 levels. The v = 1 rotational levels are populated at rotational temperatures between the radiation temperature around 4.6 micron and the gas kinetic temperature. The UV pumping efficiency increases with decreasing disc mass. The consequence is that the vibrational temperatures, which measure the relative populations between the vibrational levels, are higher than the disc gas kinetic temperatures (suprathermal population). Rotational temperatures from fundamental transitions derived using optically thick 12CO lines do not reflect the gas kinetic temperature. CO pure rotational levels with energies lower than 1000 K are populated in LTE but are sensitive to a number of vibrational levels included in the model. The 12CO pure rotational lines are highly optically thick for transition from levels up to Eupper=2000 K. (abridged)

Context. HD 172555 is a young A7 star belonging to the Beta Pictoris Moving Group that harbours a debris disc. The Spitzer IRS spectrum of the source showed mid-IR features such as silicates and glassy silica species, indicating the presence of a warm dust component with small grains, which places HD 172555 among the small group of debris discs with such properties. The IRS spectrum also shows a possible emission of SiO gas. Aims. We aim to study the dust distribution in the circumstellar disc of HD 172555 and to asses the presence of gas in the debris disc. Methods. As part of the GASPS Open Time Key Programme, we obtained Herschel-PACS photometric and spectroscopic observations of the source. We analysed PACS observations of HD 172555 and modelled the Spectral Energy Distribution (SED) with a modified blackbody and the gas emission with a two-level population model with no collisional de-excitation. Results. We report for the first time the detection of OI atomic gas emission at 63.18 microns in the HD 172555 circumstellar disc. We detect excesses due to circumstellar dust toward HD 172555 in the three photometric bands of PACS (70, 100, and 160 microns). We derive a large dust particle mass of 4.8e-4 Earth masses and an atomic oxygen mass of 2.5e-2*R^2 Earth masses, where R in AU is the separation between the star and the inner disc. Thus, most of the detected mass of the disc is in the gaseous phase.

Most of the mass in protoplanetary disks is in the form of gas. The study of the gas and its diagnostics is of fundamental importance in order to achieve a detailed description of the thermal and chemical structure of the disk. The radiation from the central star (from optical to X-ray wavelengths) and viscous accretion are the main source of energy and dominates the disk physics and chemistry in its early stages. This is the environment in which the first phases of planet formation will proceed. We investigate how stellar and disk parameters impact the fine-structure cooling lines [NeII], [ArII], [OI], [CII] and H2O rotational lines in the disk. These lines are potentially powerful diagnostics of the disk structure and their modelling permits a thorough interpretation of the observations carried out with instrumental facilities such as Spitzer and Herschel. Following Aresu et al. (2011), we computed a grid of 240 disk models, in which the X-ray luminosity, UV-excess luminosity, minimum dust grain size, dust size distribution power law and surface density distribution power law, are systematically varied. We solve self-consistently for the disk vertical hydrostatic structure in every model and apply detailed line radiative transfer to calculate line fluxes and profiles for a series of well known mid- and far-infrared cooling lines. The [OI] 63 micron line flux increases with increasing FUV luminosity when Lx < 1e30 erg/s, and with increasing X-ray luminosity when LX > 1e30 erg/s. [CII] 157 micron is mainly driven by FUV luminosity via C+ production, X-rays affect the line flux to a lesser extent. [NeII] 12.8 micron correlates with X-rays; the line profile emitted from the disk atmosphere shows a double-peaked component, caused by emission in the static disk atmosphere, next to a high velocity double-peaked component, caused by emission in the very inner rim. (abridged)

Context. Planets are thought to eventually form from the mostly gaseous (~99% of the mass) disks around young stars. The density structure and chemical composition of protoplanetary disks are affected by the incident radiation field at optical, FUV, and X-ray wavelengths, as well as by the dust properties. Aims. The effect of FUV and X-rays on the disk structure and the gas chemical composition are investigated. This work forms the basis of a second paper, which discusses the impact on diagnostic lines of, e.g., C+, O, H2O, and Ne+ observed with facilities such as Spitzer and Herschel. Methods. A grid of 240 models is computed in which the X-ray and FUV luminosity, minimum grain size, dust size distribution, and surface density distribution are varied in a systematic way. The hydrostatic structure and the thermo-chemical structure are calculated using ProDiMo. Results. The abundance structure of neutral oxygen is stable to changes in the X-ray and FUV luminosity, and the emission lines will thus be useful tracers of the disk mass and temperature. The C+ abundance distribution is sensitive to both X-rays and FUV. The radial column density profile shows two peaks, one at the inner rim and a second one at a radius r=5-10 AU. Ne+ and other heavy elements have a very strong response to X-rays, and the column density in the inner disk increases by two orders of magnitude from the lowest (LX = 1e29 erg/s) to the highest considered X-ray flux (LX = 1e32 erg/s). FUV confines the Ne+ ionized region to areas closer to the star at low X-ray luminosities (LX = 1e29 erg/s). H2O abundances are enhanced by X-rays due to higher temperatures in the inner disk and higher ionization fractions in the outer disk. The line fluxes and profiles are affected by the effects on these species, thus providing diagnostic value in the study of FUV and X-ray irradiated disks around T Tauri stars. (abridged)

Cosmic rays (CRs) control the thermal, ionization and chemical state of the dense H_2 gas regions that otherwise remain shielded from far-UV and optical stellar radiation propagating through the dusty ISM of galaxies. It is in such CR-dominated regions (CRDRs) rather than Photon-dominated regions (PDRs) of H_2 clouds where the star formation initial conditions are set, making CRs the ultimate star-formation feedback factor in galaxies, able to operate even in their most deeply dust-enshrouded environments. CR-controlled star formation initial conditions naturally set the stage for a near-invariant stellar Initial Mass Function (IMF) in galaxies as long as their average CR energy density U_{CR} permeating their molecular ISM remains within a factor of ~10 of its Galactic value. Nevertheless, in the extreme environments of the compact starbursts found in merging galaxies, where U_{CR}\sim(few)x10^{3}U_{CR,Gal}, CRs dramatically alter the initial conditions of star formation. In the resulting extreme CRDRs H_2 cloud fragmentation will produce far fewer low mass (<8 M_{sol}) stars, yielding a top-heavy stellar IMF. This will be a generic feature of CR-controlled star-formation initial conditions, lending a physical base for a bimodal IMF during galaxy formation, with a top-heavy one for compact merger-induced starbursts, and an ordinary IMF preserved for star formation in isolated gas-rich disks. In this scheme the integrated galactic IMFs (IGIMF) are expected to be strong functions of the star formation history of galaxies.

We observed a sample of 20 representative Herbig Ae/Be stars and five A-type debris discs with PACS onboard of Herschel. The observations were done in spectroscopic mode, and cover far-IR lines of [OI], [CII], CO, CH+, H2O and OH. We have a [OI]63 micron detection rate of 100% for the Herbig Ae/Be and 0% for the debris discs. [OI]145 micron is only detected in 25%, CO J=18-17 in 45% (and less for higher J transitions) of the Herbig Ae/Be stars and for [CII] 157 micron, we often found spatially variable background contamination. We show the first detection of water in a Herbig Ae disc, HD 163296, which has a settled disc. Hydroxyl is detected as well in this disc. CH+, first seen in HD 100546, is now detected for the second time in a Herbig Ae star, HD 97048. We report fluxes for each line and use the observations as line diagnostics of the gas properties. Furthermore, we look for correlations between the strength of the emission lines and stellar or disc parameters, such as stellar luminosity, UV and X-ray flux, accretion rate, PAH band strength, and flaring. We find that the stellar UV flux is the dominant excitation mechanism of [OI]63 micron, with the highest line fluxes found in those objects with a large amount of flaring and greatest PAH strength. Neither the amount of accretion nor the X-ray luminosity has an influence on the line strength. We find correlations between the line flux of [OI]63 micron and [OI]145 micron, CO J = 18-17 and [OI]6300 \AA, and between the continuum flux at 63 micron and at 1.3 mm, while we find weak correlations between the line flux of [OI]63 micron and the PAH luminosity, the line flux of CO J = 3-2, the continuum flux at 63 micron, the stellar effective temperature and the Brgamma luminosity. (Abbreviated version)

Line spectra of 68 Taurus T Tauri stars were obtained with the Herschel-PACS (Photodetector Array Camera & Spectrometer) instrument as part of the GASPS (Gas Evolution in Protoplanetary Systems) survey of protoplanetary discs. A careful examination of the line scans centred on the [OI] 63.18 microns fine-structure line unveiled a line at 63.32 micron in some of these spectra. We identify this line with a transition of ortho-water. It is detected confidently (i.e., >3 sigma) in eight sources, i.e., 24% of the sub-sample with gas-rich discs. Several statistical tests were used to search for correlations with other disc and stellar parameters such as line fluxes of [OI] 6300 Armstrong and 63.18 microns; X-ray luminosity and continuum levels at 63 microns and 850 microns. Correlations are found between the water line fluxes and the [OI] 63.18 microns line luminosity, the dust continuum, and possibly with the stellar X-ray luminosity. This is the first time that this line of warm water vapour has been detected in protoplanetary discs. We discuss its origins, in particular whether it comes from the inner disc and/or disc surface or from shocks in outflows and jets. Our analysis favours a disc origin, with the observed water vapour line produced within 2-3AU from the central stars, where the gas temperature is of the order of 500-600 K.

HD 181327 is a young F5/F6V star belonging to the Beta Pictoris moving group (12 Myr). It harbors an optically thin belt of circumstellar material at 90 AU. We aim to study the dust properties in the belt in details, and to constrain the gas-to-dust ratio. We obtained far-IR observations with the Herschel/PACS instrument, and 3.2 mm observations with the ATCA array. The geometry of the belt is constrained with newly reduced HST/NICMOS images that break the degeneracy between the disk geometry and the dust properties. We use the radiative transfer code GRaTer to compute a large grid of models, and we identify the grain models that best reproduce the Spectral Energy Distribution through a Bayesian analysis. We attempt to detect the [OI] and [CII] lines with PACS spectroscopy, providing observables to our photochemical code ProDiMo. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected in the far-IR with PACS and the disk is resolved with both PACS and ATCA. A medium integration of the gas spectral lines only provides upper limits on the line fluxes. We show that the HD 181327 dust disk consists of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an important layer of ice, for a total dust mass of 0.05 M\oplus (up to 1 mm). We discuss evidences that the grains are fluffy aggregates. The upper limits on the gas lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than 17 M\oplus. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive Kuiper Belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets.

We present detailed model fits to observations of the disc around the Herbig Ae star HD 163296. This well-studied object has an age of ~ 4 Myr, with evidence of a circumstellar disc extending out to ~ 540AU. We use the radiation thermo-chemical disc code ProDiMo to model the gas and dust in the circumstellar disc of HD 163296, and attempt to determine the disc properties by fitting to observational line and continuum data. These include new Herschel/PACS observations obtained as part of the open-time key program GASPS (Gas in Protoplanetary Systems), consisting of a detection of the [OI]63mic line and upper limits for several other far infrared lines. We complement this with continuum data and ground-based observations of the 12CO 3-2, 2-1 and 13CO J=1-0 line transitions, as well as the H2 S(1) transition. We explore the effects of stellar ultraviolet variability and dust settling on the line emission, and on the derived disc properties. Our fitting efforts lead to derived gas/dust ratios in the range 9-100, depending on the assumptions made. We note that the line fluxes are sensitive in general to the degree of dust settling in the disc, with an increase in line flux for settled models. This is most pronounced in lines which are formed in the warm gas in the inner disc, but the low excitation molecular lines are also affected. This has serious implications for attempts to derive the disc gas mass from line observations. We derive fractional PAH abundances between 0.007 and 0.04 relative to ISM levels. Using a stellar and UV excess input spectrum based on a detailed analysis of observations, we find that the all observations are consistent with the previously assumed disc geometry.

Despite its importance in the thermal-balance of the gas and in the determination of primeval planetary atmospheres, the chemistry in protoplanetary discs remains poorly constrained with only a handful of detected species. We observed the emission from disc around the Herbig Be star HD 100546 with the PACS instrument in the spectroscopic mode on board the Herschel Space Telescope as part of the Gas in Protoplanetary Systems (GASPS) programme and used archival data from the DIGIT programme to search for the rotational emission of CH+. We detected in both datasets an emission line centred at 72.16 micron that most likely corresponds to the J=5-4 rotational emission of CH+. The J=3-2 and 6-5 transitions are also detected albeit with lower confidence. Other CH+ rotational lines in the PACS observations are blended with water lines. A rotational diagram analysis shows that the CH+ gas is warm at 323 (+2320/-151) K with a mass of 3e-14-5e-12 M_Sun. We modelled the CH+ chemistry with the chemo-physical code ProDiMo using a disc density structure and grain parameters that match continuum observations and near- and mid-infrared interferometric data. The model suggests that CH+ is most abundant at the location of the disc rim at 10-13 AU from the star where the gas is warm, consistent with previous observations of hot CO gas emission.

Studying the physical conditions in circumstellar disks is a crucial step toward understanding planet formation. Of particular interest is the case of HD 100546, a Herbig Be star that presents a gap within the first 13 AU of its protoplanetary disk, that may originate in the dynamical interactions of a forming planet. We gathered a large amount of new interferometric data using the AMBER/VLTI instrument in the H- and K-bands to spatially resolve the warm inner disk and constrain its structure. Then, combining these measurements with photometric observations, we analyze the circumstellar environment of HD 100546 in the light of a passive disk model based on 3D Monte-Carlo radiative transfer. Finally, we use hydrodynamical simulations of gap formation by planets to predict the radial surface density profile of the disk and test the hypothesis of ongoing planet formation. The SED and the NIR interferometric data are adequately reproduced by our model. We show that the H- and K-band emissions are coming mostly from the inner edge of the internal dust disk, located near 0.24 AU from the star, i.e., at the dust sublimation radius in our model. We directly measure an inclination of $33^{\circ} \pm 11^{\circ}$ and a position angle of $140^{\circ} \pm 16^{\circ}$ for the inner disk. This is similar to the values found for the outer disk ($i \simeq 42^{\circ}$, $PA \simeq 145^{\circ}$), suggesting that both disks may be coplanar. We finally show that 1 to 8 Jupiter mass planets located at $\sim 8$ AU from the star would have enough time to create the gap and the required surface density jump of three orders of magnitude between the inner and outer disk. However, no information on the amount of matter left in the gap is available, which precludes us from setting precise limits on the planet mass, for now.

We present high resolution (R = 100,000) L-band spectroscopy of 11 Herbig AeBe stars with circumstellar disks. The observations were obtained with the VLT/CRIRES to detect hot water and hydroxyl radical emission lines previously detected in disks around T Tauri stars. OH emission lines are detected towards 4 disks. The OH P4.5 (1+,1-) doublet is spectrally resolved as well as the velocity profile of each component of the doublet. Its characteristic double-peak profile demonstrates that the gas is in Keplerian rotation and points to an emitting region extending out to ~ 15-30 AU. The OH, emission correlates with disk geometry as it is mostly detected towards flaring disks. None of the Herbig stars analyzed here show evidence of hot water vapor at a sensitivity similar to that of the OH lines. The non-detection of hot water vapor emission indicates that the atmosphere of disks around Herbig AeBe stars are depleted of water molecules. Assuming LTE and optically thin emission we derive a lower limit to the OH/H2O column density ratio > 1 - 25 in contrast to T Tauri disks for which the column density ratio is 0.3 -- 0.4.

Aims. We want to understand the chemistry and physics of disks on the basis of a large unbiased and statistically relevant grid of disk models. One of the main goals is to explore the diagnostic power of various gas emission lines and line ratios for deriving main disk parameters such as the gas mass. Methods. We explore the results of the DENT grid (Disk Evolution with Neat Theory) that consists of 300 000 disk models with 11 free parameters. Through a statistical analysis, we search for correlations and trends in an effort to find tools for disk diagnostic. Results. All calculated quantities like species masses, temperatures, continuum and line fluxes differ by several orders of magnitude across the entire parameter space. The broad distribution of these quantities as a function of input parameters shows the limitation of using a prototype T Tauri or Herbig Ae/Be disk model. The statistical analysis of the DENT grid shows that CO gas is rarely the dominant carbon reservoir in disks. Models with large inner radii (10 times the dust condensation radius) and/or shallow surface density gradients lack massive gas phase water reservoirs. Also, 60% of the disks have gas temperatures averaged over the oxygen mass in the range between 15 and 70 K; the average gas temperatures for CO and O differ by less than a factor two. Studying the observational diagnostics, the [CII] 158 \mum fine structure line flux is very sensitive to the stellar UV flux and presence of a UV excess and it traces the outer disk radius (Rout). In the submm, the CO low J rotational lines also trace Rout. Low [OI] 63/145 line ratios (< a few) can be explained with cool atomic O gas in the uppermost surface layers leading to self-absorption in the 63 \mum line; this occurs mostly for massive non-flaring, settled disk models without UV excess. ... abbreviated

We present new continuum and line observations, along with modelling, of the faint (6-8) Myr old T Tauri star ET Cha belonging to the eta Chamaeleontis cluster. We have acquired HERSCHEL/PACS photometric fluxes at 70 mic and 160 mic, as well as a detection of the [OI] 63 mic fine-structure line in emission, and derived upper limits for some other far-IR OI, CII, CO and o-H2O lines. The HERSCHEL data is complemented by new ANDICAM B-K photometry, new HST/COS and HST/STIS UV-observations, a non-detection of CO J=3-2 with APEX, re-analysis of a UCLES high-resolution optical spectrum showing forbidden emission lines like [OI] 6300A, [SII] 6731A and 6716A, and [NII] 6583A, and a compilation of existing broad-band photometric data. We used the thermo-chemical disk code ProDiMo and the Monte-Carlo radiative transfer code MCFOST to model the protoplanetary disk around ET Cha. Based on these models we can determine the disk dust mass Mdust = (2.E-8 - 5.E-8) Msun, whereas the total disk gas mass is found to be only little constrained, Mgas = (5.E-5 - 3.E-3) Msun. In the models, the disk extends from 0.022 AU (just outside of the co-rotation radius) to only about 10 AU. Larger disks are found to be inconsistent with the CO J=3-2 non-detection. The low velocity component of the [OI] 6300A emission line is consistent with being emitted from the inner disk. The model can also reproduce the line flux of H2 v=1-0 S(1) at 2.122 mic. An additional high-velocity component of the [OI] 6300A emission line, however, points to the existence of an additional jet/outflow of low velocity (40 - 65) km/s with mass loss rate ~1.E-9 Msun/yr. In relation to our low estimations of the disk mass, such a mass loss rate suggests a disk lifetime of only ~(0.05 - 3) Myr, substantially shorter than the cluster age. The evolutionary state of this unusual protoplanetary disk is discussed.

We examine in detail the recent proposal that extreme Cosmic-Ray-Dominated-Regions (CRDRs) characterize the ISM of galaxies during events of high-density star formation, fundamentally altering its initial conditions (Papadopoulos 2010). Solving the coupled chemical and thermal state equations for dense UV-shielded gas reveals that the large cosmic ray energy densities in such systems (U_{CR} (few)x(10^3-10^4) U_{CR,Gal}) will indeed raise the minimum temperature of this phase (where the initial conditions of star formation are set) from ~10K (as in the Milky Way) to (50-100)K. Moreover in such extreme CRDRs the gas temperature remains fully decoupled from that of the dust, with T_{kin} >> T_{dust}, even at high densities (n(H_2)~10^5--10^6 cm^{-3}), quite unlike CRDRs in the Milky Way where T_k T_{dust} when n(H_2) >= 10^5 cm^{-3}. These dramatically different star formation initial conditions will: a) boost the Jeans mass of UV-shielded gas regions by factors of ~10--100 with respect to those in quiescent or less extreme star forming systems, and b) "erase" the so-called inflection point of the effective equation of state (EOS) of molecular gas. Both these effects occur across the entire density range of typical molecular clouds, and may represent {\it a new paradigm for all high-density star formation in the Universe}, with cosmic rays as the key driving mechanism, operating efficiently even in the high dust extinction environments of extreme starbursts...

Context: T Tauri stars have X-ray luminosities ranging from L_X = 10^28-10^32 erg/s. These luminosities are similar to UV luminosities (L_UV 10^30-10^31 erg/s) and therefore X-rays are expected to affect the physics and chemistry of the upper layers of their surrounding protoplanetary disks. Aim: The effects and importance of X-rays on the chemical and hydrostatic structure of protoplanetary disks are investigated, species tracing X-ray irradiation (for L_X >= 10^29 erg/s) are identified and predictions for [OI], [CII] and [NII] fine structure line fluxes are provided. Methods: We have implemented X-ray physics and chemistry into the chemo-physical disk code ProDiMo. We include Coulomb heating and H2 ionization as heating processes and primary and secondary ionization due to X-rays in the chemistry. Results: X-rays heat up the gas causing it to expand in the optically thin surface layers. Neutral molecular species are not much affected in their abundance and spatial distribution, but charged species such as N+, OH+, H2O+ and H3O+ show enhanced abundances in the disk surface. Conclusions: Coulomb heating by X-rays changes the vertical structure of the disk, yielding temperatures of ~ 8000 K out to distances of 50 AU. The chemical structure is altered by the high electron abundance in the gas in the disk surface, causing an efficient ion-molecule chemistry. The products of this, OH+, H2O+ and H3O+, are of great interest for observations of low-mass young stellar objects with the Herschel Space Observatory. [OI] (at 63 and 145 mic) and [CII] (at 158 mic) fine structure emission are only affected for L_X > 10^30 erg/s.

We study the hydrostatic density structure of the inner disc rim around HerbigAe stars using the thermo-chemical hydrostatic code ProDiMo. We compare the Spectral Energy Distributions (SEDs) and images from our hydrostatic disc models to that from prescribed density structure discs. The 2D continuum radiative transfer in ProDiMo includes isotropic scattering. The dust temperature is set by the condition of radiative equilibrium. In the thermal-decoupled case the gas temperature is governed by the balance between various heating and cooling processes. The gas and dust interact thermally via photoelectrons, radiatively, and via gas accommodation on grain surfaces. As a result, the gas is much hotter than in the thermo-coupled case, where the gas and dust temperatures are equal, reaching a few thousands K in the upper disc layers and making the inner rim higher. A physically motivated density drop at the inner radius ("soft-edge") results in rounded inner rims, which appear ring-like in near-infrared images. The combination of lower gravity pull and hot gas beyond ~1 AU results in a disc atmosphere that reaches a height over radius ratio z/r of 0.1 while this ratio is 0.2 only in the thermo-coupled case. This puffed-up disc atmosphere intercepts larger amount of stellar radiation, which translates into enhanced continuum emission in the 3- 30 micron wavelength region from hotter grains at ~500 K. We also consider the effect of disc mass and grain size distribution on the SEDs self-consistently feeding those quantities back into the gas temperature, chemistry, and hydrostatic equilibrium computation.

Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.

Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.

In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.

The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.

We have combined the thermo-chemical disc code ProDiMo with the Monte Carlo radiative transfer code MCFOST to calculate a grid of ~300000 circumstellar disc models, systematically varying 11 stellar, disc and dust parameters including the total disc mass, several disc shape parameters and the dust-to-gas ratio. For each model, dust continuum and line radiative transfer calculations are carried out for 29 far IR, sub-mm and mm lines of [OI], [CII], 12CO and o/p-H2O under 5 inclinations. The grid allows to study the influence of the input parameters on the observables, to make statistical predictions for different types of circumstellar discs, and to find systematic trends and correlations between the parameters, the continuum fluxes, and the line fluxes. The model grid, comprising the calculated disc temperatures and chemical structures, the computed SEDs, line fluxes and profiles, will be used in particular for the data interpretation of the Herschel open time key programme GASPS. The calculated line fluxes show a strong dependence on the assumed UV excess of the central star, and on the disc flaring. The fraction of models predicting [OI] and [CII] fine-structure lines fluxes above Herschel/PACS and Spica/SAFARI detection limits are calculated as function of disc mass. The possibility of deriving the disc gas mass from line observations is discussed.

We present VLT/SINFONI integral field spectroscopy of RCW 34 along with Spitzer/IRAC photometry of the surroundings. RCW 34 consists of three different regions. A large bubble has been detected on the IRAC images in which a cluster of intermediate- and low-mass class II objects is found. At the northern edge of this bubble, an HII region is located, ionized by 3 OB stars. Intermediate mass stars (2 - 3 Msun) are detected of G- and K- spectral type. These stars are still in the pre-main sequence (PMS) phase. North of the HII region, a photon-dominated region is present, marking the edge of a dense molecular cloud traced by H2 emission. Several class 0/I objects are associated with this cloud, indicating that star formation is still taking place. The distance to RCW 34 is revised to 2.5 +- 0.2 kpc and an age estimate of 2 - 1 Myrs is derived from the properties of the PMS stars inside the HII region. The most likely scenario for the formation of the three regions is that star formation propagates from South to North. First the bubble is formed, produced by intermediate- and low-mass stars only, after that, the HII region is formed from a dense core at the edge of the molecular cloud, resulting in the expansion as a champagne flow. More recently, star formation occurred in the rest of the molecular cloud. Two different formation scenarios are possible: (a) The bubble with the cluster of low- and intermediate mass stars triggered the formation of the O star at the edge of the molecular cloud which in turn induces the current star-formation in the molecular cloud. (b) An external triggering is responsible for the star-formation propagating from South to North. [abridged]

We present a R=10,000 M-band spectrum of LLN19 (IRAS 08470-4321), a heavily embedded intermediate-mass young stellar object located in the Vela Molecular Cloud, obtained with VLT-ISAAC. The data were fitted by a 2-slab cold-hot model and a wind model. The spectrum exhibits deep broad ro-vibrational absorption lines of 12CO v=1<-0 and 13CO v=1<-0. A weak CO ice feature at 4.67 micron is also detected. Differences in velocity indicate that the warm gas is distinct from the cold millimeter emitting gas, which may be associated with the absorption by cooler gas (45K). The outflowing warm gas at 300-400K and with a mass-loss rate varying between 0.48E-7 and 4.2E-7 MSun /yr can explain most of the absorption. Several absorption lines were spectrally resolved in subsequent spectra obtained with the VLT-CRIRES instrument. Multiple absorption substructures in the high-resolution (R=100,000) spectra indicate that the mass-loss is episodic with at least two major events that occurred recently (<28 years). The discrete mass-loss events together with the large turbulent width of the gas (dv=10-12 km/s) are consistent with the predictions of the Jet-Bow shock outflow and the wide-angle wind model. The CO gas/solid column density ratio of 20-100 in the line-of-sight confirms that the circumstellar environment of LLN~19 is warm. We also derive a 12C/13C ratio of 67 +/- 3, consistent with previous measurements in local molecular clouds but not with the higher ratios found in the envelope of other young stellar objects.

The origin of Earth oceans is controversial. Earth could have acquired its water either from hydrated silicates (wet Earth scenario) or from comets (dry Earth scenario). [HDO]/[H2O] ratios are used to discriminate between the scenarios. High [HDO]/[H2O] ratios are found in Earth oceans. These high ratios are often attributed to the release of deuterium enriched cometary water ice, which was formed at low gas and dust temperatures. Observations do not show high [HDO]/[H2O] in interstellar ices. We investigate the possible formation of high [HDO]/[H2O] ratios in dense (nH> 1E6 cm^{-3}) and warm gas (T=100-1000 K) by gas-phase photochemistry in the absence of grain surface chemistry. We derive analytical solutions, taking into account the major neutral-neutral reactions for gases at T>100 K. The chemical network is dominated by photodissociation and neutral-neutral reactions. Despite the high gas temperature, deuterium fractionation occurs because of the difference in activation energy between deuteration enrichment and the back reactions. The analytical solutions were confirmed by the time-dependent chemical results in a 1E-3 MSun disc around a typical TTauri star using the photochemical code ProDiMo. The ProDiMo code includes frequency-dependent 2D dust-continuum radiative transfer, detailed non-LTE gas heating and cooling, and hydrostatic calculation of the disc structure. Both analytical and time-dependent models predict high [HDO]/[H2O] ratios in the terrestrial planet forming region (< 3 AU) of circumstellar discs. Therefore the [HDO]/[H2O] ratio may not be an unique criterion to discriminate between the different origins of water on Earth.

Methods. We use the recently developed disk code ProDiMo to calculate the physico-chemical structure of protoplanetary disks and apply the Monte-Carlo line radiative transfer code RATRAN to predict observable line profiles and fluxes. We consider a series of Herbig Ae type disk models ranging from 10^-6 M_Sun to 2.2 10^-2 M_Sun (between 0.5 and 700 AU) to discuss the dependency of the line fluxes and ratios on disk mass for otherwise fixed disk parameters. Results. We find the [CII] 157.7 mum line to originate in LTE from the surface layers of the disk, where Tg > Td . The total emission is dominated by surface area and hence depends strongly on disk outer radius. The [OI] lines can be very bright (> 10^-16 W/m^2) and form in slightly deeper and closer regions under non-LTE conditions. The high-excitation [OI] 145.5 mum line, which has a larger critical density, decreases more rapidly with disk mass than the 63.2 mum line. Therefore, the [OI] 63.2 mum/145.5 mum ratio is a promising disk mass indicator, especially as it is independent of disk outer radius for Rout > 200 AU. CO is abundant only in deeper layers A_V >~ 0.05. For too low disk masses (M_disk <~10^-4 M_Sun) the dust starts to become transparent, and CO is almost completely photo-dissociated. For masses larger than that the lines are an excellent independent tracer of disk outer radius and can break the outer radius degeneracy in the [OI] 63.2 mum/[CII]157.7 mum line ratio. Conclusions. The far-IR fine-structure lines of [CII] and [OI] observable with Herschel provide a promising tool to measure the disk gas mass, although they are mainly generated in the atomic surface layers. In spatially unresolved observations, none of these lines carry much information about the inner, possibly hot regions < 30 AU.

The spatial origin and detectability of rotational H2O emission lines from Herbig Ae type protoplanetary disks beyond 70 micron is discussed. We use the recently developed disk code ProDiMo to calculate the thermo-chemical structure of a Herbig Ae type disk and apply the non-LTE line radiative transfer code Ratran to predict water line profiles and intensity maps. The model shows three spatially distinct regions in the disk where water concentrations are high, related to different chemical pathways to form the water: (1) a big water reservoir in the deep midplane behind the inner rim, (2) a belt of cold water around the distant icy midplane beyond the snow-line r>20AU, and (3) a layer of irradiated hot water at high altitudes z/r=0.1...0.3, extending from about 1AU to 30AU, where the kinetic gas temperature ranges from 200K to 1500K. Although region 3 contains only little amounts of water vapour (~3x10^-5 M_Earth), we find this warm layer to be almost entirely responsible for the rotational water emission lines, execpt for the 3 lowest excitation lines. Thus, Herschel will probe first and foremost the conditions and radial extension of region 3, where water is predominantly formed via neutral-neutral reactions and the gas is thermally decoupled from the dust T_gas>T_dust. The observations do not allow for a determination of the snow-line, because the snow-line truncates the radial extension of region 1, whereas the lines originate from region 3. Different line transfer approximations (LTE, escape probability, Monte Carlo) are discussed. A non-LTE treatment is required in most cases, and the results obtained with the escape probability method are found to underestimate the Monte Carlo results by 2%...45%.

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.

This article summarizes a Splinter Session at the Cool Stars XV conference in St. Andrews with 3 review and 4 contributed talks. The speakers have discussed various approaches to understand the structure and evolution of the gas component in protoplanetary disks. These ranged from observational spectroscopy in the UV, infrared and millimeter, through to chemical and hydrodynamical models. The focus was on disks around low-mass stars, ranging from classical T Tauri stars to transitional disks and debris disks. Emphasis was put on water and organic molecules, the relation to planet formation, and the formation of holes and gaps in the inner regions.

We aim at understanding the formation of cloud layers in quasi-static substellar atmospheres. The time-dependent description presented in (Helling & Woitke 2006) is a kinetic model describing nucleation, growth and evaporation. It is extended to treat gravitational settling and is applied to the static-stationary case of substellar model atmospheres. From the solution for the dust moments, we determine the grain size distribution function which, together with the calculated material volume fractions, provides the basis to calculate the opacities of the composite dust grains. The cloud particles in brown dwarfs and hot giant-gas planets are found to be small in the high atmospheric layers (0.01mum), and composed of a rich mixture of all considered condensates, in particular the abundant MgSiO3[s], Mg2SiO4[s] and SiO2[s]. As the particles settle downward, they increase in size and reach several 100mum in the deepest layers. The more volatile parts of the grains evaporate and the particles stepwise purify to form composite particles of high-temperature condensates in the deeper layers, mainly Fe[s] and Al2O3[s]. The gas phase abundances of the elements involved in the dust formation process vary by orders of magnitudes throughout the atmosphere. The grain size distribution is found to be relatively broad in the upper atmospheric layers but often strongly peaked in the deeper layers. The spectral appearance of the cloud layers in the mid IR (7-20mum) is close to a grey body with only weak broad features on a few percent level, mainly caused by MgSiO3[s], and Mg2SiO4[s]. Our models predict that the gas phase depletion is much weaker as compared to phase-equilibrium calculations in the high atmospheric layers. [abridged]

Near-infrared imaging surveys of high-mass star-forming regions reveal an amazingly complex interplay between star formation and the environment (Churchwell et al. 2006; Alvarez et al. 2004). By means of near-IR spectroscopy the embedded massive young stars can be characterized and placed in the context of their birth site. However, so far spectroscopic surveys have been hopelessly incomplete, hampering any systematic study of these very young massive stars. New integral field instrumentation available at ESO has opened the possibility to take a huge step forward by obtaining a full spectral inventory of the youngest massive stellar populations in star-forming regions currently accessible. Simultaneously, the analysis of the extended emission allows the characterization of the environmental conditions. The Formation and Early Evolution of Massive Stars (FEMS) collaboration aims at setting up a large observing campaign to obtain a full census of the stellar content, ionized material, outflows and PDR's over a sample of regions that covers a large parameter space. Complementary radio, mm and infrared observations will be used for the characterization of the deeply embedded population. For the first eight regions we have obtained 40 hours of SINFONI observations. In this contribution, we present the first results on three regions that illustrate the potential of this strategy.

We observed the Herbig Ae/Be stars UX Ori, HD 34282, HD 100453, HD 101412, HD 104237 and HD 142666, and the T Tauri star HD 319139 and searched for H2 0-0 S(2) emission at 12.278 micron and H2 0-0 S(1) emission at 17.035 micron with VISIR, ESO-VLT's high-resolution MIR spectrograph. None of the sources present evidence for H2 emission. Stringent 3sigma upper limits to the integrated line fluxes and the mass of optically thin warm gas in the disks are derived. The disks contain less than a few tenths of Jupiter mass of optically thin H2 gas at 150 K at most, and less than a few Earth masses of optically thin H2 gas at 300 K and higher temperatures. We compare our results to a Chiang and Goldreich (1997, CG97) two-layer disk model. The upper limits to the disk's optically thin warm gas mass are smaller than the amount of warm gas in the interior layer of the disk, but they are much larger than the amount of molecular gas in the surface layer. We present a calculation of the expected thermal H2 emission from optically thick disks, assuming a CG97 disk structure, a gas-to-dust ratio of 100 and Tgas = Tdust. The expected H2 thermal emission fluxes from typical disks around Herbig Ae/Be stars (10^-16 to 10^-17 erg/s/cm2 at 140 pc) are much lower than the detection limits of our observations (5*10^-15 erg/s/cm2). H2 emission levels are very sensitive to departures from the thermal coupling between the molecular gas and dust. Additional sources of heating of gas in the disk's surface layer could have a major impact on the expected H2 disk emission. In the observed sources the molecular gas and dust in the surface layer have not significantly departed from thermal coupling (Tgas/Tdust< 2) and that the gas-to-dust ratio in the surface layer is very likely lower than 1000.

Near-infrared surveys of high-mass star-forming regions start to shed light onto their stellar content. A particular class of objects found in these regions, the so-called massive Young Stellar Objects (YSOs) are surrounded by dense circumstellar material. Several near- and mid-infrared diagnostic tools are used to infer the physical characteristics and geometry of this circumstellar matter. Near-infrared hydrogen emission lines provide evidence for a disk-wind. The profiles of the first overtone of the CO band-heads, originating in the inner 10 AU from the central star, are well fitted assuming a keplerian rotating disk. The mid-infrared spectral energy distribution requires the presence of a more extended envelope containing dust at a temperature of about 200 K. CRIRES observations of CO fundamental absorption lines confirm the presence of a cold envelope. We discuss the evolutionary status of these objects.

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

Dust clouds influence the atmospheric structure of brown dwarfs, and they affect the heat transfer and change the gas-phase chemistry. However, the physics of their formation and evolution is not well understood. In this letter, we predict dust signatures and propose a potential observational test of the physics of dust formation in brown dwarf atmosphere based on the spectral features of the different solid components predicted by dust formation theory. A momentum method for the formation of dirty dust grains (nucleation, growth, evaporation, drift) is used in application to a static brown dwarf atmosphere structure to compute the dust grain properties, in particular the heterogeneous grain composition and the grain size. Effective medium and Mie theory are used to compute the extinction of these spherical grains. Dust formation results in grains whose composition differs from that of grains formed at equilibrium. Our kinetic model predicts that solid amorphous SiO2[s] (silica) is one of the most abundant solid component followed by amorphous MgSiO4[s] and MgSiO3[s], while SiO2[s] is absent in equilibrium models because it is a metastable solid. Solid amorphous SiO2[s] possesses a strong broad absorption feature centered at 8.7mum, while amorphous Mg2SiO4[s]/MgSiO3[s] absorb at 9.7mum beside other absorption features at longer wavelength. Those features at lambda < 15mum are detectable in absorption if grains are small (radius < 0.2mum) in the upper atmosphere as suggested by our model. We suggest that the detection of a feature at 8.7mum in deep infrared spectra could provide evidence for non-equilibrium dust formation that yields grains composed of metastable solids in brown dwarf atmospheres. This feature will shift towards 10mum and broaden if silicates (e.g. fosterite) are much more abundant.

We performed a spectroscopic survey toward five intermediate-mass class I YSOs located in the Southern Vela molecular cloud in the L and M bands at resolving powers 600-800 up to 10,000, using the Infrared Spectrometer and Array Camera mounted on the VLT-ANTU. Lower mass companion objects were observed simultaneously in both bands. Solid H2O at 3 micron is detected in all sources, including the companion objects. CO ice at 4.67 micron is detected in a few main targets and one companion object. One object (LLN 19) shows little CO ice but strong gas-phase CO ro-vibrational lines in absorption. The CO ice profiles are different from source to source. The amount of water ice and CO ice trapped in a water-rich mantle may correlate with the flux ratio at 12 and 25 micron. The abundance of H2O-rich CO likely correlates with that of water ice. A weak feature at 3.54 mu attributed to solid CH3OH and a broad feature near 4.62 mu are observed toward LLN17, but not toward the other sources. The derived abundances of solid CH3OH and OCN- are ~10% and ~1% of the H2O ice abundance respectively. The H2O optical depths do not show an increase with envelope mass, nor do they show lower values for the companion objects compared with the main protostar. The line-of-sight CO ice abundance does not correlate with the source bolometric luminosity. Comparison of the solid CO profile toward LLN17, which shows an extremely broad CO ice feature, and that of its lower mass companion at a few thousand AU, which exhibits a narrow profile, together with the detection of OCN- toward LLN17 provide direct evidences for local thermal processing of the ice.

High-quality K-band spectra of point sources, deeply embedded in massive star-forming regions, have revealed a population of 20 young massive stars showing no photospheric absorption lines, but only emission lines. The K-band spectra exhibit one or more features commonly associated with massive Young Stellar Objects surrounded by circumstellar material: a very red color (J-K) = 2, CO bandhead emission, hydrogen emission lines (sometimes doubly peaked), and FeII and/or MgII emission lines. The CO emission comes from a relatively dense (~10^10 cm^(-3)) and hot (T ~ 2000-5000 K) region, sufficiently shielded from the intense UV radiation field of the young massive star. Modeling of the CO-first overtone emission shows that the CO gas is located within 5 AU of the star. The hydrogen emission is produced in an ionized medium exposed to UV radiation. The best geometrical configuration is a dense and neutral circumstellar disk causing the CO bandhead emission, and an ionized upper layer where the hydrogen lines are produced. We argue that the circumstellar disk is likely a remnant of the accretion via a circumstellar disk.

We report on the observation of hot water vapor (steam) in the inner AU of a young massive star located in the star-forming region IRAS 08576-4334. The water lines are detected in a medium resolution (R=10,000) K-band spectrum taken by the infrared spectrometer ISAAC mounted on the VLT-ANTU. The water vapor is at a mean temperature of 1565+/-510 K, cooler than the hot CO found in the same object, which is at 1660K and the column density is N(H2O)=(2.5 +/- 0.4)x1E18 cm-2. The profile of both H2O and CO lines is best reproduced by the emission from a Keplerian disk. To interpret the data, we also investigate the formation of molecules and especially CO and water vapor in the inner hot and dense part of disks around young high mass stars using a pseudo time-dependent gas-phase chemical model. Molecules are rapidly photodissociated but this destruction is compensated by an efficient formation due to fast neutral-neutral reactions. The ability of CO molecules to self-shield significantly enhances its abundance. Water molecules are sufficiently abundant to be detectable. The observed H2O/CO ratio is reproduced by gas at 1600K and an enhanced UV field over gas density ratio I_UV/nH=1E(-4)- 1E(-6). The simulations support the presence of CO and H2O molecules in the inner disks around young massive stars despite the strong UV radiation and show that the OH radical plays an essential role in hot gas chemistry.

We present the results of medium-resolution spectroscopy of five nearby Herbig Ae/Be stars at 4.7 micron: UX Ori, HD 34282, HD 50138, V380 Ori, HK Ori. The goal was to search for CO fundamental ro-vibrational emission. None of the targets show CO features, either in absorption nor in emission. We derive a 5sigma upper limit of < 10^-12 cm^-2 to the column density of hot CO (T ~ 1500 K) in the sources. These upper limits are considerably lower than the values of Herbig Ae/Be stars for which warm and hot CO emission has been reported.The non-detection of CO v=1-0 emission in these five targets suggest that Herbig Ae/Be stars are not a homogeneous group with respect to the structure of the gaseous disk and/or the amount of CO in the inner 50 AU of their disks.

We report on the observation of CO bandhead emission around 51 Oph (dv=2). A high resolving power (R~10,000) spectrum was obtained with the infrared spectrometer ISAAC mounted on VLT-ANTU. Modeling of the profile suggests that the hot (Tgas=2000-4000K) and dense (n_H>1e10 cm-3) molecular material as probed by the CO bandhead is located in the inner AU of a Keplerian disk viewed almost edge-on. Combined with the observation of cooler gas (Tgas=500-900K) by ISO-SWS and the lack of cold material, our data suggest that the disk around 51 Oph is essentially warm and small. We demonstrate the presence of a dust-free inner disk that extents from the inner truncation radius until the dust sublimation radius. The disk around 51 Oph may be in a rare transition state toward a small debris disk object.

We present observations of CO overtone bandhead emission toward four massive Young Stellar Objects (spectral type O6--B5). The high signal-to-noise ratio K-band spectra were obtained with VLT-ISAAC at a resolution of 30 km/s, sufficient to resolve the bandheads, but not the individual J-lines. We are able to explain the shape of the lines by assuming a simple isothermal keplerian disk model seen at different inclinations. The gas temperature ranges from 1500 to 4500 K and the CO column density is between 0.1 and 4 10^21 cm^-2. The emission probably arises within the first few astronomical units of the disk, consistent with the high gas temperature. Our results indicate that molecules can survive close to a hot star and suggest that dense (n_H > 10^10 cm^-3) inner disks may be relatively common at an advanced stage of high-mass star formation.

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 present 3-5 micron VLT-ISAAC spectroscopy searching for evidence of methanol ices in edge-on disks of young embedded stars. Examples include the disks of L1489 IRS in Taurus and CRBR 2422.8-3423 in Ophiuchus, the last of which has the highest column density of solid CO known toward a YSO. We find no unambiguous evidence for abundant methanol in the observed disks, but give strict upper limits. Several additional low-mass sources in the Serpens and Chameleon molecular clouds exhibit abundant solid methanol, although it is not clear if the ice is associated with a disk or with the envelope. These are the first detections of solid methanol in the circumstellar environments of embedded young low-mass stars providing evidence that complex molecular species previously observed only in the solid state toward high-mass star forming regions are also present near solar-type young stars. The constraints on the formation mechanisms of methanol and the chemical evolution of ices as the material is incorporated into circumstellar disks are discussed.

Medium resolution (lambda/Delta lambda = 5000-10000) VLT-ISAAC M-band spectra are presented of 39 young stellar objects in nearby low-mass star forming clouds showing the 4.67 micron stretching vibration mode of solid CO. By taking advantage of the unprecedentedly large sample, high S/N ratio and high spectral resolution, similarities in the ice profiles from source to source are identified. It is found that excellent fits to all the spectra can be obtained using a phenomenological decomposition of the CO stretching vibration profile at 4.67 micron into 3 components. All observed interstellar CO profiles can thus be uniquely described by a model depending on only 3 linear fit parameters, indicating that a maximum of 3 specific molecular environments of solid CO exist under astrophysical conditions. A simple physical model of the CO ice is presented, which shows that the 2139.9 cm-1 component is indistinguishable from pure CO ice. It is concluded, that in the majority of the observed lines of sight, 60-90% of the CO is in a nearly pure form. In the same model the 2143.7 cm-1 component can possibly be explained by the longitudinal optical (LO) component of the vibrational transition in pure crystalline CO ice which appears when the background source is linearly polarised. Furthermore the absorption band from solid 13CO at 2092 cm-1 is detected towards IRS 51 in the Rho Ophiuchi cloud complex and an isotopic ratio of 12CO/13CO=68+/-10 is derived. It is shown that all the observed solid 12CO profiles, along with the solid 13CO profile, are consistent with grains with an irregularly shaped CO ice mantle simulated by a Continuous Distribution of Ellipsoids (CDE), but inconsistent with the commonly used models of spherical grains in the Rayleigh limit. [Abstract abridged]

We present detections of the absorption band at 3.53 micron due to solid methanol toward three low-mass young stellar objects located in the Serpens and Chameleon molecular cloud complexes. The sources were observed as part of a large spectroscopic survey of ~40 protostars. This is the first detection of solid methanol in the vicinity of low mass (M <1 Msol) young stars and shows that the formation of methanol does not depend on the proximity of massive young stars. The abundances of solid methanol compared to water ice for the three sources are in the range 15-25% which is comparable to those for the most methanol-rich massive sources known. The presence of abundant methanol in the circumstellar environment of some low mass young stars has important consequences for the formation scenarios of methanol and more complex organic species near young solar-type stars.

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.

General results from a 3-5 micron spectroscopic survey of nearby low-mass young stellar objects are presented. L and M-band spectra have been obtained of \~50 low mass embedded young stars using the ISAAC spectrometer mounted on UT1-Antu at Paranal Observatory. For the first time, a consistent census of the CO, H2O ices and the minor ice species CH3OH and OCN- and warm CO gas present around young stars is obtained, using large number statistics and resolving powers of up to R=10000. The molecular structure of circumstellar CO ices, the depletion of gaseous CO onto grains in protoplanetary disks, the presence of hot gas in the inner parts of circumstellar disks and in outflows and infalls are studied. Furthermore, the importance of scattering effects for the interpretation of the spectra have been addressed.

We present direct evidence for CO freeze-out in a circumstellar disk around the edge-on class I object CRBR 2422.8-3423, observed in the M band with VLT-ISAAC at a resolving power R~10,000. The spectrum shows strong solid CO absorption, with a lower limit on the column density of 2.2E18 cm-2. The solid CO column is the highest observed so far, including high-mass protostars and background field stars. Absorption by foreground cloud material likely accounts for less than 10% percent of the total solid CO, based on the weakness of solid CO absorption toward nearby sources and the absence of gaseous C18O J=2-1 emission 30'' south. Gas-phase ro-vibrational CO absorption lines are also detected with a mean temperature of 50 +/-10 K. The average gas/solid CO ratio is ~1 along the line of sight. For an estimated inclination of 20 +/- 5 degree, the solid CO absorption originates mostly in the cold, shielded outer part of the flaring disk, consistent with the predominance of apolar solid CO in the spectrum and the non-detection of solid OCN-, a thermal/ultraviolet processing of the ice mantle. The gaseous CO comes from the warm upper layers closer to the star.

There is now strong evidence suggesting that the 12CO J = 1-0 transition, widely used to trace H2 gas, significantly underestimates its mass in metal-poor regions. In spiral disks such regions are found in large galactocentric distances where we show that any unaccounted H2 gas phase is likely to be diffuse (~5-20 cm^-3) and warmer (T(kin) ~ 50-100 K) than the cool (T(kin) ~ 15-20 K) CO-luminous one. Moreover we find that a high value of the H2 formation rate on grains, suggested by recent observational work, can compensate for the reduction of the available grain surface in the metal-poor part of typical galactic disks and thus enhance this CO-poor H2 component which may be contributing significantly to the mass and pressure of spiral disks beyond their optical radius.

We have carried out an in-depth study of the peripheral region of the molecular cloud L1204/S140, where the far ultraviolet radiation and the density are relatively low. Our observations test theories of photon-dominated regions (PDRs) in a regime that has been little explored. CII 158 micron and OI 63 micron lines are detected by ISO at all 16 positions along a 1-dimensional cut in right ascension. Emission from molecular hydrogen rotational transitions, at 28 and 17 micron, was also detected at several positions. The CII, OI, and molecular hydrogen intensities along the cut show much less spatial variation than do the rotational lines of CO and other CO isotopes. The average CII and OI intensities and their ratio are consistent with models of PDRs with low FUV radiation (Go) and density. The best-fitting model has Go about 15 and density, n about 1000 per cubic cm. Standard PDR models underpredict the intensity in the H2 rotational lines by up to an order of magnitude. This problem has also been seen in bright PDRs and attributed to factors, such as geometry and gas-grain drift, that should be much less important in the regime studied here. The fact that we see the same problem in our data suggests that more fundamental solutions, such as higher H2 formation rates, are needed. Also, in this regime of low density and small line width, the OI line is sensitive to the radiative transfer and geometry. Using the ionization structure of the models, a quantitative analysis of timescales for ambipolar diffusion in the peripheral regions of the S140 cloud is consistent with a theory of photoionization-regulated star formation. Observations of CII in other galaxies differ both from those of high Go PDRs in our galaxy and from the low Go regions we have studied. ~

The first results from a VLT-ISAAC program on L- and M-band infrared spectroscopy of deeply-embedded young stellar objects are presented. The advent of 8-m class telescopes allows high S/N spectra of low-luminosity sources to be obtained. In our first observing run, low- and medium-resolution spectra have been measured toward a dozen objects, mostly in the Vela and Chamaeleon molecular clouds. The spectra show strong absorption of H2O and CO ice, as well as weak features at `3.47' and 4.62 mu. No significant solid CH3OH feature at 3.54 mu is found, indicating that the CH3OH/H2O ice abundance is lower than toward some massive protostars. Various evolutionary diagnostics are investigated for a set of sources in Vela.

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.

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.