A&A 437, 227-234 (2005)
DOI: 10.1051/0004-6361:20042554
S. J. O'Toole 1 - S. Jordan 2 - S. Friedrich 3 - U. Heber 1
1 - Dr. Remeis-Sternwarte, Astronomisches Institut der
Universität Erlangen-Nürnberg, Sternwartstr. 7, Bamberg 96049,
Germany
2 - Astronomisches Rechen-Institut, Mönchhofstr. 12-14, 69221
Heidelberg, Germany
3 - Max-Planck-Institut für extraterrestrische Physik, Giessenbachstr.,
85748 Garching, Germany
Received 16 December 2004 / Accepted 4 March 2005
Abstract
We present polarisation measurements of sdB and sdO
stars using FORS1 on the VLT. The observations were made as
part of a project to determine whether magnetic fields in two
super-metal-rich stars can explain their extreme abundance
peculiarities. Field strengths of up to
1.5 kG range have been
measured at varying levels of significance in each of our six targets,
however no clear evidence was found between apparently normal subdwarfs and
the metal-rich objects. The origin of the magnetic fields is unknown. We
also discuss the implications of our measurements for magnetic flux
conservation in late stages of stellar evolution.
Key words: stars: subdwarfs - stars: magnetic fields
Hot subdwarfs are subluminous objects that dominate the population of faint
blue stars in our own galaxy. The subdwarf B stars (sdBs) have hydrogen-rich
atmospheres with effective temperatures below about 40 000 K
(e.g. Heber 1986). They are also typically helium-poor. Subdwarf O (sdO)
stars on the other hand cover a much larger range of atmospheric compositions
with a large spread of hydrogen and helium abundances. Their effective
temperatures are in the range between 40 000 and 90 000 K. In the
diagram they are found close to ("born-again'') post-AGB or
post-extreme horizontal branch tracks (e.g. Husfeld et al. 1989; Ströer et al. 2005). Both
groups are believed to be the progeny of 1-2% of the white dwarfs.
Recently, two subdwarf B stars (UVO 0512-08 and PG 0909+276) were
discovered to have unusually high abundances of iron-group elements, e.g.
these stars have up to 105 times more vanadium and nickel than the Sun
(Edelmann et al. 2001). While no direct evidence was found in their spectra, one
obvious possibility is that these anomalously high abundances are somehow
related to large magnetic fields, as seen in the chemically peculiar A
stars. Inspired by this discovery, we started a project to measure magnetic
fields in hot subdwarfs.
Table 1: List of targets observed. Super-metal-rich objects are marked with an asterisk.
There is only one published study of a search for magnetic fields in sdB and
sdO stars, which was carried out using the Special Astronomical Observatory
6 m telescope (Elkin 1996).
The author inferred a longitudinal magnetic field of
G (based
on three measurements of the polarisation in the He I 5875 Å line of the sdO BD +75
325 (V=9.21) and a variable field in the
sdOB Feige 66 (V=10.64) with measurements of
G
(He I 5875),
(He I 5875), and
G (H
,
He I 5875), and +460 G (H
), each
value based on one or two measurements. The uncertainties were (probably)
calculated from the deviations of successive observations without taking into
account individual errors in measurements or sytematic errors.
In this paper we present spectropolarimetric observations of sdBs and sdOs based on polarisation measurements using the Balmer lines and He I and II lines.
Circular polarisation spectra of six hot subdwarfs were obtained using the VLT/UT1 (Antu) with the FORS1 spectrograph on February 23, 2004 (UT) in visitor mode. FORS1 is a multi-mode focal reducer imager and grism spectrograph equipped with a Wollaston prism and rotatable retarder plate mosaics in the parallel beam allowing linear and circular polarimetry and spectropolarimetry. We used the 600B grism, covering the spectral range 3400-5900 Å, with a 0.5'' slit, leading to a spectral resolution of 2.8 Å. Details about our targets can be found in Table 1. Two out of six of them are the super-metal-rich sdB stars found by Edelmann et al. (2001), one is an apparently normal sdB, one is a post-EHB sdB, while the last two stars are sdOs. The two metal-rich objects were observed in the first half of the night, and subsequent bad weather meant that more "normal'' abundance targets could not be observed.
Bias frames, flat-fields and He+HgCd arc spectra were taken at the end of the night. The data were reduced using standard IRAF routines for bias subtraction, cosmic ray removal and flat fielding. After the flat-field correction, the stellar spectra were extracted from each frame by summing up all CCD rows for the e and obeams. Background sky light was averaged over 10 rows on either side of the object spectrum (giving a total of 20 rows) and subtracted.
As noted by Jordan et al. (2005),
it is very important to make sure that the spectra
are correctly wavelength calibrated; an incorrect calibration might lead to
a spurious polarisation signal.
Calibration was done independently for the spectra of each beam (e and o)
and each position of the retarder plate (
). The wavelength
accuracy of each spectrum is about 3 km s-1 or 0.05 Å at H
,
which is much lower than the spectral resolution. Stokes I spectra
(i.e. unpolarised, the sum of the e and o beams) are shown in
Fig. 1 with the continuum normalised to unity and each spectrum
offset by 0.5 (not scaled) for clarity. Note that the top two spectra are the
helium-rich sdO stars, the two middle spectra are the super-metal-rich sdBs
and those at the bottom are the "normal'' abundance sdBs.
![]() |
Figure 1: Normalised spectra of our six targets. The spectra are offset by 0.5 for clarity. |
| Open with DEXTER | |
To derive the level of circular polarisation from such a sequence, we added
the exposures taken with the same quarter-wave plate angle. Stokes V/I can
then be found using
Table 2: Effective Landé factors for the observed He I and He II lines, along with the next three strongest lines that may be observed with a different FORS1 grism.
The mean longitudinal field can be derived using the weak field approximation,
which is valid in the presence of instrumental broadening, but not in the case
of rotational broadening. The limiting value of
is dependent
on the intrinsic line width; the Balmer lines are not expected to be
affected because Stark broadening is dominant. For example, if the line
width is
15 Å (typical for the Balmer lines in hot subdwarfs),
would have to be greater than 1000 km s-1, while a narrow
line with width 3 Å would be useless if
km s-1.
None of our targets are known to be rapid rotators in any case. The
field strength is compared with the Stokes V/I spectrum using the following
equation (Landstreet 1982):
The longitudinal component of the magnetic field for each measurement was
determined by comparing the observed circular polarisation for an interval of
Å around the strongest absorption lines with the prediction of
Eq. (2). As in Aznar Cuadrado et al. (2004) and Jordan et al. (2005),
we determined
using a
-minimisation procedure. Unlike the case
of central stars of planetary nebulae, systematic errors due to the blending
of Balmer lines and the Pickering series He II lines should be
small. This is because in the stars where the hydrogen lines are strong, the
helium lines are not strong, and vice versa. In the case of strong Balmer
lines, some of the helium lines are very weak. The strongest are
He I 5875, He I 4471, He I 4921 and
He II 4686. We used the criterion that the line must be stronger
than 10% below the continuum to be measurable. Because of this, only in the
two He-rich sdO stars could other helium lines be used.
Other systematic errors besides the blending effect are difficult to estimate
in strength, but we believe that they are lower than 200 G.
Table 3: Mean longitudinal field strengths (in Gauss) measured from all lines that satisfied our 10% below the continuum criterion. No entry means the line is either not present, or does not satisfy the criterion.
Previous spectropolarimetric studies using Balmer or He II lines
have not needed to consider the effect of the Landé g-factor on the
magnetic field strength. This is because for these lines
is approximately unity. This is not the case, however, for
many of the He I lines that we can measure in our
subdwarfs. We have determined the effective Landé factors with the classical
formula (e.g. Leone et al. 2000).
Table 4: For Bz=0, 600, and 1200 G we produced 1000 simulated polarisation spectra with the same noise level as the programme stars and used the same spectral lines as in the fits of Table 3. We list the weighted means and standard deviations (from all spectral lines) of the fits at each input Bz and additionally the smallest and largest value from all fits.
The results of our analysis are shown in Table 3 with results
for each line measured, along with the weighted average field
strength. The values shown here must be intrinsic to the star (and not
e.g. interstellar or instrumental), since no polarisation effects are
visible in the continuum. The line errors are simply a product of the
test (see Aznar Cuadrado et al. 2004; and Jordan et al. 2005,
for details), while the
error in the average field strength was calculated using
,
where the
are the individual line
errors. To the untrained eye, it appears that the results for
individual lines are internally inconsistent. This will be discussed below,
as will the reality of our errors, and the statistical significance of
field strengths we have measured.
As was already demonstrated by Aznar Cuadrado et al. (2004) for white dwarfs and Jordan et al. (2005) for central stars of planetary nebulae, even a small polarisation signature, slightly lower than the noise level of the spectro-polarimetric data, can be detected when several spectral lines (and multiple observations) are used. We again investigated this with simulated data of the same noise level than in the subdwarf observations.
For this purpose we produced simulated polarisation spectra for
assumed magnetic fields between 0 G and +1200 G and steps of
200 G. To these polarisation
spectra we added Gaussian noise with the same standard deviation
as in the observations of the respective subdwarfs.
We assume Gaussian noise since we expect the CCDs
to be linear and that photon noise is the dominating limitation. The
difference between Poisson and Gaussian noise is negligible at such high S/N,
so we believe that our approximation is justified on theoretical
grounds. We have also verified this by examining the frequency distribution
of the noise.
For each of the assumed magnetic fields we made 1000 simulations and
calculated the magnetic fields of the fits from the same spectral lines we
have used for the observations of each star. In Table 4 we
list the results for 0 G, +600 G, and +1200 G and provided the mean
result, the standard deviation, and the smallest and largest value. This
should give us an impression of how realistic our statistical errors from the
analyses are. In this table we present only the mean field
strength derived from our simulations, not the value for individual
lines. It is important to note, however, that the apparently internally
inconsistent line-by-line results are also seen in our simulations. This
occurs because
in Eq. (2) is nonlinear,
meaning the noise deviations do not translate linearly into V/I deviations. Put simply, the weighted average field strengths are more
well defined.
Firstly, the mean results of the respective 1000 simulations are very
close to the predescribed magnetic fields with deviations up to 100 G
for an assumed value of 1200 G. Secondly, the standard deviations are
close to the
errors determined from the
fits.
Moreover, the smallest and largest values from all 1000 fits for
an assumed field of 0 G provides
a very conservative estimate of what magnetic field can be mimicked by pure
noise: in all cases the fitted field strengths are larger in their
absolute value than these extreme numbers; only in the case of
CPD
the derived magnetic field strength of
Bz=-885 G comes close to the simulated extreme deviation from zero
of -862 G. Therefore, we conclude that all of our results
are incompatible with Bz=0, with a somewhat smaller significance for
CPD
.
Two of the objects (UVO 0512-08 and HD 76431) are
even incompatible with Bz=600 G. This indicates that our
errors
are very realistic and that all stars very probably have a magnetic
field. Note also that our error ranges are much smaller than for the stars
studied by both Aznar Cuadrado et al. (2004) and Jordan et al. (2005),
in part because we often have
more spectral lines available to measure and because some of the objects
are much brighter. The individual objects are discussed in more detail
below, and their Stokes V/I spectra are shown in Figs. 2-7.
![]() |
Figure 2:
Observed Stokes V/I spectra for PG 0909+276 for various
hydrogen and helium lines, with value predicted by Eq. (2)
overlaid (thick line). Note that the He I 4713 Å line (second
panel from the top) has a Landé factor |
| Open with DEXTER | |
![]() |
Figure 3: Similar to Fig. 2, except for UVO 0512-08. |
| Open with DEXTER | |
![]() |
Figure 4:
Similar to Fig. 2, except for CPD
|
| Open with DEXTER | |
![]() |
Figure 5: Similar to Fig. 2, except for HD 76431. In this star the polarisation signature is evident in every line. |
| Open with DEXTER | |
![]() |
Figure 6:
Similar to Fig. 2, except for CD
|
| Open with DEXTER | |
![]() |
Figure 7: Similar to Fig. 2, except for LSE 153. |
| Open with DEXTER | |
These stars are the super-metal-rich sdOB objects discovered by
Edelmann et al. (2001), so we are particularly interested in their magnetic field
strengths. As can be seen from Table 3, we measure
for UVO 0512-08 and
for PG 0909+276. Both of these results are
statistically significant (see the previous section), and are also the two
highest field strengths we have measured. Whether or not this is physically
significant remains to be seen, since the two "normal'' sdBs we have observed
may have evolved differently (as a binary system), or have evolved well away
from the EHB. The super-metal-rich stars are both apparently single, although
with its low surface gravity and high temperature, UVO 0512-08 appears to
have also evolved away from the EHB.
CPD
is an sdB in a short-period (0.27 d) binary with an
unseen companion of low mass, possibly a brown dwarf (Edelmann et al., in
preparation). This star has the lowest measured field and the measurement with
the lowest significance, with
,
however,
the probability for an almost kilogauss field is still very
high. CPD
had the smallest number of spectral lines available
to measure polarisation. We are uncertain how binary evolution - most likely
the system has passed through a common envelope phase - might affect the
magnetic field of the star's progenitor.
HD 76431 was originally classified as a main sequence star at high
galactic latitude. Its subluminous nature was discovered by a detailed
quantitative spectral analysis by Ramspeck et al. (2001). The helium deficiency
and a peculiar metal abundance is similar to that of the sdB stars, although
the deviation from solar values are smaller than typical for these
stars. The gravity of HD 76431 is slightly higher than that of a main
sequence star, but considerably lower than for EHB stars. Comparison to
EHB evolutionary tracks suggests that it has already left the EHB and is
evolving towards the white dwarf cooling sequence as has been suggested for
some sdO stars. Therefore HD 76431 is not a typical sdB, but might be
considered as a link between the sdB and sdO stars. The magnetic field
measured (
G) has very high significance as found in
Sect. 4.1. This star has the clearest "by eye''
polarisation signature.
Both of the sdO stars we observed are helium-rich. For CD
Bauer & Husfeld (1995) found
K and
and
measured abundances; the star contains
30 time more helium then
hydrogen, and appears to be a "garden variety'' sdO (see,
e.g. Ströer et al. 2005). We measure
based on the 11
lines shown in Table 3, and find from our simulations that
this value is statistically significant.
LSE 153 is an extremely hydrogen-deficient sdO that has been suggested as a
possible successor of a R CrB star, with a mass of
0.7
(Husfeld et al. 1989). The star's very low hydrogen abundance is accompanied by both
nitrogen and carbon enhancements. The field strength we measure
(
G) was also found to be statistically significant. This
star is at a similar evolutionary stage as the central stars of planetary
nebulae studied by Jordan et al. (2005),
and has a field strength of approximately
the same order of magnitude.
In both of the helium-rich sdOs we have observed the magnetic field determined using the He II 5412 line does not match the other He II lines. We cannot see any obvious reason for this strange behaviour.
The goal of our project was to determine whether magnetic field strength is
directly connected to abundance anomalies in these stars. The magnetic field
could have an influence on the chemical composition by suppressing convection
or by influencing mass-loss or diffusion. While convection is not expected to
be important in helium-poor sdB stars, He II ionisation can initiate
convection in helium-rich sdO stars at
between 40 000 and
50 000 K (Groth et al. 1985). This may be relevant for CD
,
but
not for the other programme stars. However, no clear differences exist between
the field strengths of the "normal" and super metal-rich stars in our sample.
The origin of the fields detected here are unknown, and are found in
objects from various evolutionary channels. It is possible that they are
dynamo induced or a remnant from their progenitor. We consider a
dynamo-induced magnetic field may be possible if rapid internal rotation
occurs as suggested recently by Kawaler & Hostler (2005), who found evidence that
several of the sdB pulsators may have rapidly rotating cores. This is
contrasted with the stars' slow surface rotation (typically
km s-1, e.g. Heber et al. 2000). While the sdBs have probably followed a
similar evolutionary path, the same cannot be said for at least one of the sdO
stars. LSE 153 can be associated with the "born-again" post-AGB scenario
(Husfeld et al. 1989), while the helium-rich sdO, CD
,
is similar to
most other objects of its class.
It is interesting to estimate the strength of the magnetic fields when these
stars have evolved to white dwarfs. If we assume complete conservation of
magnetic flux the magnetic fields are amplified inversely proportional to the
square of the stellar radius. The radii of the sdB stars are
0.15-0.25
(taking the canonical mass of an EHB star to be
0.48
). Since they will evolve directly into white dwarfs, their
radii will shrink by a factor of
20, leading to field strengths of up to
500 kG, values which are apparently rarely seen. The search for rotation in
white dwarfs (Heber et al. 1997; Koester et al. 1998; Karl et al. 2005) from high resolution H
spectroscopy also resulted in constraints on field strengths. Amongst the
50 DA white dwarfs studied, only four turned out to be magnetic with Bup to 180 kG, while upper limits of 10-20 kG could be derived for the rest.
Aznar Cuadrado et al. (2004) have detected magnetic fields between 2000 and 4000 G in four
white dwarfs, while no magnetic fields on this level of accuracy were found in
the other eight programme stars.
However, since subdwarfs contribute to only a small percentage (1-2%) of
white dwarfs progenitors, and the current sample of sub-MG white dwarfs is
still small, we cannot completely rule out the existence of low-mass
(
0.5
)
magnetic white dwarfs, although the population must not be
large. It is much more difficult to work backwards to values one might
expect for the main sequence progenitors of our targets, since sdB formation
is not well defined. However, if we take a star with the same parameters as
the Sun as an example, we find field strengths in the 20-40 G range, which
seem reasonable (although currently undetectable). This also appears to rule
out any evolutionary connection with the Ap stars.
Jordan et al. (2005) have found kilogauss magnetic fields in all of the central stars of planetary nebula that they have observed by means of spectro-polarimetry. By magnetic field amplification these stars will have magnetic fields strengths between 0.35 and 2 MG if the assumption of complete conservation of magnetic flux is true. The same holds for the helium-rich sdO star LSE 153 analysed in this paper. Since about two thirds of the white dwarfs seem to have magnetic field even lower than 1 kG (Aznar Cuadrado et al. 2004) one must conclude that magnetic fields can be destroyed during the final stages of stellar evolution, although the detailed mechanisms are unknown.
There are several questions that have arisen from our detection of magnetic fields in hot subdwarfs. Firstly, are the field strengths dependent on the binary status of the star? The one star we have observed that is known to be in a binary system has the lowest and a somewhat less statistically significant measurement. Since it is only one object we cannot make any conclusions at this stage, but this may warrant further investigation. Secondly, what effect, if any, might kilogauss magnetic fields have on pulsations in sdB stars? Rotation is often considered in connection to sdB pulsators, however magnetic fields have not as yet. Asteroseismology has been used to infer a magnetic field in at least one white dwarf (the pulsating DB GD 358, see Winget et al. 1994). This field may be dynamo induced (Markiel et al. 1994). Spectropolarimetric observations of sdB pulsators, along with a theoretical investigation, would help to clarify this situation.
In this study, we have used polarisation measurements of sdB and sdO stars to
try to determine whether magnetic fields in two super-metal-rich stars can
explain their extreme abundance peculiarities. Field strengths of up to
1.5 kG range have been measured at varying levels of significance in
each of our six targets, however no clear difference was found between
apparently normal subdwarfs and the metal-rich objects. The origin of the
magnetic fields is unknown. We also considered the implications of our
measurements for magnetic flux conservation in late stages of evolution.
Unless there is a population of low-mass white dwarfs with field strengths up
to
500 kG, it is unlikely that magnetic flux is conserved in late
stages of evolution. This idea is consistent with findings for central stars
of planetary nebulae.
Acknowledgements
We would like to thank E. degl'Innocenti for his advice on Landé factors, and Svetlana Hubrig for help with the observational setup. S.J.O.T. is supported by the Deutsches Zentrum für Luft- und Raumfahrt (DLR) through grant No. 50 OR 0202.