A&A 465, 993-1002 (2007)
DOI: 10.1051/0004-6361:20066401
J. H. Groh1 - A. Damineli1 - F. Jablonski2
1 - Instituto de Astronomia, Geofísica e Ciências Atmosféricas,
Universidade de São Paulo, Rua do Matão 1226, Cidade Universitária,
05508-900 São Paulo, SP, Brasil
2 - Instituto Nacional de Pesquisas Espaciais/MCT, Avenida dos
Astronautas 1758, 12227-010 São José dos Campos, SP, Brasil
Received 14 September 2006 / Accepted 8 November 2006
Abstract
We present a digital atlas of peculiar, high-luminosity massive stars in the
near-infrared region (10 470-11 000 Å) at medium resolution (
). The spectra
are centered around He I 10 830 Å, which is formed in the wind of those stars, and is
a crucial line to obtain their physical parameters. The instrumental configuration also
sampled a rich variety of emission lines of Fe II, Mg II, C I,
N I, and Pa
.
Secure identifications for most spectral lines are given,
based on synthetic atmosphere models calculated by our group. We also propose that two
unidentified absorption features have interstellar and/or circumstellar origin.
For the strongest one (10 780 Å) an empirical calibration between E(B-V) and equivalent
width is provided. The atlas displays the spectra of massive stars organized in four
categories, namely Be stars, OBA Iape (or luminous blue variables, LBV candidates and
ex/dormant LBVs), OB supergiants and Wolf-Rayet stars. For comparison, the photospheric
spectra of non emission-line stars are presented. Selected LBVs were observed in different
epochs from 2001 to 2004, and their spectral variability reveals that some stars, such as
Car, AG Car and HR Car, suffered dramatic spectroscopic changes during this time
interval.
Key words: atlases - stars: emission-line, Be - stars: Wolf-Rayet - stars: winds, outflows - stars: early-type
In the last decades tremendous progress has been achieved in understanding how a massive star evolves along the HR diagram. A new scenario of massive star evolution was developed after the models allowed for the effects of rotation in the stellar structure and evolution (Maeder & Meynet 2000b). Also, the physical parameters of massive stars can currently be constrained with the advent of a new generation of fully-blanketed, spherical-symmetric, non-LTE radiative transfer codes (Gräfener et al. 2002; Puls et al. 2005; Hubeny & Lanz 1995; Hamann & Koesterke 1998; Hillier & Miller 1998; Pauldrach et al. 2001). Even for extremely massive objects such as Eta Carinae (Hillier et al. 2001) we now have a much better comprehension of the physical parameters of the central star by analyzing the emerged spectrum.
However, the short transitional stages of massive star evolution are still poorly understood. The radiation emitted by the photosphere of these stars interacts with the optically-thick, dense wind, which makes the spectral classification a very challenging task. In addition to the luminosity, effective temperature and surface gravity, several other stellar parameters might change the spectral morphology of these objects, notably the mass-loss rate and the wind terminal velocity. In the discussion on the characteristics of the atlas, we point out the similarities and differences in the spectral morphology of massive stars belonging to different evolutionary stages. wf2000WF2000
The aim of this paper is to extend the idea of the OB Zoo published by Walborn & Fitzpatrick (2000, hereafter WF2000)
to the near-infrared region of the spectrum.
Some of the objects
present in this work were observed at lower resolution from 1.0 to 2.5 m by
McGregor et al. (1988). In addition, peculiar massive stars were also observed
around 8500 Å and 10 000 Å by Lopes et al. (1992), and in the H and K bands
by Morris et al. (1996).
We chose to explore the "forgotten'' region around
He I 10 830 Å due to the strength of this line,
which is often the strongest one in
the near-infrared region. Pioneering
observations of early-type stars around this line were published by
Andrillat & Vreux (1979) and by Meisel et al. (1982).
Especially for objects with dense winds and high terminal
velocities, He I 10 830 Å is one of the few strong He I lines that is unblended. This line is
formed over a large region of the wind and thus it is a crucial diagnostic of the physical
properties of the wind and of the underlying photosphere. Moreover, with a relatively short
spectral coverage around 10 830 Å, we can sample a number of strong, unblended lines,
such as Pa ,
and rich emission-line spectra of Fe II, Mg II, C I,
and N I. They provide important constraints on the physical conditions found in
those stars and their winds.
We observed other peculiar objects that are not present in wf2000, but are as fascinating and intriguing as their original program stars. Meanwhile, we did not observe some objects present in wf2000, especially the faintest ones, that were beyond the limiting magnitude of our instrumental configuration. Nevertheless, we tried to cover a large range of evolutionary phases and physical parameters found in hot stars. However, it is beyond the scope of this paper to present a full evolutionary sequence of massive stars around 10 830 Å, especially considering that a significant fraction of the objects shown here does not have a well-determined evolutionary stage.
As part of the atlas, we present the spectral variability of a large sample of luminous blue variables (henceforth LBV) and LBV candidates during a time interval of 3 yr. To the best of our knowledge, this is the first work in any spectral region to present the spectroscopic variability on a timescale of a few years (which is typical for S Dor-type variations), for a large sample of LBVs and related objects. As we show along the atlas, different behaviors of variability can be noticed (radial velocity, line profile, strength).
This paper is organized as follows. In Sect. 2 we describe the observations made at the Observatório Pico dos Dias (LNA/Brazil); in Sect. 3 we present the atlas, briefly contextualizing and introducing each class of objects, and then describing the past bibliography and the spectral features of each individual object. In Sect. 4 we discuss the unidentified features present in the program stars.
Table 1: Basic data for the program stars, ordered by increasing right ascension. Projected rotational velocities are indicated for Be stars, while the presence of a circumstellar nebula is indicated for LBVs and related objects. References are: CN04 = Clark & Negueruela (2004); Cro06 = Crowther et al. (2006); Mir01 = Miroshnichenko et al. (2001); L97 = Lyubimkov et al. (1997); LL06 = Levenhagen & Leister (2006); Riv97 = Rivinius et al. (1997); R97 = Roche et al. (1997); S06 = Smith & Balona (2006); Sta03 = Stahl et al. (2003); vG92 = van Genderen et al. (1992); vG01 = van Genderen (2001); WF2000 = Walborn & Fitzpatrick (2000); Z05 = Zorec et al. (2005).
The near-infrared spectra were obtained using a HgCdTe
infrared array
(CamIV) mounted on the Coudé focus of the 1.6 m telescope at the Observatório do Pico
dos Dias/LNA (Brazil). The stars were observed on 2001 June 09-12, and again
on 2004 April 30 and 2004 July 02-04. Spectra of variable stars were also taken on 2002 November 04, and 2003 June 28.
For each object, we combined up to 35 individual exposures varying between 10 s and 60 s
to reach typically
in the continuum. We used a 600 gr/mm grating and a 1-arcsec
slit width to obtain
and 0.64 Å/pixel of dispersion, covering the
spectral range of 10 470-11 000 Å. The sky+thermal background subtraction was
accomplished by dithering the star along the slit in 7 positions. We median-combined
these 7 frames to obtain a median sky image, which was then subtracted from each original frame. Each
sky-subtracted image was then divided by a normalized flat-field image. The following data reduction,
extraction of the spectrum and analysis was made by using usual IRAF routines. The wavelengths were
calibrated using an internal ThAr emission lamp, achieving a typical rms uncertainty of 0.08 Å
(
2.2 km s-1). The telluric features were removed from the program stars' spectra by
dividing them by the spectrum of a fast-rotating hot star. The photospheric Pa
and
He I 10 912 Å absorptions were previously removed from the fast-rotating star by
interpolating the adjacent continuum with a low-order Legendre polynomial. The extracted spectra were
corrected to the heliocentric frame and normalized to the continuum intensity using a low-order
Legendre polynomial (errors in normalization are typically 1%). No flux calibration or de-reddening was done.
The atlas consists of 10 figures, containing from two to six spectra of hot luminous stars each, thus 40 peculiar objects. The spectral evolution of selected objects during 2001-2004 is shown in Figs. 16 to 18. An additional figure (Fig. 15) contains the spectra of normal early-type stars, to be used as a reference for the characteristic photospheric features. Each figure shows the air wavelength in Å versus the normalized intensity, plus an arbitrary offset, except where indicated otherwise. Some figures were re-sampled to show only the region that has strong spectral lines for a given group of stars. The peculiar objects were grouped based on their spectroscopic characteristics, and the basic data of these stars are presented in Table 1.
Many stars have been extensively studied in other spectral regions, whilst some of them
have very scarce recent bibliography. Compared to the previous published data in this
spectral region, the higher resolution and higher S/N of our spectra allowed us to identify many
faint spectral lines. They were identified by comparing the value of the oscillator strength of the
transition with the observed line intensities, and then checking for the presence of lines
of the same multiplet or ionization stage in other wavelength regions (e.g. optical).
Almost all the identifications were confirmed by comparing the observations with synthetic
spectra (Groh et al. 2007, in preparation), using the radiative transfer code CMFGEN
(Hillier & Miller 1998). However, there are still
unidentified spectral features remaining, which will be discussed in Sect. 4. In Table 2 we summarize the detected features in the spectra of
the program stars; the air wavelengths were taken from the Atomic Line List
v2.04. The line identifications are
overploted on the spectrum of the LBV AG Carinae taken in 2004 April (Fig. 1).
Table 2: Observed line transitions in the program stars.
The atlas is organized as follows: in Sect. 3.1 we present the spectra of the Be stars, in Sect. 3.2 we show the spectra of the LBVs and LBV candidates, in Sect. 3.3 we group the OB supergiants, and in Sect. 3.4 the spectra of the Wolf-Rayet stars. In each subsection we briefly introduce the general characteristics of the class, and discuss the relevant aspects and details of the observed features for each star, comparing them with another objects of our atlas when appropriated.
The 13 classical Be stars of our sample are shown in decreasing order of spectral type
in Figs. 5-7. Most stars
in this class show Fe II lines, He I 10 830 Å and Pa
in emission, often
double-peaked. The sole exception is HD 120991, which does not show Fe II lines. Some of the Be
stars also show emissions in C I and Mg II lines, that have the same morphology as the
Fe II emission.
(Obj)X Persei(Obj) (Fig. 7) is a Be/X-ray binary in our sample, such as HD 110432
(Clark et al. 2001). X Persei is the optical
counterpart of the pulsar 4U 0352+30 (Clark et al. 2001). The system has a low eccentricity
(e=0.11) and a long orbital period of 250 days (Delgado-Martí et al. 2001). Its
optical spectrum is characterized by both H
and He I 6678 Å showing significant
variability (Clark et al. 2001). In the near-infrared, its spectrum shows double-peaked emission in Fe II 10 500 Å (-100 and 195 km s-1), while He I 10 830 Å presents a weak
double-peaked profile (-98 and 13 km s-1), such as Pa
(-100 and -5 km s-1). Fe II 10 862 Å
is present too, but it is very weak. Mg II lines at 10 915 and 10 952 Å are
narrow.
(Obj) Scorpii(Obj) (Fig. 7) has been claimed to be a classic Be star based on
its optical and infrared spectroscopy history, recently undergoing a H
outburst (Fabregat et al. 2002). However, for many years it was
considered a typical B-type star. Previous observations did not show any emission
lines (Grigsby et al. 1992; Heasley & Wolff 1983).
Scorpii is in a highly eccentric binary
system, with an orbital period of 10.6 yr (Banerjee et al. 2001). Since the appearance
of emission lines coincided with the periastron passage,
binarity could be an alternative hypothesis to explain the Be phenomenon in this object (Miroshnichenko et al. 2001; Fabregat et al. 2000).
In our near-infrared spectrum
Scorpii shows double-peaked Fe II 10 500 Å and 10 862 Å
(-65 and +130 km s-1 and -45 and +80 km s-1, respectively). C I
10 683, 10 685 and 10 691 Å (blend), C I 10 707 Å, and C I 10 729 Å are
present in emission. The He I 10 830 Å line is strong, while Pa
is weak. The
Mg II lines are broad and may be double-peaked.
(Obj)HD 110432(Obj) (Fig. 7) is a Be star that is located just beyond the southern part of the
Coalsack at a distance of
pc (Rachford et al. 2001). It is one of
the brightest stars lying behind interstellar material and it has narrow
absorption resonance lines in the UV spectrum (Codina et al. 1984; Prinja & Henrichs 1987).
The optical data displays line-profile variations from night to night (Codina et al. 1984).
Torrejón & Orr (2001) classified HD 110432 as a low luminosity Be/X-ray binary.
Its near-infrared spectrum looks like
that of
Scorpii, although we have found a symmetrically double-peaked
profile in He I 10 830 Å and in Pa
,
both with peaks at -80 and 80 km s-1.
As in
Centauri, a weak feature was found at 10 541 Å and we identified
it as N I 10 539.6 Å.
(Obj) Ophiuchi(Obj) (Fig. 6) displays fast variability in the Balmer series
lines (Doazan 1970). Borisova (1992) observed variations in the He I optical lines
with a period of 0.913 day, and suggested that
Oph is a non-radial pulsator.
Its near-infrared spectrum is dominated by narrow emissions lines of Fe II, Mg II,
C I and Pa
.
The C I lines are more prominent than in the other Be stars of our
sample. The He I 10 830 Å emission is weak.
![]() |
Figure 1:
Identification of the features present in the 2004 April spectrum of the LBV
AG Carinae (full line). Thanks to the high resolution and high S/N provided by
the infrared spectrograph (CamIV) we could identify a number of weak features present in the spectrum. The dashed line
shows a non-LTE
radiative transfer model for AG Car (Groh et al. 2006, in preparation) which was used to
confirm the line identifications. Extended, strong electron-scattering wings are noticeable in
He I 10 830 Å and Pa ![]() |
Open with DEXTER |
(Obj) Cygni(Obj) (Fig. 6) is a classical Be star first detected by
Fleming (1890). This object shows strong H
emission (Neiner et al. 2005), variable on
different timescales. Outbursts are quite common in this object, such as the one detected by the
Hipparcos satellite. In our atlas its spectrum resembles that of
Centauri, except that the
C I lines are stronger, Mg II lines are weaker, and Pa
is symmetrically
double-peaked (-50 and +55 km s-1). The Fe II 10 500 Å and 10 862 Å lines are
double-peaked as well (-30 and +115 km s-1 and -36 and +77 km s-1, respectively).
(Obj)66 Ophiuchus(Obj) (Fig. 5) is a well-observed Be star. Peters (1982) analyzed the first IUE
observations and detected high velocity absorption components (from -250 to -850 km s-1) in the
resonance lines of C IV, Si IV, and Si III. Based on the UV spectrum, strong variations were detected in
its stellar-wind lines (Barker & Marlborough 1985), vanishing and appearing in intervals of several months
(Grady et al. 1987). 66 Oph has a long history of optical variability, showing H emission with
irregular variations. Cyclic V/R variability was also detected in H
during the period 1982-1993.
In our near infrared spectrum 66 Oph shows weak Fe II, C I and He I emission
lines. Pa
is almost symmetrically double-peaked (violet and red peaks at -110 and +90 km s-1, respectively), reaching up to 0.4 of the continuum intensity (i.e.,
/
).
Pa
also presents an absorption profile on the blue side that can be P Cygni absorption.
(Obj)
Centauri(Obj) (Fig. 5) is a bright and fast-rotating star
(
= 310 km s-1 ), and is a member of the Scorpio-Centaurus OB association
(de Geus et al. 1989). It rotates at 0.65 of its critical velocity (Hutchings et al. 1979) and the estimated
mass-loss rate from Si IV doublet profiles is log
yr-1 (Snow 1982).
Stefl et al. (1995) found a nearly sinusoidal variation in the radial velocity of
Si III 4552 Å with a period of
0.6424 day. In the near-infrared spectrum
Cen presents weak lines of C I and Mg II.
Fe II 10 500 Å is symmetrically double-peaked (-175 and 260 km s-1), while Fe II
10 862 Å and Fe II 10 872 Å are blended. He I 10 830 Å displays a shell feature
(central intensity of
), indicating that this object is seen nearly at the equator. Pa
shows a
complex structure, with three emission peaks.
(Obj) Centauri(Obj) (Fig. 6) is a member of the Sco-Cen OB association,
and Fleming (1890) again was the first to detect its emission-line
features. Hanuschik et al. (1995) found V/R variability in H
and Fe II lines, which he
proposed to be due to global disk oscillations. Clark et al. (1998) identified
Cen
as a radio source, which may be variable. In the near-infrared region it shows a rich
variety of Fe II lines. Fe II 10 490 Å and Fe II 10 525 Å are weak (
= 1.05),
while Fe II 10 500 Å appears double peaked (-60 and +80 km s-1). The Fe II 10 862 Å profile is double-peaked too (-88 and +30 km s-1). We found weak emission lines at 10 541 Å and 10 550 Å which we
identified as N I 10 539.6 Å and 10 549.6 Å, respectively. The blend of C I lines
is weak (
= 1.05), peaking at 10 684 Å and 10 692 Å. The He I 10 830 Å line shows a broad profile, reaching up to
= 1.15.
Pa
is strong (
= 1.25) and the peak is
slightly redshifted (+25 km s-1). The Mg II lines are present as well.
(Obj) Arae(Obj) (Fig. 5) does not show V/R variability in the optical spectrum, presenting
only a constant and moderated Balmer emission (Hanuschik et al. 1995; Mennickent 1991). In our spectrum it presents
double-peaked Fe II emission both in 10 500 Å (-105 and +142 km s-1) and 10 862 Å (-60 and +80 km s-1). As almost all Be stars of our atlas, it displays a broad emission at 10 686 Å that
we identified as a blend of C I lines (10 683 Å, 10 685 Å and 10 691 Å). An unidentified
absorption feature was found at 10 906 Å.
Arae also shows weak emission of Mg II 10 915
Å and 10 952 Å, and very weak (
= 1.05) emission of He I 10 830 Å. Pa
is strong and double-peaked (-109 and +85 km s-1).
(Obj) Pavonis(Obj) (Fig. 7) is a intermediate-rotating Be star
(
km s-1 ). This object revealed highly-variable H
emission episodes
during 1984-1987 (Mennickent 1991). Our monitoring at Observatório Pico dos Dias
have shown many episodes of fading and reappearance of H
during the last 15 yr (Damineli 2007, in preparation).
Its optical spectrum shows strong, double-peaked H emission (Sahade et al. 1988).
Chen et al. (1989) analyzed its IUE spectra and found that the profiles of the UV lines are
rotationally broadened into two groups, with rotational velocities of 170 and 210 km s-1. They
inferred two different regions in the extended stellar atmosphere, one that is rotating and
one expanding, carrying the wind out at a maximum velocity of 155 km s-1.
Only absorption lines are found in the IUE spectra of
Pav, ranging from neutral
species (e.g., C I and N I) to highly ionized ions (C IV and Si IV).
This object was detected as a
hard X-ray source by the Einstein Observatory and might have a neutron star as a companion.
The only line present in its near-infrared spectrum is He I 10 830 Å, with
= 1.45, and symmetrically double-peaked (-130 and +130 km s-1).
(Obj)HD 142983(Obj) (Fig. 5) is a V/R variable star with an extended atmosphere. It shows cyclical
variations with a period of about 10 yr. In the UV spectrum it does not show any emission lines, and
the optical H and He lines are broadened by rotation (McLaughlin 1961; Delplace & Chambon 1976). Floquet et al. (1996) showed
multiple-frequency patterns in the line-profile variations in HD 142983. They suggested that non-radial
pulsations and enhanced mass-loss rates are linked in Be stars. The near-infrared spectrum of HD 142983
shows Fe II 10 500 Å symmetrically double-peaked (-135 and +127 km s-1) and Fe II 10 862 Å
double-peaked too, but non-symmetrically (-150 and +95 km s-1). C I and Mg II lines are present,
but very weakly. He I 10 830 Å shows a shell-like absorption profile (
), indicating that
HD 142983 is seen nearly equator-on. The same unidentified feature found in
Arae at 10 906 Å was found here as well. Pa
shows a shell profile, consisting of two emission shoulders separated
by a strong absorption feature, centered at -40 km s-1.
(Obj) Monocerotis A(Obj) (Fig. 7) is the brightest member of a complex
triple system where the three components are Be stars of similar spectral
type (Maranon di Leo et al. 1994).
Mon A shows V/R variability with a
period of 12.5 yr (McLaughlin 1961). Mara et al. (1991) reported
profile variation of Mg II 4481 Å, showing episodic emissions. Its
near-infrared spectrum resembles HD 142983, except for some unidentified
features in the former. It shows Fe II 10 500 Å double-peaked (-48 and +190 km s-1) as well as
Fe II 10 862 Å (-80 and +135 km s-1). Fe II 10 525 Å is also clearly present.
We detected a broad emission around 10 689 Å which we identified as a
blend of C I (10 683, 10 685, 10 691 Å). We also detected narrow emissions at
10 714, 10 723 Å, which we related to a double-peaked emission of N I 10 718 Å (-70
and +150 km s-1, respectively). We found a broad emission at 10 740 Å and very
narrow absorptions at 10 757 Å and 10 770 Å that remain unidentified. He I 10 830 Å displays a shell absorption (
), indicating that this object is seen
nearly equator-on. Pa
is double-peaked (-70 and +110 km s-1), presenting a strong
absorption centered at 30 km s-1. The line intensity drops from
in the Doppler-shifted peaks to
in the dip.
The spectra of OBA Iape stars are presented in Figs. 8-11. We grouped peculiar objects that share common morphological properties in the same figure, although this is almost impossible in some cases. The stars in this group are all related to the LBV phase; some of them are known as bona fide LBVs (such as AG Car and HR Car), some are LBV candidates (W 243 in Westerlund 1, MWC 314) and some are believed to be in a dormant/post-LBV phase (such as He 3-519). For a recent review on the properties of these objects see Humphreys & Davidson (1994), van Genderen (2001), and Clark et al. (2005).
(Obj)HD 168625(Obj) (Fig. 8) is classified as a marginally-dormant LBV based on its last 25-year
lightcurve (Sterken et al. 1999), although a typical LBV nebula was discovered around it
(Nota et al. 1996; Robberto & Herbst 1998; Hutsemekers et al. 1994; O'Hara et al. 2003; Pasquali et al. 2002). Its optical spectrum was recently published by
Chentsov et al. (2003), while Hanson et al. (1996) showed the K-band spectrum. Our spectrum
shows weak emission lines of Fe II, Mg II and NI. He I 10 830 Å shows a deep P Cyg absorption with
a very weak emission, while Pa
is seen in absorption. As in HR Car, we also detected the 10 780 Å and the 10 792 Å absorption features. The comparison between the 2001 and 2004 datasets do not show any evidence of variations in the spectral lines, confirming its current
dormant stage.
(Obj)HD 168607(Obj) (Fig. 8) is only 62'' away from HD 168625 in the sky, and they are both believed to be members of the M 17 association, lying at 2.2 kpc wf2000. At this distance they would be separated by only 0.7 pc in projection, and regarding their similar blue spectra, it has been proposed that they could be a pair of massive twin stars (see discussion in wf2000). The 25-year lightcurve of HD 168607 presented by Sterken et al. (1999) shows the same qualitative behavior of other LBVs, such as AG Car, but on a much longer timescale. In our atlas, the near-infrared spectrum of HD 168607 resembles that of HD 168625, although the former shows more developed emission lines, suggesting a higher activity and a higher mass-loss rate for HD 168607. This objected has also shown spectroscopic variations in the last 3 yr (Fig. 18), with an increase in the absorption of He I 10 830 Å and an overall decrease in the intensity of the emission lines.
(Obj)HD 160529(Obj) (Fig. 8) is a bright A-type hypergiant star. Sterken et al. (1991),
through Strömgren uvby and near-infrared photometry, found a change in the
spectral type from A9 to B8 between 1983 and 1991, and a brightness decrease of
0.5 mag in this period. Based on this behavior, they classified HD 160529 as an
LBV star, and analyzed its lightcurve variations as well. Its spectral
variability was analyzed by Wolf et al. (1974), who found variations in the
Balmer lines profile and in the radial velocities of about 40 km s-1.
More recently Stahl et al. (2003) reported a long-term spectroscopic monitoring
and the recent lightcurve of HD 160529. They found little variation in the
spectral type around A2Ia from 1991 to 2002, mainly using He I 5876 Å as a diagnostic line. In our
atlas, the near-infrared spectrum of HD 160529 is very similar to the spectrum
of HR Car, except that the former shows stronger emission in He I 10 830 Å, with a P Cyg
profile. Moreover, Pa
is seen in absorption. From 2001 to 2004 we detected a
radial velocity variation of 55 km s-1. We also noticed an overall increase in the
strength of the absorption lines, with a new absorption feature appearing at 10 581 Å, and a drastic decrease in the Fe II emission. However, He I 10 830 Å did
not change its strength (Fig. 18).
(Obj)HR Carinae(Obj) (Fig. 8) is one of the rare bona-fide LBV stars in the Galaxy. Its variability
was first reported by Hertzsprung (Hoffleit 1940), and since then several lightcurves have been published
(see van Genderen et al. 1997). The high values of the luminosity and mass-loss rate of HR Car were first noted by
Viotti (1971), and the 1989-1990 spectral variability with the associated changes in the lightcurve
were analyzed by Hutsemekers & van Drom (1991). Their optical spectra show evidence of an expanding atmosphere with
multiple shells. The spectrum of the bipolar nebula associated with HR Car shows evidence of nitrogen
enrichment, as does the star itself. A detailed study of the optical nebula around HR Car is given in
Nota et al. (1997) and Weis et al. (1997), and the mid-infrared imaging of the nebula is presented by
Voors et al. (1997). The most recent analysis of the wind of HR Car was performed by Machado et al. (2002), based on
spectra collected in April 1999. They derived an effective temperature of 10 000 K, a luminosity of
500 000
and a mass-loss rate of
yr-1. Previously,
van Genderen (2001) derived an effective temperature of 8000 K for HR Car in 1992, while Nota et al. (1997) found 20 000 K (May 1995) and 15 000 K (January 1996). The near-infrared spectra of HR Car in June 2001 shows weak
Fe II emission (10 500 Å, 10 525 Å 10 862 Å and 10 871 Å). Pa
shows an
inverse P Cygni profile, while He I 10 830 Å is almost absent, with a very weak broad
emission. HR Car shows a rich absorption spectrum of N I, C I and Mg II. An
unidentified feature was found at 10 780 Å and another at 10 792 Å (see Sect. 4).
The data taken in 2004 show much more variability than in 2001 (Fig. 18). The inverse P Cyg profile of Pa
changed to a pure strong emission, with the equivalent width rising from 2 to 14 Å. The Fe II lines became stronger too, while the Mg II lines that were present in
absorption turned into emission. The line of He I 10 830 Å presented an intrincated profile,
with a double-peak emission superimposed on the previous broad emission.
(Obj)W 243(Obj) (Fig. 8) is a recently-discovered LBV candidate in the massive, young open
cluster Westerlund 1 (Clark & Negueruela 2004). Those authors presented an optical
spectrum characterized by a double-peaked emission of H
and
single-peaked emission in the He I lines. They also reported the strengthening
of H
during 2001-2003. The comparison with older spectra revealed that
the star evolved from B2-5Ia to A2I from 1981 to 2003, suggesting an LBV classification. A K-band spectrum obtained by Groh et al. (2006a) in July 2005 suggests further evolution
towards a hotter temperature.
In this work we obtained near-infrared data in 2004 that confirm its LBV
nature. Its spectrum is similar to other cool LBVs, showing a strong Pa
emission that resembles HR Car as it was in 2004. Another remarkable feature of
cool LBVs, namely, the presence of Mg II, N I and C I lines in absorption, are also seen in W243.
However, He I 10 830 Å emission is present and stronger than in other cool LBVs. The
P-Cygni absorption of this line is noticeable, extending up to 400 km s-1, which is much higher
than the expected wind terminal velocity for A supergiants. Further monitoring of W 243 is required to
address the origin of this high-velocity absorption.
(Obj) Carinae(Obj) (Fig. 9) is one of the most observed stars. It underwent a giant
eruption in the 1840 s that generated the bipolar-shaped Homunculus nebula (Gaviola 1950), and
another outburst in the 1890 s that probably generated another nebula, the "little Homunculus''
(Ishibashi et al. 2003). The reader is referred to Davidson & Humphreys (1997), who reviewed the system and
its environment. Since then, many works related to the nature of the central source in
Car have been
published, mainly regarding the 5.54-year cycle and the binary scenario (Damineli et al. 1997) or alternative
ones (Davidson 1997). Its near-infrared spectrum shows prominent Fe II emission at 10 500 Å, 10 525 Å,
and 10 862 Å, with a broad and a narrow component. Pa
also shows the same shape of the
Fe II lines. He I 10 830 Å is stronger, with a P Cygni profile extending up to -660 km s-1. There are other absorptions at -380, -195 and -45 km s-1. The changes in this
spectral region are dramatic over the 5.54-year cycle (Fig. 17): a few days before 2003.5 (minimum
of the spectroscopic event, see Steiner & Damineli 2004) the narrow emission components from
Fe II, He I and Pa
shrunk until they faded completely. The behavior of He I 10 830 Å is striking: the emission components faded almost completely, while the absorption component of
the P Cyg profile extended up to -1300 km s-1.
(Obj)He 3-519(Obj) (Fig. 9) is a galactic prototype of the WN11h class. Its optical and ultraviolet spectrum is nearly identical to the well-known LBV AG Carinae at minimum (i.e., during the hot phase), as discussed by Smith et al. (1994) and wf2000. Although He 3-519 does not have the typical photometric or spectroscopic behavior of a LBV, it shows a circumstellar nebula that relates it to that class. Therefore, it is believed that He 3-519 belongs to a rare class of post-LBV, pre-WN stars (Davidson et al. 1993). The near-infrared spectrum of He 3-519 is similar to that of HD 316285, but the former displays a stronger He I 10 830 Å emission and is almost free of Fe II lines. The 2001-2004 data evolution shows significant strengthening of He I 10 830 Å, but without changes in the velocity of the edge of the P Cyg absorption (Fig. 17), resembling the pattern of the variations of HD 316285.
(Obj)HD 316285(Obj) (Fig. 9) is a luminous, evolved massive star as discussed in
Hillier et al. (1998). Based on a detailed spectral analysis, those authors found evidence for nitrogen and
helium enhancement on the stellar surface, and a large mass-loss rate of
yr-1 was derived. Although there is
no evidence of photometric variability, Hillier et al. also noticed
spectroscopic variation, specially in the Fe II lines. The near-infrared
spectrum of this star resembles that of
Car, except that the nebular and
high-excitation lines
seen in the former are not present in HD 316285. Also, He I 10 830 Å does not show the additional absorption components seen in the spectrum
of
Car. Regarding the variability of HD 316285, our 3 datasets from 2001 to 2004 show significant strengthening of the emission lines (Fig. 17),
although no variations are noticed in the velocity of the wind (as seen from the P-Cyg
absorption).
(Obj)MWC 314(Obj) (Fig. 10) is a scarcely-studied emission line star. It was first
identified by Merrill (1927), while Allen (1973) reported a late-type spectral energy distribution
for this object. It was recognized as a high-luminosity B[e] star by Miroshnichenko (1996),
who estimated its stellar parameters (log
,
T=30 000 K,
)
and
derived a high reddening (AV=5.7 mag). A comparison between B[e] supergiants and
LBVs led him to propose that MWC 314 is an LBV candidate. A high resolution optical
spectrum of this star is presented by Miroshnichenko et al. (1998); its optical spectrum is
dominated by emissions of H I, He I and Fe II.
Its near-infrared spectrum
resembles AG Car and P Cygni. However, its spectroscopic evolution shows an
intriguing behavior in He I 10 830 Å, with a variable high-velocity absorption
wing extending up to 1500 km s-1 (Fig. 18). This behavior is striking because such
a high velocity is not expected for this kind of object. In Fig. 5 of Miroshnichenko et al. (1998) a high-velocity
absorption can be seen in He I 5876 Å, although errors in the continuum normalization may
affect the interpretation (Miroshnichenko, private communication). Moreover, the same qualitative
behavior appears in
Car during periastron passage, although much more dramatically in the latter
(see Fig. 17).
Additional multi-wavelenght observations
of MWC 314 using a finer temporal sampling are required to address the origin of this high-velocity
absorption component.
The equivalent width of the emission lines also changed by a factor of 1.5 from 2001 to 2004.
(Obj)AG Carinae(Obj) (Fig. 10) is the "Rosetta Stone'' of the LBV class,
lying at 6 kpc from the Sun (Humphreys et al. 1989). AG Car and its environment have
been extensively studied in the last decades. The star is surrounded by a
bipolar nebula first reported by Thackeray (1950), which was ejected from the star 104 yr ago (de Freitas Pacheco et al. 1992; Lamers et al. 2001; Nota et al. 1992). The photometric variability in the
last century was analyzed by van Genderen et al. (1997), while the recent
spectroscopic variability was studied by Leitherer et al. (1994) and Stahl et al. (2001).
The morphology and kinematics of the ejecta were analyzed by Nota et al. (1995),
and the nebular abundances were determined by Smith et al. (1997). Shore et al. (1996)
found a color excess of
E(B-V)=0.65 for AG Car based on IUE data. The optical spectra at minimum were analyzed by Smith et al. (1994),
who found an enhanced helium abundance at the surface of the star. Recently, Groh et al. (2006b)
found that AG Car is rotating very fast (up to 0.86 of its critical velocity), and
therefore it lies close to the so-called
limit (Maeder & Meynet 2000a). The AG Car
near-infrared spectrum shows prominent emission lines of Fe II and Mg II, and weak
emission of C I and N I. He I 10 667 and 10 830 Å show P Cygni profiles and
Pa
is seen in emission. The variability of the lines from 2001 (hot
phase) to 2004 (cool phase) is striking in this spectral region, specially in
He I 10 830 Å (Fig. 18). The line equivalent width dropped from 122 Å in 2001 to 0.2 Å in 2004. Also, in 2001 an additional
absorption component appeared in this line.
(Obj)P Cygni(Obj) (Fig. 10) is a famous member of the LBV class,
sharing its name to denote the characteristic shape of the spectral lines present
in those stars. P Cyg underwent a major eruption in the 1600 s, lying in a
quiescent phase since then. This eruption (and perhaps a previous one)
gave rise to complex ejecta around the star (Skinner et al. 1998; Chesneau et al. 2000; Exter et al. 2002).
The atmospheric parameters derived by
Najarro et al. (1997) revealed a helium-enriched photosphere and confirmed the
high mass-loss rate and luminosity of the star. A conference regarding the
latest P Cyg progresses was organized in 2000; we refer the reader to the
proceedings of the event (ASP Conf. Proceedings 233) for a general review.
The near-infrared spectrum of P Cygni is very similar to AG Car at minimum.
It shows strong He I 10 830 Å and Pa
emission, followed by more modest Fe
II and Mg II emission. From 2001 to 2003 the emission lines of P Cyg became
stronger, and since then they started to fade (Fig.
18).
(Obj)HD 326823(Obj) (Fig. 11) is located in the direction of the Galactic Center at a
distance of 2 kpc (McGregor et al. 1988). Sterken et al. (1995) suggested that
HD 326823 could be a LBV in a quiescent state, based on their photometric study.
(Lopes et al. 1992) presented the optical and the 1 m spectrum of this object. The
majority of the lines are double-peaked, probably indicating deviations from
spherical symmetry and suggesting the presence of a disk. The helium and nitrogen
lines are the strongest ones, while H
and H
are weakly present. Hence,
they concluded that HD 326823 is an evolved star, entering the Wolf-Rayet phase.
Borges Fernandes et al. (2001) present a high-resolution optical atlas of this star,
identifying the strongest features. The lines in the near-infrared spectrum show the same
double-peaked morphology as the optical ones. The strongest emission line seen is He I
10 830 Å, followed by Fe II and Mg II lines. Pa
and C I lines are barely seen.
Unfortunately, we observed HD 326823 only in the 2001 run, preventing us from
studying the variability of the star.
(Obj)GG Carinae(Obj) (Fig. 11) is suggested to be an eclipsing binary with a period
of 31 or 62 days (Gosset et al. 1985; Hernandez et al. 1981), although the
radial velocity curve is not well determined. The ultraviolet spectrum was
published by Brandi et al. (1987), who recognize numerous Fe II and Ni II
features. Our near infrared data from 2001 shows Fe II, Mg II, C I and He I 10 830 Å in emission. Pa
is the strongest line in our spectrum, showing a
blueshifted feature at 10 931.5 Å, that we could not identify, which may be a
component of Pa
formed in a high-velocity shell.
(Obj)HD 87643(Obj) (Fig. 11) is thought to be an evolved B[e] star (Oudmaijer et al. 1998)
and appears to be related to the LBV class, although the
link is not clear yet. This object has an optical spectrum dominated by
emission lines of Fe II and P Cygni profiles in the Balmer series, together
with low excitation forbidden lines. The spectrum is produced by a fast polar
wind combined with a slow disk wind (Oudmaijer et al. 1998). The ultraviolet
and optical spectra of HD 87643 was previously discussed by de Freitas Pacheco et al. (1985,1982),
who reported a strong spectral line variability.
HD 87643 has a bright reflection nebula, which was analyzed by Surdej et al. (1981) and
Surdej & Swings (1983).
Its 2001 near infrared spectrum shows prominent Fe II, CI lines and Pa emission, compatible with the presence of a cold wind. He I 10 830 Å
presents a P-Cygni profile, with a weak emission and a strong absorption that
goes up to -1750 km s-1, which is probably formed in the fast polar wind.
In Fig. 12 we present the OB supergiants of our sample, some of them which are related to
the LBV class. The stars are presented in decreasing order of spectral type. The spectra of the OB supergiants typically show He I 10 830 Å with a P Cyg profile and He I 10 667 Å, He I 10 913 Å and Pa
in absorption, with a variable veiling of their absorption depending on
how developed the wind is, culminating with HD 152408, that has only emission lines. The later subtypes
also present weak emission of Fe II 10 862 Å, Mg II 10 915 Å and Mg II 10 952 Å.
(Obj)HD 152408(Obj) and (Obj)HD 151804(Obj) (Fig. 12) are closely-associated
O supergiants in Sco OB1. The stellar wind of HD 152408 is stronger
despite the fact that the evolutionary
difference between both stars must be very small [see][]wf2000.
Both objects were analyzed by Crowther & Bohannan (1997) who derived
similar stellar parameters but a higher mass-loss rate and a
higher helium content in HD 152408 than in HD 151804, indicating the more evolved nature of the former.
The morphology of their near-infrared spectra agrees with the analysis of Crowther & Bohannan (1997),
with HD 152408 having strong emission of He I 10 830 Å,
He I 10913 Å and Pa ,
while HD 151804 shows the same morphology as the
hotter B supergiants stars. An unidentified feature is seen at 10 922.5 Å in
emission in HD 152408, while in HD 151804 it is in absorption.
(Obj)HD 154090(Obj) (Fig. 12) is
thought to be a normal B supergiant star, with an abundance typical of
early-type stars (Kane et al. 1981). Indeed, the near-infrared spectrum that
we present is similar to other B supergiants. However, there are
significant variations between the 2001 and the 2004 datasets, especially in He I 10 830 Å
(Fig. 18). In three years, the equivalent width of this
line reduced by a factor of two, with the P Cyg absorption almost disappearing.
These changes occurred together with radial velocity variations of 120 km s-1, and minor changes in the absorption profiles of Pa
and He I 10 913 Å.
This behavior encourages future monitoring of this star for better insight into its nature.
(Obj)HD 169454(Obj) (Fig. 12) is a member of the Sct OB3 association, located
at 1.5 kpc from the Sun (Humphreys 1978). However, the analysis of the
high-velocity interstellar gas in the direction of HD 169454 puts the star at
least at 3 kpc, which yields an absolute magnitude of
.
As a
consequence, HD 169454 may be one of the brightest B-type stars in the Galaxy. Its
optical spectrum was analyzed by Rivinius et al. (1997), who found a
variability pattern similar to that shown by HD 152236. We found that their near-infrared
spectrum and the changes from 2001 to 2004 are very similar as well (Fig. 16).
(Obj)HD 152236(Obj) (Fig. 12) is located in the Sco OB1 association, belonging
to the open cluster NGC 6231 (van Genderen et al. 1984). Photometric variability
in this object was reported by Sterken (1977) and Sterken et al. (1997),
who suggested the relationship between the cyclic charactheristics of the lightcurve of HD 152236
and the LBV class. These authors also claim an increase in the apparent visual magnitude through the
past centuries, with the star being 2 mag brighter in the 18th century. Burki et al. (1982)
analyzed the UV variability, showing the presence of discrete absorption
components (DACs). Rivinius et al. (1997) investigated the changes in the
optical spectrum revealing an intense variability during 1990-91. The
quantitative study of the spectrum by these authors also confirm that HD 152236
is a very luminous star (log
). Its near-infrared spectrum shows
He I 10 830 Å with a P Cyg profile with an absorption extending up to -450 km s-1 and He I
10 667 Å, He I 10 913 Å and Pa
in absorption. The comparison between
spectra taken in 2001 and 2004 shows radial
velocity changes of about 120 km s-1 and in the profile of He I 10 830 Å.
Also, the emission was stronger and the absorption was shallower in 2001 when compared
to the 2004 dataset (Fig. 16).
(Obj)HD 80077(Obj) (Fig. 12) is a very luminous B star (
)
with a
comparatively low mass-loss rate of
yr-1 (Carpay et al. 1991,1989)
among luminous B stars. These authors
claim that the star could be an LBV, since it has a low surface gravity and a
high luminosity. Historical observations also suggest that the mass-loss rate was
higher some decades ago (Houk 1978). The lightcurve presented by van Genderen et al. (1992)
shows that during 15 yr the star had an apparent
visual magnitude of
,
with micro-variations of about 0.05 mag. It is still
not known why HD 80077 is such a stable star, while its LBV neighbors in the HR diagram are
unstable. We observed this star only in 2001, and its near-infrared
spectrum is very similar to other luminous B stars, such as HD 152236 and HD 169454.
The only noticeable difference is the partial filling of the Pa
absorption and the presence of Mg II lines in emission. HD 80077 also displays
an unidentified absorption feature at 10 700 Å.
We have obtained spectra of 4 WN stars (Fig. 13) and 3 WC stars (Fig. 14). All of them are listed in the VIIth Catalog of Galactic Wolf-Rayet stars (van der Hucht 2001), where their properties are extensively listed. Physical properties of WR stars are thoroughly reviewed in Crowther (2007).
The spectra of the WN stars, namely (Obj)WR 6(Obj) (WN4), (Obj)WR 136(Obj) (WN6(h)),
(Obj)WR 78(Obj)
(WN7h) and (Obj)WR 22(Obj) (WN7h+O9III-V), show a very strong P Cygni profile in He I 10 830 Å blended with emission in Pa .
The WNE stars WR 6 and WR 136 show
flat-topped line emission profiles, while the later WN stars of our sample
show a more round-topped line emission profile. In WR 22 it is also noticeable
the He I 10 830 Å absorption from the O9 star.
Our sample of WC stars contains only the later subtypes ((Obj)WR 90(Obj) = WC7,
(Obj)WR 11(Obj) = WC8+O and (Obj)WR 103(Obj) = WC9d). This class of objects shows a broad, strong
emission of He I 10 830 Å with a small P Cyg absorption. There is also a
blended emission from C III + C IV centered at 10 545 Å. In the WC9d star
WR103 there is an emission of He I 10 913 Å that indicates its cool nature
compared to other WC stars. In WR 11 ( Vel)
the Pa
absorption from the secondary star (O7.5III-V) is also prominent.
In the previous section we presented unambiguous identification of most emission lines present in the spectral region covered by this work. The presence of C I lines in Be stars and LBVs is remarkable, as carbon is expected to be depleted in those objects, either due to rotational mixing or by the ejection of the outer layers. The near-infrared C I lines potentially can be used to constrain the carbon abundance in these stars.
![]() |
Figure 2:
Unidentified absorption lines
found in the LBVs HD 316285 (full), AG Car (dot-dashed) and ![]() |
Open with DEXTER |
![]() |
Figure 3: Same as Fig. 2, showing the presence of the same absorption lines in HD 160529, HD 168607, HD 168625 and W243. For clarity, the spectra were shifted up by an arbitrary offset. |
Open with DEXTER |
![]() |
Figure 4: Calibration of EW (10 780) versus E(B-V). The color excesses E(B-V) were taken from van Genderen (2001), except for W243 (Teodoro, private communication). A least-squares linear fit to the data is also shown (see Eq. (1)). |
Open with DEXTER |
However, the presence of unidentified absorption features is striking in most of the LBVs
at 10 780 Å and 10 792 Å, some of which do not have any photospheric absorption lines in
their spectra (e.g. Car). Figures 2 and 3 display a zoom around
those absorption lines. Their similar strength in stars spanning a range of spectral types,
together with their constance in variable stars such as LBVs, make it likely that the
unidentified absorption features arise in the interstellar or circumstellar medium of those
stars. The comparison between the EW of the absorption line found at 10 780 Å and the color
excess E(B-V) is displayed in Fig. 4. A correlation is clearly seen
between these quantities,
Compared to the feature at 10 780 Å, the absorption at 10 792 Å is weaker in
all objects except Car. Hence, it is difficult to perform a similar
calibration as the one made for the 10 780 Å feature. An observing campaign is underway
by our group to address this point. Nevertheless, this line is much stronger
in Eta Car than in other stars with similar interstellar reddening. It
was previously interpreted as a high-velocity (-1040 km s-1)
component of the stellar He I 10 830 Å (Smith 2002; Damineli et al. 1993). However, this line is present
in several other LBVs which do not present any high-velocity component, and it is constant
in
Car along its spectroscopic event. Therefore, we suggest that it is
formed in the circumstellar environment of
Car.
Indeed, if true, a precise calibration between EW(10 792 Å) and
E(B-V) may indicate the amount of circumstellar reddening present in
Car, which is still speculative. It would also provide constraints on
the luminosity and nature of
Car and other objects with unknown amounts of circumstellar reddening.
Acknowledgements
We are grateful to the referee Dr. J. Puls for construtive suggestions and detailed comments. J.H. Groh and A. Damineli thank Brazilian Agencies FAPESP (grants 02/11446-5 and 2005/51742-0) and CNPq (grant 200984/2004-7) for financial support. J.H. Groh acknowledges suppport from CNPq through its undergraduate research program (PIBIC). This work made extensive use of the SIMBAD database, which is operated by CDS. We thank Dr. Peter van Hoof for making available his atomic line database.
![]() |
Figure 5:
Spectra of the classical Be stars of our sample ordered by spectral type. Identification, spectral
type and projected rotational velocity, respectively, are indicated for each object. From top to bottom,
X Per, ![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
Same as Fig. 5, showing HD 120991, ![]() ![]() |
Open with DEXTER |
![]() |
Figure 7:
Same as Fig. 5, showing ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 8: Spectra of LBVs and LBV candidates. From top to bottom, HD 168625, HD 168607, HD 160529, HR Car, and W243. Spectra of variable stars also give year and month of the observation. |
Open with DEXTER |
![]() |
Figure 9:
Same as Fig. 8, showing ![]() |
Open with DEXTER |
![]() |
Figure 10: Same as Fig. 8, showing MWC 314, AG Car, and P Cyg. |
Open with DEXTER |
![]() |
Figure 11: Same as Fig. 8, showing HD 326823, GG Car, and HD 87643. |
Open with DEXTER |
![]() |
Figure 12: Spectra of the OB supergiants of our sample, showing from top to bottom HD 152408, HD 151804, HD 154090, HD 169454, HD 152236, and HD 80077. The spectra were offset by an arbitrary offset. |
Open with DEXTER |
![]() |
Figure 13: Spectra of the Wolf-Rayet stars from the WN subtype. From top to bottom, WR 6, WR 136, WR 78, and WR 22. |
Open with DEXTER |
![]() |
Figure 14: Same as Fig. 13, but for the Wolf-Rayets of the WC subtype. From top to bottom, WR 90, WR 11 and WR 103. |
Open with DEXTER |
![]() |
Figure 15: Spectra of non emission-line early-type stars, showing the typical photospheric lines present in the temperature range of our sample. |
Open with DEXTER |
![]() |
Figure 16: Spectroscopic variability of LBVs and related objects, from 2001 to 2004 (see legend). From left to right, AG Car, P Cyg and MWC 314. For each star, the upper panel presents a general view around He I 10 830 Å, while the lower panel shows a zoom around the P-Cyg absorption. |
Open with DEXTER |
![]() |
Figure 17:
Same as Fig. 16, but for He 3-519, HD 316285 and ![]() |
Open with DEXTER |
![]() |
Figure 18: Same as Fig. 16, but for HD 160529, HR Car, HD 154090, HD 168607, HD 152236 and HD 169454. Note the significant changes in radial velocity of HD 160529, HD 152236 and HD 169454. |
Open with DEXTER |