Context: The Galilean satellites are thought to have formed from a circumplanetary disk (CPD) surrounding Jupiter. When it reached a critical mass, Jupiter opened an annular gap in the solar protoplanetary disk (PPD) that might have exposed the CPD to radiation from the young Sun or from the stellar cluster in which the Solar System formed. Aims: We investigate the radiation field to which the Jovian CPD was exposed during the process of satellite formation. The resulting photoevaporation of the CPD is studied in this context to constrain possible formation scenarios for the Galilean satellites and explain architectural features of the Galilean system. Methods: We constructed a model for the stellar birth cluster to determine the intracluster far-ultraviolet (FUV) radiation field. We employed analytical annular gap profiles informed by hydrodynamical simulations to investigate a range of plausible geometries for the Jovian gap. We used the radiation thermochemical code ProDiMo to evaluate the incident radiation field in the Jovian gap and the photoevaporation of an embedded 2D axisymmetric CPD. Results: We derive the time-dependent intracluster FUV radiation field for the solar birth cluster over 10 Myr. We find that intracluster photoevaporation can cause significant truncation of the Jovian CPD. We determine steady-state truncation radii for possible CPDs, finding that the outer radius is proportional to the accretion rate $\dot{M}^{0.4}$. For CPD accretion rates $\dot M < 10^{-12} M_{\odot}$ yr$^{-1}$, photoevaporative truncation explains the lack of additional satellites outside the orbit of Callisto. For CPDs of mass $M_{\rm CPD} < 10^{-6.2 M_{\odot}}$ , photoevaporation can disperse the disk before Callisto is able to migrate into the Laplace resonance. This explains why Callisto is the only massive satellite that is excluded from the resonance.

Gas has been detected in a number of debris disks. It is likely secondary, i.e. produced by colliding solids. Here, we report ALMA Band 8 observations of neutral carbon in the CO-rich debris disk around the 15--30 Myr old A-type star HD 32297. We find that C$^0$ is located in a ring at $\sim$110 au with a FWHM of $\sim$80 au, and has a mass of $(3.5\pm0.2)\times10^{-3}$ M$_\oplus$. Naively, such a surprisingly small mass can be accumulated from CO photo-dissociation in a time as short as $\sim$10$^4$ yr. We develop a simple model for gas production and destruction in this system, properly accounting for CO self-shielding and shielding by neutral carbon, and introducing a removal mechanism for carbon gas. We find that the most likely scenario to explain both C$^0$ and CO observations, is one where the carbon gas is rapidly removed on a timescale of order a thousand years and the system maintains a very high CO production rate of $\sim$15 M$_\oplus$ Myr$^{-1}$, much higher than the rate of dust grind-down. We propose a possible scenario to meet these peculiar conditions: the capture of carbon onto dust grains, followed by rapid CO re-formation and re-release. In steady state, CO would continuously be recycled, producing a CO-rich gas ring that shows no appreciable spreading over time. This picture might be extended to explain other gas-rich debris disks.

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.

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

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

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

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.

Planets form in protoplanetary disks and inherit their chemical compositions. It is thus crucial to map the distribution and investigate the formation of simple organics, such as formaldehyde and methanol, in protoplanetary disks. We analyze ALMA observations of the nearby disk-jet system around the T Tauri star DG Tau in the o-H$_2$CO $3_{1,2}-2_{1,1}$ and CH$_3$OH $3_{-2,2}-4_{-1,4}$ E, $5_{0,5}-4_{0,4}$ A transitions at an unprecedented resolution of $\sim0.15"$, i.e., $\sim18$ au at a distance of 121 pc. The H$_2$CO emission originates from a rotating ring extending from $\sim40$ au with a peak at $\sim62$ au, i.e., at the edge of the 1.3mm dust continuum. CH$_3$OH emission is not detected down to an r.m.s. of 3 mJy/beam in the 0.162 km/s channel. Assuming an ortho-to-para ratio of 1.8-2.8 the ring- and disk-height-averaged H$_2$CO column density is $\sim0.3-4\times10^{14}$ cm$^{-2}$, while that of CH$_3$OH is $<0.04-0.7\times10^{14}$ cm$^{-2}$. In the inner $40$ au no o-H$_2$CO emission is detected with an upper limit on its beam-averaged column density of $\sim0.5-6\times10^{13}$ cm$^{-2}$. The H$_2$CO ring in the disk of DG Tau is located beyond the CO iceline (R$_{\rm CO}\sim30$ au). This suggests that the H$_2$CO abundance is enhanced in the outer disk due to formation on grain surfaces by the hydrogenation of CO ice. The emission peak at the edge of the mm dust continuum may be due to enhanced desorption of H$_2$CO in the gas phase caused by increased UV penetration and/or temperature inversion. The CH$_3$OH/H$_2$CO abundance ratio is $<1$, in agreement with disk chemistry models. The inner edge of the H$_2$CO ring coincides with the radius where the polarization of the dust continuum changes orientation, hinting at a tight link between the H$_2$CO chemistry and the dust properties in the outer disk and at the possible presence of substructures in the dust distribution.

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.

VLT instruments and ALMA have revolutionized in the past five years our view and understanding of how disks turn into planetary systems. They provide exquisite insights into non-axisymmetric structures likely closely related to ongoing planet formation processes. The following cannot be a complete review of the physical and chemical properties of disks; instead I focus on a few selected aspects. I will review our current understanding of the physical properties (e.g. solid and gas mass content, snow and ice lines) and chemical composition of planet forming disks at ages of 1-few Myr, especially in the context of the planetary systems that are forming inside them. I will highlight recent advances achieved by means of consistent multi-wavelength studies of gas AND dust in protoplanetary disks.

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.

The debris disk around $\beta$~Pictoris is known to contain gas. Previous ALMA observations revealed a CO belt at $\sim$85 au with a distinct clump, interpreted as a location of enhanced gas production. Photodissociation converts CO into C and O within $\sim$50 years. We resolve CI emission at 492 GHz using ALMA and study its spatial distribution. CI shows the same clump as seen for CO. This is surprising, as C is expected to quickly spread in azimuth. We derive a low C mass (between $5\times10^{-4}$ and $3.1\times10^{-3}$ M$_\oplus$), indicating that gas production started only recently (within $\sim$5000 years). No evidence is seen for an atomic accretion disk inwards of the CO belt, perhaps because the gas did not yet have time to spread radially. The fact that C and CO share the same asymmetry argues against a previously proposed scenario where the clump is due to an outward migrating planet trapping planetesimals in an resonance; nor can the observations be explained by an eccentric planetesimal belt secularly forced by a planet. Instead, we suggest that the dust and gas disks should be eccentric. Such a configuration, we further speculate, might be produced by a recent tidal disruption event. Assuming that the disrupted body has had a CO mass fraction of 10%, its total mass would be $\gtrsim$3 $M_\mathrm{Moon}$.

This paper investigates how the far-IR water ice features can be used to infer properties of disks around T Tauri stars and the water ice thermal history. We explore the power of future observations with SOFIA/HIRMES and SPICA's proposed far-IR instrument SAFARI. A series of detailed radiative transfer disk models around a representative T Tauri star are used to investigate how the far-IR water ice features at 45 and 63 micron change with key disk properties: disk size, grain sizes, disk dust mass, dust settling, and ice thickness. In addition, a series of models is devised to calculate the water ice emission features from warmup, direct deposit and cooldown scenarios of the water ice in disks. Photodesorption from icy grains in disk surfaces weakens the mid-IR water ice features by factors 4-5. The far-IR water ice emission features originate from small grains at the surface snow line in disks at distance of 10-100 au. Unless this reservoir is missing in disks (e.g. transitional disks with large cavities), the feature strength is not changing. Grains larger than 10 micron do not contribute to the features. Grain settling (using turbulent description) is affecting the strength of the ice features by at most 15%. The strength of the ice feature scales with the disk dust mass and water ice fraction on the grains, but saturates for dust masses larger than 1.e-4 Msun and for ice mantles that increase the dust mass by more than 50%. The various thermal histories of water ice leave an imprint on the shape of the features (crystalline/amorphous) as well as on the peak strength and position of the 45 micron feature. SOFIA/HIRMES can only detect crystalline ice features much stronger than simulated in our standard T Tauri disk model in deep exposures (1 hr). SPICA/SAFARI can detect the typical ice features in our standard T Tauri disk model in short exposures (10 min). (abbreviated)

Measurements in the infrared wavelength domain allow us to assess directly the physical state and energy balance of cool matter in space, thus enabling the detailed study of the various processes that govern the formation and early evolution of stars and planetary systems in galaxies over cosmic time. Previous infrared missions, from IRAS to Herschel, have revealed a great deal about the obscured Universe, but sensitivity has been limited because up to now it has not been possible to fly a telescope that is both large and cold. SPICA is a mission concept aimed at taking the next step in mid- and far-infrared observational capability by combining a large and cold telescope with instruments employing state-of-the-art ultra-sensitive detectors. The mission concept foresees a 2.5-meter diameter telescope cooled to below 8 K. With cooling provided by mechanical coolers instead of depending on a limited cryogen supply, the mission lifetime can extend significantly beyond the required three years. SPICA offers instrumentation with spectral resolving powers ranging from R ~50 through 11000 in the 17-230 $\mu$m domain as well as R~28.000 spectroscopy between 12 and 18 $\mu$m. Additionally SPICA will provide efficient 30-37 $\mu$m broad band mapping, and polarimetric imaging in the 100-350 $\mu$m range. SPICA will provide unprecedented spectroscopic sensitivity of ~5 x $10^{-20}$ W/m$^2$ (5$\sigma$/1hr) - at least two orders of magnitude improvement over what has been attained to date. With this exceptional leap in performance, new domains in infrared astronomy will become accessible, allowing us, for example, to unravel definitively galaxy evolution and metal production over cosmic time, to study dust formation and evolution from very early epochs onwards, and to trace the formation history of planetary systems.

Our understanding of protoplanetary disks has greatly improved over the last decade due to a wealth of data from new facilities. Unbiased dust surveys with Spitzer leave us with good constraints on the dust dispersal timescale of small grains in the terrestrial planet forming region. In the ALMA era, this can be confronted for the first time also with evolutionary timescales of mm grains in the outer disk. Gas surveys in the context of the existing multi-wavelength dust surveys will be a key in large statistical studies of disk gas evolution. Unbiased gas surveys are limited to ALMA CO submm surveys, where the quantitative interpretation is still debated. Herschel gas surveys have been largely biased, but [OI] 63 mic surveys and also accretion tracers agree qualitatively with the evolutionary timescale of small grains in the inner disk. Recent advances achieved by means of consistent multi-wavelength studies of gas AND dust in planet forming disks reveal the subtleties of the quantitative interpretation of gas surveys. Observational methods to determine disk masses e.g. from CO submm lines require the knowledge of the dust properties in the disk. Understanding not only the gas evolution, but also its chemical composition will provide crucial input for planet formation models. Kinetic chemical results give profoundly different answers than thermodynamic equilibrium in terms of the C/O ratios as well as the water ice/rock ratios. Again, dust has a key impact on the chemical evolution and composition of the gas. Grain growth for example affects freeze-out processes and strongly increases the cosmic ray induced UV field.

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.

The properties of disks around brown dwarfs and very-low mass stars (hereafter VLMOs) provide important boundary conditions on the process of planet formation and inform us about the numbers and masses of planets than can form in this regime. We use the Herschel Space Observatory PACS spectrometer to measure the continuum and [OI] 63um line emission towards 11 VLMOs with known disks in the Taurus and Chamaeleon I star-forming regions. We fit radiative transfer models to the spectral energy distributions of these sources. Additionally, we carry out a grid of radiative transfer models run in a regime that connects the luminosity of our sources with brighter T~Tauri stars. We find VLMO disks with sizes [1.3--78] au, smaller than typical T~Tauri disks, fit well the spectral energy distributions assuming disk geometry and dust properties are stellar-mass independent. Reducing the disk size increases the disk temperature and we show that VLMOs do not follow previously derived disk temperature-stellar luminosity relationships if the disk outer radius scales with stellar mass. Only 2 out of 11 sources are detected in [OI] despite a better sensitivity than was achieved for T Tauri stars, suggesting that VLMO disks are underluminous. Using thermochemical models we show that smaller disks can lead to the unexpected [OI] non-detections in our sample. The disk outer radius is an important factor in determining the gas and dust observables. Hence, spatially resolved observations with ALMA -- to establish if and how disk radii scale with stellar mass -- should be pursued further.

At early stages of stellar evolution young stars show powerful jets and/or outflows that interact with protoplanetary discs and their surroundings. Despite the scarce knowledge about the interaction of jets and/or outflows with discs, spectroscopic studies based on Herschel and ISO data suggests that gas shocked by jets and/or outflows can be traced by far-IR (FIR) emission in certain sources. We want to provide a consistent catalogue of selected atomic ([OI] and [CII]) and molecular (CO, OH, and H$_{2}$O) line fluxes observed in the FIR, separate and characterize the contribution from the jet and the disc to the observed line emission, and place the observations in an evolutionary picture. The atomic and molecular FIR (60-190 $\rm \mu m$) line emission of protoplanetary discs around 76 T Tauri stars located in Taurus are analysed. The observations were carried out within the Herschel key programme Gas in Protoplanetary Systems (GASPS). The spectra were obtained with the Photodetector Array Camera and Spectrometer (PACS). The sample is first divided in outflow and non-outflow sources according to literature tabulations. With the aid of archival stellar/disc and jet/outflow tracers and model predictions (PDRs and shocks), correlations are explored to constrain the physical mechanisms behind the observed line emission. The much higher detection rate of emission lines in outflow sources and the compatibility of line ratios with shock model predictions supports the idea of a dominant contribution from the jet/outflow to the line emission, in particular at earlier stages of the stellar evolution as the brightness of FIR lines depends in large part on the specific evolutionary stage. [Abridged Abstract]

We present ~0.4 resolution images of CO(3-2) and associated continuum emission from the gas-bearing debris disk around the nearby A star 49 Ceti, observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We analyze the ALMA visibilities in tandem with the broad-band spectral energy distribution to measure the radial surface density profiles of dust and gas emission from the system. The dust surface density decreases with radius between ~100 and 310 au, with a marginally significant enhancement of surface density at a radius of ~110 au. The SED requires an inner disk of small grains in addition to the outer disk of larger grains resolved by ALMA. The gas disk exhibits a surface density profile that increases with radius, contrary to most previous spatially resolved observations of circumstellar gas disks. While ~80% of the CO flux is well described by an axisymmetric power-law disk in Keplerian rotation about the central star, residuals at ~20% of the peak flux exhibit a departure from axisymmetry suggestive of spiral arms or a warp in the gas disk. The radial extent of the gas disk (~220 au) is smaller than that of the dust disk (~300 au), consistent with recent observations of other gas-bearing debris disks. While there are so far only three broad debris disks with well characterized radial dust profiles at millimeter wavelengths, 49 Ceti's disk shows a markedly different structure from two radially resolved gas-poor debris disks, implying that the physical processes generating and sculpting the gas and dust are fundamentally different.

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.

Protoplanetary disks undergo substantial mass-loss by photoevaporation, a mechanism which is crucial to their dynamical evolution. However, the processes regulating the gas energetics have not been well constrained by observations so far. We aim at studying the processes involved in disk photoevaporation when it is driven by far-UV photons. We present a unique Herschel survey and new ALMA observations of four externally-illuminated photoevaporating disks (a.k.a. proplyds). For the analysis of these data, we developed a 1D model of the photodissociation region (PDR) of a proplyd, based on the Meudon PDR code and computed the far infrared line emission. We successfully reproduce most of the observations and derive key physical parameters, i.e. densities at the disk surface of about $10^{6}$ cm$^{-3}$ and local gas temperatures of about 1000 K. Our modelling suggests that all studied disks are found in a transitional regime resulting from the interplay between several heating and cooling processes that we identify. These differ from those dominating in classical PDRs, i.e. grain photo-electric effect and cooling by [OI] and [CII] FIR lines. This energetic regime is associated to an equilibrium dynamical point of the photoevaporation flow: the mass-loss rate is self-regulated to set the envelope column density at a value that maintains the temperature at the disk surface around 1000 K. From our best-fit models, we estimate mass-loss rates - of the order of $10^{-7}$ $\mathrm{M}_\odot$/yr - that are in agreement with earlier spectroscopic observation of ionised gas tracers. This holds only if we assume an evaporation flow launched from the disk surface at sound speed (supercritical regime). We have identified the energetic regime regulating FUV-photoevaporation in proplyds. This regime could be implemented into models of the dynamical evolution of protoplanetary disks.

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

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.

We present H-band (1.6 {\mu}m) scattered light observations of the transitional disk RX J1615.3-3255, located in the ~1 Myr old Lupus association. From a polarized intensity image, taken with the HiCIAO instrument of the Subaru Telescope, we deduce the position angle and the inclination angle of the disk. The disk is found to extend out to 68 $\pm$ 12 AU in scattered light and no clear structure is observed. Our inner working angle of 24 AU does not allow us to detect a central decrease in intensity similar to that seen at 30 AU in the 880 {\mu}m continuum observations. We compare the observations with multiple disk models based on the Spectral Energy Distribution (SED) and submm interferometry and find that an inner rim of the outer disk at 30 AU containing small silicate grains produces a polarized intensity signal which is an order of magnitude larger than observed. We show that a model in which the small dust grains extend smoothly into the cavity found for large grains is closer to the actual H-band observations. A comparison of models with different dust size distributions suggests that the dust in the disk might have undergone significant processing compared to the interstellar medium.

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.

Gas plays a major role in the dynamical evolution of young stellar objects. Its interaction with the dust is the key to our understanding planet formation later on in the protoplanetary disc stage. Studying the gas content is a crucial step towards understanding YSO and planet formation. Such a study can be made through spectroscopic observations of emission lines in the far-infrared, where some of the most important gas coolants emit. We provide a compilation of observations of far-IR lines in 362 young stellar objects covering all evolutionary stages, from Class 0 to Class III with debris discs. In the present paper we focus on [OI] and o-H2O emission at 63 microns. We have retrieved all the available Herschel -PACS spectroscopic observations at 63 microns that used the dominant observing mode, the chop-nod technique. We provide measurements of line fluxes for the [OI] 3P1-3P2 and o-H2O 808-717 transitions at 63 microns computed using different methods. We check for spatially extended emission and further study the presence of multiple dynamical components in line emission. The final compilation consists of line and continuum fluxes at 63 microns for a total of 362 young stellar objects (YSOs). We detected [OI] line emission at 63 microns in 194 sources out of 362, and line absorption in another five sources. o-H2O was detected in 42 sources. We find evidence of extended [OI] emission in 77 sources, and detect 3sigma residual emission in 71 of them. We also looked for different components contributing to the line emission, and found evidence for multiple components in 30 sources. We explored correlations between line emission and continuum emission and find a clear correlation between WISE fluxes from 4.6 to 22 microns and [OI] line emission. We conclude that the observed emission is typically a combination of disc/envelope and jet emission.

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

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

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

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.

Protoplanetary discs are the birthplace for planets. Studying protoplanetary discs is the key to constraining theories of planet formation. By observing dust and gas in associations at different ages we can study the evolution of these discs, their clearing timescales, and their physical and geometrical properties. The stellar association Eta Cha is peculiar; some members still retain detectable amounts of gas in their discs at the late age of 7 Myr, making it one of the most interesting young stellar associations in the solar neighbourhood. We characterise the properties of dust and gas in protoplanetary and transitional discs in the Eta Cha young cluster, with special emphasis on explaining the peculiarities that lead to the observed high disc detection fraction and prominent IR excesses at an age of 7 Myr. We observed 17 members of the Eta Cha association with Herschel-PACS in photometric mode and line spectroscopic mode. A subset of members were also observed in range spectroscopic mode. The observations trace [OI] and $H2O emissions at 63.18 and 63.32 microns, respectively, as well as CO, OH, CH+ and [CII] at different wavelengths for those systems observed in range mode. The photometric observations were used to build complete spectral energy distributions (SEDs) from the optical to the far-IR. High-resolution multi-epoch optical spectra with high signal-to-noise ratios were also analysed to study the multiplicity of the sources and look for further gas (accreting) and outflow indicators.

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.

Context: The 440 Myr old main-sequence A-star Fomalhaut is surrounded by an eccentric debris belt with sharp edges. This sort of a morphology is usually attributed to planetary perturbations, but the orbit of the only planetary candidate detected so far, Fomalhaut b, is too eccentric to efficiently shape the belt. Alternative models that could account for the morphology without invoking a planet are stellar encounters and gas-dust interactions. Aims: We aim to test the possibility of gas-dust interactions as the origin of the observed morphology by putting upper limits on the total gas content of the Fomalhaut belt. Methods: We derive upper limits on the CII 158 $\mu$m and OI 63 $\mu$m emission by using non-detections from the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel Space Observatory. Line fluxes are converted into total gas mass using the non-local thermodynamic equilibrium (non-LTE) code RADEX. We consider two different cases for the elemental abundances of the gas: solar abundances and abundances similar to those observed for the gas in the $\beta$ Pictoris debris disc. Results: The gas mass is shown to be below the millimetre dust mass by a factor of at least $\sim$3 (for solar abundances) respectively $\sim$300 (for $\beta$ Pic-like abundances). Conclusions: The lack of gas co-spatial with the dust implies that gas-dust interactions cannot efficiently shape the Fomalhaut debris belt. The morphology is therefore more likely due to a yet unseen planet (Fomalhaut c) or stellar encounters.

Gas plays a major role in the dynamical evolution of protoplanetary discs. Its coupling with the dust is the key to our understanding planetary formation. Studying the gas content is therefore a crucial step towards understanding protoplanetary discs evolution. Such a study can be made through spectroscopic observations of emission lines in the far-infrared, where some of the most important gas coolants emit, such as the [OI] 3P1-3 P2 transition at 63.18 microns. We aim at characterising the gas content of protoplanetary discs in the intermediate-aged Chamaeleon II (Cha II) star forming region. We also aim at characterising the gaseous detection fractions within this age range, which is an essential step tracing gas evolution with age in different star forming regions. We obtained Herschel-PACS line scan spectroscopic observations at 63 microns of 19 Cha II Class I and II stars. The observations were used to trace [OI] and o-H2O at 63 microns. The analysis of the spatial distribution of [OI], when extended, can be used to understand the origin of the emission. We have detected [OI] emission toward seven out of the nineteen systems observed, and o-H2O emission at 63.32 microns in just one of them, Sz 61. Cha II members show a correlation between [OI] line fluxes and the continuum at 70 microns, similar to what is observed in Taurus. We analyse the extended [OI] emission towards the star DK Cha and study its dynamical footprints in the PACS Integral Field Unit (IFU). We conclude that there is a high velocity component from a jet combined with a low velocity component with an origin that may be a combination of disc, envelope and wind emission. The stacking of spectra of objects not detected individually in [OI] leads to a marginal 2.6sigma detection that may indicate the presence of gas just below our detection limits for some, if not all, of them.

We present Hubble Space Telescope STIS far-UV spectra of the edge-on disk around 49 Ceti, one of the very few debris disks showing sub-mm CO emission. Many atomic absorption lines are present in the spectra, most of which arise from circumstellar gas lying along the line-of-sight to the central star. We determined the line-of-sight CI column density, estimated the total carbon column density, and set limits on the OI column density. Surprisingly, no line-of-sight CO absorption was seen. We discuss possible explanations for this non-detection, and present preliminary estimates of the carbon abundances in the line-of-sight gas. The C/Fe ratio is much greater than the solar value, suggesting that 49 Cet harbors a volatile-rich gas disk similar to that of Beta Pictoris.

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.

Context. Debris discs are thought to be formed through the collisional grinding of planetesimals, and can be considered as the outcome of planet formation. Understanding the properties of gas and dust in debris discs can help us to comprehend the architecture of extrasolar planetary systems. Herschel Space Observatory far-infrared (IR) photometry and spectroscopy have provided a valuable dataset for the study of debris discs gas and dust composition. This paper is part of a series of papers devoted to the study of Herschel PACS observations of young stellar associations. Aims. This work aims at studying the properties of discs in the Beta Pictoris Moving Group (BPMG) through far-IR PACS observations of dust and gas. Methods. We obtained Herschel-PACS far-IR photometric observations at 70, 100 and 160 microns of 19 BPMG members, together with spectroscopic observations of four of them. Spectroscopic observations were centred at 63.18 microns and 157 microns, aiming to detect [OI] and [CII] emission. We incorporated the new far-IR observations in the SED of BPMG members and fitted modified blackbody models to better characterise the dust content. Results. We have detected far-IR excess emission toward nine BPMG members, including the first detection of an IR excess toward HD 29391.The star HD 172555, shows [OI] emission, while HD 181296, shows [CII] emission, expanding the short list of debris discs with a gas detection. No debris disc in BPMG is detected in both [OI] and [CII]. The discs show dust temperatures in the range 55 to 264 K, with low dust masses (6.6*10^{-5} MEarth to 0.2 MEarth) and radii from blackbody models in the range 3 to 82 AU. All the objects with a gas detection are early spectral type stars with a hot dust component.

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.

The measurement of isotopic ratios provides a privileged insight both into nucleosynthesis and into the mechanisms operating in stellar envelopes, such as gravitational settling. In this article, we give a few examples of how isotopic ratios can be determined from high-resolution, high-quality stellar spectra. We consider examples of the lightest elements, H and He, for which the isotopic shifts are very large and easily measurable, and examples of heavier elements for which the determination of isotopic ratios is more difficult. The presence of 6Li in the stellar atmospheres causes a subtle extra depression in the red wing of the 7Li 670.7 nm doublet which can only be detected in spectra of the highest quality. But even with the best spectra, the derived $^6$Li abundance can only be as good as the synthetic spectra used for their interpretation. It is now known that 3D non-LTE modelling of the lithium spectral line profiles is necessary to account properly for the intrinsic line asymmetry, which is produced by convective flows in the atmospheres of cool stars, and can mimic the presence of 6Li. We also discuss briefly the case of the carbon isotopic ratio in metal-poor stars, and provide a new determination of the nickel isotopic ratios in the solar atmosphere.

The far-IR range is a critical wavelength range to characterize the physical and chemical processes that transform the interstellar material into stars and planets. Objects in the earliest phases of stellar and planet evolution release most of their energy at these long wavelengths. In this contribution we briefly summarise some of the most relevant scientific advances achieved by the Herschel Space Observatory in the field. We also anticipate those that will be made possible by the large increase in sensitivity of SPICA cooled telescope. It is concluded that only through sensitive far-IR observations much beyond Herschel capabilities we will be able to constrain the mass, the energy budget and the water content of hundreds of protostars and planet-forming disks.

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

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

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

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

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

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

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

Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system's counterparts are the asteroid and Edgeworth-Kuiper belts. The DUNES survey aims at detecting extra-solar analogues to the Edgeworth-Kuiper belt around solar-type stars, putting in this way the solar system into context. The survey allows us to address some questions related to the prevalence and properties of planetesimal systems. We used {\it Herschel}/PACS to observe a sample of nearby FGK stars. Data at 100 and 160 $\mu$m were obtained, complemented in some cases with observations at 70 $\mu$m, and at 250, 350 and 500 $\mu$m using SPIRE. The observing strategy was to integrate as deep as possible at 100 $\mu$m to detect the stellar photosphere. Debris discs have been detected at a fractional luminosity level down to several times that of the Edgeworth-Kuiper belt. The incidence rate of discs around the DUNES stars is increased from a rate of $\sim$ 12.1% $\pm$ 5% before \emph{Herschel} to $\sim$ 20.2% $\pm$ 2%. A significant fraction ($\sim$ 52%) of the discs are resolved, which represents an enormous step ahead from the previously known resolved discs. Some stars are associated with faint far-IR excesses attributed to a new class of cold discs. Although it cannot be excluded that these excesses are produced by coincidental alignment of background galaxies, statistical arguments suggest that at least some of them are true debris discs. Some discs display peculiar SEDs with spectral indexes in the 70-160$\mu$m range steeper than the Rayleigh-Jeans one. An analysis of the debris disc parameters suggests that a decrease might exist of the mean black body radius from the F-type to the K-type stars. In addition, a weak trend is suggested for a correlation of disc sizes and an anticorrelation of disc temperatures with the stellar age.

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.

Proplyds are photodissociation region (PDR)-like cometary cocoons around young stars which are thought to originate through photo-evaporation of the central protoplanetary disk by external UV radiation from the nearby OB stars. This letter presents spatially resolved mid-infrared imaging and spectroscopy of the proplyd HST10 obtained with the VLT/VISIR instrument. These observations allow us to detect Polycyclic Aromatic Hydrocarbons (PAH) emission in the proplyd photodissociation region and to study the general properties of PAHs in proplyds for the first time. We find that PAHs in HST10 are mostly neutral and at least 50 times less abundant than typical values found for the diffuse ISM or the nearby Orion Bar. With such a low PAH abundance, photoelectric heating is significantly reduced. If this low abundance pertains also to the original disk material, gas heating rates could be too low to efficiently drive photoevaporation unless other processes can be identified. Alternatively, the model behind the formation of proplyds as evaporating disks may have to be revised.

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

Context. We present a method for including gas extinction of cosmic-ray-generated UV photons in chemical models of the midplane of protoplanetary disks, focusing on its implications on ice formation and chemical evolution. Aims. Our goal is to improve on chemical models by treating cosmic rays, the main source of ionization in the midplane of the disk, in a way that is consistent with current knowledge of the gas and grain environment present in those regions. We trace the effects of cosmic rays by identifying the main chemical reaction channels and also the main contributors to the gas opacity to cosmic-ray-induced UV photons. This information is crucial in implementing gas opacities for cosmic-ray-induced reactions in full 2D protoplanetary disk models. Methods. We considered time-dependent chemical models within the range 1-10 AU in the midplane of a T Tauri disk. The extinction of cosmic-ray-induced UV photons by gaseous species was included in the calculation of photorates at each timestep. We integrated the ionization and dissociation cross sections of all atoms/molecules over the cosmic-ray-induced UV emission spectrum of H_2. By analyzing the relative contribution of each gas phase species over time, we were able to identify the main contributors to the gas opacity in the midplane of protoplanetary disks. Results. At 1 AU the gas opacity contributes up to 28.2% of the total opacity, including the dust contribution. At 3-5 AU the gas contribution is 14.5% of the total opacity, and at 7-8 AU it reaches a value of 12.2%. As expected, at 10-15 AU freeze-out of species causes the gas contribution to the total opacity to be very low (6%). The main contributors to the gas opacity are CO, CO_2, S, SiO, and O_2. OH also contributes to the gas opacity, but only at 10-15 AU.

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)

Observations of the atomic and molecular line emission associated with jets and outflows emitted by young stellar objects can be used to trace the various evolutionary stages they pass through as they evolve to become main sequence stars. To understand the relevance of atomic and molecular cooling in shocks, and how accretion and ejection efficiency evolves with the source evolutionary state, we will study the far-infrared counterparts of bright optical jets associated with Class I and II sources in Taurus (T Tau, DG Tau A, DG Tau B, FS Tau A+B, and RW Aur). We have analysed Herschel/PACS observations of a number of atomic ([OI]63um, 145um, [CII]158um) and molecular (high-J CO, H2O, OH) lines, collected within the OTKP GASPS. To constrain the origin of the detected lines we have compared the FIR emission maps with the emission from optical-jets and millimetre-outflows, and the line fluxes and ratios with predictions from shock and disk models. All of the targets are associated with extended emission in the atomic lines correlated with the direction of the optical jet/mm-outflow. The atomic lines can be excited in fast dissociative J-shocks. The molecular emission, on the contrary, originates from a compact region, that is spatially and spectrally unresolved. Slow C- or J- shocks with high pre-shock densities reproduce the observed H2O and high-J CO lines; however, the disk and/or UV-heated outflow cavities may contribute to the emission. While the cooling is dominated by CO and H2O lines in Class 0 sources, [OI] becomes an important coolant as the source evolves and the environment is cleared. The cooling and mass loss rates estimated for Class II and I sources are one to four orders of magnitude lower than for Class 0 sources. This provides strong evidence to indicate that the outflow activity decreases as the source evolves.

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

Context. Time dependent gas-grain chemistry can help us understand the layered structure of species deposited onto the surface of grains during the lifetime of a protoplanetary disk. The history of trapping important quantities of carbon- and oxygen-bearing molecules onto the grains is of special significance for the formation of more complex (organic) molecules on the surface of grains. Aims. Among other processes, cosmic ray-induced UV photo-processes can lead to the efficient formation of OH. Using a more accurate treatment of cosmic ray-gas interactions for disks, we obtain an increased cosmic ray-induced UV photon flux of 3.8x10^5 photons cm^-2s^-1 for a cosmic-ray ionization rate of H2 value of 5x10^-17 s^-1 (compared to previous estimates of 10^4 photons cm^-2s^-1 based on ISM dust properties). We explore the role of the enhanced OH abundance on the gas-grain chemistry in the midplane of the disk at 10 AU, which is a plausible location for comet formation. We focus on studying the formation/destruction pathways and timescales of the dominant chemical species. Methods. We solve the chemical rate equations based on a gas-grain chemical network and correcting for the enhanced cosmic rayinduced UV field. This field is estimated from an appropriate treatment of dust properties in a protoplanetary disk, as opposed to previous estimates that assume an ISM-like grain size distribution. We also explore the chemical effects of photo-desorption of water ice into OH+H. Results. Near the end of the disk's lifetime our chemical model yields H2O, CO, CO2 and CH4 ice abundances at 10 AU (consistent with a midplane density of 10^10 cm^-3 and a temperature of 20 K) that are compatible with measurements of the chemical composition of cometary bodies for a [C/O] ratio of 0.16. Such comparison provides constraints on the physical conditions in which comets were formed.

We present a data model describing the structure of spectrophotometric datasets with spectral and temporal coordinates and associated metadata. This data model may be used to represent spectra, time series data, segments of SED (Spectral Energy Distributions) and other spectral or temporal associations.

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.

The data archives from space and ground-based telescopes present a vast opportunity for the astronomical community. We describe a classification encoding system for stellar spectra designed for archival databases that organizes the spectral data by "spectral classes." These classes are encoded into a digital format of the form TT.tt.LL.PPPP, where TT and tt refer to spectral type and subtype, LL to luminosity class, and PPPP to possible spectral peculiarities. Archive centers may wish to utilize this system to quantify classes of formerly arbitrary spectral classification strings found in classification catalogs corresponding to datasets of pointed spectroscopic observations in their holdings. The encoding system will also allow users to request archived data based on spectral class ranges, thereby streamlining an otherwise tedious data discovery process. Material in Appendix A is "normative" (part of the defined standard). Appendices B and C are "informative," meant to show how one data provider (MAST) has opted to handle some practical details.

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

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

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

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.

Depending on whom you ask, PAHs are either the smallest dust particles or the largest gas-phase molecules in space. Whether referred to as gas or dust, these PAHs can contain up to 20% of the total cosmic carbon abundance and as such also play an important role in the carbon chemistry of protoplanetary disks. The interpretation of PAH bands is often a complex procedure involving not only gas physics to determine their ionization stage and temperature, but also radiative transfer effects that can bury these bands in a strong thermal continuum from a population of larger dust particles. PAHs are most readily seen in the spectral energy distributions (SEDs) of disks around Herbig AeBe stars where they are photoprocessed by the stellar radiation field. Resolved images taken in the PAH bands confirm their origin in the flaring surfaces of circumstellar disks: if the SED is consistent with a flat disk structure (less illuminated), there is little or no evidence of PAH emission. The very low detection rates in the disks around T Tauri stars often require an overall lower abundance of PAHs in these disk surface as compared to that in molecular clouds. In this review, I will adress three aspects of PAHs in protoplanetary disks: (1) Do PAHs form in protoplanetary disks or do they originate from the precursor molecular cloud? (2) Is the presence of PAH features in SEDs a consequence of the disk structure or do PAHs in fact shape the disk structure? (3) How can we use PAHs as tracers of processes in protoplanetary disks?

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.

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

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

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

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 present the first far-IR observations of the solar-type stars delta Pav, HR 8501, 51 Peg and zeta^2 Ret, taken within the context of the DUNES Herschel Open Time Key Programme (OTKP). This project uses the PACS and SPIRE instruments with the objective of studying infrared excesses due to exo-Kuiper belts around nearby solar-type stars. The observed 100 um fluxes from delta Pav, HR 8501, and 51 Peg agree with the predicted photospheric fluxes, excluding debris disks brighter than Ldust/Lstar ~ 5 x 10^-7 (1 sigma level) around those stars. A flattened, disk-like structure with a semi-major axis of ~ 100 AU in size is detected around zeta^2 Ret. The resolved structure suggests the presence of an eccentric dust ring, which we interpret as an exo-Kuiper belt with Ldust/Lstar ~ 10^-5.

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 use of hydrodynamical simulations, the selection of atomic data, and the computation of deviations from local thermodynamical equilibrium for the analysis of the solar spectra have implied a downward revision of the solar metallicity. We are in the process of using the latest simulations computed with the CO5BOLD code to reassess the solar chemical composition. We determine the solar photospheric carbon abundance by using a radiation-hydrodynamical CO5BOLD model, and compute the departures from local thermodynamical equilibrium by using the Kiel code. We measure equivalent widths of atomic CI lines on high resolution, high signal-to-noise ratio solar atlases. Deviations from local thermodynamic equilibrium are computed in 1D with the Kiel code. Our recommended value for the solar carbon abundance, relies on 98 independent measurements of observed lines and is A(C)=8.50+-0.06, the quoted error is the sum of statistical and systematic error. Combined with our recent results for the solar oxygen and nitrogen abundances this implies a solar metallicity of Z=0.0154 and Z/X=0.0211. Our analysis implies a solar carbon abundance which is about 0.1 dex higher than what was found in previous analysis based on different 3D hydrodynamical computations. The difference is partly driven by our equivalent width measurements (we measure, on average, larger equivalent widths with respect to the other work based on a 3D model), in part it is likely due to the different properties of the hydrodynamical simulations and the spectrum synthesis code. The solar metallicity we obtain from the CO5BOLD analyses is in slightly better agreement with the constraints of helioseismology than the previous 3D abundance results. (Abridged)

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.

We present a new determination of the solar nitrogen abundance making use of 3D hydrodynamical modelling of the solar photosphere, which is more physically motivated than traditional static 1D models. We selected suitable atomic spectral lines, relying on equivalent width measurements already existing in the literature. For atmospheric modelling we used the co 5 bold 3D radiation hydrodynamics code. We investigated the influence of both deviations from local thermodynamic equilibrium (non-LTE effects) and photospheric inhomogeneities (granulation effects) on the resulting abundance. We also compared several atlases of solar flux and centre-disc intensity presently available. As a result of our analysis, the photospheric solar nitrogen abundance is A(N) = 7.86 +/- 0.12.

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.

CONTEXT: In recent years, the solar chemical abundances have been studied in considerable detail because of discrepant values of solar metallicity inferred from different indicators, i.e., on the one hand, the "sub-solar" photospheric abundances resulting from spectroscopic chemical composition analyses with the aid of 3D hydrodynamical models of the solar atmosphere, and, on the other hand, the high metallicity inferred by helioseismology. AIMS: After investigating the solar oxygen abundance using a CO5BOLD 3D hydrodynamical solar model in previous work, we undertake a similar approach studying the solar abundance of nitrogen, since this element accounts for a significant fraction of the overall solar metallicity, Z. METHOD: We used a selection of atomic spectral lines to determine the solar nitrogen abundance, relying mainly on equivalent width measurements in the literature. We investigate the influence on the abundance analysis, of both deviations from local thermodynamic equilibrium ("NLTE effects") and photospheric inhomogeneities ("granulation effects"). RESULTS: We recommend use of a solar nitrogen abundance of A(N)=7.86+-0.12 whose error bar reflects the line-to-line scatter. CONCLUSION: The solar metallicity implied by the CO5BOLD-based nitrogen and oxygen abundances is in the range 0.0145<= Z <= 0.0167. This result is a step towards reconciling photospheric abundances with helioseismic constraints on Z. Our most suitable estimates are Z=0.0156 and Z/X=0.0213.

We present a model that describes stellar infrared excesses due to heating of the interstellar (IS) dust by a hot star passing through a diffuse IS cloud. This model is applied to six lambda Bootis stars with infrared excesses. Plausible values for the IS medium (ISM) density and relative velocity between the cloud and the star yield fits to the excess emission. This result is consistent with the diffusion/accretion hypothesis that lambda Bootis stars (A- to F-type stars with large underabundances of Fe-peak elements) owe their characteristics to interactions with the ISM. This proposal invokes radiation pressure from the star to repel the IS dust and excavate a paraboloidal dust cavity in the IS cloud, while the metal-poor gas is accreted onto the stellar photosphere. However, the measurements of the infrared excesses can also be fit by planetary debris disk models. A more detailed consideration of the conditions to produce lambda Bootis characteristics indicates that the majority of infrared-excess stars within the Local Bubble probably have debris disks. Nevertheless, more distant stars may often have excesses due to heating of interstellar material such as in our model.

One of the fundamental astrobiology questions is how life has formed in our Solar System. In this context the formation and stability of abiotic organic molecules such as CH4, formic acid and amino acids, is important for understanding how organic material has formed and survived shocks and energetic particle impact from winds in the early Solar System. Shock waves have been suggested as a plausible scenario to create chondrules, small meteoritic components that have been completely molten by energetic events such as shocks and high velocity particle impacts. We study here the formation and destruction of certain gas-phase molecules such as methane and water during such shock events and compare the chemical timescales with the timescales for shocks arising from gravitational instabilities in a protosolar nebula.

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.

The A star 49 Ceti, at a distance of 61 pc, is unusual in retaining a substantial quantity of molecular gas while exhibiting dust properties similar to those of a debris disk. We present resolved observations of the disk around 49 Ceti from the Submillimeter Array in the J=2-1 rotational transition of CO with a resolution of 1.0x1.2 arcsec. The observed emission reveals an extended rotating structure viewed approximately edge-on and clear of detectable CO emission out to a distance of ~90 AU from the star. No 1.3 millimeter continuum emission is detected at a 3-sigma sensitivity of 2.1 mJy/beam. Models of disk structure and chemistry indicate that the inner disk is devoid of molecular gas, while the outer gas disk between 40 and 200 AU from the star is dominated by photochemistry from stellar and interstellar radiation. We determine parameters for a model that reproduces the basic features of the spatially resolved CO J=2-1 emission, the spectral energy distribution, and the unresolved CO J=3-2 spectrum. We investigate variations in disk chemistry and observable properties for a range of structural parameters. 49 Ceti appears to be a rare example of a system in a late stage of transition between a gas-rich protoplanetary disk and a tenuous, virtually gas-free debris disk.

We have discovered a bow shock shaped mid-infrared excess region in front of delta Velorum using 24 micron observations obtained with the Multiband Imaging Photometer for Spitzer (MIPS). The excess has been classified as a debris disk from previous infrared observations. Although the bow shock morphology was only detected in the 24 micron observations, its excess was also resolved at 70 micron. We show that the stellar heating of an ambient interstellar medium (ISM) cloud can produce the measured flux. Since delta Velorum was classified as a debris disk star previously, our discovery may call into question the same classification of other stars. We model the interaction of the star and ISM, producing images that show the same geometry and surface brightness as is observed. The modeled ISM is 15 times overdense relative to the average Local Bubble value, which is surprising considering the close proximity (24 pc) of delta Velorum. The abundance anomalies of lambda Bootis stars have been previously explained as arising from the same type of interaction of stars with the ISM. Low resolution optical spectra of delta Velorum show that it does not belong to this stellar class. The star therefore is an interesting testbed for the ISM accretion theory of the lambda Bootis phenomenon.

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

We present the results of a coronagraphic scattered-light imaging survey of six young disk candidate stars using the Hubble Space Telescope Advanced Camera for Surveys. The observations made use of the 1.8" occulting spot through the F606W (broad V) filter. Circumstellar material was imaged around HD 97048, a Herbig Ae/Be star located in the Chamaeleon I dark cloud at a distance of 180 pc. The material is seen between ~2" (360 AU) and ~4" (720 AU) from the star in all directions. A V-band azimuthally-averaged radial surface brightness profile peaks at r = 2" with a value of 19.6 +/- 0.2 mag arcsec^-2 and smoothly decreases with projected distance from the star as r^(-3.3 +/- 0.5). An integrated flux of 16.8 +/- 0.1 mag is measured between 2" and 4", corresponding to a scattered-light fractional luminosity lower limit of 8.4 x 10^-4. Filamentary structure resembling spiral arms similar to that seen in Herbig Ae/Be disks is observed. Such structure has been attributed to the influence of orbiting planets or stellar encounters. Average surface brightness upper limits are determined for the five non-detections: HD 34282, HD 139450, HD 158643, HD 159492, and HD 195627. Possible reasons for the non-detections are disks that are too faint or disks hidden by the occulter.

We present models for the HI 21 cm emission from circumstellar disks and use them to convert observed upper limits on the HI 21 cm flux to limits on the total disk gas mass. The upper limits we use come from earlier Australia Telescope Compact Array observations of the debris disk around beta Pictoris as well as fresh Giant Meterwave Radio Telescope observations of HD135344, LkCa 15 and HD163296. Our observations and models span a range of disk types, from young proto-planetary disks to old debris disks. The models self-consistently calculate the gas chemistry (H/H2 balance) and the thermal structure of UV irradiated disks. Atomic hydrogen production is dominated by UV irradiation in transition phase objects as well as debris disks, but for very young disks, HI production by stellar X-rays (which we do not account for) is important. We use a simple radiative transfer approach to convert the model disk parameters into predicted HI 21 cm line maps and spectral profiles. This allows a direct comparison of the observations to the model. We find that the HI traces the disk surface layers, and that the Hi emission, if detected, could be used to study the e ects of irradiation and evaporation, in addition to the kinematics of the disk. Our models cover massive protoplanetary disks, transition phase disks and dusty debris disks. In massive protoplanetary disks, UV produced HI constitutes less than 0.5% of the total disk mass, while X-rays clearly dominate the chemistry and thus the HI production. For the two such disks that we have observed, viz. those around LkCa 15 and HD163296, the predicted 21 cm flux is below the current detection limit. On the other hand, transition phase disks at distances of 100 pc have predicted 21 cm fluxes that are close to the detection limit.... (abstract abbreviated)

Aims: The chemistry, distribution and mass of the gas in the transitional disk around the 5 Myr old B9.5 V star HD 141569A are constrained. Methods: A quasi 2-dimensional (2D) chemistry code for photon dominated regions (PDR) is used to calculate the chemistry and gas temperatures in the disk. The calculations are performed for several gas distributions, PAH abundances and values of the total gas mass. The resulting CO J=2-1 and J=3-2 emission lines are computed with a 2D radiative transfer code and are compared to observations. Results: The CO abundance is very sensitive to the total disk mass because the disk is in a regime where self-shielding just sets in. The observed CO emission lines are best fit by a power-law gas distribution of 80 M_earth starting at 80 AU from the central star, indicating that there is some gas in the inner hole. Predictions are made for intensities of atomic fine-structure lines. [C I], which is the dominant form of carbon in large parts of the disk, is found to be a good alternative tracer of the gas mass.

We review advances in the modeling of protoplanetary disks. This review will focus on the regions of the disk beyond the dust sublimation radius, i.e. beyond 0.1 - 1 AU, depending on the stellar luminosity. We will be mostly concerned with models that aim to fit spectra of the dust continuum or gas lines, and derive physical parameters from these fits. For optically thick disks, these parameters include the accretion rate through the disk onto the star, the geometry of the disk, the dust properties, the surface chemistry and the thermal balance of the gas. For the latter we are mostly concerned with the upper layers of the disk, where the gas and dust temperature decouple and a photoevaporative flow may originate. We also briefly discuss optically thin disks, focusing mainly on the gas, not the dust. The evolution of these disks is dominated by accretion, viscous spreading, photoevaporation, and dust settling and coagulation. The density and temperature structure arising from the surface layer models provide input to models of photoevaporation, which occurs largely in the outer disk. We discuss the consequences of photoevaporation on disk evolution and planet formation.

We present a short atlas illustrating the unusual Ca {\sc ii} K-line profiles in upper main sequence stars with anomalous abundances. Slopes of the profiles for 10 cool, magnetic chemically peculiar (CP) stars change abruptly at the very core, forming a deep "nib." The nibs show the same or nearly the same radial velocity as the other atomic lines. The near wings are generally more shallow than in normal stars. In three magnetic CP stars, the K-lines are too weak to show this shape, though the nibs themselves are arguably present. The Ca {\sc ii} H-lines also show deep nibs, but the profiles are complicated by the nearby, strong H$\epsilon$ absorption. The K-line structure is nearly unchanged with phase in $\beta$ CrB and $\alpha$ Cir. Calculations, including NLTE, show that other possibilities in addition to chemical stratification may yield nib-like cores.

The structure and chemistry of protoplanetary disks depends strongly on the nature of the central star around which it has formed. The dust temperature is mainly set by the stellar luminosity, while the chemistry of the upper disk layers depends on the amount of intercepted UV and X-ray flux. We will study the differences in chemistry, thermal sturcture and line emission around Herbig Ae/Be, T Tauri stars and low mass M dwarfs. Predictions will be made for future observations with SOFIA and Herschel.

The radiation hydrodynamic code CO5BOLD has been supplemented with the time-dependent treatment of chemical reaction networks. Advection of particle densities due to the hydrodynamic flow field is also included. The radiative transfer is treated frequency-independently, i.e. grey, so far. The upgraded code has been applied to two-dimensional simulations of carbon monoxide (CO) in the non-magnetic solar photosphere and low chromosphere. For this purpose a reaction network has been constructed, taking into account the reactions which are most important for the formation and dissociation of CO under the physical conditions of the solar atmosphere. The network has been strongly reduced to 27 reactions, involving the chemical species H, H2, C, O, CO, CH, OH, and a representative metal. The resulting CO number density is highest in the cool regions of the reversed granulation pattern at mid-photospheric heights and decreases strongly above. There, the CO abundance stays close to a value of 8.3 on the usual logarithmic abundance scale with [H]=12 but is reduced in hot shock waves which are a ubiquitous phenomenon of the model atmosphere. For comparison, the corresponding equilibrium densities have been calculated, based on the reaction network but also under assumption of instantaneous chemical equilibrium by applying the Rybicki & Hummer (RH) code by Uitenbroek (2001). Owing to the short chemical timescales, the assumption holds for a large fraction of the atmosphere, in particular the photosphere. In contrast, the CO number density deviates strongly from the corresponding equilibrium value in the vicinity of chromospheric shock waves. Simulations with altered reaction network clearly show that the formation channel via hydroxide (OH) is the most important one under the conditions of the solar atmosphere.

In the past year, the Multimission Archive at the Space Telescope Science Institute (MAST) has taken major steps in making MAST's holdings available using VO-defined protocols and standards, and in implementing VO-based tools. For example, MAST has implemented the Simple Cone Search protocol, and all MAST mission searches may be returned in the VOTable format, allowing other archives to use MAST data for their VO applications. We have made many of our popular High Level Science Products available through Simple Image Access Protocol (SIAP), and are implementing the VO Simple Spectral Access Protocol (SSAP). The cross correlation of VizieR catalogs with MAST missions is now possible, and illustrates the integration of VO services into MAST. The user can easily display the results from searches within MAST using the plotting tool VOPlot. MAST also participates in the NVO registry service. Thus, the user can harvest MAST holdings simultaneously with data from many other surveys and missions through the VO DataScope Data Inventory Service.

Models for the structure of protoplanetary disks have so far been based on the assumption that the gas and the dust temperature are equal. The gas temperature, an essential ingredient in the equations of hydrostatic equilibrium of the disk, is then determined from a continuum radiative transfer calculation, in which the continuum opacity is provided by the dust. It has been long debated whether this assumption still holds in the surface layers of the disk, where the dust infrared emission features are produced. In this paper we compute the temperature of the gas in the surface layers of the disk in a self-consistent manner. The gas temperature is determined from a heating-cooling balance equation in which processes such as photoelectric heating, dissociative heating, dust-gas thermal heat exchange and line cooling are included. The abundances of the dominant cooling species such as CO, C, C+ and O are determined from a chemical network based on the atomic species H, He, C, O, S, Mg, Si, Fe (Kamp & Bertoldi 2000). The underlying disk models to our calculations are the models of Dullemond, van Zadelhoff & Natta (2002). We find that in general the dust and gas temperature are equal to withing 10% for A_V >~ 0.1, which is above the location of the `super-heated surface layer' in which the dust emission features are produced (e.g. Chiang & Goldreich 1997). High above the disk surface the gas temperature exceeds the dust temperature and can can become -- in the presence of polycyclic aromatic hydrocarbons -- as high as 600 K at a radius of 100 AU. This is a region where CO has fully dissociated, but a significant fraction of hydrogen is still in molecular form. The densities are still high enough for non-negligible H_2 emission to be produced.....(see paper for full abstract)

Interpretation of molecular line observations in tenuous circumstellar disks around young G-type stars in terms of a disk mass is difficult without a model that describes the chemical structure of these disks. This paper now discusses the chemistry in tenuous disks around young solar-type stars based on disk models that take into account the presence of a stellar chromosphere. The example of the disk around a 70 Myr old solar-type star shows that the dissociating radiation from the chromosphere is stronger than the interstellar ultraviolet radiation field up to a distance of ~400 AU from the star. Similar to other studies in this research field, it is found that, due to photodissociation, the CO-to-H_2 ratio is far from the canonical value of 10^-4 for molecular clouds. Moreover, the dust-to-gas mass ratio, as well as the dust grain size play an important role for the H_2 abundance in these disks.

We have obtained FUSE and HST STIS spectra of HR 4796A, a nearby 8 Myr old main sequence star that possesses a dusty circumstellar disk whose inclination has been constrained from high resolution near-infrared observations to be ~17 deg from edge-on. We searched for circumstellar absorption in the ground states of C II at 1036.3 A, O I at 1039.2 A, Zn II at 2026.1 A, Lyman series H2, and CO (A-X) and failed to detect any of these species. We place upper limits on the column densities and infer upper limits on the gas masses assuming that the gas is in hydrostatic equilibrium, is well-mixed, and has a temperature, Tgas ~ 65 K. Our measurements suggest that this system possesses very little molecular gas. Therefore, we infer an upper limit for the gas:dust ratio (<4.0) assuming that the gas is atomic. We measure less gas in this system than is required to form the envelope of Jupiter.

The group of lambda Bootis type stars comprises late B- to early F-type, Population I objects which are basically metal weak, in particular the Fe group elements, but with the clear exception of C, N, O and S. One of the theories to explain the abundance pattern of these stars involves circumstellar or interstellar matter around the objects. Hence, we have compiled all available data from the literature of well established members of the lambda Bootis group redward of 7000A in order to find evidence for matter around these objects. Furthermore, we present unpublished ISO as well as submillimeter continuum and CO (2-1) line measurements to complete the data set. In total, measurements for 34 (26 with data redward of 20 mu m) well established lambda Bootis stars are available. There is evidence for an infrared excesses in six stars (HD 31295, HD 74873, HD 110411, HD 125162, HD 198160/1 and HD 210111) and two are doubtful cases (HD 11413 and HD 192640) resulting in a percentage of 23% (excluding the two doubtful cases). Dust models for these objects show fractional dust luminosities comparable to the Vega-type stars and slightly higher dust temperatures. ISO-SWS spectroscopy for HD 125162 and HD 192640 resulted in the detection of pure stellar HI lines ruling out an active accretion disk (as found for several Herbig Ae/Be stars) around these objects. The submillimeter measurements gave only upper limits for the line and continuum fluxes.

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

The lambda Bootis stars are Population I, late B to early F-type stars, with moderate to extreme (up to a factor 100) surface underabundances of most Fe-peak elements and solar abundances of lighter elements (C, N, O, and S). To put constraints on the various existing theories that try to explain these peculiar stars, we investigate the observational properties of lambda Bootis stars compared to a reference sample of normal stars. Using various photometric systems and Hipparcos data, we analyze the validity of standard photometric calibrations, elemental abundances, and Galactic space motions. There crystallizes a clear picture of a homogeneous group of Population I objects found at all stages of their main-sequence evolution, with a peak at about 1 Gyr. No correlation of astrophysical parameters such as the projected rotational velocities or elemental abundances with age is found, suggesting that the a-priori unknown mechanism, which creates lambda Bootis stars, works continuously for late B to early F-type stars in all stages of main-sequence evolution. Surprisingly, the sodium abundances seem to indicate an interaction between the stars and their local environment.

The lambda Bootis stars, a group of late B to early F-type population I stars, have surface abundances that resemble the general metal depletion pattern found in the interstellar medium. Inspired by the recent result that the fundamental parameters of these peculiar stars differ in no respect from a comparison sample of normal stars, the hypothesis of an interaction between a star and a diffuse interstellar cloud is considered as a possible explanation of the peculiar abundance pattern. It is found that such a scenario is able to explain the selective accretion of interstellar gas depleted in condensable elements as well as the spectral range of the lambda Bootis phenomenon.

We determined abundances of O, Ca, Fe, Ba and Y for a sample of dusty and dust-free A stars, taken from the list of Cheng et al. (1992). Five of the stars have an infrared-excess due to circumstellar dust. Ongoing accretion from their circumstellar surroundings might have modified the abundances in the photospheres of these stars, but our results clearly show, that there is no difference in the photospheric composition of the dusty and dust-free stars. Instead all of them show the typical diffusion pattern which diminishes towards larger rotational velocities.

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

High-resolution spectra of lambda Bootis stars reveal the presence of circumstellar gas for example in the Ca K line. The example of the normal A star beta Pictoris shows, that the narrow stable absorption component in Ca K can be reproduced using appropriate disk models and a calcium underabundance in the circumstellar gas of a factor of ~30. Similar models are suggested for the group of metal-deficient lambda Bootis stars, but the observational material is still very poor.

We present preliminary results of our search for circumstellar absorption features in the Ca K lines based on high S/N observations obtained with the ESO CAT/CES system.