A&A 418, 177-184 (2004)
DOI: 10.1051/0004-6361:20034331
R. van Boekel1,2 - L. B. F. M. Waters1,3 - C. Dominik1 - C. P. Dullemond4 - A. G. G. M. Tielens5 - A. de Koter1
1 - Astronomical Institute "Anton Pannekoek'', University of
Amsterdam, Kruislaan 403,
1098 SJ Amsterdam, The Netherlands
2 -
European Southern Observatory, Karl-Schwarzschildstrasse 2,
85748 Garching bei München, Germany
3 -
Instituut voor Sterrenkunde, Katholieke Universiteit Leuven,
Celestijnenlaan 200B, 3001 Heverlee, Belgium
4 -
Max-Planck-Institut für Astrophysik, Karl-Schwarzschildstrasse 1,
Postfach 1317, 85748 Garching bei München, Germany
5 -
Kapteyn Astronomical Institute, PO Box 800, 9700 AV
Groningen, The Netherlands
Received 17 September 2003 / Accepted 23 January 2004
Abstract
We present new mid-infrared spectroscopy of the emission
from warm circumstellar dust grains in the Herbig Ae stars HD 100546,
HD 97048 and HD 104237, with a spatial resolution of
0
9. We find that the emission in the UIR bands at 8.6,
11.3 and (HD 97048 only) 12.7
m is extended in the first two
sources. The continuum emission is resolved in HD 97048 and
possibly in HD 100546. HD 104237 is not spatially resolved in our
observations. We find that the UIR emission in HD 100546 and HD 97048
is extended on a scale of (several) 100 AU, corresponding to the outer
disk scale in flaring disk models. Small carbonaceous
particles are the dominant source of opacity in the HD 97048 disk.
Key words: circumstellar matter - stars: pre-main sequence - infrared: ISM - ISM: lines and bands
Herbig Ae/Be stars (HAEBEs, Herbig 1960)
are intermediate mass
pre-main-sequence stars surrounded by material which is left from the
star formation process. A sub-group of mostly late B and A-F type
HAEBE stars (hereafter HAEs) show little or no optical extinction and
usually low mass accretion rates, as derived from
radio analysis (Skinner et al. 1993)
and the lack of significant
veiling in optical spectra. The spatial
distribution of the circumstellar material of these HAEs is a matter
of debate. While the millimeter emission observed from some stars
originates from a compact dust disk which appears to be keplerian
(Mannings & Sargent 1997), near-IR
interferometric observations seem consistent with the hottest dust
having a more spherical distribution (Millan-Gabet et al. 2001). Only recently very high
resolution interferometric measurements at 2.2 m have provided
strong evidence for non-circular symmetry on sub-AU scales,
indicating disk like geometries in the innermost regions of HAE
systems (Eisner et al. 2003).
Optical and near-infrared images in scattered light argue for
disk geometries on large scales of 100 AU
(Grady et al. 2001;
Augereau et al. 2001).
Thus, evidence is growing that HAE stars indeed have disks.
The mid-infrared spectra of HAEs are dominated by a strong continuum caused by warm grains in the immediate vicinity of the star, and solid state emission bands from silicates and carbonaceous material (Meeus et al. 2001a; Natta et al. 2001). The IR emission is believed to originate mostly from the surface layer of a passively heated dusty disk (Chiang & Goldreich 1997). Indeed, recent hydrostatic disk models are successful in explaining the spectra of HAEs in the context of such disk models, without requiring other spatial components, such as a spherical halo (Dullemond et al. 2001).
The location of the dust
and the dust opacities determine the structure and emission properties
of the disk. Direct studies of the inner disk structure require
interferometric measurements.
On larger spatial scales, where we can probe the outermost
regions of the disk using single telescopes, the blackbody temperatures
are much too low for significant 10 m emission to be produced.
Such 10
m emission could only arise from
a population of very small "super-heated'' grains (VSG),
and emission from large carbonaceous molecules
(Polycyclic Aromatic Hydrocarbons, or PAHs) that can be excited
by single photon absorption. If present in the outer disk, these
molecules offer a unique probe of the outer disk structure since
they can radiate at all stello-centric distances, irrespective
of the local blackbody temperature.
Meeus et al. (2001b) noted that based on
the IR Spectral Energy Distribution (SED), HAEs can be divided into
two main groups: group I sources that have a very strong, rising
IR excess peaking around 60 m, and group II sources displaying
a more moderate and less steeply rising IR excess.
In the disk hypothesis, group I sources show a flaring outer disk
geometry, such that the disk can intercept and re-process stellar
radiation out to large stello-centric radii. In group II sources,
the outer disk is not flaring, leading to lower temperatures. It is
therefore expected that group I sources will show more extended
emission and the disk emission in group II sources will be
concentrated close to the central star.
A critical test of disk models requires high angular resolution observations at IR and millimeter wavelengths. Mid-infrared imaging using single-dish telescopes has shown that several HAEs are extended on a scale of 0.5 to 2 arcsec, corresponding to about 100 AU (Grady et al. 2001; Fischer et al. 2000; Jayawardhana et al. 2001; Ressler & Barsony 2003).
Long slit spectroscopy is ideally suited for detecting spatially resolved structures emitting in narrow spectral regions (i.e. emission lines/bands), when observing close to the resolution limit. The profile of the spectrum in the spatial direction is measured instantaneously in the emission band and in the continuum next to it. Therefore, a relative increase in size can be detected with high confidence, even if the response of the telescope to a point source (the point spread function, PSF) is not known accurately. However, a long-slit spectrum yields spatial information only along the direction of the slit.
Here, we report on spatially and spectrally resolved mid-IR observations observations of three HAE stars, HD 97048, HD 100546 (both group I), and HD 104237 (group II). Both group I stars show prominent emission in the Unidentified InfraRed (UIR) bands (Aitken & Roche 1981; Waelkens et al. 1996), usually attributed to PAHs. HD 97048 is in a reflection nebula and has been shown to have extended mid-IR emission from UIR bands (e.g. Siebenmorgen et al. 2000) on a scale of about 10 arcsec.
Long slit infrared spectra in the 10 m atmospheric window were
taken on 27 December 2001 and 19 March 2003
with the Thermal Infrared
Multi Mode Instrument 2 (TIMMI2, Reimann et al. 1998), mounted at the 3.60 m telescope at
the ESO La Silla observatory. We used the low resolution
(
)
N band grism and a slit width of 1.2 arcsec, the
pixel scale in the spectroscopic mode of TIMMI2 is 0.45 arcsec.
To correct for the strong atmospheric and instrumental background at
10
m, we employed chopping and nodding, using a
+10
chop throw north-south, and a -10
nod
throw north-south.
This renders us insensitive to diffuse emission on scales larger than 5
.
Spectroscopic standard stars were observed
regularly and were used to correct for the non-uniform
atmospheric transmission.
A summary of the observations is given in Table 1.
Table 1:
Summary of the observations. The fourth column lists the
seeing at 0.5 m as measured by the DIMM monitor. Distances
are calculated from Hipparcos parallaxes.
The (absolute) photometric accuracy of our observations is
about 15%. In order to faciliate the comparison of the TIMMI 2
spectra and ISO-SWS spectra, simple scaling factors (1)
have been applied to the TIMMI 2 spectra, such as to give the best
match with ISO (note that the absolute accuracy of ISO is about the
same as that of our TIMMI 2 data).
![]() |
Figure 1: FWHM2 of the spatial profile as a function of wavelength of HD 100546 and the model PSF ( lower panel). The solid grey curves show polynomial fits to the FWHM2 of the PSF and the continuum (i.e. the emission outside the UIR bands). The fitted continuum width is indicated with the dashed grey curve. The vertical dashed lines indicate the wavelengths of the UIR bands. The upper panel shows the ISO SWS spectrum of HD 100546 in black and the ground based spectrum in grey. |
Open with DEXTER |
![]() |
Figure 2: FWHM2 of the spatial profile as a function of wavelength of HD 97048 and the observed psf ( lower panel). The solid grey curves show polynomial fits to the FWHM2 of the PSF and the continuum (i.e. the emission outside the UIR bands). The fitted continuum width is indicated with the dashed grey curve. The vertical dashed lines indicate the wavelengths of the UIR bands. The upper panel shows the ISO SWS spectrum of HD 97048 in black and the ground based spectrum in grey. |
Open with DEXTER |
![]() |
Figure 3: FWHM2 of the spatial profile as a function of wavelength of HD 104237 and the model PSF (lower panel). The solid grey curve shows a polynomial fit to the FWHM2 of the PSF. The vertical dashed lines indicate the wavelengths of the UIR bands. The upper panel shows the ISO SWS spectrum of HD 104237 in black and the ground based spectrum in grey. |
Open with DEXTER |
We use this method to determine the PSF during the observations of
HD 100546 and HD 104237. The method assumes that the static
abberations due to the telescope optics are independent of telescope
elevation. Unfortunately, no measurement of a calibration star is
available at airmasses comparable to the HD 100546 and HD 104237
measurements. However, the spatial profile measured in our observation
of HD 104237 agrees very well with the predicted PSF. If the spatial
extent of HD 104237 is 1
(as one may expect, based on
modeling of its SED, and the non-detection in scattered light images
(Danks et al. 2001), contrary to HD 100546
and several other group I HAe stars (Pantin et al. 2000, Grady et al. 2001) this measurement shows that our
assumptions are valid. Note that HD 97048 which, according to the SED
and the spatially resolved data in this work, should show extended
structure in scattered light images, has not yet been observed in this
fashion. The calibration stars used for the HD 97048 were taken just
before and after the science measurement, at nearly identical airmass
and under indentical seeing conditions, and we have therefore a direct
determination of the PSF.
Grady et al. (2001) (hereafter GR01) find evidence
for extended emission at 11.7 m using a narrow-band filter, which
they attribute to either UIBs or silicates. We show here that the
extended emission is in fact due to UIBs. The measured FWHM in the
GR01 image is 1
03, slightly larger than
the mean FWHM in the spatial direction
of our longslit spectrum in the passband of the 11.7
m filter, which
is 0
94. GR01 find that over the whole N-band the source is
essentially unresolved, which is consistent with our findings: if
the spatial extent that is tentatively detected here is real,
an observation at an 8 m class telescope would be needed to make
a solid detection.
Due to the atmospheric cut-off at 7.8
m we only see the
red wing of the strong 7.7
m UIR band. Even though the PSF is not
fully understood here, the very strong rise in FWHM shortward of
8.4
m indicates that also the UIR emission here is
resolved in our data. Assuming a PSF FWHM of 0
95 and an UIR
contribution of 65%, a rough estimate of the UIR FWHM
in the 7.7-7.9
m complex yields a value of 140 AU. Since
conditions did not allow measurements at wavelengths shortward
of 8.0 micron, this estimate involves extrapolation of the profile,
adding considerable uncertainty. We note however, that the derived size
is consistent with the values measured in the 8.6 and 11.3
m
bands.
Despite the strong silicate emission that has significant spectral sub-structure due to crystalline material, we see no wavelength dependence of the FWHM outside the UIR emission bands. The emission from both the crystalline and amorphous silicates must arise from a physically much smaller region than the UIR emission.
Assuming a Gaussian light distribution for both the continuum and UIR
feature emission, we deduce a FWHM of
AU for the continuum
emission at 10.5
m, and FWHM values of
,
and
AU for the 8.6 11.3 and 12.7
m UIR bands,
respectively. We estimate the contribution of the 8.6, 11.3 and
12.7
m UIR bands to the total flux to be 44%, 59% and 16%.
Previous studies already showed that the mid-IR emission in HD 97048 is extended (Prusti et al. 1994, SI01) on a scale of 5-10 arcsec, and this was attributed to emission in the UIR bands. We are insensitive to diffuse emission on these scales but with our higher spatial resolution we resolve the emission from the central source. Possibly, we have resolved the emission from a large, flared disk, whereas the extended emission seen by SI00 originates in the loose surroundings.
In our interpretation of the SED of HD 97048, namely that it is a
passive, flaring disk, the outer disk needs to have opacity both at
optical/UV wavelengths (to absorb the stellar radiation) and infrared
wavelengths (to emit the observed 10 m continuum radiation). The
spatial extent of the continuum emission implies that the grain
population dominating the continuum is not in thermal equilibrium.
The blackbody temperature at a distance of 50 AU from the central star
is about 100 K, too low for significant emission in the 10
m
region. The continuum must be dominated by very small "superheated''
grains. HD 97048 has no silicate feature and therefore it is excluded
that silicates are the source of the continuum. As the grains
producing the continuum radiation are necessarily small, they would
show a very strong emission band, peaking around 9.7
m, if
they would consist of silicate material. For a continuum opacity
source small carbonaceous particles are the most likely candidate.
Carbon has a high opacity both in the optical and at 10
m.
We note that metallic iron nano-particles also show no spectral structure
and thus may be responsible for the observed extended
emission. However, the strong UIR bands suggest a prominent population
of small carbon rich particles is likely present. Thus, carbon
appears to dominate the opacity of Herbig star disks if no silicates
are present. In sources where small silicate grains are seen in the
spectrum, carbonaceous material may also be responsible for a
significant fraction of the opacity, along with the silicates.
Also HD 100546 has excess emission around 8.2 m
compared to the PAH spectrum as observed in the ISM. It can be seen
in Fig. 1 that as we trace the spatial extent of the emission
blue-ward of the peak of the 8.6
m feature, after an initial
drop, the FWHM starts to increase again at 8.4
m (in HD 97048
this happens only at 8.1
m), indicating that the excess emission
is spatially extended on a scale of approximately 100 AU. The carrier
of this band must therefore a very small, non-thermal equilibrium
grain. Possibly, it is a processed PAH molecule, larger than
a typical ISM PAH. Our data suggest that the 8.2
m
excess emission originates predominantly from the central source
rather than the surrounding nebula (see also Sect. 5.2).
The reason we observe the excess emission at 8.2 m
to be much more extended than
the continuum in HD 100546, whereas in HD 97048 this is not obviously
so, is likely a matter of contrast. In HD 100546, the continuum is
dominated by silicates, which do not emit much at 8.2
m (and emit
significantly in the 10
m window only within about 15 AU from the
star). The 8.2
m emission in this star is fully dominated by the PAHs.
In HD 97048 the continuum emission at 8.2
m dominates the total
flux. It arises from small carbonaceous grains and the continuum
emission has a FWHM of about 100 AU, equal within errors to the
FWHM of the 8.2
m excess emission in HD 100546. It is therefore
not surprising that we observe the 8.2
m excess emission in HD 97048
not to be significantly more extended than the continuum, contrary to
HD 100546.
![]() |
(1) |
![]() |
(2) |
![]() |
(3) |
The temperature of a PAH molecule, after it has just absorbed a UV photon,
is given by (Schutte et al. 1992):
![]() |
(4) |
Table 2:
The the photon absorption rate
for a typical ISM PAH molecule (at the distance from the star where the PAH emission is observed), and the minimum number of C-atoms per molecule needed for multi-photon processes to be significant (
).
where
denotes the photon energy in eV.
The observed UV flux of HD 97048 is
erg s-1 cm-2(van Kerckhoven et al. 2002),
with an assumed distance of 180 pc
(van den Ancker et al. 1998),
this yields a UV luminosity of
erg s-1, which
corresponds to 6.54
.
The bolometric luminosities of
HD 97048 and HD 100546 are 40.7 and 36.0
,
respectively
(van den Ancker et al. 1998).
Since the spectral types of both stars (A0 and B9)
are almost identical, we assume that the spectral shape of HD 100546 equals
that of HD 97048, yielding a UV luminosity of
for HD 100546.
The photon absorbtion rates
for PAH
molecules with 50 C-atoms (typical, small ISM PAH), and the minimum number of
C-atoms needed to absorb 1 photon per second
(the cooling
time of a PAH molecule is about 1 s
(Bakes et al. 2001),
so multi-photon processes will
become significant for particles containing more than
carbon
atoms), are given in Table 2. The median energy for a
UV photon in these stars is 7.5 eV
(van Kerckhoven et al. 2002).
In both stars, multi-photon processes will become important only for PAH particles containing more than about 1000 C-atoms, at the distances
where the PAH emission is observed (note that extinction, which
is likely present in the lines of sight from the star to the outer
disk, could significantly increase this number). The temperature that such a PAH molecule attains upon absorbtion of a single 7.5 eV photon is about
280 K, which is just warm enough to emit significantly at 8 m.
Taking into account multi-photon processes, this temperature will
increase somewhat. Therefore, the carrier of the 8.2
m excess emission
could either be a relatively large PAH, in which case multi-photon
processes are important, or a smaller PAH consisting of less than about
500 C-atoms which is heated to some 400 K upon absorption of a single
UV photon. In the former case, the 8.2
m excess emission must
be confined to regions close (<1
)
to the central star, in the
latter, the carrier can be excited out to larger distances.
HD 97048 provides the opportunity to constrain the nature of
the carrier of the 8.2 m excess emission, as it is surrounded by both a
circumstellar disk and a remnant cloud. If the 8.2
m excess emission
is seen both in the central source (disk) and the surrounding cloud,
the carrier must be relatively small, and is excited by single-photon
processes. If the excess emission is seen in the central source only,
the carrier could either be large (requiring multi-photon excitation)
and/or be synthesized locally in the disk. Unfortunately there is no
straightforward way to distinguish between the carrier of the excess
emission being present both in the cloud and the disk but seen only
in the disk because of the required multi-photon excitation, and the
carrier being present in the disk only.
HD 97048 is surrounded by a reflection nebula, which has been shown to be a source PAH emission (SI01). In the previous sections we discussed the properties of the central source, which we associate with the star and its circumstellar disk. In this section we compare the disk spectrum, which we will call the "on source spectrum'', with the ISO-SWS spectrum. The ISO spectrum contains both the disk emission and emission from the diffuse cloud that surrounds the system.
Comparing the two spectra (top panel of Fig. 2), we see that the peak/continuum ratio in the PAH features is different between the on-source spectrum and the ISO SWS spectrum. The ISO spectrum shows stronger PAH emission bands than the on-source spectrum. Note that we multiplied the on-source spectrum by a factor 1.08, such that the continuum level in the ISO and TIMMI 2 spectra match, faciliating comparison of the PAH peak/continuum ratios. Since the absolute photometric accuracy of both ISO-SWS and ground based spectra is about 15%, the continuum level in the two spectra is equal within errors. Moreover, the slope of the continuum emission is indistinguishable between the on-source and ISO spectrum. This again suggests that the disk is the only source of continuum radiation in the nebula. Whereas the central source shows both continuum and PAH feature emission, emission from surrounding nebula seems to be restricted to the PAH bands. Very low level continuum emission in the nebula cannot be excluded on the basis of our data, however.
An interesting possibility for the study of objects like HD 97048, is the extraction of multiple spectra along the slit. In this way, one can study the properties of the PAHs as a function of position, and possibly detect differences between the nebular PAHs and the PAH population in the disk. Unfortunately, our data are not of sufficient quality to do this, but we emphasize the potential of such studies (see also Sect. 5.2).
To account for the
observed FWHM of the central component, as shown in Fig. 2, the
physical size of the emitting region responsible for the continuum
radiation cannot be significantly smaller than 100 AU. A population
of very small, non-thermal equilibrium grains (VSG) must be the source
of the continuum. Note that these grains cannot be silicates, since
HD 97048 neither shows the strong, broad 9.7 m emission band from
small, amorphous silicates, nor the emission bands associated with
crystalline silicates.
Figures 1 and 2 show that both group I stars are extended at 8.6,
11.3 and (HD 97048) 12.7 m, corresponding to the wavelengths of
the UIR emission bands. The red wing of the 7.7-7.9
m UIR complex
is observable from the ground and is spatially extended in our spectra
as well. Peeters et al. (2002) show that the group of isolated
HAEBE stars show similar behavior of the position and relative
strength of the UIR bands as detected by ISO. for HD 97048 however
they note that the ratio of 7.6 to 7.8
m band strength is
different from that of the other HAEBE stars. HD 97048 also differs
from HD 100546 in that the former has an extended optical reflection
nebula while the latter has not.
Recently, an extensive study of interstellar and circumstellar PAHs
has revealed the presence of systematic variations in the profiles and
peak positions of the main bands from source to source
(Peeters et al. 2002).
Specifically, the 6.2 and 7.7 m bands are distinctly
different for sources which have recently synthesized their PAHs in
their ejecta (e.g., post AGB objects and planetary nebula) from those
which are illuminating general interstellar medium materials. Two of
the sources in this study - HD 97048 and HD 100546 - were also
included in the Peeters et al. sample. The spectral characteristics of
the UIR bands in these two sources are very different from those in
sources with pure ISM materials. In particular, the 6.2 and 7.7
m
bands in the spectrum of HD 97048 are intermediate between those of of
pure interstellar and pure circumstellar sources. In fact, a linear
combination of these types of spectra provides a good fit to the
observations of HD 97048 (van Kerckhoven et al. 2002). It is tempting to
localize the PAHs with circumstellar spectral characteristics in the
flaring disk while the ISM-PAHs might reside in the envelope of this
source. Our data do not allow us to critically test this hypothesis, but
our method may do so (see Sect. 5.3).
The 6.2 and 7.7
m bands in HD 100546 are very similar to
those of circumstellar materials
(Peeters et al. 2002) and presumably
originate in the disk. In this view, while all the material associated
with these two sources derives originally from interstellar material,
the family of PAHs that made it into the disks - and that dominates
the emission spectra of these two sources - has been strongly
modified. Presently, there is no direct evidence to link this
modification to the disk environment - in principle the modification
might have occurred during the collapse phase - but we note that HD 97300
- the neighbor of HD 97048 - shows a very similar UIR spectrum
while the PAH emission in this source is known to be extended on a
scale of several 1000 AU (van Kerckhoven et al. 2002, SI01).
The processes driving the modification of PAHs around young stellar objects are not known. They might range from changes in the charge state of the emitting PAH family driven by variations in the physical conditions (e.g., illuminating FUV field or electron density), variations in the size spectrum of the PAHs driven by coagulation, chemical changes driven by the exposure to the stellar radiation field, or even a complete chemical re-formation of PAHs in these disk environment.
The continuum emission in HD 100546 arises from a much smaller region
than the PAH emission. The warm silicates are confined to the innermost
regions of the disk. The continuum emission in HD 97048 is spatially
extended on a scale of 100 AU. Together with the absence of
the spectral signature of small silicate grains, this indicates that
carbon rich very small grains are responsible for the continuum emission
in the disk. Small carbonaceous particles are an important source of
opacity in circumstellar disks and may dominate the opacity when no
small silicate grains are present.
Our data suggest that in HD 97048, the excess emission seen around
8.2 m, indicative of processed carbon rich material, is confined
to the disk. This is possibly consistent with an ISM-type PAH
contribution at scales of
1000 AU, and a modified PAH
population at
100 AU. We argue that spatially resolved
longs-lit spectroscopy using large telescopes on good 10
m sites
has strong potential for the study of (differential) dust evolution
in the outer disks of HAe stars, and - if present - their
surrounding cloud remnant.
Acknowledgements
We would like to thank the TIMMI 2 team for excellent assistance during the observations. Miska Le Louarn and Marc Sarazin are acknowledged for clarifying discussions on atmospheric turbulence.