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 |
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
7
to 18
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
of molecular gas in a disc of
400 AU radius in
Keplerian rotation around a 2.5
star, seen at an inclination of 50
.
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 ( =
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
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
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
resolution,
with the flux of 36
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
from the star viewed at an inclination of 50
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 (
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
7
to 18
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
25
4
and Dec = -70
11
41
(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
at 691 GHz and 7
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
, 0
)
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
30
at both frequencies. Pointing was performed directly prior to and after
the on-source integration providing an accuracy better than 3
.
The pointing source was 07454-7112, about 20
away from HD 100546, at the same airmass.
![]() |
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
1
at 345.796 GHz and 18
9
at 330.588 GHz, respectively. Calibration is uncertain by
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
in terms of the integrated intensity. The 13CO J=
3-2 line is marginally detected at 2
,
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
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
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
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
7 and 9
).
We obtain the following integrated intensity ratios, corrected for beam
dilution (scaled to the same beam): 12CO lines
(6-5)/(3-2) = 1.1
0.6 and (12CO 3-2)/(13CO 3-2)
2.8
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
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
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
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
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
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
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
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
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 ,
an H2 surface
density
,
and temperature
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.
![]() |
Figure 2:
a): radial temperature structure of a disc
model from Panic &
Hogerheijde (2009), with |
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
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
located at a distance R from the star is
given by
,
hence
.
The relative contribution of the annuli to the total line flux arising
from the model
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
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
.
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
,
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
,
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 (
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
.
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
,
consistent with the 400 AU and 50
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
(see Augereau et al. 2001,
and the discussion of disc size in Sect. 2). We use
the Keplerian velocity field around a 2.5
star. The dependence of the spectral profile on the stellar mass in the
range 2.0-3.0
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 ,
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
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
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
on one side of the disc and T=
70 K
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.
![]() |
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 15
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), 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
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
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 limit
on the integrated intensity of
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 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
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
star of a disc with a 400 AU radius, viewed at an inclination of 50
from face-on.
- The 12CO (6-5)/(3-2) line ratio of
1.1
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
. 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
= 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.
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
- Acke, B., & van den Ancker, M. E. 2006a, A&A, 457, 171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Acke, B., & van den Ancker, M. E. 2006b, A&A, 449, 267 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ardila, D. R., Golimowski, D. A., Krist, J. E., et al. 2007, ApJ, 665, 512 [NASA ADS] [CrossRef] [Google Scholar]
- Augereau, J. C., Lagrange, A. M., Mouillet, D., & Ménard, F. 2001, A&A, 365, 78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Beckwith, S., & Sargent, A. 1987, in Circumstellar Matter, ed. I. Appenzeller, & C. Jordan, IAU Symp., 122, 81 [Google Scholar]
- Bouwman, J., de Koter, A., Dominik, C., & Waters, L. B. F. M. 2003, A&A, 401, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Brittain, S. D., Najita, J. R., & Carr, J. S. 2009, ApJ, 702, 85 [NASA ADS] [CrossRef] [Google Scholar]
- Calvet, N., D'Alessio, P., Hartmann, L., et al. 2002, ApJ, 568, 1008 [NASA ADS] [CrossRef] [Google Scholar]
- D'Alessio, P., Merín, B., Calvet, N., Hartmann, L., & Montesinos, B. 2005, Rev. Mex. Astron. Astrofis., 41, 61 [NASA ADS] [Google Scholar]
- Dartois, E., Dutrey, A., & Guilloteau, S. 2003, A&A, 399, 773 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dent, W. R. F., Greaves, J. S., & Coulson, I. M. 2005, MNRAS, 359, 663 [NASA ADS] [CrossRef] [Google Scholar]
- Dullemond, C. P., & Dominik, C. 2005, A&A, 434, 971 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dutrey, A., Guilloteau, S., & Simon, M. 1994, A&A, 286, 149 [NASA ADS] [Google Scholar]
- Geers, V. C., Augereau, J.-C., Pontoppidan, K. M., et al. 2006, A&A, 459, 545 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Geers, V. C., van Dishoeck, E. F., Visser, R., et al. 2007, A&A, 476, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Glassgold, A. E., Najita, J., & Igea, J. 2004, ApJ, 615, 972 [NASA ADS] [CrossRef] [Google Scholar]
- Gorti, U., & Hollenbach, D. 2008, ApJ, 683, 287 [NASA ADS] [CrossRef] [Google Scholar]
- Grady, C. A., Polomski, E. F., Henning, T., et al. 2001, AJ, 122, 3396 [NASA ADS] [CrossRef] [Google Scholar]
- Greaves, J. S., Mannings, V., & Holland, W. S. 2000, Icarus, 143, 155 [NASA ADS] [CrossRef] [Google Scholar]
- Guilloteau, S., & Dutrey, A. 1998, A&A, 339, 467 [NASA ADS] [Google Scholar]
- Guilloteau, S., & Dutrey, A. 2008, Ap&SS, 313, 95 [NASA ADS] [CrossRef] [Google Scholar]
- Güsten, R., Nyman, L. Å., Schilke, P., et al. 2006, A&A, 454, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Henning, T., Burkert, A., Launhardt, R., Leinert, C., & Stecklum, B. 1998, A&A, 336, 565 [NASA ADS] [Google Scholar]
- Hernández, J., Calvet, N., Briceño, C., et al. 2007, ApJ, 671, 1784 [NASA ADS] [CrossRef] [Google Scholar]
- Hillenbrand, L. A. 2008, Phys. Scr. Vol. T, 130, 014024 [NASA ADS] [CrossRef] [Google Scholar]
- Hogerheijde, M. R., & van der Tak, F. F. S. 2000, A&A, 362, 697 [NASA ADS] [Google Scholar]
- Hollenbach, D. J., Yorke, H. W., & Johnstone, D. 2000, Protostars and Planets IV, 401 [Google Scholar]
- Isella, A., Testi, L., Natta, A., et al. 2007, A&A, 469, 213 [CrossRef] [EDP Sciences] [Google Scholar]
- 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]
- Jonkheid, B., Kamp, I., Augereau, J., & van Dishoeck, E. F. 2006, A&A, 453, 163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jonkheid, B., Dullemond, C. P., Hogerheijde, M. R., & van Dishoeck, E. F. 2007, A&A, 463, 203 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kamp, I., & Dullemond, C. P. 2004, ApJ, 615, 991 [NASA ADS] [CrossRef] [Google Scholar]
- Kamp, I., Freudling, W., Robberto, M., Chengalur, J., & Keto, E. 2008, Phys. Scr. Vol. T, 130, 014013 [Google Scholar]
- Klein, B., Philipp, S. D., Krämer, I., et al. 2006, A&A, 454, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Koerner, D. W., Sargent, A. I., & Beckwith, S. V. W. 1993, Icarus, 106, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Lecavelier des Etangs, A., Deleuil, M., Vidal-Madjar, A., et al. 2003, A&A, 407, 935 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Malfait, K., Waelkens, C., Waters, L. B. F. M., et al. 1998, A&A, 332, L25 [NASA ADS] [Google Scholar]
- Nomura, H., & Millar, T. J. 2005, A&A, 438, 923 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
- Panic, O., & Hogerheijde, M. R. 2009, A&A, 508, 707 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Panic, O., Hogerheijde, M. R., Wilner, D., & Qi, C. 2008, A&A, 491, 219 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Panic, O., Hogerheijde, M. R., Wilner, D., & Qi, C. 2009, A&A, 501, 269 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pantin, E., Waelkens, C., & Lagage, P. O. 2000, A&A, 361, L9 [NASA ADS] [Google Scholar]
- Pinte, C., Padgett, D. L., Ménard, F., et al. 2008, A&A, 489, 633 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Qi, C., Ho, P. T. P., Wilner, D. J., et al. 2004, ApJ, 616, L11 [NASA ADS] [CrossRef] [Google Scholar]
- Qi, C., Wilner, D. J., Calvet, N., et al. 2006, ApJ, 636, L157 [NASA ADS] [CrossRef] [Google Scholar]
- Quillen, A. C. 2006, ApJ, 640, 1078 [NASA ADS] [CrossRef] [Google Scholar]
- Quillen, A. C., Varnière, P., Minchev, I., & Frank, A. 2005, AJ, 129, 2481 [NASA ADS] [CrossRef] [Google Scholar]
- Raman, A., Lisanti, M., Wilner, D. J., Qi, C., & Hogerheijde, M. 2006, AJ, 131, 2290 [NASA ADS] [CrossRef] [Google Scholar]
- 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]
- Simon, M., Dutrey, A., & Guilloteau, S. 2000, ApJ, 545, 1034 [NASA ADS] [CrossRef] [Google Scholar]
- Sturm, B., Bouwman, J. Henning, Th., et al. 2010, A&A, 518, L129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tannirkulam, A., Monnier, J. D., Harries, T. J., et al. 2008, ApJ, 689, 513 [NASA ADS] [CrossRef] [Google Scholar]
- Thé, P. S., de Winter, D., & Perez, M. R. 1994, A&AS, 104, 315 [Google Scholar]
- Thi, W. F., van Dishoeck, E. F., Blake, G. A., et al. 2001, ApJ, 561, 1074 [NASA ADS] [CrossRef] [Google Scholar]
- Thi, W.-F., van Zadelhoff, G.-J., & van Dishoeck, E. F. 2004, A&A, 425, 955 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van den Ancker, M. E., de Winter, D., & Tjin A Djie, H. R. E. 1998, A&A, 330, 145 [NASA ADS] [Google Scholar]
- 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]
- 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]
- 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]
- 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]
- 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]
- Vieira, S. L. A., Pogodin, M. A., & Franco, G. A. P. 1999, A&A, 345, 559 [NASA ADS] [Google Scholar]
- Waelkens, C., Waters, L. B. F. M., de Graauw, M. S., et al. 1996, A&A, 315, L245 [NASA ADS] [Google Scholar]
- Wilner, D. J., Bourke, T. L., Wright, C. M., et al. 2003, ApJ, 596, 597 [NASA ADS] [CrossRef] [Google Scholar]
- Wilson, T. L., & Rood, R. 1994, ARA&A, 32, 191 [NASA ADS] [CrossRef] [Google Scholar]
- 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
![]() |
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 |
![]() |
Figure 2:
a): radial temperature structure of a disc
model from Panic &
Hogerheijde (2009), with |
Open with DEXTER | |
In the text |
![]() |
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.