A&A 417, 201-215 (2004)
DOI: 10.1051/0004-6361:20031769
G. Gauba - M. Parthasarathy
Indian Institute of Astrophysics, Koramangala, Bangalore 560034, India
Received 24 June 2003 / Accepted 3 December 2003
Abstract
Parameters of the circumstellar
dust shells of 15 hot post-AGB stars have been derived using a radiative transfer code (DUSTY). Combining the
optical, near and far-infrared (ISO, IRAS) data of the stars, we have
reconstructed their spectral energy distributions (SEDs) and estimated
the dust temperatures, mass loss rates, angular radii of the inner
boundary of the dust envelopes and the distances to these stars.
The mass loss rates (10
-6-10-5 yr-1)
are intermediate between stars at the tip of the AGB and the PN phase.
We have also studied the ISO spectra of 7 of these stars. Amorphous and
crystalline silicate features were observed in IRAS 14331-6435 (Hen3-1013),
IRAS 18062+2410 (SAO 85766) and IRAS 22023+5249 (LSIII +5224) indicating
oxygen-rich circumstellar dust shells. The presence of unidentified
infrared (UIR) band at 7.7
,
SiC emission at 11.5
and the "26
"
and "main 30
" features in the ISO spectrum of IRAS 17311-4924 (Hen3-1428)
suggest that the central star may be carbon-rich. The ISO spectrum of
IRAS 17423-1755 (Hen3-1475) shows a broad absorption feature at 3.1
due
to C2H2 and/or HCN which is usually detected in the circumstellar
shells of carbon-rich stars.
Key words: stars: AGB and post-AGB - stars: early-type - stars: evolution - stars: circumstellar matter - infrared: stars
In the evolution of low and intermediate mass stars (0.8-8 ),
the post-asymptotic giant branch (post-AGB) or protoplanetary nebula (PPN) phase
is a transition stage from the tip of the AGB to the planetary nebula (PN)
stage (Kwok 1993). The hot post-AGB stars form an evolutionary link between
the cooler G, F, A supergiant post-AGB stars (Parthasarathy & Pottasch 1986)
and the hotter O-B central stars of PNe (Parthasarathy 1993a).
Analysis of the UV(IUE) spectra of hot post-AGB stars
(Gauba & Parthasarathy 2003), revealed that in many cases,
the hot (OB) central stars of PPNe are partially obscured by
circumstellar dust shells. Stars on the AGB and beyond are characterised
by severe mass loss (10-8-10-3 yr-1) which
results in the formation of circumstellar envelopes. The physical mechanisms
responsible for the intensive mass loss from AGB stars are not well understood
although the most promising mechanism to date involves radiation pressure on the
dust grains (Tielens 1983). While AGB stars appear to have spherically
symmetric dust outflows (e.g. Habing & Blommaert 1993), PN tend to
have axially symmetric inner regions and spherical outer halos
(e.g. Schwarz et al. 1992). Inorder to understand the mass loss mechanisms,
wind velocities and time scales responsible for the evolution of PNe,
we need to study the circumstellar environment of the stage intermediate between
the AGB and the PN phase, i.e. the post-AGB/PPN phase. Circumstellar dust
shells of some cooler post-AGB stars (e.g. Hoogzaad et al. 2002; Hony
et al. 2003) and PNe (e.g. Siebenmorgen et al. 1994) have been modelled
to derive the dust composition, mass loss rates and dynamical ages.
Table 1: Hot post-AGB stars.
Table 2: Log of ISO observations.
As a consequence of dredge-up of byproducts of helium burning to the
surface of stars on the AGB, the oxygen-rich atmospheres of some of these stars may be
transformed into carbon-rich atmospheres (see e.g. Iben & Renzini 1983).
This change of chemistry would also be reflected in the composition of the dust
grains formed in the cirumstellar envelopes of AGB and post-AGB stars.
With the resolution and wavelength coverage of the ISO mission
(Kessler et al. 1996) the detection of prominent gas and solid state features
specific to oxygen-rich and carbon-rich chemistries became possible.
Amorphous and crystalline silicate features and crystalline water have been
reported in the ISO spectra of some AGB and post-AGB stars and
the nebulae surrounding [WC] central stars of PNe (see e.g. Waters & Molster 1999;
Hoogzaad et al. 2002). Hrivnak et al. (2000) detected the "21 " and "30
"
emission features besides the unidentified infrared (UIR) emission bands
at 3.3, 6.2, 7.7 and 11.3
in the ISO spectra of a sample of carbon-rich PPNe.
Bogdanov (2000, 2002, 2003) modelled the complete spectral energy distribution (SED) of three hot post-AGB stars, IRAS 18062+2410 (SAO 85766), IRAS 19590-1249 (LSIV-12 111) and IRAS 20462+3416 (LSII+34 26) using radiative transfer codes and derived their mass loss rates, inner radii of the dust envelopes, optical depth of the envelopes and the distances to these stars. We need to study a bigger sample of such stars to understand the evolution of the infrared spectrum as the stars evolve from the cooler post-AGB phase to the hot central stars of PNe. In this paper, we have used the radiative transfer code, DUSTY (Ivezic et al. 1999) to model the circumstellar dust shells of 15 hot post-AGB stars. Additionally, 7 stars from our list were found to have ISO spectra. We also report the analysis of the ISO spectra of these stars.
The hot post-AGB stars for this study (Table 1) were selected from the papers of Parthasarathy & Pottasch (1989), Parthasarathy (1993b) and Parthasarathy et al. (2000a). High and low resolution optical spectra have confirmed the post-AGB nature of several of these stars, e.g. IRAS 17119-5926 (Hen3-1357; Parthasarathy et al. 1993c, 1995; Bobrowsky et al. 1998), IRAS 18062+2410 (SAO 85766; Parthasarathy et al. 2000b; Mooney et al. 2002; Ryans et al. 2003), IRAS 19590-1249 (LSIV-12 111; McCausland et al. 1992; Conlon et al. 1993a, 1993b; Mooney et al. 2002; Ryans et al. 2003) and IRAS 20462+3416 (LSII+34 26; Parthasarathy 1993a; García-Lario et al. 1997a). The UV(IUE) spectra of some hot post-AGB stars listed in Table 1 showed violet shifted P-Cygni profiles of CIV and NV indicating stellar wind and post-AGB mass loss (Gauba & Parthasarathy 2003).
The ISO data archive was searched for spectra of the hot post-AGB
stars listed in Table 1. Six of the fifteen sources were found to have
only SWS (Short Wavelength Spectrometer) spectra
while one source, IRAS 17423-1755 (Hen3-1475) had both SWS and LWS
(Long Wavelength Spectrometer) spectra. The log of the observations is
given in Table 2. ISO SWS spectra have a wavelength coverage from
2.38-45.2 .
The LWS spectra extend from 43-197
.
Our objects were observed in the low resolution (AOT 01) mode of the SWS instrument (de Graauw et al. 1996). A spectrum scanned with SWS contains 12 subspectra, that each consist of two scans, one in the direction of decreasing wavelength ("up'' scan) and the other in the direction of increasing wavelength ("down'' scan). There are small regions of overlap in the wavelength between the sub-spectra. Each sub-spectrum is recorded by 12 independent detectors.
The LWS observations were carried out in LWS01 mode, covering the
full spectral range at a resolution (/
)
of
200. The characteristics of the ISO LWS instrument are
described in Clegg et al. (1996) and the calibration of the instrument
is described in Swinyard et al. (1996).
![]() |
Figure 1a:
The ISO SWS spectrum of IRAS 14331-6435 (Hen3-1013)
shows emission due to amorphous (10.8 ![]() ![]() ![]() ![]() |
![]() |
Figure 1b:
The ISO SWS spectrum of IRAS 17311-4924 (Hen3-1428)
shows emission due to the UIR band at 7.7 ![]() ![]() ![]() ![]() ![]() |
![]() |
Figure 1c:
The ISO SWS spectrum of IRAS 22023+5249 (LSIII +5224)
shows emission due to amorphous (10.8 ![]() ![]() |
In this section we describe the analysis of the ISO spectra and the modelling of the SEDs of the hot post-AGB stars.
Offline processed ISO SWS01 (OLP version 10.1) and LWS01 (OLP version 10) data were retrieved from the ISO data archive. These were further processed using ISAP (ISO Spectroscopic Analysis Package) version 2.1.
The data analysis using ISAP consisted of extensive bad data removal
primarily to minimize the effect of cosmic ray hits. All detectors
were compared to identify possible features. For each sub-band, the
flux level of the 12 detectors were shifted and brought to a mean value.
They were then averaged, using median clipping to discard points that
lay more than 2.5
from the median flux. The spectra were
averaged typically to a resolution of 300, 500, 800 and 1500
for SWS01 data taken with speed 1, 2, 3 and 4 respectively (Table 2).
Appropriate scaling factors were applied to the averaged spectra of
each sub-band to form a continous spectrum from 2.38-45.2
.
The data of subband 3E (27.5-29.0
)
are generally noisy
and unreliable (see e.g. Hrivnak et al. 2000; Hony et al. 2002).
Figures 1a-c and 2a show the SWS spectra of the hot post-AGB
stars. The SWS data of IRAS 14331-6435 (Hen3-1013), IRAS 18062+2410
(SAO 85766) and IRAS 22023+5249 (LSIII +5224) below 7
and
that of IRAS 16206-5956 (SAO 243756), IRAS 22495+5134 (LSIII +5142)
below 12.5
were noisy and have been excluded. This is also
evident from the low IRAS 12
flux (<5 Jy) of these stars.
Identification of the infrared spectral features is based on Waters & Molster
(1999), Cox (1993), Hrivnak et al. (2000), Volk et al. (2002),
Cernicharo et al. (1999) and J
rgensen et al. (2000).
![]() |
Figure 2:
ISO SWS (2.4
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 3:
Spectral energy distributions (SEDs) of the hot post-AGB stars.
![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 3: continued. |
Open with DEXTER |
Reduction of the LWS observation of IRAS 17423-1755 (Hen3-1475) consisted of
extensive bad data removal using ISAP rebinning on a fixed resolution grid of
/
.
Appropriate scaling factors were applied to
the data from different LWS detectors to form a continuous spectrum (Fig. 2).
The LWS spectrum of IRAS 17423-1755 (Hen3-1475) appears featureless.
Although the spectral shape is very reliable, the absolute flux
calibration uncertainity is 30% for the SWS at 45
(Schaeidt et al. 1996) and 10-15% for the LWS at
45
(Swinyard et al. 1998). The SWS and LWS spectra
of IRAS 17423-1755 (Hen3-1475) were scaled according to their
fluxes in the overlap region. The difference between the flux
levels of LWS and SWS in the overlap region was smaller than
30% which is acceptable within the limits of the combined
error bars. The combined SWS-LWS spectra of
IRAS 17423-1755 (Hen3-1475) was used in Fig. 3.
To re-construct the spectral energy distributions (SEDs) of the
objects, we combined the ISO data with available U, B, V, R, I, J, H, K, L, L'
and M magnitudes of the stars from literature (Table 3a). We also
searched the 2MASS (2Micron All Sky Survey) Catalog within 15
of each object for their JHK magnitudes and the Midcourse Space Experiment
(MSX) catalog within 3
of each object. The 2MASS data was
included only when it was found to be free from confusion flags. MSX fluxes were
found only for IRAS 17460-3114 (SAO 209306; Table 3b). Infrared data
(8.7, 10, 11.4, 12.6 and 19.5
)
on IRAS 18062+2410
(SAO 85766; Table 3c) was obtained from Lawrence et al. (1990).
The temperatures of the central stars (Table 4) are mainly from Gauba & Parthasarathy (2003). For IRAS 12584-4837 (Hen3-847) and IRAS 17423-1755 (Hen3-1475), we estimated the central star temperatures based on their spectral types. For the PN, IRAS 22495+5134 (LSIII +5142) we used a central star temperature of 35 000 K (Tylenda & Stasinska 1994). Gauba & Parthasarathy (2003) found that the UV(IUE) spectrum of IRAS 22023+5249 (LSIII +5224) closely resembles that of a B2-supergiant. Hence we adopted a temperature of 18 500 K for the star, corresponding to the spectral type B2I. LTE analysis of the high resolution optical spectra of IRAS 18062+2410 (SAO 85766) was carried out by Mooney et al. (2002) and Parthasarathy et al. (2000b). Recently, using non-LTE analysis, Ryans et al. (2003) reported an effective temperature of 20 750 K for IRAS 18062+2410 (SAO 85766) and 20 500 K for IRAS 19590-1249 (LSIV-12 111). In this paper, we have adopted the temperature estimates by Ryans et al. (2003).
The interstellar extinction (E(B-V)
)
in the direction
of the stars were estimated using the Diffuse Infrared Background
Experiment (DIRBE)/IRAS dust maps (Schlegel et al. 1998; Table 4).
The optical spectral types of the stars are mainly from Parthasarathy
et al. (2000a). The intrinsic B-V values, (B-V)
,
for
the optical spectral types of the stars, were taken from
Schmidt-Kaler (1982). We adopted (B-V)
for
IRAS 12584-4837 (Hen3-847) and IRAS 17423-1755 (Hen3-1475)
corresponding to
K. For IRAS 22023+5249 (LSIII +5224), we
used (
corresponding to B2I spectral type
and for the PN, IRAS 22495+5134 (LSIII +5142), we adopted
correspoding to
K.
Using the observed and intrinsic B-V values
we derived the total (interstellar plus circumstellar) extinction,
E(B-V)
(=(B-V)
)
towards these stars. Comparing
and
E(B-V)
,
it is evident that there is considerable
circumstellar extinction in most cases. Gauba & Parthasarathy (2003)
found that the circumstellar extinction law in the UV
(from
1300 Å to 3200 Å) varies linearly as
in the case of IRAS 13266-5551 (CPD-55 5588), IRAS 14331-6435 (Hen3-1013),
IRAS 16206-5956 (SAO 243756), IRAS 17311-4924 (Hen3-1428),
IRAS 18023-3409 (LSS 4634), IRAS 18062+2410 (SAO 85766),
IRAS 18371-3159 (LSE 63), IRAS 22023+5249 (LSIII +5224) and
IRAS 22495+5134 (LSIII +5142).
Since little is known about circumstellar extinction laws,
we corrected the observed optical and near infrared magnitudes
of the stars for the total extinction (E(B-V)
)
using the standard extinction laws by Rieke & Lebofsky (1985).
In particular, Rieke & Lebofsky (1985) assume Rv=3.1,
where,
Rv=A v/E(B -V). Although this is true for
interstellar extinction, it may not be strictly true
for circumstellar extinction obeying a
law
in the UV. In particular, Rv may be different from 3.1 in the case of circumstellar extinction.
The use of the radiative transfer code, DUSTY (Ivezic et al. 1999)
for modelling the circumstellar dust shells of hot post-AGB stars was
described in Gauba et al. (2003). DUSTY uses six different grain
types: "warm'' (Sil-Ow) and "cold'' (Sil-Oc) silicates from Ossenkopf
et al. (1992), silicates and graphites (Sil-Dl and grf-DL) from
Draine & Lee (1984), amorphous carbon (amC-Hn) from Hanner (1988)
and SiC (SiC-Pg) from Pégourié (1988). The central stars were
assumed to be point sources at the centers of the spherical density
distributions. The SEDs of the central stars were assumed to be
Planckian. The standard Mathis, Rumpl, Nordsieck (MRN)
(Mathis et al. 1977) power-law was used for the grain size (n(a))
distributions, i.e.
for
with q=3.5,
and
.
For each object, the
dust temperature (
)
on the inner shell boundary and the optical
depth (
)
at 0.55
(V-band) were varied. We assumed an inverse
square law (y-2) for the spherical density distribution. The shell
was assumed to extend to 1000 times its inner radius. We adopted the
fits for which the sum of squares of the deviations between the observed
and modelled fluxes (after scaling) were a minimum. Table 4 lists
the adopted input parameters. Figure 3 shows the spectral energy
distribution of the stars. DUSTY does not allow simultaneous modelling
of warm and cold dust shells. Hence, the cold dust in the case of
IRAS 12584-4837 (Hen3-847) and IRAS 17423-1755 (Hen3-1475) had to be
modelled and treated independent of the warm dust around these stars.
Having fixed
and
,
we then used the gas-dynamical mode
of the DUSTY code to derive the inner radii, r1(cm) where the dust temperatures
(
)
are specified and the mass-loss rates (
).
The radius scales in proportion to L1/2 where L is the luminosity
and the code output value corresponds to L=104
.
The
mass-loss rate scales in proportion to L3/4(
where,
the gas-to-dust mass ratio,
and the dust grain density,
g cm-3. The hot post-AGB stars discussed in this paper
have a range of core-masses (Gauba & Parthasarathy 2003).
Pottasch (1992) pointed out that the white
dwarf distribution is sharply peaked with a mean mass between 0.56 and
0.58
and central stars of PNe have core-masses which show a
peak at approximately 0.6
.
We carried out calculations for
our hot post-AGB stars with core masses of 0.565
and 0.605
corresponding to luminosities of 4500
(Schönberner 1983)
and 6300
(Blöcker 1995) respectively. Distances (d)
to the stars were derived using r1 and the ratio of the observed and
modelled fluxes at 0.55
.
(=r1/d) is the angular radii of
the inner boundary of the cold circumstellar dust envelopes.
Tables 4, 5a and b list the respective values for ,
r1, d,
and
.
All calculations were
carried out using the best fit parameters for the cold circumstellar
dust shells.
The hot post-AGB stars in this paper, except IRAS 19590-1249 (LSIV-12 111), have been described in Gauba et al. (2003) and Gauba & Parthasarathy (2003). Here, we describe, the ISO spectra and dust shell characteristics of some of these objects.
#IRAS 12584-4837 (Hen3-847)
It was found to be variable in the optical (Kazarovets et al. 2000;
de Winter et al. 2001).The Hipparcos magnitudes at maximum and
minimum are 10
52 and 10
70 respectively. Comparison of the
J, H, K magnitudes of the star from Fouque et al. (1992) and the
2MASS catalog indicate infrared variability as well. The SED of this star
shows the presence of both warm and cold circumstellar dust. We have
modelled the warm dust based on the data of Fouque et al. (1992) (775 K) and
the 2MASS catalog (700 K). The presence of warm circumstellar dust may
indicate ongoing post-AGB mass loss. The use of amorphous carbon
grains to model the warm dust and silicate grains for the cold dust,
indicates that during its evolution along the AGB, the central star
may have evolved from an oxygen-rich to a carbon-rich star.
#IRAS 14331-6435 (Hen3-1013)
The ISO SWS spectrum reveals the presence of amorphous
(10.8 )
and crystalline silicates and/or water (33.6
,
40.4
,
43.1
)
in the circumstellar environment of this star. Crystalline
silicates have been detected in the dust shells around evolved
oxygen-rich stars (see e.g. Waters et al. 1996; Waters & Molster 1999).
Waters et al. (1996) found that these emission features are more prominent
for objects with cooler dust shells (T < 300 K). Emission from
crystalline water has been reported in the 40-70
spectrum of the
Frosty Leo nebula and other cool oxygen-rich envelopes (Omont et al. 1990).
#IRAS 16206-5956 (SAO 243756)
The ISO spectrum of this star is noisy and no features could be identified.
The continuum from 12.5
to 45.2
was used in addition to the
IRAS fluxes to better constrain the model fit to the SED.
#IRAS 17311-4924 (Hen3-1428)
The broad "30
emission feature" was detected in several AGB stars and PN sources
(Forrest et al. 1981; Cox 1993). More recently it was also detected in
carbon-rich PPNs possessing the 21
emission feature (Omont et al. 1995).
Hony et al. (2002) detected this feature in the ISO SWS spectrum of
IRAS 17311-4924 (Hen3-1428) and a large sample of carbon-rich AGB stars
(C-stars), post-AGB stars and PNe. Substructure in the 30
feature was
recognised by Szczerba et al. (1999). In the spectra of several PPNs
and carbon stars, it was resolved into two components
(Hrivnak et al. 2000; Volk et al. 2000, 2002), the "26
"
and "main 30
" feature (Fig. 1b). Goebel & Moseley (1985) first
suggested that the feature is due to magnesium sulfide (MgS).
Hony et al. (2002) too identified MgS as the carrier of the "26
"
and "main 30
" emission features. However, the 30
band is never
seen in oxygen-rich sources (Forrest et al. 1979). Since the feature is seen
only in carbon-rich objects, the suggestion that its carrier is a carbonaceous material
continues to be appealing (Volk et al. 2002). We modelled the SED of
IRAS 17311-4924 (Hen3-1428) using graphite and silicon carbide (SiC).
However, we could not obtain a fit in the 30
emission region
(from
20-30
).
The 11.5
band (Fig. 1b) is attributed to SiC (Treffers & Cohen 1984)
and has been detected in carbon-rich evolved objects (see e.g., Cernicharo et al. 1989). Cox (1993) pointed out that all AGB stars where a 30
emission band is present
show the 11.5
band in emission with the notable exception of GL 3068 where
it is seen in absorption. However, the reverse is not true and some caron-rich
sources with a a prominent 11.5
emission band do not show the 30
band (V Cyg, S Cep and Y CVn). Using the DUSTY model, we
obtained a good fit to the SiC emission.
Table 3a: Photometric data on hot post-AGB stars.
Table 3b: MSX data.
Table 3c: Infrared data on IRAS 18062+2410 (SAO 85766).
We also detected the 7.7
UIR feature in the ISO
spectrum of this object. UIR features are commonly attributed to
polyaromatic hydrocarbons (PAHs) and have been detected in other
carbon-rich protoplanetary nebulae (see e.g. Beintema et al. 1996; Hrivnak et al. 2000).
#IRAS 17423-1755 (Hen3-1475)
Gauba et al. (2003) modelled the SED of this star using DUSTY. They found a warm dust
component at 1500 K in addition to the cold-dust at 100 K indicating
ongoing post-AGB mass-loss. Plotting the ISO spectrum of the star (2.38 -171
)
alongwith the photometric data from Gauba et al. (2003), near-IR data
from García-Lario et al. (1997b) and IRAS fluxes and using the DUSTY code, we
detected a second warm dust component at 1000 K.
The broad absorption feature at 3.1
seen in the ISO spectrum is due
to C2H2 and/or HCN (Ridgway et al. 1978; Cernicharo et al. 1999;
J
rgensen et al. 2000). This feature is observed in carbon stars
(Merrill & Stein 1976; Noguchi et al. 1977; Groenewegen et al. 1994).
The far-IR flux distribution of the star was modelled with amorphous carbon and
silicon carbide confirming the carbon rich nature of the circumstellar
dust shell. The LWS spectrum appears featureless.
#IRAS 18062+2410 (SAO 85766)
High resolution optical spectra (Parthasarathy et al. 2000b;
Mooney et al. 2002; Ryans et al. 2003) indicated the underabundance
of carbon in this star. The ISO spectrum shows strong emission due
to amorphous silicates at 10.8
and 17.6
in conformity
with an oxygen-rich chemistry (i.e.
)
for the central
star.
The J, H, K magnitudes from García-Lario et al. (1997b) and the 2MASS catalog indicate that the star is variable in the near-infrared. Optical variations of the star were detected by Arkhipova et al. (1999). The K, L, M fluxes by Lawrence et al. (1990) lie above the modelled SED (Fig. 3). In particular the K-band flux of Lawrence et al. (1990) lies above the K-band flux estimated by García-Lario et al. (1997b). If the K, L, M fluxes of the star by Lawrence et al. (1990) were not overestimated, the observed mismatch may be due to the variable nature of this object. Gauba & Parthasarathy (2003) reported variable circumstellar extinction which in addition to stellar pulsations may be due to a dusty torus in motion around the hot central star.
Bogdanov (2000) modelled the SED of this star using a different radiative transfer
code. They derived
K. This value is much higher than the dust temperature
(
K) derived by us from the model fit. Bogdanov (2000) had used the K, L, M fluxes by Lawrence et al. (1990) only. Besides they did not use the ISO
spectrum of the star to constrain the SED. We believe that we have a better estimate
of the physical parameters of the star especially since our model
gives a very good fit to the IRAS fluxes and ISO spectrum of the star.
Table 4: Input physical parameters for DUSTY and the adopted reddening values.
#IRAS 18379-1707 (LSS 5112)
The star was found to be variable in the near-infrared based on the J, H, K magnitudes from García-Lario et al. (1997b) and the 2MASS catalog.
#IRAS 19590-1249 (LSIV-12 111)
It was classified from low dispersion spectroscopy by Kilkenny & Pauls (1990)
as having a spectral type around B0. Based on an analysis of its high
resolution optical spectra, McCausland et al. (1992) and Conlon et al.
(1993a,b) concluded that its chemical composition and atmospheric
parameters are consistent with a post-AGB evolutionary state. The
presence of nebular emission lines of [NII], [OII] and [SII] in its
optical spectrum suggest that it may be in the early stages of PN
formation. Study of the atmospheric parameters and abundance analysis
of the star have also been carried out
by Mooney et al. (2002) and more recently a non-LTE analysis by Ryans
et al. (2003). Based on these studies it has been confirmed that the
central star of LSIV-12 111 shows severe carbon deficiency.
Using silicate grain composition for the dust envelope, we obtained a
good fit to the SED (Fig. 3). We estimated a dust temperature of 120 K
which is in good agreement with the estimate of Conlon et al. (1993b)
based on the 100
IRAS flux and the color-color diagram of
Jourdain de Muizon et al. (1990). Recently, Bogdanov (2003) modelled
the spectral energy distribution of this star using the DUSTY code.
Although they adopted different values for the extinction and the
effective temperature of the central star, we find their derived
values for the dust temperature, distance and mass loss rate
to be in good agreement with our results.
#IRAS 22023+5249 (LSIII +5224)
Amorphous (10.8 )
and crystalline (33.6
)
silicate features
were detected in the ISO spectrum of the dust shell surrounding this
star indicating the oxygen-rich nature of the central star.
#IRAS 22495+5134 (LSIII +5142)
The ISO spectrum of this PN is noisy and spectral features could not be
identified. The angular radii ()
derived using the DUSTY code
(Tables 5a and b) are in good agreement with the angular diameter of
4
reported by Tylenda & Stasinska (1994) and Acker et al. (1992)
respectively. From the Hipparcos and Tycho Catalogues (ESA 1997), the PN
has a parallax of 0.42 mas. This implies a distance of 2.38 kpc.
Assuming a core mass (
)
of 0.605
,
we derived a
distance of 2.17 kpc to the PN.
We have modelled the circumstellar dust shells of 15 hot post-AGB stars using the radiative transfer code, DUSTY and derived their dust temperatures, distances to the stars, mass loss rates and angular radii of the inner boundary of the dust envelopes (Tables 5a and b). These stars have detached dust shells (as is evident from the SEDs, Fig. 3), OB-giant or supergiant spectra and cold dust between 100-315 K, satisfying the observational properties of PPNe as defined by Kwok (1993, 2001). In addition to the cold dust, warm dust was detected in the case of IRAS 12584-4837 (Hen3-847) and IRAS 17423-1755 (Hen3-1475) indicating ongoing mass loss. From the grain types used for the model fits, we may infer the chemical composition of the circumstellar dust shells. The use of both silicate and amorphous carbon grains to model the SEDs of IRAS 12584-4837 (Hen3-847) and IRAS 17460-3114 (SAO 209306) suggests that the central stars in these two cases may have undergone a recent change from an oxygen-rich to a carbon-rich chemistry. Such hot post-AGB stars may evolve into the [WC] central stars of PNe. Recently, Waters et al. (1998) detected carbon-rich PAH features in the near-infrared and crystalline silicates in the far-infrared ISO spectra of two PNe with [WC] central stars, BD+30 3639 and He2-113.
Table 5a:
Derived stellar and dust envelope parameters for
.
Table 5b:
Derived stellar and dust envelope parameters for
.
Observational evidence (e.g. Chu et al. 1991) suggests that three winds
are involved in stripping the outer envelope of the AGB star on its
way to becoming a PN (Marten et al. 1993; Frank 1994): the
spherically symmetric AGB wind (e.g. Habing & Blommaert 1993) when the
star loses mass at rates of
10-7-10-6 yr-1 with
a wind velocity of
10 km s-1; the superwind phase when the
mass loss is thought to increase dramatically at the end of the AGB,
upto
10-5-10-3
yr-1, still with a wind velocity of
10 km s-1; once the superwind exhausts most of the AGB star's envelope,
a fast wind with mass loss rate of 10-8
yr-1 and
velocity of
1000 km s-1 develops at some point during the PPN phase.
Velocities of 1000 km s-1 and mass loss rates of
10-8
yr-1have been observed in the central stars of PNe (e.g. Gauba et al. 2001).
For our hot post-AGB stars, we derived mass loss rates of 10
-5-10-6
yr-1.
The mass loss rates (
)
scale with the
gas-to-dust mass ratio (
). We have adopted
.
For
carbon-rich AGB and post-AGB stars values between 200 and 250 are often used
(e.g. Jura 1986; Meixner et al. 1997). For the cool (F3Ib) post-AGB star,
HD 161796 (Parthasarathy & Pottasch 1986) with an oxygen-rich
circumstellar environment, Hoogzaad et al. (2002)
estimated
.
Furthermore, our models assume that the dust density
distribution falls off as y-2 in the entire circumstellar dust shell. Such an
assumption would break down in the case of episodic mass loss (Olofsson et al. 1990).
e.g. episodic mass loss may have been responsible for the rapid evolution (30-40 years) of
IRAS 17119-5926 (Hen3-1357) and IRAS 18062+2410 (SAO 85766) from B-type
post-AGB supergiants to young PNe (Parthasarathy et al. 1993c, 1995;
Bobrowsky et al. 1998; Parthasarathy et al. 2000b).
Table 6: Predicted and observed (V-J) colors.
The proper motions ()
of the stars from the
Tycho-2 Catalogue (Hog et al. 2000) have been listed in Tables 5a and b.
Using the derived distances (d) in conjunction with the proper motions
we estimated the component of the stellar space velocities of the targets
tangent to the line of sight (
).
For IRAS 17203-1534, IRAS 18062+2410 (SAO 85766) and IRAS 18371-3159
(LSE63), the large
values (Tables 5a and b) imply very high space
velocities (
;
where
is the radial velocity of a star), close to the
escape velocity from the Galaxy of 290 km s-1 near the Sun.
Mooney et al. (2002) estimated a distance of 8.1 kpc to IRAS 18062+2410
(SAO 85766) which is much greater than our estimate of
5 kpc. Such a large distance, if correct, would imply a still higher
space velocity. We believe our distance estimates to be closer to
the actual values for these stars. However, the assumption
of spherical density distributions in our models, may be an over
simplification for some of these objects.
e.g. IRAS 17423-1755 (Hen3-1475) has IRAS colors moderately close to
those of HD 233517. HD 233517 is unresolved and may have a disk
instead of a spherical outflow (see e.g. Jura 2003; Fisher et al.
2003). The predicted angular sizes of the inner radii of the dust shells
(Tables 5a and b) suggests that these objects should be easily
resolvable in the mid-IR images with large ground based telescopes.
Imaging of these objects in the IR would serve to test the basic assumptions
such as those of spherical symmetry for our models.
In Table 6 we have compared the predicted and observed (V-J) values
for these stars (
)
where,
.
The intrinsic (V-J) colors ((V-J)
), for the spectral types of
the stars are from Ducati et al. (2001). For stars with emission lines in their
optical spectra, IRAS 12584-4837 (Hen3-847), IRAS 17423-1755 (Hen3-1475),
IRAS 22023+5249 (LSIII+5224) and for the PN, IRAS 22495+5134 (LSIII+5142),
(
could not be assumed. In the case of
IRAS 14331-6435 (Hen3-1013), IRAS 17203-1534,
IRAS 18062+2410 (SAO 85766) and IRAS 18379-1707 (LSS 5112), we find
significant differences between the values of
and
.
On first sight, this would then raise a suspicion
about the adopted
values. However, we would like to point out
that the V and J magnitudes of these stars have not been recorded
simultaneously. Many of these stars are variable as evidenced from
the J, H, K data on IRAS 12584-4837 (Hen3-847), IRAS 18062+2410 (SAO 85766)
and IRAS 18379-1707 (LSS 5112). The V and J magnitudes may have be recorded
at different epochs of the variability cycle and hence it may
not be suitable to compare (
and (
.
We also studied the ISO spectra of 7 hot post-AGB stars,
IRAS 14331-6435 (Hen3-1013), IRAS 16206-5956 (SAO 243756),
IRAS 17311-4924 (Hen3-1428), IRAS 17423-1755 (Hen3-1475),
IRAS 18062+2410 (SAO 85766), IRAS 22023+5249 (LSIII +5224) and IRAS 22495+5134
(LSIII +5142). A weak amorphous silicate feature (10.8 )
alongwith
crystalline silicate features was found in the dust shells of
IRAS 14331-6435 (Hen3-1013) and IRAS 22023+5249 (LSIII +5224).
The 17.6
amorphous silicate feature was missing in these two stars.
The post-AGB star IRAS 18062+2410 (SAO 85766) did not show evidence
for the presence of crystalline silicates but strong amorphous silicate
features at 10.8
and 17.6
were detected. Volk & Kwok (1989)
predict that at dust temperatures of typically a few 100 K for
post-AGB stars, the spectrum should increase from 8
to 23
and the 10.8
and 17.6
silicate features should cease
to be observable. This appears to be consistent with the observed
spectral features and the dust temperatures of 230 K, 130 K and 120 K
for IRAS 18062+2410 (SAO 85766), IRAS 14331-6435 (Hen3-1013) and
IRAS 22023+5249 (LSIII +5224) respectively from our model fits.
The presence of silicate features in these stars indicates the
O-rich nature of the central stars. The formation of crystalline
silicates in the circumstellar shells of post-AGB stars
is still not well understood (see e.g., Waters et al. 1996).
In contrast, PAH emission at 7.7
,
the "26
" and "main 30
" features
and 11.5
SiC emission in IRAS 17311-4924 (Hen3-1428), are typical of
circumstellar dust shells around carbon-rich post-AGB stars. However, the
21
emission feature detected in several carbon-rich PPNe (Hrivnak et al. 2000) was
notably absent in the ISO spectrum of IRAS 17311-4924 (Hen3-1428).
Volk et al. (2002) pointed out that although all sources
with the 21
emission feature also display the
"26
" and "main 30
" features, the converse is not true.
The hot post-AGB star, IRAS 01005+7910 (Klochkova et al. 2002) which showed the "26"
and "main 30
" emission also did not show the 21
emission feature
(Hrivnak et al. 2000). It may be that the dust grains responsible for
the 21
emission are destroyed as the central star evolves
towards hotter temperatures. The broad absorption feature at 3.1
in
IRAS 17423-1755 (Hen3-1475) attributed to C2H2 and/or HCN indicates
that the central star may be carbon-rich.
Acknowledgements
We are grateful to the referee for helpful comments. G.G. would also like to thank Luciano Rezzolla and John Miller of the Astrophysics sector at SISSA, Trieste, Italy for providing computer facilities during her stay in Italy. The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS Instrument Teams and Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, RAL and SRON.