A&A 441, 289-302 (2005)
DOI: 10.1051/0004-6361:20047013
M. Kraus1 - M. Borges Fernandes1,2 - F. X. de Araújo2 - H. J. G. L. M. Lamers1
1 - Astronomical Institute, Utrecht University,
Princetonplein 5, NL 3584 CC Utrecht, The Netherlands
2 -
Observatório Nacional-MCT, Rua General José Cristino 77,
20921-400 São Cristovão, Rio de Janeiro, Brazil
Received 5 January 2004 / Accepted 29 April 2005
Abstract
We study the optical spectrum of the exciting B[e] star Hen 2-90
based on new high-resolution observations that cover the innermost 2
of this object whose total extent is more than 3
.
Our investigation is split in two parts, (i) a
qualitative study of the presence of the numerous emission lines and
classification of their line profiles, which indicate a circumstellar environment of
high complexity, (ii) and a quantitative analysis of numerous forbidden lines, e.g.
[O I], [O II], [O III], [S II], [S III], [Ar III],
[Cl II], [Cl III], and [N II].
We find correlation between the different ionization states of the elements and
the velocities derived from the line profiles: the highly ionized atoms have the
highest outflow velocity, while the neutral lines have the lowest.
The recent HST image of Hen 2-90 (Sahai et al. 2002, ApJ, 573, L123) reveals a bipolar,
highly ionized region, a neutral disk-like structure, and an intermediate region of
moderate ionization. This HST image covers about the same innermost regions as our
observations. When combining the velocity information with the HST image of Hen
2-90, it seems that a non-spherical stellar wind model is a good option to
explain the ionization and spatial distribution of the circumstellar material.
Such a wind might expand into the cavity formed during the AGB phase of the
star, which is still visible as a large nebula, seen
e.g. on H
plates. We modelled the forbidden lines under the assumption of
a non-spherically symmetric wind that can be split into a polar, a disk forming,
and an intermediate wind, based on the HST image. We find that in order to fit the
observed line luminosities, the mass flux, surface temperature, and terminal
wind velocities need to be latitude dependent, which might be explained in terms
of a rapidly rotating central star. A rotation speed of 75-80% of the critical velocity
was derived from the terminal velocities extracted from the observed line wings
considering the inclination of the system as suggested from the HST image.
The total mass loss rate of the star was determined to
be on the order of
yr-1.
The combination of this wind scenario and the underabundance of
C, O, and N in comparison to the solar abundance of
S, Ar, and Cl might be explained in terms of a
rapidly rotating post-AGB star.
Key words: stars: mass-loss - ISM: planetary nebulae: individual: Hen 2-90 - line: identification - methods: data analysis
Hen 2-90 (PN Sa 2-90, PN G305.1+01.4, IRAS 13064-6103) is a very interesting object
whose evolutionary stage is unclear. Its distance was estimated to be about
d = 1.5 kpc by Costa et al. (1993), while Sahai & Nyman
(2000) adopted d = 2.5 kpc. Its effective temperature was estimated
to be 51 000 K (Kaler & Jacoby 1991), and its luminosity was derived
as
(Costa et al. 1993).
The star was first classified as a planetary nebula by Webster (1966)
and Henize (1967). Later, Costa et al. (1993) and Maciel
(1993) classified it as a planetary nebula with low metal abundance and
with a central star of low mass and low luminosity. Lamers et al. (1998)
included it in the list of objects presenting the B[e] phenomenon, as a compact
planetary nebula B[e].
Sahai & Nyman (2000), based on HST images, and Guerrero et al. (2001), using ground-based images, have described the presence of a nebula bisected by a disk and with both a bipolar jet and knots, spaced uniformly. Guerrero et al. (2001) note that the dynamical stability of the jets and knots makes Hen 2-90 a unique object. These characteristics could point to a compact planetary nebula, since many planetary nebulae present a bipolar nebula. This fact implies that the asymmetries in the wind can already take place in the late stages of the AGB phase, where the jets can probably shape the spherical AGB wind in a bipolar scenario (Sahai & Trauger 1998; Imai et al. 2002; Vinkovic et al. 2004). However, the presence of a disk and jets could also be explained assuming that the system is a binary (Sahai & Nyman 2000). Guerrero et al. (2001) also noted based on near-IR colors, that Hen 2-90 might be a symbiotic system, i.e. a binary composed of a cool giant and a hot component with an accretion disk. Due to these two completely different interpretations (compact planetary nebula or symbiotic object) the real nature of Hen 2-90 remains unclear.
We study this object by means of high-resolution spectroscopy and low resolution spectrophotometry data obtained with the FEROS and B&C spectrographs, respectively, at the 1.52 m telescope in ESO (La Silla, Chile). The goals of this study are:
High and low resolution observations were obtained with the Fiber-fed
Extended Range Optical Spectrograph (FEROS) and with the Boller & Chivens (B&C) spectrograph, respectively, at the ESO 1.52 m telescope in La Silla
(Chile). FEROS is a bench-mounted Echelle spectrograph with the fibers
located at the Cassegrain focus with a spectral resolution of
R = 48 000corresponding to 2.2 pixels of 15 m and with a wavelength coverage
from 3600 Å to 9200 Å. The aperture of FEROS through which the fibers
are illuminated is 2
.
High resolution observations of Hen 2-90 (higher than the spectra previously described in the literature) were taken on June 10, 2000, with an exposure time of 4800 s. These data were used to get a spectral atlas and a detailed description of the line profiles and also to deblend some features present in the low resolution spectra. The S/N ratio in the 5500 Å region is approximately 20. FEROS has a complete automatic on-line reduction, which we adopted. Equivalent widths were measured using an IRAF task that computes the line area above the adopted continuum.
Low resolution spectra of Hen 2-90 were taken on June 9, 2000, with an
exposure time of 900 s. The slit width was 2
.
The instrumental
setup made use of grating #23 with 600 l mm-1, providing a
resolution of
4.6 Å in the range 3800-8700 Å. The
efficiency of the CCD is a function of wavelength. It increases from about 70% at 3500 Å to its maximum value of 90% at about 7000 Å and then
decreases again, reaching 75% at 8000 Å. This behaviour results in larger
flux uncertainties of about 20% especially in the very low (<5000 Å) and
very high (>7500 Å) wavelength regions of the spectrum, while in the range
in between the flux uncertainties are more on the order of 10%.
In the 5500 Å region, the S/N ratio in the
continuum is approximately 40 for the B&C spectrum (hereafter Cassegrain
spectrum). Since there is no completely line-free region in the spectrum, the
S/N derived is an upper limit. The Cassegrain spectra were reduced using standard
IRAF tasks, such as bias subtraction, flat-field normalization, and wavelength
calibration. We did flux calibrations and extinction corrections with
E(B-V) = 1.3 (Costa et al. 1993). Spectrophotometric standards from
Hamuy et al. (1994) were used for absolute flux calibration. In the
linearized spectra the line intensities were measured with the conventional
method of adjusting a Gaussian function to the profile. Another source of uncertainty
in the line intensities is the position of the underlying continuum and we
estimated the errors to be about 20% for the weakest lines (line fluxes
10 on a scale with H
)
and about 10% for the strongest lines.
![]() |
Figure 1: Extinction corrected optical spectrum of Hen 2-90, observed with the B&C spectrograph. |
Open with DEXTER |
In order to identify the lines and make a spectral atlas of the optical region, we used the line lists provided by Moore (1945), Thackeray (1967), and Landaberry et al. (2001). We have also looked up two sites on the web: NIST Atomic Spectra Database Lines Form (URL physics.nist.gov/cgi-bin/AtData/lines_form) and The Atomic Line List v2.04 (URL www.pa.uky.edu/peter/atomic/).
Figure 1 shows the extinction corrected low resolution Cassegrain spectrum,
which we can see
is dominated by emission lines superimposed on a flat continuum. No absorption line
is present. Table 4, which is only available in the online version
of the journal, shows many emission lines that we have identified. There, the observed
wavelength (Col. 1), the observed intensity
(
(
), Col. 2), the extinction corrected intensity (
(
), Col. 3), and the proposed identification (Col. 4) for each
line are given. The intensities are relative to H
with an observed
H
flux of
ergs cm-2 s-1 and an
extinction-corrected H
flux of
ergs cm-2 s-1.
The line identification given in Col. 4 of Table 4 encloses
the ionization state of the element with the proposed transition and multiplet,
as well as the rest wavelength of the transition. It is possible
that more than one ion can be allocated to a single feature. In these
cases, we give some possible alternative identifications. For some lines no
identification could be found, so they remain unidentified,
labelled as "Uid'' in the table.
Many permitted and forbidden lines were identified, most of them coming from singly or doubly ionized elements. Iron is by far the element with the richest spectrum, and [Fe III] lines are the strongest ones of this element, as cited by Guerrero et al. (2001). On the other hand, our data show an even larger number of permitted and forbidden lines of Fe II (although less intense than the [Fe III] lines), not reported by Guerrero et al. (2001).
The presence of very intense H I Balmer lines, He I,
[O II], [O III], [N II], and [S III] lines is remarkable. H
is the most intense line, followed by [O III]
5007. This spectral characteristic is very
curious, because it is different from that usually seen in a typical
planetary nebula, where the [O III]
5007 is the most intense
line. It also differs from the spectrum of a low-excitation planetary nebulae,
where the H
,
H
,
and [N II]
6548,6583 are
more intense than [O III]
5007 (Kwok 2000). Guerrero et al.
(2001) could neither confirm nor deny the presence of TiO absorption
bands, the main signature of a symbiotic system, coming from the cool component.
However, our spectrum clearly shows that the TiO bands are not present. In addition,
we also could not identify the He II lines that would come from a hot component.
![]() |
Figure 2: The [O III] and [S III] lines taken from the Hen 2-90 high-resolution (FEROS) spectrum. They are only two examples of the presence of double-peaked emission line profiles. The velocities obtained from the wings of these lines are given in Table 1. Due to the strong lines, the high velocity wings are not resolved in this figure. |
Open with DEXTER |
Identification of many Fe II permitted and forbidden lines, H I Balmer emission lines, and also [O I] lines confirms the presence of the
B[e] phenomenon, described by Lamers et al. (1998).
Concerning the line profiles, Guerrero et al. (2001) grouped the emission
lines of Hen 2-90 into three groups: lines with (1) double peaks, (2) broad single
peaks (
km s-1), and (3) narrow single peaks (
km s-1). Using our FEROS spectrum, we confirm the presence
of those groups and define a new one, (4) the lines that clearly show a "shoulder'',
an intermediate case between the double peaks and the single peak.
We list the lines falling into each profile group shortly and show an example:
(1) Double-peaked profiles are shown by the forbidden lines of O II,
O III (left panel in Fig. 2), Ar III, Fe III, S
III (right panel in Fig. 2), and Cl III with a peak separation of 40 km s-1.
Also H
shows a double-peak structure (Fig. 3).
It is interesting to note that the asymmetric double-peaked lines have
a more intense red peak than blue, indicating
that the receding material is brighter than the approaching one.
This finding agrees with Guerrero et al. (2001). However, contrary
to these authors, our [O III] lines (except of the 4363 Å line) show
double peaks with equal strength. The difference between the line profiles (and also
the line intensities) is probably due to the different slit widths used for the
observations.
![]() |
Figure 3: The H I Balmer lines taken from the Hen 2-90 high-resolution (FEROS) spectrum. These lines show an evolution from single-peaked profiles over profiles with a shoulder to double-peaked profiles with decreasing quantum number. |
Open with DEXTER |
![]() |
Figure 4:
Wings of the Balmer lines taken from the Hen 2-90 high-resolution (FEROS)
spectrum. They show an increase in velocity from H8 to H![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 5: [S II] line profiles present in the high-resolution spectrum (FEROS) of Hen 2-90. The derived wing velocities are given in Table 1. |
Open with DEXTER |
![]() |
Figure 6:
The [N II] ![]() ![]() |
Open with DEXTER |
The double-peaked forbidden lines show, in general, emission wings extending from 45
to around 100 km s-1 (except of H), indicating the presence of a low
velocity wind. H
,
however, shows wings extending up to 1800 km s-1
(Fig. 4). Whether these wings indicate an outflow velocity or are produced
due to electron scattering will be discussed in Sect. 5.4.
The velocity extent of the H wings, in our data, is higher than the
1050 km s-1 cited by Costa et al. (1993) and the 1500 km s-1found by Guerrero et al. (2001). The differences in these measurements
are due to different resolutions of the spectra and due to the different
slit widths. However, only more observations with a similar instrumental setup can
discard a spectral variation.
Another difference worth mentioning between our spectra and the data presented by
Guerrero et al. (2001) is the absence of the [S III] 6312 emission in their spectra. This line is very strong in our spectrum and shows double
peaks (see Fig. 2). A possible explanation might be that, since
this emission line arises in the outer parts of the intermediate wind (at least
according to our non-spherical wind model, see Sect. 5.3 and
Fig. 9), the slit width used by Guerrero et al. (2001) was
too narrow to detect this line.
(2) Broad single-peaked profiles are found for H,
H
,
H8 (Fig. 3), Paschen lines, and some He I lines.
The FWHM of these lines is typically
40 km s-1.
(3) Narrow single-peaked profiles are shown, especially by the forbidden lines
of Fe II and S II (Fig. 5). These lines typically
have
km s-1.
(4) Profiles with a "shoulder'' or almost a second peak, which are in the
new group identified by us. This group is represented by some He I, [N
II], and [O I] lines (Fig. 6), as well as by H
and H
(Fig. 3). These lines are asymmetrical, and in some cases their blue
and red wings show different velocities (see Table 1).
It is interesting to note that the H I Balmer lines show an evolution
from a single peak in the high order lines to a double peak in H,
having
an intermediate case with shoulder in H
and H
(Fig. 3).
Table 1 shows the velocities derived from the maximum extend of the wings of
several forbidden emission lines and some Balmer lines in our FEROS spectrum. There,
the ion identification (Col. 1), laboratory wavelength ( (Å), Col. 2), blue wing
velocity (
(km s-1), Col. 3), and red wing velocity (
(km s-1), Col. 4) for each line are given. Included are only lines whose
velocities could be determined with high accuracy and which are not blended.
Table 1: Derived maximum velocities where the blue and red wings reach the continuum. Listed are some forbidden emission lines and some Balmer lines in the FEROS spectrum whose maximum wing velocities could be determined well and which are not blended.
In summary, we can say that compared with the Guerrero et al. (2001)
data, our spectra revealed many more emission lines due to the higher signal-to-noise
ratio (see Table 4).
The fact that with the same applied reddening correction value some of the line
profiles and intensities seem to have changed over a period of 6 months needs some
further clarification. Our line profiles as shown were taken
with FEROS using the fiber technique and with an aperture of 2
.
A slit width of 2
was also used for the Cassegrain spectrum
from which we derived the line intensities.
Guerrero et al. (2001) used a slit width of 1
,
i.e.
only half of our value. They therefore observed a region very close to the star,
while our spectra cover a much larger region of the circumstellar material.
In addition, the dispersion and signal-to-noise ratio of the spectra are different,
which makes it difficult to compare individual lines. The question of line variability
can thus only be answered by subsequent observations with identical setups,
and differences between the Guerrero et al. (2001) spectrum and our data
should not be interpreted as source variability. What we can conclude, however,
is that the circumstellar medium around Hen 2-90 is not homogeneous. On the contrary,
the huge zoo of observed emission lines gives us a clear hint that the circumstellar
medium must have a rather complex structure in both density and temperature,
as indicated by the variety of line profiles, velocities, and ionization states.
![]() |
Figure 7:
Top panel: HST image of Hen 2-90 replotted with kind permission from Sahai
et al. (2002). Overplotted is a big circle that indicates
the width of the slit (which is 2
![]() ![]() |
Open with DEXTER |
Emission lines: we observe from the innermost regions in Hen 2-90 permitted emission lines, as well as forbidden emission lines. The huge amount of forbidden lines might speak in favour of a nebula nature of the emitting material. However, we cannot explain the evenly huge amount of permitted lines with the nebula approximation because both sorts of lines need completely different density conditions: permitted lines are produced in regions of high density, forbidden lines are produced in regions of low density. The fact that both types of lines are very prominent in our spectra indicates that we need at least two different regions, one with high and one with low density. This might be consistent with the picture of a high density (ring?) nebula in the equatorial region and a low density nebula in polar region.
Ionization structure: inspection of the HST image (Fig. 7) shows clearly high-ionized material in the polar direction (= jet direction in the image), i.e. in the presumably low-density region, dominated by emission of [O III], and low-ionized material at more intermediate latitudes (counted from the pole), dominated by [N II] emission. In our spectra, we found many forbidden emission lines as listed in Table 4. These lines are from ions in different ionization states and must therefore arise in regions of different ionization and, therefore, different temperature, consistent with the HST image. In addition, we found very strong emission of [O I], which means that this emission region must be neutral, also in hydrogen. The emission of the [O I] lines therefore must arise in (i) a (hydrogen) neutral region, speaking for high density material, at least close to the star, which guarantees shielding from the ionizing stellar flux, and (ii) in a medium of low density for the forbidden lines to be strong. The best location is the equatorial neutral disk-like structure visible in the HST image.
Velocity structure: our high-resolution spectra allow us also to derive the velocities from individual line profiles (see Table 1). We find that forbidden lines of different ionization stages, which can be connected with the different ionization regions seen in the HST image, have different velocities. The high-ionized lines thereby have the highest velocities, while the [O I] lines have the lowest.
The combination of all these results i.e. - (i) the co-existence of high- and low-density regions; (ii) a latitudal velocity structure having highest values in polar and lowest values in equatorial directions; (iii) a latitude-dependent ionization structure of the circumstellar material which also infers a latitudal temperature structure; and (iv) the variety of line profiles described in the previous section - seems to be more consistent with a non-spherical stellar wind scenario than with a simple constant density nebula scenario or even a piece-wise constant density nebula approximation. This conclusion was drawn solely from qualitative analysis of our observations in combination with the information that follows from the HST image. As we show in Sect. 5.3 when modeling the line luminosities of the forbidden emission lines, it turns out that saturation of the line luminosities happens for different lines at different distances from the star. This means that the forbidden lines are very sensitive to the density structure, supporting the idea of a radial density distribution rather than a constant density.
In the previous sections we have shown that the optical spectrum of Hen 2-90
contains many forbidden and permitted emission lines. Here we
analyse some of the forbidden emission lines to derive the physical parameters
of the emitting gas, such as electron temperature and density distribution, mass loss
rate of the star, and the ionization structure of the circumstellar nebula.
The stellar parameters of Hen 2-90 given by Costa et al. (1993) were
derived using the spherically symmetric nebula approximation and are therefore
relatively uncertain. Due to the lack of better parameter estimates we assume the star to be at a mean distance of 2 kpc for our model calculations, having an effective
temperature of about 50 000 K and a radius of 0.38
(Kaler & Jacoby
1991; Costa et al. 1993).
Based on the HST image of Sahai et al. (2002, see Fig. 7) we can
distinguish three major wind regions: (i) a biconical high-ionization
([O III] dominated) polar wind; (ii) a low-ionization ([N II] dominated)
wind at intermediate latitudes, which we will call the intermediate wind;
and (iii) an equatorial outflowing disk. The big white circle indicates
the slit width of 2
of our Cassegrain observations (the flux observed
comes from this region). This corresponds to an outer edge
of about 0.01 pc at a distance of 2 kpc. The HST resolved structure is embedded
into a much larger H II region with extensions of at least
(Tylenda et al. 2003). Our observations therefore contain
neither contribution from this extended nebula nor from the jets and knots, but
are concentrated on the innermost non-spherical wind structure.
The HST image also hints at the location of the different ions. In
combination with the velocity information retrieved from the line profiles
(Table 1) we can determine the terminal wind velocities of the different
regions. For the polar wind we find a mean velocity of 60 km s-1
from the [O III] and [Cl III] lines. Taking the opening angle of the cone
into account, which is about 90
,
along with the fact that the system is seen
(almost perfectly) edge-on, results in a terminal velocity of
85 km s-1. For the intermediate wind where the [S III], [S II],
and [N II] lines come from, we analogously derive a mean wind velocity of
45 km s-1; and for the disk wind, represented by the [O I] lines,
we find a mean velocity of about 35 km s-1.
The wind velocity shows a latitude dependence highest in the polar directions and
lowest in equatorial direction, and the difference is a factor of
2.5.
With the adopted distance of 2 kpc we can derive linear scales for the different
wind regions from the HST image. The polar wind is visible up to a distance of about
R* which is slightly beyond our observing
slit, which only extends to about
pc.
The outer edge of the disk is about
R*and we estimate a total opening angle of
,
based on the dark disk-like
structure seen on the HST image. The remaining volume is filled by the intermediate
wind, which also shows indication of a high degree of ionization ([O III]
dominated) at its innermost parts. We estimate that this region extends roughly to
about
O III
R*.
We concentrate our detailed analysis on the forbidden lines, which are excellent indicators of the circumstellar material for two reasons:
The line luminosity of a forbidden line is given by
![]() |
(1) |
The hydrogen density in a wind follows from the mass continuity equation
![]() |
(2) |
![]() |
(3) |
![]() |
(4) |
The electron density is given in terms of the hydrogen density and depends on the
degree of ionization in each part of the wind; in polar direction and in the inner
parts of the intermediate wind we set
.
In these
regions helium is singly ionized contributing about 10% to the total
electron density. The contribution from other elements, like O, S, N, which are also
(partly twice) ionized, to the total electron density is, however, negligible.
For the outer parts of the intermediate wind, where He is assumed to be neutral, we use
(and again neglect the small contribution from
the metals). For the disk wind where we assume that
hydrogen is neutral, we have only electrons from elements with much lower
ionization potential (i.e. lower than about 10 eV). Summing up the number
densities of all these elements (assuming solar abundances), results in a maximum
electron density in the disk of about
.
As we cannot
expect that all possible elements are completely ionized, we use an electron density of
,
which is an upper limit to the real disk
electron density. Since the
excitation of the forbidden lines depends on the electron density, the fitting
of the disk lines therefore results in a lower limit of the mass flux,
i.e. the disk mass loss rate, because a decrease in electron density needs an increase
in total density to account for the observed line luminosities.
For the temperature distribution we make use of the following equation
One additional parameter that needs to be specified is the elemental abundance. Costa et al. (1993) have found that the elements are slightly sub-solar in Hen 2-90. However, these values have been derived under the assumption of a spherically symmetric homogeneous nebula, which is far from being the case. In addition, Kraus (2005) shows that the abundances can be much higher when derived using a stellar wind rather than a homogeneous nebula. We therefore start our calculations by assuming solar abundances (taken from Grevesse & Sauval 1998) for all elements. Deviations in the calculated line luminosities compared with observations might then be due to individual deviations in the abundances.
![]() |
Figure 8:
Temperature distribution in the wind according to Eq. (5).
The starting values T0 and terminal values
![]() |
Open with DEXTER |
![]() |
Figure 9: Model results for the [S III] line arising in the outer parts of the intermediate wind. Shown is the increase in line luminosity with distance from the star. The arrow on the right side of the panel indicates the observed value. |
Open with DEXTER |
We started our modeling with the [S III] 6312 Å line
(Fig. 9). Since S III
and O II have about the same ionization potential,
it can be expected that no S III emission will come from regions
of O III emission. S III emission is therefore restricted to the
intermediate wind at distances from the star where O III, and hence S
IV, have recombined already. Therefore, the inner edge of the S III emission
is the outer edge of the O III emission, i.e.
(S III
(O III
) = 105 R* as found in Sect. 5.1.
The only free parameters are the temperature distribution
and the mass loss rate at these intermediate latitudes. Since the S III emission
region is already rather far away from the star, we assumed the temperature to be
constant and found a good fit for
K. The mass flux of
the intermediate wind has been found to be on the order of
g s-1 cm
yr-1 steradian-1.
Since S III has only one prominent forbidden line in our spectral region, this best-fit temperature had to be confirmed with the other forbidden lines from the same wind region. Other elements like N II, S II, and Cl II have more than one forbidden line, and several line ratios can be used as tracers for the electron temperature. From combining the results from all these temperature indicators, the best solution was found to be on the order of 10 000 K. The valid range of temperatures is very small, and the resulting uncertainty in mass flux is on the order of 20%. We want to stress that the dependence of the mass flux on electron temperature cannot be described with a power law, because the collision parameters have a more complicated temperature dependence.
![]() |
Figure 10: Model results for the different [O III] lines listed in the bottom panel. Shown are individual contributions from the polar and intermediate wind and their increase in line luminosity with distance from the star. The bottom panel shows the sum of both contributions, and the arrows on the right side indicate the observed value for each emission line (same line style). |
Open with DEXTER |
![]() |
Figure 11: Model results for the [O II] lines. These lines come from the intermediate wind only and from distances where O III has already recombined (see mid-panel of Fig. 10). The arrows on the right side indicate the observed value for each emission line (same line style). |
Open with DEXTER |
![]() |
Figure 12: Model results for the [O I] lines. These lines come from the disk wind only. The arrows on the right side indicate the observed value for each emission line (same line style). |
Open with DEXTER |
The next element to fit is oxygen, so we started with the lines of O III. These lines are supposed to come from 2 different regions (see Figs. 7 and 10): the inner parts of the intermediate wind and the polar wind. To get self-consistent results for all oxygen lines of different ionization state (Figs. 11 and 12), we had to decrease the O abundance from solar to 0.3 solar.
Contribution of the intermediate wind to the [O III] lines is restricted
to distances within
(O III). For the temperature
distribution we used
and
K.
Using the effective stellar temperature as the starting temperature is quite
reasonable because the star is supposed to be small, which means that the surface
temperature is about equal to the effective temperature. In fact, the choice of the
starting temperature does not really affect the behaviour of the line luminosities
because the majority of the luminosity comes from the lower density regions
where the temperature has reached its terminal value.
Once we have the contribution from the intermediate wind, the remaining line
luminosities to the [O III] lines must come from the polar wind. Here, the
free parameters are again the mass flux and the temperature distribution. We
found the best agreement by using the same temperature distribution as
for the intermediate wind and a polar mass flux of
g s-1 cm
yr-1 steradian-1.
With these fixed parameters for the intermediate wind, we found good
agreement between observed and modelled [O II] lines (Fig. 11).
The remaining O lines are the [O I] lines, which are
assumed to arise in the cool and neutral equatorial disk (Fig. 12).
The only constraint on the
modeling is the temperature distribution, which should start with
a surface temperature around 10 000 K to guarantee that H is neutral. The terminal
temperature, as well as the disk mass flux, could then be fixed by the best fit model
(see Fig. 12) to
K and
g s-1 cm
yr-1 steradian-1.
All the parameters for the three different wind regions are now fixed and summarised in Table 2. If this picture of the non-spherical wind is about correct, we should be able to fit the luminosities of the remaining forbidden lines (i.e. from S II, Cl II, Cl III, Ar III, and N II) with the same sets of parameters.
Table 2:
Best fit model parameters for the polar (p), intermediate (i), and disk wind (d). The mass fluxes are assumed to be constant within the indicated -ranges
and have an error on the order of 20%.
Table 3: Observed and modelled forbidden line luminosities. For the observed values, we took the intensities given in Table 4 and assumed a distance of 2 kpc. The fifth column contains the ratio of observed over modelled line luminosity, and in the sixth column we describe in which part(s) of the wind the line is produced. The classifications of p, i, and d mean polar wind, intermediate wind, and disk, respectively. In the upper part we gathered the lines used to derive the mass fluxes and temperature distributions in each wind zone, while the lines in the lower part have been used to confirm the models and to derive the abundance, especially in the case of N.
The modeling of the [S II] lines was rather tricky, because they arise
in two different regions: the hydrogen neutral disk and the outer parts of the
(ionized) intermediate wind, but with no clear boundary conditions.
For the intermediate wind, S II can only arise if S III has recombined.
Saturation of the [S III] line luminosities happens at about 106 R* (Fig. 9), which is set as the minimum distance for
the [S II] lines. For emission from the disk the inner edge is set to the
stellar surface. A reliable fit is found from the simultaneous fitting of both
contributions where we found that the disk contributes less than half to
the total emission and that S II has to recombine at a distance of about 6000 R*. Contribution from the intermediate wind, which is more than half
of the total line luminosities, comes from a very narrow (in radius) region, starting at
a distance of about
just beyond the region of S III line luminosity saturation and extending to the edge of our slit, which is at
.
Since the line luminosities from the intermediate wind
did not saturate at the outer edge of our slit, we expect the [S II] lines would
appear much more luminous when observed with a larger slit width.
The predicted line luminosities of the nitrogen lines were much too high, so that we
had to reduce the N abundance to
solar to achieve good fits.
In Table 3 we summarise all the observed (Col. 3) and modelled (Col. 4) line luminosities, as well as the ratio of observed over modelled luminosity (Col. 5) and the location(s) of the ions in either the polar wind (p), intermediate wind (i) or disk (d) alone, or a combination of two of them. We thereby distinguish two sets of lines: (1) those which are used to derive the mass fluxes and temperature distributions in the different wind regions are collected in the upper part, and (2) the lines that are used for confirmation of the models and for derivation of the abundance (especially of N) are given in the lower part.
Inspecting Table 3 shows that some of the modelled lines can be off by a factor of 2 or larger. Here, we give some arguments that might explain these differences.
From the results presented in Table 2 for the different wind parts, it
is obvious that the mass flux varies with latitude
being lowest at the poles and highest at the equator with
.
The total mass loss rate of the star is given by the
integral of the mass flux over the stellar surface, which is
![]() |
(6) |
With this complete wind scenario derived in the previous section we can return to
the question whether the broad wings of the H line are due to outflow or
electron scattering. From our data we derive a contribution of about 10%
of the broad wings to the total intensity, which implies an electron scattering
optical depth of about 0.09 in and above the line formation region of H
.
On the other hand, we can calculate the electron optical depth from
![]() |
(7) |
Additional confirmation of the broadening due to electron scattering rather than due to high velocity outflow is provided by the wind scenario proposed above. Our non-spherical wind model covers the complete surroundings of the star except for the highly-collimated jets, which have been found to have real outflow velocities of less than 400 km s-1 (see e.g. Sahai et al. 2002). What we found from our wing velocity measurements in all the different parts of the wind is: a terminal outflow velocity of (i) about 85 km s-1 from the [O III] lines in the polar wind, (ii) about 40 km s-1 from the [N II] and [S III] lines in the intermediate wind, and (iii) about 35 km s-1 from the [S II] and [O I] lines in the disk. We could not identify any additional wind component that might have an outflow velocity as high as 1800 km s-1.
For the model calculations we use the geometry and ionization structure based on the HST image (Fig. 7), which clearly shows a polar wind, an intermediate wind, and a disk so that the geometry used seems to be well justified.
Such a non-spherical wind model with latitude-dependent temperature, terminal wind velocity, and mass flux might be understood in terms of a rotating star. According to the van Zeipel theorem, rapid rotation will lead to polar brightening, which means that the surface temperature of the star can be much hotter on the pole compared to the equator. Our calculations show that in equatorial regions the surface temperature should not be higher than about 10 000 K to guarantee that hydrogen is neutral. In polar directions we found that the terminal temperature should not be lower than about 10 000 K, which means that the surface temperature can be much hotter. We do not have a good handle on the real polar surface temperature because the forbidden emission lines are produced far away from the star where the temperature has dropped already to its terminal value. In our calculations we used the effective temperature of 50 000 K found by Kaler & Jacoby (1991) with the Zanstra method. Cidale et al. (2001) derived a much lower effective temperature of only 32 000 K using the energy balance method. We tested our model for different polar surface temperatures in the range T0 = 15 000-50 000 K and found that, for different surface temperatures, large differences occur in the line luminosities close to the star. However, in the regions of our interest, where the lines saturate, the line luminosities are all about equal. This implies that from our observations we cannot draw any conclusion about the real polar surface temperature.
Rotation not only influences the surface temperature structure, but also the terminal wind velocities resulting in a decrease in terminal velocity from pole to equator (see e.g. Lamers & Cassinelli 1999; Kraus & Lamers 2005). This is exactly what we found from our observations, when we take the known inclination of the system into account as suggested from the HST image.
The ratio of polar (85 km s-1) to disk terminal velocity (35 km s-1) of about 2.4, derived from the forbidden emission lines, is consistent with a rotation speed of about 75-80% of the critical velocity (Kraus & Lamers 2005). Such high velocities are observed e.g. for the classical Be stars (see e.g. Porter & Rivinius 2003). Some of those stars have recently been suggested as rotating even close to break-up, i.e. with almost critical velocity (Townsend et al. 2004).
An additional effect of rotation is that it leads to a latitude-dependent mass flux, which is also the case for Hen 2-90. A mechanism like the wind compressed disk (Bjorkman & Cassinelli 1993) or the rotation-induced bistability (Lamers & Pauldrach 1991), or a combination of both, can result in the observed high density equatorial outflowing disk.
This combination of the latitude-dependent mass flux and surface temperature also explains why we see different ionization structures at different latitudes. The lower surface temperature results in less ionizing UV photons, and the higher mass flux into a higher shielding density, which in the end reduces the degree of ionization.
For our model calculations the outflowing disk is assumed to be neutral in hydrogen. Kraus & Lamers (2003) showed that in the case of supergiants, the disk created by an equatorial wind can indeed become neutral just above the stellar surface, if the equatorial mass flux is high enough to prevent ionizing stellar photons from penetrating into the high density equatorial disk wind. Although Hen 2-90 is hotter than a typical supergiant (but its effective temperature might well be only 32 000 K instead of 50 000 K as stated above), it also has a much smaller radius leading to an enormous equatorial mass flux, whose material is indeed dense enough to shield the disk material efficiently from stellar irradiation. Consequently, the disk will stay neutral even close to the stellar surface.
The terminal temperature found for the gas disk is about 6000 K. From our modeling
we found that the luminosities of the forbidden lines coming from the disk saturate at
a distance of about
R*. The total size of the disk is, however,
much larger, at least
R*. The ISO spectrum of Hen 2-90 shows
strong IR excess emission due to hot circumstellar dust that can only be located in the disk.
The evaporation temperature of dust is between 1500 K and 2000 K depending on its
composition. In the temperature region between the atomic and the dust dominated regions,
the atoms might be locked into molecules that can form (and survive) at temperatures
below about 5000 K. Our results are therefore also consistent with a circumstellar
dust disk where the dust forms far away from the star, at distances
AU. A detailed study of the ISO spectrum
to calculate the infrared emission coming from the predicted dust disk is necessary
and currently under investigation (Borges Fernandes & Lorenz-Martins, in preparation).
It would be interesting to see whether the ionization structure we found from our modeling of the forbidden emission lines might be reproduced by more sophisticated ionization structure calculations in such a non-spherically symmetric wind scenario. Unfortunately, as far as we know, no numerical codes exist that can perform the 3D radiation transfer needed for non-spherical ionization structure calculations.
Finally, we qualitatively discuss how the variety of the line profiles, as shown in Sect. 3, might be produced in our non-spherical wind scenario. Even though the central star might be rapidly rotating, we do not think that the double-peaked line profiles mirror the stellar rotation or are due to rotation at all (although we cannot ad hoc exclude a rotational contribution), because the effective temperature of about 50 000 K indicates that the star must have a line-driven wind. The acceleration of the wind to its terminal velocity happens within only a few stellar radii, leading to a radially outflowing wind.
Double-peaked lines: these lines are mainly produced in the intermediate wind (some of them with a Gaussian contribution from the polar wind), which can be regarded as a thick expanding torus. Since we see this torus edge-on, the emission lines formed within it show double-peaked profiles. Since the maximum (terminal) wind velocity of the intermediate wind is about 45 km s-1, the peak separation will be on the order of 40 to 50 km s-1, as observed.
Broad single-peaked lines:
they are found for the hydrogen Balmer lines. As discussed
already in Sect. 5.4, the broad wings of H
are due
to electron scattering with a scattered intensity of about 10%.
With increasing quantum number, the Balmer line intensities decrease and
a 10% contribution in the broad wings becomes much harder to detect above
the continuum. That's why with increasing quantum number the derived wing
velocities of the Balmer lines decrease until they adapt to the
terminal wind velocities derived from the forbidden lines (see e.g. line H8
in Table 1).
Narrow single-peaked lines: these lines are mainly formed in the disk, even though some of them can have small contributions from the intermediate wind. The lines formed in the outflowing edge-on seen disk should also show double-peaks as do the lines formed in the intermediate wind. However, since the disk outflow velocity is smaller than the velocity in the intermediate wind, the two peaks merge, leading to a single-peaked line profile.
The variety of profiles observed can be explained with our model assumptions. However, we have no explanation for those asymmetric line profiles with a strong red peak. Probably, our scenario is not completely correct, and effects like clumping in the wind might be worth investigating.
A quantitative analysis of the line profiles would be necessary, but this is difficult with our present data and beyond the scope of this paper, because the profiles are derived from the FEROS spectrum which is not flux calibrated, while the line luminosities have been derived from the low-resolution Cassegrain spectrum, which makes it impossible to model line luminosities and profiles simultaneously.
In the literature Hen 2-90 has been classified either as a (compact) planetary nebula (e.g. Costa et al. 1993; Maciel 1993; Lamers et al. 1998) or as a symbiotic object, i.e. a binary consisting of a cool giant and a hot component with an accretion disk (Sahai & Nyman 2000; Guerrero et al. 2001).
The HST image (Fig. 7, Sahai et al. 2002) clearly shows the presence of a nebula bisected by a disk, with a highly collimated and bipolar jet with several pairs of knots on both sides. This scenario could be explained assuming that Hen 2-90 is either a compact planetary nebula, where the wind asymmetries started during the AGB phase, or a binary system, where the jets are caused by an accretion disk.
During our modelling procedure we found, as also cited by Costa et al. (1993)
and Maciel (1993), a depletion of N and O by a factor of 2 and 3.3,
respectively, with respect to the electron density with S, Ar, and Cl having
normal (i.e. solar) abundances. In addition, the FEROS spectrum contains only one
single permitted carbon line, C II 6578. This line has previously been
detected by Costa et al. (1993) and Guerrero et al. (2001), but it is
very weak in the FEROS spectrum (just at the detection limit) and absent
in the Cassegrain spectrum. We can only estimate its line flux and we found
about 0.06 for the observed flux and about 0.02 for the dereddened flux in the scaling
of Table 4 (H
). The uncertainty of these line fluxes is
about 50%. No forbidden emission line of carbon was detected, even though
in the optical spectrum we would expect to observe several [C I] lines coming
from the disk. We applied our disk wind model to calculate the
line luminosities for these carbon lines, and found that these lines would only show up
in the spectrum for a carbon abundance of at least twice solar.
Inspection of the ISO spectrum of Hen 2-90 reveals emission from oxygen-rich
dust (Borges Fernandes & Lorenz-Martins, in preparation). This dust can
only be the dominant dust component in the case
during the ejection
time of the dust-forming material. We found an underabundance of O by more
than a factor of 3. Together with the fact that the dust composition should
mirror the actual wind composition (unless the star has changed its
composition from O-rich to C-rich since the ejection of the dust forming
material, which is rather unlikely),
we conclude that (i) the progenitor star that ejected the dust forming material
cannot have been carbon rich, and (ii) the actual carbon
abundance must be less than 0.6 solar to guarantee
.
Assuming a standard scenario with dredge-ups, we can understand the depletion
of C and O; however, we have no explanation for the fact that N is also
depleted. There are some B-type post-AGB stars in the halo of our Galaxy
which also show this behaviour (Moehler & Heber 1998), and recently Lennon et al. (2005)
reported on a Be star with an unexpected low N abundance. However, the origin of those
anomalies is still poorly understood. Rotation or a possible binarity could play
an important role. Unfortunately, as far as we know, no stellar evolution models exist
that follow the whole life sequence of rapidly rotating intermediate mass stars,
for which Hen 2-90 seems to be a candidate.
Another point that should be noted is that from our spectra, which cover the
complete optical wavelength range, we have identified no lines from ions with ionization
potential higher than 40 eV and consequently no He II lines. In addition,
there is no indication for the presence of TiO bands either. These two characteristics
are the main signature of a symbiotic system. Since we see the disk of Hen 2-90
almost edge-on, it might be possible that the disk hides the atmosphere of a cool giant,
where the TiO bands are formed, and also absorbs the He II emission from a hot
component. However, the fact that we don't see hints for any other line from ions
with ionization potential higher than
40 eV leads us to the conclusion that
these ions (and therefore also He II) do not exist in the wind of Hen 2-90
indicating that (i) the effective temperature of Hen 2-90 is indeed much lower than
the 50 000 K found from the Zanstra method and may well be around 32 000 K as
derived by Cidale et al. (2001) from the energy balance, and (ii) the
symbiotic system is not the favoured classification.
Besides the absence of clear indications for a classification as a symbiotic object, our results favour the conclusion of Hen 2-90 being a compact planetary nebula and two more features support this idea:
Observed features that might hint at a binary nature are the jet-like structure with blob ejections (Sahai 2002; Sahai et al. 2002) and the DENIS NIR photometry cited by Guerrero et al. (2001):
Jet-like structure and blob ejection:
the jet structure seen in Hen 2-90 is very collimated, and the blobs are found to
be ejected regularly on a 40-year timescale. This perfect alignment of the jet
axis and the blobs is unique. All PNe that are known to be binaries
show either asymmetric or point-symmetric jet and blob ejection.
In addition, to date no X-ray emission coming from an accretion disk has been
detected from the system, which might indicate a completely different formation
mechanism.
García-Segura et al. (2001) modeled the jet and knots of Hen 2-90
using the concept of magnetohydrodynamics in the wind of a single star and under
the assumption of a solar-like magnetic cycle. This model does not need a binary
for jet production, because the jet formation is not driven by an accretion disk.
A binary component might be helpful in the first instance to spin-up the star to high
rotation rates (e.g. in a merger process) and, therefore, to increase the stellar
magnetic field to the values needed for their model calculations. Such a spin-up might
also explain the high rotation velocity of the central star which we claimed in the
previous section and derived from the observed terminal velocities.
DENIS NIR photometry: Schmeja & Kimeswenger (2001) used the new DENIS survey data to probe different types of PNe. Their Fig. 1 shows a distinction between genuine PNe, symbiotic Miras, and IR-[WC] PNe on the basis of a colour-colour plot. With the published DENIS photometry data of Hen 2-90 (Schmeja & Kimeswenger 2002), and corrected for the extinction value of Costa et al. (1993), Hen 2-90 falls into the intermediate region between the genuine PNe and the IR-[WC] PNe. However, the IR-[WC] PNe are supposed to have strong PAH emission, which increases their K-band flux and shifts them off the genuine PNe region. A pollution of the K-band emission might also be present for Hen 2-90, which shows a huge zoo of emission features in the IR. Subtraction of this polluting emission might shift Hen 2-90 back into the regime of the genuine PNe. In addition, the carbon abundance found for Hen 2-90 is much too small to classify it as a carbon rich star.
From our observations, modeling, and the above discussion, we conclude that Hen 2-90 must be at least an evolved and probably rapidly rotating object, and the classification as a compact planetary nebula seems to be favourable. The question whether this compact planetary nebula is part of a binary system can only be solved if clear direct indications of a companion star are observed, which is up to now not the case but makes it worth looking at this object in much more detail.
In this paper we studied the non-spherical mass loss history of Hen 2-90 by splitting the analysis into a qualitative and a quantitative part. In the first part, we presented high and low resolution optical observations of the innermost non-spherical wind structure around Hen 2-90. The spectra contain a huge number of forbidden and permitted emission lines of different ionization states justifying the classification of Hen 2-90 as an object showing the B[e] phenomenon. The emission lines can be separated into four groups according to the different shape of their profiles: broad and narrow single-peaked lines, double-peaked lines, and lines with a well pronounced shoulder. This variety of line profiles indicates a complex structure of the circumstellar medium. There are no absorption lines present in the spectrum.
In the second part of the paper we performed a detailed analysis of the forbidden
emission lines. The wind geometry used is based on the HST image, which reveals a
non-spherical wind structure consisting of a polar wind, an outflowing disk, and an
intermediate wind in between them (see Fig. 7). In all three parts, the
ionization structure is different,
indicating a latitude dependence of the surface temperature and the mass flux.
These assumptions could be confirmed by our forbidden line analysis and might be
interpreted in terms of a rapidly rotating (75-80% critical) underlying star.
We could determine
mass fluxes of
g s-1 cm-2,
g s-1cm-2,
and
g s-1 cm-2 (i.e.
yr-1 steradian-1,
yr-1 steradian-1, and
yr-1 steradian-1)
for the polar, intermediate, and disk wind,
respectively. The surface temperature might change from 50 000 K-32 000 K (or less)
at the pole to about 10 000 K at the equator, and the terminal velocities are 85 km s-1 (polar wind), 45 km s-1 (intermediate wind) and 35 km s-1
(disk). The total mass loss rate was found to be
yr-1.
From the near absence of observable carbon lines in our spectra and from the fact that Hen 2-90 clearly shows O-rich dust in the ISO spectrum, we could restrict the C abundance to a value less than about 0.6 solar. In addition we found that the O abundance is 0.3 solar and the N abundance is about 0.5 solar leading to an enhanced N/O ratio of 5/3 with respect to the solar value, while Ar, S, and Cl could be modelled with solar abundances. The depletion of C and O follow from stellar evolution, but for the depletion of N we have no satisfying explanation, as it cannot be explained by a standard stellar evolution scenario with dredge-ups.
From our observations and modeling results we conclude that Hen 2-90 must be an evolved object with most probably a rapidly rotating central star, and we favour the interpretation of Hen 2-90 as a compact planetary nebula. Whether it is part of a binary system is still an unsolved question. More observations and better evolutionary models for rapidly rotating stars are certainly needed for better comprehension of this very curious object.
Acknowledgements
We would like to thank the anonymous referees for remarks and suggestions that have led to a significant improvement of this paper. M.K. also thanks Guillermo García-Segura for helpful discussions on PNe. M.K. was supported by the German Deutsche Forschungsgemeinschaft, DFG grant number Kr 2163/2-1 and by the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO) grant No. 614.000.310. M.B.F. acknowledges financial support from CAPES (Ph.D. scholarship). M.B.F. also acknowledges Utrecht University for the warm hospitality during his one-year stay there.
Table 4:
Emission line intensities relative to H
= 100 obtained
with B&C spectrograph.