EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 519 (September 2010) A&A, 519 (2010) A110 Full HTML
Free Access
Issue
A&A
Volume 519, September 2010
Article Number A110
Number of page(s) 9
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/200913709
Published online 20 September 2010
A&A 519, A110 (2010)

Observations of warm molecular gas and kinematics in the disc around HD 100546

O. Panic1,2 - E. F. van Dishoeck1,3 - M. R. Hogerheijde1 - A. Belloche4 - R. Güsten4 - W. Boland1,5 - A. Baryshev6

1 - Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
2 - European Southern Observatory, Karl Schwarzschild str. 2, 85748 Garching, Germany
3 - Max-Planck-Institut für extraterrestrische Physik, Postfach 1312, 85741 Garching, Germany
4 - Max-Planck Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
5 - Nederlandse Onderzoeksschool Voor Astronomie (NOVA), PO Box 9513, 2300 RA Leiden, The Netherlands
6 - SRON Netherlands Institute for Space Research, PO Box 800, 9700 AV Groningen, The Netherlands

Received 15 January 2010 / Accepted 3 May 2010

Abstract
Context. The disc around the Herbig Ae/Be star HD 100546 is one of the most extensively studied discs in the southern sky. Although there is a wealth of information about its dust content and composition, not much is known about its gas and large-scale kinematics. Many recent results have stressed the importance of studying both the gas and dust in discs. 12CO is an excellent gas tracer in the submillimetre, and the intensity ratio between lines originating from low and high rotational levels probes the gas temperature. Emerging submillimetre facilities in the Southern hemisphere allow us to characterise the gas and dust content in objects like HD 100546 better.
Aims. Our aim is to establish whether the disc is gas-rich and to study the disc temperature and kinematics.
Methods. We detected and studied the molecular gas in the disc at spatial resolution from 7 $\hbox{$.\!\!^{\prime\prime}$ }$7 to 18 $\hbox{$.\!\!^{\prime\prime}$ }$9 using the Atacama Pathfinder Experiment telescope[*]. We observed the lines 12CO J= 7-6, J= 6-5, J= 3-2, 13CO J= 3-2 and [C I] 3P2-3P1, as diagnostics of disc temperature, size, chemistry, and kinematics. We use parametric disc models that reproduce the low-J 12CO emission from Herbig Ae stars and we vary the basic disc parameters - temperature, mass, and size. With the help of a molecular excitation and radiative transfer code we fit the observed spectral line profiles.
Results. Our observations are consistent with more than 10 $^{-3}~M_{\odot}$ of molecular gas in a disc of $\approx$400 AU radius in Keplerian rotation around a 2.5 $M_{\odot}$ star, seen at an inclination of 50 $\hbox{$^\circ$ }$. The detected 12CO lines are dominated by gas at 30-70 K. Not detecting the [C I] line indicates excess ultraviolet emission above that of a B9 type model stellar atmosphere. Asymmetry in the 12CO line emission suggests that one side of the outer disc is colder by 10-20 K than the other. A plausible scenario is a warped geometry in the inner disc casting a partial shaddow on the outer disc. We exclude pointing offsets, foreground cloud absorption, and asymmetry in the disc extent as possible causes of the observed line asymmetry.
Conclusions. Efficient heating of the outer disc by the star HD 100546 ensures that low- and high-J 12CO lines are dominated by the outermost disc regions, indicating a 400 AU radius. The 12CO J= 6-5 line arises from a disc layer higher than disc midplane, and warmer by 15-20 K than the layer emitting the J= 3-2 line. The existing models of discs around Herbig Ae stars, asuming a B9.5 type model stellar atmosphere, overproduce the [CI] 3P2-3P1 line intensity from HD 100546 by an order of magnitude.

Key words: planetary systems: protoplanetary disks - stars: pre-main sequence - ISM: molecules

1 Introduction

Over the past decade our understanding of the structural and physical properties of discs around young stars has increased from basic theoretical modelling of the spectral energy distributions (SEDs) constrained by observations with no spatial information to modelling based not only on the SEDs, but also on spatially resolved dust observations, such as scattered light images and interferometry (Pinte et al. 2008; Tannirkulam et al. 2008; Panic et al. 2008). Two decades ago, the first submillimetre interferometer observations resolved the molecular gas emission spatially, and this allowed major progress in understanding the disc kinematics, structure, and chemistry (e.g., Koerner et al. 1993; Beckwith & Sargent 1987; Dutrey et al. 1994). A more recent example is a study of one of the brightest discs around a Herbig Ae star, HD 163296, shown in Isella et al. (2007). (Sub)millimetre gas and dust emission is the ideal probe of the global disc properties, such as size, mass, and radial distribution of disc material, because the bulk of the disc mass is located beyond 100 AU from the star, at temperatures of 10-50 K that dominate this part of the spectrum. Disc models that include constraints on both dust and molecular gas observations have stressed the importance of analysing the gas and dust components simultaneously, in the context of a single disc model with a three-dimensional temperature and density structure (Panic et al. 2009; Raman et al. 2006; Panic et al. 2008; Wilner et al. 2003).

Until recently, observations of rotational transitions of molecules in the submillimetre regime have been focused primarily on the low-J emission from 12CO, up to the J= 3-2 line (Greaves et al. 2000; Thi et al. 2001,2004; Qi et al. 2004; Dent et al. 2005). In two of the brightest and most studied sources, TW Hya and LkCa 15, the observations of higher-J transitions of 12CO, up to J= 6-5 ($E_{\rm k}$ = 116 K), have been compared to the low-J lines, providing estimates of the gas temperature in the intermediate-height molecular layer (van Zadelhoff et al. 2001a), crucial ingredients for chemical modelling of discs. These single-dish line spectra were fitted using simplistic disc models, deriving a temperature of 20-40 K in the 12CO line emitting layers of LkCa 15, and more than 40 K in TW Hya. Qi et al. (2006) analysed submillimetric interferometer observations of TW Hya in the context of a disc structure based on an accretion disc model (Calvet et al. 2002). Based on 12CO J= 6-5, 3-2 and 2-1 observations, they show that X-ray heating of the gas is efficient in this source, in addition to the stellar radiation field. Such diagnostics of gas heating and ionisation improve our understanding of how the gas content evolves in discs.

The emerging (sub)millimetre facilities in the Southern hemisphere like the Atacama Pathfinder EXperiment (APEX) and the Australia Telescope Compact Array (ATCA) are opening a window towards the star-forming regions of the Southern sky and are well suited to studies of circumstellar disc emission. These instruments also pave the way for future observations with the Atacama Large Millimetre/Submillimeter Array (ALMA), which will drastically improve our knowledge of disc structure and evolution (e.g. Guilloteau & Dutrey 2008). We use APEX receivers APEX-2a and CHAMP+ to observe the 12CO J= 7-6, J= 6-5, J= 3-2, 13CO J= 3-2 and [C I] 3P2-3P1 line emission towards the disc around the young intermediate-mass star HD 100546. A wealth of observations of dust in this bright disc (Malfait et al. 1998; Augereau et al. 2001; Ardila et al. 2007; Waelkens et al. 1996; Bouwman et al. 2003; Grady et al. 2001; Acke & van den Ancker 2006a) has motivated us to probe its molecular gas content and kinematics. The chosen transitions are particularly sensitive to the gas in the warm upper layers and kinematics of the outer disc. Our millimetre line observations probe the outer radius and inclination. The existing observational constraints on these parameters in the disc around HD 100546 provide an excellent basis for the analysis of our data. Our observations also provide a bridge toward even higher-J far infrared 12CO transitions to be observed with the Herschel Space Observatory. HD 100546 is a young B9V type, 2.5 $M_{\odot}$ star, classified as a Herbig Be star because of its isolation, infrared excess, and silicate emission (Malfait et al. 1998; Thé et al. 1994). With a distance of 103 $\pm$ 6 pc, measured by Hipparcos, this is one of the nearest Herbig Ae/Be stars. The age of the star is estimated to be greater than 10 Myr (van den Ancker et al. 1998). This makes the presence of circumstellar material intriguing, considering that discs are found to dissipate within 10 Myr in most young stars (e.g., Hillenbrand 2008; Hernández et al. 2007; Hollenbach et al. 2000). Based on SED modelling, Bouwman et al. (2003) postulate there is an inner hole in the disc with a 10 AU radius, possibly caused by a Jupiter-sized planet (see also Acke & van den Ancker 2006a). Direct evidence of cold disc material at larger radii is provided by ATCA observations of Wilner et al. (2003) at 89 GHz (3.4 mm) and 2 $\hbox{$^{\prime\prime}$ }$ resolution, with the flux of 36 $\pm$ 3 mJy, values consistent with the 1.3 mm observations of Henning et al. (1998). They do not detect HCO+ J= 1-0 line emission and speculate that photodissociation of its progenitor species, 12CO, in the upper disc layers or an overall gas depletion may be the reason for this. Recent spectroastrometric observations of rovibrational 12CO transitions by van der Plas et al. (2009) suggest that this molecular species is missing from the inner disc, at least up to 8 AU from the star (also see Brittain et al. 2009).

Scattered light imaging of HD 100546 reveals the disc extending up to 4 $\hbox{$^{\prime\prime}$ }$ from the star viewed at an inclination of 50 $\hbox{$^\circ$ }$ and an interesting disc structure resembling spiral arms (Augereau et al. 2001; Pantin et al. 2000; Grady et al. 2001). This structure was interpreted as the result of disc perturbation by a companion (Quillen et al. 2005) or a warped disc structure (Quillen 2006). Coronographic imaging by Augereau et al. (2001) shows steep surface brightness profiles in the environment of HD 100546 indicative of optically thin emission in the near-infrared, with surface densities as low as 10-3 g cm-2. Their images trace the emission of small dust (${<}5~\mu$m), extending out to 800 AU from the star. The authors suggest there are an optically thick disc with a 400 AU radius and an optically thin flattened halo or envelope farther from the star. The scattered light images hint at the presence of gas in the disc that supports the disc vertical structure.

In this work, we detect and study the molecular gas in the disc, in particular its kinematics and temperature, at spatial resolution from 7 $\hbox{$.\!\!^{\prime\prime}$ }$7 to 18 $\hbox{$.\!\!^{\prime\prime}$ }$9. In Sect. 2 we present our observations of 12CO, 13CO, and [C I] lines. All 12CO lines are detected, there is a tentative detection of the 13CO line, while [C I] emission is not detected. We model the spectra in Sect. 3 discussing the implications for the disc size, mass, and kinematics. We identify the regions dominating the observed lines, and derive their temperatures. Section 4 summarises our results.

2 Observations and results

The observations of 12CO J= 6-5 at 691.472 GHz and [C I] 3P2-3P1 ([C I] J= 2-1 hereafter) at 809.344 GHz towards HD 100546 at RA = 11 $\rm ^{h}33^{m}$25 $\hbox{$.\!\!^{\rm s}$ }$4 and Dec = -70 $\hbox{$^\circ$ }$11 $\hbox{$^\prime$ }$41 $\hbox{$^{\prime\prime}$ }$ (J2000) were obtained simultaneously, using the CHAMP+ heterodyne array receiver (Güsten et al. 2006) on APEX on 2008 November 11. The 7 pixels in each wavelength band are arranged in a hexagon of 6 pixels around one central pixel pointed towards the source, with beam sizes of 9 $\hbox{$^{\prime\prime}$ }$ at 691 GHz and 7 $\hbox{$.\!\!^{\prime\prime}$ }$7 at 809 GHz. The data were obtained in the position switching mode in which the telescope moved between the on-source position and a reference position at an equatorial offset of (-1000 $\hbox{$^{\prime\prime}$ }$, 0 $\hbox{$^{\prime\prime}$ }$) from the source. The backend consisting of the Fast Fourier Transform Spectrometer unit was used on all pixels, providing a spectral resolution of 0.37 MHz or 0.14-0.16 km s-1 at these frequencies, and covering a bandwidth of 1.5 GHz or 4096 channels. Main beam efficiencies are 0.40 at 691 GHz and 0.37 at 809 GHz. Calibration is uncertain by $\approx$30$\%$ at both frequencies. Pointing was performed directly prior to and after the on-source integration providing an accuracy better than 3 $\hbox{$^{\prime\prime}$ }$. The pointing source was 07454-7112, about 20 $\hbox{$^\circ$ }$ away from HD 100546, at the same airmass.

\begin{figure} \par\includegraphics[width=8.8cm,clip]{13709fg1.ps} \end{figure} Figure 1:

From top to bottom: spectra of the 12CO J= 7-6, J= 6-5, J= 3-2, 13CO J= 3-2 and CI J= 2-1 lines observed towards HD 100546 shown with black lines plotted as a histogram. The 12CO and 13CO spectra computed from our models of the line emitting layers of the disc (see Sects. 3.2.1 and 3.2.2) are shown with continuous red lines. At 100 AU from the star, our vertically isothermal models of the 12CO line emission have temperatures T100 = 60-70 K, while a lower temperature T100 = 30 K is used for the 13CO line.

Open with DEXTER

The CO J = 6-5 line was also observed on 2008 November 10 in jiggle mode and its intensity and spectral profile were found to be the same within 20%, as in the 2008 November 11 data. Since these data are noisier, the data sets taken on different dates were not combined. During this latter observation, the high band was tuned to 12CO J= 7-6 at 806.665 GHz, with the beam efficiency of 0.37. In the further analysis, we used only the high S/N 12CO J= 6-5 data that were taken on 2008 November 11. The 12CO J= 3-2 line at 345.796 GHz and the 13CO J= 3-2 line at 330.588 GHz were observed on 2005 July 27 and 28 with the APEX-2a receiver using a single pointing. The channel spacing of these data is 61 kHz or 0.05-0.06 km s-1, and the spectral resolution is 98 kHz or 0.09 km s-1 (Klein et al. 2006). The main beam efficiency of APEX at 346 GHz is 0.73, and the beam sizes are 18 $\hbox{$.\!\!^{\prime\prime}$ }$1 at 345.796 GHz and 18 $\hbox{$.\!\!^{\prime\prime}$ }$9 at 330.588 GHz, respectively. Calibration is uncertain by $\approx$20$\%$ at these frequencies. Our 12CO J= 3-2 line data were presented in Panic & Hogerheijde (2009). In all our observations, APEX forward efficiency of 0.95 was taken into account.

Table 1:   Integrated intensities of the observed spectral lines.

The data were reduced and analysed using the CLASS and STARLINK software. We detect all the observed 12CO line transitions at 7-20$\sigma$ in terms of the integrated intensity. The 13CO J= 3-2 line is marginally detected at 2$\sigma$, and the [C I] J= 2-1 line is not detected. Figure 1 shows the observed spectra, corrected for the beam efficiency and re-binned to a lower spectral resolution. A first-degree polynomial fit, except for the 0-10 km s-1 range in which the lines are emitted, is subtracted from the full spectral range for each data set. The 12CO J= 3-2 and J= 6-5 lines are detected at the highest signal-to-noise ratio and show a double-peaked profile characteristic of disc rotation. The intensities integrated over the velocity range 0-10 km s-1, over which line emission is detected, are listed in Table 1, together with the full width at half-maximum (FWHM) of the lines with sufficiently well-defined profiles.

In addition to the observations towards the source, the CHAMP+ array provides measurements at nearby offsets. This setup provides an excellent way to discern the emission from the disc, with an estimated size of at least 400 AU in radius from the surrounding material known to be present further away from the star (Augereau et al. 2001; Henning et al. 1998; Grady et al. 2001). The 400 AU lower limit to the size of the gas disc is set by the scattered light observations (Augereau et al. 2001) is close to the spatial resolution of our data. The central pixel of our CHAMP+ data probes the 12CO J= 6-5 line emission from the region of 9 $\hbox{$^{\prime\prime}$ }$ that is centered on the position of the star (460 AU radius), while the surrounding pixels probe the more distant regions (at a distance 1000-2000 AU from the star). Similarly, at the frequency of the 12CO J= 7-6 line, a small 7 $\hbox{$.\!\!^{\prime\prime}$ }$7-sized region (400 AU radius) around the star is probed with the central pixel, whereas the surrounding pixels probe the regions that are located further away. Table 2 provides an overview of the pixel positions and the corresponding fluxes integrated over a velocity range from 0 to 10 km s-1, over which the 12CO lines are firmly detected in the central pixel. Compared to the on-source fluxes, these measurements show that the 12CO emission from the surrounding material is about 20 times weaker than from the region within 400 AU from the star. The 12CO J= 3-2 line was observed with a single pointing, and the beam of 18 $\hbox{$.\!\!^{\prime\prime}$ }$1 is large enough to include any emission from regions beyond 400 AU. However, the strong resemblance in the line profile suggests that both low-J and high-J lines arise from the disc and that any contribution to the line emission by an extended low-temperature and low-velocity component is negligible.

Table 2:   The 12CO J= 6-5 and J= 7-6 integrated line intensities at offset positions.

Our data testify to the presence of warm molecular gas at distances less than roughly 400-500 AU from the star, corresponding to the smallest beam sizes in our observations (7 $\hbox{$.\!\!^{\prime\prime}$ }$7 and 9 $\hbox{$^{\prime\prime}$ }$). We obtain the following integrated intensity ratios, corrected for beam dilution (scaled to the same beam): 12CO lines (6-5)/(3-2) = 1.1 $\pm$ 0.6 and (12CO 3-2)/(13CO 3-2) $\geq$ 2.8 $\pm$ 2.2. The errors include both the rms and flux calibration uncertainty.

There is a clear asymmetry in the profile of the 12CO J= 3-2 and J= 6-5 lines, observed at a high signal-to-noise ratio and characterised by a 20-40% difference in the line intensity between the two peaks. This asymmetry is certain, given that the spectra shown in Fig. 1 have been smoothed to a lower spectral resolution and that each spectral bin contains an average of several measurements.

In some sources, asymmetries like this are explained through confusion with foreground cloud emission, as seen in 12CO lines towards IM Lup (van Kempen et al. 2007) and DL Tau (Simon et al. 2000). However, in the off-source pixels of the CHAMP+ array, we see no significant extended emission. While a slight pointing offset in the 6-5 line observations towards a symmetric source may result in a line asymmetry, the persistance of line asymmetry in our 18 $\hbox{$.\!\!^{\prime\prime}$ }$1 resolution observations of the 3-2 line indicates that the asymmetry is caused by the source properties. Furthermore, the two observations were taken years apart, so a systematic offset in both observations is highly unlikely. We therefore exclude pointing offset as a cause of the observed line asymmetry. If one side of the disc extends slightly further out than the other (e.g., on the SE side), the increase in disc surface on that side will contribute to the line flux on one side of the spectral line, and result in line asymmetry. This may be a plausible explanation for the 3-2 line, but not the observed 6-5 line with a 9 $\hbox{$^{\prime\prime}$ }$ beam size, unless the asymmetry is within 450 AU. Such an asymmetry may affect the scattered light images of HD 100546, with the disc extending further to the SE than to the NW of the star. However, this is not seen in the observations of Augereau et al. (2001). Furthermore, the additional emission from this region would dominate the frequencies closer to the line centre and not the blue-shifted peak.

Asymmetry is seen in 12CO lines from sources with a pronounced structural disc asymmetry, like HD 141569 (Dent et al. 2005). Another possibility is a temperature asymmetry in the disc. In the following section, we analyse and model the observed lines and investigate these scenarios.

3 Discussion

3.1 Line ratios as temperature and density diagnostics

3.1.1 12CO 3-2/13CO 3-2 line ratio and disc density

The submillimetre 12CO line emission from circumstellar discs is generally optically thick and arises from warm upper disc layers. These lines are particularly sensitive to the temperature of these layers, and therefore to the stellar and external illumination of the disc surface. Comparison to isotopologue line emission is often used to assess the optical thickness. If the 13CO 3-2 emission is treated as a detection, we obtain the line ratio of (12CO 3-2)/(13CO 3-2) = 2.8 $\pm$ 2.2. We make a rough estimate of the column density needed to reproduce the derived ratio using the RADEX online tool[*] (van der Tak et al. 2007). For kinetic temperatures 20-40 K, H2 densities 106-108 cm-3, and an isotopic ratio [12C]/[13C] = 77 (Wilson & Rood 1994), the (12CO 3-2)/(13CO 3-2) ratio derived above indicates the 12CO column density of roughly 1018 cm-2, and excludes lower column densities. Using a standard 12CO abundance of 10-4 with respect to H2, this translates to the gas column density of more than 1022 cm-2.

3.1.2 12CO 6-5/12CO 3-2 line ratio and disc temperature

We compared the line ratio (12CO 6-5)/(12CO 3-2) = 1.1 $\pm$ 0.6 to the calculations done using RADEX, and assumed that these two lines are dominated by the same region in the disc. This assumption is generally true up to 200-300 AU from the star (see calculations done in van Zadelhoff et al. 2001b). For H2 densities of 106-108 cm-3, and 12CO column densities of 1017-1019 cm-2, the (12CO 6-5)/(12CO 3-2) line ratio $\approx$1.1 is reproduced at temperatures of 60-100 K. If the disc is large and sufficiently cold, the assumption that the two lines arise from the same region is not necessarily true, as the 3-2 line may be dominated by colder material at radii larger than the 6-5 line. It is therefore uncertain whether the line ratio reflects the difference in the temperature between the two vertical layers or in the temperature of the radial regions emitting the lines. We explore this further in Sect. 3.2.1 using more sophisticated radiative transfer tools and parametric models of disc structure.

The observed CO 6-5/3-2 integrated intensity ratio of 1.1 $\pm$ 0.6 (observations scaled by the beam size) is higher than the ratios close to 0.5 found for the discs around the T Tauri stars LkCa 15 and TW Hya by van Zadelhoff et al. (2001a) and Qi et al. (2006). The TW Hya 6-5/3-2 ratio has been interpreted as proof that the gas temperature is higher than that of the dust in the surface layers where gas and dust are not thermally coupled. Both UV radiation and X-rays have been invoked to provide the additional gas heating (Nomura & Millar 2005; Gorti & Hollenbach 2008; Jonkheid et al. 2004; Kamp & Dullemond 2004; Glassgold et al. 2004), although the X-rays can be neglected for Herbig Ae stars (Kamp et al. 2008).

The higher ratio found for HD 100546 implies higher gas temperatures than for the aforementioned discs around T Tauri stars. This is expected based on models where most of the gas heating comes from UV radiation from the central star, since the cooler stars have less UV radiation (e.g., Woitke et al. 2009). The gas temperature is also strongly affected by the PAH abundance in the disc. PAH features are prominent in the HD 100546 mid-infrared spectrum (e.g., Malfait et al. 1998), and have been shown to be spatially extended across the disc (Geers et al. 2007). In contrast, neither of the two T Tauri discs for which the12CO J= 6-5 line has been detected show PAH emission, implying typical PAH abundances that are a factor of 10-100 lower (Geers et al. 2006).

For the specific case of Herbig Ae discs, Jonkheid et al. (2007) have computed the gas heating and chemistry, as well as the resulting CO 6-5 and 3-2 line emission, starting from a set of dust disc models developed by Dullemond & Dominik (2005). Models with decreasing disc mass from 10-1 to 10 $^{-4}~M_\odot$ and decreasing dust/gas ratios from 10-2 to 10-6 (simulating grain growth and settling) were investigated. The UV radiation field was taken to be that of a B9.5 star (see Fig. 3 of Jonkheid et al. 2006), very close to what is expected for HD 100546. An accurate treatment of the shape of the UV field at short wavelengths, <1100 Å, is very important for correct calculation of the CO photodissociation and atomic carbon photoionisation rates, since the UV intensity of a B9 star is orders of magnitude weaker in this wavelength range than the (scaled) standard interstellar radiation field. The PAH abundance is taken to follow the dust/gas ratio, with an abundance of 10-7 for a normal dust/gas ratio of 0.01.

The 12CO J= 6-5/3-2 integrated intensity ratios computed by Jonkheid et al. (2007) summed over the full extent of the disc model are remarkably close to unity for the entire range of disc parameters investigated, consistent with our observations. The absolute 12CO and 13CO intensities for model B2 (a disc with a mass of 0.01 $M_{\odot}$ with a standard gas/dust ratio of 100, appropriate for HD 100546) are also within 40% of the observed values when scaled to the same source distance and beam size, indicating good agreement between models and observations.

3.2 Parametric disc models and origin of 12CO line emission

Submillimetre 12CO line emission is analysed using two different modelling approaches in the literature. Disc physical models, like the irradiated accretion disc models of Dullemond & Dominik (2005); D'Alessio et al. (2005), are especially suitable when the emission is spatially resolved. In this way, the three-dimensional structure of the disc can be investigated, for example when transitions of different optical depths are observed (Panic et al. 2008). For spatially unresolved observations, like those presented here, simplistic models with a limited number of free parameters are more appropriate for deriving some basic constraints on disc properties based on the line spectrum (method employed in, e.g., Dartois et al. 2003; Guilloteau & Dutrey 1998; van Zadelhoff et al. 2001a; Dent et al. 2005; Dutrey et al. 1994). Here we go a step further than in Sect. 3.1 in analysing the observations in the context of such models.

Density and temperature structure. Panic & Hogerheijde (2009) show that simple power-law disc models, with a disc mass $M=0.01~M_{\odot}$, an H2 surface density  $\Sigma \propto R^{-1}$, and temperature $T=T_{100}~{\rm K}~(R/{\rm 100~AU})^{-q}$ with T100=60 K and q=0.5, are useful tools for the analysis of low-J 12CO transitions from discs around gas-rich Herbig Ae stars. These models are vertically isothermal and have a vertical density structure in hydrostatic equilibrium. In the upper panels of Fig. 2 the temperature and surface density are shown of one of these models that has an outer radius of 400 AU. The 12CO abundance with respect to H2 is assumed to be 10-4, constant throughout the disc. Their analysis uses the assumption that the low-J 12CO lines arise from the upper disc layers (optically thick) and that their fluxes are dominated by the outermost disc regions. Therefore the spectral profiles of these lines are very sensitive to the extent of the emission, with strong and narrow lines indicating a large disc and weak and relatively wider double-peaked lines characteristic of small discs.

\begin{figure} \par\includegraphics[width=16cm,clip]{13709fg2.ps} \vspace*{3.6mm} \end{figure} Figure 2:

a): radial temperature structure of a disc model from Panic & Hogerheijde (2009), with $R_{\rm out}=$ 400 AU, T100= 60 K and q= 0.5 is shown with the black solid line. The radial dependence of the line flux ${\rm d}F(R)/{\rm d}F(R)_{\rm max}$ for the 12CO 3-2 line and 6-5 line arising from this model is overplotted using the dashed and the dotted lines, respectively. b): the surface density of the same model is shown with the black full line, and line fluxes as in  a). c) and  d): same as in  a) and  b), for the model with a shallower temperature profile, q=0.2 and its corresponding line fluxes as a function of radius.

Open with DEXTER

Spatial origin of the emission. The regions of the disc probed by the 12CO J= 6-5 and J= 3-2 lines are studied in van Zadelhoff et al. (2001b). Their Fig. 6 shows the locations of the surfaces of opacity $\tau=1$ for these lines in a cross-section of their disc model. At distances up to 150 AU from the star, the two transitions trace similar depths into the discs, but at larger radii they trace different vertical layers. The low-excitation 3-2 line probes deeper and colder layers than the 6-5 line. To assess whether both 6-5 and 3-2 line fluxes are dominated by the same radial disc region and to identify this region, we used the calculations of synthetic 6-5 and 3-2 spectra from disc models as described in the previous paragraph. We used two different temperature profiles with q=0.5 and q=0.2. The radial temperature and surface density structure of the two models is shown in Fig. 2. In this exercise we use the face-on disc orientation for simplicity, and calculate the emerging 12CO J= 6-5 and J= 3-2 line intensity along lines of sight at a number of different radial locations from the star. This is done using the molecular excitation and radiative transfer code RATRAN (Hogerheijde & van der Tak 2000) and CO-H2 collision rates from the Leiden Atomic and Molecular Database[*] (LAMBDA; Schöier et al. 2005). We are primarily interested in knowing if there is a single radial region dominating the line flux more than other regions or the flux is distributed more evenly with the radius. Also, we investigate how the regions of strongest line flux spatially compare for the two 12CO transitions. The contribution to the line flux from an annulus of a fixed width ${\rm d}R$ located at a distance R from the star is given by ${\rm d}F(R)=2~\pi~I(R)~R~{\rm d}R$, hence ${\rm d}F(R)\propto I(R)~R$. The relative contribution of the annuli to the total line flux arising from the model ${\rm d}F(R)/F_{\rm total}$ is also proportional to I(R) R, allowing us to directly compare how much one radial region contributes to the line flux with respect to another region. Figure 2 shows  ${\rm d}F(R)$ for the 12CO J= 6-5 and J= 3-2 in case q=0.5 (upper panels) and q=0.2 (lower panels). For ease of comparison between the two lines, the 12CO J= 3-2 and J= 6-5 line fluxes are plotted with respect to their maximum value  ${\rm d}F_{\rm max}$. For both temperature profiles, the outermost disc region (200-400 AU) is very important, especially for the 3-2 line. For this line the strongest flux arises at the largest radii due to the efficient excitation at low temperatures found in the outer disc combined with the large surface area. This result is in line with our previous results in Panic & Hogerheijde (2009) where the 12CO J = 3-2 spectra are used as probes of disc sizes. While the 6-5 line is also sensitive to the outer disc regions, the upper panels of Fig. 2 show that our model assumptions at 50 AU matter just as much as those at 400 AU. This is explained by the excitation energy of the 12CO J = 6-5 line, what is somewhat high for the 30-40 K temperatures found in the outer disc in the model with a steeper temperature profile q=0.5. The flatter temperature distribution with T >40 K in the outer disc allows the higher energy levels of 12CO to be efficiently populated, and in this case the outermost disc regions are clearly more important than 50-100 AU, merely because of the increase in the surface area with radius. Understanding the diagnostic potential of the high-J 12CO lines in discs is already useful for the analysis of spatially unresolved APEX observations (see Sect. 3.2.1), but also shows that the future sub-arcsecond observations of high-J 12CO lines in discs with ALMA could be used to probe the radial temperature structure in discs.

Disc mass. Our initial choice of parameters, described above, roughly matches the observed 12CO J= 6-5 and J= 3-2 line fluxes, confirming the outer radius of 400 AU seen in scattered light observations. We performed test calculations for a lower value of the disc mass of 10 $^{-3}~M_{\odot}$, i.e., ten times lower densities. The modelled line fluxes are found not to vary significantly with respect to the noise levels of our observations. Thus, in the disc mass range of 10-3-10 $^{-2}~M_{\odot}$, values found typically towards circumstellar discs and temperatures similar to those assumed here, the observed 12CO J= 6-5 and J= 3-2 lines are insensitive to disc mass. If we decrease the disc mass to allow the 12CO line emission in our models to be sensitive to the density variations in the midplane (i.e., emission dominated by the midplane), the disc midplane temperatures ($\leq$20 K in the outer disc) would be the only physically appropriate assumption for the temperature in the ``low-mass'' model. While the asymmetry in the line profile could be reproduced in this way, the resulting lines would be much weaker than observed. Given these considerations, we assume the two lines to be optically thick in the disc around HD 100546 and fix the disc mass to the initial value of 10 $^{-2}~M_{\odot}$. Our model is consistent with the thermal continuum emission of dust observed toward HD 100546 by Henning et al. (1998) and Wilner et al. (2003). In the further analysis, we fit the exact line shapes and intensities, with a particular focus on the line asymmetry seen in both J= 6-5 and J= 3-2 spectra.

3.2.1 12CO J = 6-5 and J = 3-2 line fit

The characteristic temperatures and masses derived in Sect. 3.1 for the disc around HD 100546 are similar to the models used in Panic & Hogerheijde (2009). The modelling of the 12CO 3-2 line of HD 100546 in Panic & Hogerheijde (2009) yielded very rough estimates of the outer radius of 300 AU and inclination of 35 $\hbox{$^\circ$ }$, consistent with the 400 AU and 50 $\hbox{$^\circ$ }$ from the scattered light. We use their models to fit the 12CO spectra, with T100 as free parameter and varying q from its initial value of 0.5 where necessary. We fix the outer radius and inclination to the observationally constrained values 400 AU and 50 $\hbox{$^\circ$ }$ (see Augereau et al. 2001, and the discussion of disc size in Sect. 2). We use the Keplerian velocity field around a 2.5 $M_{\odot}$ star. The dependence of the spectral profile on the stellar mass in the range 2.0-3.0 $M_{\odot}$ does not affect our fit significantly.

The inner radius is assumed to be 0.6 AU, close to the dust sublimation radius. Although an inner hole of 13 AU is found in HD 100546, its presence would not affect our results because the molecular lines observed are dominated by disc regions far beyond the inner tens of AU. Because the observed 12CO lines trace warm molecular material (>20 K) and are insensitive to the colder regions deeper in the disc, it is reasonable to neglect the 12CO freeze-out.

As above, the line emission is calculated using RATRAN. The synthetic images from our disc models are convolved with the beam sizes of the corresponding observations, and the spectra are extracted towards the centre of the image. Where temperature asymmetry is modelled (see 3.2.1), pairs of synthetic images with different T100 parameter were first combined - each contributing to a given side of the disc with respect to the semi-major axis. Following this, the combined image is convolved and spectra extracted. Although negligible for the molecular line transfer, dust continuum emission is included in the calculation and subtracted from the synthetic image cubes.

In view of the result that the line emission appears to be insensitive to the density, the observed line asymmetry cannot be caused by a density asymmetry in the disc, but rather reflects an asymmetry in the disc temperature structure. To model such an asymmetry, we use the fact that the emission at the two peaks in the line profile is dominated by the two sides of the disc with respect to the minor axis. We assume that one side of the disc is colder than the other. This may stem from a combined effect of an asymmetry in disc illumination and geometry. We use different T100 parameters for the two sides of the disc and explore different radial temperature distributions with $T \propto R^{-q}$, q=0-0.5.

The best fit to the 12CO J= 3-2 spectrum, seen in Fig. 1, is obtained by assuming the radial temperature structure with q=0.5, as in Panic & Hogerheijde (2009), and  $T_{\rm 100}$ of 60 and 70 K for the two disc sides, respectively. The line shape is consistent with the assumed Keplerian velocity field and disc inclination. The temperatures of 60-70 K compare well to the theoretical predictions of temperature in the regions where these lines saturate in discs (see Fig. 6 in van Zadelhoff et al. 2001a). They show that at 100 AU from the star, the J= 6-5 and J= 3-2 lines trace the same disc layer. Therefore we assume that the 12CO J= 6-5 line is emitted at  $T_{\rm 100}$ of 60 and 70 K, as found for the 3-2 line. We find the best fit to the 12CO J= 6-5 spectrum, shown in Fig. 1, assuming T= 60 K $\times$ $(R/100~{\rm AU})^{-0.3}$ on one side of the disc and T= 70 K $\times$ $(R/100~{\rm AU})^{-0.2}$ on the other. Thus, the radial temperature distribution in the 6-5 line emitting layer is flatter than in the 3-2 layer, consistent with the result of van Zadelhoff et al. (2001b) for the outer disc regions where the 6-5 line traces warmer material higher in the disc surface than the 3-2 line.

Therefore our 12CO line observations, although spatially unresolved, provide some information on the temperature radial distribution in the 6-5 and 3-2 line emitting layers and disc temperature asymmetry. It is important to note that slightly different sets of temperature values and slopes than those found here may provide equally good fits to these observations. However, any such model requires higher outer disc temperatures to fit the 6-5 line than required for the 3-2 line, clearly showing that the two lines arise from different vertical layers in the outermost disc regions.

\begin{figure} \par\includegraphics[width=8.5cm,clip]{13709fg3.eps}\vspace{-3mm} \end{figure} Figure 3:

Illustration of a warped disc structure, in an edge-on projection. Straight light and dark grey lines show the plane of the disc and of the inner warped disc component, respectively. The white curve shows the surface region with a higher temperature, heated well by the star, and the black curve shows the cooler region partly obscured by the warp. The image of the disc with a 400 AU radius is adapted from Quillen (2006), where a detailed warped disc model is presented.

Open with DEXTER

The difference in temperature between the two sides of the disc may be explained by a warped inner disc, as illustrated in Fig. 3, with the elevated side of the inner disc intercepting a fraction of stellar light that would otherwise reach the outer disc, while the opposite side is slightly more illuminated, with the inner part of the disc shifted downwards. The possibility of an inner warp is suggested in Quillen (2006) for HD 100546, with an inner component extending up to 200 AU inclined by $\approx$15 $\hbox{$^\circ$ }$ with respect to the outer component extending beyond that radius.

The temperature asymmetry is also possible if the disc has a different thickness on the two sides. This may happen in a disc with dust settling underway, if a planet or another body embedded in the disc stirs the dust back up. In this case, the ``stirred'' part of the disc intercepts more stellar light and becomes somewhat warmer. A companion body in HD 100546 is suggested in the literature, to explain the observed inner hole, gas kinematics, and spiral arms (Quillen et al. 2005; Acke & van den Ancker 2006b).

3.2.2 13CO J = 3-2 and 12CO J = 7-6 line fit

With the low signal to noise in our 13CO J= 3-2 spectrum (marginal detection at 2$\sigma$), it is not possible to assess the spectral line shape. The disc temperatures used to fit the 12CO lines in Sect. 3.2.1 overestimate the 13CO J= 3-2 line intensity. A lower temperature of the 13CO J= 3-2 emitting layers with respect to the 12CO J= 3-2 (as discussed in Sect. 3.2.2), a higher [12C]/[13C] isotopic ratio than the standard interstellar medium value of 77, or an overall decrease in CO abundance can explain our marginal detection of the 13CO line. Freeze-out and/or selective photodissociation processes in discs are known to affect the CO abundance and isotopic ratios. Detailed modelling of spatially resolved submillimetre line observations of 12CO and isotopologues would allow us to constrain the disc structure and the amount of molecular gas better, shedding light on these processes (e.g., Panic et al. 2008). Thereby, there are a number of uncertainties involved in the interpretation of our 13CO J= 3-2 data. For these reasons we limit the analysis to a comparison of our marginal detection to one simple model. A lower temperature of the 13CO emitting layers is expected because of a different optical depth of 13CO lines with respect to the 12CO. In Fig. 1 we compare our data to model calculations done assuming low temperatures, T100 = 30 K and q=0.5, an isotopic ratio [12C]/[13C] = 77, and all remaining model parameters as in the simple axially symmetric power-law model shown in Fig. 2 and described in Sect. 3.2.1. Considering that the 13CO line emission probes colder and deeper disc layers, it is likely that temperature asymmetry derived for the upper layers from 12CO lines is negligible for 13CO. High-sensitivity observations of 13CO or C18O submillimetre lines in the future may reveal symmetric line profiles.

The 12CO J= 7-6 and J= 6-5 transitions arise from energy levels at $E_{\rm k}$ of 155 K and 116 K, respectively, and are expected to trace disc layers that are roughly the same. Figure 1 shows how our best-fit model of the 6-5 line compares to the observed 7-6 line. The observed line flux is about 50% lower than the model prediction. This may be due to a slight mispointing of the telescope, as the size of the source (800 AU in diameter) is close to the 7 $\hbox{$.\!\!^{\prime\prime}$ }$7 beam size at this frequency. The noise levels of these observations are notably higher than in the 6-5 line observations, and spectral profile not well defined. For these reasons, we do not fit the 7-6 line further but use the more reliable 6-5 line fit results in our analysis.

3.3 Implications of the [C I] J = 2-1 non-detection

Figure 1 includes the high-quality spectrum around the [C I] J= 2-1 line at 809.344 GHz. No significant feature is detected down to 0.3 K rms in a 0.27 km s-1 velocity bin, implying a 2$\sigma$ limit on the integrated intensity of ${\approx}1.0$ K km s-1 over the range 0-10 km s-1 (same width as for the detected 12CO lines). Model B2 of Jonkheid et al. (2007) predicts integrated [C I] intensities, scaled to the distance of HD 100546, around 15-20 K km s-1, whereas the model intensities are even higher in discs with significant grain growth and settling (BL model series in the aforementioned paper). Thus, while the CO data appear entirely consistent with their sophisticated UV-heated disc atmosphere models, the [C I] data are clearly discrepant by an order of magnitude.

One possible solution could be that the radiation field contains more carbon ionising photons than assumed here, shifting the chemical balance from neutral to ionised atomic carbon. The predicted [C II] line intensities for the model discs of Jonkheid et al. (2007) are very low, <0.1 K km s-1, whereas they are more than an order of magnitude higher for discs around T Tauri stars with excess UV emission (Jonkheid et al. 2004), consistent with an efficient carbon ionisation around Herbig Ae stars. Such excess of UV emission could come from the disc-star accretion boundary layer. Indeed, HD 100546 is observed to undergo significant accretion, in spite of the known (dust) gap in the inner disc (Vieira et al. 1999). Also, the FUSE UV spectrum of HD 100546 shows significant flux, both in the continuum and in lines such as O VI, in the wavelength range where atomic carbon can be ionised (Lecavelier des Etangs et al. 2003, and G. Herczeg, priv. comm.). The continuum emission should also photodissociate CO, but not the stellar UV lines unless they overlap with the discrete CO photodissociating transitions. Also, because CO is self-shielding and the lines are optically thick, the CO emission may be less affected than that of atomic carbon by any extra UV. However, detailed chemical modelling is required to test this hypothesis.

[C II] 158 $\mu$m emission has been detected by Malfait et al. (1998) with ISO-LWS and most recently with the PACS instrument on the Herschel Space Observatory by Sturm et al. (2010). The ISO fluxes are up to an order of magnitude larger than the model fluxes by Jonkheid et al. (2007), but a significant fraction of the emission observed in the large ISO beam may be caused by general Galactic background emission. More detailed modelling is needed to pin down the origin of the absence of the [C I] 809 $\mu$m line and check whether the observed [C II] flux is consistent with enhanced UV radiation. Future observations of the [C I] J= 1-0 line at 492.161 GHz (3P1-3P0), which is generally brighter than the [C I] J= 2-1 line, may provide a stronger constraint on the amount of neutral atomic carbon in the disc. Observations of the [C II] line with Herschel-HIFI can directly resolve the line profile and determine whether this line indeed comes from the disc. Overall, the [C I] flux and [C I]/[C II] line ratio could become interesting diagnostics of the UV radiation field to which the disc is exposed.

4 Conclusions

We summarise our conclusions as follows.

  • We present evidence of warm molecular gas associated with the disc around HD 100546, in the regions within 400 AU from the star, successfully separated from more extended material in our CHAMP+ observations.

  • The gas kinematics are consistent with Keplerian rotation around an 2.5 $M_{\odot}$ star of a disc with a 400 AU radius, viewed at an inclination of 50 $\hbox{$^\circ$ }$ from face-on.

  • The 12CO (6-5)/(3-2) line ratio of 1.1 $\pm$ 0.6 is higher than measured towards discs around T Tauri stars, likely due to a more efficient gas heating of the disc containing PAHs by the stronger UV radiation from the B9 star.

  • Our data testify to the significant molecular gas reservoir in the disc, consistent with the total disc masses of more than 10 $^{-3}~M_{\odot}$. We exclude the possibility of a low-density disc and optically thin 12CO emission.

  • Line asymmetry seen in 12CO J= 6-5 and J= 3-2 lines is explained by a temperature asymmetry, with one side of the disc slightly colder than the other, possibly from partial obscuration of one side by a warped inner disc or a high disc rim. We exclude radial asymmetry, midplane density asymmetry, and mispointing as possible causes

  • Our modelling shows that, owing to the efficient heating of the disc gas by the star, both low-J and high-J 12CO lines are dominated by the outermost regions of the disc, though slightly different vertical disc layers with $\Delta T$ = 15-20 K.

  • We find that in ``colder'' discs where temperatures of the emitting regions are close to 20-30 K in the outer disc - colder Herbig Ae discs and especially T Tauri discs - the high-J lines probe a larger extent of the disc, starting from as little as 50 AU.

  • The non-detection of [C I] J= 2-1 line may indicate the presence of more carbon-ionising photons than assumed in the B9 model atmosphere.
Future observations with ALMA will allow us to characterise the disc around HD 100546 and spatially resolve its kinematics and structure. In particular, these observations will allow a detailed comparison between the spatial distribution of the gas traced by the rotational transitions of 12CO and its isotopologues, and the dust traced with the millimetre continuum emission. Herschel far-infrared data can probe even higher J 12CO transitions, as well as [O I] and [C II] lines originating from the disc surface. Because it is a bright, isolated source suspected of harbouring a planet, the disc around HD 100546 is one of the prime targets to probe disc structure in the early planet-building phases.

Acknowledgements
The research of O.P. and M.R.H. is supported through a VIDI grant from The Netherlands Organisation for Scientific Research (NWO). We thank L. Kristensen for advice regarding the APEX data reduction, and C. M. Wright for useful discussions. We thank the APEX staff for assistance during the observations. Construction of CHAMP+ is a collaboration between the Max-Planck-Institut für Radioastronomie Bonn, The Netherlands Institute for Space Research (SRON), The Netherlands Research School for Astronomy (NOVA), and the Kavli Institute of Nanoscience at Delft University of Technology, with support from NWO grant 60006331010.

References

  1. Acke, B., & van den Ancker, M. E. 2006a, A&A, 457, 171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Acke, B., & van den Ancker, M. E. 2006b, A&A, 449, 267 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Ardila, D. R., Golimowski, D. A., Krist, J. E., et al. 2007, ApJ, 665, 512 [NASA ADS] [CrossRef] [Google Scholar]
  4. Augereau, J. C., Lagrange, A. M., Mouillet, D., & Ménard, F. 2001, A&A, 365, 78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Beckwith, S., & Sargent, A. 1987, in Circumstellar Matter, ed. I. Appenzeller, & C. Jordan, IAU Symp., 122, 81 [Google Scholar]
  6. Bouwman, J., de Koter, A., Dominik, C., & Waters, L. B. F. M. 2003, A&A, 401, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Brittain, S. D., Najita, J. R., & Carr, J. S. 2009, ApJ, 702, 85 [NASA ADS] [CrossRef] [Google Scholar]
  8. Calvet, N., D'Alessio, P., Hartmann, L., et al. 2002, ApJ, 568, 1008 [NASA ADS] [CrossRef] [Google Scholar]
  9. D'Alessio, P., Merín, B., Calvet, N., Hartmann, L., & Montesinos, B. 2005, Rev. Mex. Astron. Astrofis., 41, 61 [NASA ADS] [Google Scholar]
  10. Dartois, E., Dutrey, A., & Guilloteau, S. 2003, A&A, 399, 773 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Dent, W. R. F., Greaves, J. S., & Coulson, I. M. 2005, MNRAS, 359, 663 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dullemond, C. P., & Dominik, C. 2005, A&A, 434, 971 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Dutrey, A., Guilloteau, S., & Simon, M. 1994, A&A, 286, 149 [NASA ADS] [Google Scholar]
  14. Geers, V. C., Augereau, J.-C., Pontoppidan, K. M., et al. 2006, A&A, 459, 545 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Geers, V. C., van Dishoeck, E. F., Visser, R., et al. 2007, A&A, 476, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Glassgold, A. E., Najita, J., & Igea, J. 2004, ApJ, 615, 972 [NASA ADS] [CrossRef] [Google Scholar]
  17. Gorti, U., & Hollenbach, D. 2008, ApJ, 683, 287 [NASA ADS] [CrossRef] [Google Scholar]
  18. Grady, C. A., Polomski, E. F., Henning, T., et al. 2001, AJ, 122, 3396 [NASA ADS] [CrossRef] [Google Scholar]
  19. Greaves, J. S., Mannings, V., & Holland, W. S. 2000, Icarus, 143, 155 [NASA ADS] [CrossRef] [Google Scholar]
  20. Guilloteau, S., & Dutrey, A. 1998, A&A, 339, 467 [NASA ADS] [Google Scholar]
  21. Guilloteau, S., & Dutrey, A. 2008, Ap&SS, 313, 95 [NASA ADS] [CrossRef] [Google Scholar]
  22. Güsten, R., Nyman, L. Å., Schilke, P., et al. 2006, A&A, 454, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Henning, T., Burkert, A., Launhardt, R., Leinert, C., & Stecklum, B. 1998, A&A, 336, 565 [NASA ADS] [Google Scholar]
  24. Hernández, J., Calvet, N., Briceño, C., et al. 2007, ApJ, 671, 1784 [NASA ADS] [CrossRef] [Google Scholar]
  25. Hillenbrand, L. A. 2008, Phys. Scr. Vol. T, 130, 014024 [NASA ADS] [CrossRef] [Google Scholar]
  26. Hogerheijde, M. R., & van der Tak, F. F. S. 2000, A&A, 362, 697 [NASA ADS] [Google Scholar]
  27. Hollenbach, D. J., Yorke, H. W., & Johnstone, D. 2000, Protostars and Planets IV, 401 [Google Scholar]
  28. Isella, A., Testi, L., Natta, A., et al. 2007, A&A, 469, 213 [CrossRef] [EDP Sciences] [Google Scholar]
  29. Jonkheid, B., Faas, F. G. A., van Zadelhoff, G.-J., & van Dishoeck, E. F. 2004, A&A, 428, 511 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Jonkheid, B., Kamp, I., Augereau, J., & van Dishoeck, E. F. 2006, A&A, 453, 163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Jonkheid, B., Dullemond, C. P., Hogerheijde, M. R., & van Dishoeck, E. F. 2007, A&A, 463, 203 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Kamp, I., & Dullemond, C. P. 2004, ApJ, 615, 991 [NASA ADS] [CrossRef] [Google Scholar]
  33. Kamp, I., Freudling, W., Robberto, M., Chengalur, J., & Keto, E. 2008, Phys. Scr. Vol. T, 130, 014013 [NASA ADS] [CrossRef] [Google Scholar]
  34. Klein, B., Philipp, S. D., Krämer, I., et al. 2006, A&A, 454, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Koerner, D. W., Sargent, A. I., & Beckwith, S. V. W. 1993, Icarus, 106, 2 [NASA ADS] [CrossRef] [Google Scholar]
  36. Lecavelier des Etangs, A., Deleuil, M., Vidal-Madjar, A., et al. 2003, A&A, 407, 935 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Malfait, K., Waelkens, C., Waters, L. B. F. M., et al. 1998, A&A, 332, L25 [NASA ADS] [Google Scholar]
  38. Nomura, H., & Millar, T. J. 2005, A&A, 438, 923 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
  39. Panic, O., & Hogerheijde, M. R. 2009, A&A, 508, 707 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Panic, O., Hogerheijde, M. R., Wilner, D., & Qi, C. 2008, A&A, 491, 219 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Panic, O., Hogerheijde, M. R., Wilner, D., & Qi, C. 2009, A&A, 501, 269 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Pantin, E., Waelkens, C., & Lagage, P. O. 2000, A&A, 361, L9 [NASA ADS] [Google Scholar]
  43. Pinte, C., Padgett, D. L., Ménard, F., et al. 2008, A&A, 489, 633 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Qi, C., Ho, P. T. P., Wilner, D. J., et al. 2004, ApJ, 616, L11 [NASA ADS] [CrossRef] [Google Scholar]
  45. Qi, C., Wilner, D. J., Calvet, N., et al. 2006, ApJ, 636, L157 [NASA ADS] [CrossRef] [Google Scholar]
  46. Quillen, A. C. 2006, ApJ, 640, 1078 [NASA ADS] [CrossRef] [Google Scholar]
  47. Quillen, A. C., Varnière, P., Minchev, I., & Frank, A. 2005, AJ, 129, 2481 [NASA ADS] [CrossRef] [Google Scholar]
  48. Raman, A., Lisanti, M., Wilner, D. J., Qi, C., & Hogerheijde, M. 2006, AJ, 131, 2290 [NASA ADS] [CrossRef] [Google Scholar]
  49. Schöier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 369 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Simon, M., Dutrey, A., & Guilloteau, S. 2000, ApJ, 545, 1034 [NASA ADS] [CrossRef] [Google Scholar]
  51. Sturm, B., Bouwman, J. Henning, Th., et al. 2010, A&A, 518, L129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Tannirkulam, A., Monnier, J. D., Harries, T. J., et al. 2008, ApJ, 689, 513 [NASA ADS] [CrossRef] [Google Scholar]
  53. Thé, P. S., de Winter, D., & Perez, M. R. 1994, A&AS, 104, 315 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
  54. Thi, W. F., van Dishoeck, E. F., Blake, G. A., et al. 2001, ApJ, 561, 1074 [NASA ADS] [CrossRef] [Google Scholar]
  55. Thi, W.-F., van Zadelhoff, G.-J., & van Dishoeck, E. F. 2004, A&A, 425, 955 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. van den Ancker, M. E., de Winter, D., & Tjin A Djie, H. R. E. 1998, A&A, 330, 145 [NASA ADS] [Google Scholar]
  57. van der Plas, G., van den Ancker, M. E., Acke, B., et al. 2009, A&A, 500, 1137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. van der Tak, F. F. S., Black, J. H., Schöier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, A&A, 468, 627 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. van Kempen, T. A., van Dishoeck, E. F., Brinch, C., & Hogerheijde, M. R. 2007, A&A, 461, 983 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. van Zadelhoff, G.-J., van Dishoeck, E. F., Thi, W.-F., & Blake, G. A. 2001a, A&A, 377, 566 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. van Zadelhoff, G.-J., van Dishoeck, E. F., Thi, W.-F., & Blake, G. A. 2001b, A&A, 377, 566 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Vieira, S. L. A., Pogodin, M. A., & Franco, G. A. P. 1999, A&A, 345, 559 [NASA ADS] [Google Scholar]
  63. Waelkens, C., Waters, L. B. F. M., de Graauw, M. S., et al. 1996, A&A, 315, L245 [NASA ADS] [Google Scholar]
  64. Wilner, D. J., Bourke, T. L., Wright, C. M., et al. 2003, ApJ, 596, 597 [NASA ADS] [CrossRef] [Google Scholar]
  65. Wilson, T. L., & Rood, R. 1994, ARA&A, 32, 191 [NASA ADS] [CrossRef] [Google Scholar]
  66. Woitke, P., Kamp, I., & Thi, W.-F. 2009, A&A, 501, 383 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Footnotes

... telescope[*]
This publication is based on data acquired with the Atacama Pathfinder Experiment (APEX). APEX is a collaboration between the Max-Planck-Institut für Radioastronomie, the European Southern Observatory, and the Onsala Space Observatory.
... tool[*]
http://www.strw.leidenuniv.nl/ moldata/radex.html
... Database[*]
http://www.strw.leidenuniv.nl/ moldata/

All Tables

Table 1:   Integrated intensities of the observed spectral lines.

Table 2:   The 12CO J= 6-5 and J= 7-6 integrated line intensities at offset positions.

All Figures

  \begin{figure} \par\includegraphics[width=8.8cm,clip]{13709fg1.ps} \end{figure} Figure 1:

From top to bottom: spectra of the 12CO J= 7-6, J= 6-5, J= 3-2, 13CO J= 3-2 and CI J= 2-1 lines observed towards HD 100546 shown with black lines plotted as a histogram. The 12CO and 13CO spectra computed from our models of the line emitting layers of the disc (see Sects. 3.2.1 and 3.2.2) are shown with continuous red lines. At 100 AU from the star, our vertically isothermal models of the 12CO line emission have temperatures T100 = 60-70 K, while a lower temperature T100 = 30 K is used for the 13CO line.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=16cm,clip]{13709fg2.ps} \vspace*{3.6mm} \end{figure} Figure 2:

a): radial temperature structure of a disc model from Panic & Hogerheijde (2009), with $R_{\rm out}=$ 400 AU, T100= 60 K and q= 0.5 is shown with the black solid line. The radial dependence of the line flux ${\rm d}F(R)/{\rm d}F(R)_{\rm max}$ for the 12CO 3-2 line and 6-5 line arising from this model is overplotted using the dashed and the dotted lines, respectively. b): the surface density of the same model is shown with the black full line, and line fluxes as in  a). c) and  d): same as in  a) and  b), for the model with a shallower temperature profile, q=0.2 and its corresponding line fluxes as a function of radius.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=8.5cm,clip]{13709fg3.eps}\vspace{-3mm} \end{figure} Figure 3:

Illustration of a warped disc structure, in an edge-on projection. Straight light and dark grey lines show the plane of the disc and of the inner warped disc component, respectively. The white curve shows the surface region with a higher temperature, heated well by the star, and the black curve shows the cooler region partly obscured by the warp. The image of the disc with a 400 AU radius is adapted from Quillen (2006), where a detailed warped disc model is presented.

Open with DEXTER
In the text

Copyright ESO 2010

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.