
We show that in decaying hydromagnetic turbulence with initial kinetic
helicity, a weak magnetic field eventually becomes fully helical. The sign of
magnetic helicity is opposite to that of the kinetic helicity  regardless of
whether or not the initial magnetic field was helical. The magnetic field
undergoes inverse cascading with the magnetic energy decaying approximately
like t^{1/2}. This is even slower than in the fully helical case, where it
decays like t^{2/3}. In this parameter range, the product of magnetic energy
and correlation length raised to a certain power slightly larger than unity, is
approximately constant. This scaling of magnetic energy persists over long time
scales. At very late times and for domain sizes large enough to accommodate the
growing spatial scales, we expect a crossover to the t^{2/3} decay law that
is commonly observed for fully helical magnetic fields. Regardless of the
presence or absence of initial kinetic helicity, the magnetic field experiences
exponential growth during the first few turnover times, which is suggestive of
smallscale dynamo action. Our results have applications to a wide range of
experimental dynamos and astrophysical timedependent plasmas, including
primordial turbulence in the early universe.

The relation between magnetic and velocity fields in the solar atmosphere is
a topic that can now be addressed quantitatively at the statistical level. Here
we analyze kinetic energy and cross helicity spectra along with magnetic energy
and helicity spectra in two solar active regions using vector magnetograms and
Dopplergrams of NOAA~11158 and 12266. Within these active regions, we find
similar slopes for kinetic and magnetic energy spectra at intermediate
wavenumbers, where the contribution from the granulation velocity has been
removed. At wavenumbers around $0.3\Mm^{1}$, the magnetic helicity is found to
be close to its maximal value. The cross helicity spectra are found to be
within about 10\% of the maximum possible value. We also develop a twoscale
method for cross helicity spectra, which allows us to take the cancellation
from the bipolarity of active regions into account. In the quiet Sun, by
comparison, the cross helicity spectrum is found to be small.

Using a recently developed twoscale formalism to determine the magnetic
helicity spectrum (Brandenburg et al. 2017), we analyze synoptic vector
magnetograms built with data from the Vector Spectromagnetograph (VSM)
instrument on the \emph{Synoptic Optical Longterm Investigations of the Sun}
(SOLIS) telescope during January 2010July 2016. In contrast to an earlier
study using only three Carrington rotations, our analysis includes 74 synoptic
Carrington rotation maps. We recover here bihelical spectra at different phases
of solar cycle~24, where the net magnetic helicity in the majority of the data
is consistent with a largescale dynamo with helical turbulence operating in
the Sun. More than $20\%$ of the analyzed maps, however, show violations of the
expected sign rule.

We investigate the statistical properties of isotropic, stochastic, Gaussian
distributed, helical magnetic fields characterized by different shapes of the
energy spectra at large length scales and study the associated realizability
condition. We discuss smoothed magnetic fields that are commonly used when the
primordial magnetic field is constrained by observational data. We are
particularly interested in scaleinvariant magnetic fields that can be
generated during the inflationary stage by quantum fluctuations. We determine
the correlation length of such magnetic fields and relate it to the infrared
cutoff of perturbations produced during inflation. We show that this scale
determines the observational signatures of the inflationary magnetic fields on
the cosmic microwave background. At smaller scales, the scaleinvariant
spectrum changes with time. It becomes a steeper weakturbulence spectrum at
progressively larger scales. We show numerically that the critical length scale
where this happens is the turbulentdiffusive scale, which increases with the
square root of time.

Using direct numerical simulations (DNS), we study laminar and turbulent
dynamos in chiral magnetohydrodynamics (MHD) with an extended set of equations
that accounts for an additional contribution to the electric current due to the
chiral magnetic effect (CME). This quantum phenomenon originates from an
asymmetry between left and righthanded relativistic fermions in the presence
of a magnetic field and gives rise to a chiral dynamo. We show that the
magnetic field evolution proceeds in three stages: (1) a smallscale chiral
dynamo instability; (2) production of chiral magnetically driven turbulence and
excitation of a largescale dynamo instability due to a new chiral effect
(alpha_mu effect); and (3) saturation of magnetic helicity and magnetic field
growth controlled by a conservation law for the total chirality. The alpha_mu
effect becomes dominant at large fluid and magnetic Reynolds numbers and is not
related to kinetic helicity. The growth rate of the largescale magnetic field
and its characteristic scale measured in the numerical simulations agree well
with theoretical predictions based on meanfield theory. The previously
discussed twostage chiral magnetic scenario did not include stage (2) during
which the characteristic scale of magnetic field variations can increase by
many orders of magnitude. Based on the findings from numerical simulations, the
relevance of the CME and the chiral effects revealed in the relativistic plasma
of the early Universe and of protoneutron stars are discussed.

An asymmetry in the number density of left and righthanded fermions is
known to give rise to a new term in the induction equation that can result in a
smallscale instability. This is a microphysical effect characterised by a
chiral chemical potential and is mathematically similar to the alpha effect,
which is a turbulent or macrophysical effect. At high temperatures, when a
chiral asymmetry can survive for long enough, these chiral MHD dynamos can
amplify magnetic fields efficiently, which in turn drive turbulence via the
Lorentz force. While it has been demonstrated in numerical simulations that
chiral magnetically driven turbulence exists and modifies the evolution of the
plasma, the details of this process remain unclear. The goal of this paper is
to shed new light on the properties of chiral magnetically driven turbulence
using numerical simulations with the Pencil Code. We explore the generation of
turbulence 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^2,
that can be reached in chiral magnetically driven turbulence. Within the
parameter space explored in this study, Upsilon reaches a value of
approximately 0.240.29  independently of the initial chiral asymmetry and for
Pm=1. Our simulations suggest, that Upsilon decreases as a power law when
increasing Pm. While the exact scaling depends on the details of the fitting
criteria and the Reynolds number regime, an approximate result of
Upsilon(Pm)=0.3*Pm^{0.2} is reported. Using the findings from our numerical
simulations, we estimate the properties of chiral magnetically driven
turbulence in the early Universe.

High precision photometry of solarlike members of the open cluster M67 with
Kepler/K2 data has recently revealed enhanced activity for stars with a large
Rossby number, which is the ratio of rotation period to the convective turnover
time. Contrary to the well established behavior for shorter rotation periods
and smaller Rossby numbers, the chromospheric activity of the more slowly
rotating stars of M67 was found to increase with increasing Rossby number. Such
behavior has never been reported before, although it was theoretically
predicted to emerge as a consequence of antisolar differential rotation (DR)
for stars with Rossby numbers larger than that of the Sun, because in those
models the absolute value of the DR was found to exceed that for solarlike DR.
Using gyrochronological relations and an approximate age of 4 Gyr for the
members of M67, we compare with computed rotation rates using just the BV
color. The resulting rotationactivity relation is found to be compatible with
that obtained by employing the measured rotation rate. This provides additional
support for the unconventional enhancement of activity at comparatively low
rotation rates and the possible presence of antisolar differential rotation.

Recent advances in meanfield theory are reviewed and applications to the
Sun, latetype stars, accretion disks, galaxies, and the early Universe are
discussed. We focus particularly on aspects of spatiotemporal nonlocality,
which is one of the main insights that emerged from applying the testfield
method to magnetic fields of different length and timescales. We also review
the status of nonlinear quenching and the relation to magnetic helicity, which
is an important observational diagnostic of modern solar dynamo theory. Both
solar and some stellar dynamos seem to operate in an intermediate regime that
has not yet been possible to model successfully. This regime is bracketed by
antisolarlike differential rotation on one end and stellar activity cycles
belonging to the superactive stars on the other. The difficulty in modeling
this regime may be related to shortcomings in modelling solar/stellar
convection. On galactic and extragalactic length scales, the observational
constraints are still less stringent and uncertain, but recent advances both in
theory and in observations suggest that more conclusive comparisons may soon be
possible. The possibility of inversely cascading magnetic helicity throughout
all of the early Universe is particularly exciting in explaining the lower
limits of magnetic fields on cosmological length scales and the possibility of
parity breaking and finite helicity of such a field.

We consider a scaleinvariant helical magnetic field generated during
inflation. We show that, if the mean magnetic helicity density of such a field
is measured, it can be used to determine a lower bound on the duration of
inflation. Upper bounds can be used to derive constraints on the minimal
duration of inflation if one assumes that the magnetic field generated during
inflation is helical. Using threedimensional simulations, we show that an
initially scaleinvariant field develops, which is similar both with and
without magnetic helicity. In the fully helical case, however, the magnetic
field appears to have a more pronounced folded structure.

We present new simulations of decaying hydromagnetic turbulence for a
relativistic equation of state relevant to the early universe. We compare
helical and nonhelical cases either with kinetically or magnetically dominated
initial fields. Both kinetic and magnetic initial helicities lead to maximally
helical magnetic fields after some time, but with different temporal decay
laws. Both are relevant to the early universe, although no mechanisms have yet
been identified that produce magnetic helicity with strengths comparable to the
big bang nucleosynthesis limit at scales comparable to the Hubble horizon at
the electroweak phase transition. Nonhelical magnetically dominated fields
could still produce picoGauss magnetic fields under most optimistic conditions.
Only helical magnetic fields can potentially have nanoGauss strengths at scales
up to 30 kpc today.

The formation of sunspots requires the concentration of magnetic flux near
the surface. The negative magnetic pressure instability (NEMPI) might be a
possible mechanism for accomplishing this, but it has mainly been studied in
simple systems using an isothermal equation of state without a natural free
surface. We study NEMPI in a stratified Cartesian meanfield model where
turbulence effects are parameterized. We use an ideal equation of state and
include radiation transport, which establishes selfconsistently a free surface.
We use a Kramerstype opacity with adjustable exponents chosen such that the
deeper layers are approximately isentropic. No convection is therefore possible
in this model, allowing us to study NEMPI with radiation in isolation. We
restrict ourselves to twodimensional models. We use artificially enhanced
meanfield coefficients to allow NEMPI to develop, making it therefore possible
to study the reason why it is much harder to excite in the presence of
radiation. NEMPI yields moderately strong magnetic flux concentrations a
certain distance beneath the surface where optical depth is unity. The
instability is oscillatory and in the form of upward travelling waves. This
seems to be a new effect that has not been found in earlier models without
radiative transport. The horizontal wavelength is about ten times smaller than
what has been found previously in more idealized isothermal models. In our
models, NEMPI saturates at field strengths too low to explain sunspots.
Furthermore, the structures appear too narrow and too far beneath the surface
to cause significant brightness variations at the radiative surface. We
speculate that the failure to reproduce effects resembling sunspots may be
related to the neglect of convection.

We present numerical simulations of hydrodynamic overshooting convection in
local Cartesian domains. We find that a substantial fraction of the lower part
of the convection zone (CZ) is stably stratified according to the Schwarzschild
criterion while the enthalpy flux is outward directed. This occurs when the
heat conduction profile at the bottom of the CZ is smoothly varying, based
either on a Kramerslike opacity prescription as a function of temperature and
density or a static profile of a similar shape. We show that the subadiabatic
layer arises due to nonlocal energy transport by buoyantly driven downflows in
the upper parts of the CZ. Analysis of the force balance of the upflows and
downflows confirms that convection is driven by cooling at the surface. We find
that the commonly used prescription for the convective enthalpy flux being
proportional to the negative entropy gradient does not hold in the stably
stratified layers where the flux is positive. We demonstrate the existence of a
nongradient contribution to the enthalpy flux, which is estimated to be
important throughout the convective layer. A quantitative analysis of downflows
indicates a transition from a treelike structure where smaller downdrafts
merge into larger ones in the upper parts to a structure in the deeper parts
where a heightindependent number of strong downdrafts persist. This change of
flow topology occurs when a substantial subadiabatic layer is present in the
lower part of the CZ.

The presence of asymmetry between fermions of opposite handedness in plasmas
of relativistic particles can lead to exponential growth of a helical magnetic
field via a smallscale chiral dynamo instability known as the chiral magnetic
effect. Here, we show, using dimensional arguments and numerical simulations,
that this process produces through the Lorentz force chiral magnetically driven
turbulence. A k^{2} magnetic energy spectrum emerges via inverse transfer over
a certain range of wavenumbers k. The total chirality (magnetic helicity plus
normalized chiral chemical potential) is conserved in this system. Therefore,
as the helical magnetic field grows, most of the total chirality gets
transferred into magnetic helicity until the chiral magnetic effect terminates.
Quantitative results for height, slope, and extent of the spectrum are
obtained. Consequences of this effect for cosmic magnetic fields are discussed.

A turbulent dynamo in spherical geometry with an outer corona is simulated to
study the sign of magnetic helicity in the outer parts. In agreement with
earlier studies, the sign in the outer corona is found to be opposite to that
inside the dynamo. Lineofsight observations of polarized emission are
synthesized to explore the feasibility of using the local reduction of Faraday
depolarization to infer the sign of helicity of magnetic fields in the solar
corona. This approach was previously identified as an observational diagnostic
in the context of galactic magnetic fields. Based on our simulations, we show
that this method can be successful in the solar context if sufficient
statistics is gathered by using averages over ring segments in the corona
separately for the regions north and south of the solar equator.

The magnetic activity of the Sun becomes stronger and weaker over roughly an
11 year cycle, modulating the radiation and charged particle environment
experienced by the Earth as "space weather". Decades of observations from the
Mount Wilson Observatory have revealed that other stars also show regular
activity cycles in their Ca II H+K line emission, and identified two different
relationships between the length of the cycle and the rotation rate of the
star. Recent observations at higher cadence have allowed the discovery of
shorter cycles with periods between 13 yr. Some of these shorter cycles
coexist with longer cycle periods, suggesting that two underlying dynamos can
operate simultaneously. We combine these new observations with previous data,
and we show that the longer and shorter cycle periods agree remarkably well
with those expected from an earlier analysis based on the mean activity level
and the rotation period. The relative turbulent length scales associated with
the two branches of cyclic behavior suggest that a nearsurface dynamo may be
the dominant mechanism that drives cycles in more active stars, whereas a
dynamo operating in deeper layers may dominate in less active stars. However,
several examples of equally prominent long and short cycles have been found at
all levels of activity of stars younger than 2.3 Gyr. Deviations from the
expected cycle periods show no dependence on the depth of the convection zone
or on the metallicity. For some stars that exhibit longer cycles, we compute
the periods of shorter cycles that might be detected with future highcadence
observations.

Sunspots are of basic interest in the study of the Sun. Their relevance
ranges from them being an activity indicator of magnetic fields to being the
place where coronal mass ejections and flares erupt. They are therefore also an
important ingredient of space weather. Their formation, however, is still an
unresolved problem in solar physics. Observations utilize just 2D surface
information near the spot, but it is debatable how to infer deep structures and
properties from local helioseismology. For a long time, it was believed that
flux tubes rising from the bottom of the convection zone are the origin of the
bipolar sunspot structure seen on the solar surface. However, this theory has
been challenged, in particular recently by new surface observation,
helioseismic inversions, and numerical models of convective dynamos. In this
article we discuss another theoretical approach to the formation of sunspots:
the negative effective magnetic pressure instability. This is a largescale
instability, in which the total (kinetic plus magnetic) turbulent pressure can
be suppressed in the presence of a weak largescale magnetic field, leading to
a converging downflow, which eventually concentrates the magnetic field within
it. Numerical simulations of forced stratified turbulence have been able to
produce strong superequipartition flux concentrations, similar to sunspots at
the solar surface. In this framework, sunspots would only form close to the
surface due to the instability constraints on stratification and rotation.
Additionally, we present some ideas from local helioseismology, where we plan
to use the Hankel analysis to study the preemergence phase of a sunspot and to
constrain its deep structure and formation mechanism.

Recent direct numerical simulations (DNS) of largescale turbulent dynamos in
strongly stratified layers have resulted in surprisingly sharp bipolar
structures at the surface. Here we present new DNS of helically and
nonhelically forced turbulence with and without rotation and compare with
corresponding meanfield simulations (MFS) to show that these structures are a
generic outcome of a broader class of dynamos in densitystratified layers. The
MFS agree qualitatively with the DNS, but the period of oscillations tends to
be longer in the DNS. In both DNS and MFS, the sharp structures are produced by
converging flows at the surface and might be driven in nonlinear stage of
evolution by the Lorentz force associated with the largescale dynamodriven
magnetic field if the dynamo number is at least 2.5 times supercritical.

Analytic solutions of the meanfield induction equation predict a
nonoscillatory dynamo for homogeneous helical turbulence or constant alpha
effect in unbounded or periodic domains. Oscillatory dynamos are generally
thought impossible for constant alpha. We present an analytic solution for a
onedimensional bounded domain resulting in oscillatory solutions for constant
alpha, but different (Dirichlet and von Neumann or perfect conductor and
vacuum) boundary conditions on the two boundaries. We solve a second order
complex equation and superimpose two independent solutions to obey both
boundary conditions. The solution has timeindependent energy density. On one
end where the function value vanishes, the second derivative is finite, which
would not be correctly reproduced with sinelike expansion functions where a
node coincides with an inflection point. The field always migrates away from
the perfect conductor boundary toward the vacuum boundary, independently of the
sign of alpha. The obtained solution may serve as a benchmark for numerical
dynamo experiments and as a pedagogical illustration that oscillatory migratory
dynamos are possible with constant alpha.

We develop a twoscale formalism to determine global magnetic helicity
spectra in systems where the local magnetic helicity has opposite signs on both
sides of the equator, giving rise to cancelation with conventional methods. We
verify this approach using first a synthetic onedimensional magnetic field and
then twodimensional slices from a threedimensional alpha effecttype
dynamogenerated magnetic field with forced turbulence of opposite helicity
above and below the midplane of the domain. We then apply this formalism to
global solar synoptic vector magnetograms. To improve the statistics, data from
three consecutive Carrington rotations (21612163) are combined into a single
map. We find that the spectral magnetic helicity representative of the northern
hemisphere is negative at all wavenumbers and peaks at ~ 0.06 Mm^{1} (scales
around 100 Mm). There is no evidence of bihelical magnetic fields that are
found in threedimensional turbulence simulations of helicitydriven alpha
effecttype dynamos.

We perform numerical simulations of decaying hydrodynamic and
magnetohydrodynamic turbulence. We classify our timedependent solutions by
their evolutionary tracks in parametric plots between instantaneous scaling
exponents. We find distinct classes of solutions evolving along specific
trajectories toward points on a line of selfsimilar solutions. These
trajectories are determined by the underlying physics governing individual
cases, while the infrared slope of the initial conditions plays only a limited
role. In the helical case, even for a scaleinvariant initial spectrum
(inversely proportional to wavenumber k), the solution evolves along the same
trajectory as for a Batchelor spectrum (proportional to k^4.

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 meanfield simulations. Finally, the relevance
to sunspot formation is discussed.

The effects of a nongradient 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 latetype 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 Schwarzschildstably
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 largescale flows with speeds
and scales expected from MLT.

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 largescale 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 largescale magnetic field
evolution. Numerical simulations have been carried out for nonhelical and
helical magnetic field configurations. The numerical results show the
possibility of inverse transfer of energy in magnetically dominated nonhelical
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 largescale 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.

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\omega$ 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.

We use high resolution direct numerical simulations (DNS) to show that
helical turbulence can generate significant largescale fields even in the
presence of strong smallscale dynamo action. During the kinematic stage, the
unified large/smallscale dynamo grows fields with a shapeinvariant
eigenfunction, with most power peaked at small scales or large $k$, as in
\citet{SB14}. Nevertheless, the largescale 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 smallscale field, increases significantly
from the kinematic stage to saturation. This implies that the smallscale 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/smallscale dynamo in
helical turbulence.