A&A 423, 1081-1094 (2004)
DOI: 10.1051/0004-6361:20040355
R. Aznar Cuadrado1 - S. Jordan2,3 - R. Napiwotzki4 - H. M. Schmid5 - S. K. Solanki1 - G. Mathys6
1 - Max-Planck-Institut für Aeronomie, Max-Planck-Str. 2,
37191 Katlenburg-Lindau, Germany
2 -
Institut für Astronomie und Astrophysik, Eberhard-Karls-Universität
Tübingen, Sand 1, 72076 Tübingen, Germany
3 -
Astronomisches Rechen-Institut, Mönchhofstr. 12-14, 69120 Heidelberg,
Germany
4 -
Department of Physics & Astronomy, University of Leicester, University
Road, Leicester LE1 7RH, UK
5 -
Institut für Astronomie, ETH Zentrum, 8092 Zürich, Switzerland
6 -
European Southern Observatory, Casilla 19001, Santiago 19, Chile
Received 27 February 2004 / Accepted 12 May 2004
Abstract
We have detected longitudinal magnetic fields between 2 and 4 kG in three
(WD 0446-790, WD 1105-048, WD 2359-434) out
of a sample of 12 normal DA white dwarfs by using optical spectropolarimetry
done with the VLT Antu 8 m telescope equipped with FORS1. With the exception
of 40 Eri B (4 kG) these are the first positive detections of magnetic fields
in white dwarfs below 30 kG. Although suspected, it was not clear whether a
significant fraction of white dwarfs contain magnetic fields at this level.
These fields may be explained as relics from magnetic fields in the
main-sequence progenitors considerably enhanced by magnetic flux conservation
during the shrinkage of the core. A detection rate of 25% (3/12) may
indicate now for the first time that a substantial fraction of white dwarfs
have a weak magnetic field. This result, if confirmed by future observations,
would form a cornerstone for our understanding of the evolution of stellar
magnetic fields.
Key words: stars: magnetic fields - stars: individual: WD0446-790, WD1105-048, WD2359-434
According to the field amplification theory, the white dwarfs play an important role in the investigation of stellar magnetic fields. In main-sequence stars magnetic fields have been detected directly mainly for peculiar magnetic Ap and Bp stars with rather well organized fields and field strengths of the order 102-104 Gauss. For weak fields in A to O stars (B<102 G) direct magnetic field detections are still very rare (e.g. the field detections in early B stars reported by Neiner et al. 2003b,a). For sun-like stars ample evidence (coronal activity) for the presence of complicated small-scale fields exists, but direct measurements are only possible for the more active stars (Valenti & Johns-Krull 2001; Rüedi et al. 1997; Saar 1996). The contraction to a white dwarf amplifies the magnetic fields by about 4 orders of magnitude, so that weak and often undetectable magnetic fields on the main sequence become measurable during the white dwarf phase. This is supported by the known magnetic white dwarfs with megagauss fields ( B=106-109 G). Their frequency and space distribution, as well as their mass, are consistent with the widely accepted view that they are the descendents of the magnetic Ap and Bp stars (e.g. Mathys 2001). Another origin seems to be required for the magnetic degenerates with weaker fields (unless magnetic flux is lost during the contraction phase). Magnetic main-sequence stars with weaker magnetic fields have been suggested as their possible progenitor candidates (Kawka et al. 2003; Schmidt et al. 2003). The B stars on which weaker fields have been detected may be the missing stars. However, even the most sensitive observations are limited to some tens of gauss on main-sequence stars.
Thus, magnetic field amplification during stellar evolution may offer the
opportunity to investigate 1 G magnetic fields (averaged global fields)
in normal main-sequence stars with observations of
1 kG magnetic fields
during the white dwarf stage. White dwarfs with magnetic fields below 100 kG
have been either found by searching for circular polarization
(Schmidt & Smith 1994) or by looking for Zeeman splitting of narrow NLTE
line cores in the Balmer lines, particularly in H
(Koester et al. 1998). However, the splitting becomes undetectable in
intensity spectra for weak fields (<20 kG) or for objects without narrow
line core. Therefore, spectropolarimetry is the most promising technique for
successful detections of weak magnetic fields. Up to now detections of
magnetic fields below 30 kG have not been achieved, except for the very bright
white dwarf 40 Eri B (V=8.5), in which Fabrika et al. (2003) have detected
a magnetic field as low as 4 kG. The magnetic field detection limit can now
be pushed down to a few kG for many white dwarfs with spectropolarimetry using
8-10 m class telescopes.
In this paper we present and analyse VLT spectropolarimetric data of a sample
of 12 white dwarfs in a search for weak magnetic fields. In Sect. 2 the
observations and data reduction are described, while in Sect. 3 the
observational method for obtaining the Stokes parameter (V/I) is described.
Section 4 presents the method for determining weak magnetic fields analysing the
circular polarisation due to a given magnetic field. In that section we also
present the results of our analysis, along with the description of the
-minimisation procedure applied to our data. The determination of the
atmospheric and stellar parameters is presented in Sect. 5 and compared with
those found in the literature. A discussion and conclusions are presented in
Sect. 6.
With a 0.8
wide slit we obtained a (FWHM) spectral resolution of
4.5 Å. The data were recorded with a backside-illuminated thinned, AR
coated Tek.
CCD with 24
m pixels which correspond to a
pixel scale of 0.2
/pixel in spatial and 1 Å/pixel in spectral
direction.
Spectra were acquired with grism G600B, which allows observations in the
spectral range 3400-6000 Å, covering all H I Balmer lines from H
to the Balmer jump simultaneously. A reflex image from the FORS
optics affects the wavelength region from 4000 to 4100 Å which corresponds
to the blue wing of the H
line. Although the reflex shows up in the
intensity spectrum it seems to produce no spurious signal in circular
polarisation. The reflex is known to occur for this particular grism G600B / Wollaston configuration (e.g. Schmid et al. 2003).
Table 1:
Details of VLT observations. The provided
and
coordinates refer to epoch 2000 as measured in the course of the SPY project
(see Koester et al. 2001).
Spectral types,
,
and measured V magnitudes were taken from the
catalogue of McCook & Sion (1999).
The observations were split into several cycles to avoid saturation. The feasibility of circular spectropolarimetry with the required high signal-to-noise ratio using FORS1 had been demonstrated by Bagnulo et al. (2004,2002) with test measurements of one magnetic and one non-magnetic A-star. Stellar rotation may cause continuous changes in the field orientation and, therefore, in the polarisation signal. Hence, we have split the observations of an individual target into more than one polarisation measurement taken during different nights (see Sect. 3). In this way, there is a much lower probability that a candidate with a magnetic field escapes detection due to a special orientation at the time of the observations (where circular polarisation cancels). Hence, having for most objects data for more than one epoch, we are able to assess the presence of rotational modulation of a possible magnetic field.
Our sample was carefully selected on the basis of VLT-UVES spectra taken within the SPY survey (Napiwotzki et al. 2003): SPY (Supernovae type Ia Progenitor surveY) is a radial velocity search for close binary systems of two white dwarfs (double degenerates; DD). If these systems are close enough they will merge due to gravitational wave radiation and if the combined mass of these mergers exceed the Chandrasekhar limit for white dwarfs these are potential progenitors of Supernovae type Ia.
SPY was carried out with the high-resolution Echelle spectrograph UVES at the Kueyen (UT2) of VLT. With the set-up used for SPY UVES provides almost complete spectral coverage of the wavelength range from 3200 Å to 6650 Å, with a spectral resolution of 0.3 Å. For more details, please refer to Napiwotzki et al. (2003,2001). SPY observations were used to check candidates for our project for spectral peculiarities and magnetic fields strong enough to be detected in intensity spectra. Hence, our targets were selected with the criterion of not having any sign of Zeeman splitting visible in the SPY spectra, and hence no magnetic fields above a level of about 20 kG.
All our targets have strong hydrogen lines, ideal for measuring line polarisation, and no peculiarities (such as MG-magnetic fields or a bright companion or any indication of magnetic fields).
In all frames the bias level was subtracted and the frames were cleaned of cosmic ray hits. For each observing night a unique master flat-field was calculated from the median flux of all flat-fields taken that night. The stellar spectra were extracted from the flat-field corrected frames as a sum over about forty CCD rows for each beam. Sky spectra were obtained from adjacent regions (about twenty CCD rows) on the detector to the observed stars (below and above) and subtracted from the object spectra. Wavelength calibration was done using HgCd, He and Ar arc spectra, which was independently performed for wavelengths of the ordinary and extra-ordinary beams with the corresponding beams of a reference spectrum.
Stokes I was obtained as a sum over all beams, while the calculation of Stokes (V/I) is described in the following section.
In the ideal case, the polarisation information (one Stokes parameter per
exposure) is contained in the ratio, at each wavelength, of the intensities
in the two spectra (from the ordinary and extra-ordinary beams) but it is
mixed up with the system gain ratio for the pixels concerned. Alternatively,
the effect of the not well defined pixel gain can be significantly reduced
by inverting the sign of the polarisation effects in a second exposure (by
rotating the -plate by
), while leaving the gain ratios
identical (Tinbergen & Rutten 1997). Hence, Stokes V is obtained from a
differential measurement of photon counts in either the ordinary or
extra-ordinary beams, measured at two different angles of the retarder
waveplate. In this way errors from changes in the sky transparency,
atmospheric scintillation, and various instrumental effects are significantly
reduced, so that photon noise remains as the dominant error source.
Another important source of sytematic error when obtaining Stokes V could be
the wavelength calibration procedure. An incorrect calibrated wavelength scale
for each of the two analysed spectra will yield the line profiles of the two
spectra not perfectly aligned, even in the absence of a magnetic field,
leading to spurious polarization signals in each line, that change sign as
the wave plate rotates. No such spurious signals were detected during the
calibration process, confirming that FORS is a very stable instrument. Even
if spurious line signals at the noise level are present, they would be
compensated for at least partly by the combination of data taken with retarder
plate position angles
and
.
We note that wavelength
calibrations were made for each observing date separately. Because most
targets were observed for two different dates we could check that no spurious
magnetic field detection due to wavelength calibration errors are present.
We adopted the FORS1 standard observing sequence for circular polarimetry
consisting of exposures with retarder plate position angles
and
.
The number of exposures, n, is reported in Table 1 for each
object. For instance, for n=4 the sequence of position angles would go:
,
,
,
.
In order to derive the circular polarisation from a sequence of exposures, we
added up the exposures with the same quarter-wave plate position angle.
The Stokes (V/I) can be obtained as
The 12 selected white dwarfs are listed in Table 1 together with some characteristics of the observations. The Heliocentric Julian Dates correspond to the beginning of each observing sequence.
With FORS1 on the VLT we reach a noise level in circular spectropolarimetry
of about
,
where
is the continuum
intensity. This allows us to detect magnetic fields of a few kG for the
brightest white dwarfs in our sample with exposure times of about 1 h.
This is a significant improvement compared to the observations of
Schmidt & Smith (1995), who reached a detection limit of about 20 kG
(
limit) with a 2 m class telescope. By adding the signal of
several Balmer lines, which are observed simultaneously with our instrument
setup, the magnetic sensitivity can be further enhanced. However, the error
of the magnetic field determination increases for the higher series numbers,
so that a reliable analysis can be based on H
and H
only.
The theory of spectral line formation in a magnetic atmosphere shows that the
splitting of a spectral line observed in both senses of circular polarisation
is proportional to
,
the average of the component of the
magnetic field along the line of sight averaged over the visible stellar
hemisphere, i.e., the mean longitudinal magnetic field
(Babcock 1947,1958).
For field strengths below 10 kG the Zeeman splitting of the Balmer lines is
less than approximately 0.1 Å. This is well below the width of the cores
of the Balmer lines in all the stars of our sample (typically a few Å).
Therefore, we can apply the weak-field approximation
(e.g., Angel & Landstreet 1970; Landi degl'Innocenti & Landi degl'Innocenti 1973) without any loss of
accuracy. According to this approximation the measured V and I profiles
are related to
by the expression:
Equation (2) is evaluated by using the high signal-to-noise spectra
and their derivative
.
For a given
wavelength
we approximate
by the
average of
and
.
The error associated with the determination of the longitudinal field obtained
from individual Balmer lines is larger for Balmer lines at shorter wavelengths
than for lines at longer wavelengths. This is due to the combination of two
effects: while the Zeeman effect increases as lambda squared, most other
line
broadening effects depend linearly on lambda, so that the magnetic field is
better detected at longer wavelengths than at shorter wavelengths;
furthermore, the Balmer lines at shorter wavelengths are less deep, so that
is smaller. Using H
and H
simultaneously, we obtained a determination of the mean longitudinal magnetic
field that best fit the observed (V/I) (as explained in the next section).
![]() |
Figure 1: Normalised spectra of our sample of white dwarfs in the region of the Balmer series. When more than one sequence of exposures was available for a white dwarf, the average spectrum is presented. The positions of the Balmer lines are indicated. The spectra are displaced vertically by 0.7 units (or multiples thereof) for a better visualisation. |
Open with DEXTER |
For the circular polarisation spectra (V/I) flux weighted means were
calculated according to the following formula for two measurements 1
and 2:
Figure 2 shows the circular polarisation spectra (V/I) of our sample
of white dwarfs in the region of the Balmer series (no binning along the
spectrum was applied to
the data). The positions of the Balmer lines are indicated by vertical dashed
lines. In two of the averaged spectra (corresponding to stars
WD 0446-789 and WD 2359-434) a moderate S-shape circular
polarisation signature across H
and H
can be noticed. For
WD 1105-048 polarisation reversals of those lines were only present
in the observation of 9 Jan. 2003 (see Fig. 5).
![]() |
Figure 2: Circular polarisation spectra (V/I) of our sample of white dwarfs in the region of the Balmer series. The average (V/I)-spectrum is plotted for multiple observations. For all spectra a horizontal line is drawn to indicate the zero level. Dashed vertical lines represent the positions where the Balmer lines are located. Spectra are displaced with relative shifts for a better visualisation. |
Open with DEXTER |
To determine the longitudinal component of the magnetic field for each
measurement we compared the observed circular polarisation for an interval
of 20 Å around H
and H
with the prediction of
Eq. (2). The best fit for
,
the only free
parameter, was found by a
-minimisation procedure. If we assume that
no magnetic field is present, all deviations from zero polarization are due
to noise. This can be expressed by the standard deviation
over the
respective intervals around the Balmer lines. If the best fit magnetic field
is indeed close to zero, the reduced
should be automatically close to
unity. However, in the case of a finite magnetic field not all deviations from
zero would be due to noise, so that the best fit would have a reduced
smaller than 1. Following Press et al. (1986) we determined the statistical error
from the rms deviation of the observed circular polarisation from the best-fit
model. The
(68.3%) confidence range for a degree of freedom of 1 is
the interval of Bz where the deviation from the minimum is
.
Note that this is a purely statistical error and does not
account for systematic errors, e.g. from the limitation of our low-field
approximation which, however, we expect to be rather small, or from wavelength
calibration errors.
The results for each of our single observations is summarized in
Table 2, where we have listed the best fits for Bz for the
H
and H
lines as well as the weighted value
where
(
). The probable error is given by
.
When multiple observations were added up to
reduce the noise level (by means of Eq. (3)), results are labeled
as AVERAGE.
Table 2:
Magnetic fields derived from the H
and H
lines for our
sample of white dwarfs.
provides the magnetic field in units of
the
level. Detections exceeding the
levels are given in
bold.
is the standard deviation of the observed (V/I)-spectrum
obtained inthe region 4500-4700 Å. Lower limits on the detectability of
the magnetic field from the line polarisation peaks, calculated at the
level of the noise, are given in the last two columns. Multiple
observations that were averaged prior to analysis are labeled AVERAGE.
In three stars we find a significant magnetic field: WD 0446-789
(
G, see Fig. 3), WD 2359-434
(
G, see Fig. 4), and WD 1105-048 (
G). Positive and negative signs of Bz indicate opposite
magnetic polarity. For the first two stars the magnetic field is detected at
the
level individually from H
and H
(at least for
the averaged spectra), as well as from the combination of both lines. This
increases our confidence in both detections (cf. Sect. 4.3).
Although several lines of the Balmer series are available, the higher members
of the series do not contain enough V-signal to give reliable results.
However, the analysis of H
(the only line among the higher members
with its V-signal exceeding the
), confirms the positive
detections of magnetic fields in the three objects WD 0446-789,
WD 1105-048 and WD 2359-434.
In WD 0446-789 the analysis of the 30 Nov. 2002 observation around
H
yields a magnetic field of
G, corresponding to a
(instead of
for
the 28 Jan. 2003 observation), i.e. a
level detection. If we do
not use this outlier for our combined result from the two observations, we
obtain a
detection of
G.
In the case of WD 1105-048 the different values of the magnetic
field obtained for the two observations on 1 Jan. 2003 (
G) and
29 Jan. 2003 (
G, see Fig. 5) differ
significantly. This may indicate a different orientation of the magnetic
field between the two epochs due to stellar rotation.
The case of WD 0135-052, where the formal value for the combined
result from H
and H
exceeds the
error by a factor of
two, is rather uncertain since it is based on one observation only. In
addition, this object is a double-lined binary with a period of 1.556 days
and it is not clear whether the orbital variation can mimic a magnetic field.
This could just be a
level detection, which can be expected for our
sample size of field measurements by chance. However, WD 0135-052 is a good candidate for follow-up observations, although one must bear in mind
that the magnetic field would be diluted by the presence of a non-magnetic
companion.
In all other stars in our sample the best-fit magnetic field is below or close
to the error.
![]() |
Figure 3:
Circular polarisation spectra (V/I) of WD 0446-789
(average of observations from 30/11/03 and 28/01/03, thin solid line) in the
region of H![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
As in Fig. 3 but for WD 2359-434 (average of
observations from 04/11/03 and 29/11/03), where the best fit is at
![]() |
Open with DEXTER |
![]() |
Figure 5:
As in Fig. 3 but for the 09/01/03 observation of the
white dwarf WD 1105-048 where the best fit is at
![]() |
Open with DEXTER |
As an example we used the solution for the average of both
WD 0446-789 observations (
G), its noise level
(
at both H
and H
)
and calculated 1000 simulated
polarisation spectra for the interval size of
20 Å. The average
result for the magnetic field was:
G (solution for
H
), 4276
554 G (for H
only), and
G
(for H
only). The mean value of B agrees very well with the prescribed
value. The mean of the individual standard deviations of the 1000 runs are
413 G, 541 G and 667 G, for H
,
H
and H
,
respectively. This is slightly smaller than 640 G, 832 G and 1005 G from
the WD 0446-789 observation. This may be due to the fact that we
assume a homogeneous magnetic field which may not be the case in reality, or
it may indicate a small hidden source of errors in the data (e.g. a
non-Gaussian distribution of the noise).
Figure 6 shows one of the artificial polarisation measurements from
our 1000 simulations. The fit is
G in this particular example.
It demonstrates that a highly significant detection is possible at our noise
level.
![]() |
Figure 6:
Artificial polarisation spectrum with the same noise level as
WD 0446-789 plotted in Fig. 3. For this particular
example, out of the 1000 simulations, the fitting procedure results in
![]() |
Open with DEXTER |
The resulting magnetic field changed by less than 2% if an interval of
40 Å around H
and H
was used. In the case of a
smaller interval, however, larger deviations occur since not the whole range
where the predicted (V/I) differs from zero is included.
We are aware of the fact that the polarisation feature is mostly buried in
noise in the (V/I)-spectrum, with the exception of the narrow peaks on both
sides of the hydrogen line centers. Therefore, on a plot the magnetic
information is, to a large extent, invisible to the eye since in the
wings
the solution is based on an average small excess of right- and left-handed
polarisation on different side of the line core, respectively, which clearly
contributes to our analysis.
However, if one wants to obtain an idea
of what magnetic field is needed so that the predicted polarisation peaks
exceed the
level of the noise, we provide lower limits for this
detectability in Cols. 8 and 9 of Table 2.
The
was obtained as the standard deviation of the observed
(V/I)-spectrum in a region with no lines over 200 Å (i.e. 4500-4700 Å), for each single observation.
We note that in the three objects (WD 0446-789,
WD 1105-048 and WD 2359-434) with a positive detection
according to the
analysis, the peaks of the S-wave circular
polarisation signature reaches the
level of the noise so the
detection is, therefore, confirmed by this very conservative approach. For
the averaged spectra of WD 0446-789 and WD 2359-434,
the detections derived from the more reliable H
line even reach the
level.
Since all white dwarfs selected for our project are bright by white dwarf
standards, at least one model atmosphere analysis was published for each
object. This is a somewhat inhomogeneous collection of data relying on
spectra of different quality and a variaty of methods, including analysis of
UV spectra and parallax measurements. Most of the programme stars have already
been observed with the SPY project. Although our spectra have lower resolution
than the SPY spectra, which aim at the measurement of radial velocity
variations seen in the NLTE core of H,
our flux spectra have a very
high signal-to-noise ratio, allowing a very accurate determination of the
effective temperatures and gravities.
The observed line profiles are fitted with theoretical spectra from a large grid of NLTE spectra calculated with the NLTE code developed by Werner (1986). Basic assumptions are those of static, plane-parallel atmospheres in hydrostatic and radiative equilibrium. The adopted chemical composition is pure hydrogen, which is appropriate for the stars of our sample. A description of the model calculations and the adopted atomic physics is given in Napiwotzki et al. (1999), where a discussion of the impact of NLTE on the atmospheres and line profiles of DA white dwarfs is provided as well. The coolest four white dwarfs of our sample (WD 0310-688, WD 0839-327, WD 1105-048 and WD 2359-434) were analysed with a grid of LTE model spectra computed by D. Koester for the analysis of DA white dwarfs. The input physics of these model atmospheres is described in some detail in Finley et al. (1997). LTE models are more reliable below 17 000 K, because the NLTE models do not take into account convection and collision induced absorption by hydrogen quasi-molecules. On the other side, NLTE effects are small at these temperatures (Napiwotzki et al. 1999) and of no relevance for our results.
Table 3: Fitted parameters of the white dwarfs and supplementary data from literature. The three objects with positive detections of magnetic fields are labelled in bold face.
The line fits were performed with a modified version of the least-squares
algorithm of Bergeron et al. (1992) described in Napiwotzki et al. (1999). The
observed and theoretical Balmer line profiles are normalised to a linear
continuum in a consistent manner. Radial velocity offsets are corrected by
shifting the spectra to a common wavelength scale. The synthetic spectra are
convolved to the resolution of the observed spectra (4.5 Å) with a Gaussian
and interpolated to the actual parameters. The atmospheric parameters are then
determined by minimising the
value by means of a Levenberg-Marquardt
steepest descent algorithm (Press et al. 1986).
This procedure is applied simultaneously to all Balmer lines of one
observed spectrum. Formal errors can be derived from
the covariance matrix. However, due to the very high signal-to-noise ratios of
our spectra, the estimated statistical errors are extremely small (not larger
than 20 K in temperature and 0.005 dex in the log of gravity) and do not
provide realistic estimates for the whole error budget. The real external
errors for the white dwarfs investigated in this article can be estimated to
be 2.3% in T
and 0.07 dex in
(Napiwotzki et al. 1999).
![]() |
Figure 7:
Model atmosphere fit (solid line) of the 30/11/02 observation of
WD 0446-789 (solid histogram). The dashed model profile for
H![]() |
Open with DEXTER |
Our fit results are provided in columns two and three of Table 3.
The parameters are the average of the fit results for the individual spectra.
For the reasons discussed above, we do not give the formal fit errors.
WD 0135-052 was not fitted by us, because this is a well known
double-lined binary consisting of two DA stars (Saffer et al. 1988), which
requires a special treatment. White dwarf masses were computed from a
comparison of parameters derived from the fit with the grid of white dwarf
cooling sequences of Benvenuto & Althaus (1999) for an envelope hydrogen
mass of
,
and are given in column four of
Table 3. A sample fit is shown in Fig. 7.
Most Balmer lines are fitted very well, but the observed profile of H
is only poorly reproduced.
This results from the reflex in the FORS1 optics mentioned in Sect. 2. Thus we
excluded this line from our fitting procedure. However, parameter changes are
small, if H
is included anyway. Spectroscopic distances d(spec) are
determined from the absolute magnitudes computed for the given stellar
parameters using the synthetic V band fluxes computed by
Bergeron et al. (1995) and the measured V magnitudes given in our Table 1,
and are provided in column five of Table 3. Note that the error
limits for T
and
given above correspond to a 7% error in the
distance.
Table 3 is supplemented by physical parameters of our sample of
white dwarfs collected from the literature. Fundamental parameters such as
temperature and gravity determined by previous spectroscopic studies are given
in columns six and seven. The agreement is generally good. Differences are
usually within the error limits discussed above. Rotational velocities ,
if available, are shown in column eight. If parallax measurements exist, the
resulting distances d(trig) are given in column nine (H indicates
Hipparcos values, while Y are ground based data from the Yale general
catalogue of trigonometric parallaxes).
We noticed during the fitting process that our FORS data show flat Balmer
lines cores. Available UVES spectra of this white dwarf show that the core of H
is also flat (as already reported by Koester et al. 1998).
Since an accurate parallax measurement exists for WD 2359-434, we
modified our fitting procedure for this object. The gravity was fixed at a
value, which reproduced the absolute brightness of the white dwarf computed
from the parallax. The temperature was determined from a fit of the line wings
of the Balmer lines, excluding the cores. Since the absolute brightness
depends on the stellar temperature as well, a few iterations were necessary,
before we got a self-consistent solution. The derived temperature (8.66 kK)
is consistent with the result of Koester & Allard (1993), who fitted UV
spectra taken with IUE. Kepler & Nelan (1993) using the same technique,
derived a lower temperature (7.76 kK). This is possibly caused by the low
value of the surface gravity (
)
adopted during their fitting
procedure, which is much lower than plausible values for
WD 2359-434.
The reason for the flat Balmer line cores of WD 2359-434 remains a mystery. Being a white dwarf with a broad range of magnetic fields is a possible explanation since it would smear out the sigma components. However, in our polarization measurements there is no indication for such a broad range.
With the exception of the bright white dwarf 40 Eri B, for which a magnetic
field of only 4 kG had been detected (Fabrika et al. 2003),
WD 0446-789, WD 1105-048 and WD 2359-434
have the weakest magnetic fields detected so
far in white dwarfs. In previous extended searches for weak fields in white
dwarfs (e.g. Schmidt & Smith 1995) only very few objects have been found.
Only six detections have been reported for magnetic fields below 100 kG,
which are not all confirmed, and only three objects have a field weaker than
about 50 kG. Kawka et al. (2003) have reported longitudinal magnetic fields
in three stars but no significant detection was made because their
error was almost as large as the observed value itself. They concluded that
the population of white dwarfs with magnetic fields in excess of 1 MG is well
known, but that lower-field white dwarfs remained undetected.
Note, however, that our investigation is based on the averaged longitudinal component of the magnetic field, meaning that the maximum magnetic field at the white dwarf surface can be stronger, depending on the field geometry (described e.g. by offset dipoles, or more complex distributions; with the underestimate being larger for a more complex magnetic distribution) and on the orientation relative to the observer. Therefore, our results for the three objects with a positive detection are lower limits, since cancellation effects are expected.
The population of known magnetic white dwarfs presently comprises some 125
stars (Schmidt et al. 2003; Wickramasinghe & Ferrario 2000; Gänsicke et al. 2002).
Wickramasinghe & Ferrario (2000) and Jordan (2001) established that
about 3% of all white dwarfs had a magnetic field above 100 kG. Results from
the first two years of the Sloan Digital Sky Survey
(Schmidt et al. 2003; Gänsicke et al. 2002) show the total number of
known magnetic white dwarfs with
MG being 6%.
Recently, Liebert et al. (2003) have found that the incidence of magnetism at
the level of
2 MG or greater is at least
10%, or higher. They
suggest that the total fraction of magnetic WDs may be substantially higher
than 10% due to the limited spectropolarimetric analyses capable of detecting
lower field strengths down to
10 kG.
Our 3 detections out of 12 objects seem to indicate that low magnetic fields
on white dwarfs (<10 kG) are frequent while high magnetic fields are
relatively rare. However, with only three detections this hypothesis remains
insecure. If confirmed by future observations, the investigation of weak
magnetic fields in white dwarfs could form a cornerstone for the future
investigation of the properties and evolution of stellar magnetic fields.
Our sample of white dwarfs is too small to discuss in detail the dependence
of the magnetic field strength on the stellar parameters (masses and cooling
ages). However, two of our detections (WD 0446-789 and
WD 1105-048) have masses of only 0.5 .
This means that their
progenitors on the main-sequence had less than 1
(Weidemann 2000). These two stars are therefore very different from the
majority of white dwarfs with megagauss magnetic fields which tend to have
higher masses (Liebert 1988; Greenstein & Oke 1982) and, therefore, high-mass
parent stars.
Measurements of weak magnetic fields are now possible for many white dwarfs with the new large telescopes, which allow a magnetic field function (MFF, in analogy to the mass function) to be constructed in the 1-100 kG range once a sufficient number of detections have been made. Such a MFF can be compared to the corresponding function for main-sequence stars (Bychkov et al. 1997) and will provide input for answers to the following key questions on the evolution of magnetic fields in stars: Are the magnetic fields in white dwarfs just the fossil relics of magnetic main-sequence stars strengthened by contraction due to conservation (to a large extent) of magnetic flux? Or do the magnetic fields develop considerably through the final stages of stellar evolution? Are the strongly magnetic white dwarfs a distinct class of objects or do they represent a tail of the distribution of magnetic fields present in all white dwarfs? Is there a dependence between magnetic field strength and mass as found in the case of magnetic WDs with higher field strengths? Do the magnetic field strengths correlate with temperature, which would be a hint for a decay on the white dwarf cooling sequence?
Several authors have suggested that the frequency of magnetic white dwarfs may increase with decreasing effective temperature, luminosity and with increasing cooling age (e.g. Liebert et al. 2003; Valyavin & Fabrika 1998), and may decrease sharply with distance (Fabrika & Valyavin 1998).
Alternatively to the fossil origin of the magnetic field in white dwarfs,
Markiel et al. (1994) and Thomas et al. (1995) have shown that a weak
magnetic field of 1.3 kG in the variable DB star GD 358 can be
explained by an
dynamo. The magnetic field in this star has
been inferred indirectly by analyzing the g-mode oscillation spectrum taken
with the WET (Whole Earth Telescope, Winget et al. 1994). However,
according to the atmospheric parameters, the convection zone in all of our
sample stars should be too shallow to support an
dynamo.
Another long-standing problem in white dwarf research is the issue of why
metals are accreted by helium-rich white dwarfs in the range
K during the passage through an interstellar
cloud while almost no hydrogen is brought into the white dwarf atmosphere.
Illarionov & Sunyaev (1975) suggested that fields below 105 G provide a
screening mechanism to separate hydrogen and ionized species from grains in
white dwarfs accreting from the interstellar matter. Since this phenomenon
always occurs in this type of star, it is possible that at this level all
white dwarfs contain magnetic fields. Friedrich et al. (2004) have searched for
circular polarisation in one DBZ and one DBAZ which have accreted metals, but
three or four orders of magnitude less hydrogen than expected.
In one case (L745-46A) a magnetic field of 7 kG (1
error of
2 kG; 99% confidence interval of
6 kG) was found, which,
however, was based on the H
line only. For the second object (GD 40),
only an upper limit of 12 kG (99% confidence) could be derived from
polarization measurements around three spectral lines. Theoretically,
magnetic fields of 4 kG and 250 kG, respectively, would be required for
the screening by the propeller mechanism to be efficient around the
two stars.
At our signal-to-noise ratio, magnetic fields down to about 2 kG can be measured. There is still the possibility that all magnetic white dwarfs contain surface magnetic fields at the 1 kG level. To test this hypothesis, much longer exposure times would be necessary, even with the VLT.
Acknowledgements
We gratefully acknowledge useful comments on the effective Landé factor by E. Landi degl'Innocenti. We acknowledge the use of LTE model spectra computed by D. Koester. We thank the staff of the ESO VLT for carrying out the service observations. Work on magnetic white dwarfs in Tübingen is supported by the DLR grant 50 OR 0201. R.N. acknowledges support by a PPARC Advanced Fellowship. We would like to thank the referee J. D. Landstreet for valuable suggestions.