Issue |
A&A
Volume 557, September 2013
|
|
---|---|---|
Article Number | L16 | |
Number of page(s) | 4 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/201321779 | |
Published online | 12 September 2013 |
Discovery of a magnetic field in the CoRoT hybrid B-type pulsator HD 43317⋆
1 Institut d’Astrophysique et de Géophysique, Université de Liège, Allée du 6 Août 17, Bât B5c, 4000 Liège, Belgium
e-mail: maryline.briquet@ulg.ac.be
2 LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Univ. Paris Diderot, 5 place Jules Janssen, 92195 Meudon Cedex, France
3 Instituut voor Sterrenkunde, Katholieke Universiteit Leuven, Celestijnenlaan 200 D, 3001 Leuven, Belgium
Received: 26 April 2013
Accepted: 12 August 2013
Context. A promising way of testing the impact of a magnetic field on internal mixing (core overshooting, internal rotation) in main-sequence B-type stars is to perform asteroseismic studies of a sample of magnetic pulsators.
Aims. The CoRoT satellite revealed that the B3IV star HD 43317 is a hybrid SPB/β Cep-type pulsator that has a wealth of pulsational constraints on which one can perform a seismic modelling, in particular, probing the extent of its convective core and mixing processes. Moreover, indirect indicators of a magnetic field in the star were observed: rotational modulation due to chemical or temperature spots and X-ray emission. Our goal was to directly investigate the field in HD 43317 and, if it is magnetic, to characterise it.
Methods. We collected data with the Narval spectropolarimeter installed at Télescope Bernard Lyot (TBL, Pic du Midi, France) and applied the least-squares deconvolution technique to measure the circular polarisation of the light emitted from HD 43317. We modelled the longitudinal field measurements directly with a dipole.
Results. Zeeman signatures in the Stokes V profiles of HD 43317 are clearly detected and rotationally modulated, which proves that this star exhibits an oblique magnetic field. The modulation with the rotation period deduced from the CoRoT light curve is also confirmed, and we found a field strength at the poles of about 1 kG. Our result must be taken into account in future seismic modelling work of this star.
Key words: stars: magnetic field / stars: individual: HD 43317
© ESO, 2013
1. Introduction
The study of the magnetic properties of pulsating B-type stars − β Cep and slowly pulsating B (SPB) stars − is particularly interesting because when it is combined with the study of their pulsational properties, it provides a unique way of probing the impact of magnetism on the physics of non-standard mixing processes inside these stars. Comparing the amount of mixing obtained by asteroseismic investigation for a sample of magnetic β Cep stars with that of a sample of non-magnetic objects would already allow us to corroborate or disprove that magnetic fields inhibit mixing in stellar interiors. Moreover, the surface field strength obtained from modelling the magnetic measurements is an additional constraint for a comparison with the internal field strength needed to inhibit mixing in the radiative zone, as obtained by different criteria (e.g., Mathis & Zahn 2005).
Up to now, asteroseismic modelling has been performed for two magnetic β Cep stars only: β Cephei itself (Shibahashi & Aerts 2000) and V2052 Oph (Briquet et al. 2012). For V1449 Aql, which was also modelled asteroseismically (Aerts et al. 2011), a magnetic field was detected by Hubrig et al. (2011), but Shultz et al. (2012) disputed the presence of this field. Magnetic fields have also been detected in a number of other β Cep and SPB stars (e.g., Donati et al. 2001; Hubrig et al. 2006, 2009), but there is no seismic modelling available for them. Moreover, the approach presented in the previous paragraph has been accomplished only for V2052 Oph so far (Briquet et al. 2012; Neiner et al. 2012a). Asteroseismic investigations of this magnetic target, which has a polar field strength of about 400 G (Neiner et al. 2012a), showed that the stellar models that explained the observed pulsational behaviour needed no convective core overshooting (Briquet et al. 2012). This outcome is in contrast to other results of dedicated asteroseismic studies of non-magnetic β Cep stars (e.g., Briquet et al. 2007): it is usually found that convective core overshooting needs to be included in the stellar models (Aerts et al. 2010). The most plausible explanation is that the magnetic field inhibits non-standard mixing processes inside V2052 Oph. Indeed, the field strength observed in this star is higher than the critical field limit needed to inhibit mixing that was determined from theory (see Zahn 2011; Neiner et al. 2012a).
A recent study of the B3IV single star HD 43317 revealed it to be a promising target to explore the effects of magnetism on stellar interiors of main-sequence B-type stars in more detail. Indeed, an analysis of the CoRoT light curve of this star showed p and g pulsation modes of β Cep and SPB types, as well as two series of consecutive frequencies with an almost constant period spacing near 6300 seconds (Pápics et al. 2012, hereafter P12). Such an observation carries information on the extent of the convective region and mixing processes (Miglio et al. 2008), as applied to only one hybrid B-type pulsator so far (Degroote et al. 2010). Moreover, if the star is magnetic, the detected high-order g-modes have the potential to probe its internal magnetic field, as suggested by the theoretical work of Hasan et al. (2005); but this has never been applied to a massive star.
In both the photometry and spectroscopy of P12, rotational modulation connected to chemical or temperature spots were observed, providing us with a precise value of the rotation period of the star (Prot = 0.8969 days), which rotates at 50% of its critical velocity. The clear signature of rotational modulation in the He lines, as shown and discussed in P12 (see their Fig. 7), is an indirect indicator of a magnetic field in HD 43317. Indeed, it is generally believed that magnetic fields are responsible for the formation of chemical spots at the surface of chemically peculiar Bp stars (e.g., Michaud et al. 1981). Moreover, X-ray emission for our object is reported in the ROSAT all-sky survey catalogue of optically bright OB-type stars with LX ~ 8.3 × 1030 erg s-1 (Berghoefer et al. 1996). This may also point towards a magnetic field, because X-rays can be produced by shock waves in the stellar winds or by confinement of matter at the magnetic equator.
In this paper, we present the first direct detection of a magnetic field in HD 43317 by making use of the high-efficiency and high-resolution Narval spectropolarimeter installed at the Télescope Bernard Lyot (TBL) 2-m telescope (Pic du Midi, France).
2. Observations
Nineteen high-resolution (R = 65 000) spectropolarimetric Stokes V Narval observations of HD 43317 were collected in 2012 (see Table 1). The first four measurements (of 4 × 1000 s each) were obtained to investigate the presence of a magnetic field. A Zeeman signature being detected, we continued our sequence of observations to sample other rotational phases, but changed our observing strategy to consider the fact that the star is pulsating and to increase the signal-to-noise ratio (S/N). The next 15 observations consisted of 5 × 3 successive sequences (which were afterwards averaged) of 4 × 800 s each. Each sequence was shorter than 1/20th of the pulsation period to avoid that the deformations of the intensity profile by pulsations during the measurement polluted the Stokes V signal. Because the main pulsation period detected in the spectroscopy of HD 43317 is 0.735 d (P12), the exposure time is expected not be longer than 4 × 800 s at once.
Journal of Narval/TBL observations and magnetic field measurements of HD 43317.
The usual bias, flat-fields, and ThAr calibrations were obtained each night. The data reduction was performed using libre-esprit (Donati et al. 1997), the dedicated software available at TBL. The intensity spectra were then normalised to the continuum by fitting a cubic spline function. Stokes V spectra were derived by constructively combining the left- and right-hand circularly polarised light together from the four sub-exposures. The three measurements obtained successively each night were averaged so that we finally had nine magnetic measurements that covered different phases of rotation. They were normalised to the continuum intensity.
We adopted the ephemeris of HD 43317 from P12, that is, a rotation period Prot = 0.8968634 d and a reference date HJD0 = 2 455 174.57665.
We applied the least-squares deconvolution (LSD) technique (Donati et al. 1997) to the photospheric spectral lines in each échelle spectrum (λ = [3750−10 500] Å) to construct a single profile with an increased S/N (note that no additional normalisation was applied to the LSD profiles). In that way, the Zeeman signature induced by a magnetic field is clearly visible in the high S/N Stokes V profiles (Fig. 1), proving that HD 43317 is a magnetic object. To diagnose possible spurious polarisation signatures, the null profiles N were also computed by destructively combining the four subexposures. The flatness of these profiles (Fig. 1) indicates that none of our measurements has been polluted by stellar pulsations, observing conditions, or instrumental effects, even those with slightly longer individual exposure times. To compute the LSD profiles, we made use of a mask created from the Kurucz atomic database and ATLAS 9 atmospheric models of solar abundance (Kurucz 1993), for Teff = 17 000 K and log g = 4.0 (following P12), with intrinsic line depths higher than 0.1. This mask contains 199 photospheric He and metal lines of various chemical elements, together with their wavelength, depth, and Landé factor. The depths were modified so that they correspond to those of the observed spectral lines, with an interactive graphical user interface (GUI) program, written in IDL, that was kindly put at our disposal by Grunhut. The adjustment was straightforward for unblended individual lines. The relative depths of the blended lines in the original mask were kept fixed and the entire line blend depth was adjusted to the spectrum. Lines whose Landé factor is unknown or that did not appear in the spectrum were excluded. The average S/N is 5100 per 2.6 km s-1 pixels in the LSD Stokes V profiles.
Fig. 1 LSD Stokes V (top), N (middle), and I (bottom) profiles for the nine rotational phases. The number of the magnetic measurement, as listed in Table 1, is indicated on the right of each profile. Typical error bars are shown next to each I and V profile. |
3. Magnetic field
3.1. Longitudinal magnetic field
We used the LSD Stokes I and V profiles to compute the line-of-sight component of the magnetic field integrated over the visible stellar surface, that is, the longitudinal magnetic field in Gauss, given for instance by Eq. (1) of Wade et al. (2000). In this equation we used λ = 501 nm and g = 1.20. The integration limits were chosen to be high enough to cover the line width but also small enough to avoid artificially increased error bars due to noisy continuum. A range of 123 km s-1 around the line centre was adopted.
The values for the longitudinal magnetic field are reported in Table 1. They vary between about −140 and 180 G. The error bars are typically ~35 G. A phase diagram of the longitudinal field folded with the rotation period is shown in Fig. 2.
From a dipolar fit of the Bl measurements, we derived Bl,max = 193 ± 30 G, Bl,min = −115 ± 30 G, and the phase difference ΔΦ = 0.58 ± 0.02 between two successive crossings of the Bl fit at 0. The dipole fit has a reduced χ2 value of 1.3. From Eqs. (18) and (19) of Shore (1987), for instance, we then derived the obliquity angle β ∈ [70, 86] deg, assuming the inclination angle i ∈ [20, 50] deg as deduced for the star by P12, but considering the errors on the Teff and log g values.
With these values of i and β, we also deduced the polar field strength Bpol ∈ [650, 1750] G from Eq. (5) of Borra & Landstreet (1979).
Fig. 2 Longitudinal magnetic field curve of HD 43317, folded in phase with the rotation period Prot = 0.8968634 d and a reference date HJD0 = 2 455 174.57665 from P12. A sinusoidal (dipole) fit is plotted in red. |
3.2. Magnetic confinement and magnetosphere
We checked whether the magnetic field we discovered in HD 43317 is sufficient to confine the stellar wind. We used the magnetic confinement parameter, η∗, defined by ud-Doula & Owocki (2002), which depends on the stellar equatorial magnetic field strength, radius, and wind momentum. When η∗ > 1, wind material can be channelled along the field lines and be confined into a circumstellar magnetosphere.
To be conservative, we used the lower value of the polar field strength of Bpol = 650 G discovered in this work and an upper value of the mass loss of Ṁ = 10-9 M⊙ yr-1 according to Fig. 4 of Petit et al. (2013) for a B3IV star. We measured the terminal wind velocity v∞ ~ 150 km s-1 from the Si iv 1394+1403 Å and Al iii 1855+1863 Å lines available in the SWP IUE archival spectrum. Finally, we used i ∈ [20, 50] deg as derived above. We imposed that Prot ~ 0.89 d. We found that η∗ = 3400 for i = 50 deg and η∗ = 14 500 for i = 20 deg. In addition, we obtained that the Alfvén radius is RA = 7.6 R∗ for i = 50 deg and RA = 11 R∗ for i = 20 deg, and the corotating Kepler radius is RK = 2.6 R∗ for i = 50 deg and RK = 1.2 R∗ for i = 20 deg, respectively. Therefore, wind material is confined by the magnetic field (η∗ > 1) and RA > RK, that is, the magnetosphere is supported by the centrifugal force. See ud-Doula et al. (2008), Petit et al. (2012) and Neiner et al. (2012b) for more details on magnetospheres.
We did not detect emission in the Hα line of HD 43317. This means that the emission measure in the centrifugally supported magnetosphere is too low, probably because the wind is weak and rather slow. Nevertheless, we recall that HD 43317 is detected as an X-ray source in the ROSAT catalogue (Berghoefer et al. 1996) with log (LX/Lbol) = −5.94 and an estimated X-ray temperature of 0.49 keV. X-ray emission could be produced in the magnetic equator where wind particles from the two magnetic hemispheres collide.
4. Conclusions
Thanks to photometric data of unprecedented precision assembled by the CoRoT satellite complemented by ground-based spectroscopy (P12) and spectropolarimetry (presented in this paper), HD 43317 was discovered to be a single magnetic B-type hybrid pulsator with much observational information that can be used to constrain its interior properties and to study the effect of magnetism in the mixing processes that act inside the star. In future work on HD 43317, the magnetic field must be taken into account, for example to compute the oscillation periods that need to be compared with the two almost constant period spacings of high-order gravity modes observed in the star (Salmon et al., in prep.).
HD 43317 is also an interesting target in which to study physical mechanisms that take place at its stellar surface, such as the interplay between radiatively driven diffusion and magnetic field, which is expected to be responsible for the surface abundance inhomogeneities. To this end, additional spectropolarimetric data that cover the rotation cycle are required so that a Zeeman-Doppler imaging of the stellar surface can be performed, which will provide a detailed characterisation of the magnetic configuration as well as abundance mapping of the stellar surface.
Acknowledgments
We thank Evelyne Alecian for fruitful discussions. The research leading to these results has received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007–2013)/ERC grant agreement n°227224 (PROSPERITY), as well as from the Belgian Science Policy Office (Belspo, C90309: CoRoT Data Exploitation).
References
- Aerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology (Springer) [Google Scholar]
- Aerts, C., Briquet, M., Degroote, P., Thoul, A., & van Hoolst, T. 2011, A&A, 534, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Berghoefer, T. W., Schmitt, J. H. M. M., & Cassinelli, J. P. 1996, A&AS, 118, 481 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Borra, E. F., & Landstreet, J. D. 1979, ApJ, 228, 809 [NASA ADS] [CrossRef] [Google Scholar]
- Briquet, M., Morel, T., Thoul, A., et al. 2007, MNRAS, 381, 1482 [NASA ADS] [CrossRef] [Google Scholar]
- Briquet, M., Neiner, C., Aerts, C., et al. 2012, MNRAS, 427, 483 [NASA ADS] [CrossRef] [Google Scholar]
- Degroote, P., Aerts, C., Baglin, A., et al. 2010, Nature, 464, 259 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Donati, J.-F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, MNRAS, 291, 658 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Donati, J.-F., Wade, G. A., Babel, J., et al. 2001, MNRAS, 326, 1265 [NASA ADS] [CrossRef] [Google Scholar]
- Hasan, S. S., Zahn, J.-P., & Christensen-Dalsgaard, J. 2005, A&A, 444, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hubrig, S., Briquet, M., Schöller, M., et al. 2006, MNRAS, 369, L61 [NASA ADS] [CrossRef] [Google Scholar]
- Hubrig, S., Briquet, M., De Cat, P., et al. 2009, Astron. Nachr., 330, 317 [NASA ADS] [CrossRef] [Google Scholar]
- Hubrig, S., Ilyin, I., Briquet, M., et al. 2011, A&A, 531, L20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km s-1 grid. Kurucz CD-ROM No. 13 (Cambridge, Mass.: Smithsonian Astrophysical Observatory), 13 [Google Scholar]
- Mathis, S., & Zahn, J.-P. 2005, A&A, 440, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Michaud, G., Charland, Y., & Megessier, C. 1981, A&A, 103, 244 [NASA ADS] [Google Scholar]
- Miglio, A., Montalbán, J., Noels, A., & Eggenberger, P. 2008, MNRAS, 386, 1487 [NASA ADS] [CrossRef] [Google Scholar]
- Neiner, C., Alecian, E., Briquet, M., et al. 2012a, A&A, 537, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Neiner, C., Landstreet, J. D., Alecian, E., et al. 2012b, A&A, 546, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pápics, P. I., Briquet, M., Baglin, A., et al. 2012, A&A, 542, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Petit, V., Owocki, S. P., Oksala, M. E., & MiMeS Collaboration 2012, in ASP Conf. Ser. 465, eds. L. Drissen, C. Rubert, N. St-Louis, & A. F. J. Moffat, 48 [Google Scholar]
- Petit, V., Owocki, S. P., Wade, G. A., et al. 2013, MNRAS, 429, 398 [NASA ADS] [CrossRef] [Google Scholar]
- Shibahashi, H., & Aerts, C. 2000, ApJ, 531, L143 [NASA ADS] [CrossRef] [Google Scholar]
- Shore, S. N. 1987, AJ, 94, 731 [NASA ADS] [CrossRef] [Google Scholar]
- Shultz, M., Wade, G. A., Grunhut, J., et al. 2012, ApJ, 750, 2 [NASA ADS] [CrossRef] [Google Scholar]
- ud-Doula, A., & Owocki, S. P. 2002, ApJ, 576, 413 [NASA ADS] [CrossRef] [Google Scholar]
- ud-Doula, A., Owocki, S. P., & Townsend, R. H. D. 2008, MNRAS, 385, 97 [NASA ADS] [CrossRef] [Google Scholar]
- Wade, G. A., Donati, J.-F., Landstreet, J. D., & Shorlin, S. L. S. 2000, MNRAS, 313, 851 [NASA ADS] [CrossRef] [Google Scholar]
- Zahn, J.-P. 2011, in IAU Symp. 272, eds. C. Neiner, G. Wade, G. Meynet, & G. Peters, 14 [Google Scholar]
All Tables
All Figures
Fig. 1 LSD Stokes V (top), N (middle), and I (bottom) profiles for the nine rotational phases. The number of the magnetic measurement, as listed in Table 1, is indicated on the right of each profile. Typical error bars are shown next to each I and V profile. |
|
In the text |
Fig. 2 Longitudinal magnetic field curve of HD 43317, folded in phase with the rotation period Prot = 0.8968634 d and a reference date HJD0 = 2 455 174.57665 from P12. A sinusoidal (dipole) fit is plotted in red. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.