A&A 369, 1027-1047 (2001)
DOI: 10.1051/0004-6361:20010156
1 - Institut für Astronomie und Astrophysik, TU Berlin, Sekr. PN 8-1, Hardenbergstraße 36, 10623 Berlin, Germany
2 -
Research School of Astronomy and Astrophysics, Australian National University, Cotter Road, Weston ACT 2611, Australia
Received 30 October 2000 / Accepted 25 January 2001
Abstract
We present time-resolved observations of metallic emission lines of MgI,
MnI, SiI, FeI and FeII, including forbidden emission lines of [FeII],
for the six M-type Mira variables RRSco, RAql, RCar, RLeo, SScl and
RHya, which range in period from 281 to 389 days. Data is also presented for
the Balmer emission lines H,
H
,
H
and H
.
The observations were carried out in the optical wavelength region
3600-5700Å. Narrow-slit observations with dispersion of 1.53 km s-1 pixel-1 and a
3-pixel resolution of 4.6 km s-1 were made in some cases
while most observations were done with a wide slit corresponding to a
resolution of 15.3 km s-1.
The variation of line shape, flux and
velocity with phase is discussed. The data presented in this paper will be
used in subsequent papers for comparison with detailed models of the emission
from shock waves in the upper atmospheric layers of Mira variables. It is in these
important upper layers that dust formation occurs and mass loss is initiated.
Key words: line: profiles - shock waves - stars: circumstellar matter - stars: late-type - stars: variables: general
Two of the outstanding problems concerning M-type Miras are the basic mechanism
for dust formation and the mechanism for the generation of high mass loss rates
of the order of 10-6yr-1 (Habing 1996).
Theoretical works of Wood (1979), Bowen (1988), Fleischer et al.
(1992), Feuchtinger et al. (1993), Höfner & Dorfi
(1997) and Winters et al. (2000) show that the high mass loss
rates of Miras can be explained by a combination of pulsation of the star and
radiation pressure acting on dust. Pulsation generates shock waves which move
out through the atmosphere thereby increasing the mass density in the outer
parts of the atmosphere. The passage of a shock leads to a sudden increase of
the pressure, density and temperature of the gas. Behind the shock the gas
cools via the emission of radiation, which leads to conditions that are
favourable for dust formation. The shocks in the outer parts of the atmosphere
thereby trigger the formation of dust. Radiation pressure acting on this dust
leads to the high mass loss rate.
One way to reveal the thermo- and hydro-dynamical conditions in a Mira
atmospheres is to study the various emission lines which are emitted behind the
shock front and which can be observed over a substantial portion of the
pulsation period. Analysing a time-resolved series of these emission lines
offers the possibility to determine the hydrodynamical conditions in different
layers of the atmosphere influenced by the passing shock wave. In particular,
the hydrodynamical conditions of the outer, dust-producing layers of the
atmosphere can be studied with the help of metal emission lines which appear
late in the pulsation cycle when the shock wave has reached these layers.
M-type Miras show a variety of metal emission lines in all regions of
wavelength. The most substantial observations in the optical wavelength region
are the long-term observation of o Ceti by Joy (1954) and the
observation of several Miras by Merrill (1940, 1945,
1946a, 1946b, 1947a, 1947b). In the
near-IR wavelength region (
m), Gillet et al. 1985a, 1985b,
1985c) observed various metal emission lines of SCar and o Ceti
and interpreted the variation of the lines using a shock model.
Using ISO, Aoki et al. (1998) found a number of high intensity emission
lines of ionized metals in the mid-IR spectral region, although the
stars examined were semi-regular and irregular variables rather than Miras.
Finally, in the
UV spectral region, Wood & Karovska (2000) observed emission lines in M-type Miras and
suggested shock waves as an explanation for their variations.
Around maximum light, the hydrogen Balmer emission lines are prominent in the
optical spectra. Deutsch & Merrill (1959) were the first to interpret
these lines as due to shock waves and the theoretical work of Gorbatski
(1961) verified their interpretation. A number of observational
studies of Balmer emission lines have been made over the years, many
concentrating on the mean radial velocity (Merrill 1940, 1945,1946a,1946b,1947a,1947b; Joy 1947, 1954; Gillet et al.
1985a,1985b,1985c). Quantitive data on Balmer line shapes, widths and line
fluxes, as well as theoretical shock model estimates for the shock speed and
temperature and line fluxes, are presented in Fox et al. (1984,
1985) and Gillet et al. (1983, 1985c).
The emphasis of our work is on various metal emission lines of FeI, FeII,
MgI, MnI and SiI, especially the lines which appear at late phases and
which therefore should originate from the outer parts of the
atmosphere. For two of the stars in our sample (RCar and RLeo) we
observed forbidden lines [FeII] around minimum light.
These lines were observed by Joy (1954) in o Ceti as well as by Merrill (1940,
1947b) in the M-type Miras UOri,
RLeo, RHya and the S-type Mira Cyg. These lines
must occur only very briefly in each cycle or perhaps only in occasional cycles
as we did not see the lines in R Leo during the minimum of March 1999 whereas
the lines were strong during the minimum of December 1999. Similarly, we did
not see the lines during the minima of o Ceti (which we also observed,
but do not present in this paper since no metal emission lines were detected)
or R Hya although the lines
have been reported in the past in these stars. We also note that although Joy
(1954) reports the lines in o Ceti, Merrill (1940)
explicitly notes that he did not see them.
Existing spectral observations have mostly concentrated on the mean radial velocities of the metal emission lines. There is very little published quantitive data on the shapes, widths and fluxes of metal emission lines. Here we present these quantities for a sample of six M-type Miras, namely RAql, RRSco, RCar, RLeo, SScl and RHya, which range in period from 281 days to 389 days. These Miras have been observed at multiple phases of the pulsation cycle. Because of this phase coverage, the data shows the history of the shock as it emerges through the deep photosphere (before maximum light) and then moves out through the atmosphere.
As well as the metal line observations, we also present similar data for the
hydrogen emission lines H,
H
,
H
,
H
.
These
lines are most prominent around maximum light and they provide information
about early history of the specific shock wave that generated the metal lines
we observed.
Star | P[days] | Sp. | ![]() |
RR Sco | 281 | M6e-M9e | -37a |
R Aql | 284 | M5e-M9e | 29b |
R Car | 309 | M4e-M8e | 20c |
R Leo | 310 | M6e-M9.5e | 6b |
S Scl | 363 | M3e-M9e(Tc) | (12) |
R Hya | 389 | M6e-M9eS(Tc) | -11b |
The observations presented here were taken from December 1998 to
December 1999. The observed line spectra were obtained with the coudé
echelle spectrograph and 81 cm camera of the Mount Stromlo Observatory
1.88 m telescope. A Site 24k CCD was used as detector, allowing
coverage of the spectral region 3600-5700Å in each exposure: the
corresponding spectral range in each echelle order was 70-60Å.
The echelle grating gave a dispersion of 0.0204 Å pixel-1 at
4000Å, equivalent to 1.53 km s-1 pixel-1. Maximum resolution
observations were done with a 300
m (1.2
)
slit, giving a 3-pixel
resolution of 4.6 km s-1. For spectophotometric observations, a 1500
m (6
)
slit was used, giving a resolution of 15.3 km s-1.
Because of the high slit losses experienced with the 300
m slit, and the
fact that the emission lines are mostly considerably broader than
15.3 km s-1, a large fraction of all observations were in fact done using
the wide slit.
The periods, spectral types and stellar center-of-mass velocities of the Miras
in our sample are listed in Table 1. The periods and spectral
types are taken from the index to variable stars of the Variable Star Network
(VSNET) (Nogami et al.
1997). The heliocentric stellar center-of-mass velocity
was
determined from the midpoint of the circumstellar molecular emission of CO or
SiO given in the papers cited below the table. This method of deriving stellar
center-of-mass velocities has already been described by Reid & Dickinson
(1976) and the velocities should be correct to 1-2 km s-1. For
SScl it was not possible to find a direct estimate of
in the
literature. We adopted the value given in parentheses: it was obtained by
choosing
to give the same velocity (relative to the stellar
center-of-mass) for various emission lines as found for stars with known
.
This estimate is probably correct to 4 km s-1.
We observed the Balmer lines H
(4340.46Å), H
(4101.73Å), H
(3889.05Å) and H
(3835.38Å) near
maximum visible light when possible in order to get the shock velocity deep in
the atmosphere. Between maximum and minimum light we searched for emission
lines of the metals MgI, SiI, MnI, FeI and FeII. Table 2
lists the multiplet number M, wavelength
,
upper and lower energy
levels of the transition
and
,
and the
Einstein-coefficients
for the metal emission lines observed. This
data was obtained from the NIST
database
(Fuhr et al. 1988). Note that the detected lines include the forbidden lines
of the multiplets [FeII]6F, [FeII]7F and [FeII]21F.
Ion | M | ![]() |
![]() |
![]() |
![]() |
SiI | 2 | 4102.95 | 39760 | 15394 | 9.6+4 |
MnI | 2 | 4030.75 | 24802 | 0 | 1.7+7 |
MgI | 1 | 4571.10 | 21870 | 0 | 2.2+2 |
3 | 3829.32 | 47957 | 21850 | 8.9+7 | |
3 | 3832.35 | 47957 | 21870 | 6.7+7 | |
3 | 3838.29 | 47957 | 21911 | 4.5+6 | |
FeI | 2 | 4461.65 | 23111 | 704 | 2.9+4 |
2 | 4375.93 | 22846 | 0 | 2.9+4 | |
3 | 4206.70![]() |
24181 | 416 | 7.2+3 | |
3 | 4216.18 | 23711 | 0 | 1.8+4 | |
3 | 4291.46![]() |
23711 | 416 | 4.1+3 | |
42 | 4307.90 | 35768 | 12560 | 3.4+7 | |
42 | 4202.03 | 35768 | 11976 | 8.2+6 | |
73 | 3852.57 | 32499 | 17550 | 2.9+6 | |
648 | 4374.49![]() |
49477 | 26624 | 4.9+5 | |
828 | 4427.30 | 22997 | 416 | 3.4+4 | |
FeII | 38 | 4583.84 | 44449 | 22637 | 3.8+5 |
6F | 4457.95 | 22810 | 385 | 2.9-1 | |
7F | 4359.33 | 23318 | 385 | 1.1+0 | |
7F | 4287.39 | 23318 | 0 | 1.5+0 | |
21F | 4276.83 | 25805 | 2430 | 6.5-1 | |
21F | 4243.97 | 25429 | 1872 | 9.0-1 |
For data reduction, the Image Reduction and Analysis Facility (IRAF) of the National Optical Astronomy Observatories was used. The wavelength calibration of the spectra was obtained from thorium-argon arcs taken during each run. The standard stars HR 718, HR 3454 or HR 7596 were observed on photometric nights to allow calibration of the spectra to absolute fluxes. The calibrated fluxes from multiple observations of standard stars during each night were compared to estimate the errors in the fluxes tabulated in Tables 4-8: these comparisons indicate that the flux errors should be less than 5%. This error is much less than the typical variation of flux throughout the pulsation cycle (see Sect. 4).
In the case of non-photometic nights, we had to estimate the absolute fluxes
from the spectra by assuming that the flux in the V band at 5100Å (a
quasi-continuum point) varies in the same way as the light curves of The
American Association Of Variable Star Observers
(AAVSO) (Mattei et al.
1980). With this method we were able to estimate the continuum flux at
5100Å by scaling the non-photometric flux relative to the calibrated flux
from a photometric night. We were then able to calibrate the flux at other
wavelengths using the relative spectral response of the spectrograph. The line
fluxes obtained in this way are marked with a colon in Tables 4-8.
In order to derive the emission line fluxes, it is necessary to have an estimate of the continuum level. This was estimated by eye: initially, the whole order containing the line was examined to get an idea of the overall continuum shape, then a wavelength interval on either side of the line approximately equal to the line width was used for the final estimate. In the case of P-Cygni line profiles, the associated absorption was ignored in the continuum estimate. For strong emission lines (ex. Balmer lines), the continuum level is not a significant contributor to line flux errors, but in the case of weak emission lines in these stars which exhibit many absorption lines (ex. FeI 4307.90Å at phase 0.63 in RLeo, Fig. 16), the continuum level can be a major contributor to the flux error. In the worst cases, we estimate that the flux could be wrong by a factor of 2.
The phases for each observation of each Mira variable were determined from the AAVSO maxima: phase zero corresponds to visible maximum. Table 3 lists the Julian dates and the corresponding phases for the variables observed.
JDa | Phase | |||||
RR Sco | R Aql | R Car | R Leo | S Scl | R Hya | |
1153 | - | - | -0.16 | - | - | - |
1169 | - | - | - | - | 0.02 | - |
1206 | - | - | -0.09 | 0.43 | - | - |
1268 | - | 0.01 | 0.11 | 0.63 | - | - |
1357 | 0.29 | 0.32 | 0.39 | - | 0.52 | 0.08 |
1375 | - | 0.39 | 0.45 | - | 0.58 | 0.12 |
1412 | 0.48 | 0.50 | 0.57 | - | 0.67 | 0.21 |
1440 | 0.58 | - | 0.66 | - | 0.75 | - |
1464 | - | 0.69 | - | - | - | - |
1542 | - | - | - | 1.49 | - | 0.53 |
Observed line profiles for the Balmer lines H,
H
,
H
and
H
are presented in Figs. 1-5 for RRSco, RAql,
RCar, SScl and RHya (because of its position in the sky and the phase of
its light curve, RLeo was not observable near maximum when the Balmer lines
were in emission, but we did observe it during two different minima). The line
![]() |
Figure 1:
Left: Line profiles of H![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 2:
Left: Line profiles of H![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 3:
Left: Line profiles of H![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
Left: Line profiles of H![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 5:
Left: Line profiles of H![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
From left to right: Line profiles of SiI 4102.94Å (M 2) as a
function of phase in RRSco, RAql, RCar, SScl and RHya
(plotted as in the left panel of Fig. 1). The feature on the left of each plot
is the edge of the nearby H![]() |
Open with DEXTER |
![]() |
Figure 7:
From left to right: Narrow (300![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 8: From left to right: Line profiles of MgI 4571.10Å (M 1) as a function of phase in RRSco, SScl, RLeo and RHya (plotted as in the left panel of Fig. 1). Table 5 lists the maximum flux level in each line |
Open with DEXTER |
![]() |
Figure 9:
Left: Line profiles of MgI 4571.10Å (M 1) and MgI
3829.32Å, 3832.35Å and 3838.29Å (all M 3) as a function of
phase in RAql (plotted as in the left
panel of Fig. 1).
Table 5 lists the maximum flux level in each line.
Right: Profiles of the MgI multiplet 3 lines observed with a narrow slit
(300![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 10: Same as Fig. 9 but for RCar |
Open with DEXTER |
![]() |
Figure 11: Line profiles of FeI 4427.30Å (M 828), 4375.93Å (M 2), 4307.90Å, 4202.03Å (both M 42) and 3852.57Å (M 73) as a function of phase in RRSco (plotted as in the left panel of Fig. 1). Tables 6 and 7 list the maximum flux level in each line |
Open with DEXTER |
![]() |
Figure 12: Same as Fig. 11 but for R Car, and with FeI 4461.65Å (M 2) replacing FeI 4427.30Å (M 828) in the left panel |
Open with DEXTER |
![]() |
Figure 13:
Left: Line profiles of FeI 4307.90Å, 4202.03Å (both M 42)
and 3852.57Å (M 73) as a function of phase in SScl (plotted as in the left
panel of Fig. 1). Tables 6 and 7
list the maximum flux level in each line. Right: The FeI 4202.03Å line
profile observed with a narrow (300![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 14: Left: Same as Fig. 13 but for RHya, and including FeI 4307.90Å (M 42) and 3852.57Å (M 73) in the right panel |
Open with DEXTER |
![]() |
Figure 15: Left: Same as Fig. 13 but for RAql, and including FeI 3852.57Å (M 73) in the right panel |
Open with DEXTER |
![]() |
Figure 16:
Left: Line profiles of FeI 4307.90Å (M 42) as a function of
phase in RLeo. Right: Line profiless of FeI 4461.65Å,
4375.93Å (both M 2), 4427.30Å (M 828), 4202.03Å (M 42) and
4216.18Å (M 3) at phase ![]() |
Open with DEXTER |
![]() |
Figure 17:
Left: Line profile of FeII 4583.84Å (M 38) at phase
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 18:
Left: Line profiles of MnI 4030.75Å (M 2), FeII
4583.84Å (M 38) and [FeII]6F 4457.95Å at phase ![]() ![]() ![]() |
Open with DEXTER |
profiles are plotted against velocity rather than wavelength,
because we want
to analyse the atmospheric kinematics associated with pulsation. Zero velocity
in the plots corresponds to the center-of-mass velocity of the star and a
positive velocity indicates a motion outward from the center of the star, and
visa versa. Most of the observations in the plots (left panels) were taken
with a wide slit, but we also present selected observations taken with a narrow
slit (right panels) to give a comparison. Clearly, since the Balmer lines have
intrinsic full widths of 80 km s-1, the instrumental broadening of
the wide slit (15.3 km s-1) does not influence the line profiles
significantly.
The main result seen in the presented spectra is the change in the
overall line shape with phase. The full base width (hereinafter FBW) of the
hydrogen lines is around 80 km s-1 when they first appear. With increasing phase,
the hydrogen lines become weaker and more narrow (FBW60-70 km s-1).
The emission lines are generally centered around
an outward velocity of
10 km s-1.
Previous studies of the hydrogen
lines in Miras (e.g. Fox et al. 1984) have shown similar results.
A noticeable feature of the hydrogen lines is that the strengths of
superimposed absorption features diminish with advancing phase. Joy
(1947) identified absorption lines of TiO in H,
of FeI and
VI in H
and of TiI in H
.
A decrease in absorption with
increasing phase is just as expected in a model where the shock producing the
emission lines moves above more of the atmosphere as time progresses. The
strength of the absorption features also seems to vary from star to star; for
instance, superimposed absorption in H
is stronger in RAql than in SScl at
comparable phases (Fig. 2
and
Fig. 4
). Fox et al.
(1984) noted a similar variation of overlying absorption line strength
from star to star and suggested that this could be due to a cycle-to-cycle
variation in shock strength: at fainter maxima weaker shocks are formed and
they are not able to propagate out through as much overlying absorbing
material, leading to stronger absorption.
A remarkable result is the appearance of faint Balmer lines at the very early
phase
in RRSco (Fig. 1). These emission lines
are dominated by strong overlying absorption, as expected. They could also be
evidence for an unusually strong shock emerging from the star in the observed
pulsation cycle.
Various metal emission lines can be observed in Mira variables around maximum light and during the post-maximum phases. In the spectral range of our observation, we observed the following emission lines, which are ordered by the mutiplet number M: SiI (M2: 4102.94Å), MgI (M1: 4571.10Å), MgI (M3: 3829.32Å, 3832.35Å and 3838.29Å), FeI (M2: 4461.65Å and 4375.93Å), FeI (M3: 4216.18Å), FeI (M42: 4307.90Å and 4202.03Å), FeI (M73: 3852.57Å), FeI (M828: 4427.30Å), FeII (M38: 4583.84Å) as well as forbidden lines [FeII]6F (4457.95Å), [FeII]7F (4359.33Å and 4287.39Å) and [FeII]21F (4276.83Å and 4243.97Å) and finially MnI (M2: 4030.75Å). We discuss these emission lines in the order of their appearance during the pulsation cycle.
The SiI line of the multiplet 2 at 4102.94Å appears at the same phases
as the hydrogen lines, namely around maximum light, and varies in strength in
the same way as the Balmer lines. The variation in appearance with phase is presented
in Fig. 6 for RRSco, RAql, RCar, SScl and RHya. For
comparison of wide and narrow slit spectra, additional observations, taken with
a narrow slit, are shown in Fig. 7 for one phase in each of RCar
and SScl. We measured a FBW of the silicon line of around 60 km s-1when it first appears, somewhat less than the FBW of the Balmer lines. As with
the Balmer emission lines, the SiI emission line can be seen at the
remarkably early phase
in RRSco (Fig. 6, left panel).
The observed MgI emission lines (M1: 4571.10Å; M3: 3829.32Å, 3832.35Å and 3838.29Å) are presented in Figs. 8-10. Additional narrow slit observations are shown in the same figures for the multiplet 3 lines, and in Fig. 7 for the 4571.10Å line. Note that in some of the spectra, the 4571.10Å line appears totally in absorption.
The lines of the MgI multiplet 3 at 3829.32Å, 3832.35Å and
3838.29Å appear around the phase of maximum when the hydrogen lines are
dominant. These lines are high excitation lines (lower state excitation
potential 2.7 eV) and the Einstein-coefficient is rather large (see
Table 2). The emission lines were detected in RAql and RCar
(Figs. 9 and 10) where we observed the early
post-maximum phases in detail. We measured a FBW of
30-40 km s-1. It is obvious that these lines are always
blueshifted. For RAql (Fig. 9) at phase
,
the emission
lines show equal velocity shifts of
+10 km s-1, as can also be seen
in the panel of the narrow slit observation in Fig. 9 (all emission
line velocities quoted in this paper are measured at half maximum height). For
RCar (Fig. 10), the velocity shift of the lines at 3829.32Å
and 3832.35Å is about +20 km s-1 when they first appear
(
): the weak 3838.29Å line appears to be strongly affected by
overlying absorption causing it to show a smaller velocity shift of only
+10 km s-1. At phase
,
all multiplet 3 lines are equally
blueshifted by
+10 km s-1. For RHya, we could not detect the
multiplet 3 MgI lines at phases
and 0.12 although one would
expect them to be detectable.
At later phases, the lines of the multiplet 3 disappear and the
MgI 4571.10Å (M1) emission line appears. The lower level of this
emission line is the ground state and the Einstein-coefficient is rather low
(see Table 2). Because it is a ground state line, it
is prominent in absorption as well as emission. Before and around maximum
light, two absorption components appear, especially in RCar
(Fig. 10,
-0.11), presumably
corresponding to material behind, and in front of, the emerging shock front.
The majority of the observed MgI 4571.10Å
emission lines are affected
by obvious overlying absorption
on the negative velocity side, although this absorption decreases as the
phase advances.
The FBW of the MgI 4571.10Å emission line is generally around
40-50 km s-1 when it first appears and it is blueshifted by
+10 km s-1 with respect to the center-of-mass velocity. There is a
general decrease in the blueshift towards minimum light; this is probably due
to a combination of a real slowing of the shock and a decrease in overlying
absorption on the red edge. In RAql (Fig. 9), the velocity shift around
the minimum phase (
)
seems to increase again, possibly due to
some asymmetry in the upper atmosphere being transversed by the shock.
RLeo is an interesting case to study, since there are two minima observed
which appear rather different. Around the first minimum, the MgI 4571Å
line was very faint and narrow ( km s-1) with the velocity
shift decreasing from +5 km s-1 to +2 km s-1 for the phases
to
.
At the next minimum, the emission line was bright
and showed a FBW of
40 km s-1 and the line was centered at a
velocity shift of
+3 km s-1. The AAVSO light curve shows that the
light maximum preceding the first minimum was much fainter than the maximum
preceding the second minimum. A bright preceding maximum obviously leads to a
much wider and stronger MgI 4571.10Å emission line. This is consistent
with the hypothesis that a stronger (visually brighter) maximum is associated
with a stronger shock.
The development of the MgI 4571Å line in RHya (Fig. 8, right panel), the star with the
longest period (389 d) of our sample, is obviously different. The emission
line appears at a relatively early phase ()
and it has a sharp low
velocity edge indicative of strong overlying absorption. Clearly, velocities
derived from this line need to be treated with caution, at least until very
advanced phases.
Out of the variety of neutral iron emission lines in Mira variables, we studied the following: (M2) 4461.65Å and 4375.93Å, (M3) 4216.18Å, (M42) 4307.90Å and 4202.03Å, (M73) 3852.57Å and (M828) 4427.30Å. Figs. 11-16 show the observed FeI lines for all stars of our sample. For SScl, RHya and RAql we present narrow slit as well as broad slit observations for some of the neutral iron lines.
It is obvious that several of the FeI emission lines, namely FeI 4307.90Å, 4202.03Å and 3852.57Å, appear brightly during the post-maximum phase of every star of our sample. Table 2 shows that these are all high excitation lines (lower level excitation potential 1.5-2.2 eV) with large Einstein-coefficents. It has been suggested that the great strength of all three lines is the result of fluorescence, the lines being pumped by the MgII 2795, 2802Å lines (Thackeray 1937; Willson 1976).
The 3852.57Å line appears around maximum visible light with a very
narrow FBW of 25-30 km s-1. At later phases it widens up to a
FBW of
40-60 km s-1. It appears that this is due to
overlying absorption at early phases removing the negative velocity side of the
line (see RCar, RHya and RAql in Figs. 12, 14 and
15).
The 4202.03Å line also appears near maximum light (see RCar, SScl and RAql in Figs. 12, 13 and 15). At early phases, it shows an obvious inverted P-Cygni profile caused by overlying absorption of infalling material.
The 4307.90Å line does not appear as prominently as the
previous line for phases
0.2 even though both lines have identical upper
states (see Table 2). This is almost certainly due to more overlying
absorption affecting this line: around maximum light in RCar,
SScl and RAql (Figs. 12, 13 and 15), the
line is totally in absorption while the 4202.03Å line shows distinct
emission. Multiple absorption components are clearly visible cutting into the
emission line except at the latest phases. This may be
partly self-absorption by pre- and post-shock material as the 4307.90Å
has the largest Einstein-coefficent of the two lines. The large
Einstein-coefficent may also explain the fact that the 4307.90Å line
persists the longest and has the largest flux (see Table 6) of
the two lines.
All other FeI emission lines observed are low excitation lines (see
Table 2).
The 4375.93Å line of the multiplet 2 appears around the phase of maximum
light (see RRSco, RCar and RLeo in Figs. 11, 12,
16), which is an exceptional early appearance for a low
excitation line. The line shows an obvious inverted P Cygni profile with a
steep red wing: the absorption component decreases towards later phases.
The 4461.65Å line was only detected in emission in RCar at the phase
and in RLeo at
(Figs. 12 and
16). Although it belongs to the same multiplet (2) as the 4375.93Å line,
it is not as prominent at early phases as the latter.
The remainder of the observed FeI low excitation lines only appear in some stars
and no regularity can be seen. The FeI 4427Å emission line was observed
in RRSco and RLeo at just one phase each (see Figs. 11 and
16). Finally, the low excitation line at 4216.18Å was observed in
RLeo at the phase .
The permitted FeII 4583.84Å emission line was observed in three of the
stars of our sample, namely RRSco, RCar and RLeo
(Figs. 17 and 18). This is a high excitation
line (lower level excitation potential 2.8 eV) which appears towards the
minimum light of the star, at phase
in RRSco,
in
RCar and
in RLeo. The FBW varies from
35 km s-1 in
RRSco and RCar to
45 km s-1 in RLeo. The velocity shifts of
the FeII line are around +2 km s-1 in these stars.
In RLeo and RCar, forbidden [FeII] lines can be clearly seen at the same
phase as the permitted FeII line at 4583.84Å
(Figs. 17 and 18). In RLeo at phase
,
following the bright maximum, we found the transitions
[FeII]6F (4457.95Å), [FeII]21F (4276.83Å and 4243.97Å)
and [FeII]7F (4359.33Å and 4287.39Å). No forbidden lines could
be seen in the previous cycle following a fainter maximum: the cycle-to-cycle
variability of the forbidden lines is consistent with the suggestion noted
earlier that stronger shocks are associated with brighter maxima.
In RCar at phase ,
we identified only the two [FeII]7F
emission lines at 4359.33Å and 4287.39Å. Unfortunately, the
observation for RCar at phase
was taken with a narrow slit in
poor weather conditions, so the signal to noise of these faint lines is low.
Other forbidden lines were probably present but undetectable. In RRSco we
could not identify any of the forbidden lines, although we could clearly detect
the permitted FeII line (Fig. 17).
The FBW of the [FeII]7F emission lines in RCar is
21-26 km s-1. The velocity shift was measured at 0-2 km s-1,
which is the same as for the permitted FeII 4583.84Å emission. In
RLeo, all the [FeII] emission lines have a FBW of
40 km s-1.
The velocity shifts are approximately 0-2 km s-1 for the 7F and
the 21F emission lines, similar to the velocity shift of permitted
FeII 4583.84Å line.
The manganese line at 4030.75Å is a ground state line. It was observed
in just one star, namely RLeo, at phase (Fig. 18). The MnI emission line shows a FBW of
28 km s-1 and has a large velocity shift of 10 km s-1, unlike
the FeII and [FeII] lines at the same phase. As with the neutral Fe lines
of low excitation (e.g. FeI 4375.93Å), the sharper red edge of the line
indicates absorption by infalling, neutral Mn atoms in the layers above the
shock.
Total line fluxes
(ergcm-2s-1) for the observed emission lines
are given in Tables 4 to 8, along with the peak
flux
(ergcm-2s-1Å-1) in the line plots.
The peak flux can be used in conjunction with Figs. 1-18
(except for the narrow slit observations) to get intensities through the line
profiles. Note that no correction for superimposed absorption lines has been
made in calculating total line fluxes.
H![]() |
H![]() |
H![]() |
H![]() |
SiI ( M2; 4102Å) | |||||||
Star | Phase |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
RR Sco |
0.29 | 8.2-12 | 3.2-12 | 1.1-11 | 3.7-12 | 5.3-12 | 1.6-12 | 4.3-12 | 1.5-12 | 8.5-13 | 7.2-13 |
0.58 | 5.0-13 | 2.7-13 | 3.5-12 | 1.2-12 | - | - | - | - | 1.2-12 | 6.8-13 | |
R Aql | 0.01 | 1.4-10 | 8.3-11 | 3.5-10 | 1.6-10 | 1.2-10 | 5.8-11 | 1.0-10 | 5.2-11 | 1.9-11 | 1.1-11 |
0.32 | 3.0-12 | 1.6-12 | 3.2-12 | 1.0-12 | 1.1-12 | 4.9-13 | 1.2-12 | 2.8-13 | - | - | |
0.39 | 1.5-12: | 1.0-12: | 1.9-12: | 6.2-13: | - | - | - | - | - | - | |
R Car | -0.16 | 6.1-11: | 2.0-11: | 1.6-10: | 4.7-11: | - | - | - | - | 1.9-11: | 6.9-12: |
-0.09 | 1.1-9: | 7.8-10: | 1.8-9: | 9.7-10: | 3.1-10: | 1.9-10: | 1.9-10: | 1.2-10: | 1.2-10: | 5.2-11: | |
0.11 | 6.8-10 | 2.7-10 | 1.0-9 | 3.7-10 | 5.0-10 | 1.9-10 | 4.1-10 | 1.6-10 | 4.3-11 | 1.2-11 | |
0.39 | 3.7-12: | 2.2-12: | - | - | - | - | - | - | - | - | |
S Scl | 0.02 | 7.0-11 | 3.7-11 | 3.3-10 | 1.4-10 | 8.0-11 | 4.1-11 | 5.0-11 | 3.5-11 | 1.4-11 | 6.8-12 |
R Hya | 0.08 | 3.0-10: | 1.3-10: | 2.7-10: | 1.0-10: | 1.7-10: | 6.0-11: | 8.8-11: | 3.0-11: | 1.8-11: | 6.8-12 |
0.12 | 2.1-10: | 9.1-11: | 1.1-10: | 4.3-11: | 9.5-11: | 3.4-11: | 3.2-11: | 1.0-11: | 2.2-11: | 4.9-12: | |
0.21 | 3.2-11 | 1.7-11 | 2.2-11 | 8.5-12 | 1.3-11 | 4.3-12 | 1.2-11 | 5.4-12 | 2.4-12 | 9.8-13 |
MgI ( M1; 4571Å) | MgI ( M3; 3829Å) | MgI ( M3; 3832Å) | MgI ( M3; 3838Å) | MnI ( M2; 4030Å) | |||||||
Star | Phase |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
RR Sco |
0.29 | 1.3-12 | 4.5-13 | - | - | - | - | - | - | - | - |
0.48 | 2.9-13 | 6.6-14 | - | - | - | - | - | - | - | - | |
R Aql | 0.01 | - | - | 1.1-11 | 2.8-12 | 1.4-11 | 3.4-12 | 1.1-11 | 2.8-12 | - | - |
0.32 | 2.0-12 | 4.8-13 | - | - | - | - | - | - | - | - | |
0.39 | 2.6-12: | 7.5-13: | - | - | - | - | - | - | - | - | |
0.50 | 6.5-13 | 1.7-13 | - | - | - | - | - | - | - | - | |
R Car | -0.09 | - | - | 3.6-11: | 1.0-11: | 3.4-11: | 5.7-12: | 1.9-11: | 3.7-12: | - | - |
0.11 | - | - | 2.7-11 | 9.6-12 | 3.2-11 | 7.2-12 | 3.0-11 | 1.5-12 | - | - | |
0.39 | 1.8-11: | 3.7-12: | - | - | - | - | - | - | - | - | |
0.45 | 6.7-12: | 1.5-12: | - | - | - | - | - | - | - | - | |
0.57 | 1.0-12 | 3.1-13 | - | - | - | - | - | - | - | - | |
0.66 | 4.0-13 | 1.4-13 | - | - | - | - | - | - | - | - | |
R Leo | 0.43 | 6.6-12: | 1.4-12: | - | - | - | - | - | - | - | - |
0.63 | 9.5-13 | 1.8-13 | - | - | - | - | - | - | - | - | |
1.49 | 7.4-12 | 2.0-12 | - | - | - | - | - | - | 8.8-13 | 2.6-13 | |
S Scl | 0.52 | 2.6-13 | 6.8-14 | - | - | - | - | - | - | - | - |
0.58 | 1.5-13 | 4.9-14 | - | - | - | - | - | - | - | - | |
R Hya | 0.12 | 9.4-12: | 3.2-12: | - | - | - | - | - | - | - | - |
0.21 | 8.4-12 | 3.1-12 | - | - | - | - | - | - | - | - | |
0.53 | 1.4-12: | 4.7-13: | - | - | - | - | - | - | - | - |
FeI ( M2; 4461Å) | FeI ( M2; 4375Å) | FeI ( M3; 4216Å) | FeI ( M42; 4307Å) | FeI ( M42; 4202Å) | |||||||
Star | Phase |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
RR Sco |
0.29 | - | - | 3.6-13 | 1.2-13 | - | - | 1.4-12 | 7.1-13 | 3.9-12 | 7.4-13 |
0.48 | - | - | - | - | - | - | 1.3-13 | 4.7-14 | - | - | |
R Aql | 0.01 | - | - | - | - | - | - | - | - | 9.9-12 | 2.3-12 |
0.32 | - | - | - | - | - | - | 9.7-13 | 5.0-13 | 2.2-12 | 5.1-13 | |
0.39 | - | - | - | - | - | - | 1.4-12: | 5.6-13: | 2.0-12: | 4.3-13: | |
0.50 | - | - | - | - | - | - | 3.5-13 | 1.1-13 | - | - | |
R Car | -0.09 | - | - | 1.6-11: | 1.6-12: | - | - | - | - | 2.9-11: | 6.5-12: |
0.11 | - | - | 2.4-11 | 4.4-12 | - | - | - | - | 4.5-11 | 1.0-11 | |
0.39 | 1.6-12: | 2.7-13: | 1.4-12: | 1.7-13: | - | - | 1.2-11: | 5.6-12: | 9.5-12: | 3.5-12: | |
0.45 | - | - | - | - | - | - | 2.1-12: | 7.6-13: | 1.8-12: | 6.2-13: | |
0.57 | - | - | - | - | - | - | 7.3-13 | 2.1-13 | - | - | |
0.66 | - | - | - | - | - | - | 2.6-13 | 6.6-14 | - | - | |
R Leo | 0.43 | - | - | - | - | - | - | 1.5-12: | 2.8-13: | - | - |
0.63 | - | - | - | - | - | - | 4.1-13 | 1.0-13 | - | - | |
1.49 | 3.6-13 | 7.9-14 | 4.0-13 | 8.4-14 | 2.5-13 | 5.3-14 | 3.7-12 | 1.5-12 | 2.6-12 | 1.0-12 | |
S Scl | 0.02 | - | - | - | - | - | - | - | - | 6.5-12 | 1.7-12 |
0.52 | - | - | - | - | - | - | 2.6-13 | 8.4-14 | - | - | |
0.58 | - | - | - | - | - | - | 1.3-13 | 3.9-14 | - | - | |
R Hya | 0.08 | - | - | - | - | - | - | 8.0-12: | 3.1-12: | 4.3-11: | 1.1-11: |
0.12 | - | - | - | - | - | - | 9.8-12: | 2.8-12: | 1.9-11: | 5.2-12: | |
0.21 | - | - | - | - | - | - | 5.5-12 | 2.5-12 | 1.1-11 | 2.9-12 | |
0.53 | - | - | - | - | - | - | 2.5-12: | 7.8-13: | - | - |
FeI ( M73; 3852Å) | FeI ( M828; 4427Å) | FeII ( M38; 4583Å) | |||||
Star | Phase |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
RR Sco |
0.29 | 4.8-12 | 1.2-12 | 4.2-13 | 9.2-14 | 2.3-13 | 4.7-14 |
R Aql | 0.01 | 1.5-11 | 3.5-12 | - | - | - | - |
0.32 | 3.2-12 | 9.7-13 | - | - | - | - | |
R Car | -0.09 | 2.7-11: | 3.9-12: | - | - | - | - |
0.11 | 1.1-10 | 2.9-11 | - | - | - | - | |
0.39 | - | - | - | - | 7.7-13: | 5.8-13: | |
R Leo | 1.49 | - | - | 4.2-13 | 9.1-14 | 2.4-13 | 7.3-14 |
S Scl | 0.02 | 8.5-12 | 1.9-12 | - | - | - | - |
R Hya | 0.08 | 3.6-11: | 6.8-12: | - | - | - | - |
0.12 | 1.9-11: | 4.2-12: | - | - | - | - | |
0.21 | 8.3-12 | 3.7-12 | - | - | - | - |
[FeII]6F | [FeII]7F | [FeII]7F | [FeII]21F | [FeII]21F | |||||||
( 4457Å) | ( 4359Å) | ( 4287Å) | ( 4276Å) | ( 4243Å) | |||||||
Star | Phase |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
R Car |
0.39 | - | - | 1.2-12: | 2.5-13: | 2.7-12: | 5.5-13: | - | - | - | - |
R Leo | 1.49 | 2.3-13 | 9.0-14 | 5.4-13 | 1.4-13 | 7.8-13 | 3.0-13 | 2.6-13 | 7.9-14 | 5.4-13 | 1.4-13 |
![]() |
Figure 19:
Left: Absolute line fluxes plotted against phase in RLeo. Note
the change in scale of the phase for we plotted the fluxes of the minimum of
the next cycle in the same panel. Note also that several line fluxes are coincident at
phase ![]() ![]() ![]() |
Open with DEXTER |
Figures 19-22 show the variation with phase
of all fluxes measured in each star of our sample. For the Balmer lines, both
the data given in this paper and the data of Fox et al. (1984) are used
in the plots. The typical variation with phase can be seen best in RCar
(Fig. 22), for which a separate plot is shown for the Balmer
lines. The Balmer lines H
and H
appear first, roughly between
phases -0.4 and -0.2 (see the Fox et al. 1984 data for RCar, SScl,
RRSco and RAql as well as the present data for RCar - Figs. 20,
21 and 22).
The lines then strengthen rapidly, reaching maximum flux around phase -0.1. At this phase,
the other Balmer lines are observed too although they tend to
reach maximum flux around phase +0.1 when the amount of overlying absorption is
reduced. The Balmer lines finally vanish about phase
.
Fox et al. (1984) noted that the Balmer line fluxes can vary enormously
from cycle to cycle (see o Ceti in their data), with brighter maxima yielding
brighter Balmer lines. A comparison of their fluxes with ours for RCar also
shows this effect. The Balmer lines we observed in RCar are roughly ten
times brighter than the lines observed by Fox et al: at the same time, the AAVSO
light curves show that the visual maximum during our observations (
)
was about one magnitude brighter than the maximum (
)
during
which Fox et al. made their observations.
The SiI 4102Å (M2) high
excitation line appears in all stars near maximum light. The
line flux is always less than in the nearby Balmer line H.
In general,
the variation with phase and the phase of the peak flux is
similar to that of the hydrogen lines, although the peak flux is a factor of
10 less than in the hydrogen lines.
The high excitation lines of multiplet 3 of MgI (3829.32Å, 3832.35Å and
3838.29Å) appear near maximum light and were observed in RAql and
RCar (Fig. 21 right panel, and Fig. 22
left panel). The flux in these lines is about a factor of 30 less than in the
Balmer lines. The lines disappear at
0.2-0.3 and are
replaced by the low excitation line of MgI (M1) at 4571.10Å. The
latter appears around
0.2-0.3 and was observed in
every star of our sample, disappearing finally at the late phase
0.6. It appears exceptionally early (
)
in RHya, as already noted in Sect. 3.2.2.
The FeI lines show a rather complex variation with phase. In general, one can
divide them into two groups: type A, those appearing early in the cycle (
maximum); and type B, those appearing around
0.3-0.4. The first group consists of 4202.03Å (M42), 3852.57Å
(M73) (both high excitation lines) and the 4375.93Å (M2) (a low
excitation line). Variation of flux with phase in RHya, RAql and RCar
can be seen Fig. 19 right panel,
Fig. 21 right panel, and Fig. 22 left
panel). These lines vary less rapidly than the Balmer or SiI lines and reach
their maximum at the later phase of
0.1.
The type B FeI emission lines 4427.29Å (M828), 4461.65Å (M2),
4216.18Å (M3) and 4307.90Å (M42) appear around
0.3-0.4 and are usually fainter than FeI emission
lines discussed above. Apart from the 4307.90Å line, these emission
lines are all of low excitation. As noted in Sect. 3.2.3, the
4307.90Å line has the same upper level as the 4202.03Å line of the
same multiplet and its late appearance and low flux relative to the
4202.03Å line is almost certainly due to overlying
absorption.
However, the flux in the 4307.90Å line at late phases is
usually greater than in other late emission lines of FeI 4427.29Å,
4461.65Å and 4216.18Å. Its flux and variation with phase is quite
similar to the MgI 4571.10Å line.
FeII lines were observed at late phases in a number of stars. Line fluxes of
the permitted FeII 4583.84Å (M38) emission line are shown for RRSco,
RCar and RLeo in Figs. 21, 22
and 19 (left panels). We note that this high excitation
line appears around 0.3-0.5 and is one of the
faintest lines observed.
Fluxes for the forbidden lines [FeII]7F at 4359.33Å and
4287.40Å are plotted in Fig. 22
(left panel) for RCar at phase
and in
Fig. 19 (left panel) for RLeo at
.
The
line flux in RCar is similar to that of the
permitted FeII emission line at 4583.84Å. In RLeo, however, the
forbidden line flux is
3 times higher than the permitted line flux, at
least for
where the forbidden lines were observed. At phase
in RLeo, we also measured the line fluxes of the [FeII]6F and
[FeII]21F emission lines. These lines are weaker than the [FeII]7F
lines: their flux is comparable to the faint, late-appearing FeI lines.
It is worth noting the much higher metal emission line fluxes during the second
observed minimum of RLeo. Comparing the line fluxes of MgI 4571.10Å
and FeI 4307.90Å, we see that the MgI line flux is stronger by a
factor of 3 in the second cycle whereas the FeI line flux is stronger
by a factor of almost 10. Finally, in Fig. 23 we summarize the
overall duration and flux variation with phase of emission lines in M-type
Miras.
Another property of the emission line spectra that is useful for comparison with theoretical models of Mira atmospheres is the velocity of the emission lines. This line velocity should represent the velocity of the emission region associated with the shock wave passing through the Mira atmosphere.
![]() |
Figure 20:
Absolute line fluxes plotted against phase in SScl. Note the nearly
equal fluxes for H![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
In Fig. 24, we plot the velocities (relative the stellar
center-of-mass) of representative emission lines against the phase of the
pulsation cycle. The velocities are defined at half-height of the line
profiles and are obtained from all stars in the sample. When overlying
absorption is dominant (for example, at early phases of the H
line), no
velocity measurement was made.
The plot clearly shows that when the shock emerges from deep in the
photosphere, the post-shock emission region has a measured outward velocity of
10-12 km s-1. This result is in good agreement with the
velocity measured in the infrared (Hinkle 1978; Hinkle et al.
1982, 1984) for deep pulsating layers when converted to
center-of-mass velocities (Wood 1987). As phase advances, the emission
line velocity decreases until it becomes essentially zero around minimum light.
Note that the material whose velocity is measured by the emission lines is
associated with the near-shock zone: this is quite different from the deeper,
infall zone whose velocity is measured near minimum light by the infrared
spectra.
![]() |
Figure 21:
Left: Absolute line fluxes plotted against phase in
RRSco. Note that at phase ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 22:
Left: Absolute line fluxes plotted against phase in RCar.
Note that there are coincident points for the phases
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
At first sight, it is surprising that in the average Mira the apparent velocity
of the emission lines approaches zero half a cycle after the shock emerges from
the deep photosphere. Since the lines are still in emission, the shock must
still be propagating outward so that the post-shock material, from which the
line emission presumably originates, should show a positive outward velocity.
(We note in passing that since pulsation in the outer layers can be quite
irregular (e.g. Bessell et al. 1996), in an individual cycle the shock
could be stalled or even reversed by infalling material from a previous cycle,
but, on average, shocks must progress outward.) By the time of minimum light,
the shock is far above the photosphere (defined for the present discussion to
be at optical depth one). For example, in a typical fundamental mode Mira
model (e.g. see the models of Bessell et al. 1996 or Hofmann et al.
1998), the photosphere is at 240
while the shock is at
420
.
In a simple geometric model for emission from such a
system, neglecting absorption above the photosphere, the emission lines would
have square profiles going from +
to
,
where
would be
5 km s-1: only the emission with the most
negative velocity is hidden behind the star. We would therefore expect the
line emission to be centered close to velocity zero, as observed. Detailed
models for the transfer of line photons originating from the shock are needed
to make quantitative estimates of shock velocities high in the Mira atmosphere.
![]() |
Figure 23: Simplified time sequence for the observed emission lines and their total fluxes (in arbitary units) |
Open with DEXTER |
![]() |
Figure 24: The velocities of representative emission lines plotted against the phase of the pulsation cycle |
Open with DEXTER |
A sample of six M-type Miras has been observed with a phase coverage of about
one pulsation period. We obtained time-resolved, high-dispersion emission line
spectra and line fluxes in the optical spectral region for the Balmer lines
(H,
H
,
H
,
H
)
and for metal lines of MgI, MnI,
SiI, FeI, FeII and [FeII]. The variation of line shape, velocity and
flux with phase has been discussed.
Evidence was presented that brighter maxima are associated with stronger shock
waves. In RCar, the Balmer line fluxes we observed around maximum were much
stronger than the fluxes observed by Fox et al. (1984), and the maximum
during which we observed the Balmer lines was 1 mag. brighter than the
maximum during which Fox et al. observed. In RLeo, the second minimum during
which we observed showed much stronger metal emission lines than the first, and
the second minimum followed a brighter maximum than the first. Forbidden
[FeII] lines were only observed during the second minimum of RLeo and after
the relatively bright maximum of RCar. This suggests that their erratic
appearance may be because they require an exceptionally bright maximum for
production.
A proper interpretation of the observations presented in this paper will require the help of detailed theoretical modelling. Our future work will involve self-consistent shock models for the production of the observed emission lines and radiative transfer calculations in moving atmospheres under conditions of non-LTE, to model the transfer of the emission line photons through the Mira atmosphere.
Acknowledgements
This work was supported by the Deutscher Akademischer Austauschdienst, DAAD with a HSPIII grant in 1999, project number D/99/15656. In 2000 this work was supported by the Deutsche Forschungsgemeinschaft, DFG project number SE 420/21-1.