A&A 416, 1107-1115 (2004)
DOI: 10.1051/0004-6361:20031385
T. Kipper 1 - V. G. Klochkova2 - K. Annuk1 - A. Hirv1 - I. Kolka1 - L. Leedjärv1 - A. Puss1 - P. Skoda3 - M. Slechta3
1 - Tartu Observatory, Tõravere 61602, Estonia
2 -
Special Astrophysical Observatory, Nizhnij Arkhyz 369167,
Russia
Isaac Newton Institute of Chile, SAO Branch
3 -
Astronomical Institute CAS, Ondrejov 25165, Czech Republic
Received 6 June 2003 / Accepted 25 August 2003
Abstract
Spectroscopic observations of the peculiar variable V838 Mon
during the period from the second light outburst until the fast dimming
are presented. We describe high resolution
high S/Nspectra obtained a day before the second light maximum and low resolution
spectra covering the whole period. The temporal
run of intensities and radial velocities of various lines
is presented. Using Na I D
IS lines we determine the reddening distance of V838 Mon d>3.1 kpc, and
kinematic distance d>4 kpc. We estimate that V838 Mon is slightly metal
deficient but otherwise has a quite solar-like chemical composition except
for enhanced abundances of Li, Ba and La.
Key words: stars: variables: general - stars: novae, cataclysmic variables - stars: individual: V838 Mon
The peculiar object V838 Mon (RA 704
04
9, Dec -3$^$50
51
1 (2000)) was
photographically discovered by Brown (Brown et al. 2002) on
January 06, 2002 when its brightness was
.
Brown noted that nothing appeared at the same position on
December 22, 2001. The precise position
of the object was measured by Oksanen &
Solonen (Brown et al. 2002). Some
authors have estimated the brightness of the
progenitor. The results are quite different. Munari et al.
(2002a) found using
SERC plates B=16.10 and
.
Kimeswenger et al.
(2002)
estimated R=14.56 and B=15.87 using SuperCOSMOS frames. Goranskii et al.
(2002) found
using the Digitized Sky Survey and SAI Crimean Station plates
and
.
V838 Mon first showed a nova-like outburst with the usual decline, but this was followed by a second outburst after about one month, and then again after one month a third, somewhat weaker, one. After the third outburst the brightness declined quite rapidly and according to Henden (2002) reached V=16.04 on August 22, 2002. The compiled light curve of V838 Mon is presented in Fig. 1. The dates of our spectral observations are also indicated in this figure. The available photometry of V838 Mon has been discussed by Munari et al. (2002a) and Kimeswenger et al. (2002). Here we present and discuss our spectroscopic observations.
A set of high resolution (
)
spectra was obtained
with the Nasmyth Echelle Spectrometer
(Panchuk et al. 1999; Panchuk et al. 2002)
of Russian 6 m telescope on
February 04/05, 2002 (JD 2 452 310).
The rapidly rotating A1V star HR 4687 was also
exposed for reduction purposes. The spectra cover
without any gaps.
These spectra were reduced using the NOAO astronomical
data analysis facility IRAF. The spectral orders were extracted using
HR 4687 exposures. The continuum was placed by fitting low order
functions through
manually indicated points in every order.
The spectrum of HR 4687 was used for determining the sensitivity
function of orders. These high resolution spectra are discussed in the Sect. 3.2.
![]() |
Figure 1: V838 Mon lightcurve in V. Filled circles are the data from Kimeswenger et al. (2002), triangles and diamonds - Munari et al. (2002a). Some measurements, noted by asterisks, are supplemented from various sources (IAUC 7785, 7786, 7822, and 7889). The open circles are the quick-look data from AAVSO just for showing the fast weakening of the star. The last plotted point V=14.98 was measured by Henden (2002). The arrow indicates the time of our high resolution spectral observation and the dates of the lower resolution spectra are indicated by the vertical dashes. |
Open with DEXTER |
On the same night the low resolution (
)
observations were
started with the Cassegrain spectrograph of the Tartu Observatory 1.5 m
telescope. The spectra were recorded with the cryogenically cooled
CCD camera Orbis-1. This limited the wavelength coverage of a single
exposure to about 400 Å. During various nights different numbers of
exposures were obtained for the total wavelength region
and
,
mostly depending on
weather conditions. The nights when at least one exposure was secured
are indicated in Fig. 1. These spectra
were reduced using the ESO data
analysis package MIDAS. Sections 3.3-3.5 are based on these low
resolution spectra. For the dates
for which we have no spectra observed at the Tartu
Observatory, we used for analysing the temporal changes of Si II
,
Fe I
,
and H
lines
the spectra obtained by Buil (2002) and spectra obtained with the
Ondrejov 2 m telescope coude spectrograph (
). The
Ondrejov coude spectra were reduced using IRAF.
Two additional spectra were obtained on March 08 and 18 with the
Ondrejov 2 m telescope and echelle spectrograph HEROS (
).
These spectra
were reduced by the HEROS standard package written by O. Stahl in
Heidelberg. These spectra were obtained in two channels - blue,
covering
,
and red, covering
.
The red spectra of March 18
(JD 2 452 352) were used to
estimate the chemical composition. The blue spectra are too noisy for that
purpose.
Already the first spectra
(January 8 and 9) obtained prior to our observations were not
typical of Fe- or He/N-type classical novae early in their evolution
(Wagner et al. 2002; Della Valle & Iijima 2002).
The spectra cover
(Wagner et al. 2002)
and
(Della Valle & Iijima 2002).
These spectra showed numerous
absorption lines (Fe II, O I, Ca II) and four
prominent emission lines with P Cygni profiles. Three of these were
identified with Ba II and the fourth was the Na I D doublet.
H
was also weakly in emission. No absorption bands of TiO and C2 were seen. These features allowed Della Valle & Iijima to suggest that
the object could be a peculiar "slow'' nova or a post-AGB star in the
course of the He-flash.
The near infrared spectra (range 1.03-2.52 micrometers) obtained on January 12 by Geballe et al. (2002) showed shallow water bands at 1.4 and 1.9 micrometers, deep CO bands longward of 2.29 micrometers, hydrogen Brackett series lines and numerous strong metal lines (Mg I, Si I, Fe I) in absorption. Many of the metal lines had P Cygni profiles indicating outflow velocities of several hundred km s-1.
Zwitter & Munari (2002)
obtained the first high resolution
(
)
spectra (
)
on January 26. They found that the spectrum resembles a
heavily reddened cool K giant at heliocentric velocity +53 km s-1.
Numerous absorption lines of Ca I, Ca II, S I,
Si I, Na I, Ni I, Y II, Zr II, C I,
Cr I, Gd II, Fe I, and Ti I showed P Cygni
profiles. They did not find Ba II lines.
The average terminal
heliocentric velocity from the P Cygni absorption components was -458 km s-1, with the line core at -252 km s-1.
The H line showed narrow absorption at +54 km s-1 and weak
emission at +72 km s-1. Interstellar Na I and K I lines were observed at velocities of +37.1 and +64.9 km s-1 and
corresponded to
E(B-V)=0.56 and 0.24, respectively.
They also found a diffuse
interstellar band (DIB) at
with the equivalent width 0.29 Å. From the reddening Zwitter & Munari suggested a distance of V838 Mon greater than 3 kpc.
Iijima & Della Valle (2002) also found in the spectra obtained a week later that numerous ionized metal lines (Fe II, Ti II, Cr II, Mg II, Si II) showed a P Cygni profile. These lines emerged between Feb. 1 and 2.
Morrison et al. (2002) obtained high dispersion (
)
yellow region spectra about the same time as our spectra were observed.
They found a few lines with P Cygni profiles belonging to Fe II,
Na I, Si II, and H
.
As indicated in Fig. 1, our high
resolution spectra were obtained about a day before
V838 Mon reached its light maximum.
Almost all metallic lines visible in our spectra showed P Cygni profiles.
The lines identified belong to Na I, Mg I,
Sc II, Ti II,
Cr II, Fe I, Fe II, Sm II, Ba II.
Some of the line profiles are illustrated in Figs. 2 and 3. Faint
absorption lines (
)
of Fe I and Ca I were also found. Their identification is rather
uncertain. Some of those faint Fe I lines also show P Cygni
type emission components (Fig. 4).
![]() |
Figure 2: The profiles of some lines in the spectrum of V838 Mon on February 04, 2002 (JD 2 452 310) on a velocity scale. |
Open with DEXTER |
![]() |
Figure 3:
Profile of the H![]() |
Open with DEXTER |
![]() |
Figure 4: Faint Fe I lines at 5565.71, 5569.63 and 5572.85 Å. |
Open with DEXTER |
The metallic lines with P Cygni profiles
have several
components both in absorption and emission.
The wavelengths of the components were determined by
manualy placing the cursor when using IRAF's task "splot''.
The mean heliocentric
velocities of absorption components are
and
km s-1. In the weaker lines additional components
at velocities
and
km s-1 could be found.
In emission two components with heliocentric velocities
and
km s-1 are visible. If these emission peaks
are interpreted as belonging to approaching and receding parts of the
expanding shell, the expansion velocity would be around 35 km s-1 and
the systemic velocity
km s-1.
For the faintest (
)
Fe I and Ca I lines
only the absorption components at
km s-1
are visible. For slightly stronger lines
(
)
this is accompanied
by the absorption at
and emission at
km s-1.
As noted before, Zwitter & Munari (2002) estimated systemic velocity +53 km s-1 on January 26. Munari et al. (2002a) found for the same time using La II absorption lines the mean heliocentric velocity +54.4 km s-1. In our spectra no La II lines are visible as stronger La II fall outside of our observed spectral region.
The terminal velocity derived from strong Fe II 42 lines is -388 km s-1 and -240 km s-1 from H line. These
velocities are lower than those found by Munari et al. (2002a) from
observations two weeks earlier. Munari et al. found that the reduction
of terminal velocities continued also in February and March. These
velocities are much lower than in classical novae. Only in the slowest
classical nova HR Del and in RR Pic the absorption components of H I
and metallic lines were blue-shifted by about -600 to -700 km s-1in the early stage and about -70 to -200 km s-1 at the end of the
pre-maximum stage (Iijima et al. 1998).
The Na I D lines also show a P Cygni profile with three strong interstellar (IS) components (Fig. 11). The terminal velocity from D2line is -250 km s-1. Three separate Gaussians could be fitted to the IS components. The results are given in Table 1. Without Gaussian decomposition the shorter wavelength IS components correspond to the heliocentric velocity +36.1 km s-1. Our velocities are close to the ones found by Zwitter & Munari (2002) - +37.1 and +64.9 km s-1. Kolev et al. (2002) found velocities +37.3 and +63.9 km s-1. They suggested that the bluer component could be of circumstellar origin together with unusually strong DIBs (see also Sect. 3.2.3).
Using the calibration of Munari & Zwitter (1997) we found the reddening from the equivalent widths of IS lines (Table 1). The total reddening is E(B-V)=0.49 and the total visual extinction AV=1.5. This extinction corresponds to a distance of 3.1 kpc if the extinction maps by Neckel & Klare (1980) are used. The other estimates of E(B-V) span from 0.7 (Kimeswenger et al. 2002) to 0.9 (Munari et al. 2002b). The latter value was derived assuming that the apparent hot companion star (or just the star in the same position) is a genuine B3 dwarf. The distance to that star is then 10.5 kpc. V838 Mon is placed in a relatively low reddening region and the relevant field of the Neckel & Klare (1980) calibration is quite large and apparently overestimates the reddening towards V838 Mon. The dust distribution maps by Schlegel et al. (1998) allow us to determine the upper limit of E(B-V)=0.89 in that direction. This coincides with the Munari et al. (2002b) value, therefore fixing the upper limit of the distance at 10.5 kpc.
In the Leiden/Dwingeloo H I survey (Hartmann & Burton 1997) in the direction of V838 Mon three complexes of clouds with the mean heliocentric velocities of 67.7, 42.8, and 27.1 km s-1 (or 50.3, 25.3, and 9.7 km s-1 relative to LSR) could be identified. These velocities are close to the ones found from the Na I D IS lines.
Assuming that the observed radial velocities reflect the Galaxy rotation and
the peculiar velocities of the clouds do not exceed 7 km s-1, the
distances of these clouds would be kpc,
kpc, and
kpc. These estimates, however, being very crude, still allow us
to place V838 Mon not much closer than 4 kpc. Wisniewski et al. (2003) used
Na I D IS line radial velocities together with the velocity map of
the outer Galaxy by Brand & Blitz (1993) and found
kpc as
the lower limit of the distance to V838 Mon.
These kinematic distances are consistent with the distances estimated
from the reddening.
On the other hand, these distance estimates are in sharp disagreement with the
ones found using the light echo expansion measurements.
Munari et al. (2002a)
estimated in this way the distance of
pc and
Kimeswenger et al.
(2002) found the distance of 640 to 680 pc
relying on the assumption of spherical geometry.
However, Bond et al. (2002a)
used imaging of the V838 Mon light echo in polarized light with
the Hubble Space
Telescope on May 20, 2002 and
found the lower limit of the distance 2.5 kpc if
the polarized ring was illuminated by
the light emitted in the March 11 outburst.
If that light was emitted on February 7, when the star was the brightest, the
distance would be 3.6 kpc. The latest estimates from the light echo
expansion by Bond et al. (2002b) give the distance between 3 and 7 kpc,
again in the same range as the kinematic distances and the distances found
from the interstellar reddening.
Using the polarimetric imaging with HST Bond et al. (2003) even set
the lower limit to the distance at 6 kpc.
Table 1: Heliocentric velocities, equivalent widths, and corresponding reddenings E(B-V) of interstellar Na I D lines in V838 Mon spectra.
The spectra show also the presence of diffuse interstellar
bands (DIB) at 5780.37, 5796.96, 5849.80 and 8620.79 Å (mean rest wavelengths from
Galazutdinov et al. 2000). The observed DIB wavelengths, equivalent
widths, and corresponding reddenings are given in Table 2.
The reddenings
are estimated from equivalent widths using the calibration data
by Jenniskens & Désert (1994) and Josafatsson & Snow (1987).
These reddenings turn out to be much larger than those obtained from
Na I IS lines. However, the reddenings derived from
the visual region DIBs should
be considered very inaccurate as their extremely different strength ratios
indicate (Kreowski & Schmidt 1997).
The velocities, however, are close to the velocity of the most redshifted Na I D IS line velocity and also to the assumed
systemic velocity. Therefore these DIBs do not originate from the same
region as the blue IS components of the Na I D doublet as proposed by
Kolev et al. (2002).
On our low resolution spectra a medium intensity DIB at
is
visible. Munari (1999) has
shown that the equivalent width of this DIB correlates tightly with the reddening. Our estimate of E(B-V) from
this DIB is very close to the one from the Na I D IS lines.
Table 2:
Observed wavelengths,
heliocentric velocities, and equivalent widths of DIBs in the spectrum
of V838 Mon, and
corresponding reddenings E(B-V). For the visual region DIBs
the first column of E(B-V) corresponds
to Jenniskens & Désert (1994) calibration and the
second to Josafatsson & Snow (1987) data. For
DIB
the calibration by Munari (1999) was used.
Our sequence of low resolution spectra of V838 Mon covers the period from 4th of February to 4th of April 2002 representing the brightness evolution from the rising branch of of the second outburst to the beginning of the deep light weakening after the third outburst (Fig. 1).
The spectral evolution of V838 Mon after the second outburst is best
characterized by the continuous decrease of excitation with a
temporary increase
during the third light maximum and possible expulsion of a new amount
of material during the minimum between the second and third maxima.
We have found that the intensity pattern of typical lines belonging to both
neutral and ionized species has been remarkably correlated with the
lightcurve in parallel to the main trend of dropping excitation.
This trend is illustrated in Fig. 5 where the spectra
in the H region are depicted. The vertical shift of individual spectra
in this figure is
proportional to the time which runs downwards. One could clearly see the
gradual weakening of Fe II and increase of Fe I lines. Also,
the emission in Fe II and H
decreases rather fast.
![]() |
Figure 5:
The temporal change of the V838 Mon spectrum in the vicinity
of the H![]() |
Open with DEXTER |
Just before the second outburst, the
H line showed strong and very broad
wings extending more than
1500 km s-1 (Fig. 6). These wings disappeared
10 days later (JD 2 452 320). In the spectra obtained
by Buil (2002) these
wings are slightly visible until JD 2 452 318
(below we use for marking the dates shorter expressions used as
tick labels in figures e.g. JD..318 instead of JD 2 452 318).
According to Buil's observations
the intensity of the broad
wings changes together with the central intensity of the
H
line (Fig. 7).
Osiwa
a et al. (2003) have suggested Raman
scattering of Ly
photons by neutral hydrogen as an explanation
of these wings.
![]() |
Figure 6:
Some representative profiles of H![]() |
Open with DEXTER |
![]() |
Figure 7:
The temporal change of H![]() |
Open with DEXTER |
At the time when H
broad emission wings were visible, the absorption
component had a radial velocity of about -180 km s-1 and total
width about 150 km s-1. When the broad wings disappeared, the absorption
velocity increased to -200 km s-1 and the width to 300 km s-1.
Starting from JD..326 a new absorption component starts to form and is
clearly visible 5 days later at velocity -60 km s-1. This is the
time when the visual brightness of the star had its minimum between the second
and third outburst. At the same time the H
emission was slightly
stronger than on the neighbouring days.
The radial velocity of the
H
main absorption trough at -200 km s-1did not change during the whole
following period covered by our observations,
although the intensity of the trough
was decreasing.
In contrast to that, the originally developed
P Cygni type profiles of ionized species (Si II, Fe II),
while weakening, also decreases the blueshift of the absorption trough.
On JD..320 a local lightcurve bump starts to develop. This bump is not visible in B, but is slightly noticeable in V, and well developed in I (Bond et al. 2003).
Simultaneously with the H
absorption component at
a radial velocity -60 km s-1, Fe I lines start to strengthen (Fig. 8). Until that time
Fe I absorption lines were extremely weak and the radial velocity
estimates for the
line in Fig. 8
are very uncertain but not
in conflict with estimates from high resolution spectra for JD..310
(see Sect. 3.2.1).
The new absorption component developed at larger blueshift than the
disappearing one. During the periods of deepest absorption a few percent
emission peak developed in Fe I
.
Other Fe I lines
behave similarly.
The strengthening of Fe I lines is reversed during
the third light maximum and continues when the star weakens again. During
all that time the blueshift of Fe I lines decreases at a mean rate 1 km s-1 d-1. One of the best measured cases of the Fe I
line is depicted in Fig. 9.
The intensity of the emission and absorption components of the
low excitation Ba II
line (
eV) increases
considerably between JD..313 and JD..316. The absorption component
continues to grow for
several days after that, and then stays at nearly the same level until the end of
our observations.
![]() |
Figure 8:
The run of the residual intensities and corresponding
radialvelocities of absorption component of the Fe I
![]() |
Open with DEXTER |
![]() |
Figure 9:
The change with time of the radial velocity
of the absorption component of Fe I
![]() |
Open with DEXTER |
This means that a new layer of mass was expelled around JD 2 452 320.
Further evidence supporting that suggestion is given in Fig. 10,
where the
temporal run of the residual intensity and the
radial velocity of two absorption
components of Si II
line is presented. One of these
components, which had its
maximum strength slightly after the second light maximum,
disappeared and a new one started to form. The velocity change of this new
component is close to the Fe I line velocities. The line strength,
however, changes in the
opposite direction. This is what is expected if the line
strength is mainly determined by excitation conditions.
The low excitation strong lines such as Na I D lines and Ca II
behave somewhat differently - even at the end of March
high terminal velocities up to -350 km s-1 were observed. Before
the second light outburst these lines showed terminal velocities up
to -500 km s-1 (Munari et al. 2002a; Kolev et al.
2002) but
highly blueshifted wings
of Na I D lines disappeared during several days before the
second outburst. A possible explanation for that is that the sudden rise
of colour temperature at the first moments of the second brightening
caused much higher excitation in the whole volume of gas outflow.
During the minimum between the second and third outbursts, Na I D lines develop stronger blue absorption line wings
again. At the third
outburst these wings disappeared again (Fig. 11).
![]() |
Figure 10:
The run of the residual intensities and corresponding
radial velocities of two absorption components of Si II
![]() |
Open with DEXTER |
![]() |
Figure 11: Profiles of Na I D2 line for various dates. Note, that at the light maxima on February, 04 and March, 08 the profiles are virtually identical. The spectral resolution of the Feb., 04 spectrum is much higher than for the other spectra. |
Open with DEXTER |
Munari et al. (2002a)
have interpreted the broad-band colour variability of V838 Mon in terms of changing blackbody temperature having maxima at
the peaks of lightcurve. This is in accord with our findings on spectral lines.
We have made a crude quantitative estimate of the excitation temperature in the
absorption line forming region by comparing our spectra of V838 Mon in the
blue
and near-IR
regions
with a set of standard spectra of supergiants ranging from A0 I to K0 I. These
stars were observed at the Tartu Observatory in the blue region with the same
resolution as V838 Mon. The near-IR reference spectra of comparable resolution
were obtained from the database described by Cenarro et al. (2001).
Our spectral classification
results for different dates are presented in Table 3. In the fourth column of
this Table the corresponding effective temperatures of supergiants
(Drilling & Landolt 2000) are presented.
In the fifth column the blackbody
temperatures by Munari et al. (2002a) are given.
The temperatures from the
spectral classification are higher by
up to 2000 K. It is possible to remove this
difference by assuming larger IS reddening when reducing broad band data.
Munari et al. (2002a)
preferred
E(B-V)=0.5 but larger values up to 0.9
were later suggested
by themselves (Munari et al. 2002b). The two last colums of
Table 3
show the observed colours (B-V) (Goranskii et al. 2002) and
unreddened colours (B-V)0 for the corresponding spectral class. The
resulting colour excess is close to 1.0.
Table 3: Spectral classification of V838 Mon.
On February, 14
(JD..320) we noticed that a few percent emission appeared
in He I
line. A closer inspection revealed that the
weak emission line was present between 14th and 25th of February 2002.
The emission peak of the line was blueshifted by 150 km s-1
with respect
to other emission lines which had nearly constant heliocentric radial velocity
around +100 km s-1. This line of highest excitation
(
eV)
was probably not directly connected to the temperature
rise at the beginning of the second outburst, otherwise it should have been
visible already on JD..310. Combining our data with the similar ones
by Wisniewski et al. (2003) we confirm that
the He I
line was absent up to JD..315. We tried to check the presence
of He I
line. This line is severely blended
with an Fe I line. Up to JD..313 the blend behaved as other neutral
iron lines. Starting from JD..316 a clear blend of two lines appeared and
the presence of extra emission was obvious. Later the intensity of the
neutral iron lines started to grow and the contribution from He I could not been monitored. In summary, we are convinced that He I emission was present from Feb., 10
(JD..316) until the end of the month. The origin
of these lines should be different from the other lines considered. One
could suspect that they formed near an another hotter star. The B3V companion found by Desidera & Munari (2002) could be such a star if it
really is the companion to V838 Mon.
In the spectra obtained in March the P Cygni emission components were significantly weakened and therefore there was some hope that the abundance analysis using spectrum synthesizing could be performed. The red echelle spectra obtained in Ondrejov on March 18, 2002 are of reasonable quality for that.
According to Munari et al. (2002a)
the energy distribution in the continuum
corresponded to 4600 K on March 08 and to somewhat
lower temperature on March 18.
If the mass of the star is around 1 and the distance is 6 kpc, the absolute
bolometric magnitude with this temperature (4500 K) will be
and the surface gravity
.
With these parameters we started our
experiments.
A small grid of atmospheric models with solar abundances
was computed for temperatures
K and gravities
.
The model atmospheres were computed with an improved version
of the MARCS program (Gustafsson et al. 1975)
which was modified in order to
include molecular and metallic line opacities in an OS approximation.
Computed spectra with the starting temperature show a strong molecular (CN)
spectrum which is not observed. This forced us to raise
to 4750 K, assuming solar abundances of C and N and
.
This temperature is close to the mean of the colour temperature
(Munari et al. 2002a)
and the excitation temperature derived from spectral
classification.
To compute the CN spectrum we used the data from the SCAN-CN
tape (Jørgensen & Larsson 1990).
For metallic lines we used the line-list
by Bell (1976) with some oscillator strengths corrected according
to Thevenin (1989, 1990).
The simultaneous fit of the weak and strong Fe I lines
allows us to fix
km s-1. The large uncertainty
here is due to the small number of strong lines. The abundances from lines
stronger than
suffer from this uncertainty.
The observed lines are much
wider than the instrumental profile. We attribute this extra broadening
to the macroturbulence of 75 km s-1.
The great width of the lines makes the blending large. We could not
therefore use the many lines that are clearly present but hopelessly blended.
The use of the ionization equilibrium of iron allowed to confirm the
surface gravity
.
The resulting abundances are given in Table 4. We stress that the indicated uncertainties show only the differences in abundances found from different lines. The atmospheric parameters were not firmly established and the static LTE atmospheric models are very sketchy for this star, which on March, 18 still showed P Cygni-type lines - clear indicators of mass outflow. Therefore these abundances should be considered only as very crude estimates.
From Table 4 one finds that the abundances of most elements
are slightly less than the solar ones,
.
At the same time the lithium
,
barium
,
and lanthanum
abundances are clearly
enhanced.
The abundances of Li and Ba are found
from strong lines assuming
km s-1. With
higher microturbulent velocity these overabundances are less.
These abundances do not resemble the ones found in the "late He-flash''
objects V4334 Sgr and FG Sge. This is an argument for not including
V838 Mon in that type of object.
Table 4: Abundances in V838 Mon on March, 18 2002.
Several scenarios have been proposed to explain the outburst of V838 Mon. The spectral evolution and mainly the absence of any nebular phase excludes the simplest scenario of classical nova. The nondetection of X-ray luminosity a year after the outburst with Chandra (Orio et al. 2003) is also strong evidence against nova-like thermonuclear runaway as a mechanism for the outburst.
The scenario of V838 Mon being a post-AGB object having a late helium flash was also discussed (Munari et al. 2002a) and rejected, mainly due to too fast evolution and the lack of strong carbon overabundance. The abundances found by us also rule out this model.
Two stars, V4332 Sgr in the Galaxy and a red variable in Andromeda
galaxy (RV M31),
were found to behave similarly to V838 Mon (Munari et al. 2002a;
Kimeswenger et al. 2002). Kimeswenger et al. rule out RV M31 as a
possible twin to V838 Mon on the basis that V838 Mon had the maximum
luminosity much lower than RV M31. In the present paper we
offer some arguments favouring
the large distance for V838 Mon, and therefore it could have been as luminous
as RV M31 which peaked at
(Rich et al. 1989).
If the distance of V838 Mon is greater than 6 kpc its absolute
bolometric magnitude peaked at
.
In both cases (RV M31 and V4332 Sgr) there is no good explanation. For RV M31 Iben & Tutukov (1992) proposed a model where such a cold outburst of a short period CV is possible. But this would then led to the complete destruction of the WD. In case of RV M31 at least two outbursts have been recorded (Sharov 1990). Sharov, however, later retracted his claim due to the low quality of the photographic material used (Sharov 1993). In the case of V838 Mon there have been previous outbursts which left the circumstellar envelope visible in the light echo. The nondetection of X-rays is also an argument against this scenario in the case of V838 Mon.
Soker & Tylenda (2003) proposed a scenario of merging of two main
sequence stars of masses 1.5 and about 0.5 .
This kind
of merging also destroys the star and therefore is not applicable for
recurrent outbursts. But they also suggest the merger of a zero-age horizontal
branch star with a low mass main sequence star. Then the dust-halo
illuminated by light-echo is a result of earlier high mass loss on RGB.
However,
in his recent study of the light echo of V838 Mon, Tylenda (2004)
found that the light echo is of interstellar origin rather than produced by
dust generated by
mass loss in the past.
If nature produces one of the merging scenarios, V838 Mon should
be a triple system.
The progenitor was estimated to be a hot
star (Munari et al.
2002a; Kimeswenger et al. 2002; Goranskii et al.
2002). After the outburst
faded, this star is still there (Munari et al. 2002b).
Although there is the
possibility that this hot star lies in the same direction but
farther away than V838 Mon, the brief existence of emission in the
He I
line seems to support the idea of
physically bound multiple star.
Acknowledgements
The financial support to this work was provided by Estonian Science Foundation grant 5003. This study was also supported by grant 205/02/0445 of Grant Agency of the Czech Republic and by project 02-02-16085a of the Russian Foundation for Basic Research. We thank the referee, J. Landstreet, for the detailed and constructive remarks.