Gas in the protoplanetary disc of HD 169142: Herschel’s view ^†^†thanks: Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.

G. Meeus Dep. de Física Teórica, Fac. de Ciencias, UAM Campus Cantoblanco, 28049 Madrid, Spain
gwendolyn.meeus@uam.es C. Pinte School of Physics, University of Exeter, Stocker Road, Exeter EX4 4QL, United Kingdom Université Joseph Fourier - Grenoble 1/CNRS, Laboratoire d’Astrophysique de Grenoble (LAOG) UMR 5571, BP 53, 38041 Grenoble Cedex 09, France P. Woitke UK Astronomy Technology Centre, Royal Observatory, Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, UK SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory Edinburgh, UK Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK School of Physics & Astronomy, University of St. Andrews, North Haugh, St. Andrews KY16 9SS, UK B. Montesinos LAEX, Depto. Astrofísica, Centro de Astrobiología (INTA-CSIC), P.O. Box 78, E-28691 Villanueva de la Cañada, Spain I. Mendigutía LAEX, Depto. Astrofísica, Centro de Astrobiología (INTA-CSIC), P.O. Box 78, E-28691 Villanueva de la Cañada, Spain P. Riviere-Marichalar LAEX, Depto. Astrofísica, Centro de Astrobiología (INTA-CSIC), P.O. Box 78, E-28691 Villanueva de la Cañada, Spain C. Eiroa Dep. de Física Teórica, Fac. de Ciencias, UAM Campus Cantoblanco, 28049 Madrid, Spain
gwendolyn.meeus@uam.es G.S. Mathews Institute for Astronomy (IfA), University of Hawaii, 2680 Woodlawn Dr., Honolulu, HI 96822, USA B. Vandenbussche Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, Leuven, Belgium C.D. Howard SOFIA-USRA, NASA Ames Research Center, MS 211-3, Moffett Field, CA 94035, USA A. Roberge NASA Goddard Space Flight Center, Exoplanets & Stellar Astrophysics lab, Code 667, Greenbelt, MD 20771, USA G. Sandell SOFIA-USRA, NASA Ames Research Center, MS 211-3, Moffett Field, CA 94035, USA G. Duchêne Université Joseph Fourier - Grenoble 1/CNRS, Laboratoire d’Astrophysique de Grenoble (LAOG) UMR 5571, BP 53, 38041 Grenoble Cedex 09, France Astronomy Department, University of California, Berkeley CA 94720-3411 USA F. Ménard Université Joseph Fourier - Grenoble 1/CNRS, Laboratoire d’Astrophysique de Grenoble (LAOG) UMR 5571, BP 53, 38041 Grenoble Cedex 09, France C.A. Grady Eureka ScientiÞc, 2452 Delmer, Suite 100, Oakland, CA 96002, USA NASA Goddard Space Flight Center, Exoplanets & Stellar Astrophysics lab, Code 667, Greenbelt, MD 20771, USA W.R.F. Dent ALMA, Avda Apoquindo 3846, Piso 19, Edificio Alsacia, Las Condes, Santiago, Chile European Southern Observatory, Alonso de Córdova 3107, Vitacura, Santiago, Chile I. Kamp Kapteyn Astronomical Institute, P.O. Box 800, 9700 AV Groningen, The Netherlands J.C. Augereau Université Joseph Fourier - Grenoble 1/CNRS, Laboratoire d’Astrophysique de Grenoble (LAOG) UMR 5571, BP 53, 38041 Grenoble Cedex 09, France W.F. Thi SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory Edinburgh, UK Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK Université Joseph Fourier - Grenoble 1/CNRS, Laboratoire d’Astrophysique de Grenoble (LAOG) UMR 5571, BP 53, 38041 Grenoble Cedex 09, France I. Tilling SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory Edinburgh, UK Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK J.M. Alacid Spanish Virtual Observatory, Unidad de Archivo de Datos, Depto. Astrofísica, Centro de Astrobiología (INTA-CSIC), P.O. Box 78, E-28691 Villanueva de la Cañada, Spain S. Andrews Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA, USA D.R. Ardila NASA Herschel Science Center, California Institute of Technology, Pasadena, USA G. Aresu Kapteyn Astronomical Institute, P.O. Box 800, 9700 AV Groningen, The Netherlands D. Barrado Calar Alto Observatory, Centro Astronómico Hispano-Alemán, c/ Jesús Durbán Remón, 2-2, 04004 Almería, Spain LAEX, Depto. Astrofísica, Centro de Astrobiología (INTA-CSIC), P.O. Box 78, E-28691 Villanueva de la Cañada, Spain S. Brittain Department of Physics & Astronomy, Clemson University, Clemson SC 29634-0978 USA D.R. Ciardi NASA Exoplanet Science Institute/Caltech 770 South Wilson Avenue, Mail Code: 100-22, Pasadena, CA USA 91125 W. Danchi NASA Goddard Space Flight Center, Exoplanets & Stellar Astrophysics lab, Code 667, Greenbelt, MD 20771, USA D. Fedele Dep. de Física Teórica, Fac. de Ciencias, UAM Campus Cantoblanco, 28049 Madrid, Spain
gwendolyn.meeus@uam.es Max Planck Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany; Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA I. de Gregorio-Monsalvo ALMA, Avda Apoquindo 3846, Piso 19, Edificio Alsacia, Las Condes, Santiago, Chile European Southern Observatory, Alonso de Córdova 3107, Vitacura, Santiago, Chile A. Heras Research and Scientific Support Department of ESA, ESTEC/SRE-C, Postbus 299, 2200 AG Noordwijk, The Netherlands N. Huelamo LAEX, Depto. Astrofísica, Centro de Astrobiología (INTA-CSIC), P.O. Box 78, E-28691 Villanueva de la Cañada, Spain A. Krivov Astrophysikalisches Institut und Universitätssternwarte, Friedrich-Schiller-Universität, Schillergäßchen 2-3, 07745 Jena, Germany J. Lebreton Université Joseph Fourier - Grenoble 1/CNRS, Laboratoire d’Astrophysique de Grenoble (LAOG) UMR 5571, BP 53, 38041 Grenoble Cedex 09, France R. Liseau Department of Radio and Space Science, Chalmers University of Technology, Onsala Space Observatory, 439 92 Onsala, Sweden C. Martin-Zaidi Université Joseph Fourier - Grenoble 1/CNRS, Laboratoire d’Astrophysique de Grenoble (LAOG) UMR 5571, BP 53, 38041 Grenoble Cedex 09, France A. Mora ESA-ESAC Gaia SOC, P.O. Box 78, 28691 Villanueva de la Cañada, Madrid, Spain M. Morales-Calderon Spitzer Science Center, California Institute of Technology, 1200 E California Blvd, 91125 Pasadena, USA. H. Nomura Department of Astronomy, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan E. Pantin CEA/IRFU/SAP, AIM UMR 7158, 91191 Gif-sur-Yvette, France I. Pascucci Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA N. Phillips SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory Edinburgh, UK Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK L. Podio Kapteyn Astronomical Institute, P.O. Box 800, 9700 AV Groningen, The Netherlands D.R. Poelman School of Physics & Astronomy, University of St. Andrews, North Haugh, St. Andrews KY16 9SS, UK S. Ramsay European Southern Observatory, Karl-Schwarzschild-Strasse, 2, 85748 Garching bei München, Germany. B. Riaz Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA K. Rice SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory Edinburgh, UK Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK E. Solano Spanish Virtual Observatory, Unidad de Archivo de Datos, Depto. Astrofísica, Centro de Astrobiología (INTA-CSIC), P.O. Box 78, E-28691 Villanueva de la Cañada, Spain H. Walker The Rutherford Appleton Laboratory, Chilton, Didcot, OX11 OQL, UK G.J. White The Rutherford Appleton Laboratory, Chilton, Didcot, OX11 OQL, UK Department of Physics & Astronomy, The Open University, Milton Keynes MK7 6AA, UK and The Rutherford Appleton Laboratory, Chilton, Didcot, OX11 OQL, UK J.P. Williams Institute for Astronomy (IfA), University of Hawaii, 2680 Woodlawn Dr., Honolulu, HI 96822, USA G. Wright UK Astronomy Technology Centre, Royal Observatory, Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, UK

Received 2010; accepted

Key Words.:

Stars: pre-main sequence; (Stars:) planetary systems; (Stars:) circumstellar matter; protoplanetary discs; Infrared: planetary systems

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 $μ$ m, [CII] 157.7 $μ$ m, CO 72.8 and 90.2 $μ$ m, and o-H $_{2}$ O 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 / 13}$ CO 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 $\pm 2.0 \times$ 10 $^{- 3}$ M $_{⊙}$ is still present in this disc, in agreement with earlier CO observations.

1 Introduction

Giant gas-planets form in protoplanetary discs within the first 10 Myr after protostar formation. Therefore the amount of gas that is present in a disc at a given time is very important, as it determines whether these planets can still be formed. The study of the dust is equally important, as it witnesses the first steps of planet formation. A lot is known about the dust in protoplanetary discs thanks to the Space Observatories ISO and Spitzer : different degrees of dust processing were observed with no clear relation to stellar properties, while grain growth could be related to the disc structure (e.g. Sicilia-Aguilar et al. sicilia2007 (), Meeus et al. meeus2009 ()). Gas is more difficult to observe, as the spectral lines are not very strong. Most gas studies are based on CO lines in the near-IR (warm gas) or in the mm, where CO freeze-out onto grains is a complicating factor (e.g. Brittain et al. brittain2007 (), Dent et al. dent2005 ()). However, as gas lines are stronger in the far-IR, we expect this field to experience major breakthroughs in the coming years based on data obtained with the Herschel Space Observatory (Pilbratt et al. pilbratt2010 ()), which provides sensitive far-IR photometry and spectroscopy.

The source HD 169142 is a young - age 6 $_{- 3}^{+ 6}$ Myr (Grady et al. grady2007 ()) - Herbig Ae star with an IR to millimetre excess attributed to a circumstellar (CS) disc. Submillimeter Array (SMA) observations show a disc in Keplerian rotation with radius $r$ = 235 AU and inclination $i \sim$ 13°(Raman et al. raman2006 ()) and a total gas mass of 0.6-3.0 $\times 10^{- 2}$ M $_{⊙}$ (Panić et al. panic2008 ()). For a Herbig Ae, the star is unusual as it has a small near-IR excess (e.g. Dominik et al. dominik2003 ()). Furthermore, SWS/ISO spectra revealed that the silicate 10 $μ$ m feature, detected in the majority of the Herbig Ae/Be stars (only 8 out of 53 lack the feature, Juhász et al., juhasz2010 ()) is lacking in HD 169142 (Meeus et al. meeus2001 ()). This absence can be explained if the silicate grains are either too large or too cold to emit at 10 $μ$ m (Meeus et al. meeus2002 ()). On the other hand, features of polycyclic aromatic hydrocarbons (PAHs), which can be excited by UV photons, were clearly detected with ISO and Spitzer. Based on the near to far-IR ratio, Grady et al. (grady2007 ()) suggested that the inner region has already cleared some material, and that the inner and outer disc are not coupled. However, this might also be explained by a reduced opacity due to grain growth.

In this paper we show the first Herschel observations of HD 169142, and use radiative transfer and chemical models to constrain the gas properties of the outer disc.

2 Herschel PACS observations

Our observations are part of the Herschel Open Time Key Project GASPS (P.I. Dent, see Mathews and Thi, this issue, and Dent et al. dent2010 ()). We obtained PACS (Poglitsch et al. poglitsch2010 ()) photometry (obsid 1342183656, Point Source, 400 s) and spectroscopy (obsid 1342186309, PacsRangeSpec, 5150 s and obsid 1342186310, PacsLineSpec, 1669 s). The spectroscopic data were reduced with the developer build version 3.0.1212 of the Herschel Interactive Processing Environment (HIPE; Ott ott2010 ()), using standard tasks provided in HIPE. In order to conserve the best signal and not to introduce additional noise, we only extracted the central spaxel and corrected for the flux loss with an aperture correction. Furthermore, we applied another correction factor as supplied by the PACS team (1.3 in the blue, and 1.1 in the red), to obtain a flux calibration accuracy of 40%. More details on the data reduction can be found in Mathews et al. (mathews2010 ()). In Fig. 1 we show our only detection, the fine structure line [OI] 63.2 $μ$ m. We also searched for other gas lines that were targeted in our spectroscopy, but did not detect: [OI] 145.5 $μ$ m, [CII] 157.7 $μ$ m, CO 72.8 and 90.2 $μ$ m, nor o-H $_{2}$ O 78.7 and 179.5 $μ$ m. In Table 1 we list the measured line properties. We note that in the LWS/ISO spectrum of HD 169142, the [CII] 157.7 $μ$ m line was detected, but LWS uses a larger beam ( $\sim$ 40″) than the size of the central spaxel of PACS (9 $\aas@@fstack′′$ 4 by 9 $\aas@@fstack′′$ 4). Therefore we also analysed the other 24 spaxels of the IFU and found a tentative detection in spaxels 13 and 14, which are to the West of HD 169142. The cause of this emission (if confirmed) is currently unknown.

Figure 1: Observed [OI] 63 $μ$ m emission line (black) together with our line fit (red; parameters listed in Table 1).

Species	Transition	$λ_{c}$	Continuum	Line Flux
		( $μ$ m)	(Jy)	(10 $^{- 18}$ W/m $^{2}$ )
[OI]	$^{3} P_{1} \to^{3} P_{2}$	63.18	18.89 (0.13)	71.7 (3.8)
[OI]	$^{3} P_{0} \to^{3} P_{1}$	145.53	13.76 (0.03)	$< 10.4$
[CII]	$^{2} P_{3 / 2} \to^{2} P_{1 / 2}$	157.74	14.48 (0.03)	$< 6.4$
o-H $_{2}$ O	2 $_{12} \to 1_{01}$	179.53	13.12 (0.05)	$< 8.8$
o-H $_{2}$ O	4 $_{23} \to 3_{12}$	78.74	17.82 (0.32)	$< 10.6$
CO	36 $\to$ 35	72.84	17.92 (0.04)	$< 15.5$
CO	29 $\to$ 28	90.16	17.00 (0.05)	$< 10.6$
$^{12}$ CO	2 $\to$ 1	1300.40	0.169 (0.005)	0.093 (0.004)
$^{13}$ CO	2 $\to$ 1	1360.22	0.169 (0.005)	0.048 (0.004)

Table 1: Summary of the PACS targeted lines and the CO lines that were observed with the SMA (Raman et al. raman2006 (), Panić et al. panic2008 ()). We list the central wavelength, continuum and line flux or upper limits (three sigma) in the case of a non-detection. Between brackets we give the one sigma statistical error for the detections.

Photometric observations were obtained in Point Source Mode for both the blue (70 $μ$ m) and red (160 $μ$ m) bands, and reduced using HIPE v2.3.4 with the Point Source Observations pipeline. The source is spatially unresolved, with FWHM 5″ at 70 $μ$ m and 11″ at 160 $μ$ m. We derived photometry using an aperture of 21″, and applied an aperture correction from the PACS PhotChopNod Release Note (Feb. 22, 2010). This gives a flux of 27.35 $\pm$ 0.03 Jy at 70 $μ$ m and 17.39 $\pm$ 0.05 Jy at 160 $μ$ m, with a flux calibration uncertainty of 5% in the blue, and 10% in the red. Applying the flux uncertainty of 40% for the spectroscopy, our photometry is consistent with the continuum fluxes measured from the spectral scans. Given the smaller error in the photometry, however, we give preference to those fluxes, and expect to reduce the spectroscopic flux uncertainties in the future.

3 Determination of stellar parameters

Because previous studies list conflicting effective temperatures, we re-estimated the fundamental parameters. Photometry excluding the contribution from the cool companion (at a distance of 9 $\aas@@fstack′′$ 3; Grady et al. grady2007 ()) was collected from Cutri et al. (cutri2003 ()), Zacharias et al. (zacharias2004 ()) and Sylvester et al. (sylvester1996 ()). We also analysed optical spectra, taken with CES/CAT with a resolution of $\sim$ 60000. Temperatures were selected by simultaneously fitting the lines from neutral and singly ionised species. The synthetic stellar spectra were computed using the ATLAS9 and SYNTHE codes by Kurucz (kurucz1993 ()). $χ^{2}$ tests show that models between $T_{e f f} = 7500$ K, $log g_{*} = 4.0$ , [Fe/H]= $- 0.50$ and $T_{e f f} = 7800$ K, $log g_{*} = 4.1$ , [Fe/H]= $- 0.25$ match our data.

Furthermore, the International Ultraviolet Explorer (IUE) obtained five spectra of HD 169142: one between 1200–1900 Å (SW) and four between 1900–3200 Å (LW). The SW spectrum is unusable for measurements below 1650 Å, but suggests emission lines of OI (1304 Å), CII (1335 Å) and CIV (1550 Å) - unfortunately, most are dominated by bad pixels. The LW spectra do not show any variability (over a period of 4.5 months), match the stellar photosphere, and do not show emission features. Overall, the IUE data (obtained from the INES INTA archive) do not give evidence for a UV excess.

4 Analysis

4.1 Dust and continuum modelling

Far-IR lines emerging from a CS disc are affected by the stellar UV irradiation, disc mass and geometry, dust size and composition, as well as PAH abundance. In order to interpret the lines observed with Herschel, it is crucial to first obtain a solid knowledge of the dust structure in the disc, based on as many observations as possible, like SED, scattered light images and visibilities. Each of these observations provide complementary views of the disc structure and the dust properties. The disc model is calculated with the Monte Carlo radiative transfer code MCFOST (Pinte et al. pinte2006 (); pinte2009 ()), as outlined below.

We consider an axisymmetric, slightly flared density structure with a Gaussian vertical profile, assuming power-laws for the surface density and scale height. We assume homogeneous and spherical dust grains (Mie theory), with sizes distributed according to the power-law $f (a) \propto a^{- 3.5}$ between $a_{m i n}$ and $a_{m a x}$ , where $a$ is the grain radius. The dust is assumed to be well-mixed with the gas, i.e. the dust/gas ratio is constant throughout the disc. The star is reproduced by a uniformly radiating sphere with previously determined parameters $T_{e f f} = 7800$ K, $log g_{*} = 4.1$ , [Fe/H] = -0.25, and R $_{*}$ = 1.6 R $_{⊙}$ . Parameters are adjusted to simultaneously fit the SED, the Spitzer/IRS spectrum (Sloan et al. sloan2005 ()), the $1.1 μ$ m HST image (Grady et al. grady2007 ()) and the 1.3 mm SMA visibilities (Panić et al. panic2008 ()). The scattered light image mainly constrains the flaring index to a low value (around 1.0) and the mm visibilities indicate a surface density varying as $r^{- 1}$ .

We reproduce the 1.3 mm emission with a dust mass of $1.5 \times 10^{- 4}$ M $_{⊙}$ with a grain size distribution¹¹1Note that for these size distribution parameters, only 31.5 % of the dust mass is in grains $\leq 1$ mm, the rest (68.5%) is in larger grains. between $a_{m i n} = 0.03 μ$ m and $a_{m a x} = 1$ cm. The weak emission around 10 $μ$ m implies that there is a discontinuity in the disc $τ = 1$ surface, possibly due to a shadowed area or a gap. In this paper, we explore a solution with a gap. For the dust composition, we use a mixture of 70 % silicates (Draine et al. draine2003 ()) and 30 % amorphous carbon (ACAR sample; Zubko et al. zubko1996 ()), and calculate the effective optical index with the Bruggeman mixing rule. The strong PAH features around 10 $μ$ m can be reproduced with a low PAH abundance due to the low continuum emission in that region. We modelled the emission with a single grain-size (54 carbon atoms, positively ionised) and assume that the PAHs are uniformly distributed in the outer disc. We did not try to reproduce in detail the various bands observed with IRS, but rather constrain their abundance. We obtain M $_{P A H}$ /M $_{d u s t}$ = $5 \times 10^{- 4}$ , corresponding to $f_{P A H} = 0.03$ (where $f_{P A H}$ is the uniform PAH abundance relative to a standard ISM abundance of $10^{- 6.52}$ PAH particles/H-nucleus, $m_{P A H} = 667$ amu, and gas/dust = 100).

Figure 2: Best fit with MCFOST to the observed data. Shown is the photometry obtained from the literature (black), Spitzer IRS spectrum (red), PACS photometric observations (blue crosses), and PACS continua derived from the spectroscopic observations (red crosses). The green line represents the best MCFOST model, in red the atmosphere model. The blue dot-dashed and pink dashed lines are the contributions (scattered light and thermal emission) from the inner and outer disc, respectively.

Parameter	Inner Disc	Outer Disc
$r_{i n}$ (AU)	0.1	20
$r_{o u t}$ (AU)	5	235
surface dens. exp $ϵ$	-1.0	-1.0
flaring exponent $β$	1.05	1.00
ref. scale height $h_{0}$	0.07 AU @ 1 AU	12.5 AU @ 100 AU
M $_{d u s t}$ (M $_{⊙}$ )	$2 \times 10^{- 9}$	$1.5 \times 10^{- 4}$
$f_{P A H}$	0	0.03

Table 2: Parameters of the best fitting dust model.

In Fig. 2 we show the best fitting model on top of the SED. Our model provides good constraints on the outer disc (mid-, far-IR and mm emission): dust mass, scale height, flaring index and surface density profile, as well as dust properties and amount of PAHs. However, the inner disc (near-IR emission) remains poorly constrained. We tried to improve the fit by moving the inner radius closer to the star ( $r_{i n} = 0.06$ AU), but then the temperature is too high (2400 K) for dust, even carbon, to survive. Because we lack simultaneous visible and near-IR photometry as well as spatially resolved near-IR data, we did not further improve the fit to the inner disc. However, as most of the mid-IR and mm lines originate from the outer disc, this is not critical for the present paper.

4.2 Gas and line modelling

	model #1	model #2	model #3	observed
gas/dust	1	100	33	–
$f_{P A H}$	0.01	0.0055	0.0087	–
$f_{U V}$	0.005	0.0	0.0	–
Line	Line Fluxes $[10^{- 18} W / m^{2}]$
[OI] 63.2 $μ$ m	154	71.6	71.6	71.7
[OI] 145.5 $μ$ m	5.17	10.1	7.01	$<$ 10.4
[CII] 157.7 $μ$ m	4.58	0.04	0.06	$<$ 6.4
o-H $_{2}$ O 179.5 $μ$ m	5.66	5.15	1.76	$<$ 8.8
$^{12}$ CO 2 $\to$ 1	0.060	0.092	0.093	0.093
$^{13}$ CO 2 $\to$ 1	0.011	0.059	0.048	0.048
$^{12}$ CO $/^{13}$ CO	5.69	1.55	1.92	1.94

Table 3: Predicted line fluxes for different models, compared with the observed fluxes. The best fit is obtained by model #3. All models in this table have turbulent line broadening

v_{t u r b}

= 0.15 km/s.

Figure 3: Spatial origin of the various gas emission-lines. The two vertical white dashed lines indicate 15% and 85% of the radially cumulative face-on line flux, respectively, i.e. 70% of the line flux originate from within the two white dashed lines. The red lines mark the cells that contribute most to the line flux in their vertical column.

The best model reproducing the continuum observations is now fed into the gas thermo-chemical code ProDiMo (Woitke et al. woitke2009 ()) to calculate the chemical and gas temperature structure in the disc and to predict the line fluxes, following the pipeline described in Woitke et al. (woitke2010 (); see their Fig. 1). This final modelling step has additional free parameters such as the dust/gas ratio. We computed models for the following parameters: gas/dust mass ratio $\in [1, 100]$ , turbulent broadening with $v_{t u r b} \in [0, 0.15]$ km/s, PAH abundance $f_{P A H} \in [0, 0.06]$ and stellar UV excess $f_{U V} = L_{U V} / L_{⋆} \in [0, 0.005]$ (for details, see Woitke et al. woitke2010 ()). Due to the stellar UV irradiation, the models generally result in disc surface-layers where the gas is warmer than the dust. Table 3 shows the calculated line fluxes for a few selected models. An important result from our modelling effort is that the gas heating by PAHs plays a central role for the line flux predictions, in our case useful for the temperature sensitive $[O I] 63.2 μ$ m line. For each selection of the parameters dust/gas and $f_{U V}$ , we tuned the PAH abundance until a fit with the observed $[O I] 63.2 μ$ m line flux was obtained, if possible. However, models with UV excess $f_{U V} = 0.005$ result in a much too hot gas and hence too strong gas emission lines, even for $f_{P A H} = 0$ , and even if we decrease gas/dust $\to 1$ , as in model #1. For the standard ratio $g a s / d u s t = 100$ , we found model #2, which fits the $[O I] 63.2 μ$ m and $^{12}$ CO line fluxes, but the $^{13}$ CO line is slightly too high (2.5 sigma). A better fit is obtained with $g a s / d u s t = 33$ (model #3), which simultaneously fits all three detected lines, and agrees with all other line upper limits. However, we emphasise that apart from the uncertainties in the observations, systematic uncertainties in the physical description of various chemical and radiative processes in the models render a proper gas mass determination difficult. The observational uncertainty in the $^{13}$ CO line flux - which appears to be a dominant gas mass tracer - translates to a ratio range $g a s / d u s t = 22$ to 50. An increase of the [OI] 63.2 $μ$ m flux by 40% (to account for the calibration uncertainty and match the photometry) would indicate a higher gas temperature. In our modelling, this is best interpreted by a larger PAH abundance ( $f_{P A H}$ = 0.02), as the other more relevant parameter, $f_{U V}$ , is constrained by the non-detection of the [CII] 157.7 $μ$ m line.

The spatial distribution of the $^{12 / 13}$ CO molecules predicted by our chemical model is consistent with the interferometric observations by Panić et al. (panic2008 ()). Figure 3 shows a detailed analysis of the spatial origin of the various emission lines, as derived from model #3. The [OI] 63.2 $μ$ m and 145.5 $μ$ m lines originate from a radial disc region extending from the beginning of the outer disc at 20 AU to about 75 AU and 60 AU, respectively. Both lines are optically thick ( $τ_{l i n e} \approx 10 - 100$ ) and come from relative heights $z / r \approx 0.2 - 0.3$ . The [CII] 157.7 $μ$ m line is optically thin and extraordinarily weak in this model, due to the lack of stellar UV photons $< 110$ nm that are capable of ionising carbon ( $f_{U V} = 0$ ). The $^{12}$ CO and $^{13}$ CO $J = 2 \to 1$ lines probe the conditions in the outermost disc regions $43 - 185$ AU and $52 - 190$ AU, respectively, and are extremely optically thick ( $τ_{l i n e} \approx 100 - 5000$ ). Whereas the $^{12}$ CO line comes from the PDR-like CO surface $z / r \approx 0.4$ , the more transparent $^{13}$ CO line probes slightly deeper layers $z / r \approx 0.2$ . These findings are quite stable and robust among all calculated models.

5 Conclusion

We presented the first PACS observations of HD 169142 and showed the unique capability of Herschel to obtain an independent gas mass determination by simultaneously modelling the atomic fine-structure lines that are temperature sensitive (due to their high excitation energies), when combined with ground-based observations of $^{12 / 13}$ CO-lines. We showed that the observations are consistent with a simple model of a disc hosting a gap, with a continuous structure between 20 and 200 AU in which the gas and dust are well-mixed, but stress that this solution is not unique. We determined the location of the emitting species, and constrained the gas/dust ratio to be $\sim$ 22-50. We also derived that the UV excess emission, if present, must be lower than $f_{U V}$ = 0.005. We derive a gas mass between 3.0 and 7.0 $\times$ 10 $^{- 3}$ M $_{⊙}$ , in agreement with the gas mass determination based on the CO lines alone (Panić et al. panic2008 ()). Despite indications that the HD 169142 disc is transitional, it still is a gas-rich disc. Furthermore, although there is dusty material close to the central star, the UV excess induced by accretion is extremely weak, if at all present. We will apply a similar modelling strategy for future GASPS observations, with the aim to study gas dissipation in a wide range of protoplanetary discs.

Acknowledgements.

We thank the PACS instrument team for their dedicated support, O. Panić for a discussion on her data and A. Juhász for providing the Spitzer/IRS data. C. Eiroa, G. Meeus, I. Mendigutía and B. Montesinos are partly supported by AYA 2008-01727; J.M. Alacid and E. Solano by AYA 2008-02156 and I. d. Gregorio by AYA 2008-06189 (Spanish grants). C. Pinte acknowledges funding from the EC 7

^{t h}

FP as a Marie Curie Intra-European Fellow (PIEF-GA-2008-220891). J.-C. Augereau, G. Duchêne, J. Lebreton, C. Martin-Zaidi, F. Ménard and C. Pinte acknowledge PNPS, CNES and ANR (contract ANR-07-BLAN-0221) for financial support. D.R. Ardila, S.D. Brittain, B. Danchi, C.A. Grady, I. Pascucci, B. Riaz, A. Roberge, G. Sandell and C.D. Howards acknowledge NASA/JPL for funding support. PACS has been developed by a consortium of institutes led by MPE (Germany) and including UVIE (Austria); KUL, CSL, IMEC (Belgium); CEA, OAMP (France); MPIA (Germany); IFSI, OAP/AOT, OAA/CAISMI, LENS, SISSA (Italy); IAC (Spain). This development has been supported by the funding agencies BMVIT (Austria), ESA-PRODEX (Belgium), CEA/CNES (France), DLR (Germany), ASI (Italy), and CICT/MCT (Spain).

References

(2007) Brittain, S.D., Simon, Th. , Najita, J.R., and Rettig, T.W. 2007, ApJ, 659, 685
(2003) Cutri, R.M., Skrutskie, M.F., van Dyk, S. et al. 2003, The IRSA 2MASS All-Sky Point Source Catalog, NASA/IPAC Infrared Science Archive
(2005) Dent, W.R.F., Greaves, J.S., Coulson, I.M. 2005, MNRAS 359, 663
(2010) Dent, W.R.F. et al. 2010, PASP
(2003) Dominik, C., Dullemond, C.P., Waters, L.B.F.M. et al. 2003, A&A 398, 607
(2003) Draine, B. T. 2003, ApJ, 598, 1017
(2007) Grady, C. A., Schneider, G., Hamaguchi, K., et al. 2007, ApJ, 665, 1391
(submitted) Juhász, A., Bouwman, J. Henning, Th. et al. ApJ, submitted
(1993) Kurucz, R.L. 1993, CD-ROM No. 13. Cambridge, MA, SAO
(2010) Mathews, G.S. et al., this issue
(2001) Meeus, G., Waters, L.B.F.M., Bouwman, J. et al. 2001, A&A 365, 476
(2002) Meeus, G., Bouwman, J., Dominik, C., et al. 2002, A&A 392, 1039
(2009) Meeus, G., Juhász, A., Henning, Th. et al. 2009, A&A, 497, 379
(2010) Ott, S. 2010, in ASP Conference Series, Astronomical Data Analysis Software and Systems XIX, Y. Mizumoto, K.-I. Morita, and M. Ohishi, eds., in press
(2008) Panić, O., Hogerheijde, M.R., Wilner, D., and Qi, C. 2008, A&A 491, 219
(2010) Pilbratt, G. et al. 2010, this issue
(2006) Pinte, C., Ménard, F., Duchêne, G., & Bastien, P. 2006, A&A, 459, 797
(2009) Pinte, C., Harries, T. J., Min, M., et al. 2009, A&A, 498, 967
(2010) Poglitsch et al. 2010, this issue
(2006) Raman, A., Lisanti, M., Wilner, et al. 2006, AJ, 131, 2290
(2007) Sicilia-Aguilar, A., Hartmann, L. W., Watson, D. et al., ApJ 659, 1637
(2005) Sloan G.C., Keller, L.D., Forrest, W.J. et al. 2005, ApJ632, 956
(1996) Sylvester, R.J., Skinner, C.J., Barlow, M.J. et al. 1996, MNRAS, 279, 915
(2010) Woitke, P., Pinte, C., Tilling, I., et al. 2010, MNRAS, in press
(2009) Woitke, P., Kamp, I. Thi, W.-F. et al. 2009, A&A 501, 383
(2004) Zacharias, N., Monet, D.G., Levine, S.E. et al. 2004, BAAS Vol. 36, p.1418
(1996) Zubko, V.G., Mennella V., Colangeli L. & Bussoletti E. 1996, MNRAS 282, 1321

Comments 0

Request Comment

You are adding the first comment!

How to quickly get a good reply:

Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.

The feedback must be of minimum 40 characters and the title a minimum of 5 characters

Cancel

Loading ...

Gas in the protoplanetary disc of HD 169142: Herschel’s view ††thanks: Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.