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 ... ~540 AU. 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 63 μm 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 an upper limit for the H2 0–0 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.
Aims: We want to understand the chemistry and physics of discs on the basis of a large unbiased and statistically relevant grid of disc models. One of the main goals is to explore the diagnostic ...power of various gas emission lines and line ratios for deriving main disc parameters such as the gas mass. Methods: We explored the results of the DENT grid (Disk Evolution with Neat Theory) that consists of 300 000 disc models with 11 free parameters. Through a statistical analysis, we searched for correlations and trends in an effort to find tools for disc 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 disc model. The statistical analysis of the DENT grid shows that CO gas is rarely the dominant carbon reservoir in discs. 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 discs 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. Our study of the observational diagnostics shows that the C ii 158 μm fine structure line flux is very sensitive to the stellar UV flux and presence of a UV excess, and that it traces the outer disc radius (Rout). In the submm, the CO low J rotational lines also trace Rout. Low O i 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 μm line. This occurs mostly for massive non-flaring, settled disc models without UV excess. A combination of the O i 63 line and low J CO lines correlates with several disc properties, such as the average O i gas temperature in discs, the outer disc radius, and the UV excess. Conclusions: The O i 63/CO 2-1 line ratio is a powerful diagnostic for breaking disc modelling degeneracies. A combination of the O i 63 μm flux and the O i 63/CO 2-1 line ratio can be used for Mgas ≤ 10-3 M⊙ to obtain an order of magnitude estimate for the disc gas mass purely from gas observations. The previously used conversion of a CO submm line flux alone generally leads to larger uncertainties.
Planets are formed in disks around young stars. With an age of similar to 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 (12) CO 3-2 and (CO)-C-13 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 mu m. The other lines that were observed, OI at 145 mu m and CII at 157 mu m, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming (CO)-C-12/(CO)-C-13 = 69 suggests a dust mass for grains with radius < 1 mm of similar to 1.9 x 10(-4) M-circle dot (total solid mass of 3 x 10(-3) M-circle dot) and a gas mass of (0.5-5) x 10(-3) M-circle dot. The gas-to-dust mass may be lower than the standard interstellar value of 100.
We have combined the thermo-chemical disc code ProDiMo with the Monte Carlo radiative transfer code MCFOST to calculate a grid of ∼300 000 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-infrared, sub-mm and mm lines of Oi, Cii, 12CO and o/p-H2O under five inclinations. The grid allows us 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 temperature and chemical structures, the computed spectral energy distributions, line fluxes and profiles, will be used in particular for the data interpretation of the HERSCHEL open time-key program GASPS. The calculated line fluxes show a strong dependence on the assumed ultraviolet 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 is calculated as a function of disc mass. The possibility of deriving the disc gas mass from line observations is discussed.
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 program 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 μm, CII 157.7 μm, CO 72.8 and 90.2 μm, and o-H2O 78.7 and 179.5 μm. We only detect the OI 63.2 μm 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 × 10-3 $M_{\odot}$ 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 approximate to 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 gimel similar to 60-190 mu m using the Photodetector Array Camera and Spectrometer instrument on the Herschel Space Observatory.Results. We detect OI 63 mu m 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 mu m. Extensive models are presented in associated papers.