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.

Little is known about the interaction between atmospheres and crusts of exoplanets so far, but future space missions and ground-based instruments are expected to detect molecular features in the spectra of hot rocky exoplanets. We aim to understand the composition of the gas in an exoplanet atmosphere which is in equilibrium with a planetary crust. Methods. The molecular composition of the gas above a surface made of a mixture of solid and liquid materials is determined by assuming phase equilibrium for given pressure, temperature and element abundances. We study total element abundances that represent different parts of the Earth's crust (Continental Crust, Bulk Silicate Earth, Mid Oceanic Ridge Basalt), CI chondrites and abundances measured in polluted white dwarfs. For temperatures between 600K and 3500K, the near-crust atmospheres of all considered total element abundances are mainly composed of H2O, CO2, SO2 and in some cases of O2 and H2. For temperatures 500K, only N2-rich or CH4-rich atmospheres remain. For 3500K, the atmospheric gas is mainly composed of atoms (O, Na, Mg, Fe), metal oxides (SiO, NaO, MgO, CaO, AlO, FeO) and some metal hydroxides (KOH, NaOH). The inclusion of phyllosilicates as potential condensed species is crucial for lower temperatures, as they can remove water from the gas phase below about 700K and inhibit the presence of liquid water. Measurements of the atmospheric composition could, in principle, characterise the rock composition of exoplanet crusts. H2O, O2 and CH4 are natural products from the outgassing of different kinds of rocks that had time to equilibrate. These are discussed as biomarkers, but do emerge naturally as result of the thermodynamic interaction between crust and atmosphere. Only the simultaneous detection of all three molecules might be a sufficient biosignature, as it is inconsistent with chemical equilibrium.

The atmospheres of ultra-hot Jupiters are commonly considered to be at thermochemical equilibrium. We aim to provide disequilibrium chemistry maps for a global understanding of the chemistry in HAT-P-7b's atmosphere and assess the importance of disequilibrium chemistry on UHJs. We apply a hierarchical modelling approach utilising 97 1D atmospheric profiles from 3D GCM of HAT-P-7b. For each 1D profile, we evaluate our kinetic cloud formation model consistently with the local gas-phase composition in chemical equilibrium. We then evaluate quenching results from a zeroth-order approximation in comparison to a kinetic gas-phase approach. We find that the zeroth-order approach of estimating quenching points agrees well with the full gas-kinetic modeling results. Chemical disequilibrium has the greatest effect on the nightside and morning abundance of species such as H, H$_2$O, CH$_4$, CO$_2$, HCN, and all C$_n$H$_m$ molecules; heavier C$_n$H$_m$ molecules are more affected by disequilibrium processes. CO abundance, however, is affected only marginally. While dayside abundances also notably change, those around the evening terminator of HAT-P-7b are the least affected by disequilibrium processes. The latter finding may partially explain the consistency of observed transmission spectra of UHJs with atmospheres in thermochemical equilibrium. Photochemistry only negligibly affects molecular abundances and quenching levels. In general, the quenching points of HAT-P-7b's atmosphere are at much lower pressures in comparison to the cooler hot-jupiters. We propose several avenues to look for the effect of disequilibrium processes on UHJs that are, in general, based on abundance and opacity measurements at different local times. It remains a challenge to completely disentangle this from the chemical effects of clouds and that of a primordial non-solar abundance.

The precipitation of cloud particles in brown dwarf and exoplanet atmospheres establishes an ongoing downward flux of condensable elements. To understand the efficiency of cloud formation, it is therefore crucial to quantify the replenishment mechanism that is able to compensate for the element losses in the upper atmosphere, to keep the extrasolar weather cycle running. In this paper, we introduce a new cloud formation model by combining the cloud particle moment method of Helling & Woitke with a diffusive mixing approach, taking into account turbulent mixing and gas-kinetic diffusion for both gas and cloud particles. The equations are of diffusion-reaction type and are solved time-dependently for a prescribed 1D atmospheric structure, until the model has relaxed toward a time-independent solution. In comparison to our previous models, the new hot Jupiter model results (Teff=2000K, log g=3) show fewer but larger cloud particles which are more concentrated towards the cloud base. The abundances of condensable elements in the gas phase are featured by a steep decline above the cloud base, followed by a shallower, monotonous decrease towards a plateau, the level of which depends on temperature. The chemical composition of the cloud particles also differs significantly from our previous models. Due to the condensation of specific condensates like Mg2SiO4[s] in deeper layers, certain elements, such as Mg, are almost entirely removed from the gas phase early. This leads to unusual (and non-solar) element ratios in higher atmospheric layers, which then favours the formation of SiO[s] and SiO2[s], for example, rather than MgSiO3[s]. Such condensates are not expected in phase-equilibrium models that start from solar abundances. Above the main silicate cloud layer, which is enriched with iron and metal oxides, we find a second cloud layer made of Na2S[s] particles in cooler models (Teff<1500K).

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.

Ultra-hot Jupiters have recently attracted interest from observers and theoreticians alike, as they provide observationally accessible test cases. We apply a hierarchical modelling approach as a virtual laboratory to study cloud formation and gas-phase chemistry. We utilise 97 vertical 1D profiles of a 3D GCM for HAT-P-7b to evaluate our kinetic cloud formation model consistently with the local equilibrium gas-phase composition. The day/night temperature difference on HAT-P-7b (~ 2500K) causes clouds to form on the nightside (dominated by H2/He) while the dayside (dominated by H/He) retains cloud-free equatorial regions. The cloud particles vary in composition and size throughout the vertical extension of the cloud, but also globally. TiO2[s]/Al2O3[s]/CaTiO3[s]-particles of cm-sized radii occur in the higher dayside-latitudes, resulting in a dayside dominated by gas-phase opacity. The opacity on the nightside, however, is dominated by 0.01 ... 0.1 mum particles made of a material mix dominated by silicates. The gas pressure at which the atmosphere becomes optically thick is ~1d-4 bar in cloudy regions, and ~0.1 bar in cloud-free regions. HAT-P-7b features strong morning/evening terminator asymmetries, providing an example of patchy clouds and azimuthally-inhomogeneous chemistry. The large temperature differences result in an increasing geometrical extension from the night- to the dayside. The chemcial equilibrium H2O abundance at the terminator changes by < 1 dex with altitude and < 0.3 dex (a factor of 2) across the terminator for a given pressure, indicating that H2O abundances derived from transmission spectra can be representative of the well-mixed metallicity at P > 10 bar. We suggest the atmospheric C/O as an important tool to trace the presence and location of clouds in exoplanet atmospheres. Phase curve variability of HAT-P-7b is unlikely to be caused by dayside clouds.

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.

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.

WASP-18b is an utra-hot Jupiter with a temperature difference of upto 2500K between day and night. Such giant planets begin to emerge as planetary laboratory for understanding cloud formation and gas chemistry in well-tested parameter regimes in order to better understand planetary mass loss and for linking observed element ratios to planet formation and evolution. We aim to understand where clouds form, their interaction with the gas phase chemistry through depletion and enrichment, the ionisation of the atmospheric gas and the possible emergence of an ionosphere on ultra-hot Jupiters. We utilize 1D profiles from a 3D atmosphere simulations for WASP-18b as input for kinetic cloud formation and gas-phase chemical equilibrium calculations. We solve our kinetic cloud formation model for these 1D profiles that sample the atmosphere of WASP-18b at 16 different locations along the equator and in the mid-latitudes and derive consistently the gas-phase composition. The dayside of WASP-18b emerges as completely cloud-free due to the very high atmospheric temperatures. In contrast, the nightside is covered in geometrically extended and chemically heterogeneous clouds with disperse particle size distributions. The atmospheric C/O increases to $>0.7$ and the enrichment of the atmospheric gas with cloud particles is $\rho_{\rm d}/\rho_{\rm gas}>10^{-3}$. The clouds that form at the limbs appear located farther inside the atmosphere and they are the least extended. Not all day-night terminator regions form clouds. The gas-phase is dominated by H$_2$, CO, SiO, H$_2$O, H$_2$S, CH$_4$, SiS. In addition, the dayside has a substantial degree of ionisation due to ions like Na$^+$, K$^+$, Ca$^+$, Fe$^+$. Al$^+$ and Ti$^+$ are the most abundant of their element classes. We find that WASP-18b, as one example for ultra-hot Jupiters, develops an ionosphere on the dayside.

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.

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.

We introduce a fast and versatile computer code, GGchem, to determine the chemical composition of gases in thermo-chemical equilibrium down to 100 K, with or without equilibrium condensation. We review the data for molecular equilibrium constants, kp(T), from several sources and discuss which functional fits are most suitable for low temperatures. We benchmark our results against another chemical equilibrium code. We collect Gibbs free energies, dG(T), for about 200 solid and liquid species from the NIST-JANAF database and the geophysical database SUPCRTBL. We discuss the condensation sequence of the elements with solar abundances in phase equilibrium down to 100 K. Once the major magnesium silicates Mg2SiO4[s] and MgSiO3[s] have formed, the dust/gas mass ratio jumps to a value of about 0.0045 which is significantly lower than the often assumed value of 0.01. Silicate condensation is found to increase the carbon/oxygen ratio (C/O) in the gas from its solar value of ~0.55 up to ~0.71, and, by the additional intake of water and hydroxyl into the solid matrix, the formation of phyllosilicates at temperatures below ~400 K increases the gaseous C/O further to about 0.83. Metallic tungsten (W) is the first condensate found to become thermodynamically stable around 1600 - 2200 K (depending on pressure), several hundreds of Kelvin before subsequent materials like zirconium dioxide (ZrO2) or corundum (Al2O3) can condense. We briefly discuss whether tungsten, despite its low abundance of ~2.E-7 times the silicon abundance, could provide the first seed particles for astrophysical dust formation. The GGchem code is publicly available at https://github.com/pw31/GGchem.

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.

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.

Recent observations indicate potentially carbon-rich exoplanet atmospheres. Spectral fitting methods for brown dwarfs and exoplanets have invoked the C/O ratio as additional parameter but carbon-rich cloud formation modeling is a challenge for the models applied. The determination of the habitable zone for exoplanets requires the treatment of cloud formation in chemically different regimes. Disk models show that carbon-rich or near-carbon-rich niches may emerge and cool carbon planets may trace these particular stages of planetary evolution. We extend our kinetic cloud formation model by including carbon seed formation and the formation of C[s], TiC[s], SiC[s], KCl[s], and MgS[s] by gas-surface reactions. The seed formation efficiency is lower in carbon-rich atmospheres than in oxygen-rich gases due to carbon being a very effective growth species. The consequence is that less particles will make up a cloud for C/O_0>1. The cloud particles will be smaller in size than in an oxygen-rich atmosphere. An increasing initial C/O ratio does not revert this trend because a much greater abundance of condensible gas species exists in a rich carbon environment. Cloud particles are generally made of a mix of materials: carbon dominates if C/O_0>1 and silicates dominate if C/O_0<1. 80-90% carbon is reached only in extreme cases where C/O_0=3.0 or 10.0. Carbon-rich atmospheres would form clouds that are made of particles of height-dependent mixed compositions, sizes and numbers. The remaining gas-phase is far less depleted than in an oxygen-rich atmosphere. Typical tracer molecules are HCN and C2H2 in combination with a featureless, smooth continuum due to a carbonaceous cloud cover, unless the cloud particles become crystalline.

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.

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

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.

Observations of exoplanet atmospheres have revealed the presence of cloud particles in their atmospheres. 3D modelling of cloud formation in atmospheres of extrasolar planets coupled to the atmospheric dynamics has long been a challenge. We investigate the thermo-hydrodynamic properties of cloud formation processes in the atmospheres of hot Jupiter exoplanets. We simulate the dynamic atmosphere of HD 189733b with a 3D model that couples 3D radiative-hydrodynamics with a kinetic, microphysical mineral cloud formation module designed for RHD/GCM exoplanet atmosphere simulations. Our simulation includes the feedback effects of cloud advection and settling, gas phase element advection and depletion/replenishment and the radiative effects of cloud opacity. We model the cloud particles as a mix of mineral materials which change in size and composition as they travel through atmospheric thermo-chemical environments. All local cloud properties such as number density, grain size and material composition are time-dependently calculated. Gas phase element depletion as a result of cloud formation is included in the model. In-situ \textit{effective medium theory} and Mie theory is applied to calculate the wavelength dependent opacity of the cloud component. We present a 3D cloud structure of a chemically complex, gaseous atmosphere of the hot Jupiter HD 189733b. Mean cloud particle sizes are typically sub-micron (0.01-0.5 $\mu$m) at pressures less than 1 bar with hotter equatorial regions containing the smallest grains. Denser cloud structures occur near terminator regions and deeper ($\sim$ 1 bar) atmospheric layers. Silicate materials such as MgSiO$_{3}$[s] are found to be abundant at mid-high latitudes, while TiO$_{2}$[s] and SiO$_{2}$[s] dominate the equatorial regions.

Discs of dusty debris around main-sequence star indicate fragmentation of orbiting planetesimals, and for a few A-type stars, a gas component is also seen that may come from collisionally-released volatiles. Here we find the sixth example of a CO-hosting disc, around the 30Myr old A0-star HD 32297. Two more of these CO-hosting stars, HD 21997 and 49 Cet, have also been imaged in dust with SCUBA-2 within the SONS project. A census of 27 A-type debris hosts within 125 pc now shows 7/16 detections of carbon-bearing gas within the 5-50 Myr epoch, with no detections in 11 older systems. Such a prolonged period of high fragmentation rates corresponds quite well to the epoch when most of the Earth was assembled from planetesimal collisions. Recent models propose that collisional products can be spatially asymmetric if they originate at one location in the disc, with CO particularly exhibiting this behaviour as it can photodissociate in less than an orbital period. Of the six CO-hosting systems, only beta Pic is in clear support of this hypothesis. However, radiative transfer modelling with the ProDiMo code shows that the CO is also hard to explain in a proto-planetary disc context.

The presence or absence of ice in protoplanetary disks is of great importance for the formation of planets. By enhancing the solid surface density and increasing the sticking efficiency, ice catalyzes the rapid formation of planetesimals and decreases the time scale for giant planet core accretion. Aims: In this paper we analyse the composition of the outer disk around the Herbig star HD~142527. We focus on the composition of the water ice, but also analyse the abundances of previously proposed minerals. Methods: We present new Herschel far infrared spectra and a re-reduction of archival data from the Infrared Space Observatory (ISO). We model the disk using full 3D radiative transfer to obtain the disk structure. Also, we use an optically thin analysis of the outer disk spectrum to obtain firm constraints on the composition of the dust component. Results: The water ice in the disk around HD~142527 contains a large reservoir of crystalline water ice. We determine the local abundance of water ice in the outer disk (i.e. beyond 130\,AU). The re-reduced ISO spectrum differs significantly from that previously published, but matches the new Herschel spectrum at their common wavelength range. In particular, we do not detect any significant contribution from carbonates or hydrous silicates, in contrast to earlier claims. Conclusions: The amount of water ice detected in the outer disk requires $\sim80\,$\% of the oxygen atoms. This is comparable to the water ice abundance in the outer solar system, in comets and in dense interstellar clouds. The water ice is highly crystalline while the temperatures where we detect it are too low to crystallize the water on relevant time scales. We discuss the implications of this finding.

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 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.

In protoplanetary disks micron-size dust grains coagulate to form larger structures with complex shapes and compositions. The coagulation process changes the absorption and scattering properties of particles in the disk in significant ways. To properly interpret observations of protoplanetary disks and to place these observations in the context of the first steps of planet formation, it is crucial to understand the optical properties of these complex structures. We derive the optical properties of dust aggregates using detailed computations of aggregate structures and compare these computa- tionally demanding results with approximate methods that are cheaper to compute in practice. In this way we wish to understand the merits and problems of approximate methods and define the context in which they can or cannot be used to analyze observations of objects where significant grain growth is taking place. For the detailed computations we used the discrete dipole approximation (DDA), a method able to compute the interaction of light with a complexly shaped, inhomogeneous particle. We compared the results to those obtained using spherical and irregular, homogeneous and inhomogeneous particles. While no approximate method properly reproduces all characteristics of large dust aggregates, the thermal properties of dust can be analyzed using irregularly shaped, porous, inhomogeneous grains. The asymmetry of the scattering phase function is a good indicator of aggregate size, while the degree of polarization is probably determined by the size of the constituent particles. Optical properties derived from aggregates significantly differ from the most frequently used standard ("astronomical silicate" in spherical grains). We outline a computationally fast and relatively accurate method that can be used for a multiwavelength analysis of aggregate dust in protoplanetary disks.

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.

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.

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..

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.

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]

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.

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)

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)

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.

(abridged) We calculate near-infrared thermal emission spectra using a doubling-adding radiative transfer code, which includes scattering by clouds and haze. Initial temperature profiles and cloud optical depths are taken from the drift-phoenix brown dwarf model. As is well known, cloud particles change the spectrum compared to when clouds are ignored. The clouds reduce fluxes in the near-infrared spectrum and make it redder than for the clear sky case. We also confirm that not including scattering in the spectral calculations can result in errors on the spectra of many tens of percent, both in magnitude and in variations with wavelength. This is especially apparent for particles that are larger than the wavelength and only have little iron in them. Scattering particles will show deeper absorption features than absorbing (e.g. iron) particles and particle size will also affect the calculated infrared colours. Large particles also tend to be strongly forward-scattering, and we show that assuming isotropic scattering in this case also leads to very large errors in the spectrum. Thus, care must be taken in the choice of radiative transfer method for heat balance or spectral calculations when clouds are present in the atmosphere. Besides the choice of radiative transfer method, the type of particles that are predicted by models will change conclusions about e.g. infrared colours and trace gas abundances. As a result, knowledge of the scattering properties of the clouds is essential when deriving temperature profiles or gas abundances from direct infrared observations of exoplanets or brown dwarfs and from secondary eclipse measurements of transiting exoplanets, since scattering clouds will change the depth of gas absorption features, among other things. Thus, ignoring the presence of clouds can yield retrieved properties that differ significantly from the real atmospheric properties.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Solving the continuum radiative transfer equation in high opacity media requires sophisticated numerical tools. In order to test the reliability of such tools, we present a benchmark of radiative transfer codes in a 2D disc configuration. We test the accuracy of seven independently developed radiative transfer codes by comparing the temperature structures, spectral energy distributions, scattered light images, and linear polarisation maps that each model predicts for a variety of disc opacities and viewing angles. The test cases have been chosen to be numerically challenging, with midplane optical depths up 10^6, a sharp density transition at the inner edge and complex scattering matrices. We also review recent progress in the implementation of the Monte Carlo method that allow an efficient solution to these kinds of problems and discuss the advantages and limitations of Monte Carlo codes compared to those of discrete ordinate codes. For each of the test cases, the predicted results from the radiative transfer codes are within good agreement. The results indicate that these codes can be confidently used to interpret present and future observations of protoplanetary discs.

The atmospheres of substellar objects contain clouds of oxides, iron, silicates, and other refractory condensates. Water clouds are expected in the coolest objects. The opacity of these `dust' clouds strongly affects both the atmospheric temperature-pressure profile and the emergent flux. Thus any attempt to model the spectra of these atmospheres must incorporate a cloud model. However the diversity of cloud models in atmospheric simulations is large and it is not always clear how the underlying physics of the various models compare. Likewise the observational consequences of different modeling approaches can be masked by other model differences, making objective comparisons challenging. In order to clarify the current state of the modeling approaches, this paper compares five different cloud models in two sets of tests. Test case 1 tests the dust cloud models for a prescribed L, L--T, and T-dwarf atmospheric (temperature T, pressure p, convective velocity vconv)-structures. Test case 2 compares complete model atmosphere results for given (effective temperature Teff, surface gravity log g). All models agree on the global cloud structure but differ in opacity-relevant details like grain size, amount of dust, dust and gas-phase composition. Comparisons of synthetic photometric fluxes translate into an modelling uncertainty in apparent magnitudes for our L-dwarf (T-dwarf) test case of 0.25 < \Delta m < 0.875 (0.1 < \Delta m M 1.375) taking into account the 2MASS, the UKIRT WFCAM, the Spitzer IRAC, and VLT VISIR filters with UKIRT WFCAM being the most challenging for the models. (abr.)

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]

We aim to understand cloud formation in substellar objects. We combined the non-equilibrium, stationary cloud model of Helling, Woitke & Thi (2008; seed formation, growth, evaporation, gravitational settling, element conservation) with the general-purpose model atmosphere code PHOENIX (radiative transfer, hydrostatic equilibrium, mixing length theory, chemical equilibrium) in order to consistently calculate cloud formation and radiative transfer with their feedback on convection and gas phase depletion. We calculate the complete 1D model atmosphere structure and the chemical details of the cloud layers. The DRIFT-PHOENIX models enable the first stellar atmosphere simulation that is based on the actual cloud formation process. The resulting (T,p) profiles differ considerably from the previous limiting PHOENIX cases DUSTY and COND. A tentative comparison with observations demonstrates that the determination of effective temperatures based on simple cloud models has to be applied with care. Based on our new models, we suggest a mean Teff=1800K for the L-dwarf twin-binary system DENIS J0205-1159 which is up to 500K hotter than suggested in the literature. We show transition spectra for gas-giant planets which form dust clouds in their atmospheres and evaluate photometric fluxes for a WASP-1 type system.

We aim to understand the correlation between cloud formation and alkali line formation in substellar atmospheres.We perform line profile calculations for Na I and K I based on the coupling of our kinetic model for the formation and composition of dust grains with 1D radiative transfer calculations in atmosphere models for brown dwarfs and giant gas planets. The Na I and K I line profiles sensibly depend on the way clouds are treated in substellar atmosphere simulations. The kinetic dust formation model results in the highest pseudo-continuum compared to the limiting cases.

A test case comparison is presented for different dust cloud model approaches applied in brown dwarfs and giant gas planets. We aim to achieve more transparency in evaluating the uncertainty inherent to theoretical modelling. We show in how far model results for characteristic dust quantities vary due to different assumptions. We also demonstrate differences in the spectral energy distributions resulting from our individual cloud modelling in 1D substellar atmosphere simulations

New dynamical models for dust-driven winds of oxygen-rich AGB stars are presented which include frequency-dependent Monte Carlo radiative transfer by means of a sparse opacity distribution technique and a time-dependent treatment of the nucleation, growth and evaporation of inhomogeneous dust grains composed of a mixture of Mg2SiO4, SiO2, Al2O3, TiO2, and solid Fe. The frequency-dependent treatment of radiative transfer reveals that the gas is cold close to the star (700-900 K at 1.5-2 R*) which facilitates the nucleation process. The dust temperatures are strongly material-dependent, with differences as large as 1000 K for different pure materials, which has an important influence on the dust formation sequence. Two dust layers are formed in the dynamical models: almost pure glassy Al2O3 close to the star (r > 1.5 R*) and the more opaque Fe-poor Mg-Fe-silicates further out. Solid Fe or Fe-rich silicates are found to be the only condensates that can efficiently absorb the stellar light in the near IR. Consequently, they play a crucial role for the wind driving mechanism and act as thermostat. Only small amounts of Fe can be incorporated into the grains, because otherwise the grains get too hot. Thus, the models reveal almost no mass loss, and no dust shells.

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.

New axisymmetric (2D) models for dust-driven winds of C-stars are presented which include hydrodynamics with radiation pressure on dust, equilibrium chemistry and time-dependent dust formation with coupled grey Monte Carlo radiative transfer. Considering the most simple case without stellar pulsation (hydrostatic inner boundary condition) these models reveal a more complex picture of the dust formation and wind acceleration as compared to earlier published spherically symmetric (1D) models. The so-called exterior $\kappa$-mechanism causes radial oscillations with short phases of active dust formation between longer phases without appreciable dust formation, just like in the 1D models. However, in 2D geometry, the oscillations can be out-of-phase at different places above the stellar atmosphere which result in the formation of dust arcs or smaller caps that only occupy a certain fraction of the total solid angle. These dust structures are accelerated outward by radiation pressure, expanding radially and tangentially like mushroom clouds, while dust-poor matter is falling back towards the star at other places. A highly dynamical and turbulent dust formation zone is created in this way, which again leads to inhomogeneous dust production. Further away from the star, flow instabilities (e.g. Rayleigh-Taylor) have time to fragment the outward moving arcs and shells to produce numerous small-scale cloud-like sub-structures.

The large scale structure of a brown dwarf atmosphere is determined by an interplay of convection, radiation, dust formation, and gravitational settling, which possibly provides an explanation for the observed variability. The result is an element depletion of the dust forming regions and an element enrichment of the dust evaporating sites. The formation of dust cloud structures in substellar atmospheres is demonstrated based on a consistent theoretical description of dust formation and destruction, gravitational settling, and element depletion including the effect of convective overshoot. Since the viscosity is small in brown dwarf atmospheres, the convection creates a turbulent environment with fluctuations of all thermodynamic state variables on a wide range of spatial scales. Hence, the classical turbulent closure problem needs to be tackled in connection with dust formation in substellar objects, because a complete three-dimensional and time-dependent solution of the model equations is simply not possible. Structure formation may be seeded on the smallest scales, when chemical processes are involved. In order to understand the interaction of turbulence and dust formation, we have performed investigations of the smallest scale regimes in 1D and in 2D in order to identify the governing processes of the unresolved scale regime.

This paper reports on computational evidence for the formation of cloud-like dust structures around C-rich AGB stars. This spatio-temporal structure formation process is caused by a radiative/thermal instability of dust forming gases as identified by Woitke et al.(2000). Our 2D (axisymmetric) models combine a time-dependent description of the dust formation process according to Gail & Sedlmayr (1988) with detailed, frequency-dependent continuum radiative transfer by means of a Monte Carlo method (Niccolini et al.2003) in an otherwise static medium (v=0). These models show that the formation of dust behind already condensed regions, which shield the stellar radiation field, is strongly favoured. In the shadow of these clouds, the temperature decreases by several hundred Kelvin which triggers the subsequent formation of dust and ensures its thermal stability. Considering an initially dust-free gas with small density inhomogeneities, we find that finger-like dust structures develop which are cooler than the surroundings and point towards the centre of the radiant emission, similar to the cometary knots observed in planetary nebulae and star formation regions. Compared to a spherical symmetric reference model, the clumpy dust distribution has little effect on the spectral energy distribution, but dominates the optical appearance in near IR monochromatic images.

Dust formation in brown dwarf atmospheres is studied by utilising a model for driven turbulence in the mesoscopic scale regime. We apply a pseudo-spectral method where waves are created and superimposed within a limited wavenumber interval. The turbulent kinetic energy distribution follows the Kolmogoroff spectrum which is assumed to be the most likely value. Such superimposed, stochastic waves may occur in a convectively active environment. They cause nucleation fronts and nucleation events and thereby initiate the dust formation process which continues until all condensible material is consumed. Small disturbances are found to have a large impact on the dust forming system. An initially dust-hostile region, which may originally be optically thin, becomes optically thick in a patchy way showing considerable variations in the dust properties during the formation process. The dust appears in lanes and curls as a result of the interaction with waves, i.e. turbulence, which form larger and larger structures with time. Aiming on a physical understanding of the variability of brown dwarfs, related to structure formation in substellar atmospheres, we work out first necessary criteria for small-scale closure models to be applied in macroscopic simulations of dust forming astrophysical systems.

We present benchmark problems and solutions for the continuum radiative transfer (RT) in a 2D disk configuration. The reliability of three Monte-Carlo and two grid-based codes is tested by comparing their results for a set of well-defined cases which differ for optical depth and viewing angle. For all the configurations, the overall shape of the resulting temperature and spectral energy distribution is well reproduced. The solutions we provide can be used for the verification of other RT codes.We also point out the advantages and disadvantages of the various numerical techniques applied to solve the RT problem.