A&A 452, 969-974 (2006)
DOI: 10.1051/0004-6361:20054635
R. S. Schnerr1 - E. Verdugo2 - H. F. Henrichs1 - C. Neiner3,4
1 - Astronomical Institute "Anton Pannekoek'', University of Amsterdam,
Kruislaan 403, 1098 SJ Amsterdam, The Netherlands
2 - European Space Astronomy Centre (ESAC), Research & Scientific
Support Department of ESA, Villafranca del Castillo, PO Box 50727, 28080 Madrid, Spain
3 - Institute for astronomy, KU Leuven, 3001 Leuven, Belgium
4 - GEPI, UMR 8111 du CNRS, Observatoire de Paris-Meudon, 92195 Meudon Cedex, France
Received 2 December 2005 / Accepted 7 March 2006
Abstract
We report on attempts to measure the magnetic field of the pulsating B star
Eridani with the Musicos spectropolarimeter attached to the 2 m telescope at the Pic du Midi,
France. This object is one of the most extensively studied stars for pulsation modes, and the existence of a magnetic field was suggested from the inequality of the frequency separations of a triplet in the stars' oscillation spectrum. We show that the inferred 5-10 kG field was not present during our observations, which cover about one year. We discuss the influence of the strong pulsations on the analysis of the magnetic field strength and set an upper limit to the effective longitudinal field strength and to the field strength for a dipolar configuration. We also find that the observed wind line variability is caused by the pulsations.
Key words: stars: magnetic fields - stars: early-type - stars: individual:
Eri - stars: oscillations - stars: activity - line: profiles
The B2III star
Eridani (HD 29248, V = 3.93) is known to show radial velocity variations for more than a century (Frost & Adams 1903). It was found to be a multi-mode non-radial pulsator belonging to the class of
Cephei variables with a main frequency of 5.76 c d-1, identified as a
,
p1 mode.
Handler et al. (2004) and Jerzykiewicz et al. (2005) detected two independent low frequency, high-order g modes, demonstrating that the star also belongs to the class of Slowly Pulsating B (SPB) stars.
The star has the richest known oscillation spectrum of all
Cephei stars. From a very extensive campaign (see Aerts et al. 2004; Handler et al. 2004; De Ridder et al. 2004; Jerzykiewicz et al. 2005, hereafter Papers I, II, III and IV), 34 photometric and 20 spectroscopic frequencies were detected, corresponding to 14 different pulsation frequencies.
Among these 14, 12 are high-frequency modes, out of which 9 form three triplets, which are slightly asymmetric. The symmetric part is attributed to the effect of stellar rotation, whereas the asymmetric parts could be due to higher order rotational effects or to a magnetic field. From line profile modeling Smith (1983) derived
km s-1 for the projected equatorial velocity. This value is consistent with a rotation period of 30-60 days as derived from the modeling of the splitting of the strongest triplet around 5.64 c d-1, consisting of
modes (Dziembowski & Jerzykiewicz 2003, hereafter DJ; Paper I; Paper II).
DJ found that the asymmetry of this triplet, as measured from data of van Hoof (1961), could only partly be explained by the quadratic effects of rotation (see, e.g., Saio 1981) and suggested that a strong magnetic dipole field of the order of 5-10 kG could explain this discrepancy. In a more recent analysis (Paper I) the asymmetry was found to be a factor of 2 smaller than before, and Pamyatnykh et al. (2004) could entirely account for the asymmetry in terms of quadratic rotational effects. However, in Papers III and IV the asymmetry was again found to be larger, and a second and third triplet were detected around 6.24 c d-1 and 7.91 c d-1. More advanced modeling of the different separations and asymmetries in all three triplets is still needed.
It is clear that an observational limit for the magnetic field strength will constrain such models, but until now no magnetic measurements of this star are available. The specific prediction by DJ motivated us to observe this star with the Musicos spectropolarimeter at the Télescope Bernard Lyot (TBL, Pic du Midi, France) to search for the presence of a magnetic field, as the detection limit of this instrument is of the order of 100 G, which is far below the predicted value. An additional argument to search for a magnetic field was the observed stellar wind variability as recorded 25 years earlier by the International Ultraviolet Explorer (IUE) satellite.
This paper describes our attempts to measure the magnetic field of Eri from spectropolarimetric observations. We show the UV line variability as observed by the IUE satellite (Sect. 2.1), describe how we interpret the signatures in our spectropolarimetric measurements as entirely due to the strong pulsations in this star rather than due to a magnetic field (Sect. 2.3) and set an upper limit to the field strength (Sect. 3).
Specific behaviour of variable stellar wind lines belongs to the well-known indirect indicators of a magnetic field in early-type stars (see Henrichs et al. 2003, for a review). In Fig. 1 we show the UV line profile changes in
Eri, as recorded with the IUE satellite in 1979.
The epochs of the observations and the radial velocity variations due to pulsation, which have been used to correct the spectra, are given in Table 1.
![]() |
Figure 1:
Ultraviolet C IV line profile variability of ![]() ![]() ![]() ![]() |
Open with DEXTER |
Table 1: Epochs of the IUE observations with applied radial velocity corrections and continuum ratios used to normalise the spectra (see text). The calculation of the radial velocities is based on Paper III.
The temporal variance measures the ratio of the observed to the expected variability, and is very similar for the C IV wind profiles in all magnetic B stars. We calculated the expected variability as a function of wavelength with a noise model as described by Henrichs et al. (1994). The temporal variance spectrum of
Eri (lower panel in Fig. 1) is very similar to that of a magnetic oblique rotator, such as
Cas (a B2IV star, see Neiner et al. 2003a). Although only 10 high-resolution IUE spectra of
Eri exist, nine of which were taken within one day and the other one month later, the peaks around +100 km s-1 in both doublet members are significant at the 3
level.
![]() |
Figure 2: UV continuum flux near C IV measured with IUE (dots - typical errors are 0.002) and the normalised U-band flux (line) deduced from the photometric pulsation analysis in Papers I and III; the typical error range is indicated by the dashed lines. The timescales of the variability observed in the UV and in the U-band are very similar. |
Open with DEXTER |
Although this variability is suggestive of the presence of a magnetic field in Eri, most of the variation is observed over one day. This period is consistent with known pulsation modes, which have periods between 0.126-2.312 days, rather than with the rotation period, expected to be of the order of months (see also Sect. 3.1).
Before the temporal variance is constructed, the 10 spectra are normalised such that the equivalent widths summed over the wavelength bands [1465, 1510] Å and [1575, 1605] Å are equal to their average value. These regions were selected to be free from stellar wind affected lines. The normalisation is necessary because of the intensity variations due to the pulsations (see, e.g., Porri et al. 1994, Paper I). In Fig. 2 we show the inverse normalisation constants, which can be considered as a measure of the UV flux for the 9 spectra of Feb. 23/24, 1979. The main timescale seen in the light curve is the same as found in several optical bands in Paper I, and corresponds to that of the strongest pulsation mode detected in this star.
Table 2:
Journal of TBL observations, with epochs, exposure times and a comparison (O-C) between the
measured radial velocities (
)
and predicted radial velocities (
)
based on the ephemeris and amplitudes of
Paper III. Note that these values correspond reasonably well, provided that a constant difference of
km s-1 for the system velocity is taken into account.
The magnetic field measurements were carried out with the Musicos spectropolarimeter attached to the 2 m TBL at the Pic du Midi, France. We obtained 32 spectra of Eri between 8 Feb. 2003 and 14 Feb. 2004 (see Table 2) from which circularly polarised (Stokes V) and unpolarised spectra (Stokes I) are calculated.
The technique to carry out
high-precision magnetic measurements with this instrument is extensively described by
Donati et al. (1997) and Wade et al. (2000). Each set of four subexposures was taken in the usual
-plate position sequence q1, q3, q3, q1, corresponding to
angles.
We used the dedicated ESpRIT data reduction package (Donati et al. 1997) for the optimal extraction of the échelle spectra and to obtain Stokes I and V spectra.
We also calculate a Null polarisation, called Stokes N, which represent the pollution by non-magnetic effects and should be null for a perfect measurement.
The package includes a Least-Squares Deconvolution (LSD) routine to calculate a normalised average Stokes I line profile and corresponding Stokes V and N line profiles of all available spectral lines (we used 107-108 lines). If a magnetic field is present it will result in a typical Zeeman signature in the average Stokes V profile, from which the effective longitudinal component of the stellar magnetic field can be determined (see Sect. 3.2).
Stokes V and N spectra were calculated using the standard equations:
![]() |
(1) |
Although no significant magnetic fields are detected (see Table 3), in several cases the average Stokes V profiles (see Figs. 3 and 4) show a significant signature which could be interpreted as due to a magnetic field. However, significant signatures are also visible in the corresponding Stokes N profile, which makes a magnetic interpretation questionable, especially because of the large changes in the line profile due to pulsations during the 4 subexposures.
![]() |
Figure 3:
LSD results for ![]() |
Open with DEXTER |
We now consider the origin of the signatures. In an ideal instrumental setup, pulsations would not create signatures in Stokes V (or N) since for each spectrum both right-handed circular and left-handed circular polarisation spectra are recorded simultaneously, and different line shapes between the four subexposures of one magnetic field measurement should cancel exactly (see Eq. (2)). The practical reason why they do not cancel is that the two spectra of one subexposure partially follow a different light path and are recorded on different pixels of the CCD. This is why usually at least 2 (and often 4) subexposures are used. As a consequence the two spectra may have a different intensity level and slight differences in wavelength calibration, on average typically several hundred m s-1 (Donati et al. 1997; Semel et al. 1993). Such inevitable inaccuracies cause the different line shapes between the different subexposures to appear in the resulting Stokes V spectrum. For stars that do not show strong changes in line shape this problem does not occur, because the differences between the two spectra of one subexposure are corrected by the next subexposure where the opposite circular polarization state is recorded through the same light path, on the same pixels and with the same wavelength calibration. In the following section we closely examine the effect of these differences.
![]() |
Figure 4: Stokes I, N and V profiles of observations 1 ( left), 2 ( middle) and 6 ( right). These profiles were selected to illustrate the effect of the pulsation. In the Stokes I plots ( top), the solid lines represent the model fits. The dashed lines are the four exposures that were used for one magnetic field measurement. In the plots of the Stokes N ( middle) and V ( bottom, in ), the dashed line represents the observations, the thin solid line the model, and the thick solid line the corrected observations after subtracting the model to correct for the effect of the pulsations. Note that the corrected Stokes N profiles are all consistent with zero, which indicates that the assumption of a constant velocity shift between the two beams is sufficient to explain the signatures in Stokes N for pulsation-affected profiles. The limits of [-0.04%, +0.04%] adopted in Sect. 3.3 as the maximum amplitude of any undetected, broad, magnetic polarisation signatures, are indicated by the dashed lines in the bottom left plot. |
Open with DEXTER |
Table 3:
Results of the TBL magnetic field measurements. The signal to noise ratio per pixel ( S/N) was measured around 550 nm in the Stokes V spectrum (order 108). The HJD was calculated halfway the four subexposures used for each magnetic measurement. Measurements shown are the effective magnetic field as measured from the Stokes V profiles before and after correcting for the pulsations (
and
respectively) and similarly for the values for N. The last column gives the best-fit velocity shift between the two beams.
To investigate the effect of pulsations on the signatures in the Stokes V profile we developed a model that predicts Stokes V and N signatures from the Stokes I profiles of all four subexposures. For each subexposure we first fit the average line profile, resulting from the LSD method, with the following function:
To model the error in the wavelength calibration between the two beams of one subexposure, we characterise the typical wavelength shift between the two beams with a velocity parameter
.
From the resulting 8 line profiles (4 subexposures
2 beams) we calculate Stokes V and N profiles, with a fixed value for
to be determined by the fit. It is important to note that such an offset in wavelength calibration has a different effect than the presence of a magnetic field. A magnetic field would cause the spectrum in one circular polarisation state to be shifted relative to the other due to the Zeeman splitting; this shift is of opposite sign in the q1 and q3 exposures due to the switching of the two polarisation states between the subsequent subexposures. In our case the shift is the same in both q1 and q3 subexposures because this parameter is related to the beams themselves.
We determine the parameter
by a minimum
fit of the calculated Stokes N profiles from our model to the measured Stokes N profiles. The results can be found in Table 3. To check whether these values are realistic, we have extracted several ThAr exposures in the same way as we do for our science exposures. Since ThAr exposures are used to wavelength calibrate the spectra, we would expect similar inaccuracies between the two beams in these exposures as for our science exposures. We indeed find shifts in velocity between spectral lines in the two beams, varying from
to +1 km s-1, explaining the average shifts found.
In Fig. 4 we show the Stokes I, N and V profiles of the three observations which were subjected to the strongest pulsations. It is clear that the features in the Stokes N profile can be fairly well reproduced with our simple model, which contains only one free parameter. The signatures in the Stokes V that are predicted by our model have a somewhat smaller amplitude and appear to be slightly broader than the observed profiles, but the overall shape is very similar. The broadness of the modeled profiles is probably due to the fact that we use only one parameter,
,
to characterise the shift of the average line profile, while in reality there is a distribution for all the different lines. Furthermore, at some phases the line profiles show extended wings which are not represented in our model, and resulted in slightly broader model fits.
With this modeled contribution to the Stokes V profile we can quantitatively determine the spurious effect on the magnetic field determination as will be done in Sect. 3.2.
In Table 3 we show the effective magnetic field strength as measured from the observations both before and after subtracting the modeled signature of the pulsations. The integration ranges used are set at two times the width of the fitted line profile as determined from Eq. (3) (see also Table 3). In general, the lower and upper limits of the integration are different due to the varying shape of the line profile. The difference in
and
with and without the correction for pulsations is quite small. This is because the signatures created by the pulsations are almost symmetric, and hence do not much influence the first moments that determine these quantities.
No magnetic field has been detected. However, for some observations significant signatures in Stokes V are found. This could indicate that although the effective longitudinal magnetic field strength is zero, there is still evidence for the presence of a magnetic field. Using our simple model, we have shown that we are able to model the signatures found in Stokes N very well, and at the same time predict Stokes V profiles that are very similar in shape to the measured ones. From this we conclude that the profiles we detect in Stokes N and V are the result of the combined effects of the pulsations and inaccuracies in wavelength calibration that were not removed by our imperfect modeling of these effects.
![]() |
Figure 5:
For a given magnetic field strength, the amplitude of the Stokes V signature depends mainly on ![]() |
Open with DEXTER |
![]() |
Figure 6:
Results of model calculations of Stokes V profiles for a star with a dipolar magnetic field with a polar strength of 300 G. The left plot shows profiles for an inclination of ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Constraining the magnetic field of a star from Stokes V profiles is not straightforward. A Stokes V profile is only sensitive to the line-of-sight component of the field and a light-intensity weighted average over the visible stellar disc is involved, similar to the formation of a line profile. Since one also has to account for the rotational Doppler shifts of the lines, it will be clear that Stokes V profiles can be very different for stars with the same (polar) field strength
and
,
even for a simple dipolar configuration.
To determine the constraints on the polar field strength from our spectropolarimetry we used a model that calculates Stokes V profiles from the Zeeman splitting of an absorption line. In this simple model a spectral line at rest wavelength
is split into two Zeeman components with wavelength
and
,
where
is the typical wavelength shift corresponding to the local line-of-sight component of the magnetic field (Mathys 1989). To determine the final profile we integrate over the visible stellar disc, using a limb darkening constant
(Gray 1992),
km s-1, and an intrinsic line width of 10 km s-1. This high value for
is required to reproduce the average line profile, which is broadened by pulsations and the averaging process. We checked this model by reproducing Stokes V profiles for
Cep which has a polar field strength of 360 G (Donati et al. 2001; Henrichs et al. 2000).
From this model we find that for a given field strength and
the amplitude of the Stokes V profile mainly depends on the angle between the rotation axis projected onto the plane of the sky and the magnetic axis (the angle
in Fig. 5). Although the shape of the profile and
depend on whether the magnetic axis is pointing towards or away from us (maximum and minimum
)
or lies in the plane of the sky (
), the amplitude of the profile is rather independent of this. Example profiles and the dependence of the maximum amplitude on
are shown in Fig. 6. The maximum amplitude approximately scales with
.
![]() |
Figure 7:
Upper limit to the magnetic field strength at the magnetic pole as a function of the angle ![]() |
Open with DEXTER |
For our observations (except for Nr. 3, which has a very low S/N), magnetic polarisation signatures in Stokes V with an amplitude larger than approximately 0.04% would have been detected (see Fig. 4). With our model it is possible to constrain the strength of a dipolar magnetic field for a given .
For a maximum amplitude of 0.04%, we find that the upper limit on the polar field strength is:
(see Fig. 7). To hide the predicted field of
kG,
would have to be smaller than about 3.5
.
So the angle between the magnetic field axis and the projected rotation axis would have to be smaller than this 3.5
,
for all observations.
Generally (unless the rotation axis lies exactly in the plane of the sky)
depends on the rotational phase. Our observations of February 2004 cover a period of 10 days, which, with an estimated period of
Eri of 1-2 months, corresponds to 1/3 to 1/6 of the full rotation period. Since there are two magnetic extrema every rotation period, the inclination of the rotation axis would have to be smaller than
10
to allow
to be
over this whole rotation phase, with both
(the angle between the rotation and the magnetic axis) and
(the rotational phase) fine-tuned to minimise
.
It seems very improbable to have all these parameters conspire to hide a magnetic signature.
In the absence of a magnetic field, we can conclude that the observed UV variability is due to the strong pulsations in this star, which is supported by the short timescale of the variability. However, the asymmetry of the splitting of the pulsation triplet around 5.26 c d-1 remains unexplained. In view of the discovery of two more triplets with different splittings and asymmetries, more sophisticated modeling of this star and all three triplets is required before further conclusions can be drawn on the relation between the stellar rotation and the splitting, and asymmetry, of the triplets.
Acknowledgements
E.V. and H.F.H. would like to thank W. Dziembowski for inspiring discussions onEri during the Mmabatho meeting on magnetic fields in South Africa, November 2002 when this project was initiated. Most of the TBL observations were taken in service mode. Without this efficient observing mode this project could not have been done. In particular we acknowledge M. Aurière and F. Paletou for observing. We are also indebted to the capable TBL staff for assisting with the observations, G. Handler for providing radial velocity information, and M. Smith, the referee, for his constructive comments. This research has been partly based on INES data from the IUE satellite and made use of the Simbad and ADS databases maintained at CDS, Strasbourg, France.