A&A 431, 615-622 (2005)
DOI: 10.1051/0004-6361:20041325
C. Aerts1,2 - K. Kolenberg1,3
1 - Instituut voor Sterrenkunde, Katholieke Universiteit Leuven,
Celestijnenlaan 200 B, 3001 Leuven, Belgium
2 -
Department of Astrophysics, University of Nijmegen, PO Box 9010, 6500 GL
Nijmegen, the Netherlands
3 - Institut für Astronomie,
Universität Wien, Türkenschanzstrasse 17, 1180 Wien, Austria
Received 20 May 2004 / Accepted 20 October 2004
Abstract
We have gathered and analysed multicolour Geneva photometry and
high-resolution spectroscopy of the stars HD 121190 (B9V) and HD 106419
(B9III) whose short-term periodic variability had become evident from their
HIPPARCOS data. We find three significant frequencies for HD 121190 in the
Geneva data: 2.6831, 2.6199 and 2.4713 c d-1 and classify the star as a
slowly pulsating B star, the coolest single star of that class known to
date. Its amplitude ratios are compatible with low-degree gravity modes. The
spectra reveal low-amplitude variability and
km s-1 which implies that the star rotates at 26% of its critical velocity.
For HD 106419 we are unable to confirm the period found from the HIPPARCOS
photometry in our ground-based data; rather we find a marginally significant
frequency of 0.8986 c d-1. This star is probably an evolved slowly
pulsating B star with complex variations. HD 106409 is also a
moderate rotator as the spectra reveal
km s-1,
which is also 26% of its critical velocity.
Key words: stars: oscillations - stars: variables: general - stars: individual: HD 121190 - stars: individual: HD 106419 - stars: early-type - line: profiles
One of the remarkable by-products of the ESA HIPPARCOS mission was the discovery of a large amount of periodically variable B stars (Waelkens et al. 1998). Among these, the large number of newly discovered slowly pulsating B stars (hereafter abbreviated as SPBs) was quite unexpected. Large follow-up studies of these stars were initiated by Aerts et al. (1999) and Mathias et al. (2001). This resulted in a sample of some 25 well-known members of this class of multiperiodic gravity-mode oscillators, the statistical properties of which were summarized by De Cat (2002) and seem to fulfill all expectations from theoretical studies (Pamyatnykh 1999).
Among the known SPBs, the double-lined spectroscopic binary HD 123515 stands out as the coolest one with an effective temperature of 12 000 K and spectral type B9IV (De Cat et al. 2000; De Cat & Aerts 2002). All other members are at least 1 500 K hotter (De Cat 2002). Given that most of the SPBs known nowadays were discovered from the HIPPARCOS unbiased sample (Waelkens et al. 1998), we must conclude that the lower SPB instability strip is less populated than its upper counterpart and one may question if single B8/9 SPBs exist at all. A similar question, albeit from a somewhat different perspective, was already tackled by Baade (1989a,b) many years before the results of HIPPARCOS became available. He performed an extensive spectroscopic search for line-profile variability among 22 late-B main-sequence stars and subgiants and came to the same null result: no B8/9 IV/V line-profile variable stars were found. B8/9 stars are also absent in the list of Be line-profile variables analysed and interpreted recently in terms of non-radial oscillations by Rivinius et al. (2003).
![]() |
Figure 1:
The selected SiII
![]() |
Open with DEXTER |
As a continuation of the quest for oscillating late-B stars, several stars
caught our attention as they were listed in the HIPPARCOS catalogue of periodic
variables with periods of only a few hours. The current paper concerns two of
these stars of spectral type B9. HD 121190 (B9V,
)
was already
classified as a new candidate SPB by Waelkens et al. (1998), although the two
classification schemes used by these authors did not give consistent results for
this star. Waelkens et al. (1998) found one period of 0.38 d in the HIPPARCOS
photometry for this star. All other confirmed SPBs have longer periods (De Cat
& Aerts 2002; De Cat 2002). Moreover, theoretical models typically predict
periods above half a day (Pamyatnykh 1999), although this lower limit is not
strict as it depends largely on metallicity and rotational effects, the latter
of which are usually neglected in theoretical SPB excitaton calculations so
far. A second star that struck us is HD 106419 (B9III,
). The
HIPPARCOS team assigned an intrinsic period of 0.11 d to it. If confirmed, such
a short period would imply very severe deviation from current theoretical model
predictions.
The current paper describes the results of our follow-up campaigns dedicated to HD 106419 and HD 121190. In Sect. 2 we discuss the data that were gathered with the goal to confirm their periodic variability. Section 3 deals with the determination of the basic physical parameters of these stars and their position with respect to the SPB instability strip. A detailed frequency analysis is described in Sect. 4. Finally, we discuss our findings in Sect. 5.
We included HD 106419 and HD 121190 in our long-term photometric monitoring programme of pulsating B stars of Leuven University. In this framework we obtained Geneva 7-colour high-precision photometry for the two stars during one campaign with a 3-week duration and one campaign of 9 weeks within one season. The integration times were typically 5 min for HD 106419 and 3 min for HD 121190 resulting in a precision of about 5 mmag per measurement. In total we obtained 165 measurements for HD 106419 with a time span of 104 days and in 145 datapoints for HD 121190 with a 115 days time span. For HD 121190 we had already 5 very old Geneva measurements available, which result in an overall time span of 6922 days. The old measurements turned out to be fully compatible with the new ones, so we kept them in our analysis.
Given the confirmed photometric variability of the stars (see Sect. 4), we obtained additional spectroscopic data with the FEROS spectrograph attached to the ESO 2.2 m telescope. The data were taken during 5 consecutive nights 2 years after the photometric campaigns. The main goal of this spectroscopic study was to exclude binarity, to check for evidence of line-profile variability, as well as to derive an accurate estimate of the rotational velocity of the stars.
The line-profile variability of the slow rotators among the SPBs is by far best
studied from the SiII doublet at
Å (see, e.g., Aerts et al. 1999).
However, for our two target stars, this doublet turned out to be
blended due to their relatively high rotation rate compared to other SPBs
and so we focused on the much better
isolated SiII
Å line. In order not to miss any line-profile
variability with short periods, we decided to integrate during 15 min for
both stars. This resulted in a signal-to-noise ratio between 80 and 150 for
HD 106419 and between 150 and 300 for HD 121190 in the region near the
SiII
Å line. We show all the data of this line, for both
target stars, in Fig. 1.
Table 1: Logbook of the dedicated measurements of HD 106419 and HD 121190. The Heliocentric Julian Dates (HJD) are given in days with respect to 2 450 000 and N denotes the number of measurements.
A logbook of all the data is provided in Table 1.
We used the method outlined by Künzli et al. (1997) to estimate the mean effective temperature and gravity from the average Geneva colours. The results are provided in Table 2 and are in very good agreement with the spectral types B9III (HD 106419) and B9V (HD 121190) listed in Simbad. We point out that the error estimates for the temperature and gravity resulting from the code by Künzli et al. are internal errors only, in the sense that they result from interpolation in the grids of standard stars. They assume the average Geneva colours and the reference tables to be error-free. The true uncertainties are typically a factor two larger and it is these more realistic error estimates that we provide in Table 2.
Subsequently, we used the photometric calibration outlined in Heynderickx et al. (1994) together with interpolation in the evolutionary models published by Schaller et al. (1992) to estimate the mass, radius, luminosity and bolometric magnitude of each star. The result of the physical parameters obtained in this way are summarised in Table 2 and are used to place the two stars in the HR diagram in Fig. 2. The SPB instability strip calculated by Pamyatnykh (1999) is also indicated, as well as some evolutionary models. The instability strip is not considered beyond the TAMS. It can be seen that HD 106419 is clearly an evolved star. HD 121190, on the other hand, is a main-sequence star situated just at the lower limit of the instability strip, which was computed for solar metallicity and ignoring core overshooting as well as rotation.
Table 2:
The physical parameters of the two programme stars. The mass and
radius are given in solar units;
and
are given in
km s-1.
The thermal broadening of metallic lines in B9 stars amounts to only a few
km s-1 and is negligible compared to the measured width of the lines. For
this reason, and keeping in mind the variable nature of the stars, we estimated
their
by averaging the full-width-at-half-maximum of all the
SiII
Å lines shown in Fig. 1, i.e.,
.
This resulted in 78 and 118 km s-1for respectively HD 106419 and HD 121190. The projected rotation velocity of
HD 106419 equals the mean value found for 70 B9-B9.5 III stars (Abt et al.
2002). The mean
-value found by Abt et al. (2002) for 145 B9-B9.5 V
stars amounts to 134 km s-1 which is only slightly higher than
HD 121190's value. Both our target stars are therefore "normal'' rotators
among the B-type stars.
Taking the values for the mass and radius given in Table 2 we
determined the critical rotation velocity for both stars from the
well-known expression valid within the Roche approximation for a centrally
condensed star:
(e.g. Tassoul 1978).
It is found that both HD 106419 and HD 121190 rotate at at least 26% of
their critical velocity. At this low rate, the equatorial radius of the star is
only a factor 1.03 larger than the polar radius such that the oblateness of the
star due to the centrifugal force is negligible.
All but one of the confirmed SPBs rotate much slower than 26% of their critical
velocity (De Cat 2002), HD 1976 being a notable exception with 32% (Mathias et al. 2001). From this perspective, both our target stars rotate much faster than
the average SPB.
![]() |
Figure 2:
The position of the two stars (crosses) with respect to the SPB
instability strip which is taken from Pamyatnykh (1999) and is based on OPAL
G93/21 opacities (Iglesias & Rogers 1996) for a composition
X, Z =0.70, 0.02.
Evolutionary tracks for 3, 4, 5 and 6 ![]() |
Open with DEXTER |
We searched for frequencies in the data with the method outlined in Scargle
(1982). The accuracy of the frequencies was calculated as
where the proportionality constant is of order
1, A is the amplitude of the frequency f, N is the number of data points,
T is the total time base and
is the error estimate of the individual
measurements, which is 5 mmag for the Geneva photometry (Cuypers 1987).
After
having derived the frequencies by a prewhitening process, we determined the
signal-to-noise (S/N) ratio in the periodogram of the residuals by averaging the
amplitude of all peaks for the frequency range [0,20] c d-1. We use the
4 S/N value as a criterion to judge a frequency to be significant or not,
following Breger et al. (1993).
![]() |
Figure 3: The periodograms for the HIPPARCOS ( upper) and Geneva U ( lower) data of HD 106419. The dashed horizontal line indicates the 4 S/N ratio determined in the way explained in the text. |
Open with DEXTER |
In the Geneva U data, for which the periodogram is shown in the lower part of
Fig. 3, we do find one marginally significant frequency:
c d-1. The amplitude of this frequency amounts to
0.0051
0.0013 mag and its variance reduction is 15%. The value of this
frequency is typical for the gravity modes expected to occur in SPBs, unlike the
high frequency determined by the HIPPARCOS team. Nevertheless, it is only
slightly above the 4 S/N level. Moreover, we do not find any clear evidence of
fU, or any other frequency, in the data for the other 6 Geneva filters. From
this we cannot but consider fU to be a candidate frequency only as long as we
have no independent evidence of its reality. We show the phase diagrams for the
HIPPARCOS and Geneva U data, for respectively
and fU, in
Fig. 4. These phase diagrams are typical for multiperiodic
stars with low amplitudes (see De Cat & Aerts 2002, for additional examples).
The SiII
Å profiles shown in Fig. 1, while rather
noisy, show small-amplitude variability of the line center. This is independent
proof that the star is variable, but we have too few spectra to do an
independent frequency search for these data.
Our conclusion for this star is that it is probably an evolved SPB with a complex low-amplitude frequency pattern, covering a wide range in frequency. This suggestion needs further confirmation, however.
![]() |
Figure 4: The phase diagrams for the Geneva U ( left) and HIPPARCOS ( right) data folded with different frequencies, of which the values are given in the text. The open dots are the data and the full line is a least-squares fit for the considered frequency. |
Open with DEXTER |
The Geneva data for HD 121190 are completely in line with the HIPPARCOS
results. Moreover, they clearly point towards a multiperiodic star. In
Fig. 5 we summarize the results for the Geneva U filter, in
which the amplitudes are largest. We find evidence for three significant
frequencies:
c d-1,
c d-1 and
c d-1. After
prewhitening with
f1, f2, f3, the frequency 2.8452 c d-1 occurs, but
it does no longer fulfill the 4 S/N criterion. The frequency f2 corresponds
to the one found in the HIPPARCOS photometry. The phase plots for the U filter
for subsequent stages of prewhitening are shown in the upper panel of
Fig. 6. In that same plot we also show the phase diagram for f1 for the colours U-B, B-V as well as the HIPPARCOS data folded with
f2. The amplitudes and phases in the 7 filters for the significant
frequencies are provided in Table 3, as well as the variance
reduction for a triperiodic least-squares fit with
f1, f2, f3 fixed.
![]() |
Figure 5: The periodograms for the Geneva U data of HD 121190 after subsequent stages of prewhitening with the frequencies provided in the text. The amplitudes are expressed in magnitudes. The dashed horizontal line indicates the 4 S/N ratio determined in the way explained in the text. |
Open with DEXTER |
From the amplitudes listed in Table 3 we derive the ratios with respect to the U filter shown in Fig. 7. Although the error bars are large, it can be seen that a clear decrease of the amplitude from U to B occurs. This result is typically found from observations for all confirmed SPBs (De Cat & Aerts 2002; De Cat et al. 2004). It is also what is expected from theoretical computations (Townsend 2002; Dupret et al. 2003). The observed decrease in amplitude from U to G is similar for all three modes, although the main mode's decrease is somewhat smaller. In general, this points towards modes with similar degree.
A more quantitative statement on the degree of the oscillation modes can in
principle be derived from the confrontation of the observed ratios with those
predicted from model calculations. Townsend (2002) has provided a complete
systematic study of the photometric behaviour expected for SPBs in the
approximation of non-rotating stars. In his Fig. 9, he provides the amplitude
and phase differences predicted by models for all unstable
modes in the SPB instability strip as well as the observed values found by De Cat & Aerts (2002). For mode identification in specific stars, however, one has
to consider only the ratios and phase differences for models with appropriate
physical parameters. The top rows in Figs. 7 and 8 of Townsend (2002) are
suitable to identify the three detected modes of HD 121190. First of all we
point out that none of the
and 4 modes are found to be unstable for the
appropriate mass. This result is confirmed by us on the basis of the entirely
independently written evolution and pulsation codes CLÉS and MAD (see Dupret
et al. 2003). This leaves us with the possibility of
or 2, if
we take the theoretical excitation predictions based on the
mechanism
applied to models with a solar mixture for granted - see, however, Pamyatnykh
et al. (2004) and Ausseloos et al. (2004) for an example where such
predictions fail, albeit for a pressure mode in the much hotter
Cep
star
Eri. The confrontation of our observed amplitude ratios and very
small phase differences with the theoretical predictions by Townsend (2002)
leads to the conclusion that the dominant mode must be
.
We are unable
to discriminate between
and 2 for the two additional modes.
![]() |
Figure 6: Upper panels: the phase diagrams for the Geneva U data of HD 121190 for the three significant frequencies f1, f2, f3 after subsequent stages of prewhitening. Lower panels: left and middle: the Geneva colours U-B and B-V folded with the main frequency; right: the HIPPARCOS data folded with f2. The open dots are the data and the full line is a least-squares fit for the indicated frequency. |
Open with DEXTER |
An important remark that has to be made regarding the achieved mode
identification is the fact that the Coriolis force plays a significant role in
the behaviour of gravity modes of low frequency whenever the
pulsational frequency in the corotating frame becomes larger than half the
rotational frequency (Townsend 2003, and references therein). In such a
situation, neglect of the Coriolis force in mode identification from amplitude
ratios is inappropriate (see Fig. 5 in Townsend 2003). We have, unfortunately,
no means to evaluate the importance of the Coriolis force in the treatment of
HD 121190's oscillations because we have no information about its inclination
angle i. In any case, from the radius estimate listed in Table 2 and
the measured ,
we deduce
c d-1, which is
about half of the measured oscillation frequencies. It is therefore quite likely
that the frequencies observed for HD 121190 are significantly shifted with
respect to those in a non-rotating star due to the Coriolis force, unless we
see the star almost equator-on. In the latter case our mode identification is
still accurate, but for smaller inclination angles our results have to be
treated with caution.
![]() |
Figure 7:
Amplitude ratios AX/AU for the three detected modes of HD 121190
with their errors, where AX stands for any of the amplitudes in the 6 Geneva
filters B1, B, B2, V1, V, G. The symbols are: ![]() ![]() ![]() |
Open with DEXTER |
Table 3:
Results of harmonic fits to the Geneva lightcurves of HD 121190. A
stands for the amplitude, expressed in millimag, and
for the phase,
expressed in 2
radians.
The adopted reference epoch for
corresponds to HJD 2 450 000.0.
![]() |
Figure 8: The phase diagram of the radial velocity of HD 121190 derived from the profiles shown in Fig. 1. The open dots are the observations and the full line is a fit for the frequency f1. |
Open with DEXTER |
Due to the moderate
rotation of HD 121190 it is not obvious from the lines shown
in Fig. 1 whether profile variability occurs and if so, whether it is
consistent with the photometric variability. In order to check that, we
computed line diagnostics for the 40 profiles shown in Fig. 1. We
considered the equivalent width, the radial velocity and the first three
velocity moments in the definition by Aerts et al. (1992). The time span of the
spectroscopic data is only 5 days and the frequency accuracy is only at
0.05 c d-1 level. The periodograms of the radial velocity and the first
and second moment show a maximum near 0.7 c d-1 and all of its one-day
aliases and so they are consistent with the star being a line-profile variable
with frequency f1 found in the multicolour photometry (f1 and f2 are
hardly resolved in the spectroscopic data). A phase plot of the radial velocity
variation for f1 is shown in Fig. 8. We find an average radial
velocity of 21.40.5 km s-1 which is entirely compatible with the
much less accurate value of 23 km s-1listed in Simbad.
The amplitude of the fit to the radial velocity
variation for f1 shown in Fig. 8 is 2.8
1.0 km s-1. The
peak-to-peak variation in the radial velocity amounts to
18 km s-1.
No
obvious periodicity is found in the equivalent width and the
full-width-at-half-maximum.
We conclude to have found clear evidence for multiperiodic variations in the 7-colour Geneva photometry of HD 121190. Three significant frequencies are present. Each of them has amplitude ratios compatible with low-degree non-radial gravity-mode oscillations. The spectra are consistent with this photometric variability. In view of its physical parameters we therefore classify HD 121190 as the coolest single and second most rapidly rotating SPB known to date.
We have provided clear observational evidence for the occurrence of line-profile variability in two stars of spectral type B9. For one of these, HD 121190, we were able to detect three significant frequencies near 2.6 c d-1 in the multicolour photometry. The variability is consistent with high-order low-degree gravity-mode oscillations. For HD 106419 we were unable to unravel the complex low-amplitude frequency spectrum.
During the last half century, there has been quite some confusion and debate in
the literature about the existence of a specific group of variables with
spectral types between B7 V-III and A2 V-II and periods between 2 and 8 h,
baptised "Maia stars'' by Struve (1955). We point out that these hypothetical
stars would be partly situated within the SPB instability strip, extending
towards the Scuti strip. Despite large search campaigns (see, e.g.,
Scholz et al. 1998, and references therein) clear detection of short-period
variability for these stellar types remains ambiguous for most candidates. At
the time when Struve (1955) made his suggestion, SPBs were still unknown.
Applying the Struve (1955) criteria, HD 121190 would be classified as a Maia
star. However, we object this type of classification. Indeed, ever since the
observational studies by Waelkens (1991) and Waelkens et al. (1998) on the one
hand, and the theoretical works by Dziembowski et al. (1993) and Gautschy &
Saio (1993) on the other hand, it does not make sense any longer to speak of
Maia stars for spectral types B. There is one clear physical mechanism that
explains such oscillating stars:
driven, non-radial, gravity-mode
oscillations, and so these objects should be termed SPBs. The short periods are
sometimes claimed to be problematic from a theoretical viewpoint (e.g. Kallinger et al. 2002) but we would like to offer a different view. From
their extensive study, Scholz et al. (1998) conclude that only 3 stars out of
the 15 for which long-term observational studies are available are confirmed as
periodic variables:
UMi (A2III),
CrB (A0V) and ET And
(B9IV). These three stars have a
of respectively 165, 112, and
80 km s-1. Hence they rotate much faster than the average SPB, just as
our two target stars. As already discussed in Sect. 4.2, such rotation
velocities imply that the effects of the Coriolis force may come into play and
that this force introduces significant frequency shifts for the low-frequency
gravity modes (Townsend 2003). We believe that such shifts offer the correct
explanation for the relatively high observed frequency values for the modes that
we have detected in HD 121190 and that Scholz et al. (1998) found for the
three stars mentioned above. We also recall that theoretical instability strips
so far have always been provided for non-rotating stars (e.g. Pamyatnykh 1999
and Fig. 2). The strip shown in Fig. 2 would no longer be
placed in exactly the same position, and the periods of the oscillations would
have a much broader allowed range if the Coriolis force would be taken into
account in an accurate way.
As for the observational search for pulsating late-B to early-A type stars, we
warn that large high-precision multicolour photometric and/or high-resolution
spectroscopic campaigns have to be undertaken before making any firm decisions
about the presence or absence of low-amplitude multiperiodic oscillations, as
also already emphasised by Scholz et al. (1998). One should not make firm
conclusions based on short observing runs. Rather, dedicated campaigns of the
type described in Aerts et al. (1999), De Cat et al. (2000), Mathias et al. (2001) and De Cat & Aerts (2002) are needed to study such oscillations. The
same is true for the gravity modes in the much cooler Doradus stars
(see, e.g., Poretti et al. 2002; Henry & Fekel 2003; Aerts et al. 2004;
Mathias et al. 2004 for some recent examples of such long-term campaigns).
Acknowledgements
We are much indebted to the referee, Dr. R. Townsend, for having pointed out some serious shortcomings in an earlier version of our paper. We thank our colleagues from the Institute of Astronomy of Leuven University who contributed to the gathering of the photometric data: Peter De Cat, Joris De Ridder and Gwendolyn Meeus.