Dynamo action in simulations of penetrative solar convection with an imposed tachocline

We summarize new and continuing three-dimensional spherical shell simulations of dynamo action by convection allowed to penetrate downward into a tachocline of rotational shear. The inclusion of an imposed tachocline allows us to examine several processes believed to be essential in the operation of the global solar dynamo, including differential rotation, magnetic pumping, and the stretching and organization of fields within the tachocline. In the stably stratified core, our simulations reveal that strong axisymmetric magnetic fields (of ∼ 3000 G strength) can be built, and that those fields generally exhibit a striking antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. In the convection zone above, fluctuating fields dominate over weaker mean fields. New calculations indicate that the tendency toward toroidal fields of antisymmetric parity is relatively insensitive to initial magnetic field configurations; they also reveal that on decade-long timescales, the magnetic fields can briefly enter (and subsequently emerge from) states of symmetric parity.We have not yet observed any overall reversals of the field polarity, nor systematic latitudinal propagation. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamo action in an ambient field

The origin of global magnetic fields in celestial bodies is generally ascribed to dynamo action by fluid motions in their electrically conducting interiors. Some objects – e.g. close‐in extra‐solar planets or the moons of some giant planets – are embedded in ambient magnetic fields which modify the generation of the internal field in these bodies. Recently, the feedback of the magnetospheric field by Chapman‐Ferraro currents in the magnetopause onto the interior dynamo has been proposed to explain the observed weakness of the intrinsic magnetic field of planet Mercury. We study a simplified mean‐field dynamo model which allows us to analytically address various issues like positive and negative feedback situations, stationary versus time‐dependent solutions, and the stability of weak and strong field branches. We discuss the influence of the response function on the solutions when the external field depends on the strength of the intrinsic field like in the situation of the feedback dynamo of Mercury. We find that the feedback mechanism works only for a narrow range of dynamo numbers in the case of Mercury which makes him unique in our solar system. We conclude with some implications for extra‐solar planets (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

High-precision magnetic field measurements of Ap and Bp stars

In this paper we describe a new approach for measuring the mean longitudinal magnetic field and net linear polarization of Ap and Bp stars. As was demonstrated by Wade et al., least-squares deconvolution (LSD; Donati et al.) provides a powerful technique for detecting weak Stokes V , Q and U Zeeman signatures in stellar spectral lines. These signatures have the potential to apply strong new constraints to models of stellar magnetic field structure. Here we point out two important uses of LSD Stokes profiles. First, they can provide very precise determinations of the mean longitudinal magnetic field. In particular, this method allows one frequently to obtain 1 σ error bars better than 50 G, and smaller than 20 G in some cases. This method is applicable to both broad- and sharp-lined stars, with both weak and strong magnetic fields, and effectively redefines the quality standard of longitudinal field determinations. Secondly, LSD profiles can in some cases provide a measure of the net linear polarization, a quantity analogous to the broad-band linear polarization recently used to derive detailed magnetic field models for a few stars (e.g. Leroy et al.). In this paper we report new high-precision measurements of the longitudinal fields of 14 magnetic Ap/Bp stars, as well as net linear polarization measurements for four of these stars, derived from LSD profiles.     

On non-axisymmetric magnetic equilibria in stars

In previous work, stable approximately axisymmetric equilibrium configurations for magnetic stars were found by numerical simulation. Here, I investigate the conditions under which more complex, non-axisymmetric configurations can form. I present numerical simulations of the formation of stable equilibria from turbulent initial conditions and demonstrate the existence of non-axisymmetric equilibria consisting of twisted flux tubes lying horizontally below the surface of the star, meandering around the star in random patterns. Whether such a non-axisymmetric equilibrium or a simple axisymmetric equilibrium forms depends on the radial profile of the strength of the initial magnetic field. The results could explain observations of non-dipolar fields on stars such as the B0.2 main-sequence star τ Sco or the pulsar 1E 1207.4-5209. The secular evolution of these equilibria due to Ohmic and buoyancy processes is also examined.     

Stability of toroidal magnetic fields in the solar tachocline and beneath

Stability of toroidal magnetic field in a stellar radiation zone is considered for the cases of uniform and differential rotation. In the rigidly rotating radiative core shortly below the tachocline, the critical magnetic field for instability is about 600 G. The unstable disturbances for slightly supercritical fields have short radial scales ∼1 Mm. Radial mixing produced by the instability is estimated to conclude that the internal field of the sun can exceed the critical value of 600 G only marginally. Otherwise, the mixing is too strong and not compatible with the observed lithium abundance. Analysis of joint instability of differential rotation and toroidal field leads to the conclusion that axisymmetric models of the laminar solar tachocline are stable to nonaxisymmetric disturbances. The question of whether sun-like stars can posses tachoclines is addressed with positive answer for stars with rotation periods shorter than about two months. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamo action and the period gap in cataclysmic variables

The conjecture is presented that the gap in the distribution of the orbital periods of cataclysmic variables is related to a particular kind of hydromagnetic dynamo, called an interface dynamo, operating near the base of the convective envelope of their secondary components. Such a dynamo is characterized by the spatial separation of the regions where differential rotation and the α effect operate. Unlike conventional dynamos, the linear growth rate of an interface dynamo becomes negative for highly supercritical dynamo numbers, leading to the disappearance of the dynamo action. If such a result, from linear theory, is confirmed by non-linear calculations, it may provide a physical basis for the so-called disrupted magnetic braking hypothesis, invoked to explain the existence of the period gap by several evolutionary models of cataclysmic variables.     

Secular variability of the longitudinal magnetic field of the Ap star γ Equ

We present an analysis of the secular variability of the longitudinal magnetic field B _e in the roAp star γ Equ (HD 201601). Measurements of the stellar magnetic field B _e were mostly compiled from the literature, and we appended also our 33 new B _e measurements which were obtained with the 1-m optical telescope of the Special Astrophysical Observatory (Russia). All the available data cover the time period of 58 yr, and include both phases of the maximum and minimum B _e. We determined that the period of the long-term magnetic B _e variations equals  91.1 ± 3.6 yr  , with   B _e(max) =+577 ± 31 G  and   B _e(min) =−1101 ± 31 G  .     

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.     

The magnetic field and rotational period of 53 Camelopardalis

The magnetic Ap star 53 Cam has been described as the best object for testing theories of elemental diffusion in the presence of a stellar magnetic field. It is thus vital to determine the rotational period of this very important star as accurately as possible. An improved period (8.026 81±0.000 04 d) is derived here by combining new measurements of the effective magnetic field with previously published values.     

How can α 2-dynamos generate axisymmetric magnetic fields?

Inspired by a recently observed axisymmetric field in a fully convective star we investigate the influence of an anisotropic diffusivity on the dynamo. We find that with reasonable assumptions for the anisotropy of the diffusivity and the α -effect the preference of axisymmetric modes is achieved. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Chemical spots and oscillatory diffusion modes in magnetic stars

The stars of the middle main sequence often have spot‐like chemical structures at their surfaces. We consider diffusion caused by electric currents and argue that such current‐driven diffusion can form chemical inhomogeneities in a plasma. The considered mechanism can contribute to a formation of element spots in Hg‐Mn and Ap‐stars. Due to the Hall effect, diffusion in the presence of electric currents can be accompanied by the propagation of a particular type of magnetohydrodynamic modes in which only the impurity number density oscillates. Such modes exist if the magnetic pressure is much greater than the gas pressure and can be the reason for variations of the abundance peculiarities in stars. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamo coefficients from local simulations of the turbulent ISM

Observations in polarized emission reveal the existence of large‐scale coherent magnetic fields in a wide range of spiral galaxies. Radio‐polarization data show that these fields are strongly inclined towards the radial direction, with pitch angles up to 35° and thus cannot be explained by differential rotation alone. Global dynamo models describe the generation of the radial magnetic field from the underlying turbulence via the so called α ‐effect. However, these global models still rely on crude assumptions about the small‐scale turbulence. To overcome these restrictions we perform fully dynamical MHD simulations of interstellar turbulence driven by supernova explosions. From our simulations we extract profiles of the contributing diagonal elements of the dynamo α ‐tensor as functions of galactic height. We also measure the coefficients describing vertical pumping and find that the ratio between these two effects has been overestimated in earlier analytical work, where dynamo action seemed impossible. In contradiction to these models based on isolated remnants we always find the pumping to be directed inward. In addition we observe that depends on whether clustering in terms of superbubbles is taken into account. Finally, we apply a test field method to derive a quantitative measure of the turbulent magnetic diffusivity which we determine to be ∼2 kpckms^–1. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dynamo saturation in direct simulations of the multi‐phase turbulent interstellar medium

The ordered magnetic field observed via polarised synchrotron emission in nearby disc galaxies can be explained by a mean‐field dynamo operating in the diffuse interstellar medium (ISM). Additionally, vertical‐flux initial conditions are potentially able to influence this dynamo via the occurrence of the magnetorotational instability (MRI). We aim to study the influence of various initial field configurations on the saturated state of the mean‐field dynamo. This is motivated by the observation that different saturation behaviour was previously obtained for different supernova rates. We perform direct numerical simulations (DNS) of three‐dimensional local boxes of the vertically stratified, turbulent interstellar medium, employing shearing‐periodic boundary conditions horizontally. Unlike in our previous work, we also impose a vertical seed magnetic field. We run the simulations until the growth of the magnetic energy becomes negligible. We furthermore perform simulations of equivalent 1D dynamo models, with an algebraic quenching mechanism for the dynamo coefficients. We compare the saturation of the magnetic field in the DNS with the algebraic quenching of a mean‐field dynamo. The final magnetic field strength found in the direct simulation is in excellent agreement with a quenched α) dynamo. For supernova rates representative of the Milky Way, field losses via a Galactic wind are likely responsible for saturation. We conclude that the relative strength of the turbulent and regular magnetic fields in spiral galaxies may depend on the galaxy's star formation rate. We propose that a mean field approach with algebraic quenching may serve as a simple sub‐grid scale model for galaxy evolution simulations including a prescribed feedback from magnetic fields. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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)     

A coupled model of magnetic flux generation and transport in stars

We present a combined model for magnetic field generation and transport in cool stars with outer convection zones. The mean toroidal magnetic field, which is generated by a cyclic thin-layer α Ω dynamo at the bottom of the convection zone is taken to determine the emergence probability of magnetic flux tubes in the photosphere. Following the nonlinear rise of the unstable thin flux tubes, emergence latitudes and tilt angles of bipolar magnetic regions are determined. These quantities are put into a surface flux transport model, which simulates the surface evolution of magnetic flux under the effects of large-scale flows and turbulent diffusion. First results are discussed for the case of the Sun and for more rapidly rotating solar-type stars. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Alleviation of catastrophic quenching in solar dynamo model with nonlocal alpha‐effect

The nonlocal alpha‐effect of Babcock‐Leighton type is not prone to the catastrophic quenching due to conservation of magnetic helicity. This is shown with a dynamo model, which jointly applies the nonlocal alpha‐effect, the diamagnetic pumping, and dynamical equation for the magnetic alpha‐effect. The same model shows catastrophic quenching when the alpha‐effect is changed to its local formulation. The nonlocal model shows a preferred excitation of magnetic fields of dipolar symmetry, which oscillate with a period of about ten years and have a toroidal‐to‐polar fields ratio of about a thousand (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Probability distribution functions for the Sun's magnetic field

Magnetoconvection structures the Sun's magnetic field cover a vast range of scales, down to the magnetic diffusion scale that is orders of magnitude smaller than the resolution of current telescopes. The statistical properties of this structuring are governed by probability density functions (PDFs) for the flux densities and by the angular distribution functions for the orientations of the field vectors. The magnetic structuring on sub‐pixel scales greatly affects the field properties averaged over the resolution element due to the non‐linear relation between polarization and magnetic field. Here we use a Hinode SOT/SP data set for the quiet Sun disk center to explore the complex behavior of the 6301–6302 Å Stokes line profile system and identify the observables that allow us to determine the distribution functions in the most robust and least model dependent way. The angular distribution is found to be strongly peaked around the vertical direction for large flux densities but widens with decreasing flux density to become isotropic in the limit of zero flux density. The noise‐corrected PDFs for the vertical, horizontal, and total flux densities all have a narrowly peaked maximum at zero flux density that can be fitted with a stretched exponential, while the extended wings decline quadratically (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)