A&A 396, 131-141 (2002)
D. N. Monin - S. N. Fabrika - G. G. Valyavin
Special Astrophysical Observatory of Russian AS, Nizhnij Arkhyz 369167, Russia
Received 26 June 2000 / Accepted 30 August 2002
The first results of a systematic search for magnetic fields in the brightest upper main sequence (MS) stars are presented. The main goal is to survey the stars with about the same detection limit and to improve existing statistics of their magnetism. The target list contains 57 upper MS stars and represents well B0.5-F9 stars. High-resolution Zeeman spectra were obtained for 30 stars of the list. The accuracy of the magnetic field measurements ranges from 20 to 300 G depending mainly on spectral class. In the majority of studied stars we did not detect magnetic fields. In some stars we suspect the presence of a weak magnetic field. These are the best candidates for more extensive studies. A particular case is the star Dra where we probably detected the global magnetic field. The longitudinal field strength is G. Further observations of this star are needed to confirm the detection and to ascertain if the magnetic field is variable with the period of rotation.
Key words: stars: magnetic fields - surveys
Magnetic fields have been discovered in many types of stars having different evolutionary status: pre-MS stars, MS stars, subdwarfs and white dwarfs. Often they are definitely and repeatedly detected in chemically peculiar (CP) MS stars as well. Longitudinal magnetic fields of MS stars range from 102 to 104 G. Some white dwarfs show strong magnetic fields in the range 104-109 G (Schmidt & Smith 1995; McCook & Sion 1999; Fabrika & Valyavin 1999). There is evidence that the field evolves during a star's life. It seems to appear when a star completes a significant part of its MS lifetime (Hubrig et al. 2000), and the field strength changes across the MS (Glagolevskij 1987). Valyavin & Fabrika (1999) found that the frequency of magnetic white dwarfs begins with 2% among hot (young) degenerate stars and reaches 20% among cool ones. This means that their initial fields are amplified with the time. However the nature of these initial fields is still unclear. The field in white dwarfs can be produced in the process of collapse of their progenitors, MS stars which possess surface fields of 1-104 G, by simple flux conservation. This hypothesis is based mainly on the coincidence of magnetic fluxes. To study this subject the general magnetic properties of both MS stars and white dwarfs should be analyzed and compared.
The first statistical investigation of magnetic fields in white dwarfs was done by Angel et al. (1981). Their method has been improved by Fabrika et al. (1997) and Fabrika & Valyavin (1999). They introduced the magnetic field function (MFF) - a normalized distribution of the star's frequencies versus magnetic field strength calculated per one Gauss. Now we have a powerful tool to study the incidence of magnetism in MS stars.
We need to know how frequent magnetic stars are, how many of them are weakly magnetic, whether they are a peculiar class of stars or they are just a magnetic tail in a single distribution of magnetic fields. Magnetic fields and surface chemical peculiarities in Ap stars are known to be closely related, but the mechanism of this relation is unclear. It is still unknown what the minimum magnetic field strength is that can introduce the mechanism and how it depends on other parameters of the stars. In the study of general magnetic properties one may hope to find magnetic field evolution both on the MS (Hubrig et al. 2000), and after that in subsequent evolution to the white dwarf stage.
After the first observations by Babcock (1958), many authors have published results of their studies of magnetic fields in upper MS stars. These results have been compiled by Romanyuk (2000). The measurements are not homogeneous both in methods (see for review Landstreet 1980), accuracies and selection criteria. The magnetic field measurements have so far been made mainly in stars that are likely to possess magnetic fields, sharp lined CP stars. Furthermore, authors do not often publish null results. In the next section we estimate the incidence of magnetism among MS stars using the data available.
The most extensive study of normal MS stars made with high precision is by Landstreet (1982). He searched for longitudinal magnetic fields in 31 O-F upper MS stars and 5 Am stars. No magnetic fields were detected with a median accuracy of about 65 G. Both a Balmer-line Zeeman analyzer and a two-channel coude line-profile scanner equipped with a polarization modulator were used in this study. Only 13 out of 31 stars were observed with the first device, both slow and fast rotating stars. The second device was applied to obtain helium- and iron-line measurements for slow rotating stars. Thus, the full sample cannot be considered to be unbiased.
To apply statistical methods one needs an unbiased (selection-free) sample of stars. It is obvious that the number of selection criteria must be diminished as much as possible. It seemed appropriate to undertake a systematic magnetic survey of upper MS stars. We started such a survey using the same observing method and data reduction for all samples. In this paper we report the first results. Special attention was paid to selection criteria. In Sect. 3 we describe the selection criteria, observational limitations, and the list of program stars. In Sect. 4 we present the steps of data reduction. Using the method of weighted coefficients along with the continuum reconstruction procedure we have obtained a high accuracy of magnetic field measurements.
The first attempt to construct the surface magnetic field distribution of Ap stars was made by Preston (1971). Only extremely slowly rotating Ap stars were chosen in his study. Estimates of mean surface fields were obtained from measurements of magnetically resolved or broadened absorption lines. The distribution has a maximum between 2 and 3 kG and a tail of stronger fields extended beyond 10 kG. The number of stars with a field below 2 kG is decreased. The behavior below 2 kG was explained by observational limits. In the low-field region the measurements become highly uncertain.
In their extensive study, Mathys et al. (1997) constructed the distribution of Ap stars using surface magnetic fields derived from measurements of magnetic splitting of the FeII 6149.2 line. This distribution has almost the same form as that obtained by Preston (1971). It has a maximum between 3 and 5 kG. Mathys et al. (1997) pointed out an insufficient number of stars with a field less than 3 kG. The lowest field found by them is 2.7 kG. This could not be explained by observational limitations only because the resolution limit corresponded to 1.7 kG.
It can be concluded from the analysis of both these distributions that the magnetic star frequency has a maximum at 3-5 kG in the range of surface magnetic fields from a couple of kG to a few tens of kG. Nevertheless, Bohlender & Landstreet (1990a) and Landstreet (1992) have found a continuous distribution of Ap stars in rms longitudinal magnetic fields strongly skewed toward small field values down to the limit of detectability and a long high-field tail. Note that the last authors used estimates of longitudinal fields derived by another method, the Balmer-line Zeeman analyzer technique.
In our opinion such an inconsistency is not surprising because of the direct limitations and specific selection criteria. To enhance the probability of detection, targets are isolated using chemical peculiarities, narrow lines, etc. There is a relation between magnetic field strength and a chemical peculiarity value (North 1980; Cramer & Maeder 1980). Up to the surface field kG, stronger fields are accompanied by higher peculiarities. The stars with a spectral (and photometric) peculiarity are usually used as targets to find magnetic fields. The data on stars with strong fields is incomplete, but they are dramatically incomplete for stars having weak fields.
Nevertheless, we can roughly estimate the frequency of magnetic stars. There are 161 magnetic MS stars that are brighter than 7 5 and whose spectral types are in the range B5-F2 (Romanyuk 2000). We consider MS stars of luminosity class V only. According to the Hipparcos data there are 3236 B5-F2 MS stars brighter than 7 5. Therefore, at least 5% of upper MS stars seem to be magnetic. It is obvious that this number is only a lower limit because not all 3236 stars have been studied for the presence of the field. On the other hand, Egret & Jaschek (1981) estimated the frequency of CP stars using the catalogue of stellar groups. They have found that the fraction of Ap stars (Si, Cr, Eu, Sr and Hg-Mn stars) is about 11% among B7-A3 stars. One may consider this number to be an upper limit of the frequency of magnetic stars because not every Ap star being observed for the presence of magnetic field is detected as magnetic.
Estimating the frequency of magnetic stars the question has to be answered: what is the magnetic field strength above which a star can be considered as magnetic? This limiting magnetic field value depends directly on the current observational possibilities. By improving observational techniques and decreasing we find magnetic stars to be more frequent. The MFF is free of this limitation because both positive and negative detections are used in estimations of frequencies (Fabrika et al. 1997). The first MFF for MS stars was presented by Bychkov et al. (1997). In logarithmic coordinates it has been found to be almost linear in the range of surface fields G, the spectral index is . Fabrika & Valyavin (1998, 1999) have found that the distribution of white dwarfs (the surface fields 104 - 109 G) is also linear, but it has the spectral index . Both the distributions are represented well by a power law, but their spectral indices are quite different. This fact may be related to evolution effects. However, the difference in the spectral slopes may also be related to insufficient accuracy of the distribution of MS stars. Bychkov et al. (1997) noted that the distribution can have a different behavior in the low-field part where the surface field is less than a few kG. Furthermore, Bychkov et al. (1997) pointed out that constant term could not be determined correctly because the frequency of magnetic stars is not known well. The authors used only available published data. These data are greatly distorted by observational selection.
To investigate the magnetic properties of MS stars as a whole, it is necessary to use an unbiased sample of stars with a high accuracy of measurements and to cover as many stars as possible. To get better statistics, one should observe many stars, rather than to study some of them in detail. There is a probability of null detection of the effective magnetic field in a single observation because of stellar rotation. To take into account the probability of null detection, a missing factor should be introduced (Schmidt & Smith 1995; Fabrika & Valyavin 1999). The upper MS stars show large-scale, rather dipolar magnetic fields. In dipolar geometry the missing factor depends on the accuracy of measurements and the interval of effective magnetic fields considered. Therefore, it can be calculated. This allowed us to choose a "single observation'' strategy.
There are two obvious ways to select targets - volume limited and magnitude limited. The apparent density of hot stars is low compared to cool stars. The nearest early B-type MS star is 31 pc from the Sun, and 108 F-type MS stars occupy this volume. Apart from great differences in sample sizes for cool and hot stars, the volume-limited selection also introduces differences in accuracy of magnetic field measurements. Because target magnitudes spread over a wide range, the difference of magnitudes can reach . This cannot be compensated by adjustment of exposure times.
Another way is to observe all stars up to the limiting stellar magnitude. It increases the apparent density of hot stars with respect to that of cool stars. Hot and cool stars are selected inside volumes that differ in size. This may introduce some bias, if the stars were not well mixed in space after their birth. They may have some from the common properties host.
The obvious observational selection which we cannot avoid in any selecting scheme is the different accuracy of magnetic field measurements in stars of different spectral classes. For fast rotating A-type stars with few measurable lines in their spectra the accuracy will be essentially worse. A natural way to solve the problem is to increase the number of stars involved in the study in order to test different spectral types separately. From the above arguments, the magnitude-limited selecting scheme seems more appropriate.
We formed our list using the Bright Star Catalogue (BSC, Hoffleit & Jaschek 1982). The criteria for stars to be included in the sample are as follows:
Our sample represents well MS stars in spectral classes B0.5-F9. The upper spectral type limit B0.5 is defined by the availability of stars brighter than 4 0. The northern stars cooler than F9 are mostly giants, and only 4 of them have luminosity class IV. They do not satisfy our criteria. There are 13 B stars, 36 A stars (including 6 Ap and 3 Am) and 8 F stars in the list (including one of the coolest CP star CrB). All stars may be considered as belonging to the upper MS stars where one expects global magnetic fields to be present and the convection to be not so strong as to destroy the field. The fraction of Ap stars in the list among stars in the same range of spectral types is 16%. The star UMa is probably not an Ap star as reported in BSC. It is classified as having a doubtful nature in The General catalogue of Ap and Am stars (Renson 1991). Our spectra of UMa do not show any signs suggesting it is an Ap. Exclusion of UMa reduces the fraction of Ap stars to 14%. It is close to the frequencies found by North & Cramer (1981), Egret & Jaschek (1981) and Schneider (1994).
The distributions over the sky of the selected stars, stellar magnitude,
distance and their projected
rotational velocity are shown in Fig. 1.
|Figure 1: Distributions of the target stars versus a) stellar magnitude, b) over the sky, c) versus distance and d) rotational velocity.|
|Open with DEXTER|
Analyzing all published data we find that 4 out of 57 stars of the list
have well established magnetic fields ranging between 0.1 and 1.5 kG
( UMa, Aur, CVn, CrB).
This suggests that at least 7% of upper MS stars are magnetic.
On the other hand,
24 stars of the list have already been tested for the presence
of a field. Therefore, the frequency is less than 17%.
It is an upper limit because the stars already studied were
the most promising, in expectation of the field.
Both the upper and lower limits are rough, but they
do not contradict the estimates obtained in Sect. 2.
One may expect the real frequency
to be closer to the lower limit (7%), because
every star among these 57 has been studied spectroscopically.
Probably no other strongly magnetic Ap stars can be found among our targets.
|19356||Per||2.12||B8V||65||0||320||Borra & Landstreet (1973)|
|40312||Aur||2.62||A0pSi||49||18||Borra & Landstreet (1980)|
|58715||CMi||2.90||B8Ve||276||200||Barker et al. (1985)|
|60179||Gem||1.98||A1V||14||0||37||Glagolevskij et al. (1989)|
|112185||UMa||1.77||A0pCr||38||18||Wade et al. (2000)|
|112413||CVn||2.90||A0pSiEuHg||29||55||Wade et al. (2000)|
|137909||CrB||3.68||F0p||19||21||Wade et al. (2000)|
|210418||Peg||3.53||A2Vp||117||170||Bohlender & Landstreet (1990a)|
The magnetic survey was carried out on the coude echelle spectrograph CEGS of the 1 m telescope of the Special Astrophysical Observatory. The instrumental configuration has been described by Musaev (1996). In our observations, incoming light passes through the Zeeman analyzer which contains rotational first quarter-wave plate, a calcite block and second quarter-wave plate. To improve the polarization measurement accuracy the analyzer can operate with two fixed positions of the first quarter-wave plate making an angle of . We used two different orientations of the first quarter-wave plate to record opposite circular polarizations on the same part of a CCD in two adjacent exposures. In this case many instrumental effects can be taken into account (Donati et al. 1990a). Left- and right-circularly polarized (LCP and RCP) Zeeman spectra are obtained simultaneously on the Wright Instruments CCD ( pixels). With the entrance slit set to the instrument configuration yields a resolving power of about R=40 000over a spectral range from 3700 to 9000 Å.
The spectra were taken between March 1996 and June 2000. The total number of stars observed in this period is 30. Exposure times were chosen to be longer for fainter stars so that the S/Nratio in their spectra was nearly the same as in the spectra of brighter stars.
The data reduction was carried out using the image processing package MIDAS with some modifications for polarimetric spectra (Monin 1999). The general steps are as follows: definition of spectral orders, fitting and removal of the background scattered light, cosmic ray hit removal, extraction of the spectral orders in opposite polarizations, normalization to the continuum, determination of a dispersion law, line identification, magnetic shift determination. The background is approximated by a smoothing spline interpolation between interorder points. The extraction of orders is carried out using the optimal extraction algorithm based on that described by Mukai (1990).
The transmission of LCP and RCP components through the spectrograph are slightly different. Therefore, a spurious broadband polarization signal arises. Its value changes both along the spectral order and from one order to another. As was mentioned above, to reduce most of the instrumental effects we make two exposures with different orientations of the first quarter-wave plate. Nevertheless, a small difference between LCP and RCP fluxes averaged in two exposures may exist if the exposures are obtained under slightly different conditions. Thus we need to normalize spectra or define well the continuum shape. Precise magnetic shift measurements require an accurate continuum shape reconstruction. The continuum fitting is made using an algorithm based on a smoothing-and-clipping high-frequency filter (Shergin et al. 1996). We use a convolution with a Gaussian profile for the smoothing. The algorithm takes the noise level into account.
Broad absorption lines make it difficult to determine the underlying continuum. The full width of hydrogen lines may reach tens of Angströms and line wings may spread far beyond the spectral order. The line wings can appear even in adjacent orders. To restore the true continuum flux in the region containing a broad line, we interpolate the continuum flux in several adjacent orders. A low degree (2nd or 3rd order) polynomial interpolation is applied. We use the same fitting parameters for both polarizations. It is difficult to reconstruct the continuum flux near the Balmer and Paschen jumps. Beyond the Balmer jump the spectra' quality is not sufficient for our goals. In the red part (longer than 7000 Å) almost all orders contain strong telluric bands. These orders cannot be used in the reconstruction procedure. In all other spectral regions the algorithm works well even in the case of a complicated instrumental response function.
A spectrum of scattered solar light is used as a reference spectrum for wavelength calibration. The spectrum is corrected for both Earth's rotation and orbital motion. To reach high accuracy in the determination of the dispersion law we identify at least 20 lines uniformly distributed over the image. We use the line list from Pierce & Breckinridge (1973). The dispersion relation can be expressed as: , where m is the absolute order number and x is the pixel position along the direction of the main dispersion. f(x,1/m) is a polynomial function. The standard deviation in individual line positions was typically about 0.003 Å. In the worst cases it did not exceed 0.007 Å.
For spectral line identification we use the line lists formed by the database VALD (Piskunov et al. 1995; Kupka et al. 1999; Ryabchikova at al. 1999). This database gives also the latest effective Lande factors of lines of many chemical elements (including rare-earth ones).
Magnetic shifts are measured using two different methods: the
center-of-gravity, the classical method developed for photographic
technique, and the Fourier cross-correlation method.
In the first case, we obtain the displacement of the center of gravity
of spectral lines in opposite circular polarizations:
In fact, the accuracy of magnetic shift measurements depends on many parameters. The most important are (in order of importance) the number of lines measured, the spectra quality (signal-to-noise ratio S/N), the line intensity and width (), and the method used for magnetic shift measurement. Less important are the line shape, blending, accuracy of Lande factor of individual lines, instrumental polarization and depolarization.
Taking into account the accuracy of individual line measurements
improves considerably the resulting accuracy.
Donati et al. (1997) showed that extremally small Zeeman signatures
could be detected from the analysis of many spectral lines.
They suggested analysis of a mean Stokes V profile, reproduced by
the Stokes V profiles of all selected lines. The optimal weighting
is central depth of the line.
Instead of averaging individual Stokes V profiles, we calculated the
magnetic shift for each individual line using Eq. (1).
The measurements of individual lines were averaged
or S/N as a weight w:
The central depth of hydrogen lines does not vary significantly from line to line. Thus, it is better to weight magnetic shifts of hydrogen lines by the S/N ratio, w = S/N. H has a maximum weight. The throughput of the spectrograph and detector decreases toward blue orders. Hydrogen lines which are close to the Balmer jump give a smaller contribution to the weighted average value. Pashen lines also have lower weight because of spectra quality degradation in orders with wavelengths longer than about 7000 Å.
Using weighted values we improve the precision of magnetic field estimates. In some cases, for example, when the number of measurable spectral lines is small, it allows us to reduce the uncertainty by a factor of 2.
In spectra of rapidly rotating stars the magnetic shifts are measured in a central part of the line profile (about Å). Line wings are less sensitive to the magnetic shift because of a reducing of the profile gradient. Therefore, taking a broader band around the line center, we do not improve the accuracy. In broad lines even small errors in the continuum restoration may introduce a spurious polarization signal. This signal becomes sufficient when we analyze a considerable part of the line profile. On the other hand, we get worse accuracies if we take the band too narrowly. The line center is essentially flattened by fast rotation. We have found an optimal range, Å, which provides the best accuracy.
We checked our measurements for spurious polarization using telluric lines. When the number of identified telluric lines is 100-300 we can reach an extremely high precision of about a few G. We concluded that no spurious polarization signals due to the influence of uncorrected instrumental effects are present.
We have observed two well-known magnetic stars, 53 Cam and UMa. 53 Cam is an Ap star whose longitudinal magnetic field varies between -5 and +5 kG, the rotational period is (Hill 1998). It is fainter ( ) than our targets. Nevertheless, we have obtained good quality spectra in the region 5000-6500 Å in different phases of the rotational period. Two estimates were made near the positive extreme while two others were obtained when the longitudinal magnetic field reversed its sign and became negative. About 30 unblended lines with a particularly simple Zeeman pattern were selected. Most of them are Fe I lines (70%). All measurements were phased according to the ephemeris of Hill et al. (1998). Our measurements agree well with published data (see Fig. 2).
One of the goals of our survey is to search for magnetic fields with extremely high accuracy. In Fig. 3 we present the comparison of our measurements of the weak-field star UMa with those obtained by other investigators: the Balmer-line Zeeman analyzer measurements reported by Bohlender & Landstreet (1990b), and LSD measurements by Wade et al. (2000). UMa has the weakest global magnetic field, its strength varies between approximately -100 and +100 G while the star rotates. Most of the time the longitudinal magnetic field of UMa is below the limit of detectability which can be reached by the polarimetric instruments made so far. The measurements were phased according to the ephemeris of Provin (1953). Our estimates, and G, agree well with the previous determinations (see Fig. 3), and a very good internal consistency between the two independent measurements exists. We note that internal consistency of individual measurements in those stars which have been repeatedly observed is always less than an error estimated using Eq. (3). Therefore, we suppose that the accuracy of our magnetic field measurements may be even better than we indicate. The magnetic shift in the H core of UMa corresponds to a slightly lower magnetic field, +89 and +80 G. Such a difference may be explained in terms of the concentration of Fe and Cr in the magnetic poles (Rice & Wehlay 1991). Our measurements were obtained near a phase when the positive magnetic pole was seen.
We can conclude that our data reduction software behaves properly.
measurements of known magnetic stars are in agreement
with those obtained by other authors. We are able to detect both
strong and very weak stellar magnetic fields.
|Figure 2: Longitudinal magnetic field of 53 Cam against rotational phase. Open circles - the data from Hill et al. (1998). Open square - the data from Borra & Landstreet (1977). Open diamonds - the data by Wade et al. (2000) obtained from LSD profiles. Filled circles - our measurements. The fitting curve by Hill et al. (1998) is presented. The error bars that correspond to one standard deviation are shown.|
|Open with DEXTER|
|Figure 3: Longitudinal magnetic field of UMa against rotational phase. Open circles - the data from Bohlender & Landstreet (1990b). Crosses - the data from Donati et al. (1990b). Open diamonds are measurements obtained from LSD profiles by Wade et al. (2000). Filled circles - our measurements. The error bars that correspond to one standard deviation are shown.|
|Open with DEXTER|
|Figure 4: The measurement uncertainties of against the number of measured lines. The solid line is the best fit of the points.|
|Open with DEXTER|
We have measured longitudinal magnetic fields in 30 of the 57 isolated stars. Results of the measurements are presented in Table 2. The columns give in consecutive order the name or HD number, heliocentric Julian date of the middle of an observation, the longitudinal magnetic field with corresponding errors evaluated according to Eq. (3), and the number n of spectral lines used. For some stars the results for hydrogen lines are presented separately. Both the center-of-gravity ( ) and the Fourier cross-correlation ( ) measurements are shown. Each value was obtained using two exposures with a different orientation of the first quarter-wave plate of the analyzer.
The accuracy varies from star to star. The main factor affecting the precision is the number of lines measured n. In Fig. 4 we show the measurement uncertainties against n, plotted from the data in Table 2. When the number of lines is less than 10 or more than 200 the accuracy is roughly independent of n. It appears that the smallest error which can be reached is 10-15 G. It is almost independent of the brightness of a star, because we have taken spectra of about the same quality for both bright and faint stars. A residual instrumental effect that cannot be removed even by the differential measurements lead to the precision limitation at the level of 10 G. When the number of measurable lines is small n < 10 (most of the lines have low contrast and only hydrogen lines can be analyzed), the precision tends to be around 200-300 G. Spectra' quality (the S/N ratio) drops rapidly from H to higher Balmer lines. Thus H has a weight which is considerably greater than those of other hydrogen lines, and inclusion of other hydrogen lines does not change the result significantly. This explains the behavior of the errors in Fig. 4 at small n.
For photon noise limited measurements, the standard error is given by . From Fig. 4 we find when n ranges from 10 to 100. The observed slope is steeper than expected. It is caused by the fact that contribution of cool stars increases when n grows. The cool stars, whose spectra are rich in lines, are usually slow rotators, and they have higher line contrast. The field measurement error depends on the line contrast. This results in a higher accuracy than expected from the square root dependency.
In the majority of stars in our list we did not detect magnetic fields.
However in some stars one may suspect a positive detection.
They are the star Leo with distinct features of chemical
anomalies of the Am type (Savanov 1985),
the star Cyg showing the presence of variable Be-like emission
the fast rotating star Aql,
and one of the coolest star in the list, Dra.
We intend to discuss the most interesting stars
in our next paper. Here we discuss in detail two stars, Vir and
Dra. These stars are very similar in temperature, brightness and
but while the first star shows a null result,
there is evidence of a weak field in the second one.
|Name||JD +||All lines||Hydrogen lines|
|2 450 000||, G||n||, G||n||, G||n||, G||n|
The first attempt to detect a magnetic field in this star was made by Gray (1984) by searching for Zeeman broadening in spectral lines. He detected no magnetic field. Andreetta & Giampapa (1995) concluded that only 7% of the surface of Vir may be covered by magnetic regions. They estimated the fractional area covered by active, magnetic regions.
We used more than 400 lines for the measurements. Equivalent width of the identified lines ranges from 20 to 300 mÅ. The list of chemical elements having the most numerous lines is presented in Table 3. It contains the number of lines and magnetic shifts derived from lines of the given species. Apart from those elements we used also spectral lines of Ti, Mn, Sc, Co, Cr. The magnetic field obtained from all the lines is null with an accuracy of G (see Table 2). Both the center-of-gravity and the Fourier cross-correlation methods give the same result.
A large number of lines in the spectra of Vir allowed us to
check possible systematic errors of our magnetic field estimates.
We studied the magnetic shifts versus various
spectral line parameters such as Lande factor, wavelength,
equivalent width, and equivalent width error.
In all the cases we did not find any correlation, the correlation
coefficients being less than 0.05.
We conclude that there are no notable systematic errors in our data.
Dra is a spectroscopic binary. Its orbital period is 280.55 days. The angular separation between components is (Tomkin et al. 1987; Schoeller et al. 1998). Tomkin et al. (1987) reported elements of both spectroscopic and interferometric orbits. At the time of our observations spectral lines of the primary and secondary components were unresolved (the orbital phase 0.4). The radial velocity in our composite spectra km s-1 is close to the systemic velocity km s-1. Because of the high orbital inclination eclipses in the binary do not occur. This means that spectra near the phase 0.4 consist of lines of both components and the positions of lines of both stars coincide. Spectral classes of the stars were determined by Tomkin et al. (1987) as F7V for the primary and G8V (or K2V) for the secondary.
Gray (1984) tried to detect the field in Dra searching for Zeeman broadening in 16 spectral lines. The magnetic field was not detected in that study. We cannot compare our data with the results by Gray (1984) because he did not publish an accuracy.
We have identified 124 telluric lines between 5800 and 9000 Å. The central depth ranges from 0.05 to 0.9. We used these lines to check the magnetic field measurements. Telluric lines are very narrow with steep profile gradients. No shifts in the telluric lines were found. Accepting the Lande factor equal to 1 we found the formal magnetic shift in telluric lines G and G. This fact along with others allows us to exclude any systematic instrumental effects.
The spectrum has numerous sharp lines (
We identified more than 500 lines
in spectral region 4000-8700 Å.
Most of them are Fe lines (59%). Fe I lines are much more
numerous than Fe II lines (4%).
In Table 4 we present results of measurements
by the center-of-gravity method
of those elements that have significant number of lines in the spectra of
Dra. It corresponds to the second
measurement according to Table 2, JD 2 450 326.60.
Ti (9% of all identified lines), Ni (8%), Cr (5%), Ca (3%) and Mn (3%) lines do not show reliable magnetic shift. Fe I lines give the main contribution to the magnetic field estimates, and . Using Fe II lines we found a higher field strength, G and . The difference between two average values -44 G for the Fe I lines and -147 G for the Fe II lines is about 3 error bars. This difference could be, in principle, a statistical effect. On the other hand, the primary and secondary lines in Dra were unresolved in our observations. The secondary has an effective temperature of 4900 K which is less than temperature of the primary by 1200 K (Tomkin et al. 1987). Using speckle interferometric data Schoeller et al. (1998) reported that the stellar magnitude of the components differs by about 2 in the red (6500-7500 Å). Some spectral lines are formed almost completely in the hotter star, but the cooler component becomes noticeable in others. The Fe II lines are formed mainly in the atmosphere of primary. Therefore, the difference in the field strength mentioned above may also be due to depolarization by the lines of the secondary star.
To check this hypothesis we found a line ratio
(secondary over primary) in the composite spectrum. We used line center
intensities from two models,
The models were taken from the VALD database. We found the line ratio for
each line identified in the spectrum. The flux difference of two components
(by about a factor of 6 in the red and 9 in the blue part of the spectra)
was taken into account.
is small it means that the line is formed mainly in the
lines are equally due to primary and secondary.
If the primary star is magnetic and secondary is non-magnetic or it has
a longitudinal field of opposite polarity, we may expect a decrease of the
field strength derived from individual lines with increasing of .
It may even change its sign at some if the field of the secondary is strong enough.
In Fig. 5 the longitudinal magnetic field of individual
lines derived from the second pair of spectra (JD 2 450 326.60)
is plotted against .
|Figure 5: The longitudinal magnetic field of Dra against the secondary to primary line ratio for all identified lines. Open symbols - Fe II lines. Mean values of were calculated inside the range marked by horizontal bars. error bars are shown.|
|Open with DEXTER|
We also apply another method to detect the field in the spectra
of Dra. Spectral line magnetic splitting
depends on the Lande factor (
line wavelength. We can introduce a reduced magnetic shift
We obtained the magnetic field of about G from measurements of the H core ( Å from the line center). The same field strength was registered by both the center-of-gravity and the cross-correlation methods. This fact gives independent confirmation of the magnetic field detection.
There are several independent indications that Dra is very
likely to possess a weak global magnetic field of about a few tens G.
This is the result of direct measurements of line shifts in
two pairs of spectra, the behavior of the reduced magnetic shift with
effective Lande factors, the different values of the field in the Fe I and
Fe II lines which are distributed in the same spectral region, the firm
negative detection in telluric lines also distributed in the same
spectral region and the magnetic shift in H.
A more detailed study of this star is beyond the scope of the present
Further time-resolved spectropolarimetric observations are needed
to confirm the positive detection and to ascertain
whether the longitudinal magnetic field varies with the star's rotation.
The weak field in Dra (as well as possibly in some
other stars of our list) would be very important in finding a connection
between the upper MS stars which have large-scale organized magnetic
fields and G, K, M dwarfs having complex magnetic field structures, in
studying a connection between weak magnetic fields in the upper MS stars
and weak chemical peculiarities in their atmospheres.
|Figure 6: The reduced magnetic shift in the spectrum of Dra binned in 4 intervals of the effective Lande factors. The error bars correspond to one standard deviation. The best linear fit forced through the origin ( ) is shown.|
|Open with DEXTER|
The first results of the systematic search for magnetic fields in the brightest MS stars show that the accuracy of our magnetic field measurements is sufficient to fulfill the main goal - to survey an unbiased sample of upper MS stars. The accuracy of the measurements ranges from 20 to 300 G, depending mainly on spectral type. The spectra in this study were taken under practically identical conditions using the same spectrograph and detector. We are continuing the survey observations. The observations of all targets of the list will help us to derive the selection-free distribution of magnetic fields of upper MS stars.
In some stars we suspect the presence of a weak magnetic field. These suspected stars are the best candidates for more extensive studies and additional observations of these stars are needed. A particular case is the star Dra, where we have a positive detection of a magnetic field. Further observations of Dra are needed to confirm the detection and to ascertain if the field varies with stellar rotation.
The authors thank E. A. Barsukova for assistance in performing the observations, J. D. Landstreet and I. I. Romanyuk for useful discussions. We thank the referee, J.-F. Donati, for his helpful critical comments. The VALD database operated in Vienna, Austria and the Simbad database operated at the CDS, Strasbourg, France were used. The work was supported by a grant of the Russian AS for young researchers (Project N. 8) and partly by the RFBR grant N. 01-02-16808.