We observed a sample of 20 representative Herbig Ae/Be stars and 5 A-type debris discs with PACS onboard Herschel, as part of the GAS in Protoplanetary Systems (GASPS) project. The observations were ...done in spectroscopic mode, and cover the far-infrared lines of O i, C ii, CO, CH+, H2O, and OH. We have a O i 63 μm detection rate of 100% for the Herbig Ae/Be and 0% for the debris discs. The O i 145 μm line is only detected in 25% and CO J = 18-17 in 45% (and fewer cases for higher J transitions) of the Herbig Ae/Be stars, while for C ii 157 μm, we often find 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. First seen in HD 100546, CH+ emission 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 either the stellar or disc parameters, such as stellar luminosity, ultraviolet and X-ray flux, accretion rate, polycyclic aromatic hydrocarbon (PAH) band strength, and flaring. We find that the stellar ultraviolet flux is the dominant excitation mechanism of O i 63 μm, with the highest line fluxes being found in objects with a large amount of flaring and among the largest PAH strengths. 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 O i 63 μm and O i 145 μm, CO J = 18-17 and O i 6300 Å, and between the continuum flux at 63 μm and at 1.3 mm, while we find weak correlations between the line flux of O i 63 μm and the PAH luminosity, the line flux of CO J = 3-2, the continuum flux at 63 μm, the stellar effective temperature, and the Brγ luminosity. Finally, we use a combination of theO i 63 μm and 12CO J = 2-1 line fluxes to obtain order of magnitude estimates of the disc gas masses, in agreement with the values that we find from detailed modelling of two Herbig Ae/Be stars, HD 163296 and HD 169142. Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.
The gas- and dust dissipation processes in disks around young stars remain uncertain despite numerous studies. At the distance of similar99-116 pc, HD 141569A is one of the nearest Herbig Ae stars ...that is surrounded by a tenuous disk, probably in transition between a massive primordial disk and a debris disk. Atomic and molecular gases have been found in the structured 5-Myr old HD 141569A disk, making HD 141569A the perfect object within which to directly study the gaseous atomic and molecular component. We wish to constrain the gas and dust mass in the disk around HD 141569A. We constrained simultaneously the silicate dust grain, PAH, and gas mass in a similar5-Myr old Herbig Ae disk. The disk-averaged gas-to-dust-mass is most likely around 100, which is the assumed value at the disk formation despite the uncertainties due to disagreements between the different gas tracers. If the disk was originally massive, the gas and the dust would have dissipated at the same rate.
Context. Studying the physical conditions in circumstellar disks is a crucial step toward understanding planet formation and disk evolution. Of particular interest is the case of HD 100546, a Herbig ...Be star that presents a gap within the first 13 AU of its protoplanetary disk, a gap that may originate in the dynamical interactions of a forming planet with its hosting disk. Aims. We seek a more detailed understanding of the structure of the circumstellar environment of HD 100546 and refine our previous disk model that is composed of a tenuous inner disk, a gap and a massive outer disk (see Benisty et al. 2010, A&A, 511, A75). We also investigate whether planetary formation processes can explain the complex density structure observed in the disk. Methods. We gathered a large amount of new interferometric data using the AMBER/VLTI instrument in the H- and K-bands to spatially resolve the warm inner disk and constrain its structure. Then, combining these measurements with photometric observations, we analyze the circumstellar environment of HD 100546 in the light of a passive disk model based on 3D Monte-Carlo radiative transfer. Finally, we use hydrodynamical simulations of gap formation by planets to predict the radial surface density profile of the disk and test the hypothesis of ongoing planet formation. Results. The SED (spectral energy distribution) from the UV to the millimeter range, and the NIR (near-infrared) interferometric data are adequately reproduced by our model. We show that the H- and K-band emissions are coming mostly from the inner edge of the internal dust disk, located near 0.24 AU from the star, i.e., at the dust sublimation radius in our model. At such a short distance, the survival of hot (silicate) dust requires the presence of micron-sized grains, heated at ~1750 K. We directly measure an inclination of 33° ± 11° and a position angle of 140° ± 16° for the inner disk. This is similar to the values found for the outer disk (i ≃ 42°, PA ≃ 145°), suggesting that both disks may be coplanar. We finally show that 1 to 8 Jupiter mass planets located at ~8 AU from the star would have enough time to create the gap and the required surface density jump of three orders of magnitude between the inner and outer disk. However, no information on the amount of matter left in the gap is available, which precludes us from setting precise limits on the planet mass, for now.
Context. Constraining the gas and dust disk structure of transition disks, particularly in the inner dust cavity, is a crucial step toward understanding the link between them and planet formation. HD ...135344B is an accreting (pre-)transition disk that displays the CO 4.7 μm emission extending tens of AU inside its 30 AU dust cavity. Aims. We constrain HD 135344B’s disk structure from multi-instrument gas and dust observations. Methods. 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 μm, Herschel/PACS O i 63 μm, 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. Results. We found a disk model able to describe the current gas and dust observations simultaneously. This disk has the following structure. (1) To simultaneously 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 the gas inside the cavity must increase with radius to fit the CO ro-vibrational line profile, a small gap of a few AU in the gas distribution is compatible with current data, and a large gap of tens of AU in the gas does not appear likely. (4) The gas-to-dust ratio inside the cavity is >100 to account for the 870 μm continuum upper limit and the CO P(10) line flux. (5) The gas-to-dust ratio in the outer disk (30 <R< 200 AU) is <10 to simultaneously describe the O i 63 μm line flux and the CO P(10) line profile. (6) In the outer disk, most of the gas and dust mass should be located in the midplane, and a significant fraction of the dust should be 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 cavity echoes the effect of a migrating jovian planet in the disk structure. The low gas mass (a few Jupiter masses) throughout the HD 135344B disk supports the idea that it is an evolved disk that has already lost a large portion of its mass.
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.
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 the 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 μm 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\hbox{$^{+2320}_{-151}$}+2320-151 K with a mass of ~3 × 10-14−5 × 10-12 M⊙. 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, which is consistent with previous observations of hot CO gas emission.
ABSTRACT We describe a large-scale far-infrared line and continuum survey of protoplanetary disk through to young debris disk systems carried out using the ACS instrument on the Herschel Space ...Observatory. This Open Time Key program, known as GASPS (Gas Survey of Protoplanetary Systems), targeted ∼250 young stars in narrow wavelength regions covering the OI fine structure line at 63 μm the brightest far-infrared line in such objects. A subset of the brightest targets were also surveyed in OI145 μm, CII at 157 μm, as well as several transitions of H2O and high-excitation CO lines at selected wavelengths between 78 and 180 μm. Additionally, GASPS included continuum photometry at 70, 100 and 160 μm, around the peak of the dust emission. The targets were SED Class II-III T Tauri stars and debris disks from seven nearby young associations, along with a comparable sample of isolated Herbig AeBe stars. The aim was to study the global gas and dust content in a wide sample of circumstellar disks, combining the results with models in a systematic way. In this overview paper we review the scientific aims, target selection and observing strategy of the program. We summarise some of the initial results, showing line identifications, listing the detections, and giving a first statistical study of line detectability. The OI line at 63 μm was the brightest line seen in almost all objects, by a factor of ∼10. Overall OI63 μm detection rates were 49%, with 100% of HAeBe stars and 43% of T Tauri stars detected. A comparison with published disk dust masses (derived mainly from sub-mm continuum, assuming standard values of the mm mass opacity) shows a dust mass threshold for OI63 μm detection of ∼10-5 M⊙. Normalising to a distance of 140 pc, 84% of objects with dust masses ≥10-5 M⊙ can be detected in this line in the present survey; 32% of those of mass 10-6-10-5 M⊙, and only a very small number of unusual objects with lower masses can be detected. This is consistent with models with a moderate UV excess and disk flaring. For a given disk mass, OI detectability is lower for M stars compared with earlier spectral types. Both the continuum and line emission was, in most systems, spatially and spectrally unresolved and centred on the star, suggesting that emission in most cases was from the disk. Approximately 10 objects showed resolved emission, most likely from outflows. In the GASPS sample, OI detection rates in T Tauri associations in the 0.3-4 Myr age range were ∼50%. For each association in the 5-20 Myr age range, ∼2 stars remain detectable in OI63 μm, and no systems were detected in associations with age >20 Myr. Comparing with the total number of young stars in each association, and assuming a ISM-like gas/dust ratio, this indicates that ∼18% of stars retain a gas-rich disk of total mass ∼1 MJupiter for 1-4 Myr, 1-7% keep such disks for 5-10 Myr, but none are detected beyond 10-20 Myr. The brightest OI objects from GASPS were also observed in OI145 μm, CII157 μm and CO J = 18 - 17, with detection rates of 20-40%. Detection of the CII line was not correlated with disk mass, suggesting it arises more commonly from a compact remnant envelope.
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.
Context. Molecular hydrogen (H2) gas is the most abundant molecule in the circumstellar (CS) environments of young stars. It is thus a key element in our understanding of the evolution of pre-main ...sequence stars and their environments towards the main sequence. Aims. At the present time, little is known about the gas as compared to the dust in the environments of young stars. We thus observed molecular hydrogen around a sample of pre-main sequence stars in order to better characterize their gaseous CS environments. Methods. The FUSE (Far Ultraviolet Spectroscopic Explorer) spectral range offers access to electronic transitions of H2. We analyzed the FUSE spectra of a sample of Herbig Ae/Be stars (HAeBes) covering a broad spectral range (from F4 to B2), including the main-sequence A5 star β Pictoris. To better diagnose the origin of the detected molecular gas and its excitation conditions, we used a model of a photodissociation region. Results. Our analysis demonstrates that the excitation of H2 is clearly different around most of the HAeBes compared to the interstellar medium. Moreover, the characteristics of H2 around Herbig Ae and Be stars give evidence for different excitation mechanisms. For the most massive stars of our sample (B8 to B2 type), the excitation diagrams are reproduced well by a model of photodissociation regions (PDR). Our results favor an interpretation in terms of large CS envelopes, remnants of the molecular clouds in which the stars were formed. On the other hand, the group of Ae stars (later than B9 type) known to possess disks is more inhomogeneous. In most cases, when CS H2 is detected, the lines of sight do not pass through the disks. The excitation conditions of H2 around Ae stars cannot be reproduced by PDR models and correspond to warm and/or hot excited media very close to the stars. In addition, no clear correlation has been found between the ages of the stars and the amount of circumstellar H2. Our results suggest structural differences between Herbig Ae and Be star environments. Herbig Be stars do evolve faster than Ae stars, and consequently, most Herbig Be stars are younger than Ae ones at the time we observe them. It is thus more likely to find remnants of their parent cloud around them.
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.