Source author record

Axel Brandenburg

Axel Brandenburg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

154works

26topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Big Bang Nucleosynthesis Limits and Relic Gravitational Waves Detection Prospects

We revisit the big bang nucleosynthesis (BBN) limits on primordial magnetic fields and/or turbulent motions accounting for the decaying nature of turbulent sources between the time of generation and BBN. This leads to larger estimates for the gravitational wave (GW) signal than previously expected. We address the detection prospects through space-based interferometers (for GWs generated around the electroweak energy scale) as well as pulsar timing arrays and astrometric missions (for GWs generated around the quantum chromodynamics energy scale).

preprint2022arXiv

Chirality in Astrophysics

Chirality, or handedness, enters astrophysics in three distinct ways. Magnetic field and vortex lines tend to be helical and have a systematic twist in the northern and southern hemispheres of a star or a galaxy. Helicity is here driven by external factors. Chirality can also enter at the microphysical level and can then be traced back to the parity-breaking weak force. Finally, chirality can arise spontaneously, but this requires not only the presence of an instability, but also the action of nonlinearity. Examples can be found both in magnetohydrodynamics and in astrobiology, where homochirality among biomolecules probably got established at the origin of life. In this review, all three types of chirality production will be explored and compared.

preprint2022arXiv

Collision fluctuations of lucky droplets with superdroplets

It was previously shown that the superdroplet algorithm for modeling the collision-coalescence process can faithfully represent mean droplet growth in turbulent clouds. But an open question is how accurately the superdroplet algorithm accounts for fluctuations in the collisional aggregation process. Such fluctuations are particularly important in dilute suspensions. Even in the absence of turbulence, Poisson fluctuations of collision times in dilute suspensions may result in substantial variations in the growth process, resulting in a broad distribution of growth times to reach a certain droplet size. We quantify the accuracy of the superdroplet algorithm in describing the fluctuating growth history of a larger droplet that settles under the effect of gravity in a quiescent fluid and collides with a dilute suspension of smaller droplets that were initially randomly distributed in space ('lucky droplet model'). We assess the effect of fluctuations upon the growth history of the lucky droplet and compute the distribution of cumulative collision times. The latter is shown to be sensitive enough to detect the subtle increase of fluctuations associated with collisions between multiple lucky droplets. The superdroplet algorithm incorporates fluctuations in two distinct ways: through the random spatial distribution of superdroplets and through the Monte Carlo collision algorithm involved. Using specifically designed numerical experiments, we show that both on their own give an accurate representation of fluctuations. We conclude that the superdroplet algorithm can faithfully represent fluctuations in the coagulation of droplets driven by gravity.

preprint2022arXiv

Compressible test-field method and its application to shear dynamos

In this study we present a compressible test-field method (CTFM) for computing $α$ effect and turbulent magnetic diffusivity tensors, as well as those relevant for mean ponderomotive force and mass source, applied to the full MHD equations. We describe the theoretical background of the method, and compare it to the quasi-kinematic test-field method, and to the previously studied variant working in simplified MHD (SMHD). We present several test cases using velocity and magnetic fields of the Roberts geometry, and also compare with the imposed-field method. We show that, for moderate imposed field strengths, the nonlinear CTFM (nCTFM) gives results in agreement with the imposed-field method. Comparison of different flavors of the nCTFM in the shear dynamo case also agree up to equipartition field strengths. Some deviations between the CTFM and SMHD variants exist. As a relevant physical application, we study non-helically forced shear flows, which exhibit large-scale dynamo action, and present a re-analysis of low Reynolds number, moderate shear systems, where we previously neglected the pressure gradient in the momentum equation, and found no coherent shear-current effect. Another key difference is that in the earlier study we used magnetic forcing to mimic small-scale dynamo action, while here it is self-consistently driven by purely kinetic forcing. The kinematic CTFM with general validity forms the core of our analysis. We still find no coherent shear-current effect, but do recover strong large-scale dynamo action that, according to our analysis, is driven through the incoherent effects.

preprint2022arXiv

Dynamo effect in unstirred self-gravitating turbulence

In many astrophysical environments, self-gravity can generate kinetic energy, which, in principle, is available for driving dynamo action. Using direct numerical simulations, we show that in unstirred self-gravitating subsonic turbulence with helicity and a magnetic Prandtl number of unity, there is a critical magnetic Reynolds number of about 25 above which the work done against the Lorentz force exceeds the Ohmic dissipation. The collapse itself drives predominantly irrotational motions that cannot be responsible for dynamo action. We find that, with a weak magnetic field, one-third of the work done by the gravitational force goes into compressional heating and the remaining two-thirds go first into kinetic energy of the turbulence before a fraction of it is converted further into magnetic and finally thermal energies. Close to the collapse, however, these fractions change toward 1/4 and 3/4 for compressional heating and kinetic energy, respectively. When the magnetic field is strong, the compressional heating fraction is unchanged. Out of the remaining kinetic energy, one quarter goes directly into magnetic energy via work against the Lorentz force. The fraction of vortical motions diminishes in favor of compressive motions that are almost exclusively driven by the Jeans instability. For an initially uniform magnetic field, field amplification at scales larger than those of the initial turbulence are driven by tangling.

preprint2022arXiv

Dynamo instabilities in plasmas with inhomogeneous chiral chemical potential

We study the dynamics of magnetic fields in chiral magnetohydrodynamics, which takes into account the effects of an additional electric current related to the chiral magnetic effect in high energy plasmas. We perform direct numerical simulations, considering weak seed magnetic fields and inhomogeneities of the chiral chemical potential mu_5 with a zero mean. We demonstrate that a small-scale chiral dynamo can occur in such plasmas if fluctuations of mu_5 are correlated on length scales that are much larger than the scale on which the dynamo growth rate reaches its maximum. Magnetic fluctuations grow by many orders of magnitude due to the small-scale chiral dynamo instability. Once the nonlinear backreaction of the generated magnetic field on fluctuations of mu_5 sets in, the ratio of these scales decreases and the dynamo saturates. When magnetic fluctuations grow sufficiently to drive turbulence via the Lorentz force before reaching maximum field strength, an additional mean-field dynamo phase is identified. The mean magnetic field grows on a scale that is larger than the integral scale of turbulence after the amplification of the fluctuating component saturates. The growth rate of the mean magnetic field is caused by a magnetic alpha effect that is proportional to the current helicity. With the onset of turbulence, the power spectrum of mu_5 develops a universal k^(-1) scaling independently of its initial shape, while the magnetic energy spectrum approaches a k^(-3) scaling.

preprint2022arXiv

Evolution of primordial magnetic fields during large-scale structure formation

Primordial magnetic fields could explain the large-scale magnetic fields present in the Universe. Inflation and phase transitions in the early Universe could give rise to such fields with unique characteristics. We investigate the magneto-hydrodynamic evolution of these magnetogenesis scenarios with cosmological simulations. We evolve inflation-generated magnetic fields either as (i) uniform (homogeneous) or as (ii) scale-invariant stochastic fields, and phase transition-generated ones either as (iii) helical or as (iv) non-helical fields from the radiation-dominated epoch. We find that the final distribution of magnetic fields in the simulated cosmic web shows a dependence on the initial strength and the topology of the seed field. Thus, the observed field configuration retains information on the initial conditions at the moment of the field generation. If detected, primordial magnetic field observations would open a new window for indirect probes of the early universe. The differences between the competing models are revealed on the scale of galaxy clusters, bridges, as well as filaments and voids. The distinctive spectral evolution of different seed fields produces imprints on the correlation length today. We discuss how the differences between rotation measures from highly ionized regions can potentially be probed with forthcoming surveys.

preprint2022arXiv

Polarization of gravitational waves from helical MHD turbulent sources

We use direct numerical simulations of decaying primordial hydromagnetic turbulence with helicity to compute the resulting gravitational wave (GW) production and its degree of circular polarization. We find a clear dependence of the polarization of the resulting GWs on the fractional helicity of the turbulent source and we show that driven magnetic fields produce GWs more efficiently than magnetic fields that are initially present, leading to larger spectral amplitudes. The helicity does not have a huge impact on the maximum spectral amplitude in any of the two types of turbulence considered. However, the GW spectrum at wave numbers away from the peak becomes smaller for larger values of the magnetic fractional helicity. The degree of circular polarization approaches zero at frequencies below the peak, and reaches its maximum at the peak. At higher frequencies, it stays finite if the magnetic field is initially present, and it approaches zero if it is driven. We predict that the spectral peak of the GW signal can be detected by LISA if the turbulent energy density is at least $\sim 3\%$ of the radiation energy density, and the characteristic scale is a hundredth of the horizon at the electroweak scale. We show that the resulting GW polarization is unlikely to be detectable by the anisotropies induced by our proper motion in the dipole response function of LISA. Such signals can, however, be detectable by cross-correlating data from the LISA-Taiji network for turbulent energy densities of $\sim 5\%$, and fractional helicity of 0.5 to 1. Second-generation space-based GW detectors, such as BBO and DECIGO, would allow for the detection of a larger range of the GW spectrum and smaller amplitudes of the magnetic field.

preprint2022arXiv

Production of a chiral magnetic anomaly with emerging turbulence and mean-field dynamo action

In relativistic magnetized plasmas, asymmetry in the number densities of left- and right-handed fermions, i.e., a non-zero chiral chemical potential mu_5, leads to an electric current along the magnetic field. This causes a chiral dynamo instability for a uniform mu_5, but our simulations reveal a dynamo even for fluctuating mu_5 with zero mean. It produces magnetically-dominated turbulence and generates mean magnetic fields via the magnetic alpha effect. Eventually, a universal scale-invariant k^(-1) spectrum of mu_5 and a k^(-3) magnetic spectrum are formed independently of the initial condition.

preprint2021arXiv

Circular Polarization of Gravitational Waves from Early-Universe Helical Turbulence

We perform direct numerical simulations to compute the net circular polarization of gravitational waves from helical (chiral) turbulent sources in the early Universe for a variety of initial conditions, including driven (stationary) and decaying turbulence. We investigate the resulting gravitational wave signal assuming different turbulent geneses such as magnetically or kinetically driven cases. Under realistic physical conditions in the early Universe we compute numerically the wave number-dependent polarization degree of the gravitational waves. We find that the spectral polarization degree strongly depends on the initial conditions. The peak of the spectral polarization degree occurs at twice the typical wavenumber of the source, as expected, and for fully helical decaying turbulence, it reaches its maximum of nearly 100\% {\it only} at the peak. We determine the temporal evolution of the turbulent sources as well as the resulting gravitational waves, showing that the dominant contribution to their spectral energy density happens shortly after the activation of the source. Only through an artificially prolonged decay of the turbulence can further increase of the gravitational wave amplitude be achieved. We estimate the detection prospects for the net polarization, arguing that its detection contains {\it clean} information (including the generation mechanisms, time, and strength) about the sources of possible parity violations in the early Universe.

preprint2021arXiv

On the existence of shear-current effects in magnetized burgulence

The possibility of explaining shear flow dynamos by a magnetic shear-current (MSC) effect is examined via numerical simulations. Our primary diagnostics is the determination of the turbulent magnetic diffusivity tensor $\boldsymbolη$. In our setup, a negative sign of its component $η_{yx}$ is necessary for coherent dynamo action by the SC effect. To be able to measure turbulent transport coefficients from systems with magnetic background turbulence, we present an extension of the test-field method (TFM), applicable to our setup where the pressure gradient is dropped from the momentum equation: the nonlinear TFM (NLTFM). Our momentum equation is related to Burgers' equation and the resulting flows are referred to as magnetized burgulence. We use both stochastic kinetic and magnetic forcings to mimic cases without and with simultaneous small-scale dynamo action (SSD). When we force only kinetically, negative $η_{yx}$ are obtained with exponential growth in both the radial and azimuthal mean magnetic field components. Using isotropic magnetic forcing, the field growth is no longer exponential, while NLTFM yields positive $η_{yx}$. By employing an alternative forcing from which wavevectors having small components are removed, the exponential growth is recovered, but the NLTFM results do not change significantly. Analyzing the dynamo excitation conditions for the coherent SC and incoherent $α$ and SC effects shows that the incoherent effects are the main drivers of the dynamo in the majority of cases. We find no evidence for MSC-effect-driven dynamos in our simulations.

preprint2020arXiv

Application of a helicity proxy to edge-on galaxies

We study the prospects of detecting magnetic helicity in galaxies by observing the dust polarization of the edge-on galaxy NGC 891. Our numerical results of mean-field dynamo calculations show that there should be a large-scale component of the rotationally invariant parity-odd B polarization that we predict to be negative in the first and third quadrants, and positive in the second and fourth quadrants. The large-scale parity-even E polarization is predicted to be negative near the axis and positive further away in the outskirts. These properties are shown to be mostly a consequence of the magnetic field being azimuthal and the polarized intensity being maximum at the center of the galaxy and are not a signature of magnetic helicity.

preprint2020arXiv

Condensational and collisional growth of cloud droplets in a turbulent environment

We investigate the effect of turbulence on the combined condensational and collisional growth of cloud droplets by means of high resolution direct numerical simulations of turbulence and a superparticle approximation for droplet dynamics and collisions. The droplets are subject to turbulence as well as gravity, and their collision and coalescence efficiencies are taken to be unity. We solve the thermodynamic equations governing temperature, water-vapor mixing ratio, and the resulting supersaturation fields together with the Navier-Stokes equation. We find that the droplet-size distribution broadens with increasing Reynolds number and/or mean energy dissipation rate. Turbulence affects the condensational growth directly through supersaturation fluctuations, and it influences collisional growth indirectly through condensation. Our simulations show for the first time that, in the absence of the mean updraft cooling, supersaturation fluctuation-induced broadening of droplet-size distributions enhances the collisional growth. This is contrary to classical (non-turbulent) condensational growth, which leads to a growing mean droplet size, but a narrower droplet-size distribution. Our findings, instead, show that condensational growth facilitates collisional growth by broadening the size distribution in the tails at an early stage of rain formation. With increasing Reynolds numbers, evaporation becomes stronger. This counteracts the broadening effect due to condensation at late stages of rain formation. Our conclusions are consistent with results of laboratory experiments and field observations, and show that supersaturation fluctuations are important for precipitation.

preprint2020arXiv

Dynamo in weakly collisional nonmagnetized plasmas impeded by Landau damping of magnetic fields

We perform fully kinetic simulations of flows known to produce dynamo in magnetohydrodynamics (MHD), considering scenarios with low Reynolds number and high magnetic Prandtl number, relevant for galaxy cluster scale fluctuation dynamos. We find that Landau damping on the electrons leads to a rapid decay of magnetic perturbations, impeding the dynamo. This collisionless damping process operates on spatial scales where electrons are nonmagnetized, reducing the range of scales where the magnetic field grows in high magnetic Prandtl number fluctuation dynamos. When electrons are not magnetized down to the resistive scale, the magnetic energy spectrum is expected to be limited by the scale corresponding to magnetic Landau damping or, if smaller, the electron gyroradius scale, instead of the resistive scale. In simulations we thus observe decaying magnetic fields where resistive MHD would predict a dynamo.

preprint2020arXiv

Hall cascade with fractional magnetic helicity in neutron star crusts

The ohmic decay of magnetic fields in the crusts of neutron stars is generally believed to be governed by Hall drift which leads to what is known as a Hall cascade. Here we show that helical and fractionally helical magnetic fields undergo strong inverse cascading like in magnetohydrodynamics (MHD), but the magnetic energy decays more slowly with time $t$: $\propto t^{-2/5}$ instead of $\propto t^{-2/3}$ in MHD. Even for a nonhelical magnetic field there is a certain degree of inverse cascading for sufficiently strong magnetic fields. The inertial range scaling with wavenumber $k$ is compatible with earlier findings for the forced Hall cascade, i.e., proportional to $k^{-7/3}$, but in the decaying cases, the subinertial range spectrum steepens to a novel $k^5$ slope instead of the $k^4$ slope in MHD. The energy of the large-scale magnetic field can increase quadratically in time through inverse cascading. For helical fields, the energy dissipation is found to be inversely proportional to the large-scale magnetic field and proportional to the fifth power of the root-mean square (rms) magnetic field. For neutron star conditions with an rms magnetic field of a few times $10^{14}\,$G, the large-scale magnetic field might only be $10^{11}\,$G, while still producing magnetic dissipation of $10^{33}\,$erg$\,$s$^{-1}$ for thousands of years, which could manifest itself through X-ray emission. Finally, it is shown that the conclusions from local unstratified models agree rather well with those from stratified models with boundaries.

preprint2020arXiv

Hemispheric handedness in the Galactic synchrotron polarization foreground

The large-scale magnetic field of the Milky Way is thought to be created by an alpha-Omega dynamo, which implies that it should have opposite handedness North and South of the Galactic midplane. Here we attempt to detect a variation in handedness using polarization data from the Wilkinson Microwave Anisotropy Probe. Previous analyzes of the parity-even and parity-odd parts of linear polarization of the global dust and synchrotron emission have focused on quadratic correlations in spectral space of, and between, these two components. Here, by contrast, we analyze the parity-odd polarization itself and show that it has, on average, opposite signs in Northern and Southern Galactic hemispheres. Comparison with a Galactic mean-field dynamo model shows broad qualitative agreement and reveals that the sign of the observed hemispheric dependence of the azimuthally averaged parity-odd polarization is not determined by the sign of alpha, but by the sense of differential rotation.

preprint2020arXiv

Magnetic helicity dissipation and production in an ideal MHD code

We study a turbulent helical dynamo in a periodic domain by solving the ideal magnetohydrodynamic (MHD) equations with the FLASH code using the divergence-cleaning eight-wave method and compare our results with direct numerical simulations (DNS) using the Pencil Code. At low resolution, FLASH reproduces the DNS results qualitatively by developing the large-scale magnetic field expected from DNS, but at higher resolution, no large-scale magnetic field is obtained. In all those cases in which a large-scale magnetic field is generated, the ideal MHD results yield too little power at small scales. As a consequence, the small-scale current helicity is too small compared with that of the DNS. The resulting net current helicity has then always the wrong sign, and its statistical average also does not approach zero at late times, as expected from the DNS. Our results have implications for astrophysical dynamo simulations of stellar and galactic magnetism using ideal MHD codes.

preprint2020arXiv

On the measurement of handedness in Fermi Large Area Telescope data

A handedness in the arrival directions of high-energy photons from outside our Galaxy can be related to the helicity of an intergalactic magnetic field. Previous estimates by arXiv:1310.4826 and arXiv:1412.3171 showed a hint of a signal present in the photons observed by the Fermi Large Area Telescope (LAT). An update on the measurement of handedness in Fermi-LAT data is presented using more than 10 years of observations. Simulations are performed to study the uncertainty of the measurements, taking into account the structure of the exposure caused by the energy-dependent instrument response and its observing profile, as well as the background from the interstellar medium. The simulations are required to accurately estimate the uncertainty and to show that previously the uncertainty was significantly underestimated. The apparent signal in the earlier analysis of {\em Fermi}-LAT data is rendered non-significant.

preprint2020arXiv

Piecewise quadratic growth during the 2019 novel coronavirus epidemic

The temporal growth in the number of deaths in the COVID-19 epidemic is subexponential. Here we show that a piecewise quadratic law provides an excellent fit during the thirty days after the first three fatalities on January 20 and later since the end of March 2020. There is also a brief intermediate period of exponential growth. During the second quadratic growth phase, the characteristic time of the growth is about eight times shorter than in the beginning, which can be understood as the occurrence of separate hotspots. Quadratic behavior can be motivated by peripheral growth when further spreading occurs only on the outskirts of an infected region. We also study numerical solutions of a simple epidemic model, where the spatial extend of the system is taken into account. To model the delayed onset outside China together with the early one in China within a single model with minimal assumptions, we adopt an initial condition of several hotspots, of which one reaches saturation much earlier than the others. At each site, quadratic growth commences when the local number of infections has reached a certain saturation level. The total number of deaths does then indeed follow a piecewise quadratic behavior.

preprint2020arXiv

Primordial magnetic helicity evolution with a homogeneous magnetic field from inflation

Motivated by a scenario of magnetogenesis in which a homogeneous magnetic field is generated during inflation, we study the magnetohydrodynamic evolution of the primordial plasma motions for two kinds of initial conditions -- (i) a spatially homogeneous field with an unlimited correlation length, and (ii) a zero flux scale-invariant statistically homogeneous magnetic field. In both cases, we apply, for a short initial time interval, monochromatic forcing at a certain wave number so that the correlation length is finite, but much smaller than the typical length scale of turbulence. In particular, we investigate the decay of nonhelical and helical hydromagnetic turbulence. We show that, in the presence of a homogeneous magnetic field, the decay of helical and nonhelical small-scale fields can occur rapidly. This is a special property of a system with a perfectly homogeneous magnetic field, which is sometimes considered as a local approximation to a slowly varying background field. It can never change and acts as an imposed magnetic field. This is in a sharp contrast to the case of a statistically homogeneous magnetic field, where we recover familiar decay properties: a much slower decay of magnetic energy and a faster growth of the correlation length, especially in the case with magnetic felicity. The result suggests that a homogeneous magnetic field, if generated during inflation, should persist under the influence of small-scale fields and could be the origin of the large-scale magnetic field in the Universe.

preprint2020arXiv

The turbulent stress spectrum in the inertial and subinertial ranges

For velocity and magnetic fields, the turbulent pressure and, more generally, the squared fields such as the components of the turbulent stress tensor, play important roles in astrophysics. For both one and three dimensions, we derive the equations relating the energy spectra of the fields to the spectra of their squares. We solve the resulting integrals numerically and show that for turbulent energy spectra of Kolmogorov type, the spectral slope of the stress spectrum is also of Kolmogorov type. For shallower turbulence spectra, the slope of the stress spectrum quickly approaches that of white noise, regardless of how blue the spectrum of the field is. For fully helical fields, the stress spectrum is elevated by about a factor of two in the subinertial range, while that in the inertial range remains unchanged. We discuss possible implications for understanding the spectrum of primordial gravitational waves from causally generated magnetic fields during cosmological phase transitions in the early universe. We also discuss potential diagnostic applications to the interstellar medium, where polarization and scintillation measurements characterize the square of the magnetic field.

preprint2020arXiv

Turbulent viscosity and effective magnetic Prandtl number from simulations of isotropically forced turbulence

(abridged) Context: Turbulent diffusion of large-scale flows and magnetic fields play major roles in many astrophysical systems. Aims: Our goal is to compute turbulent viscosity and magnetic diffusivity, relevant for diffusing large-scale flows and magnetic fields, respectively, and their ratio, the turbulent magnetic Prandtl number, ${\rm Pm}_{\rm t}$, for isotropically forced homogeneous turbulence. Methods: We use simulations of forced turbulence in fully periodic cubes composed of isothermal gas with an imposed large-scale sinusoidal shear flow. Turbulent viscosity is computed either from the resulting Reynolds stress or from the decay rate of the large-scale flow. Turbulent magnetic diffusivity is computed using the test-field method. The scale dependence of the coefficients is studied by varying the wavenumber of the imposed sinusoidal shear and test fields. Results: We find that turbulent viscosity and magnetic diffusivity are in general of the same order of magnitude. Furthermore, the turbulent viscosity depends on the fluid Reynolds number (${\rm Re}$) and scale separation ratio of turbulence. The scale dependence of the turbulent viscosity is found to be well approximated by a Lorentzian. The results for the turbulent transport coefficients appear to converge at sufficiently high values of ${\rm Re}$ and the scale separation ratio. However, a weak decreasing trend is found even at the largest values of ${\rm Re}$. The turbulent magnetic Prandtl number converges to a value that is slightly below unity for large ${\rm Re}$ whereas for small ${\rm Re}$, we find values between 0.5 and 0.6. Conclusions: The turbulent magnetic diffusivity is in general consistently higher than the turbulent viscosity. The actual value of ${\rm Pm}_{\rm t}$ found from the simulations ($\approx0.9\ldots0.95$) at large ${\rm Re}$ and scale separation ratio is higher than any of the analytic predictions.

preprint2019arXiv

The nature of mean-field generation in three classes of optimal dynamos

In recent years, several optimal dynamos have been discovered. They minimize the magnetic energy dissipation or, equivalently, maximize the growth rate at a fixed magnetic Reynolds number. In the optimal dynamo of Willis (2012, Phys. Rev. Lett. 109, 251101), we find mean-field dynamo action for planar averages. One component of the magnetic field grows exponentially while the other decays in an oscillatory fashion near onset. This behavior is different from that of an alpha^2 dynamo, where the two non-vanishing components of the planar averages are coupled and have the same growth rate. For the Willis dynamo, we find that the mean field is excited by a negative turbulent magnetic diffusivity, which has a non-uniform spatial profile near onset. The temporal oscillations in the decaying component are caused by the corresponding component of the diffusivity tensor being complex when the mean field is decaying and, in this way, time-dependent. The growing mean field can be modeled by a negative magnetic diffusivity combined with a positive magnetic hyperdiffusivity. In two other classes of optimal dynamos of Chen et al. (2015, J. Fluid Mech. 783, 23), we find, to some extent, similar mean-field dynamo actions. When the magnetic boundary conditions are mixed, the two components of the planar averaged field grow at different rates when the dynamo is 15% supercritical. When the mean magnetic field satisfies homogeneous boundary conditions (where the magnetic field is tangential to the boundary), mean-field dynamo action is found for one-dimensional averages, but not for planar averages. Despite having different spatial profiles, both dynamos show negative turbulent magnetic diffusivities. Our finding suggests negative turbulent magnetic diffusivities may support a broader class of dynamos than previously thought, including these three optimal dynamos.

preprint2018arXiv

Energetics of turbulence generated by chiral MHD dynamos

An asymmetry in the number density of left- and right-handed fermions is known to give rise to a new term in the induction equation that can result in a dynamo instability. At high temperatures, when a chiral asymmetry can survive for long enough, this chiral dynamo instability can amplify magnetic fields efficiently, which in turn drive turbulence via the Lorentz force. While it has been demonstrated in numerical simulations that this chiral magnetically driven turbulence exists and strongly affects the dynamics of the magnetic field, the details of this process remain unclear. The goal of this paper is to analyse the energetics of chiral magnetically driven turbulence and its effect on the generation and dynamics of the magnetic field using direct numerical simulations. We study these effects for different initial conditions, including a variation of the initial chiral chemical potential and the magnetic Prandtl number, Pm. In particular, we determine the ratio of kinetic to magnetic energy, Upsilon, in chiral magnetically driven turbulence. Within the parameter space explored in this study, Upsilon reaches a value of approximately 0.064-0.074 -independently of the initial chiral asymmetry and for Pm=1. Our simulations suggest, that Upsilon decreases as a power law when increasing Pm by decreasing the viscosity. While the exact scaling depends on the details of the fitting criteria and the Reynolds number regime, an approximate result of Upsilon(Pm)=0.1 Pm^{-0.4} is reported. Using the findings from our numerical simulations, we analyse the energetics of chiral magnetically driven turbulence in the early Universe.

preprint2018arXiv

Magnetism in the Early Universe

Blazar observations point toward the possible presence of magnetic fields over intergalactic scales of the order of up to $\sim1\,$Mpc, with strengths of at least $\sim10^{-16}\,$G. Understanding the origin of these large-scale magnetic fields is a challenge for modern astrophysics. Here we discuss the cosmological scenario, focussing on the following questions: (i) How and when was this magnetic field generated? (ii) How does it evolve during the expansion of the universe? (iii) Are the amplitude and statistical properties of this field such that they can explain the strengths and correlation lengths of observed magnetic fields? We also discuss the possibility of observing primordial turbulence through direct detection of stochastic gravitational waves in the mHz range accessible to LISA.

preprint2016arXiv

A unified large/small-scale dynamo in helical turbulence

We use high resolution direct numerical simulations (DNS) to show that helical turbulence can generate significant large-scale fields even in the presence of strong small-scale dynamo action. During the kinematic stage, the unified large/small-scale dynamo grows fields with a shape-invariant eigenfunction, with most power peaked at small scales or large $k$, as in \citet{SB14}. Nevertheless, the large-scale field can be clearly detected as an excess power at small $k$ in the negatively polarized component of the energy spectrum for a forcing with positively polarized waves. Its strength $bar{B}$, relative to the total rms field $B_{rms}$, decreases with increasing magnetic Reynolds number, $R_m$. However, as the Lorentz force becomes important, the field generated by the unified dynamo orders itself by saturating on successively larger scales. The magnetic integral scale for the positively polarized waves, characterizing the small-scale field, increases significantly from the kinematic stage to saturation. This implies that the small-scale field becomes as coherent as possible for a given forcing scale, which averts the $R_m$-dependent quenching of $\bar{B}/B_rms$. These results are obtained for $1024^3$ DNS with magnetic Prandtl numbers of $P_m=0.1$ and $10$. For $P_m=0.1$, $\bar{B}/B_{rms}$ grows from about $0.04$ to about $0.4$ at saturation, aided in the final stages by helicity dissipation. For $P_m=10$, $\bar{B}/B_{rms}$ grows from much less than 0.01 to values of the order the $0.2$. Our results confirm that there is a unified large/small-scale dynamo in helical turbulence.

preprint2016arXiv

Bipolar region formation in stratified two-layer turbulence

This work presents an extensive study of the previously discovered formation of bipolar flux concentrations in a two-layer model. We interpret the formation process in terms of negative effective magnetic pressure instability (NEMPI), which is a possible mechanism to explain the origin of sunspots. In our simulations, we use a Cartesian domain of isothermal stratified gas that is divided into two layers. In the lower layer, turbulence is forced with transverse nonhelical random waves, whereas in the upper layer no flow is induced. A weak uniform magnetic field is imposed in the entire domain at all times. In this study we vary the stratification by changing the gravitational acceleration, magnetic Reynolds number, strength of the imposed magnetic field, and size of the domain to investigate their influence on the formation process. Bipolar magnetic structure formation takes place over a large range of parameters. The magnetic structures become more intense for higher stratification until the density contrast becomes around $100$ across the turbulent layer. For the Reynolds numbers considered, magnetic flux concentrations are generated at magnetic Prandtl number between 0.1 and 1. The magnetic field in bipolar regions increases with higher imposed field strength until the field becomes comparable to the equipartition field strength of the turbulence. A larger horizontal extent enables the flux concentrations to become stronger and more coherent. The size of the bipolar structures turns out to be independent of the domain size. In the case of bipolar region formation, we find an exponential growth of the large-scale magnetic field, which is indicative of a hydromagnetic instability. Additionally, the flux concentrations are correlated with strong large-scale downward and converging flows. These findings imply that NEMPI is responsible for magnetic flux concentrations.

preprint2016arXiv

Evolution of Magnetic Helicity and Energy Spectra of Solar Active Regions

We adopt an isotropic representation of the Fourier-transformed two-point correlation tensor of the magnetic field to estimate the magnetic energy and helicity spectra as well as current helicity spectra of two individual active regions (NOAA 11158 and NOAA 11515) and the change of the spectral indices during their development as well as during the solar cycle. The departure of the spectral indices of magnetic energy and current helicity from 5/3 are analyzed, and it is found that it is lower than the spectral index of the magnetic energy spectrum. Furthermore, the fractional magnetic helicity tends to increase when the scale of the energy-carrying magnetic structures increases. The magnetic helicity of NOAA 11515 violates the expected hemispheric sign rule, which is interpreted as an effect of enhanced field strengths at scales larger than 30-60Mm with opposite signs of helicity. This is consistent with the general cycle dependence, which shows that around the solar maximum the magnetic energy and helicity spectra are steeper, emphasizing the large-scale field.

preprint2016arXiv

High-wavenumber solar $f$-mode strengthening prior to active region formation

We report a systematic strengthening of the local solar surface or fundamental $f$-mode $1$-$2$ days prior to the emergence of an active region (AR) in the same (corotating) location. Except for a possibly related increase in the kurtosis of the magnetic field, no indication can be seen in the magnetograms at that time. Our study is motivated by earlier numerical findings of Singh et al. (2014) which showed that, in the presence of a nonuniform magnetic field that is concentrated a few scale heights below the surface, the $f$-mode fans out in the diagnostic $kω$ diagram at high wavenumbers. Here we explore this possibility using data from the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory and show for six isolated ARs, 11130, 11158, 11242, 11105, 11072, and 11768, that at large latitudinal wavenumbers (corresponding to horizontal scales of around 3000 km), the $f$-mode displays strengthening about two days prior to AR formation and thus provides a new precursor for AR formation. Furthermore, we study two ARs, 12051 and 11678, apart from a magnetically quiet patch lying next to AR~12529, to demonstrate the challenges in extracting such a precursor signal when a newly forming AR emerges in a patch that lies in close proximity of one or several already existing ARs which are expected to pollute neighboring patches. We then discuss plausible procedures for extracting precursor signals from regions with crowded environments. The idea that the $f$-mode is perturbed days before any visible magnetic activity occurs at the surface can be important in constraining dynamo models aimed at understanding the global magnetic activity of the Sun.

preprint2016arXiv

Influence of a coronal envelope as a free boundary to global convective dynamo simulations

We explore the effects of an outer stably stratified coronal envelope on rotating turbulent convection, differential rotation, and large-scale dynamo action in spherical wedge models of the Sun. We solve the compressible magnetohydrodynamic equations in a two-layer model with unstable stratification below the surface, representing the convection zone, and a stably stratified coronal envelope above. The interface represents a free surface. We compare our model to models that have no coronal envelope. The presence of a coronal envelope is found to modify the Reynolds stress and the $Λ$ effect resulting in a weaker and non-cylindrical differential rotation. This is related to the reduced latitudinal temperature variations that are caused by and dependent on the angular velocity. Some simulations develop a near-surface shear layer that we can relate to a sign change in the meridional Reynolds stress term in the thermal wind balance equation. Furthermore, the presence of a free surface changes the magnetic field evolution since the toroidal field is concentrated closer to the surface. In all simulations, however, the migration direction of the mean magnetic field can be explained by the Parker--Yoshimura rule, which is consistent with earlier findings. A realistic treatment of the upper boundary in spherical dynamo simulations is crucial for the dynamics of the flow and magnetic field evolution.

preprint2016arXiv

Large-scale flow generation by inhomogeneous helicity

The effect of kinetic helicity (velocity--vorticity correlation) on turbulent momentum transport is investigated. The turbulent kinetic helicity (pseudoscalar) enters the Reynolds stress (mirrorsymmetric tensor) expression in the form of a helicity gradient as the coupling coefficient for the mean vorticity and/or the angular velocity (axial vector), which suggests the possibility of mean-flow generation in the presence of inhomogeneous helicity. This inhomogeneous helicity effect, which was previously confirmed at the level of a turbulence- or closure-model simulation, is examined with the aid of direct numerical simulations of rotating turbulence with non-uniform helicity sustained by an external forcing. The numerical simulations show that the spatial distribution of the Reynolds stress is in agreement with the helicity-related term coupled with the angular velocity, and that a large-scale flow is generated in the direction of angular velocity. Such a large-scale flow is not induced in the case of homogeneous turbulent helicity. This result confirms the validity of the inhomogeneous helicity effect in large-scale flow generation and suggests that a vortex dynamo is possible even in incompressible turbulence where there is no baroclinicity effect.

preprint2016arXiv

Magnetic concentrations in stratified turbulence: the negative effective magnetic pressure instability

In the presence of strong density stratification, hydromagnetic turbulence attains qualitatively new properties: the formation of magnetic flux concentrations. We review here the theoretical foundations of this mechanism in terms of what is now called the negative effective magnetic pressure instability. We also present direct numerical simulations of forced turbulence in strongly stratified layers and discuss the qualitative and quantitative similarities with corresponding mean-field simulations. Finally, the relevance to sunspot formation is discussed.

preprint2016arXiv

Multiple dynamo modes as a mechanism for long-term solar activity variations

Solar magnetic activity shows both smooth secular changes, such as the Grand Modern Maximum, and quite abrupt drops that are denoted as Grand Minima. Direct numerical simulations (DNS) of convection driven dynamos offer one way of examining the mechanisms behind these events. In this work, we analyze a solution of a solar-like DNS that has been evolved for roughly 80 magnetic cycles of 4.9 years, during which epochs of irregular behavior are detected. The emphasis of our analysis is to find physical causes for such behavior. The DNS employed is a semi-global (wedge) magnetoconvection model. For data analysis we use Ensemble Empirical Mode Decomposition (EEMD) and phase dispersion ($D^2$) methods. A special property of the DNS is the existence of multiple dynamo modes at different depths and latitudes. The dominant mode is solar-like. This mode is accompanied by a higher frequency mode near the surface and a low-frequency mode in the bottom of the convection zone. The overall behavior of the dynamo solution is very complex exhibiting variable cycle lengths, epochs of disturbed and even ceased surface activity, and strong short-term hemispherical asymmetries. Surprisingly, the most prominent suppressed surface activity epoch is actually a global magnetic energy maximum. We interpret the overall irregular behavior to be due to the interplay of the different dynamo modes showing different equatorial symmetries, especially the smoother part of the irregular variations being related to the variations of the mode strengths, evolving with different and variable cycle lengths. The abrupt low activity epoch in the dominant dynamo mode near the surface is related to a strong maximum of the bottom toroidal field strength, which causes abrupt disturbances especially in the differential rotation profile via the suppression of the Reynolds stresses.

preprint2016arXiv

Stellar mixing length theory with entropy rain

The effects of a non-gradient flux term originating from the motion of convective elements with entropy perturbations of either sign are investigated and incorporated into a modified version of stellar mixing length theory (MLT). Such a term, first studied by Deardorff in the meteorological context, might represent the effects of cold intense downdrafts caused by the rapid cooling in the granulation layer at the top of the convection zone of late-type stars. Such intense downdrafts were first seen in the strongly stratified simulations of Stein & Nordlund in the late 1980s. These downdrafts transport heat nonlocally, a phenomenon referred to as entropy rain. Moreover, the Deardorff term can cause upward enthalpy transport even in a weakly Schwarzschild-stably stratified layer. In that case, no giant cell convection would be excited. This is interest in view of recent observations, which could be explained if the dominant flow structures were of small scale even at larger depths. To study this possibility, three distinct flow structures are examined: one in which convective structures have similar size and mutual separation at all depths, one in which the separation increases with depth, but their size is still unchanged, and one in which both size and separation increase with depth, which is the standard flow structure. It is concluded that the third possibility with fewer and thicker downdrafts in deeper layers remains most plausible, but it may be unable to explain the suspected absence of large-scale flows with speeds and scales expected from MLT.

preprint2016arXiv

The evolution of primordial magnetic field since its generation

We study the evolution of primordial magnetic fields in an expanding cosmic plasma. For this purpose we present a comprehensive theoretical model to consider the evolution of MHD turbulence that can be used over a wide range of physical conditions, including cosmological and astrophysical applications. We model different types of decaying cosmic MHD turbulence in the expanding universe and characterize the large-scale magnetic fields in such a medium. Direct numerical simulations of freely decaying MHD turbulence are performed for different magnetogenesis scenarios: magnetic fields generated during cosmic inflation as well as electroweak and QCD phase transitions in the early universe. Magnetic fields and fluid motions are strongly coupled due to the high Reynolds number in the early universe. Hence, we abandon the simple adiabatic dilution model to estimate magnetic field amplitudes in the expanding universe and include turbulent mixing effects on the large-scale magnetic field evolution. Numerical simulations have been carried out for non-helical and helical magnetic field configurations. The numerical results show the possibility of inverse transfer of energy in magnetically dominated non-helical MHD turbulence. On the other hand, decay properties of helical turbulence depend on whether the turbulent magnetic field is in a weakly or a fully helical state. Our results show that primordial magnetic fields can be considered as a seed for the observed large-scale magnetic fields in galaxies and clusters. Bounds on the magnetic field strength are obtained and are consistent with the upper and lower limits set by observations of extragalactic magnetic fields.

preprint2016arXiv

Turbulent reconnection of magnetic bipoles in stratified turbulence

We consider strongly stratified forced turbulence in a plane-parallel layer with helicity and corresponding large-scale dynamo action in the lower part and non-helical turbulence in the upper. The magnetic field is found to develop strongly concentrated bipolar structures near the surface. They form elongated bands with a sharp interface between opposite polarities. Unlike earlier experiments with imposed magnetic field, the inclusion of rotation does not strongly suppress the formation of these structures. We perform a systematic numerical study of this phenomenon by varying magnetic Reynolds number, scale separation ratio, and Coriolis number. We focus on the formation of a current sheet between bipolar regions where reconnection of oppositely oriented field lines occurs. We determine the reconnection rate by measuring either the inflow velocity in the vicinity of the current sheet or by measuring the electric field in the reconnection region. We demonstrate that for large Lundquist numbers, S>10^3, the reconnection rate is nearly independent of S in agreement with results of recent numerical simulations performed by other groups in simpler settings.

preprint2016arXiv

VIEWPOINT: A New Twist in Simulating Solar Flares

Simulations show for the first time how the magnetic fields that produce solar flares can extend out of the Sun by acquiring a twist.

preprint2015arXiv

Bipolar magnetic spots from dynamos in stratified spherical shell turbulence

Recent work by Mitra et al. (2014) has shown that in strongly stratified forced two-layer turbulence with helicity and corresponding large-scale dynamo action in the lower layer, a magnetic field occurs in the upper layer in the form of sharply bounded bipolar magnetic spots. Here we extend this model to spherical wedge geometry covering the northern hemisphere up to 75° latitude and an azimuthal extent of 180°. The kinetic helicity and therefore also the large-scale magnetic field are strongest at low latitudes. For moderately strong stratification, several bipolar spots form that fill eventually the full longitudinal extent. At early times, the polarity of spots reflects the orientation of the underlying azimuthal field, as expected from Ω-shaped flux loops. At late times their tilt changes such that there is a radial field of opposite orientation at different latitudes separated by about 10°. Our model demonstrates the spontaneous formation of spots of sizes much larger than the pressure scale height. Their tendency to produce filling factors close to unity is argued to be reminiscent of highly active stars. We confirm that strong stratification and strong scale separation are essential ingredients behind magnetic spot formation, which appears to be associated with downflows at larger depths.

preprint2015arXiv

Hysteresis between distinct modes of turbulent dynamos

Nonlinear mean-field models of the solar dynamo show long-term variability, which may be relevant to different states of activity inferred from long-term radiocarbon data. This paper is aimed to probe the dynamo hysteresis predicted by the recent mean-field models of Kitchatinov \& Olemskoy (2010) with direct numerical simulations. We perform three-dimensional simulations of large-scale dynamos in a shearing box with helically forced turbulence. As initial condition, we either take a weak random magnetic field or we start from a snapshot of an earlier simulation. Two quasi-stable states are found to coexist in a certain range of parameters close to the onset of the large-scale dynamo. The simulations converge to one of these states depending on the initial conditions. When either the fractional helicity or the magnetic Prandtl number is increased between successive runs above the critical value for onset of the dynamo, the field strength jumps to a finite value. However, when the fractional helicity or the magnetic Prandtl number is then decreased again, the field strength stays at a similar value (strong field branch) even below the original onset. We also observe intermittent decaying phases away from the strong field branch close to the point where large-scale dynamo action is just possible. The dynamo hysteresis seen previously in mean-field models is thus reproduced by 3D simulations. Its possible relation to distinct modes of solar activity such as grand minima is discussed.

preprint2015arXiv

Is the small-scale magnetic field correlated with the dynamo cycle?

The small-scale magnetic field is ubiquitous at the solar surface---even at high latitudes. From observations we know that this field is uncorrelated (or perhaps even weakly anticorrelated) with the global sunspot cycle. Our aim is to explore the origin, and particularly the cycle dependence of such a phenomenon using three-dimensional dynamo simulations. We adopt a simple model of a turbulent dynamo in a shearing box driven by helically forced turbulence. Depending on the dynamo parameters, large-scale (global) and small-scale (local) dynamos can be excited independently in this model. Based on simulations in different parameter regimes, we find that, when only the large-scale dynamo is operating in the system, the small-scale magnetic field generated through shredding and tangling of the large-scale magnetic field is positively correlated with the global magnetic cycle. However, when both dynamos are operating, the small-scale field is produced from both the small-scale dynamo and the tangling of the large-scale field. In this situation, when the large-scale field is weaker than the equipartition value of the turbulence, the small-scale field is almost uncorrelated with the large-scale magnetic cycle. On the other hand, when the large-scale field is stronger than the equipartition value, we observe an anticorrelation between the small-scale field and the large-scale magnetic cycle. This anticorrelation can be interpreted as a suppression of the small-scale dynamo. Based on our studies we conclude that the observed small-scale magnetic field in the Sun is generated by the combined mechanisms of small-scale dynamo and tangling of the large-scale field.

preprint2015arXiv

Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

preprint2015arXiv

Magnetically controlled stellar differential rotation near the transition from solar to anti-solar profiles

Late-type stars rotate differentially owing to anisotropic turbulence in their outer convection zones. The rotation is called solar-like (SL) when the equator rotates fastest and anti-solar (AS) otherwise. Hydrodynamic simulations show a transition from SL to AS rotation as the influence of rotation on convection is reduced, but the opposite transition occurs at a different point in the parameter space. The system is bistable, i.e., SL and AS rotation profiles can both be stable. We study the effect of a dynamo-generated magnetic field on the large-scale flows, particularly on the possibility of bistable behavior of differential rotation. We solve the hydromagnetic equations numerically in a rotating spherical shell for a set of different radiative conductivities controlling the relative importance of convection. In agreement with earlier findings, our models display SL rotation profiles when the rotational influence on convection is strong and a transition to AS when the rotational influence decreases. We find that dynamo-generated magnetic fields help to produce SL differential rotation compared to the hydrodynamic simulations. We do not observe any bistable states of differential rotation. In the AS cases we get coherent single-cell meridional circulation, whereas in SL cases we get multi-cellular patterns. In both cases, we obtain poleward circulation near the surface with a magnitude close to that observed in the Sun. Moreover, both differential rotation and meridional circulation have significant magnetic cycle-related variations that are similar in strength to those of the Sun. Purely hydrodynamic simulations of differential rotation and meridional circulation are shown to be of limited relevance as magnetic fields, self-consistently generated by dynamo action, significantly affect the flows.

preprint2015arXiv

Negative magnetic eddy diffusivities from test-field method and multiscale stability theory

The generation of large-scale magnetic field in the kinematic regime in the absence of an alpha-effect is investigated by following two different approaches, namely the test-field method and multiscale stability theory relying on the homogenisation technique. We show analytically that the former, applied for the evaluation of magnetic eddy diffusivities, yields results that fully agree with the latter. Our computations of the magnetic eddy diffusivity tensor for the specific instances of the parity-invariant flow-IV of G.O. Roberts and the modified Taylor-Green flow in a suitable range of parameter values confirm the findings of previous studies, and also explain some of their apparent contradictions. The two flows have large symmetry groups; this is used to considerably simplify the eddy diffusivity tensor. Finally, a new analytic result is presented: upon expressing the eddy diffusivity tensor in terms of solutions to auxiliary problems for the adjoint operator, we derive relations between magnetic eddy diffusivity tensors that arise for opposite small-scale flows v(x) and -v(x).

preprint2015arXiv

Properties of $p$- and $f$-modes in hydromagnetic turbulence

With the ultimate aim of using the fundamental or $f$-mode to study helioseismic aspects of turbulence-generated magnetic flux concentrations, we use randomly forced hydromagnetic simulations of a piecewise isothermal layer in two dimensions with reflecting boundaries at top and bottom. We compute numerically diagnostic wavenumber-frequency diagrams of the vertical velocity at the interface between the denser gas below and the less dense gas above. For an Alfvén-to-sound speed ratio of about 0.1, a 5% frequency increase of the $f$-mode can be measured when $k_xH_{\rm p}=3$-$4$, where $k_x$ is the horizontal wavenumber and $H_{\rm p}$ is the pressure scale height at the surface. Since the solar radius is about 2000 times larger than $H_{\rm p}$, the corresponding spherical harmonic degree would be 6000-8000. For weaker fields, a $k_x$-dependent frequency decrease by the turbulent motions becomes dominant. For vertical magnetic fields, the frequency is enhanced for $k_xH_{\rm p}\approx4$, but decreased relative to its nonmagnetic value for $k_xH_{\rm p}\approx9$.

preprint2014arXiv

Evolution of magnetic field generated by the Kelvin-Helmholtz instability

The Kelvin-Helmholtz instability in an ionized plasma is considered with a focus on the generation of magnetic field via the Biermann battery mechanisms. The problem is studied through direct numerical simulations of two counter-directed flows in 2D geometry. The simulations demonstrate the formation of eddies and their further interaction resulting in a large single vortex. At early stages, the generated magnetic field evolves due to the baroclinic term in the induction equation, revealing significantly different structures from the vorticity field, despite the fact that magnetic field and vorticity obey identical equations. At later times, the magnetic field exhibits complex behavior and continues to grow even after a hydrodynamic vortex has developed.

preprint2014arXiv

Fanning out of the $f$-mode in presence of nonuniform magnetic fields

We show that in the presence of a harmonically varying magnetic field the fundamental or $f$-mode in a stratified layer is altered in such a way that it fans out in the diagnostic $kω$ diagram, but with mode power also within the fan. In our simulations, the surface is defined by a temperature and density jump in a piecewise isothermal layer. Unlike our previous work (Singh et al. 2014) where a uniform magnetic field was considered, we employ here a nonuniform magnetic field together with hydromagnetic turbulence at length scales much smaller than those of the magnetic fields. The expansion of the $f$-mode is stronger for fields confined to the layer below the surface. In some of those cases, the $kω$ diagram also reveals a new class of low frequency vertical stripes at multiples of twice the horizontal wavenumber of the background magnetic field. We argue that the study of the $f$-mode expansion might be a new and sensitive tool to determining subsurface magnetic fields with longitudinal periodicity.

preprint2014arXiv

Faraday signature of magnetic helicity from reduced depolarization

Using one-dimensional models, we show that a helical magnetic field with an appropriate sign of helicity can compensate the Faraday depolarization resulting from the superposition of Faraday-rotated polarization planes from a spatially extended source. For radio emission from a helical magnetic field, the polarization as a function of the square of the wavelength becomes asymmetric with respect to zero. Mathematically speaking, the resulting emission occurs then either at observable or at unobservable (imaginary) wavelengths. We demonstrate that rotation measure (RM) synthesis allows for the reconstruction of the underlying Faraday dispersion function in the former case, but not in the latter. The presence of positive magnetic helicity can thus be detected by observing positive RM in highly polarized regions in the sky and negative RM in weakly polarized regions. Conversely, negative magnetic helicity can be detected by observing negative RM in highly polarized regions and positive RM in weakly polarized regions. The simultaneous presence of two magnetic constituents with opposite signs of helicity is shown to possess signatures that can be quantified through polarization peaks at specific wavelengths and the gradient of the phase of the Faraday dispersion function. Similar polarization peaks can tentatively also be identified for the bi-helical magnetic fields that are generated self-consistently by a dynamo from helically forced turbulence, even though the magnetic energy spectrum is then continuous. Finally, we discuss the possibility of detecting magnetic fields with helical and non-helical properties in external galaxies using the Square Kilometre Array.

preprint2014arXiv

Hydraulic effects in a radiative atmosphere with ionization

In a paper of 1978, Eugene Parker postulated the need for hydraulic downward motion to explain magnetic flux concentrations at the solar surface. A similar process has recently also been seen in simplified (e.g., isothermal) models of flux concentrations from the negative effective magnetic pressure instability. We study the effects of partial ionization near the radiative surface on the formation of such magnetic flux concentrations. We first obtain one-dimensional (1D) equilibrium solutions using either a Kramers-like opacity or the ${\rm H}^{-}$ opacity. The resulting atmospheres are then used as initial conditions in two-dimensional (2D) models where flows are driven by an imposed gradient force resembling a localized negative pressure in the form of a blob. To isolate the effects of partial ionization and radiation, we ignore turbulence and convection. In 1D models, due to partial ionization, an unstable stratification forms always near the surface. We show that the extrema in the specific entropy profiles correspond to the extrema in degree of ionization. In the 2D models without partial ionization, flux concentrations form close to the height where the blob is placed. In models with partial ionization, such flux concentrations form at the surface much above the blob. This is due to the corresponding unstable layer in specific entropy. With ${\rm H}^{-}$ opacity, flux concentrations are weaker due to the stably stratified deeper parts. We demonstrate that, together with density stratification, the imposed source of negative pressure drives the formation of flux concentrations. We find that the inclusion of partial ionization affects entropy profiles causing the strong flux concentrations to form closer to the surface. We speculate that turbulence is needed to limit the strength of flux concentrations and homogenize the specific entropy to a more nearly marginal stratification.

preprint2014arXiv

Magnetic flux concentrations from dynamo-generated fields

The mean-field theory of magnetized stellar convection gives rise to the two possibility of distinct instabilities: the large-scale dynamo instability, operating in the bulk of the convection zone, and a negative effective magnetic pressure instability (NEMPI) operating in the strongly stratified surface layers. The latter might be important in connection with magnetic spot formation, but the growth rate of NEMPI is suppressed with increasing rotation rates, although recent direct numerical simulations (DNS) have shown a subsequent increase in the growth rate. We examine quantitatively whether this increase in the growth rate of NEMPI can be explained by an alpha squared mean-field dynamo, and whether both NEMPI and the dynamo instability can operate at the same time. We use both DNS and mean-field simulations (MFS) to solve the underlying equations numerically either with or without an imposed horizontal field. We use the test-field method to compute relevant dynamo coefficients. DNS show that magnetic flux concentrations are still possible up to rotation rates above which the large-scale dynamo effect produces mean magnetic fields. The resulting DNS growth rates are quantitatively well reproduced with MFS. As expected, for weak or vanishing rotation, the growth rate of NEMPI increases with increasing gravity, but there is a correction term for strong gravity and large turbulent magnetic diffusivity. Magnetic flux concentrations are still possible for rotation rates above which dynamo action takes over. For the solar rotation rate, the corresponding turbulent turnover time is about 5 hours, with dynamo action commencing in the layers beneath.

preprint2014arXiv

Magnetic helicity and energy spectra of a solar active region

We compute for the first time magnetic helicity and energy spectra of the solar active region NOAA 11158 during 11-15 February 2011 at 20^o southern heliographic latitude using observational photospheric vector magnetograms. We adopt the isotropic representation of the Fourier-transformed two-point correlation tensor of the magnetic field. The sign of magnetic helicity turns out to be predominantly positive at all wavenumbers. This sign is consistent with what is theoretically expected for the southern hemisphere. The magnetic helicity normalized to its theoretical maximum value, here referred to as relative helicity, is around 4% and strongest at intermediate wavenumbers of k ~ 0.4 Mm^{-1}, corresponding to a scale of 2pi/k ~ 16 Mm. The same sign and a similar value are also found for the relative current helicity evaluated in real space based on the vertical components of magnetic field and current density. The modulus of the magnetic helicity spectrum shows a k^{-11/3} power law at large wavenumbers, which implies a k^{-5/3} spectrum for the modulus of the current helicity. A k^{-5/3} spectrum is also obtained for the magnetic energy. The energy spectra evaluated separately from the horizontal and vertical fields agree for wavenumbers below 3 Mm^{-1}, corresponding to scales above 2 Mm. This gives some justification to our assumption of isotropy and places limits resulting from possible instrumental artefacts at small scales.

preprint2014arXiv

Magnetic Prandtl number dependence of kinetic to magnetic dissipation ratio

Using direct numerical simulations of three-dimensional hydromagnetic turbulence, either with helical or non-helical forcing, we show that the ratio of kinetic-to-magnetic energy dissipation always increases with the magnetic Prandtl number, i.e., the ratio of kinematic viscosity to magnetic diffusivity. This dependence can always be approximated by a power law, but the exponent is not the same in all cases. For non-helical turbulence, the exponent is around 1/3, while for helical turbulence it is between 0.6 and 2/3. In the statistically steady state, the rate of the energy conversion from kinetic into magnetic by the dynamo must be equal to the Joule dissipation rate. We emphasize that for both small-scale and large-scale dynamos, the efficiency of energy conversion depends sensitively on the magnetic Prandtl number, and thus on the microphysical dissipation process. To understand this behavior, we also study shell models of turbulence and one-dimensional passive and active scalar models. We conclude that the magnetic Prandtl number dependence is qualitatively best reproduced in the one-dimensional model as a result of dissipation via localized Alfven kinks.

preprint2014arXiv

Mean-field dynamo action from delayed transport

We analyze the nature of dynamo action that produces horizontally averaged magnetic fields in two particular flows that were studied by Roberts (1972, Phil. Trans. R. Soc. A 271, 411), namely his flows II and III. They have zero kinetic helicity either pointwise (flow II), or on average (flow III). Using direct numerical simulations, we determine the onset conditions for dynamo action at moderate values of the magnetic Reynolds number. Using the test-field method, we show that the turbulent magnetic diffusivity is then positive for both flows. However, we demonstrate that for both flows large-scale dynamo action occurs through delayed transport. Mathematically speaking, the magnetic field at earlier times contributes to the electromotive force through the off-diagonal components of the alpha tensor such that a zero mean magnetic field becomes unstable to dynamo action. This represents a qualitatively new mean-field dynamo mechanism not previously described.

preprint2014arXiv

Near-polytropic simulations with a radiative surface

Studies of solar and stellar convection often employ simple polytropic setups using the diffusion approximation instead of solving the proper radiative transfer equation. This allows one to control separately the polytropic index of the hydrostatic reference solution, the temperature contrast between top and bottom, and the Rayleigh and Peclet numbers. We extend such studies by including radiative transfer in the gray approximation using a Kramers-like opacity with freely adjustable coefficients. We study the properties of such models and compare them with results from the diffusion approximation. We use the Pencil Code, which is a high-order finite difference code where radiation is treated using the method of long characteristics. The source function is given by the Planck function. The opacity is written as kappa=kappa_0 rho^a T^b, where b is varied from -3.5 to +5, and kappa_0 is varied by four orders of magnitude. We consider sets of one dimensional models and perform a comparison with the diffusion approximation. Except for the case where b=5, we find one-dimensional hydrostatic equilibria with a nearly polytropic stratification and a polytropic index close to n=(3-b)/(1+a), covering both convectively stable (n>3/2) and unstable (n<3/2) cases. For b=3 and a=-1, the value of n is undefined a priori and the actual value of $n$ depends then on the depth of the domain. For large values of $κ_0$, the thermal adjustment time becomes long, the Peclet and Rayleigh numbers become large, and the temperature contrast increases and is thus no longer an independent input parameter, unless the Stefan Boltzmann constant is considered adjustable. Proper radiative transfer with Kramers-like opacities provides a useful tool for studying stratified layers with a radiative surface in ways that are more physical than what is possible with polytropic models using the diffusion approximation.

preprint2014arXiv

Nonhelical inverse transfer of a decaying turbulent magnetic field

In the presence of magnetic helicity, inverse transfer from small to large scales is well known in magnetohydrodynamic (MHD) turbulence and has applications in astrophysics, cosmology, and fusion plasmas. Using high resolution direct numerical simulations of magnetically dominated self-similarly decaying MHD turbulence, we report a similar inverse transfer even in the absence of magnetic helicity. We compute for the first time spectral energy transfer rates to show that this inverse transfer is about half as strong as with helicity, but in both cases the magnetic gain at large scales results from velocity at similar scales interacting with smaller-scale magnetic fields. This suggests that both inverse transfers are a consequence of a universal mechanisms for magnetically dominated turbulence. Possible explanations include inverse cascading of the mean squared vector potential associated with local near two-dimensionality and the shallower $k^2$ subinertial range spectrum of kinetic energy forcing the magnetic field with a $k^4$ subinertial range to attain larger-scale coherence. The inertial range shows a clear $k^{-2}$ spectrum and is the first example of fully isotropic magnetically dominated MHD turbulence exhibiting weak turbulence scaling.

preprint2014arXiv

On the cause of solar-like equatorward migration in global convective dynamo simulations

We present results from four convectively-driven stellar dynamo simulations in spherical wedge geometry. All of these simulations produce cyclic and migrating mean magnetic fields. Through detailed comparisons we show that the migration direction can be explained by an $αΩ$ dynamo wave following the Parker--Yoshimura rule. We conclude that the equatorward migration in this and previous work is due to a positive (negative) $α$ effect in the northern (southern) hemisphere and a negative radial gradient of $Ω$ outside the inner tangent cylinder of these models. This idea is supported by a strong correlation between negative radial shear and toroidal field strength in the region of equatorward propagation.

preprint2014arXiv

Particle energization through time-periodic helical magnetic fields

We solve for the motion of charged particles in a helical time-periodic ABC (Arnold-Beltrami-Childress) magnetic field. The magnetic field lines of a stationary ABC field with coefficients $A=B=C=1$ are chaotic, and we show that the motion of a charged particle in such a field is also chaotic at late times with positive Lyapunov exponent. We further show that in time-periodic ABC fields, the kinetic energy of a charged particle can increase indefinitely with time. At late times the mean kinetic energy grows as a power law in time with an exponent that approaches unity. For an initial distribution of particles, whose kinetic energy is uniformly distributed within some interval, the PDF of kinetic energy is, at late times, close to a Gaussian but with steeper tails.

preprint2014arXiv

Quenching and anisotropy of hydromagnetic turbulent transport

Hydromagnetic turbulence affects the evolution of large-scale magnetic fields through mean-field effects like turbulent diffusion and the $α$ effect. For stronger fields, these effects are usually suppressed or quenched, and additional anisotropies are introduced. Using different variants of the test-field method, we determine the quenching of the turbulent transport coefficients for the forced Roberts flow, isotropically forced non-helical turbulence, and rotating thermal convection. We see significant quenching only when the mean magnetic field is larger than the equipartition value of the turbulence. Expressing the magnetic field in terms of the equipartition value of the {\it quenched} flows, we obtain for the quenching exponents of the turbulent magnetic diffusivity about 1.3, 1.1, and 1.3 for Roberts flow, forced turbulence, and convection, respectively. However, when the magnetic field is expressed in terms of the equipartition value of the unquenched flows these quenching exponents become about 4, 1.5, and 2.3, respectively. For the $α$ effect, the exponent is about 1.3 for the Roberts flow and 2 for convection in the first case, but 4 and 3, respectively, in the second. In convection, the quenching of turbulent pumping follows the same power law as turbulent diffusion, while for the coefficient describing the $\bf Ω\times \bf J$ effect nearly the same quenching exponent is obtained as for $α$. For forced turbulence, turbulent diffusion proportional to the second derivative along the mean magnetic field is quenched much less, especially for larger values of the magnetic Reynolds number. However, we find that in corresponding axisymmetric mean-field dynamos with dominant toroidal field the quenched diffusion coefficients are the same for the poloidal and toroidal field constituents.

preprint2014arXiv

Simulations of galactic dynamos

We review our current understanding of galactic dynamo theory, paying particular attention to numerical simulations both of the mean-field equations and the original three-dimensional equations relevant to describing the magnetic field evolution for a turbulent flow. We emphasize the theoretical difficulties in explaining non-axisymmetric magnetic fields in galaxies and discuss the observational basis for such results in terms of rotation measure analysis. Next, we discuss nonlinear theory, the role of magnetic helicity conservation and magnetic helicity fluxes. This leads to the possibility that galactic magnetic fields may be bi-helical, with opposite signs of helicity and large and small length scales. We discuss their observational signatures and close by discussing the possibilities of explaining the origin of primordial magnetic fields.

preprint2014arXiv

Superflare occurrence and energies on G-, K- and M-type dwarfs

Kepler data from G-, K- and M-type stars are used to study conditions that lead to superflares with energies above $10^{34} {\rm erg}$. From the 117,661 stars included, 380 show superflares with a total of 1690 such events. We study whether parameters, like effective temperature or the rotation rate, have any effect on the superflare occurrence rate or energy. With increasing effective temperature we observe a decrease in the superflare rate, which is analogous to the previous findings of a decrease in dynamo activity with increasing effective temperature. For slowly rotating stars, we find a quadratic increase of the mean occurrence rate with the rotation rate up to a critical point, after which the rate decreases linearly. Motivated by standard dynamo theory, we study the behavior of the relative starspot coverage, approximated as the relative brightness variation. For faster rotating stars, an increased fraction of stars shows higher spot coverage, which leads to higher superflare rates. A turbulent dynamo is used to study the dependence of the Ohmic dissipation as a proxy of the flare energy on the differential rotation or shear rate. The resulting statistics of the dissipation energy as a function of dynamo number is similar to the observed flare statistics as a function of the inverse Rossby number and shows similarly strong fluctuations. This supports the idea that superflares might well be possible for solar-type G stars.

preprint2014arXiv

Surface flux concentrations and spherical alpha-square dynamo

In the presence of strong density stratification, turbulence can lead to a large-scale instability of a horizontal magnetic field if its strength is in a suitable range (within a few percent of the turbulent equipartition value). This instability is related to a suppression of the turbulent pressure so that the turbulence contribution to the mean magnetic pressure becomes negative. This results in the excitation of a negative effective magnetic pressure instability (NEMPI). This instability has so far only been studied for an imposed magnetic field. We want to know how NEMPI works when the mean magnetic field is generated self-consistently by an $α^2$ dynamo, whether it is affected by global spherical geometry, and whether it can influence the properties of the dynamo itself. We adopt the mean-field approach which has previously been shown to provide a realistic description of NEMPI in direct numerical simulations. We assume axisymmetry and solve the mean-field equations with the Pencil-Code for an adiabatic stratification at a total density contrast in the radial direction of approximately 4 orders of magnitude. NEMPI is found to work when the dynamo-generated field is about 4% of the equipartition value, which is achieved through strong $α$ quenching. This instability is excited in the top 5% of the outer radius provided the density contrast across this top layer is at least 10. NEMPI is found to occur at lower latitudes when the mean magnetic field is stronger. For weaker fields, NEMPI can make the dynamo oscillatory with poleward migration. NEMPI is a viable mechanism for producing magnetic flux concentrations in a strongly stratified spherical shell in which a magnetic field is generated by a strongly quenched $α$ effect dynamo.

preprint2014arXiv

Traces of large-scale dynamo action in the kinematic stage

Using direct numerical simulations (DNS) we verify that in the kinematic regime, a turbulent helical dynamo grows in such a way that the magnetic energy spectrum remains to high precision shape-invariant, i.e., at each wavenumber $k$ the spectrum grows with the same growth rate. Signatures of large-scale dynamo action can be identified through the excess of magnetic energy at small $k$, of one of the two oppositely polarized constituents. Also a suitably defined planar average of the magnetic field can be chosen such that its rms value isolates the strength of the mean field. However, these different means of analysis suggest that the strength of the large-scale field diminishes with increasing magnetic Reynolds number ${\rm Re}_{\rm M}$ like ${\rm Re}_{\rm M}^{-1/2}$ for intermediate values and like ${\rm Re}_{\rm M}^{-3/4}$ for larger ones. Both an analysis from the Kazantsev model including helicity and the DNS show that this arises due to the magnetic energy spectrum still peaking at resistive scales, even when helicity is present. As expected, the amplitude of the large-scale field increases with increasing fractional helicity, enabling us to determine the onset of large scale dynamo action and distinguishing it from that of the small-scale dynamo. Our DNS show that, contrary to earlier results for smaller scale separation (only 1.5 instead of now 4), the small-scale dynamo can still be excited at magnetic Prandtl numbers of 0.1 and only moderate values of the magnetic Reynolds numbers ($\sim 160$).

preprint2013arXiv

A mean field dynamo from negative eddy diffusivity

Using direct numerical simulations, we verify that "flow IV" of Roberts (1972) exhibits dynamo action dominated by horizontally averaged large-scale magnetic field. With the test-field method we compute the turbulent magnetic diffusivity and find that it is negative and overcomes the molecular diffusivity, thus explaining quantitatively the large-scale dynamo for magnetic Reynolds numbers above $\approx8$. As expected for a dynamo of this type, but contrary to $α$-effect dynamos, the two horizontal field components grow independently of each other and have arbitrary amplitude ratios and phase differences. Small length scales of the mean magnetic field are shown to be stabilized by the turbulent magnetic diffusivity becoming positive at larger wavenumbers. Oscillatory decaying or growing solutions have also been found in certain wavenumber intervals and sufficiently large values of the magnetic Reynolds number. For magnetic Reynolds numbers below $\approx0.5$ the turbulent magnetic diffusivity is confirmed to be positive, as expected for all incompressible flows. Earlier claims of a dynamo driven by a modified Taylor-Green flow through negative eddy diffusivity could not be confirmed.

preprint2013arXiv

Azimuthal dynamo wave in spherical shell convection

We report the finding of an azimuthal dynamo wave of a low-order (m=1) mode in direct numerical simulations (DNS) of turbulent convection in spherical shells. Such waves are predicted by mean field dynamo theory and have been obtained previously in mean-field models. Observational results both from photometry and Doppler imaging have revealed persistent drifts of spots for several rapidly rotating stars, but, although an azimuthal dynamo wave has been proposed as a possible mechanism responsible for this behavior, it has been judged as unlikely, as practical evidence for such waves from DNS has been lacking. The large-scale magnetic field in our DNS, which is due to self-consistent dynamo action, is dominated by a retrograde m=1 mode. Its pattern speed is nearly independent of latitude and does not reflect the speed of the differential rotation at any depth. The extrema of magnetic m=1 structures coincide reasonably with the maxima of m=2 structures of the temperature. These results provide direct support for the observed drifts being due to an azimuthal dynamo wave.

preprint2013arXiv

Bipolar magnetic structures driven by stratified turbulence with a coronal envelope

We report the spontaneous formation of bipolar magnetic structures in direct numerical simulations of stratified forced turbulence with an outer coronal envelope. The turbulence is forced with transverse random waves only in the lower (turbulent) part of the domain. Our initial magnetic field is either uniform in the entire domain or confined to the turbulent layer. After about 1--2 turbulent diffusion times, a bipolar magnetic region of vertical field develops with two coherent circular structures that live during one turbulent diffusion time, and then decay during 0.5 turbulent diffusion times. The resulting magnetic field strengths inside the bipolar region are comparable to the equipartition value with respect to the turbulent kinetic energy. The bipolar magnetic region forms a loop-like structure in the upper coronal layer. We associate the magnetic structure formation with the negative effective magnetic pressure instability in the two-layer model.

preprint2013arXiv

Can planetesimals form by collisional fusion?

As a test bed for the growth of protoplanetary bodies in a turbulent circumstellar disk we examine the fate of a boulder using direct numerical simulations of particle seeded gas flowing around it. We provide an accurate description of the flow by imposing no-slip and non-penetrating boundary conditions on the boulder surface using the immersed boundary method pioneered by Peskin (2002). Advected by the turbulent disk flow, the dust grains collide with the boulder and we compute the probability density function (PDF) of the normal component of the collisional velocity. Through this examination of the statistics of collisional velocities we test the recently developed concept of collisional fusion which provides a physical basis for a range of collisional velocities exhibiting perfect sticking. A boulder can then grow sufficiently rapidly to settle into a Keplerian orbit on disk evolution time scales.

preprint2013arXiv

Coronal influence on dynamos

We report on turbulent dynamo simulations in a spherical wedge with an outer coronal layer. We apply a two-layer model where the lower layer represents the convection zone and the upper layer the solar corona. This setup is used to study the coronal influence on the dynamo action beneath the surface. Increasing the radial coronal extent gradually to three times the solar radius and changing the magnetic Reynolds number, we find that dynamo action benefits from the additional coronal extent in terms of higher magnetic energy in the saturated stage. The flux of magnetic helicity can play an important role in this context.

preprint2013arXiv

Cyclic magnetic activity due to turbulent convection in spherical wedge geometry

We report on simulations of turbulent, rotating, stratified, magnetohydrodynamic convection in spherical wedge geometry. An initially small-scale, random, weak-amplitude magnetic field is amplified by several orders of magnitude in the course of the simulation to form oscillatory large-scale fields in the saturated state of the dynamo. The differential rotation is solar-like (fast equator), but neither coherent meridional poleward circulation nor near-surface shear layer develop in these runs. In addition to a poleward branch of magnetic activity beyond 50 degrees latitude, we find for the first time a pronounced equatorward branch at around 20 degrees latitude, reminiscent of the solar cycle.

preprint2013arXiv

Effects of enhanced stratification on equatorward dynamo wave propagation

We present results from simulations of rotating magnetized turbulent convection in spherical wedge geometry representing parts of the latitudinal and longitudinal extents of a star. Here we consider a set of runs for which the density stratification is varied, keeping the Reynolds and Coriolis numbers at similar values. In the case of weak stratification, we find quasi-steady dynamo solutions for moderate rotation and oscillatory ones with poleward migration of activity belts for more rapid rotation. For stronger stratification, the growth rate tends to become smaller. Furthermore, a transition from quasi-steady to oscillatory dynamos is found as the Coriolis number is increased, but now there is an equatorward migrating branch near the equator. The breakpoint where this happens corresponds to a rotation rate that is about 3-7 times the solar value. The phase relation of the magnetic field is such that the toroidal field lags behind the radial field by about $π/2$, which can be explained by an oscillatory $α^2$ dynamo caused by the sign change of the $α$-effect about the equator. We test the domain size dependence of our results for a rapidly rotating run with equatorward migration by varying the longitudinal extent of our wedge. The energy of the axisymmetric mean magnetic field decreases as the domain size increases and we find that an $m=1$ mode is excited for a full $2π$ azimuthal extent, reminiscent of the field configurations deduced from observations of rapidly rotating late-type stars.

preprint2013arXiv

Evolution of Primordial Magnetic Fields from Phase Transitions

We consider the evolution of primordial magnetic fields generated during cosmological, electroweak or QCD, phase transitions. We assume that the magnetic field generation can be described as an injection of magnetic energy to cosmological plasma at a given scale determined by the moment of magnetic field generation. A high Reynolds number ensures strong coupling between magnetic field and fluid motions. The subsequent evolution of the magnetic field is governed by decaying hydromagnetic turbulence. Both our numerical simulations and a phenomenological description allow us to recover "universal" laws for the decay of magnetic energy and the growth of magnetic correlation length in the turbulent (low viscosity) regime. In particular, we show that during the radiation dominated epoch, energy and correlation length of non-helical magnetic fields scale as conformal time to the powers -1/2 and +1/2, respectively. For helical magnetic fields, energy and correlation length scale as conformal time to the powers -1/3 and +2/3, respectively. The universal decay law of the magnetic field implies that the strength of magnetic field generated during the QCD phase transition could reach $\sim 10^{-9}$\,G with the present day correlation length $\sim 50$ kpc. The fields generated at the electroweak phase transition could be as strong as $\sim 10^{-10}$ G with correlation lengths reaching $\sim 0.3$\,kpc. These values of the magnetic fields are consistent with the lower bounds of the extragalactic magnetic fields.

preprint2013arXiv

Non-linear and chaotic dynamo regimes

An update is given on the current status of solar and stellar dynamos. At present, it is still unclear why stellar cycle frequencies increase with rotation frequency in such a way that their ratio increases with stellar activity. The small-scale dynamo is expected to operate in spite of a small value of the magnetic Prandtl number in stars. Whether or not the global magnetic activity in stars is a shallow or deeply rooted phenomenon is another open question. Progress in demonstrating the presence and importance of magnetic helicity fluxes in dynamos is briefly reviewed, and finally the role of nonlocality is emphasized in modeling stellar dynamos using the mean-field approach. On the other hand, direct numerical simulations have now come to the point where the models show solar-like equatorward migration that can be compared with observations and that need to be understood theoretically.

preprint2013arXiv

Scaling and intermittency in incoherent α-shear dynamo

We consider mean-field dynamo models with fluctuating αeffect, both with and without shear. The αeffect is chosen to be Gaussian white noise with zero mean and given covariance. We show analytically that the mean magnetic field does not grow, but, in an infinitely large domain, the mean-squared magnetic field shows exponential growth of the fastest growing mode at a rate proportional to the shear rate, which agrees with earlier numerical results of Yousef et al (2008) and recent analytical treatment by Heinemann et al (2011) who use a method different from ours. In the absence of shear, an incoherent α^2 dynamo may also be possible. We further show by explicit calculation of the growth rate of third and fourth order moments of the magnetic field that the probability density function of the mean magnetic field generated by this dynamo is non-Gaussian.

preprint2013arXiv

Self-assembly of shallow magnetic spots through strongly stratified turbulence

Recent studies have demonstrated that in fully developed turbulence, the effective magnetic pressure of a large-scale field (non-turbulent plus turbulent contributions) can become negative. In the presence of strongly stratified turbulence, this was shown to lead to a large-scale instability that produces spontaneous magnetic flux concentrations. Furthermore, using a horizontal magnetic field, elongated flux concentrations with a strength of a few per cent of the equipartition value were found. Here we show that a uniform vertical magnetic field leads to circular magnetic spots of equipartition field strengths. This could represent a minimalistic model of sunspot formation and highlights the importance of two critical ingredients: turbulence and strong stratification. Radiation, ionization, and supergranulation may be important for realistic simulations, but are not critical at the level of a minimalistic model of magnetic spot formation.

preprint2013arXiv

Spoke-like differential rotation in a convective dynamo with a coronal envelope

We report on the results of four convective dynamo simulations with an outer coronal layer. The magnetic field is self-consistently generated by the convective motions beneath the surface. Above the convection zone, we include a polytropic layer that extends to 1.6 solar radii. The temperature increases in this region to $\approx8$ times the value at the surface, corresponding to $\approx1.2$ times the value at the bottom of the spherical shell. We associate this region with the solar corona. We find solar-like differential rotation with radial contours of constant rotation rate, together with a near-surface shear layer. This non-cylindrical rotation profile is caused by a non-zero latitudinal entropy gradient that offsets the Taylor--Proudman balance through the baroclinic term. The meridional circulation is multi-cellular with a solar-like poleward flow near the surface at low latitudes. In most of the cases, the mean magnetic field is oscillatory with equatorward migration in two cases. In other cases, the equatorward migration is overlaid by stationary or even poleward migrating mean fields.

preprint2013arXiv

The kinetic helicity needed to drive large-scale dynamos

Magnetic field generation on scales large compared with the scale of the turbulent eddies is known to be possible via the so-called $α$ effect when the turbulence is helical and if the domain is large enough for the $α$ effect to dominate over turbulent diffusion. Using three-dimensional turbulence simulations, we show that the energy of the resulting mean magnetic field of the saturated state increases linearly with the product of normalized helicity and the ratio of domain scale to eddy scale, provided this product exceeds a critical value of around unity. This implies that large-scale dynamo action commences when the normalized helicity is larger than the inverse scale ratio. Our results show that the emergence of small-scale dynamo action does not have any noticeable effect on the large-scale dynamo. Recent findings by Pietarila Graham et al. (2012, Phys. Rev. E85, 066406) of a smaller minimal helicity may be an artifact due to the onset of small-scale dynamo action at large magnetic Reynolds numbers. However, the onset of large-scale dynamo action is difficult to establish when the kinetic helicity is small. Instead of random forcing, they used an ABC-flow with time-dependent phases. We show that such dynamos saturate prematurely in a way that is reminiscent of inhomogeneous dynamos with internal magnetic helicity fluxes. Furthermore, even for very low fractional helicities, such dynamos display large-scale fields that change direction, which is uncharacteristic of turbulent dynamos.

preprint2012arXiv

Active region formation through the negative effective magnetic pressure instability

The negative effective magnetic pressure instability operates on scales encompassing many turbulent eddies and is here discussed in connection with the formation of active regions near the surface layers of the Sun. This instability is related to the negative contribution of turbulence to the mean magnetic pressure that causes the formation of large-scale magnetic structures. For an isothermal layer, direct numerical simulations and mean-field simulations of this phenomenon are shown to agree in many details in that their onset occurs at the same depth. This depth increases with increasing field strength, such that the maximum growth rate of this instability is independent of the field strength, provided the magnetic structures are fully contained within the domain. A linear stability analysis is shown to support this finding. The instability also leads to a redistribution of turbulent intensity and gas pressure that could provide direct observational signatures.

preprint2012arXiv

Breakdown of chiral symmetry during saturation of the Tayler instability

We study spontaneous breakdown of chiral symmetry during the nonlinear evolution of the Tayler instability. We start with an initial steady state of zero helicity. Within linearized perturbation calculations, helical perturbations of this initial state have the same growth rate for either sign of helicity. Direct numerical simulations (DNS) of the fully nonlinear equations, however, show that an infinitesimal excess of one sign of helicity in the initial perturbation gives rise to a saturated helical state. We further show that this symmetry breaking can be described by weakly nonlinear finite--amplitude equations with undetermined coefficients which can be deduced solely from symmetry consideration. By fitting solutions of the amplitude equations to data from DNS we further determine the coefficients of the amplitude equations.

preprint2012arXiv

Catastrophic quenching in alpha-Omega dynamos revisited

At large magnetic Reynolds numbers, magnetic helicity evolution plays an important role in astrophysical large-scale dynamos. The recognition of this fact led to the development of the dynamical alpha quenching formalism, which predicts catastrophically low mean fields in open systems. Here we show that in oscillatory alpha-Omega dynamos this formalism predicts an unphysical magnetic helicity transfer between scales. An alternative technique is proposed where this artifact is removed by using the evolution equation for the magnetic helicity of the total field in the shearing--advective gauge. In the traditional dynamical alpha quenching formalism, this can be described as an additional magnetic helicity flux of small-scale fields that does not appear in homogeneous alpha-squared dynamos. In alpha-Omega dynamos, the alternative formalism is shown to lead to larger saturation fields than previously obtained with the traditional formalism.

preprint2012arXiv

Coherent structures and the saturation of a nonlinear dynamo

Eulerian and Lagrangian tools are used to detect coherent structures in the velocity and magnetic fields of a mean--field dynamo, produced by direct numerical simulations of the three--dimensional compressible magnetohydrodynamic equations with an isotropic helical forcing and moderate Reynolds number. Two distinct stages of the dynamo are studied, the kinematic stage, where a seed magnetic field undergoes exponential growth, and the saturated regime. It is shown that the Lagrangian analysis detects structures with greater detail, besides providing information on the chaotic mixing properties of the flow and the magnetic fields. The traditional way of detecting Lagrangian coherent structures using finite--time Lyapunov exponents is compared with a recently developed method called function M. The latter is shown to produce clearer pictures which readily permit the identification of hyperbolic regions in the magnetic field, where chaotic transport/dispersion of magnetic field lines is highly enhanced.

preprint2012arXiv

Cosmic ray current-driven turbulence and mean-field dynamo effect

We show that an alpha effect is driven by the cosmic ray Bell instability exciting left-right asymmetric turbulence. Alfven waves of a preferred polarization have maximally helical motion, because the transverse motion of each mode is parallel to its curl. We show how large-scale Alfven modes, when rendered unstable by cosmic ray streaming, can create new net flux over any finite region, in the direction of the original large-scale field. We perform direct numerical simulations (DNS) of an MHD fluid with a forced cosmic ray current and use the test-field method to determine the alpha effect and the turbulent magnetic diffusivity. As follows from DNS, the dynamics of the instability has the following stages: (i) in the early stage, the small-scale Bell instability that results in a production of small-scale turbulence is excited; (ii) in the intermediate stage, there is formation of larger-scale magnetic structures; (iii) finally, quasi-stationary large-scale turbulence is formed at a growth rate that is comparable to that expected from the dynamo instability, but its amplitude over much longer timescales remains unclear. The results of DNS are in good agreement with the theoretical estimates. It is suggested that this dynamo is what gives weakly magnetized relativistic shocks such as those from gamma ray bursts a macroscopic correlation length. It may also be important for large-scale magnetic field amplification associated with cosmic ray production and diffusive shock acceleration in supernova remnants (SNR) and blast waves from gamma ray bursts. Magnetic field amplification by Bell turbulence in SNR is found to be significant, but it is limited owing to the finite time available to the super-Alfvenicly expanding remnant. The effectiveness of the mechanisms is shown to be dependent on the shock velocity.

preprint2012arXiv

Current status of turbulent dynamo theory: From large-scale to small-scale dynamos

Several recent advances in turbulent dynamo theory are reviewed. High resolution simulations of small-scale and large-scale dynamo action in periodic domains are compared with each other and contrasted with similar results at low magnetic Prandtl numbers. It is argued that all the different cases show similarities at intermediate length scales. On the other hand, in the presence of helicity of the turbulence, power develops on large scales, which is not present in non-helical small-scale turbulent dynamos. At small length scales, differences occur in connection with the dissipation cutoff scales associated with the respective value of the magnetic Prandtl number. These differences are found to be independent of whether or not there is large-scale dynamo action. However, large-scale dynamos in homogeneous systems are shown to suffer from resistive slow-down even at intermediate length scales. The results from simulations are connected to mean field theory and its applications. Recent work on helicity fluxes to alleviate large-scale dynamo quenching, shear dynamos, nonlocal effects and magnetic structures from strong density stratification are highlighted. Several insights which arise from analytic considerations of small-scale dynamos are discussed.

preprint2012arXiv

Data assimilation for stratified convection

We show how the 3DVAR data assimilation methodology can be used in the astrophysical context of a two-dimensional convection flow. We study the way this variational approach finds best estimates of the current state of the flow from a weighted average of model states and observations. We use numerical simulations to generate synthetic observations of a vertical two-dimensional slice of the outer part of the solar convection zone for varying noise levels and implement 3DVAR when the covariance matrices are scalar. Our simulation results demonstrate the capability of 3DVAR to produce error estimates of system states between up to tree orders of magnitude below the original noise level present in the observations. This work exemplifies the importance of applying data assimilation techniques in simulations of the stratified convection.

preprint2012arXiv

Dynamics of Saturated Energy Condensation in Two-Dimensional Turbulence

In two-dimensional forced Navier-Stokes turbulence, energy cascades to the largest scales in the system to form a pair of coherent vortices known as the Bose condensate. We show, both numerically and analytically, that the energy condensation saturates and the system reaches a statistically stationary state. The time scale of saturation is inversely proportional to the viscosity and the saturation energy level is determined by both the viscosity and the force. We further show that, without sufficient resolution to resolve the small-scale enstrophy spectrum, numerical simulations can give a spurious result for the saturation energy level. We also find that the movement of the condensate is similar to the motion of an inertial particle with an effective drag force. Furthermore, we show that the profile of the saturated coherent vortices can be described by a Gaussian core with exponential wings.

preprint2012arXiv

Ejections of magnetic structures above a spherical wedge driven by a convective dynamo with differential rotation

We combine a convectively driven dynamo in a spherical shell with a nearly isothermal density-stratified cooling layer that mimics some aspects of a stellar corona to study the emergence and ejections of magnetic field structures. This approach is an extension of earlier models, where forced turbulence simulations were employed to generate magnetic fields. A spherical wedge is used which consists of a convection zone and an extended coronal region to $\approx1.5$ times the radius of the sphere. The wedge contains a quarter of the azimuthal extent of the sphere and $150\degr$ in latitude. The magnetic field is self-consistently generated by the turbulent motions due to convection beneath the surface. Magnetic fields are found to emerge at the surface and are ejected to the coronal part of the domain. These ejections occur at irregular intervals and are weaker than in earlier work. We tentatively associate these events with coronal mass ejections on the Sun, even though our model of the solar atmosphere is rather simplistic.

preprint2012arXiv

Evolution of inflation-generated magnetic field through phase transitions

We study the evolution of an inflation-generated magnetic field, due to its coupling to fluid motions, during cosmological phase transitions. We find that the magnetic field stays almost unchanged on large scales, while on small scales the spectrum is modified in such a way that power at small scales becomes progressively suppressed. We also show that the magnetic field generates turbulent motions in the initially turbulence-free plasma. On large scales, the slope of the resulting kinetic energy spectrum is consistent with that of white noise.

preprint2012arXiv

Kinetic helicity decay in linearly forced turbulence

The decay of kinetic helicity is studied in numerical models of forced turbulence using either an externally imposed forcing function as an inhomogeneous term in the equations or, alternatively, a term linear in the velocity giving rise to a linear instability. The externally imposed forcing function injects energy at the largest scales, giving rise to a turbulent inertial range with nearly constant energy flux while for linearly forced turbulence the spectral energy is maximum near the dissipation wavenumber. Kinetic helicity is injected once a statistically steady state is reached, but it is found to decay on a turbulent time scale regardless of the nature of the forcing and the value of the Reynolds number.

preprint2012arXiv

Lagrangian chaos in an ABC--forced nonlinear dynamo

The Lagrangian properties of the velocity field in a magnetized fluid are studied using three-dimensional simulations of a helical magnetohydrodynamic dynamo. We compute the attracting and repelling Lagrangian coherent structures, which are dynamic lines and surfaces in the velocity field that delineate particle transport in flows with chaotic streamlines and act as transport barriers. Two dynamo regimes are explored, one with a robust coherent mean magnetic field and one with intermittent bursts of magnetic energy. The Lagrangian coherent structures and the statistics of finite--time Lyapunov exponents indicate that the stirring/mixing properties of the velocity field decay as a linear function of the magnetic energy. The relevance of this study for the solar dynamo problem is discussed.

preprint2012arXiv

Magnetic Fields from QCD Phase Transitions

We study the evolution of QCD phase transition-generated magnetic fields in freely decaying MHD turbulence of the expanding Universe. We consider a magnetic field generation model that starts from basic non-perturbative QCD theory and predicts stochastic magnetic fields with an amplitude of the order of 0.02 $μ$G and small magnetic helicity. We employ direct numerical simulations to model the MHD turbulence decay and identify two different regimes: "weakly helical" turbulence regime, when magnetic helicity increases during decay, and "fully helical" turbulence, when maximal magnetic helicity is reached and an inverse cascade develops. The results of our analysis show that in the most optimistic scenario the magnetic correlation length in the comoving frame can reach 10 kpc with the amplitude of the effective magnetic field being 0.007 nG. We demonstrate that the considered model of magneto-genesis can provide the seed magnetic field for galaxies and clusters.

preprint2012arXiv

Magnetic helicity fluxes and their effect on stellar dynamos

Magnetic helicity fluxes in turbulently driven alpha^2 dynamos are studied to demonstrate their ability to alleviate catastrophic quenching. A one-dimensional mean-field formalism is used to achieve magnetic Reynolds numbers of the order of 10^5. We study both diffusive magnetic helicity fluxes through the mid-plane as well as those resulting from the recently proposed alternate dynamic quenching formalism. By adding shear we make a parameter scan for the critical values of the shear and forcing parameters for which dynamo action occurs. For this $αΩ$ dynamo we find that the preferred mode is antisymmetric about the mid-plane. This is also verified in 3-D direct numerical simulations.

preprint2012arXiv

Magnetic twist: a source and property of space weather

We present evidence for finite magnetic helicity density in the heliosphere and numerical models thereof, and relate it to the magnetic field properties of the dynamo in the solar convection zone. We use simulations and solar wind data to compute magnetic helicity either directly from the simulations, or indirectly using time series of the skew-symmetric components of the magnetic correlation tensor. We find that the solar dynamo produces negative magnetic helicity at small scales and positive at large scales. However, in the heliosphere these properties are reversed and the magnetic helicity is now positive at small scales and negative at large scales. We explain this by the fact that a negative diffusive magnetic helicity flux corresponds to a positive gradient of magnetic helicity, which leads to a change of sign from negative to positive values at some radius in the northern hemisphere.

preprint2012arXiv

Nonuniformity effects in the negative effective magnetic pressure instability

Using direct numerical simulations (DNS) and mean-field simulations (MFS), the effects of non-uniformity of the magnetic field on the suppression of the turbulent pressure is investigated. This suppression of turbulent pressure can lead to an instability which, in turn, makes the mean magnetic field even more non-uniform. This large-scale instability is caused by a resulting negative contribution of small-scale turbulence to the effective (mean-field) magnetic pressure. We show that enhanced mean current density increases the suppression of the turbulent pressure. The instability leads to magnetic flux concentrations in which the normalized effective mean-field pressure is reduced to a certain value at all depths within a structure.

preprint2012arXiv

Solar-like differential rotation and equatorward migration in a convective dynamo with a coronal envelope

We present results of convective turbulent dynamo simulations including a coronal layer in a spherical wedge. We find an equatorward migration of the radial and azimuthal fields similar to the behavior of sunspots during the solar cycle. The migration of the field coexist with a spoke-like differential rotation and anti-solar (clockwise) meridional circulation. Even though the migration extends over the whole convection zone, the mechanism causing this is not yet fully understood.

preprint2012arXiv

The negative effective magnetic pressure in stratified forced turbulence

To understand the basic mechanism of the formation of magnetic flux concentrations, we determine by direct numerical simulations the turbulence contributions to the mean magnetic pressure in a strongly stratified isothermal layer with large plasma beta, where a weak uniform horizontal mean magnetic field is applied. The negative contribution of turbulence to the effective mean magnetic pressure is determined for strongly stratified forced turbulence over a range of values of magnetic Reynolds and Prandtl numbers. Small-scale dynamo action is shown to reduce the negative effect of turbulence on the effective mean magnetic pressure. However, the turbulence coefficients describing the negative effective magnetic pressure phenomenon are found to be converged for magnetic Reynolds numbers between 60 and 600, which is the largest value considered here. In all these models the turbulent intensity is arranged to be nearly independent of height, so the kinetic energy density decreases with height due to the decrease in density. In a second series of numerical experiments, the turbulent intensity increases with height such that the turbulent kinetic energy density is nearly independent of height. Turbulent magnetic diffusivity and turbulent pumping velocity are determined with the test-field method for both cases. The vertical profile of the turbulent magnetic diffusivity is found to agree with what is expected based on simple mixing length expressions. Turbulent pumping is shown to be down the gradient of turbulent magnetic diffusivity, but it is twice as large as expected. Corresponding numerical mean-field models are used to show that a large-scale instability can occur in both cases, provided the degree of scale separation is large enough and hence the turbulent magnetic diffusivity small enough.

preprint2012arXiv

Topological constraints on magnetic field relaxation

Magnetic field relaxation is determined by both the field's geometry and its topology. For relaxation processes, however, it turns out that its topology is a much more stringent constraint. As quantifier for the topology we use magnetic helicity and test whether it is a stronger condition than the linking of field lines. Further, we search for evidence of other topological invariants, which give rise to further restrictions in the field's relaxation. We find that magnetic helicity is the sole determinant in most cases. Nevertheless, we see evidence for restrictions not captured through magnetic helicity.

preprint2012arXiv

Turbulent dynamo with advective magnetic helicity flux

Many astrophysical bodies harbor magnetic fields that are thought to be sustained by a dynamo process. However, it has been argued that the production of large-scale magnetic fields by mean-field dynamo action is strongly suppressed at large magnetic Reynolds numbers owing to the conservation of magnetic helicity. This phenomenon is known as {\it catastrophic quenching}. Advection of magnetic fields by stellar and galactic winds toward the outer boundaries and away from the dynamo is expected to alleviate such quenching. Here we explore the relative roles played by advective and turbulent--diffusive fluxes of magnetic helicity in the dynamo. In particular, we study how the dynamo is affected by advection. We do this by performing direct numerical simulations of a turbulent dynamo of $α^2$ type driven by forced turbulence in a Cartesian domain in the presence of a flow away from the equator where helicity changes sign. Our results indicate that in the presence of advection, the dynamo, otherwise stationary, becomes oscillatory. We confirm an earlier result for turbulent--diffusive magnetic helicity fluxes that for small magnetic Reynolds numbers ($\Rm\lesssim 100...200$, based on the wavenumber of the energy-carrying eddies) the magnetic helicity flux scales less strongly with magnetic Reynolds number ($\Rm^{-1/2}$) than the term describing magnetic helicity destruction by resistivity ($\Rm^{-1}$). Our new results now suggest that for larger $\Rm$ the former becomes approximately independent of $\Rm$, while the latter falls off more slowly. We show for the first time that both for weak and stronger winds, the magnetic helicity flux term becomes comparable to the resistive term for $\Rm\gtrsim 1000$, which is necessary for alleviating catastrophic quenching.

preprint2011arXiv

A model of driven and decaying magnetic turbulence in a cylinder

Using mean-field theory, we compute the evolution of the magnetic field in a cylinder with outer perfectly conducting boundaries, an imposed axial magnetic and electric field. The thus injected magnetic helicity in the system can be redistributed by magnetic helicity fluxes down the gradient of the local current helicity of the small-scale magnetic field. A weak reversal of the axial magnetic field is found to be a consequence of the magnetic helicity flux in the system. Such fluxes are known to alleviate so-called catastrophic quenching of the α-effect in astrophysical applications. Application to the reversed field pinch in plasma confinement devices is discussed.

preprint2011arXiv

Alleviating alpha quenching by solar wind and meridional flow

We study the ability of magnetic helicity expulsion to alleviate catastrophic $α$-quenching in mean field dynamos in two--dimensional spherical wedge domains. Motivated by the physical state of the outer regions of the Sun, we consider $α^2Ω$ mean field models with a dynamical $α$ quenching. We include two mechanisms which have the potential to facilitate helicity expulsion, namely advection by a mean flow ("solar wind") and meridional circulation. We find that a wind alone can prevent catastrophic quenching, with the field saturating at finite amplitude. In certain parameter ranges, the presence of a large-scale meridional circulation can reinforce this alleviation. However, the saturated field strengths are typically below the equipartition field strength. We discuss possible mechanisms that might increase the saturated field.

preprint2011arXiv

Angular momentum transport in convectively unstable shear flows

Angular momentum transport owing to hydrodynamic turbulent convection is studied using local three dimensional numerical simulations employing the shearing box approximation. We determine the turbulent viscosity from non-rotating runs over a range of values of the shear parameter and use a simple analytical model in order to extract the non-diffusive contribution (Lambda-effect) to the stress in runs where rotation is included. Our results suggest that the turbulent viscosity is of the order of the mixing length estimate and weakly affected by rotation. The Lambda-effect is non-zero and a factor of 2-4 smaller than the turbulent viscosity in the slow rotation regime. We demonstrate that for Keplerian shear, the angular momentum transport can change sign and be outward when the rotation period is greater than the turnover time, i.e. when the Coriolis number is below unity. This result seems to be relatively independent of the value of the Rayleigh number.

preprint2011arXiv

Astrophysical turbulence modeling

The role of turbulence in various astrophysical settings is reviewed. Among the differences to laboratory and atmospheric turbulence we highlight the ubiquitous presence of magnetic fields that are generally produced and maintained by dynamo action. The extreme temperature and density contrasts and stratifications are emphasized in connection with turbulence in the interstellar medium and in stars with outer convection zones, respectively. In many cases turbulence plays an essential role in facilitating enhanced transport of mass, momentum, energy, and magnetic fields in terms of the corresponding coarse-grained mean fields. Those transport properties are usually strongly modified by anisotropies and often completely new effects emerge in such a description that have no correspondence in terms of the original (non coarse-grained) fields.

preprint2011arXiv

Chandrasekhar-Kendall functions in astrophysical dynamos

Some of the contributions of Chandrasekhar to the field of magnetohydrodynamics are highlighted. Particular emphasis is placed on the Chandrasekhar-Kendall functions that allow a decomposition of a vector field into right- and left-handed contributions. Magnetic energy spectra of both contributions are shown for a new set of helically forced simulations at resolutions higher than what has been available so far. For a forcing function with positive helicity, these simulations show a forward cascade of the right-handed contributions to the magnetic field and nonlocal inverse transfer for the left-handed contributions. The speed of inverse transfer is shown to decrease with increasing value of the magnetic Reynolds number.

preprint2011arXiv

Coronal ejections from convective spherical shell dynamos

We present a three-dimensional model of rotating convection combined with a simplified corona in spherical coordinates. The motions in the convection zone generate a large-scale magnetic field that gets sporadically ejected into the outer layers above. Our model corona is approximately isothermal, but it includes density stratification due to gravity.

preprint2011arXiv

Cycles and cycle modulations

Some selected concepts for the solar activity cycle are briefly reviewed. Cycle modulations through a stochastic alpha effect are being identified with limited scale separation ratios. Three-dimensional turbulence simulations with helicity and shear are compared at two different scale separation ratios. In both cases the level of fluctuations shows relatively little variation with the dynamo cycle. Prospects for a shallow origin of sunspots are discussed in terms of the negative effective magnetic pressure instability. Tilt angles of bipolar active regions are discussed as a consequence of shear rather than the Coriolis force.

preprint2011arXiv

Decay of helical and non-helical magnetic knots

We present calculations of the relaxation of magnetic field structures that have the shape of particular knots and links. A set of helical magnetic flux configurations is considered, which we call $n$-foil knots of which the trefoil knot is the most primitive member. We also consider two nonhelical knots; namely, the Borromean rings as well as a single interlocked flux rope that also serves as the logo of the Inter-University Centre for Astronomy and Astrophysics in Pune, India. The field decay characteristics of both configurations is investigated and compared with previous calculations of helical and nonhelical triple-ring configurations. Unlike earlier nonhelical configurations, the present ones cannot trivially be reduced via flux annihilation to a single ring. For the $n$-foil knots the decay is described by power laws that range form $t^{-2/3}$ to $t^{-1/3}$, which can be as slow as the $t^{-1/3}$ behavior for helical triple-ring structures that were seen in earlier work. The two nonhelical configurations decay like $t^{-1}$, which is somewhat slower than the previously obtained $t^{-3/2}$ behavior in the decay of interlocked rings with zero magnetic helicity. We attribute the difference to the creation of local structures that contain magnetic helicity which inhibits the field decay due to the existence of a lower bound imposed by the realizability condition. We show that net magnetic helicity can be produced resistively as a result of a slight imbalance between mutually canceling helical pieces as they are being driven apart. We speculate that higher order topological invariants beyond magnetic helicity may also be responsible for slowing down the decay of the two more complicated nonhelical structures mentioned above.

preprint2011arXiv

Decay of trefoil and other magnetic knots

Two setups with interlocked magnetic flux tubes are used to study the evolution of magnetic energy and helicity on magnetohydrodynamical (MHD) systems like plasmas. In one setup the initial helicity is zero while in the other it is finite. To see if it is the actual linking or merely the helicity content that influences the dynamics of the system we also consider a setup with unlinked field lines as well as a field configuration in the shape of a trefoil knot. For helical systems the decay of magnetic energy is slowed down by the helicity which decays slowly. It turns out that it is the helicity content, rather than the actual linking, that is significant for the dynamics.

preprint2011arXiv

Detection of negative effective magnetic pressure instability in turbulence simulations

We present the first demonstration of the negative effective magnetic pressure instability in direct numerical simulations of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta. By the action of this instability, an initially uniform horizontal magnetic field forms flux concentrations whose scale is large compared to the turbulent scale. We further show that the magnetic energy of these large-scale structures is only weakly dependent on the magnetic Reynolds number, provided its value is large enough for the instability to be excited. Our results support earlier mean-field calculations and analytic work which identified this instability. Applications to the formation of active regions in the Sun are discussed.

preprint2011arXiv

Influence of Magnetic Helicity in MHD

Observations have shown that the Sun's magnetic field has helical structures. The helicity content in magnetic field configurations is a crucial constraint on the dynamical evolution of the system. Since helicity is connected with the number of links we investigate configurations with interlocked magnetic flux rings and one with unlinked rings. It turns out that it is not the linking of the tubes which affects the magnetic field decay, but the content of magnetic helicity.

preprint2011arXiv

Lagrangian coherent structures in nonlinear dynamos

Turbulence and chaos play a fundamental role in stellar convective zones through the transportof particles, energy and momentum, and in fast dynamos, through the stretching, twisting and folding of magnetic flux tubes. A particularly revealing way to describe turbulent motions is through the analysis of Lagrangian coherent structures (LCS), which are material lines or surfaces that act as transport barriers in the fluid. We report the detection of Lagrangian coherent structures in helical MHD dynamo simulations with scale separation. In an ABC--flow, two dynamo regimes, a propagating coherent mean--field regime and an intermittent regime, are identified as the magnetic diffusivity is varied. The sharp contrast between the chaotic tangle of attracting and repelling LCS in both regimes permits a unique analysis of the impact of the magnetic field on the velocity field. In a second example, LCS reveal the link between the level of chaotic mixing of the velocity field and the saturation of a large--scale dynamo when the magnetic field exceeds the equipartition value.

preprint2011arXiv

Magnetic helicity fluxes in alpha Omega dynamos

In turbulent dynamos the production of large-scale magnetic fields is accompanied by a separation of magnetic helicity in scale. The large- and small-scale parts increase in magnitude. The small-scale part can eventually work against the dynamo and quench it, especially at high magnetic Reynolds numbers. A one-dimensional mean-field model of a dynamo is presented where diffusive magnetic helicity fluxes within the domain are important. It turns out that this effect helps to alleviate the quenching. Here we show that internal magnetic helicity fluxes, even within one hemisphere, can be important for alleviating catastrophic quenching.

preprint2011arXiv

Magnetic helicity transport in the advective gauge family

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e. when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential.

preprint2011arXiv

Mean-field diffusivities in passive scalar and magnetic transport in irrotational flows

Certain aspects of the mean-field theory of turbulent passive scalar transport and of mean-field electrodynamics are considered with particular emphasis on aspects of compressible fluids. It is demonstrated that the total mean-field diffusivity for passive scalar transport in a compressible flow may well be smaller than the molecular diffusivity. This is in full analogy to an old finding regarding the magnetic mean-field diffusivity in an electrically conducting turbulently moving compressible fluid. These phenomena occur if the irrotational part of the motion dominates the vortical part, the Pèclet or magnetic Reynolds number is not too large, and, in addition, the variation of the flow pattern is slow. For both the passive scalar and the magnetic cases several further analytical results on mean-field diffusivities and related quantities found within the second-order correlation approximation are presented, as well as numerical results obtained by the test-field method, which applies independently of this approximation. Particular attention is paid to non-local and non-instantaneous connections between the turbulence-caused terms and the mean fields. Two examples of irrotational flows, in which interesting phenomena in the above sense occur, are investigated in detail. In particular, it is demonstrated that the decay of a mean scalar in a compressible fluid under the influence of these flows can be much slower than without any flow, and can be strongly influenced by the so-called memory effect, that is, the fact that the relevant mean-field coefficients depend on the decay rates themselves.

preprint2011arXiv

Nonlinear small-scale dynamos at low magnetic Prandtl numbers

Saturated small-scale dynamo solutions driven by isotropic non-helical turbulence are presented at low magnetic Prandtl numbers Pm down to 0.01. For Pm < 0.1, most of the energy is dissipated via Joule heat and, in agreement with earlier results for helical large-scale dynamos, kinetic energy dissipation is shown to diminish proportional to Pm^{1/2} down to values of 0.1. In agreement with earlier work, there is, in addition to a short Golitsyn k^{-11/3} spectrum near the resistive scale also some evidence for a short k^{-1} spectrum on larger scales. The rms magnetic field strength of the small-scale dynamo is found to depend only weakly on the value of Pm and decreases by about a factor of 2 as Pm is decreased from 1 to 0.01. The possibility of dynamo action at Pm=0.1 in the nonlinear regime is argued to be a consequence of a suppression of the bottleneck seen in the kinetic energy spectrum in the absence of a dynamo and, more generally, a suppression of kinetic energy near the dissipation wavenumber.

preprint2011arXiv

Scale-dependence of magnetic helicity in the solar wind

We determine the magnetic helicity, along with the magnetic energy, at high latitudes using data from the Ulysses mission. The data set spans the time period from 1993 to 1996. The basic assumption of the analysis is that the solar wind is homogeneous. Because the solar wind speed is high, we follow the approach first pioneered by Matthaeus et al. (1982, Phys. Rev. Lett. 48, 1256) by which, under the assumption of spatial homogeneity, one can use Fourier transforms of the magnetic field time series to construct one-dimensional spectra of the magnetic energy and magnetic helicity under the assumption that the Taylor frozen-in-flow hypothesis is valid. That is a well-satisfied assumption for the data used in this study. The magnetic helicity derives from the skew-symmetric terms of the three-dimensional magnetic correlation tensor, while the symmetric terms of the tensor are used to determine the magnetic energy spectrum. Our results show a sign change of magnetic helicity at wavenumber k~2 AU^{-1} (or frequency nu~2 uHz) at distances below 2.8 AU and at k~30 AU^{-1} (or nu~25 uHz) at larger distances. At small scales the magnetic helicity is positive at northern heliographic latitudes and negative at southern latitudes. The positive magnetic helicity at small scales is argued to be the result of turbulent diffusion reversing the sign relative to what is seen at small scales at the solar surface. Furthermore, the magnetic helicity declines toward solar minimum in 1996. The magnetic helicity flux integrated separately over one hemisphere amounts to about 10^{45} Mx^2/cycle at large scales and to a 3 times lower value at smaller scales.

preprint2011arXiv

Shear-driven instabilities in Hall-MHD plasmas

The large-scale dynamics of plasmas is well described within the framework of magnetohydrodynamics (MHD). However, whenever the ion density of the plasma becomes sufficiently low, the Hall effect is likely to become important. The role of the Hall effect has been studied in several astrophysical plasma processes, such as magnetic reconnection, magnetic dynamo, MHD turbulence or MHD instabilities. In particular, the development of small-scale instabilities is essential to understand the transport properties in a number of astrophysical plasmas. The magneto-rotational instability, which takes place in differentially rotating accretion disks embedded in relatively weak magnetic fields, is just one example. The influence of the large-scale velocity flows on small-scale instabilities is often approximated by a linear shear flow. In this paper we quantitatively study the role of the Hall effect on plasmas embedded in large-scale shear flows. More precisely, we show that a new instability develops as long as the Hall effect is present, which we therefore term as the Hall magneto-shear instability. As a particular case, we recover the so-called magneto-rotational instability and quantitatively assess the role of the Hall effect on its development and evolution.

preprint2011arXiv

Spontaneous chiral parity breaking by hydromagnetic buoyancy

Evidence for a parity-breaking nature of the magnetic buoyancy instability in a stably stratified gas is reported. In the absence of rotation, no helicity is produced, but the non-helical state is found to be unstable to small helical perturbations during the development of the instability. The parity-breaking nature of an instability in magnetohydrodynamics appears to be the first of its kind and is similar to chiral symmetry breaking in biochemistry. Applications to the production of mean fields in galaxy clusters are being discussed.

preprint2011arXiv

Spontaneous chiral symmetry breaking in the Tayler instability

The chiral symmetry breaking properties of the Tayler instability are discussed. Effective amplitude equations are determined in one case. This model has three free parameters that are determined numerically. Comparison with chiral symmetry breaking in biochemistry is made.

preprint2011arXiv

Spontaneous formation of flux concentrations in a stratified layer

The negative effective magnetic pressure instability discovered recently in direct numerical simulations (DNS) may play a crucial role in the formation of sunspots and active regions in the Sun and stars. This instability is caused by a negative contribution of turbulence to the effective mean Lorentz force (the sum of turbulent and non-turbulent contributions) and results in formation of large-scale inhomogeneous magnetic structures from initial uniform magnetic field. Earlier investigations of this instability in DNS of stably stratified, externally forced, isothermal hydromagnetic turbulence in the regime of large plasma beta are now extended into the regime of larger scale separation ratios where the number of turbulent eddies in the computational domain is about 30. Strong spontaneous formation of large-scale magnetic structures is seen even without performing any spatial averaging. These structures encompass many turbulent eddies. The characteristic time of the instability is comparable to the turbulent diffusion time, L^2/eta_t, where eta_t is the turbulent diffusivity and L is the scale of the domain. DNS are used to confirm that the effective magnetic pressure does indeed become negative for magnetic field strengths below the equipartition field. The dependence of the effective magnetic pressure on the field strength is characterized by fit parameters that seem to show convergence for larger values of the magnetic Reynolds number.

preprint2011arXiv

The fratricide of alpha-Omega dynamos by their alpha-squared siblings

Context. Helically forced magneto-hydrodynamic shearing-sheet turbulence can support different large-scale dynamo modes, although the α - Ω mode is generally expected to dominate because it is the fastest growing. In an α - Ω dynamo, most of the field amplification is produced by the shear. As differential rotation is an ubiquitous source of shear in astrophysics, such dynamos are believed to be the source of most astrophysical large-scale magnetic fields. Aims. We study the stability of oscillatory migratory α - Ω type dynamos in turbulence simulations. Methods. We use shearing-sheet simulations of hydromagnetic turbulence that is helically forced at a wavenumber that is about three times larger than the lowest wavenumber in the domain so that both α - Ω and α2 dynamo action is possible. Results. After initial dominance and saturation, the α - Ω mode is found to be destroyed by an orthogonal α2 mode sustained by the helical turbulence alone. We show that there are at least two processes through which this transition can occur. Conclusions. The fratricide of α - Ω dynamos by its α2 sibling is discussed in the context of grand minima of solar and stellar activity. However, the genesis of α - Ω dynamos from an α2 dynamo has not yet been found.

preprint2011arXiv

Vorticity from irrotationally forced flow

In the interstellar medium the turbulence is believed to be forced mostly through supernova explosions. In a first approximation these flows can be written as a gradient of a potential being thus devoid of vorticity. There are several mechanisms that could lead to vorticity generation, like viscosity and baroclinic terms, rotation, shear and magnetic fields, but it is not clear how effective they are, neither is it clear whether the vorticity is essential in determining the turbulent diffusion acting in the ISM. Here we present a study of the role of rotation, shear and baroclinicity in the generation of vorticity in the ISM.

preprint2011arXiv

Vorticity production through rotation, shear and baroclinicity

In the absence of rotation and shear, and under the assumption of constant temperature or specific entropy, purely potential forcing by localized expansion waves is known to produce irrotational flows that have no vorticity. Here we study the production of vorticity under idealized conditions when there is rotation, shear, or baroclinicity, to address the problem of vorticity generation in the interstellar medium in a systematic fashion. We use three-dimensional periodic box numerical simulations to investigate the various effects in isolation. We find that for slow rotation, vorticity production in an isothermal gas is small in the sense that the ratio of the root-mean-square values of vorticity and velocity is small compared with the wavenumber of the energy-carrying motions. For Coriolis numbers above a certain level, vorticity production saturates at a value where the aforementioned ratio becomes comparable with the wavenumber of the energy-carrying motions. Shear also raises the vorticity production, but no saturation is found. When the assumption of isothermality is dropped, there is significant vorticity production by the baroclinic term once the turbulence becomes supersonic. In galaxies, shear and rotation are estimated to be insufficient to produce significant amounts of vorticity, leaving therefore only the baroclinic term as the most favorable candidate. We also demonstrate vorticity production visually as a result of colliding shock fronts.

preprint2010arXiv

3D MHD simulations of subsurface convection in OB stars

During their main sequence evolution, massive stars can develop convective regions very close to their surface. These regions are caused by an opacity peak associated with iron ionization. Cantiello et al. (2009) found a possible connection between the presence of sub-photospheric convective motions and small scale stochastic velocities in the photosphere of early-type stars. This supports a physical mechanism where microturbulence is caused by waves that are triggered by subsurface convection zones. They further suggest that clumping in the inner parts of the winds of OB stars could be related to subsurface convection, and that the convective layers may also be responsible for stochastic excitation of non-radial pulsations. Furthermore, magnetic fields produced in the iron convection zone could appear at the surface of such massive stars. Therefore subsurface convection could be responsible for the occurrence of observable phenomena such as line profile variability and discrete absorption components. These phenomena have been observed for decades, but still evade a clear theoretical explanation. Here we present preliminary results from 3D MHD simulations of such subsurface convection.

preprint2010arXiv

Alpha effect due to buoyancy instability of a magnetic layer

A strong toroidal field can exist in form of a magnetic layer in the overshoot region below the solar convection zone. This motivates a more detailed study of the magnetic buoyancy instability with rotation. We calculate the alpha effect due to helical motions caused by a disintegrating magnetic layer in a rotating density-stratified system with angular velocity Omega making an angle theta with the vertical. We also study the dependence of the alpha effect on theta and the strength of the initial magnetic field. We carry out three-dimensional hydromagnetic simulations in Cartesian geometry. A turbulent EMF due to the correlations of the small scale velocity and magnetic field is generated. We use the test-field method to calculate the transport coefficients of the inhomogeneous turbulence produced by the layer. We show that the growth rate of the instability and the twist of the magnetic field vary monotonically with the ratio of thermal conductivity to magnetic diffusivity. The resulting alpha effect is inhomogeneous and increases with the strength of the initial magnetic field. It is thus an example of an "anti-quenched" alpha effect. The alpha effect is nonlocal, requiring around 8--16 Fourier modes to reconstruct the actual EMF based on the actual mean field.

preprint2010arXiv

Autocatalyses

Autocatalysis is a fundamental concept, used in a wide range of domains. From the most general definition of autocatalysis, that is a process in which a chemical compound is able to catalyze its own formation, several different systems can be described. We detail the different categories of autocatalyses, and compare them on the basis of their mechanistic, kinetic, and dynamic properties. It is shown how autocatalytic patterns can be generated by different systems of chemical reactions. The notion of autocatalysis covering a large variety of mechanistic realisations with very similar behaviors, it is proposed that the key signature of autocatalysis is its kinetic pattern expressed in a mathematical form.

preprint2010arXiv

Calibrating passive scalar transport in shear-flow turbulence

The turbulent diffusivity tensor is determined for linear shear flow turbulence using numerical simulations. For moderately strong shear, the diagonal components are found to increase quadratically with Peclet and Reynolds numbers below about 10 and then become constant. The diffusivity tensor is found to have components proportional to the symmetric and antisymmetric parts of the velocity gradient matrix, as well as products of these. All components decrease with the wave number of the mean field in a Lorentzian fashion. The components of the diffusivity tensor are found not to depend significantly on the presence of helicity in the turbulence. The signs of the leading terms in the expression for the diffusion tensor are found to be in good agreement with estimates based on a simple closure assumption.

preprint2010arXiv

Can catastrophic quenching be alleviated by separating shear and alpha effect?

The small-scale magnetic helicity produced as a by-product of the large-scale dynamo is believed to play a major role in dynamo saturation. In a mean-field model the generation of small-scale magnetic helicity can be modelled by using the dynamical quenching formalism. Catastrophic quenching refers to a decrease of the saturation field strength with increasing Reynolds number. It has been suggested that catastrophic quenching only affects the region of non-zero helical turbulence (i.e. where the kinematic alpha operates) and that it is possible to alleviate catastrophic quenching by separating the region of strong shear from the alpha layer. We perform a systematic study of a simple axisymmetric two-layer alpha-omega dynamo in a spherical shell for Reynolds numbers in the range 1 < Rm < 10^5. In the framework of dynamical quenching we show that this may not be the case, suggesting that magnetic helicity fluxes would be necessary.

preprint2010arXiv

Dynamo generated field emergence through recurrent plasmoid ejections

Magnetic buoyancy is believed to drive the transport of magnetic flux tubes from the convection zone to the surface of the Sun. The magnetic fields form twisted loop-like structures in the solar atmosphere. In this paper we use helical forcing to produce a large-scale dynamo-generated magnetic field, which rises even without magnetic buoyancy. A two layer system is used as computational domain where the upper part represents the solar atmosphere. Here, the evolution of the magnetic field is solved with the stress--and--relax method. Below this region a magnetic field is produced by a helical forcing function in the momentum equation, which leads to dynamo action. We find twisted magnetic fields emerging frequently to the outer layer, forming arch-like structures. In addition, recurrent plasmoid ejections can be found by looking at space--time diagrams of the magnetic field. Recent simulations in spherical coordinates show similar results.

preprint2010arXiv

How can vorticity be produced in irrotationally forced flows?

A spherical hydrodynamical expansion flow can be described as the gradient of a potential. In that case no vorticity should be produced, but several additional mechanisms can drive its production. Here we analyze the effects of baroclinicity, rotation and shear in the case of a viscous fluid. Those flows resemble what happens in the interstellar medium. In fact in this astrophysical environment supernovae explosion are the dominant flows and, in a first approximation, they can be seen as spherical. One of the main difference is that in our numerical study we examine only weakly supersonic flows, while supernovae explosions are strongly supersonic.

preprint2010arXiv

Magnetic field decay of three interlocked flux rings with zero linking number

The resistive decay of chains of three interlocked magnetic flux rings is considered. Depending on the relative orientation of the magnetic field in the three rings, the late-time decay can be either fast or slow. Thus, the qualitative degree of tangledness is less important than the actual value of the linking number or, equivalently, the net magnetic helicity. Our results do not suggest that invariants of higher order than that of the magnetic helicity need to be considered to characterize the decay of the field.

preprint2010arXiv

Magnetic helicity flux in the presence of shear

Magnetic helicity has risen to be a major player in dynamo theory, with the helicity of the small-scale field being linked to the dynamo saturation process for the large-scale field. It is a nearly conserved quantity, which allows its evolution equation to be written in terms of production and flux terms. The flux term can be decomposed in a variety of fashions. One particular contribution that has been expected to play a significant role in dynamos in the presence of mean shear was isolated by Vishniac & Cho (2001, ApJ 550, 752). Magnetic helicity fluxes are explicitly gauge dependent however, and the correlations that have come to be called the Vishniac-Cho flux were determined in the Coulomb gauge, which turns out to be fraught with complications in shearing systems. While the fluxes of small-scale helicity are explicitly gauge dependent, their divergences can be gauge independent. We use this property to investigate magnetic helicity fluxes of small-scale field through direct numerical simulations in a shearing-box system and find that in a numerically usable gauge the divergence of the small-scale helicity flux vanishes, while the divergence of the Vishniac-Cho flux remains finite. We attribute this seeming contradiction to the existence of horizontal fluxes of small-scale magnetic helicity with finite divergences even in our shearing-periodic domain.

preprint2010arXiv

Magnetic helicity fluxes in interface and flux transport dynamos

Dynamos in the Sun and other bodies tend to produce magnetic fields that possess magnetic helicity of opposite sign at large and small scales, respectively. The build-up of magnetic helicity at small scales provides an important saturation mechanism. In order to understand the nature of the solar dynamo we need to understand the details of the saturation mechanism in spherical geometry. In particular, we want to understand the effects of magnetic helicity fluxes from turbulence and meridional circulation. We consider a model with just radial shear confined to a thin layer (tachocline) at the bottom of the convection zone. The kinetic alpha owing to helical turbulence is assumed to be localized in a region above the convection zone. The dynamical quenching formalism is used to describe the build-up of mean magnetic helicity in the model, which results in a magnetic alpha effect that feeds back on the kinetic alpha effect. In some cases we compare with results obtained using a simple algebraic alpha quenching formula. In agreement with earlier findings, the magnetic alpha effect in the dynamical alpha quenching formalism has the opposite sign compared with the kinetic alpha effect and leads to a catastrophic decrease of the saturation field strength with increasing magnetic Reynolds numbers. However, at high latitudes this quenching effect can lead to secondary dynamo waves that propagate poleward due to the opposite sign of alpha. Magnetic helicity fluxes both from turbulent mixing and from meridional circulation alleviate catastrophic quenching.

preprint2010arXiv

Numerical simulations of the decay of primordial magnetic turbulence

We perform direct numerical simulations of forced and freely decaying 3D magnetohydrodynamic turbulence in order to model magnetic field evolution during cosmological phase transitions in the early Universe. Our approach assumes the existence of a magnetic field generated either by a process during inflation or shortly thereafter, or by bubble collisions during a phase transition. We show that the final configuration of the magnetic field depends on the initial conditions, while the velocity field is nearly independent of initial conditions.

preprint2010arXiv

Oscillatory migrating magnetic fields in helical turbulence in spherical domains

We present direct numerical simulations of the equations of compressible magnetohydrodynamics in a wedge-shaped spherical shell, without shear, but with random helical forcing which has negative (positive) helicity in the northern (southern) hemisphere. We find a large-scale magnetic field that is nearly uniform in the azimuthal direction and approximately antisymmetric about the equator. Furthermore, the large-scale field in each hemisphere oscillates on nearly dynamical time scales with reversals of polarity and equatorward migration. Corresponding mean-field models also show similar migratory oscillations with a frequency that is nearly independent of the magnetic Reynolds number. This mechanism may be relevant for understanding equatorward migration seen in the solar dynamo.

preprint2010arXiv

Plasmoid ejections driven by dynamo action underneath a spherical surface

We present a unified three-dimensional model of the convection zone and upper atmosphere of the Sun in spherical geometry. In this model, magnetic fields, generated by a helically forced dynamo in the convection zone, emerge without the assistance of magnetic buoyancy. We use an isothermal equation of state with gravity and density stratification. Recurrent plasmoid ejections, which rise through the outer atmosphere, is observed. In addition, the current helicity of the small--scale field is transported outwards and form large structures like magnetic clouds.

preprint2010arXiv

Recurrent flux emergence from dynamo-generated fields

We investigate the emergence of a large-scale magnetic field. This field is dynamo-generated by turbulence driven with a helical forcing function. Twisted arcade-like field structures are found to emerge in the exterior above the turbulence zone. Time series of the magnetic field structure show recurrent plasmoid ejections.

preprint2010arXiv

Shear-driven and diffusive helicity fluxes in alpha-Omega dynamos

We present nonlinear mean-field alpha-Omega dynamo simulations in spherical geometry with simplified profiles of kinematic alpha effect and shear. We take magnetic helicity evolution into account by solving a dynamical equation for the magnetic alpha effect. This gives a consistent description of the quenching mechanism in mean-field dynamo models. The main goal of this work is to explore the effects of this quenching mechanism in solar-like geometry, and in particular to investigate the role of magnetic helicity fluxes, specifically diffusive and Vishniac-Cho (VC) fluxes, at large magnetic Reynolds numbers (Rm). For models with negative radial shear or positive latitudinal shear, the magnetic alpha effect has predominantly negative (positive) sign in the northern (southern) hemisphere. In the absence of fluxes, we find that the magnetic energy follows an Rm^-1 dependence, as found in previous works. This catastrophic quenching is alleviated in models with diffusive magnetic helicity fluxes resulting in magnetic fields comparable to the equipartition value even for Rm=10^7. On the other hand, models with a shear-driven Vishniac-Cho flux show an increase of the amplitude of the magnetic field with respect to models without fluxes, but only for Rm<10^4. This is mainly a consequence of assuming a vacuum outside the Sun which cannot support a significant VC flux across the boundary. However, in contrast with the diffusive flux, the VC flux modifies the distribution of the magnetic field. In addition, if an ill-determined scaling factor in the expression for the VC flux is large enough, subcritical dynamo action is possible that is driven by the action of shear and the divergence of current helicity flux.

preprint2010arXiv

Simulations of astrophysical dynamos

Numerical aspects of dynamos in periodic domains are discussed. Modifications of the solutions by numerically motivated alterations of the equations are being reviewed using the examples of magnetic hyperdiffusion and artificial diffusion when advancing the magnetic field in its Euler potential representation. The importance of using integral kernel formulations in mean-field dynamo theory is emphasized in cases where the dynamo growth rate becomes comparable with the inverse turnover time. Finally, the significance of microscopic magnetic Prandtl number in controlling the conversion from kinetic to magnetic energy is highlighted.

preprint2010arXiv

Turbulence and magnetic spots at the surface of hot massive stars

Hot luminous stars show a variety of phenomena in their photospheres and in their winds which still lack clear physical explanations at this time. Among these phenomena are non-thermal line broadening, line profile variability (LPVs), discrete absorption components (DACs), wind clumping and stochastically excited pulsations. Cantiello et al. (2009) argued that a convection zone close to the surface of hot, massive stars, could be responsible for some of these phenomena. This convective zone is caused by a peak in the opacity due to iron recombination and for this reason is referred as the "iron convection zone" (FeCZ). 3D MHD simulations are used to explore the possible effects of such subsurface convection on the surface properties of hot, massive stars. We argue that turbulence and localized magnetic spots at the surface are the likely consequence of subsurface convection in early type stars.

preprint2010arXiv

Turbulent magnetic pressure instability in stratified turbulence

A reduction of total mean turbulent pressure due to the presence of magnetic fields was previously shown to be a measurable effect in direct numerical simulations. However, in the studied parameter regime the formation of large-scale structures, as anticipated from earlier mean-field simulations, was not found. An analysis of the relevant mean-field parameter dependency and the parameter domain of interest is conducted in order to clarify this apparent discrepancy.

preprint2009arXiv

Equatorial magnetic helicity flux in simulations with different gauges

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small-scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo-Lorenz gauge, where the speed of light is replaced by the sound speed, and the `resistive gauge' in which the Laplacian of the magnetic vector potential acts as resistive term. We find that, in the statistically steady state, the time-averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges.

preprint2009arXiv

From fibril to diffuse fields during dynamo saturation

The degree of intermittency of the magnetic field of a large-scale dynamo is considered. Based on simulations it is argued that there is a tendency for the field to become more diffuse and non-intermittent as the dynamo saturates. The simulations are idealized in that the turbulence is strongly helical and shear is strong, so the tendency for the field to become more diffuse is somewhat exaggerated. Earlier results concerning the effects of magnetic buoyancy are discussed. It is emphasized that the resulting magnetic buoyancy is weak compared with the stronger effects of simultaneous downward pumping. These findings are used to support the notion that the solar dynamo might operate in a distributed fashion where the near-surface shear layer could play an important role.

preprint2009arXiv

Homochirality and the need of energy

The mechanisms for explaining how a stable asymmetric chemical system can be formed from a symmetric chemical system, in the absence of any asymmetric influence other than statistical fluctuations, have been developed during the last decades, focusing on the non-linear kinetic aspects. Besides the absolute necessity of self-amplification processes, the importance of energetic aspects is often underestimated. Going down to the most fundamental aspects, the distinction between a single object -- that can be intrinsically asymmetric -- and a collection of objects -- whose racemic state is the more stable one -- must be emphasized. A system of strongly interacting objects can be described as one single object retaining its individuality and a single asymmetry; weakly or non-interacting objects keep their own individuality, and are prone to racemize towards the equilibrium state. In the presence of energy fluxes, systems can be maintained in an asymmetric non-equilibrium steady-state. Such dynamical systems can retain their asymmetry for times longer than their racemization time.

preprint2009arXiv

Large-scale dynamos at low magnetic Prandtl numbers

Using direct simulations of hydromagnetic turbulence driven by random polarized waves it is shown that dynamo action is possible over a wide range of magnetic Prandtl numbers from 10^-3 to 1. Triply periodic boundary conditions are being used. In the final saturated state the resulting magnetic field has a large-scale component of Beltrami type. For the kinematic phase, growth rates have been determined for magnetic Prandtl numbers between 0.01 and 1, but only the case with the smallest magnetic Prandtl number shows large-scale magnetic fields. It is less organized than in the nonlinear stage. For small magnetic Prandtl numbers the growth rates are comparable to those calculated from an alpha squared mean-field dynamo. In the linear regime the magnetic helicity spectrum has a short inertial range compatible with a -5/3 power law, while in the nonlinear regime it is the current helicity whose spectrum may be compatible with such a law. In the saturated case, the spectral magnetic energy in the inertial range is in slight excess over the spectral kinetic energy, although for small magnetic Prandtl numbers the magnetic energy spectrum reaches its resistive cut off wavenumber more quickly. The viscous energy dissipation declines with the square root of the magnetic Prandtl number, which implies that most of the energy is dissipated via Joule heat.

preprint2009arXiv

Large-scale magnetic flux concentrations from turbulent stresses

In this study we provide the first numerical demonstration of the effects of turbulence on the mean Lorentz force and the resulting formation of large-scale magnetic structures. Using three-dimensional direct numerical simulations (DNS) of forced turbulence we show that an imposed mean magnetic field leads to a decrease of the turbulent hydromagnetic pressure and tension. This phenomenon is quantified by determining the relevant functions that relate the sum of the turbulent Reynolds and Maxwell stresses with the Maxwell stress of the mean magnetic field. Using such a parameterization, we show by means of two-dimensional and three-dimensional mean-field numerical modelling that an isentropic density stratified layer becomes unstable in the presence of a uniform imposed magnetic field. This large-scale instability results in the formation of loop-like magnetic structures which are concentrated at the top of the stratified layer. In three dimensions these structures resemble the appearance of bipolar magnetic regions in the Sun. The results of DNS and mean-field numerical modelling are in good agreement with theoretical predictions. We discuss our model in the context of a distributed solar dynamo where active regions and sunspots might be rather shallow phenomena.

preprint2009arXiv

Shear-driven magnetic buoyancy oscillations

The effects of uniform horizontal shear on a stably stratified layer of gas is studied. The system is initially destabilized by a magnetically buoyant flux tube pointing in the cross-stream direction. The shear amplifies the initial field to Lundquist numbers of about 200-400, but then its value drops to about 100-300, depending on the value of the sub-adiabatic gradient. The larger values correspond to cases where the stratification is strongly stable and nearly isothermal. At the end of the runs the magnetic field is nearly axisymmetric, i.e. uniform in the streamwise direction. In view of Cowling's theorem the sustainment of the field remains a puzzle and may be due to subtle numerical effects that have not yet been identified in detail. In the final state the strength of the magnetic field decreases with height in such a way that the field is expected to be unstable. Low amplitude oscillations are seen in the vertical velocity even at late times, suggesting that they might be persistent.

preprint2009arXiv

Turbulent diffusion and galactic magnetism

Using the test-field method for nearly irrotational turbulence driven by spherical expansion waves it is shown that the turbulent magnetic diffusivity increases with magnetic Reynolds numbers. Its value levels off at several times the rms velocity of the turbulence multiplied by the typical radius of the expansion waves. This result is discussed in the context of the galactic mean-field dynamo.

preprint2009arXiv

Turbulent dynamos in spherical shell segments of varying geometrical extent

We use three-dimensional direct numerical simulations of the helically forced magnetohydrodynamic equations in spherical shell segments in order to study the effects of changes in the geometrical shape and size of the domain on the growth and saturation of large-scale magnetic fields. We inject kinetic energy along with kinetic helicity in spherical domains via helical forcing using Chandrasekhar-Kendall functions. We take perfect conductor boundary conditions for the magnetic field to ensure that no magnetic helicity escapes the domain boundaries. We find dynamo action giving rise to magnetic fields at scales larger than the characteristic scale of the forcing. The magnetic energy exceeds the kinetic energy over dissipative time scales, similar to that seen earlier in Cartesian simulations in periodic boxes. As we increase the size of the domain in the azimuthal direction we find that the nonlinearly saturated magnetic field organizes itself in long-lived cellular structures with aspect ratios close to unity. These structures tile the domain along the azimuthal direction, thus resulting in very small longitudinally averaged magnetic fields for large domain sizes. The scales of these structures are determined by the smallest scales of the domain, which in our simulations is usually the radial scale. We also find that increasing the meridional extent of the domains produces little qualitative change, except a marginal increase in the large-scale field. We obtain qualitatively similar results in Cartesian domains with similar aspect ratios.

preprint2008arXiv

Mean-field effects in the Galloway-Proctor flow

The coefficients defining the mean electromotive force in a Galloway-Proctor flow are determined. This flow shows a two-dimensional pattern and is helical. The pattern wobbles in its plane. Apart from one exception a circular motion of the flow pattern is assumed. This corresponds to one of the cases considered recently by Courvoisier, Hughes and Tobias (2006, Phys. Rev. Lett., 96, 034503). An analytic theory of the alpha effect and related effects in this flow is developed within the second-order correlation approximation and a corresponding fourth-order approximation. In the validity range of these approximations there is an alpha effect but no gamma effect, or pumping effect. Numerical results obtained with the test-field method, which are independent of these approximations, confirm the results for alpha and show that gamma is in general nonzero. Both alpha and gamma show a complex dependency on the magnetic Reynolds number and other parameters that define the flow, that is, amplitude and frequency of the wobbling motion. Some results for the magnetic diffusivity eta_t and a related quantity are given, too. Finally a result for alpha in the case of a randomly varying flow without the aforementioned circular motion is presented. This flow may be a more appropriate model for studying the alpha effect and related effects in flows that are statistical isotropic in a plane.

preprint2008arXiv

Nonlinear aspects of astrobiological research

Several aspects of mathematical astrobiology are discussed. It is argued that around the time of the origin of life the handedness of biomolecules must have established itself through an instability. Possible pathways of producing a certain handedness include mechanisms involving either autocatalysis or, alternatively, epimerization as governing effects. Concepts for establishing hereditary information are discussed in terms of the theory of hypercycles. Instabilities toward parasites and possible remedies by invoking spatial extent are reviewed. Finally, some effects of early life are discussed that contributed to modifying and regulating atmosphere and climate of the Earth, and that could have contributed to the highly oxidized state of its crust.

preprint2006arXiv

Magnetic helicity density and its flux in weakly inhomogeneous turbulence

A gauge invariant and hence physically meaningful definition of magnetic helicity density for random fields is proposed, using the Gauss linking formula, as the density of correlated field line linkages. This definition is applied to the random small scale field in weakly inhomogeneous turbulence, whose correlation length is small compared with the scale on which the turbulence varies. For inhomogeneous systems, with or without boundaries, our technique then allows one to study the local magnetic helicity density evolution in a gauge independent fashion, which was not possible earlier. This evolution equation is governed by local sources (owing to the mean field) and by the divergence of a magnetic helicity flux density. The role of magnetic helicity fluxes in alleviating catastrophic quenching of mean field dynamos is discussed.

preprint2005arXiv

Spatial dynamics of homochiralization

The emergence and spreading of chirality on the early Earth is considered by studying a set of reaction-diffusion equations based on a polymerization model. It is found that effective mixing of the early oceans is necessary to reach the present homochiral state. The possibility of introducing mass extinctions and modifying the emergence rate of life is discussed.

preprint2001arXiv

Solitons in the noisy Burgers equation

We investigate numerically the coupled diffusion-advective type field equations originating from the canonical phase space approach to the noisy Burgers equation or the equivalent Kardar-Parisi-Zhang equation in one spatial dimension. The equations support stable right hand and left hand solitons and in the low viscosity limit a long-lived soliton pair excitation. We find that two identical pair excitations scatter transparently subject to a size dependent phase shift and that identical solitons scatter on a static soliton transparently without a phase shift. The soliton pair excitation and the scattering configurations are interpreted in terms of growing step and nucleation events in the interface growth profile. In the asymmetrical case the soliton scattering modes are unstable presumably toward multi soliton production and extended diffusive modes, signalling the general non-integrability of the coupled field equations. Finally, we have shown that growing steps perform anomalous random walk with dynamic exponent z=3/2 and that the nucleation of a tip is stochastically suppressed with respect to plateau formation.

preprint2000arXiv

The inverse cascade and nonlinear alpha-effect in simulations of isotropic helical hydromagnetic turbulence

A numerical model of isotropic homogeneous turbulence with helical forcing is investigated. The resulting flow, which is essentially the prototype of the alpha^2 dynamo of mean-field dynamo theory, produces strong dynamo action with an additional large scale field on the scale of the box (at wavenumber k=1; forcing is at k=5). This large scale field is nearly force-free and exceeds the equipartition value. As the magnetic Reynolds number R_m increases, the saturation field strength and the growth rate of the dynamo increase. However, the time it takes to built up the large scale field from equipartition to its final super-equipartition value increases with magnetic Reynolds number. The large scale field generation can be identified as being due to nonlocal interactions originating from the forcing scale, which is characteristic of the alpha-effect. Both alpha and turbulent magnetic diffusivity eta_t are determined simultaneously using numerical experiments where the mean-field is modified artificially. Both quantities are quenched in a R_m-dependent fashion. The evolution of the energy of the mean field matches that predicted by an alpha^2 dynamo model with similar alpha and eta_t quenchings. For this model an analytic solution is given which matches the results of the simulations. The simulations are numerically robust in that the shape of the spectrum at large scales is unchanged when changing the resolution from 30^3 to 120^3 meshpoints, or when increasing the magnetic Prandtl number (viscosity/magnetic diffusivity) from 1 to 100. Increasing the forcing wavenumber to 30 (i.e. increasing the scale separation) makes the inverse cascade effect more pronounced, although it remains otherwise qualitatively unchanged.

preprint2000arXiv

The response of a turbulent accretion disc to an imposed epicyclic shearing motion

We excite an epicyclic motion, whose amplitude depends on the vertical position, $z$, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays we can obtain information about the interaction between the warp and the disc turbulence. A high amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, $α_{\rm v}$, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, $α_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $α_{\rm v}$ is approximately equal to $α_{\rm h}$ and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness ($\sim 0.01$) of $α_{\rm v}$ and $α_{\rm h}$ we conclude that for $β= p_{\rm gas}/p_{\rm mag} \sim 10$ the timescale for diffusion or damping of a warp is much shorter than the usual viscous timescale. Finally, we review the astrophysical implications.

preprint1997arXiv

Robustness of truncated alpha-Omega dynamos with a dynamic alpha

In a recent work (Covas et al. 1997), the behaviour and the robustness of truncated $αΩ$ dynamos with a dynamic $α$ were studied with respect to a number of changes in the driving term of the dynamic $α$ equation, which was considered previously by Schmalz and Stix (1991) to be of the form $\sim A_ϕB_ϕ$. Here we review and extend our previous work and consider the effect of adding a quadratic quenching term of the form $α|B|^2$. We find that, as before, such a change can have significant effects on the dynamics of the related truncated systems. We also find intervals of (negative) dynamo numbers, in the system considered by Schmalz and Stix (1991), for which there is sensitivity with respect to small changes in the dynamo number and the initial conditions, similar to what was found in our previous work. This latter behaviour may be of importance in producing the intermittent type of behaviour observed in the sun.

preprint1996arXiv

Large-scale magnetic fields from hydromagnetic turbulence in the very early universe

We investigate hydromagnetic turbulence of primordial magnetic fields using magnetohydrodynamics (MHD) in an expanding universe. We present the basic, covariant MHD equations, find solutions for MHD waves in the early universe, and investigate the equations numerically for random magnetic fields in two spatial dimensions. We find the formation of magnetic structures at larger and larger scales as time goes on. In three dimensions we use a cascade (shell) model, that has been rather successful in the study of certain aspects of hydrodynamic turbulence. Using such a model we find that after ${\cal O}(10^9)$ times the initial time the scale of the magnetic field fluctuation (in the comoving frame) has increased by 4-5 orders of magnitude as a consequence of an inverse cascade effect (i.e. transfer of energy from smaller to larger scales). Thus {\it at large scales} primordial magnetic fields are considerably stronger than expected from considerations which do not take into account the effects of MHD turbulence.

preprint1995arXiv

An Incoherent $α-Ω$ Dynamo in Accretion Disks

We use the mean-field dynamo equations to show that an incoherent alpha effect in mirror-symmetric turbulence in a shearing flow can generate a large scale, coherent magnetic field. We illustrate this effect with simulations of a few simple systems. In accretion disks, this process can lead to axisymmetric magnetic domains whose radial and vertical dimensions will be comparable to the disk height. This process may be responsible for observations of dynamo activity seen in simulations of dynamo-generated turbulence involving, for example, the Balbus-Hawley instability. In this case the magnetic field strength will saturate at $\sim (h/r)^2$ times the ambient pressure in real accretion disks. The resultant dimensionless viscosity will be of the same order. In numerical simulations the azimuthal extent of the simulated annulus should be substituted for $r$. We compare the predictions of this model to numerical simulations previously reported by Brandenburg et al. (1995). In a radiation pressure dominated environment this estimate for viscosity should be reduced by a factor of $(P_{gas}/P_{radiation})^6$ due to magnetic buoyancy.

Axel Brandenburg

What is connected

Connect this record

See the researcher in context

Building this map preview

154 published item(s)

Big Bang Nucleosynthesis Limits and Relic Gravitational Waves Detection Prospects

Chirality in Astrophysics

Collision fluctuations of lucky droplets with superdroplets

Compressible test-field method and its application to shear dynamos

Dynamo effect in unstirred self-gravitating turbulence

Dynamo instabilities in plasmas with inhomogeneous chiral chemical potential

Evolution of primordial magnetic fields during large-scale structure formation

Polarization of gravitational waves from helical MHD turbulent sources

Production of a chiral magnetic anomaly with emerging turbulence and mean-field dynamo action

Circular Polarization of Gravitational Waves from Early-Universe Helical Turbulence

On the existence of shear-current effects in magnetized burgulence

Application of a helicity proxy to edge-on galaxies

Condensational and collisional growth of cloud droplets in a turbulent environment

Dynamo in weakly collisional nonmagnetized plasmas impeded by Landau damping of magnetic fields

Hall cascade with fractional magnetic helicity in neutron star crusts

Hemispheric handedness in the Galactic synchrotron polarization foreground

Magnetic helicity dissipation and production in an ideal MHD code

On the measurement of handedness in Fermi Large Area Telescope data

Piecewise quadratic growth during the 2019 novel coronavirus epidemic

Primordial magnetic helicity evolution with a homogeneous magnetic field from inflation

The turbulent stress spectrum in the inertial and subinertial ranges

Turbulent viscosity and effective magnetic Prandtl number from simulations of isotropically forced turbulence

The nature of mean-field generation in three classes of optimal dynamos

Energetics of turbulence generated by chiral MHD dynamos

Magnetism in the Early Universe

A unified large/small-scale dynamo in helical turbulence

Bipolar region formation in stratified two-layer turbulence

Evolution of Magnetic Helicity and Energy Spectra of Solar Active Regions

High-wavenumber solar $f$-mode strengthening prior to active region formation

Influence of a coronal envelope as a free boundary to global convective dynamo simulations

Large-scale flow generation by inhomogeneous helicity

Magnetic concentrations in stratified turbulence: the negative effective magnetic pressure instability

Multiple dynamo modes as a mechanism for long-term solar activity variations

Stellar mixing length theory with entropy rain

The evolution of primordial magnetic field since its generation

Turbulent reconnection of magnetic bipoles in stratified turbulence

VIEWPOINT: A New Twist in Simulating Solar Flares

Bipolar magnetic spots from dynamos in stratified spherical shell turbulence

Hysteresis between distinct modes of turbulent dynamos

Is the small-scale magnetic field correlated with the dynamo cycle?

Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics

Magnetically controlled stellar differential rotation near the transition from solar to anti-solar profiles

Negative magnetic eddy diffusivities from test-field method and multiscale stability theory

Properties of $p$- and $f$-modes in hydromagnetic turbulence

Evolution of magnetic field generated by the Kelvin-Helmholtz instability

Fanning out of the $f$-mode in presence of nonuniform magnetic fields

Faraday signature of magnetic helicity from reduced depolarization

Hydraulic effects in a radiative atmosphere with ionization

Magnetic flux concentrations from dynamo-generated fields

Magnetic helicity and energy spectra of a solar active region

Magnetic Prandtl number dependence of kinetic to magnetic dissipation ratio

Mean-field dynamo action from delayed transport

Near-polytropic simulations with a radiative surface

Nonhelical inverse transfer of a decaying turbulent magnetic field

On the cause of solar-like equatorward migration in global convective dynamo simulations

Particle energization through time-periodic helical magnetic fields

Quenching and anisotropy of hydromagnetic turbulent transport

Simulations of galactic dynamos

Superflare occurrence and energies on G-, K- and M-type dwarfs

Surface flux concentrations and spherical alpha-square dynamo

Traces of large-scale dynamo action in the kinematic stage

A mean field dynamo from negative eddy diffusivity

Azimuthal dynamo wave in spherical shell convection

Bipolar magnetic structures driven by stratified turbulence with a coronal envelope

Can planetesimals form by collisional fusion?

Coronal influence on dynamos

Cyclic magnetic activity due to turbulent convection in spherical wedge geometry

Effects of enhanced stratification on equatorward dynamo wave propagation

Evolution of Primordial Magnetic Fields from Phase Transitions

Non-linear and chaotic dynamo regimes

Scaling and intermittency in incoherent α-shear dynamo

Self-assembly of shallow magnetic spots through strongly stratified turbulence

Spoke-like differential rotation in a convective dynamo with a coronal envelope

The kinetic helicity needed to drive large-scale dynamos