Axisymmetric magnetic fields in stars: relative strengths of poloidal and toroidal components

In this third paper in a series on stable magnetic equilibria in stars, I look at the stability of axisymmetric field configurations and, in particular, the relative strengths of the toroidal and poloidal components. Both toroidal and poloidal fields are unstable on their own, and stability is achieved by adding the two together in some ratio. I use Tayler's stability conditions for toroidal fields and other analytic tools to predict the range of stable ratios and then check these predictions by running numerical simulations. If the energy in the poloidal component as a fraction of the total magnetic energy is written as E_p / E , it is found that the stability condition is a ( E / U ) < E_p / E ≲ 0.8 where E /U is the ratio of magnetic to gravitational energy in the star and a is some dimensionless factor whose value is of order 10 in a main-sequence star and of order 10³ in a neutron star. In other words, whilst the poloidal component cannot be significantly stronger than the toroidal, the toroidal field can be very much stronger than the poloidal–given that in realistic stars we expect E / U < 10⁻⁶. The implications of this result are discussed in various contexts such as the emission of gravitational waves by neutron stars, free precession and a 'hidden' energy source for magnetars.     

Discovery of magnetic fields in three He variable Bp stars with He and Si spots

It is essential for the understanding of stellar structure models of high mass stars to explain why constant stars, nonpulsating chemically peculiar hot Bp stars and pulsating stars co‐exist in the slowly pulsating B stars and β Cephei instability strips. We have conducted a search for magnetic fields in the four Bp stars HD55522, HD105382, HD131120, and HD138769 which previously have been wrongly identified as slowly pulsating B stars. A recent study of these stars using the Doppler Imaging technique revealed that the elements He and Si are inhomogeneously distributed on the stellar surface, causing the periodic variability. Using FORS 1 in spectropolarimetric mode at the VLT, we have acquired circular polarisation spectra to test the presence of a magnetic field in these stars. A variable magnetic field is clearly detected in HD55522 and HD105382, but no evidence for the existence of a magnetic field was found in HD131120. The presence of a magnetic field in HD138769 is suggested by one measurement at 3σ level. We discuss the occurrence of magnetic B stars among the confirmed pulsating B stars and find strong magnetic fields of order kG and oscillations to be mutually exclusive. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Line profile variability of OB stars: Pulsation,rotation, clumps and magnetic fields

The variability of line profiles in spectra of bright OB stars has been studied.We obtain more than 1000 high quality spectra of δ Ori A, λ Ori A, α Cam, 19 Cep, ι Her, ρ Leo and other target stars. We revealed the line profile microvariability of small amplitude (0.5–3% in the adjacent continuum units) for all observed stars. For most stars only cyclic components of the line profile variability (LPV) at the time scales from hours to days were detected. These components seem to be connected with both the non-radial pulsations (NRP) and rotation line profile modulation. In the spectra of δ Ori A and λ Ori A we found the evidences of the stochastic LPV in spectra, probably connected with the small clumps in the stellar wind. On the basis of recent observations we discuss the origin of the magnetic field of early-type stars. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Destabilization of hydrodynamically stable rotation laws by azimuthal magnetic fields

We consider the effect of toroidal magnetic fields on hydrodynamically stable Taylor–Couette differential rotation flows. For current-free magnetic fields a non-axisymmetric   m = 1  magnetorotational instability arises when the magnetic Reynolds number exceeds   O (100)  . We then consider how this 'azimuthal magnetorotational instability' (AMRI) is modified if the magnetic field is not current-free, but also has an associated electric current throughout the fluid. This gives rise to current-driven Tayler instabilities (TIs) that exist even without any differential rotation at all. The interaction of the AMRI and the TI is then considered when both electric currents and differential rotation are present simultaneously. The magnetic Prandtl number Pm turns out to be crucial in this case. Large Pm have a destabilizing influence, and lead to a smooth transition between the AMRI and the TI. In contrast, small Pm have a stabilizing influence, with a broad stable zone separating the AMRI and the TI. In this region the differential rotation is acting to stabilize the TIs, with possible astrophysical applications (Ap stars). The growth rates of both the AMRI and the TI are largely independent of Pm , with the TI acting on the time-scale of a single rotation period, and the AMRI slightly slower, but still on the basic rotational time-scale. The azimuthal drift time-scale is ∼20 rotations, and may thus be a (flip-flop) time-scale of stellar activity between the rotation period and the diffusion time.     

A search for magnetic fields in cool sdB stars

Hot cluster horizontal branch (HB) stars and field subdwarf B (sdB) stars are core helium burning stars that exhibit abundance anomalies that are believed to be due to atomic diffusion. Diffusion can be effective in these stars because they are slowly rotating. In particular, the slow rotation of the hot HB stars (T_eff > 11000 K), which show abundance anomalies, contrasts with the fast rotation of the cool HB stars, where the observed abundances are consistent with those of red giants belonging to the same cluster. The reason why sdB stars and hot HB stars are rotating slowly is unknown. In order to assess the possible role of magnetic fields on abundances and rotation, we investigated the occurrence of such fields in sdB stars with T_eff < 30 000 K, whose temperatures overlap with those of the hot HB stars. We conclude that large‐scale organised magnetic fields of kG order are not generally present in these stars but at the achieved accuracy, the possibility that they have fields of a few hundred Gauss remains open. We report the marginal detection of such a field in SB 290; further observations are needed to confirm it (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the time dependence of differential rotation in young late-type stars

A model for the angular momentum transfer within the convection zone of a rapidly rotating star is introduced and applied to the analysis of recent observations of temporal fluctuations of the differential rotation on the young late-type stars AB Doradus (AB Dor) and LQ Hydrae (LQ Hya). Under the hypothesis that the mean magnetic field produced by the stellar dynamo rules the angular momentum exchanges and that the angular velocity depends only on the distance s from the rotation axis and the time, the minimum azimuthal Maxwell stress  | B_sB _φ|  , averaged over the convection zone, is found to range from ∼0.04 to  ∼0.14 T²  . If the poloidal mean magnetic field   B _s   is of the order of 0.01 T, as indicated by the Zeeman–Doppler imaging maps of those stars, then the azimuthal mean field   B _φ  can reach an intensity of several teslas, which significantly exceeds equipartition with the turbulent kinetic energy. Such strong fields can account also for the orbital period modulation observed in cataclysmic variables and RS Canum Venaticorum systems with a main-sequence secondary component. Moreover, the model allows us to compute the kinetic energy dissipation rate during the maintenance of the differential rotation. Only in the case of the largest surface shear observed on LQ Hya may the dissipated power exceed the stellar luminosity, but the lack of a sufficient statistic on the occurrence of such episodes of large shear does not allow us to estimate their impact on the energy budget of the convection zone.     

A new numerical scheme for structures of rotating magnetic stars

We have developed a new numerical scheme for obtaining structures of rapidly rotating stars with strong magnetic fields. In our scheme, both poloidal and toroidal magnetic fields can be treated for stars with compressibility and infinite conductivity. By introducing the vector potential and its integral representation, we can treat the boundary condition for the magnetic fields across the surface properly. We show structures and distributions of magnetic fields as well as the distributions of the currents of rotating magnetic polytropic stars with polytropic index   N = 1.5  . The shapes of magnetic stars are oblate as long as the magnetic vector potential decreases as 1/ r when   r →∞  . For extremely strong magnetic fields, equilibrium configurations can be of toroidal shapes.     

Differential rotation of cool active stars

The surface differential rotation of active solar‐type stars can be investigated by means of Doppler and Zeeman‐Doppler Imaging, both techniques enabling one to estimate the short‐term temporal evolution of photospheric structures (cools spots or magnetic regions). After describing the main modeling tools recently developed to guarantee a precise analysis of differential rotation in this framework, we detail the main results obtained for a small number of active G and K fast rotating stars. We evoke in particular some preliminary trends that can be derived from this sample, bearing the promise that major advances in this field will be achieved with the new generation of spectropolarimeters (ESPaDOnS/CFHT, NARVAL/TBL). (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic fields and rotation in white dwarfs and neutron stars

Recent spectropolarimetric observations of Ap and Bp stars with improved sensitivity have suggested that most Ap and Bp stars are magnetic with dipolar fields of at least a few hundred gauss. These new estimates suggest that the range of magnetic fluxes found for the majority of magnetic white dwarfs is similar to that of main-sequence Ap–Bp stars, thus strengthening the empirical evidence for an evolutionary link between magnetism on the main sequence and magnetism in white dwarfs. We draw parallels between the magnetic white dwarfs and the magnetic neutron stars and argue that the observed range of magnetic fields in isolated neutron stars  ( B_p ∼ 10¹¹–10¹⁵ G)  could also be explained if their mainly O-type progenitors have effective dipolar fields in the range of a few gauss to a few kilogauss, assuming approximate magnetic flux conservation with the upper limit being consistent with the recent measurement of a field of   B_p ∼ 1100 G  for θ Orion C.
In the magnetic field–rotation diagram, the magnetic white dwarfs can be divided into three groups of different origin: a significant group of strongly magnetized slow rotators  ( P _rot∼ 50 –100 yr)  that have originated from single-star evolution, a group of strongly magnetized fast rotators  ( P _rot∼ 700 s)  , typified by EUVE J0317–853, that have originated from a merger, and a group of modest rotators ( P _rot∼ hours–days) of mixed origin (single-star and CV-type binary evolution). We propose that the neutron stars may similarly divide into distinct classes at birth , and suggest that the magnetars may be the counterparts of the slowly rotating high-field magnetic white dwarfs.     

Dynamo-generated magnetic fields at the surface of a massive star

Spruit has shown that an astrophysical dynamo can operate in the non-convective material of a differentially rotating star as a result of a particular instability in the magnetic field (the Tayler instability). By assuming that the dynamo operates in a state of marginal instability, Spruit has obtained formulae which predict the equilibrium strengths of azimuthal and radial field components in terms of local physical quantities. Here, we apply Spruit's formulae to our previously published models of rotating massive stars in order to estimate Tayler dynamo field strengths. There are no free parameters in Spruit's formulae. In our models of 10- and  50-M_⊙  stars on the zero-age main sequence, we find internal azimuthal fields of up to 1 MG, and internal radial components of a few kG. Evolved models contain weaker fields. In order to obtain estimates of the field strength at the stellar surface, we examine the conditions under which the Tayler dynamo fields are subject to magnetic buoyancy. We find that conditions for Tayler instability overlap with those for buoyancy at intermediate to high magnetic latitudes. This suggests that fields emerge at the surface of a massive star between magnetic latitudes of about 45° and the poles. We attempt to estimate the strength of the field which emerges at the surface of a massive star. Although these estimates are very rough, we find that the surface field strengths overlap with values which have been reported recently for line-of-sight fields in several O and B stars.     

Evolution of magnetic fields in stars across the upper main sequence: II. Observed distribution of the magnetic field geometry

We re‐discuss the evolutionary state of upper main sequence magnetic stars using a sample of Ap and Bp stars with accurate Hipparcos parallaxes and definitely determined longitudinal magnetic fields. We confirm our previous results obtained from the study of Ap and Bp stars with accurate measurements of the mean magnetic field modulus and mean quadratic magnetic fields that magnetic stars of mass M < 3 M_⊙ are concentrated towards the centre of the main‐sequence band. In contrast, stars with masses M > 3 M_⊙ seem to be concentrated closer to the ZAMS. The study of a few known members of nearby open clusters with accurate Hipparcos parallaxes confirms these conclusions. Stronger magnetic fields tend to be found in hotter, younger and more massive stars, as well as in stars with shorter rotation periods. The longest rotation periods are found only in stars which spent already more than 40% of their main sequence life, in the mass domain between 1.8 and 3 M_⊙ and with log g values ranging from 3.80 to 4.13. No evidence is found for any loss of angular momentum during the main‐sequence life. The magnetic flux remains constant over the stellar life time on the main sequence. An excess of stars with large obliquities β is detected in both higher and lower mass stars. It is quite possible that the angle β becomes close to 0. in slower rotating stars of mass M > 3 M_⊙ too, analog to the behaviour of angles β in slowly rotating stars of M < 3 M_⊙. The obliquity angle distribution as inferred from the distribution of r ‐values appears random at the time magnetic stars become observable on the H‐R diagram. After quite a short time spent on the main sequence, the obliquity angle β tends to reach values close to either 90. or 0. for M < 3 M_⊙. The evolution of the obliquity angle β seems to be somewhat different for low and high mass stars. While we find a strong hint for an increase of β with the elapsed time on the main sequence for stars with M > 3 M_⊙, no similar trend is found for stars with M < 3 M_⊙. However, the predominance of high values of β at advanced ages in these stars is notable. As the physics governing the processes taking place in magnetised atmospheres remains poorly understood, magnetic field properties have to be considered in the framework of dynamo or fossil field theories. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The influence of Hall drift on the magnetic fields of neutron stars

We consider the evolution of magnetic fields under the influence of Hall drift and Ohmic decay. The governing equation is solved numerically, in a spherical shell with   r _i / r _o = 0.75  . Starting with simple free-decay modes as initial conditions, we then consider the subsequent evolution. The Hall effect induces so-called helicoidal oscillations, in which energy is redistributed among the different modes. We find that the amplitude of these oscillations can be quite substantial, with some of the higher harmonics becoming comparable with the original field. Nevertheless, this transfer of energy to the higher harmonics is not sufficient to accelerate significantly the decay of the original field, at least not at the   R _B = O (100)  parameter values accessible to us, where this Hall parameter   R _B   measures the ratio of the Ohmic time-scale to the Hall time-scale. We do find clear evidence though of increasingly fine structures developing for increasingly large   R _B   , suggesting that perhaps this Hall-induced cascade to ever-shorter length-scales is eventually sufficiently vigorous to enhance the decay of the original field. Finally, the implications for the evolution of neutron star magnetic fields are discussed.     

Large-scale magnetic field of the G8 dwarf ξ Bootis A

We investigate the magnetic geometry of the active G8 dwarf ξ Bootis A (ξ Boo A), from spectropolarimetric observations obtained in 2003 with the MuSiCoS échelle spectropolarimeter at the Télescope Bernard Lyot (Observatoire du Pic du Midi, France). We repeatedly detect a photospheric magnetic field, with periodic variations consistent with rotational modulation. Circularly polarized (Stokes V) line profiles present a systematic asymmetry, showing up as an excess in amplitude and area of the blue lobe of the profiles. Direct modelling of Stokes V profiles suggests that the global magnetic field is composed of two main components, with an inclined dipole and a large-scale toroidal field. We derive a dipole intensity of about 40 G, with an inclination of 35° of the dipole with respect to the rotation axis. The toroidal field strength is of the order of 120 G. A noticeable evolution of the field geometry is observed over the 40 nights of our observation window and results in an increase in field strength and dipole inclination.
This study is the first step of a long-term monitoring of ξ Boo A and other active solar-type stars, with the aim of investigating secular fluctuations of stellar magnetic geometries induced by activity cycles.