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 imposed 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 non-helical 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 helicity. 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.

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 much more slowly ($\propto t^{-2/5}$ instead of $\propto t^{-2/3}$ in MHD). Even for a nonhelical magnetic field there is inverse cascading with an increase of magnetic energy at the largest possible scale in the system. The strength of this effect depends on the strength of the magnetic field and disappears for weak fields. The inertial range scaling is compatible with earlier findings for the forced Hall cascade, i.e., proportional to $k^{-7/3}$, but in the decaying case, the subinertial range spectrum steepens to a novel $k^5$ slope instead of $k^4$, as in MHD. The energy dissipation is found to be proportional to the magnetic field to the 5th power for helical fields. For neutron star conditions with an rms magnetic field of $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 would manifest itself through X-ray emission.

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.

The possibility of explaining shear flow dynamos in terms of a magnetic shear-current effect (SC) is examined. A competing explanation is the incoherent $\alpha$-shear dynamo effect. Our primary diagnostics is the turbulent magnetic diffusivity tensor, in particular the off-diagonal diffusivity tensor component $\eta_{yx}$, a systematically negative sign of which would imply coherent dynamo action through SC. To be able to measure turbulent transport coefficients from systems with strong magnetic fluctuations, we present an extension of the test-field method that is capable of such measurements in the case where the pressure gradient term is dropped from the MHD equations, which is a nonlinear TFM (NLTFM). The hydrodynamic equation is related to Burger's equation and the resulting flows are referred to as magnetized burgulence. We use both kinetic and magnetic forcings, to mimic cases without and with simultaneous small-scale dynamo action. When the simplified MHD flow is forced kinetically, negative $\eta_{yx}$ values are obtained with exponential growth in both the radial and azimuthal magnetic field components. Using stochastic monochromatic magnetic forcing, the exponential growth is no longer seen, and NLTFM yields positive values of $\eta_{yx}$. In an attempt to recover the exponential growth seen in the kinetically forced case, we also employ an alternative forcing function with lowest $k$ vectors being removed - in this case the exponential growth is recovered, but the NLTFM results do not change dramatically. Our analysis of the dynamo numbers for the coherent SC and incoherent $\alpha$ and SC effects show that the incoherent effects are the main drivers of the dynamo instability in majority of cases. The coherent SC effect is mildly enhancing the dynamo action in the kinetically forced cases, while we find no evidence for magnetic SC in our simulations.

Stellar winds are an integral part of the underlying dynamo, the motor of stellar activity. The wind controls the star's angular momentum loss, which depends on the magnetic field geometry which, in turn, varies significantly in time and latitude. Here we study basic properties of a self-consistent model that includes simple representations of both the global stellar dynamo in a spherical shell and the exterior in which the wind accelerates and becomes supersonic. We numerically solve an axisymmetric mean-field model for the induction, momentum, and continuity equations using an isothermal equation of state. The model allows for the simultaneous generation of a mean magnetic field and the development of a Parker wind. The resulting flow is transonic at the critical point, which we arrange to be between the inner and outer radii of the model. The boundary conditions are assumed to be such that the magnetic field is antisymmetric about the equator, i.e., dipolar. At the solar rotation rate, the dynamo is oscillatory and of $\alpha^2$ type. In most of the domain, the magnetic field corresponds to that of a split monopole. The magnetic energy flux is largest between the stellar surface and the critical point. At rapid rotation of up to 50 times the solar value, most of the magnetic field is lost along the axis within the inner tangential cylinder of the model. The model reveals unexpected features that are not generally anticipated from models that are designed to reproduce the solar wind: highly variable angular momentum flux fluxes even from just an $\alpha^2$ dynamo in the star. For rapid rotation, magnetic fields are ejected mostly along the axis, where the wind speed is reduced.

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.

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.

We discuss selected aspects regarding the magnetic field evolution of solar-type stars. Most of the stars with activity cycles are in the range where the normalized chromospheric Calcium emission increases linearly with the inverse Rossby number. For Rossby numbers below about a quarter of the solar value, the activity saturates and no cycles have been found. For Rossby numbers above the solar value, again no activity cycles have been found, but now the activity goes up again for a major fraction of the stars. Rapidly rotating stars show nonaxisymmetric large-scale magnetic fields, but there is disagreement between models and observations regarding the actual value of the Rossby number where this happens. We also discuss the prospects of detecting the sign of magnetic helicity using various linear polarization techniques both at the stellar surface using the parity-odd contribution to linear polarization and above the surface using Faraday rotation.

We perform direct numerical simulations of magnetohydrodynamic turbulence in the early universe and numerically compute the resulting stochastic background of gravitational waves and relic magnetic fields. These simulations do not make the simplifying assumptions of earlier analytic calculations. If the turbulence is assumed to have an energy-carrying scale that is about a hundredth of the Hubble scale at the time of generation, as expected in a first-order phase transition, the peak of gravitational wave power will be in the mHz frequency range for a signal produced at the electroweak scale. The efficiency of GW production varies significantly with how the turbulence is driven. Detectability of turbulence at the electroweak scale by the planned Laser Interferometer Space Antenna requires anywhere from 0.1 to 10 percent of the thermal plasma energy density to be in plasma motions or magnetic fields, depending on the model of the driving process. At frequencies below the peak value, our results predict a new universal form of the spectrum with more power than previously thought. The enhanced low-frequency tail makes turbulence at significantly higher energy scales detectable.

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.

(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.

The alpha effect is believed to play a key role in the generation of the solar magnetic field. A fundamental test for its significance in the solar dynamo is to look for magnetic helicity of opposite signs in the two hemispheres, and at small and large scales. However, measuring magnetic helicity is compromised by the inability to fully infer the magnetic field vector from observations of solar spectra, caused by what is known as the "pi ambiguity" of spectropolarimetric observations. We decompose linear polarisation into parity-even and parity-odd E and B polarisations, which are not affected by the "pi ambiguity". Furthermore, we study whether the correlations of spatial Fourier spectra of B and parity-even quantities such as E or temperature T are a robust proxy for magnetic helicity of solar magnetic fields. We analyse polarisation measurements of active regions observed by the Helioseismic and Magnetic Imager on board the Solar Dynamics observatory. Theory predicts the magnetic helicity of active regions to have, statistically, opposite signs in the two hemispheres. We then compute the parity-odd E B and T B correlations, and test for systematic preference of their sign based on the hemisphere of the active regions. We find that: (i) E B and T B correlations are a reliable proxy for magnetic helicity, when computed from linear polarisation measurements away from spectral line cores, and (ii) E polarisation reverses its sign close to the line core. Our analysis reveals Faraday rotation to not have a significant influence on the computed parity-odd correlations. The EB decomposition of linear polarisation appears to be a good proxy for magnetic helicity independent of the "pi ambiguity". This allows us to routinely infer magnetic helicity directly from polarisation measurements.

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.

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.

Explicit radiation hydrodynamic simulations of the atmospheres of massive stars and of convection in accretion discs around white dwarfs suffer from prohibitively short time steps due to radiation. This constraint is related to the cooling time rather than the radiative pressure, which also becomes important in hot stars and discs. We show that the radiative time step constraint is governed by the minimum of the sum of the optically thick and thin contributions rather than the smaller one of the two. In simulations with the Pencil Code, their weighting fractions are found empirically. In three-dimensional convective accretion disc simulations, the Deardorff term is found to be the main contributor to the enthalpy flux rather than the superadiabatic gradient. We conclude with a discussion of how the radiative time step problem could be mitigated in certain types of investigations.

If the alpha effect plays a role in the generation of the Sun's magnetic field, the field should show evidence of magnetic helicity of opposite signs at large and small length scales. Measuring this faces two challenges: (i) in weak-field regions, horizontal field measurements are unreliable because of the pi ambiguity, and (ii) one needs a truly global approach to computing helicity spectra in the case where one expects a sign reversal across the equator at all wavenumbers. Here we develop such a method using spin-2 spherical harmonics to decompose the linear polarization in terms of the parity-even and parity-odd E and B polarizations, respectively. Using simple one- and two-dimensional models, we show that the product of the spectral decompositions of E and B, taken at spherical harmonic degrees that are shifted by one, can act as a proxy of the global magnetic helicity with a sign that represents that in the northern hemisphere. We then apply this method to the analysis of solar synoptic vector magnetograms, from which we extract a pseudo-polarization corresponding to a "pi-ambiguated" magnetic field, i.e., a magnetic field vector that has no arrow. We find a negative sign of the global EB helicity proxy at spherical harmonic degrees of around 6. This could indicate a positive magnetic helicity at large length scales, but the spectrum fails to capture clear evidence of the well-known negative magnetic helicity at smaller scales. This method might also be applicable to stellar and Galactic polarization data.

Hydromagnetic turbulence produced during phase transitions in the early universe can be a powerful source of stochastic gravitational waves (GWs). GWs can be modelled by the linearised spatial part of the Einstein equations sourced by the Reynolds and Maxwell stresses. We have implemented two different GW solvers into the {\sc Pencil Code}---a code which uses a third order timestep and sixth order finite differences. Using direct numerical integration of the GW equations, we study the appearance of a numerical degradation of the GW amplitude at the highest wavenumbers, which depends on the length of the timestep---even when the Courant--Friedrichs--Lewy condition is ten times below the stability limit. This degradation leads to a numerical error, which is found to scale with the third power of the timestep. A similar degradation is not seen in the magnetic and velocity fields. To mitigate numerical degradation effects, we alternatively use the exact solution of the GW equations under the assumption that the source is constant between subsequent timesteps. This allows us to use a much longer timestep, which cuts the computational cost by a factor of about ten.

We compute magnetic helicity and energy spectra from about 2485 patches of about 100 megameters (Mm) side length on the solar surface using data from Hinode during 2006--2017. An extensive database is assembled where we list magnetic energy and helicity, large- and small-scale magnetic helicity, mean current helicity density, fractional magnetic helicity, and correlation length along with the Hinode map identification number (MapID), as well as Carrington latitude and longitude for each MapID. While there are departures from the hemispheric sign rule for magnetic and current helicities, the weak trend reported here is in agreement with the previous results. This is argued to be a physical effect associated with the dominance of individual active regions that contribute more strongly in the better resolved Hinode maps. In comparison with earlier work, the typical correlation length is found to be 6-8 Mm, while the length scale relating magnetic and current helicity to each other is found to be around 1.4 Mm.

We study the effects of ambipolar diffusion (AD) on hydromagnetic turbulence. In most of the cases, we use the single fluid approximation where the drift velocity between charged and neutral particles is proportional to the Lorentz force. In two cases we also compare with the corresponding two-fluid model, where ionization and recombination are included in the continuity and momentum equations for the neutral and charged species. We consider the regime of large magnetic Prandtl number, relevant to the interstellar medium. The magnetic field properties are found to be well represented by the single fluid approximation. We quantify the effects of AD on total and spectral kinetic and magnetic energies, the Ohmic and AD dissipation rates, the statistics of the magnetic field, the current density, and the linear polarization as measured by the rotationally invariant E and B mode polarizations. We show that the kurtosis of the magnetic field decreases with increasing AD. The E mode polarization changes its skewness from positive values for small AD to negative ones for large AD. Even when AD is weak, changes in AD have a marked effect on the skewness and kurtosis of E, and only a weak effect on those of B. These results open the possibility of employing E and B mode polarizations as diagnostic tools for characterizing turbulent properties of the interstellar medium.

We show that at large magnetic Prandtl numbers, the Lorentz force does work on the flow at small scales and drives fluid motions, whose energy is dissipated viscously. This situation is opposite to that in a normal dynamo, where the flow does work against the Lorentz force. We compute the spectral conversion rates between kinetic and magnetic energies for several magnetic Prandtl numbers and show that normal (forward) dynamo action occurs on large scales over a progressively narrower range of wavenumbers as the magnetic Prandtl number is increased. At higher wavenumbers, reversed dynamo action occurs, i.e., magnetic energy is converted back into kinetic energy at small scales. We demonstrate this in both direct numerical simulations forced by volume stirring and in large eddy simulations of solar convectively driven small-scale dynamos. Low density plasmas such as stellar coronae tend to have large magnetic Prandtl numbers, i.e., the viscosity is large compared with the magnetic diffusivity. The regime in which viscous dissipation dominates over resistive dissipation for large magnetic Prandtl numbers was also previously found in large eddy simulations of the solar corona, i.e., our findings are a more fundamental property of MHD that is not just restricted to dynamos. Viscous energy dissipation is a consequence of positive Lorentz force work, which may partly correspond to particle acceleration in close-to-collisionless plasmas. This is, however, not modeled in the MHD approximation employed. By contrast, resistive energy dissipation on current sheets is expected to be unimportant in stellar coronae.

Large-scale magnetic fields in stars and galaxies are thought to arise by mean-field dynamo action due to the combined influence of both helical turbulence and shear. Those systems are also highly conducting and the turbulence therein leads to a fluctuation (or small-scale) dynamo which more rapidly amplifies magnetic field fluctuations on the eddy scales and smaller. Will this then interfere with and suppress the mean (or large-scale) field growth? Using direct numerical simulations of helical turbulence (with and without shear), we identify a novel quasi-kinematic large-scale dynamo which operates as the small-scale dynamo saturates. Thus both dynamos operate efficiently, one after the other, and lead to the generation of significant large-scale fields.

Recent numerical work in helioseismology has shown that a periodically varying subsurface magnetic field leads to a fanning of the $f$-mode, which emerges from the density jump at the surface. In an attempt to model a more realistic situation, we now modulate this periodic variation with an envelope, giving thus more emphasis on localized bipolar magnetic structures in the middle of the domain. Some notable findings are: (i) compared to the purely hydrodynamic case, the strength of the $f$-mode is significantly larger at high horizontal wavenumbers $k$, but the fanning is weaker for the localized subsurface magnetic field concentrations investigated here than the periodic ones studied earlier; (ii) when the strength of the magnetic field is enhanced at a fixed depth below the surface, the fanning of the $f$-mode in the $k\omega$ diagram increases proportionally in such a way that the normalized $f$-mode strengths remain nearly the same in different such cases; (iii) the unstable Bloch modes reported previously in case of harmonically varying magnetic fields are now completely absent when more realistic localized magnetic field concentrations are imposed beneath the surface, thus suggesting that the Bloch modes are unlikely to be supported during most phases of the solar cycle; (iv) the $f$-mode strength appears to depend also on the depth of magnetic field concentrations such that it shows a relative decrement when the maximum of the magnetic field is moved to a deeper layer. We argue that detections of $f$-mode perturbations such as those being explored here could be effective tracers of solar magnetic fields below the photosphere before these are directly detectable as visible manifestations in terms of active regions or sunspots.

The chiral magnetic effect (CME) is a quantum relativistic effect that describes the appearance of an additional electric current along a magnetic field. It is caused by an asymmetry between the number densities of left- and right-handed fermions, which can be maintained at high energies when the chirality flipping rate can be neglected, for example in the early Universe. The inclusion of the CME in the Maxwell equations leads to a modified set of MHD equations. We discuss how the CME is implemented in the PENCIL CODE. The CME plays a key role in the evolution of magnetic fields since it results in a dynamo effect associated with an additional term in the induction equation. This term is formally similar to the $\alpha$ effect in classical mean-field MHD. However, the chiral dynamo can operate without turbulence and is associated with small spatial scales that can be, in the case of the early Universe, orders of magnitude below the Hubble radius. A chiral $\alpha_\mu$ effect has also been identified in mean-field theory. It occurs in the presence of turbulence but is not related to kinetic helicity. Depending on the plasma parameters, chiral dynamo instabilities can amplify magnetic fields over many orders of magnitude. These instabilities can affect the propagation of MHD waves, which is demonstrated by our DNS. We also study the coupling between the evolution of the chiral chemical potential and the ordinary chemical potential, which is proportional to the sum of the number densities of left- and right-handed fermions. An important consequence of this coupling is the emergence of chiral magnetic waves (CMWs). We confirm numerically that linear CMWs and MHD waves are not interacting. Our simulations suggest that the chemical potential has only a minor effect on the non-linear evolution of the chiral dynamo.

We study the evolution of kinetic and magnetic energy spectra in magnetohydrodynamic flows in the presence of strong cross helicity. For forced turbulence, we find weak inverse transfer of kinetic energy toward the smallest wavenumber. This is plausibly explained by the finiteness of scale separation between the injection wavenumber and the smallest wavenumber of the domain, which here is a factor of 15. In the decaying case, there is a slight increase at the smallest wavenumber, which is probably explained by the dominance of kinetic energy over magnetic energy at the smallest wavenumbers. Within a range of wavenumbers covering almost an order of magnitude the decay is purely exponential, which is argued to be a consequence of a suppression of nonlinearity due to the presence of strong cross helicity.

Much work on turbulent three-dimensional dynamos has been done using triply periodic domains, in which there are no magnetic helicity fluxes. Here we present simulations where the turbulent intensity is still nearly homogeneous, but now there is a perfect conductor boundary condition on one end and a vertical field or pseudo-vacuum condition on the other. This leads to migratory dynamo waves. Good agreement with a corresponding analytically solvable alpha^2 dynamo is found. Magnetic helicity fluxes are studied in both types of models. It is found that at moderate magnetic Reynolds numbers, most of the magnetic helicity losses occur at large scales. Whether this changes at even larger magnetic Reynolds numbers, as required for alleviating the catastrophic dynamo quenching problem, remains still unclear.

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 cross-over 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 small-scale dynamo action. Our results have applications to a wide range of experimental dynamos and astrophysical time-dependent plasmas, including primordial turbulence in the early universe.

We test the sensitivity of hydrodynamic and magnetohydrodynamic turbulent convection simulations with respect to Mach number, thermal and magnetic boundary conditions, and the centrifugal force. We find that varying the luminosity, which also controls the Mach number, has only a minor effect on the large-scale dynamics. A similar conclusion can also be drawn from the comparison of two formulations of the lower magnetic boundary condition with either vanishing electric field or current density. The centrifugal force has an effect on the solutions, but only if its magnitude with respect to acceleration due to gravity is by two orders of magnitude greater than in the Sun. Finally, we find that the parameterisation of the photospheric physics, either by an explicit cooling term or enhanced radiative diffusion, is more important than the thermal boundary condition. In particular, runs with cooling tend to lead to more anisotropic convection and stronger deviations from the Taylor-Proudman state. In summary, the fully compressible approach taken here with the Pencil Code is found to be valid, while still allowing the disparate timescales to be taken into account.

We consider the effect of a subadiabatic layer at the base of the convection zone on convection itself and the associated large-scale dynamos in spherical wedge geometry. We use a heat conduction prescription based on the Kramers opacity law which allows the depth of the convection zone to dynamically adapt to changes in the physical characteristics such as rotation rate and magnetic fields. We find that the convective heat transport is strongly concentrated toward the equatorial and polar regions in the cases without a substantial radiative layer below the convection zone. The presence of a stable layer below the convection zone significantly reduces the anisotropy of radial enthalpy transport. Furthermore, the dynamo solutions are sensitive to subtle changes in the convection zone structure. We find that the kinetic helicity changes sign in the deeper parts of the convection zone at high latitudes in all runs. This region expands progressively toward the equator in runs with a thicker stably stratified layer.

Solar and stellar dynamos shed small-scale and large-scale magnetic helicity of opposite signs. However, solar wind observations and simulations have shown that some distance above the dynamo both the small-scale and large-scale magnetic helicities have reversed signs. With realistic simulations of the solar corona above an active region now being available, we have access to the magnetic field and current density along coronal loops. We show that a sign reversal in the horizontal averages of the magnetic helicity occurs when the local maximum of the plasma beta drops below unity and the field becomes nearly fully force free. Hence, this reversal is expected to occur well within the solar corona and would not directly be accessible to in-situ measurements with the Parker Solar Probe or SolarOrbiter. We also show that the reversal is associated with subtle changes in the relative dominance of structures with positive and negative magnetic helicity.

The emergence of dipolar magnetic features on the solar surface is an idealization. Most of the magnetic flux emergence occurs in complex multipolar regions. Here, we show that the surface pattern of magnetic structures alone can reveal the sign of the underlying magnetic helicity in the nearly force-free coronal regions above. The sign of the magnetic helicity can be predicted to good accuracy by considering the three-dimensional position vectors of three spots on the sphere ordered by their relative strengths at the surface and compute from them the skew product. This product, which is a pseudoscalar, is shown to be a good proxy for the sign of the coronal magnetic helicity.

The full-sky Planck polarization data at 850um revealed unexpected properties of the E and B mode power spectra of dust emission in the interstellar medium (ISM). The positive cross-correlation between the total dust intensity, T, with the B modes has raised new questions about the physical mechanisms that affect dust polarization, such as the Galactic magnetic-field structure. This is key both to better understanding ISM dynamics and to accurately describing Galactic foregrounds to the polarization of the Cosmic Microwave Background (CMB). In this theoretical paper we investigate the possibility that the observed cross-correlations in the dust polarization power spectra, and specifically between T and B, can be related to a parity-odd quantity in the ISM such as the magnetic helicity. We produce synthetic dust polarization data, derived from 3D analytical toy models of density structures and helical magnetic fields, to compare with the E and B modes of observations. Focusing on the observed T-B correlation, we propose a new line of interpretation of the Planck observations based on a large-scale helical component of the Galactic magnetic field in the solar neighborhood. Our analysis shows that: I) the sign of magnetic helicity does not affect E and B modes for isotropic magnetic-field configurations; II) helical magnetic fields threading interstellar filaments cannot reproduce the Planck results; III) a weak helical left-handed magnetic field structure in the solar neighborhood may explain the T-B correlation seen in the Planck data. This work suggests a new perspective for the interpretation of the dust polarization power spectra, which strongly supports the imprint of a large-scale structure of the Galactic magnetic field in the solar neighborhood.

Gradient- and curl-type or E- and B-type polarizations have been routinely analyzed to study the physics contributing to the cosmic microwave background polarization and galactic foregrounds. They characterize the parity-even and parity-odd properties of the underlying physical mechanisms, for example hydromagnetic turbulence in the case of dust polarization. Here we study spectral correlation functions characterizing the parity-even and parity-odd parts of linear polarization for homogeneous and inhomogeneous turbulence to show that only the inhomogeneous helical case can give rise to a parity-odd polarization signal. We also study nonhelical turbulence and suggest that a strong nonvanishing (here negative) skewness of the E polarization is responsible for an enhanced ratio of the EE to the BB (quadratic) correlation in both helical and nonhelical cases. This could explain the enhanced EE/BB ratio observed recently for dust polarization. We close with a preliminary assessment of using linear polarization of the Sun to characterize its helical turbulence without being subjected to the pi ambiguity that magnetic inversion techniques have to address.

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.

We develop a mean-field theory of compressibility effects in turbulent magnetohydrodynamics and passive scalar transport using the quasi-linear approximation and the spectral $\tau$-approach. We find that compressibility decreases the $\alpha$ effect and the turbulent magnetic diffusivity both at small and large magnetic Reynolds numbers, Rm. Similarly, compressibility decreases the turbulent diffusivity for passive scalars both at small and large P\'eclet numbers, Pe. On the other hand, compressibility does not affect the effective pumping velocity of the magnetic field for large Rm, but it decreases it for small Rm. Density stratification causes turbulent pumping of passive scalars, but it is found to become weaker with increasing compressibility. No such pumping effect exists for magnetic fields. However, compressibility results in a new passive scalar pumping effect from regions of low to high turbulent intensity both for small and large P\'eclet numbers. It can be interpreted as compressible turbophoresis of noninertial particles and gaseous admixtures, while the classical turbophoresis effect exists only for inertial particles and causes them to be pumped to regions with lower turbulent intensity.

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.

The formation of sunspots and starspots is not yet fully understood and is therefore one of the major open problems in solar and stellar physics. Magnetic flux concentrations can be produced by the negative effective magnetic pressure instability (NEMPI). This instability is strongly suppressed by rotation. However, the presence of an outer coronal envelope was previously found to strengthen the flux concentrations and make them more prominent. It also allows for the formation of bipolar regions (BRs). We want to know whether the presence of an outer coronal envelope also changes the excitation conditions and the rotational dependence of NEMPI. We use direct numerical simulations and mean-field simulations. We adopt a simple two-layer model of turbulence that mimics the jump between the convective turbulent and coronal layers below and above the surface of a star, respectively. The computational domain is Cartesian and located at a certain latitude of a rotating sphere. We investigate the effects of rotation on NEMPI by changing the Coriolis number, the latitude, and the box resolution. Rotation has a strong impact on the process of BR formation. Even rather slow rotation is found to suppress their formation. However, increasing the imposed magnetic field strength also makes the structures stronger and alleviates the rotational suppression somewhat. The presence of a coronal layer itself does not significantly alleviate the effects of rotational suppression.

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 scale-invariant 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 scale-invariant spectrum changes with time. It becomes a steeper weak-turbulence spectrum at progressively larger scales. We show numerically that the critical length scale where this happens is the turbulent-diffusive scale, which increases with the square root of time.

We develop a formalism that treats the calculation of solar kinetic energy and cross helicity spectra in an equal manner to that of magnetic energy and helicity spectra. The magnetic helicity spectrum is shown to be equal to the vertical part of the current helicity spectrum divided by the square of the wavenumber. For the cross helicity, we apply the recently developed two-scale approach globally over an entire active region to account for the sign change between the two polarities. Using vector magnetograms and Dopplergrams of NOAA 11158 and 12266, we show that kinetic and magnetic energy spectra have similar slopes 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. Using the two-scale method for NOAA 12266, the global cross helicity spectrum is found to be particularly steep, similarly to what has previously been found in theoretical models of spot generation. In the quiet Sun, by comparison, the cross helicity spectrum is found to be small.

Recent advances in mean-field theory are reviewed and applications to the Sun, late-type stars, accretion disks, galaxies, and the early Universe are discussed. We focus particularly on aspects of spatio-temporal nonlocality, which is one of the main insights that emerged from applying the test-field 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 antisolar-like 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.

Small-scale dynamos are expected to operate in all astrophysical fluids that are turbulent and electrically conducting, for example the interstellar medium, stellar interiors, and accretion disks, where they may also be affected by or competing with large-scale dynamos. However, the possibility of small-scale dynamos being excited at small and intermediate ratios of viscosity to magnetic diffusivity (the magnetic Prandtl number) has been debated, and the possibility of them depending on the large-scale forcing wavenumber has been raised. Here we show, using four values of the forcing wavenumber, that the small-scale dynamo does not depend on the scale-separation between the size of the simulation domain and the integral scale of the turbulence, i.e., the forcing scale. Moreover, the spectral bottleneck in turbulence, which has been implied as being responsible for raising the excitation conditions of small-scale dynamos, is found to be invariant under changing the forcing wavenumber. However, when forcing at the lowest few wavenumbers, the effective forcing wavenumber that enters in the definition of the magnetic Reynolds number is found to be about twice the minimum wavenumber of the domain. Our work is relevant to future studies of small-scale dynamos, of which several applications are being discussed.

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 right-handed 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 small-scale chiral dynamo instability; (2) production of chiral magnetically driven turbulence and excitation of a large-scale 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 large-scale magnetic field and its characteristic scale measured in the numerical simulations agree well with theoretical predictions based on mean-field theory. The previously discussed two-stage 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 proto-neutron stars are discussed.

We seek to understand the transition from nearly axisymmetric configurations at solar rotation rates to nonaxisymmetric configurations for rapid rotation using 3D numerical simulations of turbulent convection and considering rotation rates between 1 and 30 times the solar value. We find a transition from axi- to nonaxisymmetric solutions at around 1.8 times the solar rotation rate. This transition coincides with a change in the rotation profile from antisolar- to solar-like differential rotation with a faster equator and slow poles. In the solar-like rotation regime, the field configuration consists of an axisymmetric oscillatory field accompanied by an m=1 azimuthal mode (two active longitudes), which also shows temporal variability. At slow (rapid) rotation, the axisymmetric (nonaxisymmetric) mode dominates. The axisymmetric mode produces latitudinal dynamo waves with polarity reversals, while the nonaxisymmetric mode often exhibits a drift in the rotating reference frame and the strength of the active longitudes changes cyclically over time between the different hemispheres. Most of the obtained dynamo solutions exhibit cyclic variability either caused by latitudinal or azimuthal dynamo waves. In an activity-period diagram, the cycle lengths normalized by the rotation period form two different populations as a function of rotation rate or magnetic activity level. The slowly rotating axisymmetric population lies close to what is called the inactive branch in observations, while the rapidly rotating models are close to the superactive branch with a declining cycle to rotation frequency ratio with increasing rotation rate. We can successfully reproduce the transition from axi- to nonaxisymmetric dynamo solutions for high rotation rates, but high-resolution simulations are required to limit the effect of rotational quenching of convection at rotation rates above 20 times the solar value.

Using a recently developed two-scale 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 Long-term Investigations of the Sun} (SOLIS) telescope during January 2010-July 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 large-scale dynamo with helical turbulence operating in the Sun. More than $20\%$ of the analyzed maps, however, show violations of the expected sign rule.

Small-scale dynamo action is often held responsible for the generation of quiet-Sun magnetic fields. We aim to determine the excitation conditions and saturation level of small-scale dynamos in non-rotating turbulent convection at low magnetic Prandtl numbers. We use high resolution direct numerical simulations of weakly stratified turbulent convection. We find that the critical magnetic Reynolds number for dynamo excitation increases as the magnetic Prandtl number is decreased, which might suggest that small-scale dynamo action is not automatically evident in bodies with small magnetic Prandtl numbers as the Sun. As a function of the magnetic Reynolds number (${\rm Rm}$), the growth rate of the dynamo is consistent with an ${\rm Rm}^{1/2}$ scaling. No evidence for a logarithmic increase of the growth rate with ${\rm Rm}$ is found.

High precision photometry of solar-like 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 solar-like 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 B-V color. The resulting rotation--activity 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.

To explain the large-scale magnetic field of the Sun and other bodies, mean-field dynamo theory is commonly applied where one solves the averaged equations for the mean magnetic field. However, the standard approach breaks down when the scale of the turbulent eddies becomes comparable to the scale of the variations of the mean magnetic field. Models showing sharp magnetic field structures have therefore been regarded as unreliable. Our aim is to look for new effects that occur when we relax the restrictions of the standard approach, which becomes particularly important at the bottom of the convection zone where the size of the turbulent eddies is comparable to the depth of the convection zone itself. We approximate the underlying integro-differential equation by a partial differential equation corresponding to a reaction-diffusion type equation for the mean electromotive force, making an approach that is nonlocal in space and time feasible under conditions where spherical geometry and nonlinearity are included. In agreement with earlier findings, spatio-temporal nonlocality lowers the excitation conditions of the dynamo. Sharp structures are now found to be absent. However, in the surface layers the field remains similar to before.

We investigate dynamo action in global compressible solar-like convective dynamos in the framework of mean-field theory. We simulate a solar-type star in a wedge-shaped spherical shell, where the interplay between convection and rotation self-consistently drives a large-scale dynamo. To analyze the dynamo mechanism we apply the test-field method for azimuthally ($\phi$) averaged fields to determine the 27 turbulent transport coefficients of the electromotive force, of which six are related to the $\alpha$ tensor. This method has previously been used either in simulations in Cartesian coordinates or in the geodynamo context and is applied here for the first time to fully compressible simulations of solar-like dynamos. We find that the $\phi\phi$-component of the $\alpha$ tensor does not follow the profile expected from that of kinetic helicity. The turbulent pumping velocities significantly alter the effective mean flows acting on the magnetic field and therefore challenge the flux transport dynamo concept. All coefficients are significantly affected by dynamically important magnetic fields. Quenching as well as enhancement are being observed. This leads to a modulation of the coefficients with the activity cycle. The temporal variations are found to be comparable to the time-averaged values and seem to be responsible for a nonlinear feedback on the magnetic field generation. Furthermore, we quantify the validity of the Parker-Yoshimura rule for the equatorward propagation of the mean magnetic field in the present case.

Context: Convectively-driven flows play a crucial role in the dynamo processes that are responsible for producing magnetic activity in stars and planets. It is still not fully understood why many astrophysical magnetic fields have a significant large-scale component. Aims: Our aim is to investigate the dynamo properties of compressible convection in a rapidly rotating Cartesian domain, focusing upon a parameter regime in which the underlying hydrodynamic flow is known to be unstable to a large-scale vortex instability. Methods: The governing equations of three-dimensional nonlinear magnetohydrodynamics (MHD) are solved numerically. Different numerical schemes are compared and we propose a possible benchmark case for other similar codes. Results: In keeping with previous related studies, we find that convection in this parameter regime can drive a large-scale dynamo. The components of the mean horizontal magnetic field oscillate, leading to a continuous overall rotation of the mean field. Whilst the large-scale vortex instability dominates the early evolution of the system, it is suppressed by the magnetic field and makes a negligible contribution to the mean electromotive force that is responsible for driving the large-scale dynamo. The cycle period of the dynamo is comparable to the ohmic decay time, with longer cycles for dynamos in convective systems that are closer to onset. In these particular simulations, large-scale dynamo action is found only when vertical magnetic field boundary conditions are adopted at the upper and lower boundaries. Strongly modulated large-scale dynamos are found at higher Rayleigh numbers, with periods of reduced activity ("grand minima"-like events) occurring during transient phases in which the large-scale vortex temporarily re-establishes itself, before being suppressed again by the magnetic field.

We consider a scale-invariant 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 three-dimensional simulations, we show that an initially scale-invariant 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 magnetohydrodynamic (MHD) description of plasmas with relativistic particles necessarily includes an additional new field, the chiral chemical potential associated with the axial charge (i.e., the number difference between right- and left-handed relativistic fermions). This chiral chemical potential gives rise to a contribution to the electric current density of the plasma (\emph{chiral magnetic effect}). We present a self-consistent treatment of the \emph{chiral MHD equations}, which include the back-reaction of the magnetic field on a chiral chemical potential and its interaction with the plasma velocity field. A number of novel phenomena are exhibited. First, we show that the chiral magnetic effect decreases the frequency of the Alfv\'{e}n wave for incompressible flows, increases the frequencies of the Alfv\'{e}n wave and of the fast magnetosonic wave for compressible flows, and decreases the frequency of the slow magnetosonic wave. Second, we show that, in addition to the well-known laminar chiral dynamo effect, which is not related to fluid motions, there is a dynamo caused by the joint action of velocity shear and chiral magnetic effect. In the presence of turbulence with vanishing mean kinetic helicity, the derived mean-field chiral MHD equations describe turbulent large-scale dynamos caused by the chiral alpha effect, which is dominant for large fluid and magnetic Reynolds numbers. The chiral alpha effect is due to an interaction of the chiral magnetic effect and fluctuations of the small-scale current produced by tangling magnetic fluctuations (which are generated by tangling of the large-scale magnetic field by sheared velocity fluctuations). These dynamo effects may have interesting consequences in the dynamics of the early universe, neutron stars, and the quark--gluon plasma.

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 mean-field 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 Kramers-type 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 two-dimensional models. We use artificially enhanced mean-field 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.

Using numerical simulations of forced turbulence, it is shown that for magnetic Reynolds numbers larger than unity, i.e., beyond the regime of quasilinear theory, the turbulent magnetic diffusivity attains an additional negative contribution that is quadratic in the kinetic helicity. In particular, for large magnetic Reynolds numbers, the turbulent magnetic diffusivity without helicity is about twice the value with helicity. Such a contribution was not previously anticipated, but, as we discuss, it turns out to be important when accurate estimates of the turbulent magnetic diffusivity are needed.

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 Kramers-like 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 non-gradient 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 tree-like structure where smaller downdrafts merge into larger ones in the upper parts to a structure in the deeper parts where a height-independent 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 small-scale 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. Line-of-sight 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 1-3 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 near-surface 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 high-cadence observations.

Small-scale dynamos are ubiquitous in a broad range of turbulent flows with large-scale shear, ranging from solar and galactic magnetism to accretion disks, cosmology and structure formation. Using high-resolution direct numerical simulations we show that in non-helically forced turbulence with zero mean magnetic field, large-scale shear supports small-scale dynamo action, i.e., the dynamo growth rate increases with shear and shear enhances or even produces turbulence, which, in turn, further increases the dynamo growth rate. When the production rates of turbulent kinetic energy due to shear and forcing are comparable, we find scalings for the growth rate $\gamma$ of the small-scale dynamo and the turbulent velocity $u_{\rm rms}$ with shear rate $S$ that are independent of the magnetic Prandtl number: $\gamma \propto |S|$ and $u_{\rm rms} \propto |S|^{2/3}$. For large fluid and magnetic Reynolds numbers, $\gamma$, normalized by its shear-free value, depends only on shear. Having compensated for shear-induced effects on turbulent velocity, we find that the normalized growth rate of the small-scale dynamo exhibits the scaling, $\widetilde{\gamma}\propto |S|^{2/3}$, arising solely from the induction equation for a given velocity field.

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 large-scale instability, in which the total (kinetic plus magnetic) turbulent pressure can be suppressed in the presence of a weak large-scale 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 super-equipartition 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 pre-emergence phase of a sunspot and to constrain its deep structure and formation mechanism.

Recent direct numerical simulations (DNS) of large-scale turbulent dynamos in strongly stratified layers have resulted in surprisingly sharp bipolar structures at the surface. Here we present new DNS of helically and non-helically forced turbulence with and without rotation and compare with corresponding mean-field simulations (MFS) to show that these structures are a generic outcome of a broader class of dynamos in density-stratified 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 large-scale dynamo-driven magnetic field if the dynamo number is at least 2.5 times supercritical.

Analytic solutions of the mean-field 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 one-dimensional 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 time-independent energy density. On one end where the function value vanishes, the second derivative is finite, which would not be correctly reproduced with sine-like 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 two-scale 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 one-dimensional magnetic field and then two-dimensional slices from a three-dimensional alpha effect-type dynamo-generated 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 (2161--2163) 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 three-dimensional turbulence simulations of helicity-driven alpha effect-type dynamos.

We perform numerical simulations of decaying hydrodynamic and magnetohydrodynamic turbulence. We classify our time-dependent 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 self-similar 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 scale-invariant 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 mean-field simulations. Finally, the relevance to sunspot formation is discussed.

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

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.

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.

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.

(abidged) Context: Stellar convection zones are characterized by vigorous high-Reynolds number turbulence at low Prandtl numbers. Aims: We study the dynamo and differential rotation regimes at varying levels of viscous, thermal, and magnetic diffusion. Methods: We perform three-dimensional simulations of stratified fully compressible magnetohydrodynamic convection in rotating spherical wedges at various thermal and magnetic Prandtl numbers (from 0.25 to 2 and 5, respectively). Results: We find that the rotation profiles for high thermal diffusivity show a monotonically increasing angular velocity from the bottom of the convection zone to the top and from the poles toward the equator. For sufficiently rapid rotation, a region of negative radial shear develops at mid-latitudes as the thermal diffusivity is decreased. This coincides with a change in the dynamo mode from poleward propagating activity belts to equatorward propagating ones. Furthermore, the cyclic solutions disappear at the highest magnetic Reynolds numbers. The total magnetic energy increases with the magnetic Reynolds number in the range studied here ($5-151$), but the energies of the mean magnetic fields level off at high magnetic Reynolds numbers. The differential rotation is strongly affected by the magnetic fields and almost vanishes at the highest magnetic Reynolds numbers. In some of our most turbulent cases we find that two regimes are possible where either differential rotation is strong and mean magnetic fields relatively weak or vice versa. Conclusions: Our simulations indicate a strong non-linear feedback of magnetic fields on differential rotation, leading to qualitative changes in the behaviors of large-scale dynamos at high magnetic Reynolds numbers. Furthermore, we do not find indications of the simulations approaching an asymptotic regime where the results would be independent of diffusion coefficients.

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.

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.

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.

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.

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

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.

A brief summary of the various observations and constraints that underlie solar dynamo research are presented. The arguments that indicate that the solar dynamo is an alpha-omega dynamo of the Babcock-Leighton type are then shortly reviewed. The main open questions that remain are concerned with the subsurface dynamics, including why sunspots emerge at preferred latitudes as seen in the familiar butterfly wings, why the cycle is about 11 years long, and why the sunspot groups emerge tilted with respect to the equator (Joy's law). Next, we turn to magnetic helicity, whose conservation property has been identified with the decline of large-scale magnetic fields found in direct numerical simulations at large magnetic Reynolds numbers. However, magnetic helicity fluxes through the solar surface can alleviate this problem and connect theory with observations, as will be discussed.

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.

We study the connection between spherical wedge and full spherical shell geometries using simple mean-field $\alpha^2$ dynamos. We solve the equations for a one-dimensional time-dependent mean-field dynamo to examine the effects of varying the polar angle $\theta_0$ between the latitudinal boundaries and the poles in spherical coordinates. We investigate the effects of turbulent magnetic diffusivity and $\alpha$ effect profiles as well as different latitudinal boundary conditions to isolate parameter regimes where oscillatory solutions are found. Finally, we add shear along with a damping term mimicking radial gradients to study the resulting dynamo regimes. We find that the commonly used perfect conductor boundary condition leads to oscillatory $\alpha^2$ dynamo solutions only if the wedge boundary is at least one degree away from the poles. Other boundary conditions always produce stationary solutions. By varying the profile of the turbulent magnetic diffusivity alone, oscillatory solutions are achieved with models extending to the poles, but the magnetic field is strongly concentrated near the poles and the oscillation period is very long. By introducing radial shear and a damping term mimicking radial gradients, we again see oscillatory dynamos, and the direction of drift follows the Parker--Yoshimura rule. Oscillatory solutions in the weak shear regime are found only in the wedge case with $\theta_0 = 1^\circ$ and perfect conductor boundaries. A reduced $\alpha$ effect near the poles with a turbulent diffusivity concentrated toward the equator yields oscillatory dynamos with equatorward migration and reproduces best the solutions in spherical wedges.

(abridged) Context: The mechanisms that cause the formation of sunspots are still unclear. Aims: We study the self-organisation of initially uniform sub-equipartition magnetic fields by highly stratified turbulent convection. Methods: We perform simulations of magnetoconvection in Cartesian domains that are $8.5$-$24$ Mm deep and $34$-$96$ Mm wide. We impose either a vertical or a horizontal uniform magnetic field in a convection-driven turbulent flow. Results: We find that super-equipartition magnetic flux concentrations are formed near the surface with domain depths of $12.5$ and $24$ Mm. The size of the concentrations increases as the box size increases and the largest structures ($20$ Mm horizontally) are obtained in the 24 Mm deep models. The field strength in the concentrations is in the range of $3$-$5$ kG. The concentrations grow approximately linearly in time. The effective magnetic pressure measured in the simulations is positive near the surface and negative in the bulk of the convection zone. Its derivative with respect to the mean magnetic field, however, is positive in the majority of the domain, which is unfavourable for the negative effective magnetic pressure instability (NEMPI). Furthermore, we find that magnetic flux is concentrated in regions of converging flow corresponding to large-scale supergranulation convection pattern. Conclusions: The linear growth of large-scale flux concentrations implies that their dominant formation process is tangling of the large-scale field rather than an instability. One plausible mechanism explaining both the linear growth and the concentrate on of the flux in the regions of converging flow pattern is flux expulsion. Possible reasons for the absence of NEMPI are that the derivative of the effective magnetic pressure with respect to the mean magnetic field has an unfavourable sign and that there may not be sufficient scale separation.

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.

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.

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).

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{\deg} latitude and an azimuthal extent of 180{\deg}. 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 {\Omega}-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{\deg}. 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.

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.

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\'en-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$.

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.

In light of new results, the one-dimensional mean-field dynamo model of Brandenburg & Kapyla (2007) with dynamical quenching and a nonlocal Babcock-Leighton alpha effect is re-examined for the solar dynamo. We extend the one-dimensional model to include the effects of turbulent downward pumping (Kitchatinov & Olemskoy 2011), and to combine dynamical quenching with shear. We use both the conventional dynamical quenching model of Kleeorin & Ruzmaikin (1982) and the alternate one of Hubbard & Brandenburg (2011), and confirm that with varying levels of non-locality in the alpha effect, and possibly shear as well, the saturation field strength can be independent of the magnetic Reynolds number.

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.

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 $\alpha\Omega$ dynamo wave following the Parker--Yoshimura rule. We conclude that the equatorward migration in this and previous work is due to a positive (negative) $\alpha$ effect in the northern (southern) hemisphere and a negative radial gradient of $\Omega$ 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.

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.

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.

(abridged) Context: Solar-like differential rotation is characterized by a rapidly rotating equator and slower poles. However, theoretical models and numerical simulations can result in a slower equator and faster poles when the rotation is slow. Aims: We study the critical rotational influence under which differential rotation flips from solar-like (fast equator, slow poles) to an anti-solar one (slow equator, fast poles). We estimate the non-diffusive ($\Lambda$ effect) and diffusive (turbulent viscosity) contributions to the Reynolds stress. Methods: We present the results of three-dimensional numerical simulations of mildly turbulent convection in spherical wedge geometry. Here we apply a fully compressible setup which would suffer from a prohibitive time step constraint if the real solar luminosity was used. We regulate the convective velocities by varying the amount of heat transported by thermal conduction, turbulent diffusion, and resolved convection. Results: Increasing the efficiency of resolved convection leads to a reduction of the rotational influence on the flow and a sharp transition from solar-like to anti-solar differential rotation for Coriolis numbers around 1.3. We confirm the recent finding of a large-scale flow bistability: contrasted with running the models from an initial condition with unprescribed differential rotation, the initialization of the model with certain kind of rotation profile sustains the solution over a wider parameter range. Conclusions: Our results may have implications for real stars that start their lives as rapid rotators implying solar-like rotation in the early main-sequence evolution. As they slow down, they might be able to retain solar-like rotation for lower Coriolis numbers before switching to anti-solar rotation. This could partially explain the puzzling findings of anti-solar rotation profiles for models in the solar parameter regime.

We perform direct numerical simulations of the equations of magnetohydrodynamics with external random forcing and in the presence of gravity. The domain is divided into two parts: a lower layer where the forcing is helical and an upper layer where the helicity of the forcing is zero with a smooth transition in between. At early times, a large-scale helical dynamo develops in the bottom layer. At later times the dynamo saturates, but the vertical magnetic field continues to develop and rises to form dynamic bipolar structures at the top, which later disappear and reappear. Some of the structures look similar to $\delta$ spots observed in the Sun. This is the first example of magnetic flux concentrations, owing to strong density stratification, from self-consistent dynamo simulations that generate bipolar, super-equipartition strength, magnetic structures whose energy density can exceeds the turbulent kinetic energy by even a factor of ten.

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$).

Hydromagnetic turbulence affects the evolution of large-scale magnetic fields through mean-field effects like turbulent diffusion and the $\alpha$ 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 $\alpha$ 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 \Omega \times \bf J$ effect nearly the same quenching exponent is obtained as for $\alpha$. 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.

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.

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.

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

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

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.

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 $\alpha^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 $\alpha$ 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 $\alpha$ effect dynamo.

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.

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.

We calculate the mean electromotive force in plane Couette flows of a nonrotating conducting fluid under the influence of a large-scale magnetic field for driven turbulence. A vertical stratification of the turbulence intensity results in an alpha effect owing to the presence of horizontal shear. Here we discuss the possibility of an experimental determination of the components of the alpha tensor using both quasilinear theory and nonlinear numerical simulations. For magnetic Prandtl numbers of the order of unity, we find that in the high-conductivity limit the alpha effect in the direction of the flow clearly exceeds the component in spanwise direction. In this limit, alpha runs linearly with the magnetic Reynolds number Rm while in the low-conductivity limit it runs with the product Rm*Re, where Re is the kinetic Reynolds number so that for given Rm the alpha effect grows with decreasing magnetic Prandtl number. For the small magnetic Prandtl numbers of liquid metals, a common value for the horizontal elements of the alpha tensor appears, which makes it unimportant whether the alpha effect is measured in the spanwise or streamwise directions. The resulting effect should lead to an observable voltage in both directions of about 0.5 mV for magnetic fields of 1 kgauss and velocity fluctuations of about 1 m/s in a channel of 50 cm height (independent of its width).

Strongly stratified hydromagnetic turbulence has recently been identified as a candidate for explaining the spontaneous formation of magnetic flux concentrations by the negative effective magnetic pressure instability (NEMPI). Much of this work has been done for isothermal layers for which the density scale height is constant throughout. We now study the validity of earlier conclusions about the size and growth rate of magnetic structures in the case of polytropic layers, which scale height decreases sharply towards the surface. To allow for a continuous transition from isothermal to polytropic layers, we employ a generalization of the exponential function known as the q-exponential. Now, the top of the polytropic layer shifts with the polytropic index such that the scale height at some reference height is always the same. We use both mean-field and direct numerical simulations of forced stratified turbulence to determine the resulting flux concentrations. Magnetic structures begin to form at a depth where magnetic field strength is about 3-4% the local equipartition field strength with respect to the turbulent kinetic energy. Unlike the isothermal case where stronger fields can give rise to magnetic flux concentrations at larger depths, in the polytropic one the growth rates decreases for structures deeper down. For vertical fields, magnetic structures of super-equipartition strengths are formed because such fields survive downward advection, unlike NEMPI under horizontal magnetic fields. The horizontal cross-section of such structures is approximately circular. Results based on isothermal models can be applied locally to polytropic layers. For vertical fields, magnetic flux concentrations of super-equipartition strengths form, which supports suggestions that sunspot formation might be a shallow phenomenon.

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.

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.

Strongly stratified hydromagnetic turbulence has previously been found to produce magnetic flux concentrations if the domain is large enough compared with the size of turbulent eddies. Mean-field simulations (MFS) using parameterizations of the Reynolds and Maxwell stresses show a negative effective magnetic pressure instability and have been able to reproduce many aspects of direct numerical simulations (DNS) regarding the growth rate of this large-scale instability, shape of the resulting magnetic structures, and their height as a function of magnetic field strength. Unlike the case of an imposed horizontal field, for a vertical one, magnetic flux concentrations of equipartition strength with the turbulence can be reached. This results in magnetic spots that are reminiscent of sunspots. Here we want to find out under what conditions magnetic flux concentrations with vertical field occur and what their internal structure is. We use a combination of MFS, DNS, and implicit large-eddy simulations to characterize the resulting magnetic flux concentrations in forced isothermal turbulence with an imposed vertical magnetic field. We confirm earlier results that in the kinematic stage of the large-scale instability the horizontal wavelength of structures is about 10 times the density scale height. At later times, even larger structures are being produced in a fashion similar to inverse spectral transfer in helically driven turbulence. Using turbulence simulations, we find that magnetic flux concentrations occur for different values of the Mach number between 0.1 and 0.7. DNS and MFS show magnetic flux tubes with mean-field energies comparable to the turbulent kinetic energy. The resulting vertical magnetic flux tubes are being confined by downflows along the tubes and corresponding inflow from the sides, which keep the field concentrated.

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.

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 $\pi/2$, which can be explained by an oscillatory $\alpha^2$ dynamo caused by the sign change of the $\alpha$-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\pi$ azimuthal extent, reminiscent of the field configurations deduced from observations of rapidly rotating late-type stars.

We present results from compressible Cartesian convection simulations with and without imposed shear. In the former case the dynamo is expected to be of $\alpha^2\varOmega$ type which is generally expected to be relevant for the Sun, whereas the latter case refers to $\alpha^2$ dynamos which are more likely to occur in more rapidly rotating stars whose differential rotation is small. We perform a parameter study where the shear flow and the rotational influence are varied to probe the relative importance of both types of dynamos. Oscillatory solutions are preferred both in the kinematic and saturated regimes when the negative ratio of shear to rotation rates, $q\equiv -S/\varOmega$, is between 1.5 and 2, i.e., when shear and rotation are of comparable strengths. Other regions of oscillatory solutions are found with small values of $q$, i.e., when shear is weak in comparison to rotation, and in the regime of large negative $q$s, when shear is very strong in comparison to rotation. However, exceptions to these rules also appear so that for a given ratio of shear to rotation, solutions are non-oscillatory for small and large shear, but oscillatory in the intermediate range. Changing the boundary conditions from vertical field to perfect conductor ones changes the dynamo mode from oscillatory to quasi-steady. Furthermore, in many cases an oscillatory solution exists only in the kinematic regime whereas in the nonlinear stage the mean fields are stationary. However, the cases with rotation and no shear are always oscillatory in the parameter range studied here and the dynamo mode does not depend on the magnetic boundary conditions. The strengths of total and large-scale components of the magnetic field in the saturated state, however, are sensitive to the chosen boundary conditions.

Magnetic field generation on scales large compared with the scale of the turbulent eddies is known to be possible via the so-called $\alpha$ effect when the turbulence is helical and if the domain is large enough for the $\alpha$ 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.

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.

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.

In astrophysics, turbulent diffusion is often used in place of microphysical diffusion to avoid resolving the small scales. However, we expect this approach to break down when time and length scales of the turbulence become comparable with other relevant time and length scales in the system. Turbulent diffusion has previously been applied to the magneto-rotational instability (MRI), but no quantitative comparison of growth rates at different turbulent intensities has been performed. We investigate to what extent turbulent diffusion can be used to model the effects of small-scale turbulence on the kinematic growth rates of the MRI, and how this depends on angular velocity and magnetic field strength. We use direct numerical simulations in three-dimensional shearing boxes with periodic boundary conditions in the spanwise direction and additional random plane-wave volume forcing to drive a turbulent flow at a given length scale. We estimate the turbulent diffusivity using a mixing length formula and compare with results obtained with the test-field method. It turns out that the concept of turbulent diffusion is remarkably accurate in describing the effect of turbulence on the growth rate of the MRI. No noticeable breakdown of turbulent diffusion has been found, even when time and length scales of the turbulence become comparable with those imposed by the MRI itself. On the other hand, quenching of turbulent magnetic diffusivity by the magnetic field is found to be absent. Turbulence reduces the growth rate of the MRI in a way that is the same as microphysical magnetic diffusion.

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.

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.

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.

Recent progress in astrophysical hydromagnetic turbulence is being reviewed. The physical ideas behind the now widely accepted Goldreich--Sridhar model and its extension to compressible magnetohydrodynamic turbulence are introduced. Implications for cosmic ray diffusion and acceleration is being discussed. Dynamo-generated magnetic fields with and without helicity are contrasted against each other. Certain turbulent transport processes are being modified and often suppressed by anisotropy and inhomogeneities of the turbulence, while others are being produced by such properties, which can lead to new large-scale instabilities of the turbulent medium. Applications of various such processes to astrophysical systems are being considered.

We consider mean-field dynamo models with fluctuating \alpha effect, both with and without shear. The \alpha 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 \alpha^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.

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.

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

In a strongly stratified turbulent layer, a uniform horizontal magnetic field can become unstable to spontaneously form local flux concentrations due to a negative contribution of turbulence to the large-scale (mean-field) magnetic pressure. This mechanism, called the negative effective magnetic pressure instability (NEMPI), is of interest in connection with dynamo scenarios where most of the magnetic field resides in the bulk of the convection zone, and not at the bottom. Recent work using the mean-field hydromagnetic equations has shown that NEMPI becomes suppressed at rather low rotation rates with Coriolis numbers as low as 0.1.}{Here we extend these earlier investigations by studying the effects of rotation both on the development of NEMPI and on the effective magnetic pressure. We also quantify the kinetic helicity from direct numerical simulations (DNS) and compare with earlier work.}{To calculate the rotational effect on the effective magnetic pressure we consider both DNS and analytical studies using the $\tau$ approach. To study the effects of rotation on the development of NEMPI we use both DNS and mean-field calculations of the 3D hydromagnetic equations in a Cartesian domain.}{We find that the growth rates of NEMPI from earlier mean-field calculations are well reproduced with DNS, provided the Coriolis number is below about 0.06. In that case, kinetic and magnetic helicities are found to be weak. For faster rotation, dynamo action becomes possible. However, there is an intermediate range of rotation rates where dynamo action on its own is not yet possible, but the rotational suppression of NEMPI is being alleviated.}{Production of magnetic flux concentrations through the suppression of turbulent pressure appears to be possible only in the upper-most layers of the Sun, where the convective turnover time is less than 2 hours.}

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.

Energetic nonthermal particles (cosmic rays, CRs) are accelerated in supernova remnants, relativistic jets and other astrophysical objects. The CR energy density is typically comparable with that of the thermal components and magnetic fields. In this review we discuss mechanisms of magnetic field amplification due to instabilities induced by CRs. We derive CR kinetic and magnetohydrodynamic equations that govern cosmic plasma systems comprising the thermal background plasma, comic rays and fluctuating magnetic fields to study CR-driven instabilities. Both resonant and non-resonant instabilities are reviewed, including the Bell short-wavelength instability, and the firehose instability. Special attention is paid to the longwavelength instabilities driven by the CR current and pressure gradient. The helicity production by the CR current-driven instabilities is discussed in connection with the dynamo mechanisms of cosmic magnetic field amplification.

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.

A central quantity in mean-field magnetohydrodynamics is the mean electromotive force EMF, which in general depends on the mean magnetic field. It may however have a part independent of the mean magnetic field. Here we study an example of a rotating conducting body of turbulent fluid with non-zero cross-helicity, in which a contribution to the EMF proportional to the angular velocity occurs (Yoshizawa 1990). If the forcing is helical, it also leads to an alpha effect, and large-scale magnetic fields can be generated. For not too rapid rotation, the field configuration is such that Yoshizawa's contribution to the EMF is considerably reduced compared to the case without alpha effect. In that case, large-scale flows are also found to be generated.

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.

The surface layers of the Sun are strongly stratified. In the presence of turbulence with a weak mean magnetic field, a large-scale instability resulting in the formation of non-uniform magnetic structures, can be excited over the scale of many turbulent eddies or convection cells. This instability is caused by a negative contribution of turbulence to the effective (mean-field) magnetic pressure and has previously been discussed in connection with the formation of active regions and perhaps sunspots. We want to understand the effects of rotation on this instability in both two and three dimensions. We use mean-field magnetohydrodynamics in a parameter regime in which the properties of the negative effective magnetic pressure instability have previously been found to be in agreement with those of direct numerical simulations. We find that the instability is suppressed already for relatively slow rotation with Coriolis numbers (i.e. inverse Rossby numbers) around 0.2. The suppression is strongest at the equator. In the nonlinear regime, we find traveling wave solutions with propagation in the prograde direction at the equator with additional poleward migration away from the equator. The prograde rotation of the magnetic pattern near the equator is argued to be a possible explanation for the faster rotation speed of magnetic tracers found on the Sun. In the bulk of the domain, kinetic and current helicities are negative in the northern hemisphere and positive in the southern.

Using simulations of slowly rotating stratified turbulence, we show that the alpha effect responsible for the generation of astrophysical magnetic fields is proportional to the logarithmic gradient of kinetic energy density rather than that of momentum, as was previously thought. This result is in agreement with a new analytic theory developed in this paper for large Reynolds numbers. Thus, the contribution of density stratification is less important than that of turbulent velocity. The alpha effect and other turbulent transport coefficients are determined by means of the test-field method. In addition to forced turbulence, we also investigate supernova-driven turbulence and stellar convection. In some cases (intermediate rotation rate for forced turbulence, convection with intermediate temperature stratification, and supernova-driven turbulence) we find that the contribution of density stratification might be even less important than suggested by the analytic theory.

We present preliminary results from high resolution magneto-convection simulations where we find the formation of flux concentrations from an initially uniform magnetic field. We compute the effective magnetic pressure but find that the concentrations appear also in places where it is positive. The structures appear in roughly ten convective turnover times and live close to a turbulent diffusion time. The time scales are compatible with the negative effective magnetic pressure instability (NEMPI), although structure formation is not restricted to regions where the effective magnetic pressure is negative.

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 $\alpha^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.

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.

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.

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.

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.

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

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.

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.

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.

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.

We investigate the effects of weakly and strongly stratified turbulent convection on the mean effective Lorentz force, and especially on the mean effective magnetic pressure. Earlier studies with isotropically forced non-stratified and stratified turbulence have shown that the contribution of the turbulence to the mean magnetic pressure is negative for mean horizontal magnetic fields that are smaller than the equipartition strength, so that the effective mean magnetic pressure that takes into account the turbulence effects, can be negative. Compared with earlier cases of forced turbulence with an isothermal equation of state, we find that the turbulence effect is similar to or even stronger in the present case of turbulent convection. This is argued to be due to the anisotropy of turbulence in the vertical direction. Another important difference compared with earlier studies is the presence of an evolution equation for the specific entropy. Mean-field modelling with entropy evolution indicates that the negative effective magnetic pressure can still lead to a large-scale instability which forms local flux concentrations, even though the specific entropy evolution tends to have a stabilizing effect when applied to a stably stratified (e.g., isothermal) layer. It is argued that this large-scale instability could be important for the formation of solar large-scale magnetic structures such as active regions.

We report on simulations of mildly turbulent convection in spherical wedge geometry with varying density stratification. We vary the density contrast within the convection zone by a factor of 20 and study the influence of rotation on the solutions. We demonstrate that the size of convective cells decreases and the anisotropy of turbulence increases as the stratification is increased. Differential rotation is found to change from anti-solar (slow equator) to solar-like (fast equator) at roughly the same Coriolis number for all stratifications. The largest stratification runs, however, are sensitive to changes of the Reynolds number. Evidence for a near-surface shear layer is found in runs with strong stratification and large Reynolds numbers.

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.

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.

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.

We study the role of viscosity and the effects of a magnetic field on a rotating, self-gravitating fluid, using Newtonian theory and adopting the ideal magnetohydrodynamic approximation. Our results confirm that viscosity can generate vorticity in inhomogeneous environments, while the magnetic tension can produce vorticity even in the absence of fluid pressure and density gradients. Linearizing our equations around an Einstein-de Sitter cosmology, we find that viscosity adds to the diluting effect of the universal expansion. Typically, however, the dissipative viscous effects are confined to relatively small scales. We also identify the characteristic length bellow which the viscous dissipation is strong and beyond which viscosity is essentially negligible. In contrast, magnetism seems to favor cosmic rotation. The magnetic presence is found to slow down the standard decay-rate of linear vortices, thus leading to universes with more residual rotation than generally anticipated.

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.

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.

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.

We investigate the mean electromotive force in the kinematic framework, that is, ignoring the back-reaction of the magnetic field on the fluid velocity, under the assumption of axisymmetric turbulence determined by the presence of either rotation, density stratification, or both. We use an analogous approach for the mean passive scalar flux. As an alternative to convection, we consider forced turbulence in an isothermal layer. When using standard ansatzes, the mean magnetic transport is then determined by nine, and the mean passive scalar transport by four coefficients. We give results for all these transport coefficients. We use the test-field method and the test-scalar method, where transport coefficients are determined by solving sets of equations with properly chosen mean magnetic fields or mean scalars. These methods are adapted to mean fields which may depend on all three space coordinates. We find the anisotropy of turbulent diffusion to be moderate in spite of rapid rotation or strong density stratification. Contributions to the mean electromotive force determined by the symmetric part of the gradient tensor of the mean magnetic field, which were ignored in several earlier investigations, turn out to be important. In stratified rotating turbulence, the $\alpha$ effect is strongly anisotropic, suppressed along the rotation axis on large length scales, but strongly enhanced at intermediate length scales. Also the $\OO\times\meanJJ$ effect is enhanced at intermediate length scales. The turbulent passive scalar diffusivity is typically almost twice as large as the turbulent magnetic diffusivity. Both magnetic and passive scalar diffusion are slightly enhanced along the rotation axis, but decreased if there is gravity.

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.

As was demonstrated in earlier studies, turbulence can result in a negative contribution to the effective mean magnetic pressure, which, in turn, can cause a large-scale instability. In this study, hydromagnetic mean-field modelling is performed for an isothermally stratified layer in the presence of a horizontal magnetic field. The negative effective magnetic pressure instability (NEMPI) is comprehensively investigated. It is shown that, if the effect of turbulence on the mean magnetic tension force vanishes, which is consistent with results from direct numerical simulations of forced turbulence, the fastest growing eigenmodes of NEMPI are two-dimensional. The growth rate is found to depend on a parameter beta_* characterizing the turbulent contribution of the effective mean magnetic pressure for moderately strong mean magnetic fields. A fit formula is proposed that gives the growth rate as a function of turbulent kinematic viscosity, turbulent magnetic diffusivity, the density scale height, and the parameter beta_*. The strength of the imposed magnetic field does not explicitly enter provided the location of the vertical boundaries are chosen such that the maximum of the eigenmode of NEMPI fits into the domain. The formation of sunspots and solar active regions is discussed as possible applications of NEMPI.

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.

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 $\alpha\Omega$ dynamo we find that the preferred mode is antisymmetric about the mid-plane. This is also verified in 3-D direct numerical simulations.

Small levels of turbulence can be present in stellar radiative interiors due to, e.g., instability of rotational shear. In this paper we estimate turbulent transport coefficients for stably stratified rotating stellar radiation zones. Stable stratification induces strong anisotropy with a very small ratio of radial-to-horizontal turbulence intensities. Angular momentum is transported mainly due to the correlation between azimuthal and radial turbulent motions induced by the Coriolis force. This non-diffusive transport known as the Lambda-effect has outward direction in radius and is much more efficient compared to the effect of radial eddy viscosity. Chemical species are transported by small radial diffusion only. This result is confirmed using direct numerical simulations combined with the test-scalar method. As a consequence of the non-diffusive transport of angular momentum, the estimated characteristic time of rotational coupling (< 100 Myr) between radiative core and convective envelope in young solar-type stars is much shorter compared to the time-scale of Lithium depletion (~1 Gyr).

Direct numerical simulations of isotropically forced homogeneous stationary turbulence with an imposed passive scalar concentration gradient are compared with an analytical closure model which provides evolution equations for the mean passive scalar flux and variance. Triple correlations of fluctuations appearing in these equations are described in terms of relaxation terms proportional to the quadratic correlations. Three methods are used to extract the relaxation timescales tau_i from direct numerical simulations. Firstly, we insert the closure ansatz into our equations, assume stationarity, and solve for tau_i. Secondly, we use only the closure ansatz itself and obtain tau_i from the ratio of quadratic and triple correlations. Thirdly we remove the imposed passive scalar gradient and fit an exponential decay law to the solution. We vary the Reynolds (Re) and P\'eclet (Pe) numbers while keeping their ratio at unity and the degree of scale separation and find for large Re fair correspondence between the different methods. The ratio of the turbulent relaxation time of passive scalar flux to the turnover time of turbulent eddies is of the order of three, which is in remarkable agreement with earlier work. Finally we make an effort to extract the relaxation timescales relevant for the viscous and diffusive effects. We find two regimes which are valid for small and large Re, respectively, but the dependence of the parameters on scale separation suggests that they are not universal.

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.

The phenomenon of turbulent thermal diffusion in temperature-stratified turbulence causing a non-diffusive turbulent flux of inertial and non-inertial particles in the direction of the turbulent heat flux is found using direct numerical simulations (DNS). In simulations with and without gravity, this phenomenon is found to cause a peak in the particle number density around the minimum of the mean fluid temperature for Stokes numbers less than 1, where the Stokes number is the ratio of particle Stokes time to turbulent Kolmogorov time at the viscous scale. Turbulent thermal diffusion causes the formation of large-scale inhomogeneities in the spatial distribution of inertial particles. The strength of this effect is maximum for Stokes numbers around unity, and decreases again for larger values. The dynamics of inertial particles is studied using Lagrangian modelling in forced temperature-stratified turbulence, whereas non-inertial particles and the fluid are described using DNS in an Eulerian framework.

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.

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.

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\`eclet 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.

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.

When scale separation in space and time is poor, the alpha effect and turbulent diffusivity have to be replaced by integral kernels. Earlier work in computing these kernels using the test-field method is now generalized to the case in which both spatial and temporal scale separations are poor. The approximate form of the kernel is such that it can be treated in a straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean-field equations are solved for oscillatory alpha-shear dynamos as well as alpha^2 dynamos in which alpha is antisymmetric about the equator, making this dynamo also oscillatory. In both cases, the critical values of the dynamo number is lowered by the fact that the dynamo is oscillatory.

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.

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.

We determine the timescales associated with turbulent diffusion and isotropization in closure models using anisotropically forced and freely decaying turbulence simulations and to study the applicability of these models. We compare the results from anisotropically forced three-dimensional numerical simulations with the predictions of the closure models and obtain the turbulent timescales mentioned above as functions of the Reynolds number. In a second set of simulations, turning the forcing off enables us to study the validity of the closures in freely decaying turbulence. Both types of experiments suggest that the timescale of turbulent diffusion converges to a constant value at higher Reynolds numbers. Furthermore, the relative importance of isotropization is found to be about 2.5 times larger at higher Reynolds numbers than in the more viscous regime.

Using different analytical methods (the quasi-linear approach, the path-integral technique and tau-relaxation approximation) we develop a comprehensive mean-field theory for a pumping effect of the mean magnetic field in homogeneous non-rotating helical turbulence with imposed large-scale shear. The effective pumping velocity is proportional to the product of alpha effect and large-scale vorticity associated with the shear, and causes a separation of the toroidal and poloidal components of the mean magnetic field along the direction of the mean vorticity. We also perform direct numerical simulations of sheared turbulence in different ranges of hydrodynamic and magnetic Reynolds numbers and use a kinematic test-field method to determine the effective pumping velocity. The results of the numerical simulations are in agreement with the theoretical predictions.

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.

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.

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.

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.

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.

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.

We extend earlier models of turbulent dynamos with an upper, nearly force-free exterior to spherical geometry, and study how flux emerges from lower layers to the upper ones without being driven by magnetic buoyancy. We also study how this affects the possibility of plasmoid ejection. A spherical wedge is used that includes northern and southern hemispheres up to mid-latitudes and a certain range in longitude of the Sun. In radius, we cover both the region that corresponds to the convection zone in the Sun and the immediate exterior up to twice the radius of the Sun. Turbulence is driven with a helical forcing function in the interior, where the sign changes at the equator between the two hemispheres. An oscillatory large-scale dynamo with equatorward migration is found to operate in the turbulence zone. Plasmoid ejections occur in regular intervals, similar to what is seen in earlier Cartesian models. These plasmoid ejections are tentatively associated with coronal mass ejections. The magnetic helicity is found to change sign outside the turbulence zone, which is in agreement with recent findings for the solar wind.

We simulate the magnetic feature tracking (MFT) speed using advective-diffusive transport models in both one and two dimensions. By depositing magnetic bipolar regions at different latitudes at the Sun's surface and following their evolution for a prescribed meridional circulation and magnetic diffusivity profiles, we derive the MFT speed as a function of latitude. We find that in a one dimensional surface-transport model the simulated MFT speed at the surface is always the same as the meridional flow-speed used as input to the model, but is different in a two-dimensional transport model in the meridional (r,theta) plane. The difference depends on the value of the magnetic diffusivity and on the radial gradient of the latitudinal velocity. We have confirmed our results with two different codes in spherical and Cartesian coordinates.

Context. Earlier work has suggested that large-scale dynamos can reach and maintain equipartition field strengths on a dynamical time scale only if magnetic helicity of the fluctuating field can be shed from the domain through open boundaries. Aims. Our aim is to test this scenario in convection-driven dynamos by comparing results for open and closed boundary conditions. Methods. Three-dimensional numerical simulations of turbulent compressible convection with shear and rotation are used to study the effects of boundary conditions on the excitation and saturation of large-scale dynamos. Open (vertical-field) and closed (perfect-conductor) boundary conditions are used for the magnetic field. The shear flow is such that the contours of shear are vertical, crossing the outer surface, and are thus ideally suited for driving a shear-induced magnetic helicity flux. Results. We find that for given shear and rotation rate, the growth rate of the magnetic field is larger if open boundary conditions are used. The growth rate first increases for small magnetic Reynolds number, Rm, but then levels off at an approximately constant value for intermediate values of Rm. For large enough Rm, a small-scale dynamo is excited and the growth rate of the field in this regime increases as Rm^(1/2). Regarding the nonlinear regime, the saturation level of the energy of the total magnetic field is independent of Rm when open boundaries are used. In the case of perfect conductor boundaries, the saturation level first increases as a function of Rm, but then decreases proportional to Rm^(-1) for Rm > 30, indicative of catastrophic quenching. These results suggest that the shear-induced magnetic helicity flux is efficient in alleviating catastrophic quenching when open boundaries are used. The horizontally averaged mean field is still weakly decreasing as a function of Rm even for open boundaries.

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.

(abridged) Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Here we perform a systematic comparison between Cartesian and spherical geometries as a function of the rotation rate. Aims. We extend the earlier studies by using spherical wedges to obtain turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime. In particular, we want to clarify whether the sharp equatorial profile of the horizontal Reynolds stress found in earlier Cartesian models is reproduced in spherical models. Methods. We employ direct numerical simulations of turbulent convection. In order to reach as high spatial resolution as possible in the spherical runs, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong `banana cells'. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor--Proudman balance.

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.

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.

Using simulations of helically driven turbulence, it is shown that the ratio of kinetic to magnetic energy dissipation scales with the magnetic Prandtl number in power law fashion with an exponent of approximately 0.6. Over six orders of magnitude in the magnetic Prandtl number the magnetic field is found to be sustained by large-scale dynamo action of alpha-squared type. This work extends a similar finding for small magnetic Prandtl numbers to the regime of large magnetic Prandtl numbers. At large magnetic Prandtl numbers, most of the energy is dissipated viscously, lowering thus the amount of magnetic energy dissipation, which means that simulations can be performed at magnetic Reynolds numbers that are large compared to the usual limits imposed by a given resolution. This is analogous to an earlier finding that at small magnetic Prandtl numbers, most of the energy is dissipated resistively, lowering the amount of kinetic energy dissipation, so simulations can then be performed at much larger fluid Reynolds numbers than otherwise. The decrease in magnetic energy dissipation at large magnetic Prandtl numbers is discussed in the context of underluminous accretion found in some quasars.

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.

Convectively driven dynamos with rotation generating magnetic fields on scales large compared with the scale of the turbulent eddies are being reviewed. It is argued that such fields can be understood as the result of an alpha effect. Simulations in Cartesian domains show that such large-scale magnetic fields saturate on a time scale compatible with the resistive one, suggesting that the magnitude of the alpha effect is here still constrained by approximate magnetic helicity conservation. It is argued that, in the absence of shear and/or any other known large-scale dynamo effects, these simulations prove the existence of turbulent alpha^2-type dynamos. Finally, recent results are discussed in the context of solar and stellar 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.

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.

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.

Context. Helically forced magneto-hydrodynamic shearing-sheet turbulence can support different large-scale dynamo modes, although the {\alpha} - {\Omega} mode is generally expected to dominate because it is the fastest growing. In an {\alpha} - {\Omega} 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 {\alpha} - {\Omega} 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 {\alpha} - {\Omega} and {\alpha}2 dynamo action is possible. Results. After initial dominance and saturation, the {\alpha} - {\Omega} mode is found to be destroyed by an orthogonal {\alpha}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 {\alpha} - {\Omega} dynamos by its {\alpha}2 sibling is discussed in the context of grand minima of solar and stellar activity. However, the genesis of {\alpha} - {\Omega} dynamos from an {\alpha}2 dynamo has not yet been found.

We study the ability of magnetic helicity expulsion to alleviate catastrophic $\alpha$-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 $\alpha^2\Omega$ mean field models with a dynamical $\alpha$ 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.

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.

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.

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.

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.

In a density-stratified turbulent medium the cross helicity <u'.B'> is considered as a result of the interaction of the velocity fluctuations and a large-scale magnetic field. By means of a quasilinear theory and by numerical simulations we find the cross helicity and the mean vertical magnetic field anti-correlated. In the high-conductivity limit the ratio of the helicity and the mean magnetic field equals the ratio of the magnetic eddy diffusivity and the (known) density scale height. The result can be used to predict that the cross helicity at the solar surface exceeds the value of 1 Gauss km/s. Its sign is anti-correlated with that of the radial mean magnetic field. Alternatively, we can use our result to determine the value of the turbulent magnetic diffusivity from observations of the cross helicity.

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.

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.

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.

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.

While the rising flux tube paradigm is an elegant theory, its basic assumptions, thin flux tubes at the bottom of the convection zone with field strengths two orders of magnitude above equipartition, remain numerically unverified at best. As such, in recent years the idea of a formation of sunspots near the top of the convection zone has generated some interest. The presence of turbulence can strongly enhance diffusive transport mechanisms, leading to an effective transport coefficient formalism in the mean-field formulation. The question is what happens to these coefficients when the turbulence becomes anisotropic due to a strong large-scale mean magnetic field. It has been noted in the past that this anisotropy can also lead to highly non-diffusive behaviour. In the present work we investigate the formation of large-scale magnetic structures as a result of a negative contribution of turbulence to the large-scale effective magnetic pressure in the presence of stratification. In direct numerical simulations of forced turbulence in a stratified box, we verify the existence of this effect. This phenomenon can cause formation of large-scale magnetic structures even from initially uniform large-scale magnetic field.

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.

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.

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.

Direct numerical simulations and mean-field theory are used to model reactive front propagation in a turbulent medium. In the mean-field approach, memory effects of turbulent diffusion are taken into account to estimate the front speed in cases when the Damkohler number is large. This effect is found to saturate the front speed to values comparable with the speed of the turbulent motions. By comparing with direct numerical simulations, it is found that the effective correlation time is much shorter than for non-reacting flows. The nonlinearity of the reaction term is found to make the front speed slightly faster.

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.

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.

Twisted magnetic fields are frequently seen to emerge above the visible surface of the Sun. This emergence is usually associated with the rise of buoyant magnetic flux structures. Here we ask how magnetic fields from a turbulent large-scale dynamo appear above the surface if there is no magnetic buoyancy. The computational domain is split into two parts. In the lower part, which we refer to as the turbulence zone, the flow is driven by an assumed helical forcing function leading to dynamo action. Above this region, which we refer to as the exterior, a nearly force-free magnetic field is computed at each time step using the stress-and-relax method. Twisted arcade-like field structures are found to emerge in the exterior above the turbulence zone. Strong current sheets tend to form above the neutral line, where the vertical field component vanishes. Time series of the magnetic field structure show recurrent plasmoid ejections. The degree to which the exterior field is force free is estimated as the ratio of the dot product of current density and magnetic field strength to their respective rms values. This ratio reaches values of up to 95% in the exterior. A weak outward flow is driven by the residual Lorentz force.

Using three-dimensional convection simulations it is shown that a sinusoidal variation of horizontal shear leads to a kinematic \alpha effect with a similar sinusoidal variation. The effect exists even for weak stratification and arises owing to the inhomogeneity of turbulence and the presence of impenetrable vertical boundaries. This system produces large-scale magnetic fields that also show a sinusoidal variation in the streamwise direction. It is argued that earlier investigations overlooked these phenomena partly because of the use of horizontal averaging and also because measurements of \alpha using an imposed field combined with long time averages give erroneous results. It is demonstrated that in such cases the actual horizontally averaged mean field becomes non-uniform. The turbulent magnetic diffusion term resulting from such non-uniform fields can then no longer be neglected and begins to balance the \alpha effect.

Self-consistent convective dynamo simulations in wedge-shaped spherical shells are presented. Differential rotation is generated by the interaction of convection with rotation. Equatorward acceleration and dynamo action are obtained only for sufficiently rapid rotation. The angular velocity tends to be constant along cylinders. Oscillatory large-scale fields are found to migrate in the poleward direction. Comparison with earlier simulations in full spherical shells and Cartesian domains is made.

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.

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.

A high-order scheme for direct numerical simulations of turbulent combustion is discussed. Its implementation in the massively parallel and publicly available Pencil Code is validated with the focus on hydrogen combustion. Ignition delay times (0D) and laminar flame velocities (1D) are calculated and compared with results from the commercially available Chemkin code. The scheme is verified to be fifth order in space. Upon doubling the resolution, a 32-fold increase in the accuracy of the flame front is demonstrated. Finally, also turbulent and spherical flame front velocities are calculated and the implementation of the non-reflecting so-called Navier-Stokes Characteristic Boundary Condition is validated in all three directions.

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.

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.

The predictive power of mean-field theory is emphasized by comparing theory with simulations under controlled conditions. The recently developed test-field method is used to extract turbulent transport coefficients both in kinematic as well as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic method is provided by magnetic buoyancy-driven flows that produce an alpha effect and turbulent diffusion.

We present the results of simulations of forced turbulence in a slab where the mean kinetic helicity has a maximum near the mid-plane, generating gradients of magnetic helicity of both large and small-scale fields. We also study systems that have poorly conducting buffer zones away from the midplane in order to assess the effects of boundaries. The dynamical alpha quenching phenomenology requires that the magnetic helicity in the small-scale fields approaches a nearly static, gauge independent state. To stress-test this steady state condition we choose a system with a uniform sign of kinetic helicity, so that the total magnetic helicity can reach a steady state value only through fluxes through the boundary, which are themselves suppressed by the velocity boundary conditions. Even with such a set up, the small-scale magnetic helicity is found to reach a steady state. In agreement with earlier work, the magnetic helicity fluxes of small-scale fields are found to be turbulently diffusive. By comparing results with and without halos, we show that artificial constraints on magnetic helicity at the boundary do not have a significant impact on the evolution of the magnetic helicity, except that "softer" (halo) boundary conditions give a lower energy of the saturated mean magnetic field.

Aims: The test-field method for computing turbulent transport coefficients from simulations of hydromagnetic flows is extended to the regime with a magnetohydrodynamic (MHD) background. Methods: A generalized set of test equations is derived using both the induction equation and a modified momentum equation. By employing an additional set of auxiliary equations, we derive linear equations describing the response of the system to a set of prescribed test fields. Purely magnetic and MHD backgrounds are emulated by applying an electromotive force in the induction equation analogously to the ponderomotive force in the momentum equation. Both forces are chosen to have Roberts flow-like geometry. Results: Examples with an MHD background are studied where the previously used quasi-kinematic test-field method breaks down. In cases with homogeneous mean fields it is shown that the generalized test-field method produces the same results as the imposed-field method, where the field-aligned component of the actual electromotive force from the simulation is used. Furthermore, results for the turbulent diffusivity tensor are given, which are inaccessible to the imposed-field method. For MHD backgrounds, new mean-field effects are found that depend on the occurrence of cross-correlations between magnetic and velocity fluctuations. For strong imposed fields, $\alpha$ is found to be quenched proportional to the fourth power of the field strength, regardless of the type of background studied.

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.

In the mean-field theory of magnetic fields, turbulent transport, i.e. the turbulent electromotive force, is described by a combination of the alpha effect and turbulent magnetic diffusion, which are usually assumed to be proportional respectively to the mean field and its spatial derivatives. For a passive scalar there is just turbulent diffusion, where the mean flux of concentration depends on the gradient of the mean concentration. However, these proportionalities are approximations that are valid only if the mean field or the mean concentration vary slowly in time. Examples are presented where turbulent transport possesses memory, i.e. where it depends crucially on the past history of the mean field. Such effects are captured by replacing turbulent transport coefficients with time integral kernels, resulting in transport coefficients that depend effectively on the frequency or the growth rate of the mean field itself. In this paper we perform numerical experiments to find the characteristic timescale (or memory length) of this effect as well as simple analytical models of the integral kernels in the case of passive scalar concentrations and kinematic dynamos. The integral kernels can then be used to find self-consistent growth or decay rates of the mean fields. In mean-field dynamos the growth rates and cycle periods based on steady state values of alpha effect and turbulent diffusivity can be quite different from the actual values.

Estimates for the nonlinear alpha effect in helical turbulence with an applied magnetic field are presented using two different approaches: the imposed-field method where the electromotive force owing to the applied field is used, and the test-field method where separate evolution equations are solved for a set of different test fields. Both approaches agree for stronger fields, but there are apparent discrepancies for weaker fields that can be explained by the influence of dynamo-generated magnetic fields on the scale of the domain that are referred to as meso-scale magnetic fields. Examples are discussed where these meso-scale fields can lead to both drastically overestimated and underestimated values of alpha compared with the kinematic case. It is demonstrated that the kinematic value can be recovered by resetting the fluctuating magnetic field to zero in regular time intervals. It is concluded that this is the preferred technique both for the imposed-field and the test-field methods.

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.

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.

In mean-field magnetohydrodynamics the mean electromotive force due to velocity and magnetic field fluctuations plays a crucial role. In general it consists of two parts, one independent of and another one proportional to the mean magnetic field. The first part may be nonzero only in the presence of mhd turbulence, maintained, e.g., by small-scale dynamo action. It corresponds to a battery, which lets a mean magnetic field grow from zero to a finite value. The second part, which covers, e.g., the alpha effect, is important for large-scale dynamos. Only a few examples of the aforementioned first part of mean electromotive force have been discussed so far. It is shown that a mean electromotive force proportional to the mean fluid velocity, but independent of the mean magnetic field, may occur in an originally homogeneous isotropic mhd turbulence if there are nonzero correlations of velocity and electric current fluctuations or, what is equivalent, of vorticity and magnetic field fluctuations. This goes beyond the Yoshizawa effect, which consists in the occurrence of mean electromotive forces proportional to the mean vorticity or to the angular velocity defining the Coriolis force in a rotating frame and depends on the cross-helicity defined by the velocity and magnetic field fluctuations. Contributions to the mean electromotive force due to inhomogeneity of the turbulence are also considered. Possible consequences of the above and related findings for the generation of magnetic fields in cosmic bodies are discussed.

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.

The role of magnetic helicity in astrophysical large-scale dynamos is reviewed and compared with cases where there is no energy supply and an initial magnetic field can only decay. In both cases magnetic energy tends to get redistributed to larger scales. Depending on the efficiency of magnetic helicity fluxes, the decay of a helical field can speed up. Likewise, the saturation of a helical dynamo can speed up through magnetic helicity fluxes. The astrophysical importance of these processes is reviewed in the context of the solar dynamo and an estimated upper limit for the magnetic helicity flux of 10^46 Mx^2/cycle is given.

Using two- and three-dimensional hydromagnetic simulations for a range of different flows, including laminar and turbulent ones, it is shown that solutions expressing the field in terms of Euler potentials (EP) are in general incorrect if the EP are evolved with an artificial diffusion term. In three dimensions, standard methods using the magnetic vector potential are found to permit dynamo action when the EP give decaying solutions. With an imposed field, the EP method yields excessive power at small scales. This effect is more exaggerated in the dynamic case, suggesting an unrealistically reduced feedback from the Lorentz force. The EP approach agrees with standard methods only at early times when magnetic diffusivity did not have time to act. It is demonstrated that the usage of EP with even a small artificial magnetic diffusivity does not converge to a proper solution of hydromagnetic turbulence. The source of this disagreement is not connected with magnetic helicity or the three-dimensionality of the magnetic field, but is simply due to the fact that the nonlinear representation of the magnetic field in terms of EP that depend on the same coordinates is incompatible with the linear diffusion operator in the induction equation.

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.

The existence of large-scale dynamos in rigidly rotating turbulent convection without shear is studied using three-dimensional numerical simulations of penetrative rotating compressible convection. We demonstrate that rotating convection in a Cartesian domain can drive a large-scale dynamo even in the absence of shear. The large-scale field contains a significant fraction of the total field in the saturated state. The simulation results are compared with one-dimensional mean-field dynamo models where turbulent transport coefficients, as determined using the test field method, are used. The reason for the absence of large-scale dynamo action in earlier studies is shown to be due to the rotation being too slow: whereas the alpha-effect can change sign, its magnitude stays approximately constant as a function of rotation, and the turbulent diffusivity decreases monotonically with increasing rotation. Only when rotation is rapid enough a large-scale dynamo can be excited. The one-dimensional mean-field model with dynamo coefficients from the test field results predicts reasonably well the dynamo excitation in the direct simulations. This result further validates the test field procedure and reinforces the interpretation that the observed dynamo is driven by a turbulent alpha-effect. This result demonstrates the existence of an alpha-effect and an alpha^2-dynamo with natural forcing.

Dynamo action owing to helically forced turbulence and large-scale shear is studied using direct numerical simulations. The resulting magnetic field displays propagating wave-like behavior. This behavior can be modelled in terms of an \alpha\Omega dynamo. In most cases super-equipartition fields are generated. By varying the fraction of helicity of the turbulence the regeneration of poloidal fields via the helicity effect (corresponding to the \alpha-effect) is regulated. The saturation level of the magnetic field in the numerical models is consistent with a linear dependence on the ratio of the fractional helicities of the small and large-scale fields, as predicted by a simple nonlinear mean-field model. As the magnetic Reynolds number (Rm) based on the wavenumber of the energy-carrying eddies is increased from 1 to 180, the cycle frequency of the large-scale field is found to decrease by a factor of about 6 in cases where the turbulence is fully helical. This is interpreted in terms of the turbulent magnetic diffusivity, which is found to be only weakly dependent on Rm.

(abridged) Aims: To study turbulent transport coefficients that describe the evolution of large-scale magnetic fields in turbulent convection. Methods: We use the test field method together with 3D numerical simulations of turbulent convection with shear and rotation to compute turbulent transport coefficients describing the evolution of large-scale magnetic fields in mean-field theory in the kinematic regime. 1D mean-field models are used with the derived turbulent transport coefficients to compare with direct simulations. Results: The alpha-effect increases monotonically as rotation increases. Turbulent diffusivity, eta_t, is proportional to the square of the turbulent vertical velocity. Whereas eta_t decreases approximately inversely proportional to the wavenumber of the field, the alpha-effect and turbulent pumping show a more complex behaviour. In the presence of shear and no rotation a small alpha-effect is induced which does not seem to show any consistent trend as a function of shear. If the shear is large enough, this small alpha is able to excite a dynamo in the mean-field model. The coefficient responsible for driving the shear-current effect shows several sign changes as a function of depth but is also able to contribute to dynamo action in the mean-field model. The growth rates in these cases are well below those in direct simulations suggesting that an incoherent alpha-shear dynamo may also act in them. If both rotation and shear are present, the alpha-effect is more pronounced. The combination of the shear-current and Omega x J-effects is also stronger than in the case of shear only, but subdominant to the alpha-shear dynamo. The results of direct simulations are consistent with mean-field models where all of these effects are taken into account without the need to invoke incoherent effects.

The generation of mean magnetic fields is studied for a simple non-helical flow where a net cross helicity of either sign can emerge. This flow, which is also known as the Archontis flow, is a generalization of the Arnold--Beltrami--Childress flow, but with the cosine terms omitted. The presence of cross helicity leads to a mean-field dynamo effect that is known as the Yoshizawa effect. Direct numerical simulations of such flows demonstrate the presence of magnetic fields on scales larger than the scale of the flow. Contrary to earlier expectations, the Yoshizawa effect is found to be proportional to the mean magnetic field and can therefore lead to its exponential instead of just linear amplification for magnetic Reynolds numbers that exceed a certain critical value. Unlike $\alpha$ effect dynamos, it is found that the Yoshizawa effect is not noticeably constrained by the presence of a conservation law. It is argued that this is due to the presence of a forcing term in the momentum equation which leads to a nonzero correlation with the magnetic field. Finally, the application to energy convergence in solar wind turbulence is discussed.

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.

The turbulent diffusion tensor describing the evolution of the mean concentration of a passive scalar is investigated for non-helically forced turbulence in the presence of rotation or a magnetic field. With rotation, the Coriolis force causes a sideways deflection of the flux of mean concentration. Within the magnetohydrodynamics approximation there is no analogous effect from the magnetic field because the effects on the flow do not depend on the sign of the field. rotation and magnetic fields tend to suppress turbulent transport, but this suppression is weaker in the direction along the magnetic field. Turbulent transport along the rotation axis is not strongly affected by rotation, except on shorter length scales, i.e., when the scale of the variation of the mean field becomes comparable with the scale of the energy-carrying eddies. These results are discussed in the context of anisotropic convective energy transport in the Sun.

Selected topics in solar dynamo theory are being highlighted. The possible relevance of the near-surface shear layer is discussed. The role of turbulent downward pumping is mentioned in connection with earlier concerns that a dynamo-generated magnetic field would be rapidly lost from the convection zone by magnetic buoyancy. It is argued that shear-mediated small-scale magnetic helicity fluxes are responsible for the success of some of the recent large-scale dynamo simulations. These fluxes help in disposing of excess small-scale magnetic helicity. This small-scale magnetic helicity, in turn, is generated in response to the production of an overall tilt in each Parker loop. Some preliminary calculations of this helicity flux are presented for a system with uniform shear. In the Sun the effects of magnetic helicity fluxes may be seen in coronal mass ejections shedding large amounts of magnetic helicity.

Using mean-field models with a dynamical quenching formalism we show that in finite domains magnetic helicity fluxes associated with small-scale magnetic fields are able to alleviate catastrophic quenching. We consider fluxes that result either from advection by a mean flow, the turbulent mixing down the gradient of mean small-scale magnetic helicity concentration, or the explicit removal which may be associated with the effects of coronal mass ejections in the Sun. In the absence of shear, all the small-scale magnetic helicity fluxes are found to be equally strong both for large-scale and small-scale fields. In the presence of shear there is also an additional magnetic helicity flux associated with the mean field, but this flux does not alleviate catastrophic quenching. Outside the dynamo-active region there are neither sources nor sinks of magnetic helicity, so in a steady state this flux must be constant. It is shown that unphysical behavior emerges if the small-scale magnetic helicity flux is forced to vanish within the computational domain.

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.

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.

The axino is a promising candidate for dark matter in the Universe. It is electrically and color neutral, very weakly interacting, and could be - as assumed in this study - the lightest supersymmetric particle, which is stable for unbroken R-parity. In supersymmetric extensions of the standard model, in which the strong CP problem is solved via the Peccei-Quinn mechanism, the axino arises naturally as the fermionic superpartner of the axion. We compute the thermal production rate of axinos in supersymmetric QCD. Using hard thermal loop resummation, we obtain a finite result in a gauge-invariant way, which takes into account Debye screening in the hot quark-gluon-squark-gluino plasma. The relic axino abundance from thermal scatterings after inflation is evaluated. We find that thermally produced axinos could provide the dominant part of cold dark matter, for example, for an axino mass of 100 keV and a reheating temperature of 10^6 GeV.

Recent analytical and computational advances in the theory of large-scale dynamos are reviewed. The importance of the magnetic helicity constraint is apparent even without invoking mean-field theory. The tau approximation yields expressions that show how the magnetic helicity gets incorporated into mean-field theory. The test-field method allows an accurate numerical determination of turbulent transport coefficients in linear and nonlinear regimes. Finally, some critical views on the solar dynamo are being offered and targets for future research are highlighted.

Simulations of stochastically forced shear-flow turbulence in a shearing-periodic domain are used to study the spontaneous generation of large-scale flow patterns in the direction perpendicular to the plane of the shear. Based on an analysis of the resulting large-scale velocity correlations it is argued that the mechanism behind this phenomenon could be the mean-vorticity dynamo effect pioneered by Elperin, Kleeorin, and Rogachevskii in 2003 (Phys. Rev. E 68, 016311). This effect is based on the anisotropy of the eddy viscosity tensor. One of its components may be able to replenish cross-stream mean flows by acting upon the streamwise component of the mean flow. Shear, in turn, closes the loop by acting upon the cross-stream mean flow to produce stronger streamwise mean flows. The diagonal component of the eddy viscosity is found to be of the order of the rms turbulent velocity divided by the wavenumber of the energy-carrying eddies.

We study the dependence of turbulent transport coefficients, such as the components of the $\alpha$ tensor ($\alpha_{ij}$) and the turbulent magnetic diffusivity tensor ($\eta_{ij}$), on shear and magnetic Reynolds number in the presence of helical forcing. We use three-dimensional direct numerical simulations with periodic boundary conditions and measure the turbulent transport coefficients using the kinematic test field method. In all cases the magnetic Prandtl number is taken as unity. We find that with increasing shear the diagonal components of $\alpha_{ij}$ quench, whereas those of $\eta_{ij}$ increase. The antisymmetric parts of both tensors increase with increasing shear. We also propose a simple expression for the turbulent pumping velocity (or $\gamma$ effect). This pumping velocity is proportional to the kinetic helicity of the turbulence and the vorticity of the mean flow. For negative helicity, i.e. for a positive trace of $\alpha_{ij}$, it points in the direction of the mean vorticity, i.e. perpendicular to the plane of the shear flow. Our simulations support this expression for low shear and magnetic Reynolds number. The transport coefficients depend on the wavenumber of the mean flow in a Lorentzian fashion, just as for non-shearing turbulence.

We present local numerical models of accretion disk turbulence driven by the magnetorotational instability with varying shear rate. The resulting turbulent stresses are compared with predictions of a closure model in which triple correlations are modelled in terms of quadratic correlations. This local model uses five nondimensional parameters to describe the properties of the flow. We attempt to determine these closure parameters for our simulations and find that the model does produce qualitatively correct behaviour. In addition, we present results concerning the shear rate dependency of the magnetic to kinetic energy ratio. We find both the turbulent stress ratio and the total stress to be strongly dependent on the shear rate.

Using direct simulations, weakly nonlinear theory and nonlinear mean-field theory, it is shown that the quenched velocity field of a saturated nonlinear dynamo can itself act as a kinematic dynamo. The flow is driven by a forcing function that would produce a Roberts flow in the absence of a magnetic field. This result confirms an analogous finding by F. Cattaneo & S. M. Tobias (arXiv:0809.1801) for the more complicated case of turbulent convection, suggesting that this may be a common property of nonlinear dynamos; see also the talk given also online at the Kavli Institute for Theoretical Physics (http://online.kitp.ucsb.edu/online/dynamo_c08/cattaneo). It is argued that this property can be used to test nonlinear mean-field dynamo theories.

The effect of a dynamo-generated mean magnetic field of Beltrami type on the mean electromotive force is studied. In the absence of the mean magnetic field the turbulence is assumed to be homogeneous and isotropic, but it becomes inhomogeneous and anisotropic with this field. Using the testfield method the dependence of the alpha and turbulent diffusivity tensors on the magnetic Reynolds number Rm is determined for magnetic fields that have reached approximate equipartition with the velocity field. The tensor components are characterized by a pseudoscalar alpha and a scalar turbulent magnetic diffusivity etat. Increasing Rm from 2 to 600 reduces etat by a factor ~5, suggesting that the quenching of etat is, in contrast to the 2-dimensional case, only weakly dependent on Rm. Over the same range of Rm, however, alpha is reduced by a factor ~14, which can qualitatively be explained by a corresponding increase of a magnetic contribution to the alpha effect with opposite sign. The level of fluctuations of alpha and etat is only 10% and 20% of the respective kinematic reference values.

(abridged) Aims: Three-dimensional numerical simulations of penetrative compressible convection with uniform horizontal shear are used to study dynamo action and the generation of large-scale magnetic fields. Methods: We consider cases where the magnetic Reynolds number is either marginal or moderately supercritical with respect to small-scale dynamo action in the absence of shear and rotation. The effects of magnetic helicity fluxes are studied by comparing results for the magnetic field with open and closed boundaries. Results: Without shear no large-scale dynamos are found even if the ingredients necessary for the alpha-effect (rotation and stratification) are present in the system. When uniform horizontal shear is added, a large-scale magnetic field develops, provided the boundaries are open. In this case the mean magnetic field contains a significant fraction of the total field. For those runs where the magnetic Reynolds number is between 60 and 250, small-scale dynamo is expected to be excited, but the field distribution is found to be similar to cases with smaller magnetic Reynolds number where the small-scale dynamo is not excited. In the case of closed (perfectly conducting) boundaries, magnetic helicity fluxes are suppressed and no large-scale fields are found. Similarly, poor large-scale field development is seen when vertical shear is used in combination with periodic boundary conditions in the horizontal directions. If, however, open (normal-field) boundary conditions are used in the x-direction, a large-scale field develops. These results support the notion that shear not only helps to generate the field, but it also plays a crucial role in driving magnetic helicity fluxes out of the system along the isocontours of shear, thereby allowing efficient dynamo action.

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.

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.

Aspects of turbulence in protostellar accretion discs are being reviewed. The emergence of dead zones due to poor ionization and alternatives to the magneto-rotational instability are discussed. The coupling between dust and gas in protostellar accretion discs is explained and turbulent drag is compared with laminar drag in the Stokes and Epstein regimes. Finally, the significance of magnetic field generation in turbulent discs is emphasized in connection with driving outflows and with star-disc coupling.

The role of shear in alleviating catastrophic quenching by shedding small-scale magnetic helicity through fluxes along contours of constant shear is discussed. The level of quenching of the dynamo effect depends on the quenched value of the turbulent magnetic diffusivity. Earlier estimates that might have suffered from the force-free degeneracy of Beltrami fields are now confirmed for shear flows where this degeneracy is lifted. For a dynamo that is saturated near equipartition field strength those estimates result in a 5-fold decrease of the magnetic diffusivity as the magnetic Reynolds number based on the wavenumber of the energy-carrying eddies is increased from 2 to 600. Finally, the role of shear in driving turbulence and large-scale fields by the magneto-rotational instability is emphasized. New simulations are presented and the 3pi/4 phase shift between poloidal and toroidal fields is confirmed. It is suggested that this phase shift might be a useful diagnostic tool in identifying mean-field dynamo action in simulations and to distinguish this from other scenarios invoking magnetic buoyancy as a means to explain migration away from the midplane.

Aims: We determine the components of the $\Lambda$-effect tensor that quantifies the contributions to the turbulent momentum transport even for uniform rotation. Methods: Three-dimensional numerical simulations are used to study turbulent transport in triply periodic cubes under the influence of rotation and anisotropic forcing. Comparison is made with analytical results obtained via the so-called minimal tau-approximation. Results: In the case where the turbulence intensity in the vertical direction dominates, the vertical stress is always negative. This situation is expected to occur in stellar convection zones. The horizontal component of the stress is weaker and exhibits a maximum at latitude 30 degrees - regardless of how rapid the rotation is. The minimal tau-approximation captures many of the qualitative features of the numerical results, provided the relaxation time tau is close to the turnover time, i.e. the Strouhal number is of order unity.

We analyze the mass distribution of cores formed in an isothermal, magnetized, turbulent, and self-gravitating nearly critical molecular cloud model. Cores are identified at two density threshold levels. Our main results are that the presence of self-gravity modifies the slopes of the core mass function (CMF) at the high mass end. At low thresholds, the slope is shallower than the one predicted by pure turbulent fragmentation. The shallowness of the slope is due to the effects of core coalescence and gas accretion. Most importantly, the slope of the CMF at the high mass end steepens when cores are selected at higher density thresholds, or alternatively, if the CMF is fitted with a log-normal function, the width of the lognormal distribution decreases with increasing threshold. This is due to the fact that gravity plays a more important role in denser structures selected at higher density threshold and leads to the conclusion that the role of gravity is essential in generating a CMF that bears more resemblance with the IMF when cores are selected with an increasing density threshold in the observations.

The turbulent magnetic diffusivity tensor is determined in the presence of rotation or shear. The question is addressed whether dynamo action from the shear-current effect can explain large-scale magnetic field generation found in simulations with shear. For this purpose a set of evolution equations for the response to imposed test fields is solved with turbulent and mean motions calculated from the momentum and continuity equations. The corresponding results for the electromotive force are used to calculate turbulent transport coefficients. The diagonal components of the turbulent magnetic diffusivity tensor are found to be very close together, but their values increase slightly with increasing shear and decrease with increasing rotation rate. In the presence of shear, the sign of the two off-diagonal components of the turbulent magnetic diffusion tensor is the same and opposite to the sign of the shear. This implies that dynamo action from the shear--current effect is impossible, except perhaps for high magnetic Reynolds numbers. However, even though there is no alpha effect on the average, the components of the alpha tensor display Gaussian fluctuations around zero. These fluctuations are strong enough to drive an incoherent alpha--shear dynamo. The incoherent shear--current effect, on the other hand, is found to be subdominant.

In an attempt to model the accretion onto a neutron star in low-mass X-ray binaries, we present two-dimensional hydrodynamical models of the gas flow in close vicinity of the stellar surface. First we consider a gas pressure dominated case, assuming that the star is non-rotating. For the stellar mass we take $M_{\rm star}=1.4 \times 10^{-2} \msun$ and for the gas temperature $T=5 \times 10^6$ K. Our results are qualitatively different in the case of a realistic neutron star mass and a realistic gas temperature of $T\simeq 10^8$ K, when the radiation pressure dominates. We show that to get the stationary solution in a latter case, the star most probably has to rotate with the considerable velocity.

We introduce on/off intermittency into a mean field dynamo model by imposing stochastic fluctuations in either the alpha effect or through the inclusion of a fluctuating electromotive force. Sufficiently strong small scale fluctuations with time scales of the order of 0.3-3 years can produce long term variations in the system on time scales of the order of hundreds of years. However, global suppression of magnetic activity in both hemispheres at once was not observed. The variation of the magnetic field does not resemble that of the sunspot number, but is more reminiscent of the 10 Be record. The interpretation of our results focuses attention on the connection between the level of magnetic activity and the sunspot number, an issue that must be elucidated if long term solar effects are to be well understood.

Aims: To determine alpha effect and turbulent magnetic diffusivity for mean magnetic fields with profiles of different length scale from simulations of isotropic turbulence, and to relate these results to nonlocal formulations in which alpha and the turbulent magnetic diffusivity correspond to integral kernels. Methods: A set of evolution equations for magnetic fields is solved which gives the response to imposed test fields, that is, mean magnetic fields with various wavenumbers. Both an imposed fully helical steady flow consisting of a pattern of screw-like motions (Roberts flow) and time-dependent statistically steady isotropic turbulence are considered. In the latter case the aforementioned evolution equations are solved simultaneously with the momentum and continuity equations. The corresponding results for the electromotive force are used to calculate alpha and magnetic diffusivity tensors. Results: For both the Roberts flow under the second--order correlation approximation and isotropic turbulence alpha and turbulent magnetic diffusivity are largest at large scales and these values diminish toward smaller scales. In both cases the alpha effect and turbulent diffusion kernels are well approximated by exponentials, corresponding to Lorentzian profiles in Fourier space. For isotropic turbulence the turbulent diffusion kernel is half as wide as the alpha effect kernel. For the Roberts flow beyond the second--order correlation approximation the turbulent diffusion kernel becomes negative at large scales.

A simple explicit example of a Roberts-type dynamo is given in which the alpha-effect of mean-field electrodynamics exists in spite of point-wise vanishing kinetic helicity of the fluid flow. In this way it is shown that alpha-effect dynamos do not necessarily require non-zero kinetic helicity. A mean-field theory of Roberts-type dynamos is established within the framework of the second-order correlation approximation. In addition numerical solutions of the original dynamo equations are given, that are independent of any approximation of that kind. Both theory and numerical results demonstrate the possibility of dynamo action in the absence of kinetic helicity.

Using numerical simulations at moderate magnetic Reynolds numbers up to 220 it is shown that in the kinematic regime, isotropic helical turbulence leads to an alpha effect and a turbulent diffusivity whose values are independent of the magnetic Reynolds number, $\Rm$, provided $\Rm$ exceeds unity. These turbulent coefficients are also consistent with expectations from the first order smoothing approximation. For small values of $\Rm$, alpha and turbulent diffusivity are proportional to $\Rm$. Over finite time intervals meaningful values of alpha and turbulent diffusivity can be obtained even when there is small-scale dynamo action that produces strong magnetic fluctuations. This suggests that small-scale dynamo-generated fields do not make a correlated contribution to the mean electromotive force.

The motivation for considering distributed large scale dynamos in the solar context is reviewed in connection with the magnetic helicity constraint. Preliminary accounts of 3-dimensional direct numerical simulations (in spherical shell segments) and simulations of 2-dimensional mean field models (in spherical shells) are presented. Interesting similarities as well as some differences are noted.

Homogeneous anisotropic turbulence simulations are used to determine off-diagonal components of the Reynolds stress tensor and its parameterization in terms of turbulent viscosity and Lambda-effect. The turbulence is forced in an anisotropic fashion by enhancing the strength of the forcing in the vertical direction. The Coriolis force is included with a rotation axis inclined relative to the vertical direction. The system studied here is significantly simpler than that of turbulent stratified convection which has often been used to study Reynolds stresses. Certain puzzling features of the results for convection, such as sign changes or highly concentrated latitude distributions, are not present in the simpler system considered here.

We evaluate the gravitino production rate in supersymmetric QCD at high temperature to leading order in the gauge coupling. The result, which is obtained by using the resummed gluon propagator, depends logarithmically on the gluon plasma mass. As a byproduct, a new result for the axion production rate in a QED plasma is obtained. The implicatons for the cosmological dark matter problem are briefly discussed, in particular the intriguing possibility that gravitinos are the dominant part of cold dark matter.

Magnetic helicity effects are discussed in laboratory and astrophysical settings. First, dynamo action in Taylor-Green flows is discussed for different boundary conditions. However, because of the lack of scale separation with respect to the container, no large scale field is being produced and there is no resistively slow saturation phase as otherwise expected. Second, the build-up of a large scale field is demonstrated in a simulation where a localized magnetic eddy produces field on a larger scale if the eddy possesses a swirl. Such a set-up might be realizable experimentally through coils. Finally, new emerging issues regarding the connection between magnetic helicity and the solar dynamo are discussed. It is demonstrated that dynamos with a non-local (Babcock-Leighton type) alpha effect can also be catastrophically quenched, unless there are magnetic helicity fluxes.

The validity of a closure called the minimal tau approximation (MTA), is tested in the context of dynamo theory, wherein triple correlations are assumed to provide relaxation of the turbulent electromotive force. Under MTA, the alpha effect in mean field dynamo theory becomes proportional to a relaxation time scale multiplied by the difference between kinetic and current helicities. It is shown that the value of the relaxation time is positive and, in units of the turnover time at the forcing wavenumber, it is of the order of unity. It is quenched by the magnetic field -- roughly independently of the magnetic Reynolds number. However, this independence becomes uncertain at large magnetic Reynolds number. Kinetic and current helicities are shown to be dominated by large scale properties of the flow.

Some basic properties of the solar convection zone are considered and the use of helioseismology as an observational tool to determine its depth and internal angular velocity is discussed. Aspects of solar magnetism are described and explained in the framework of dynamo theory. The main focus is on mean field theories for the Sun's magnetic field and its differential rotation.

It is argued that accretion discs in young stellar objects may have hot coronae that are heated by magnetic reconnection. This is a consequence of the magneto-rotational instability driving turbulence in the disc. Magnetic reconnection away from the midplane leads to heating of the corona which, in turn, contributes to driving disc winds.

The backreaction of the Lorentz force on the alpha-effect is studied in the limit of small magnetic and fluid Reynolds numbers, using the first order smoothing approximation (FOSA) to solve both the induction and momentum equations. Both steady and time dependent forcings are considered. In the low Reynolds number limit, the velocity and magnetic fields can be expressed explicitly in terms of the forcing function. The nonlinear alpha-effect is then shown to be expressible in several equivalent forms in agreement with formalisms that are used in various closure schemes. On the one hand, one can express alpha completely in terms of the helical properties of the velocity field as in traditional FOSA, or, alternatively, as the sum of two terms, a so-called kinetic alpha-effect and an oppositely signed term proportional to the helical part of the small scale magnetic field. These results hold for both steady and time dependent forcing at arbitary strength of the mean field (abridged).

The thermal instability with a piecewise power law cooling function is investigated using one- and three-dimensional simulations with periodic and shearing-periodic boundary conditions in the presence of constant thermal diffusion and kinematic viscosity coefficients. Consistent with earlier findings, the flow behavior depends on the average density, <rho>. When $\bra{\rho}$ is in the range 1-5 x 10^{-24} g cm^{-3} the system is unstable and segregates into cool and warm phases with temperatures of roughly 100 and 10^4 K, respectively. However, in all cases the resulting average pressure <p> is independent of <rho> and just a little above the minimum value. For a constant heating rate of 0.015 erg g^{-1} s^{-1}, the mean pressure is around 24 x 10^{-14} dyn (corresponding to p/k_B ~ K cm^{-3}). Cool patches tend to coalesce into bigger ones. In all cases investigated there is no sustained turbulence, which is in agreement with earlier results. Simulations in which turbulence is driven by a body force show that when rms velocities of between 10 and 30 km/s are obtained, the resulting dissipation rates rates are comparable to the thermal energy input rate. The resulting mean pressures are then about 30 x 10^{-14} dyn, corresponding to p/k_B ~ 2170 K cm^{-3}. This is comparable to the value expected for the Galaxy. Differential rotation tends to make the flow two-dimensional, that is, uniform in the streamwise direction, but this does not lead to instability.

Large scale dynamo-generated fields are a combination of interlocked poloidal and toroidal fields. Such fields possess magnetic helicity that needs to be regenerated and destroyed during each cycle. A number of numerical experiments now suggests that stars may do this by shedding magnetic helicity. In addition to plain bulk motions, a favorite mechanism involves magnetic helicity flux along lines of constant rotation. We also know that the sun does shed the required amount of magnetic helicity mostly in the form of coronal mass ejections. Solar-like stars without cycles do not face such strong constraints imposed by magnetic helicity evolution and may not display coronal activity to that same extent. I discuss the evidence leading to this line of argument. In particular, I discuss simulations showing the generation of strong mean toroidal fields provided the outer boundary condition is left open so as to allow magnetic helicity to escape. Control experiments with closed boundaries do not produce strong mean fields.

The outer surface layers of the sun show a clear deceleration at low latitudes. This is generally thought to be the result of a strong dominance of vertical turbulent motions associated with strong downdrafts. This strong negative radial shear should not only contribute to amplifying the toroidal field locally and to expelling magnetic helicity, but it may also be responsible for producing a strong prograde pattern speed in the supergranulation layer. Using simulations of rotating stratified convection in cartesian boxes located at low latitudes around the equator it is shown that in the surface layers patterns move in the prograde direction on top of a retrograde mean background flow. These patterns may also be associated with magnetic tracers and even sunspot proper motions that are known to be prograde relative to the much slower surface plasma.

The macroscopic behaviour of cosmic rays in turbulent magnetic fields is discussed. An implementation of anisotropic diffusion of cosmic rays with respect to the magnetic field in a non-conservative, high-order, finite-difference magnetohydrodynamic code is discussed. It is shown that the standard implementation fails near singular X-points of the magnetic field, which are common if the field is random. A modification to the diffusion model for cosmic rays is described and the resulting telegraph equation (implemented by solving a dynamic equation for the diffusive flux of cosmic rays) is used; it is argued that this modification may better describe the physics of cosmic ray diffusion. The present model reproduces several processes important for the propagation and local confinement of cosmic rays, including spreading perpendicular to the local large-scale magnetic field, controlled by the random-to-total magnetic field ratio, and the balance between cosmic ray pressure and magnetic tension. Cosmic ray diffusion is discussed in the context of a random magnetic field produced by turbulent dynamo action. It is argued that energy equipartition between cosmic rays and other constituents of the interstellar medium do not necessarily imply that cosmic rays play a significant role in the balance of forces.

In an attempt to determine the outer scale of turbulence driven by localized sources, such as supernova explosions in the interstellar medium, we consider a forcing function given by the gradient of gaussian profiles localized at random positions. Different coherence times of the forcing function are considered. In order to isolate the effects specific to the nature of the forcing function we consider the case of a polytropic equation of state and restrict ourselves to forcing amplitudes such that the flow remains subsonic. When the coherence time is short, the outer scale agrees with the half-width of the gaussian. Longer coherence times can cause extra power at large scales, but this would not yield power law behavior at scales larger than that of the expansion waves. At scales smaller than the scale of the expansion waves the spectrum is close to power law with a spectral exponent of -2. The resulting flow is virtually free of vorticity. Viscous driving of vorticity turns out to be weak and self-amplification through the nonlinear term is found to be insignificant. No evidence for small scale dynamo action is found in cases where the magnetic induction equation is solved simultaneously with the other equations.

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.

Direct and large eddy simulations of hydrodynamic and hydromagnetic turbulence have been performed in an attempt to isolate artifacts from real and possibly asymptotic features in the energy spectra. It is shown that in a hydrodynamic turbulence simulation with a Smagorinsky subgrid scale model using 512^3 meshpoints two important features of the 4096^3 simulation on the Earth simulator (Kaneda et al. 2003, Phys. Fluids 15, L21) are reproduced: a k^{-0.1} correction to the inertial range with a k^{-5/3} Kolmogorov slope and the form of the bottleneck just before the dissipative subrange. Furthermore, it is shown that, while a Smagorinsky-type model for the induction equation causes an artificial and unacceptable reduction in the dynamo efficiency, hyper-resistivity yields good agreement with direct simulations. In the large-scale part of the inertial range, an excess of the spectral magnetic energy over the spectral kinetic energy is confirmed. However, a trend towards spectral equipartition at smaller scales in the inertial range can be identified. With magnetic fields, no explicit bottleneck effect is seen.

Aim: To study turbulence in the Orion Molecular Cloud (OMC1) by comparing observed and simulated characteristics of the gas motions. Method: Using a dataset of vibrationally excited H2 emission in OMC1 containing radial velocity and brightness which covers scales from 70AU to 30000AU, we present the transversal structure functions and the scaling of the structure functions with their order. These are compared with the predictions of two-dimensional projections of simulations of supersonic hydrodynamic turbulence. Results: The structure functions of OMC1 are not well represented by power laws, but show clear deviations below 2000AU. However, using the technique of extended self-similarity, power laws are recovered at scales down to 160AU. The scaling of the higher order structure functions with order deviates from the standard scaling for supersonic turbulence. This is explained as a selection effect of preferentially observing the shocked part of the gas and the scaling can be reproduced using line-of-sight integrated velocity data from subsets of supersonic turbulence simulations. These subsets select regions of strong flow convergence and high density associated with shock structure. Deviations of the structure functions in OMC1 from power laws cannot however be reproduced in simulations and remains an outstanding issue.

We review advances in the theoretical and computational studies of disk winds, jets and outflows including: the connection between accretion and jets, the launch of jets from magnetized disks, the coupled evolution of jets and disks, the interaction of magnetized young stellar objects with their surrounding disks and the relevance to outflows, and finally, the link between jet formation and gravitational collapse. We also address the predictions that the theory makes about jet kinematics, collimation, and rotation, that have recently been confirmed by high spatial and spectral resolution observations. Disk winds have a universal character that may account for jets and outflows during the formation of massive stars as well as brown dwarfs.

Some common properties of helical magnetic fields in decaying and driven turbulence are discussed. These include mainly the inverse cascade that produces fields on progressively larger scales. Magnetic helicity also restricts the evolution of the large scale field: the field decays less rapidly than a non-helical field, but it also saturates more slowly, i.e. on a resistive time scale if there are no magnetic helicity fluxes. The former effect is utilized in primordial field scenarios, while the latter is important for successfully explaining astrophysical dynamos that saturate faster than resistively. Dynamo action is argued to be important not only in the galactic dynamo, but also in accretion discs in active galactic nuclei and around protostars, both of which contribute to producing a strong enough seed magnetic field. Although primordial magnetic fields may be too weak to compete with these astrophysical mechanisms, such fields could perhaps still be important in producing polarization effects in the cosmic background radiation.

Magnetohydrodynamic simulations of fully convective, rotating spheres with volume heating near the center and cooling at the surface are presented. The dynamo-generated magnetic field saturates at equipartition field strength near the surface. In the interior, the field is dominated by small-scale structures, but outside the sphere by the global scale. Azimuthal averages of the field reveal a large-scale field of smaller amplitude also inside the star. The internal angular velocity shows some tendency to be constant along cylinders and is ``anti-solar'' (fastest at the poles and slowest at the equator).

Aims: To quantify the transient growth of nonaxisymmetric perturbations in unstratified magnetized and stratified non-magnetized rotating linear shear flows in the shearing sheet approximation of accretion disc flows. Method: The Rayleigh quotient in modal approaches for the linearized equations (with time-dependent wavenumber) and the amplitudes from direct shearing sheet simulations using a finite difference code are compared. Results: Both approaches agree in their predicted growth behavior. The magneto-rotational instability for axisymmetric and non-axisymmetric perturbations is shown to have the same dependence of the (instantaneous) growth rate on the wavenumber along the magnetic field, but in the nonaxisymmetric case the growth is only transient. However, a meaningful dependence of the Rayleigh quotient on the radial wavenumber is obtained. While in the magnetized case the total amplification factor can be several orders of magnitude, it is only of order ten or less in the nonmagnetic case. Stratification is shown to have a stabilizing effect. In the present case of shearing-periodic boundaries the (local) strato-rotational instability seems to be absent.

Nonlinear behaviour of galactic dynamos is studied, allowing for magnetic helicity removal by the galactic fountain flow. A suitable advection speed is estimated, and a one-dimensional mean-field dynamo model with dynamic alpha-effect is explored. It is shown that the galactic fountain flow is efficient in removing magnetic helicity from galactic discs. This alleviates the constraint on the galactic mean-field dynamo resulting from magnetic helicity conservation and thereby allows the mean magnetic field to saturate at a strength comparable to equipartition with the turbulent kinetic energy.

A two-dimensional mean field dynamo model is solved where magnetic helicity conservation is fully included. The model has a negative radial velocity gradient giving rise to equatorward migration of magnetic activity patterns. In addition the model develops longitudinal variability with activity patches travelling in longitude. These patches may be associated with active longitudes.

Arguments are presented in favor of the idea that the solar dynamo may operate not just at the bottom of the convection zone, i.e. in the tachocline, but it may operate in a more distributed fashion in the entire convection zone. The near-surface shear layer is likely to play an important role in this scenario.

The location of the solar dynamo is discussed in the context of new insights into the theory of nonlinear turbulent dynamos. It is argued that, from a dynamo-theoretic point of view, the bottom of the convection zone is not a likely location and that the solar dynamo may be distributed over the convection zone. The near surface shear layer produces not only east-west field alignment, but it also helps the dynamo disposing of its excess small scale magnetic helicity.

Aims: To investigate the stability and angular momentum transport by the strato-rotational instability in the nonlinear regime. Methods: The hydrodynamic compressible equations are solved in a cartesian box in which the outer cylinder is embedded. Gravity along the rotation axis leads to density stratification. No-slip boundary conditions are used in the radial direction, while free-slip conditions are used on the two ends of the cylinders. Results: The strato-rotational instability is confirmed and the Reynolds stress is shown to transport angular momentum away from the axis. However, the growth rate decreases with increasing Reynolds number. This, as well as the presence of boundaries renders this instability less relevant for astrophysical applications.

An efficient algorithm for calculating radiative transfer on massively parallel computers using domain decomposition is presented. The integral formulation of the transfer equation is used to divide the problem into a local but compute-intensive part for calculating the intensity and optical depth integrals, and a nonlocal part for communicating the intensity between adjacent processors. The waiting time of idle processors during the nonlocal communication part does not have a severe impact on the scaling. The wall clock time thus scales nearly linearly with the inverse number of processors.

Star-disc coupling is considered in numerical models where the stellar field is not an imposed perfect dipole, but instead a more irregular self-adjusting dynamo-generated field. Using axisymmetric simulations of the hydromagnetic mean-field equations, it is shown that the resulting stellar field configuration is more complex, but significantly better suited for driving a stellar wind. In agreement with recent findings by a number of people, star-disc coupling is less efficient in braking the star than previously thought. Moreover, stellar wind braking becomes equally important. In contrast to a perfect stellar dipole field, dynamo-generated stellar fields favor field-aligned accretion with considerably higher velocity at low latitudes, where the field is weaker and originating in the disc. Accretion is no longer nearly periodic (as it is in the case of a stellar dipole), but it is more irregular and episodic.

Accretion disc turbulence is investigated in the framework of the shearing box approximation. The turbulence is either driven by the magneto-rotational instability or, in the non-magnetic case, by an explicit and artificial forcing term in the momentum equation. Unlike the magnetic case, where most of the dissipation occurs in the disc corona, in the forced hydrodynamic case most of the dissipation occurs near the midplane. In the hydrodynamic case evidence is presented for the stochastic excitation of epicycles. When the vertical and radial epicyclic frequencies are different (modeling the properties around rotating black holes), the beat frequency between these two frequencies appear to show up as a peak in the temporal power spectrum in some cases. Finally, the full turbulent resistivity tensor is determined and it is found that, if the turbulence is driven by a forcing term, the signs of its off-diagonal components are such that this effect would not be capable of dynamo action by the shear--current effect.

The validity of a closure called the minimal tau approximation (MTA), is tested in the context of dynamo theory, wherein triple correlations are assumed to provide relaxation of the turbulent electromotive force. Under MTA, the alpha effect in mean field dynamo theory becomes proportional to a relaxation time scale multiplied by the difference between kinetic and current helicities. It is shown that the value of the relaxation time is positive and, in units of the turnover time at the forcing wavenumber, it is of the order of unity. It is quenched by the magnetic field -- roughly independently of the magnetic Reynolds number. However, this independence becomes uncertain at large magnetic Reynolds number. Kinetic and current helicities are shown to be dominated by large scale properties of the flow.

Numerical simulations of turbulent stratified convection are used to study models with approximately the same convective flux, but different radiative fluxes. As the radiative flux is decreased, for constant convective flux: the entropy jump at the top of the convection zone becomes steeper, the temperature fluctuations increase and the velocity fluctuations decrease in magnitude, and the distance that low entropy fluid from the surface can penetrate increases. Velocity and temperature fluctuations follow mixing length scaling laws.

Several one and two dimensional mean field models are analyzed where the effects of current helicity fluxes and boundaries are included within the framework of the dynamical quenching model. In contrast to the case with periodic boundary conditions, the final saturation energy of the mean field decreases inversely proportional to the magnetic Reynolds number. If a nondimensional scaling factor in the current helicity flux exceeds a certain critical value, the dynamo can operate even without kinetic helicity, i.e. it is based only on shear and current helicity fluxes, as first suggested by Vishniac & Cho (2001, ApJ 550, 752). Only above this threshold is the current helicity flux also able to alleviate catastrophic quenching. The fact that certain turbulence simulations have now shown apparently non-resistively limited mean field saturation amplitudes may be suggestive of the current helicity flux having exceeded this critical value. Even below this critical value the field still reaches appreciable strength at the end of the kinematic phase, which is in qualitative agreement with dynamos in periodic domains. However, for large magnetic Reynolds numbers the field undergoes subsequent variations on a resistive time scale when, for long periods, the field can be extremely weak.

We study the evolution of growth and decay laws for the magnetic field coherence length xi, energy E_M and magnetic helicity H in freely decaying 3D MHD turbulence. We show that with certain assumptions, self-similarity of the magnetic power spectrum alone implies that xi~t^{1/2}. This in turn implies that magnetic helicity decays as H~t^{-2s}, where s=(xi_{diff}/xi_H)^2, in terms of xi_{diff}, the diffusion length scale, and xi_H, a length scale defined from the helicity power spectrum. The relative magnetic helicity remains constant, implying that the magnetic energy decays as E_M~t^{-1/2-2s}. The parameter s is inversely proportional to the magnetic Reynolds number Re_M, which is constant in the self-similar regime.

The current understanding of astrophysical magnetic fields is reviewed, focusing on their generation and maintenance by turbulence. In the astrophysical context this generation is usually explained by a self-excited dynamo, which involves flows that can amplify a weak 'seed' magnetic field exponentially fast. Particular emphasis is placed on the nonlinear saturation of the dynamo. Analytic and numerical results are discussed both for small scale dynamos, which are completely isotropic, and for large scale dynamos, where some form of parity breaking is crucial. Central to the discussion of large scale dynamos is the so-called alpha effect which explains the generation of a mean field if the turbulence lacks mirror symmetry, i.e. if the flow has kinetic helicity. Large scale dynamos produce small scale helical fields as a waste product that quench the large scale dynamo and hence the alpha effect. With this in mind, the microscopic theory of the alpha effect is revisited in full detail and recent results for the loss of helical magnetic fields are reviewed.

The axino and the gravitino are well-motivated candidates for the lightest supersymmetric particle (LSP) and also for cold dark matter in the Universe. Assuming that a charged slepton is the next-to-lightest supersymmetric particle (NLSP), we show how the NLSP decays can be used to probe the axino LSP scenario in hadronic axion models as well as the gravitino LSP scenario at the Large Hadron Collider and the International Linear Collider. We show how one can identify experimentally the scenario realized in nature. In the case of the axino LSP, the NLSP decays will allow one to estimate the value of the axino mass and the Peccei-Quinn scale.

We study numerically the dependence of the critical magnetic Reynolds number Rmc for the turbulent small-scale dynamo on the hydrodynamic Reynolds number Re. The turbulence is statistically homogeneous, isotropic, and mirror--symmetric. We are interested in the regime of low magnetic Prandtl number Pm=Rm/Re<1, which is relevant for stellar convective zones, protostellar disks, and laboratory liquid-metal experiments. The two asymptotic possibilities are Rmc->const as Re->infinity (a small-scale dynamo exists at low Pm) or Rmc/Re=Pmc->const as Re->infinity (no small-scale dynamo exists at low Pm). Results obtained in two independent sets of simulations of MHD turbulence using grid and spectral codes are brought together and found to be in quantitative agreement. We find that at currently accessible resolutions, Rmc grows with Re with no sign of approaching a constant limit. We reach the maximum values of Rmc~500 for Re~3000. By comparing simulations with Laplacian viscosity, fourth-, sixth-, and eighth-order hyperviscosity and Smagorinsky large-eddy viscosity, we find that Rmc is not sensitive to the particular form of the viscous cutoff. This work represents a significant extension of the studies previously published by Schekochihin et al. 2004, PRL 92, 054502 and Haugen et al. 2004, PRE, 70, 016308 and the first detailed scan of the numerically accessible part of the stability curve Rmc(Re).

The excitation of internal gravity waves by penetrative convective plumes is investigated using 2-D direct simulations of compressible convection. The wave generation is quantitatively studied from the linear response of the radiative zone to the plumes penetration, using projections onto the g-modes solutions of the associated linear eigenvalue problem for the perturbations. This allows an accurate determination of both the spectrum and amplitudes of the stochastically excited modes. Using time-frequency diagrams of the mode amplitudes, we then show that the lifetime of a mode is around twice its period and that during times of significant excitation up to 40% of the total kinetic energy may be contained into g-modes.

Arguments for and against the widely accepted picture of a solar dynamo being seated in the tachocline are reviewed and alternative ideas concerning dynamos operating in the bulk of the convection zone, or perhaps even in the near-surface shear layer, are discussed. Based on the angular velocities of magnetic tracers it is argued that the observations are compatible with a distributed dynamo that may be strongly shaped by the near-surface shear layer. Direct simulations of dynamo action in a slab with turbulence and shear are presented to discuss filling factor and tilt angles of bipolar regions in such a model.

Small scale turbulent dynamo action in compressible transonic turbulence is discussed. It is shown that the critical value of the magnetic Reynolds number displays a bimodal behavior and changes from a typical value of 35 for small Mach numbers to about 80 for larger Mach numbers. The transition between the two regimes is relatively sharp. The direct simulations are then compared with simulations where shocks are captured using a shock viscosity that becomes large at locations where there are shocks. In the presence of shear it is shown that large scale dynamo action is possible.

Two different computational approaches to the magnetorotational instability (MRI) are pursued: the shearing box approach which is suited for local simulations and the embedding box approach whereby a Taylor Couette flow is embedded in a box so that numerical problems with the coordinate singularity are avoided. New shearing box simulations are presented and differences between regular and hyperviscosity are discussed. Preliminary simulations of spherical nonlinear Taylor Couette flow in an embedding box are presented and the effects of an axial field on the background flow are studied.

Magnetic helicity is nearly conserved and its evolution equation provides a dynamical feedback on the alpha effect that is distinct from the conventional algebraic alpha quenching. The seriousness of this dynamical alpha quenching is particularly evident in the case of closed or periodic boxes. The explicit connection with catastrophic alpha quenching is reviewed and the alleviating effects of magnetic and current helicity fluxes are discussed.

Three closely related stumbling blocks of solar mean field dynamo theory are discussed: how dominant are the small scale fields, how is the alpha effect quenched, and whether magnetic and current helicity fluxes alleviate the quenching? It is shown that even at the largest currently available resolution there is no clear evidence of power law scaling of the magnetic and kinetic energy spectra in turbulence. However, using subgrid scale modeling, some indications of asymptotic equipartition can be found. The frequently used first order smoothing approach to calculate the alpha effect and other transport coefficients is contrasted with the superior minimal tau approximation. The possibility of catastrophic alpha quenching is discussed as a result of magnetic helicity conservation. Magnetic and current helicity fluxes are shown to alleviate catastrophic quenching in the presence of shear. Evidence for strong large scale dynamo action, even in the absence of helicity in the forcing, is presented.

Numerical turbulence with hyperviscosity is studied and compared with direct simulations using ordinary viscosity and data from wind tunnel experiments. It is shown that the inertial range scaling is similar in all three cases. Furthermore, the bottleneck effect is approximately equally broad (about one order of magnitude) in these cases and only its height is increased in the hyperviscous case--presumably as a consequence of the steeper decent of the spectrum in the hyperviscous subrange. The mean normalized dissipation rate is found to be in agreement with both wind tunnel experiments and direct simulations. The structure function exponents agree with the She-Leveque model. Decaying turbulence with hyperviscosity still gives the usual t^{-1.25} decay law for the kinetic energy, and also the bottleneck effect is still present and about equally strong.

Non-helical hydromagnetic turbulence with an externally imposed magnetic field is investigated using direct numerical simulations. It is shown that the imposed magnetic field lowers the spectral magnetic energy in the inertial range. This is explained by a suppression of the small scale dynamo. At large scales, however, the spectral magnetic energy increases with increasing imposed field strength for moderately strong fields, and decreases only slightly for even stronger fields. The presence of Alfven waves is explicitly confirmed by monitoring the evolution of magnetic field and velocity at one point. The frequency omega agrees with vA k1, where vA is the Alfven speed and k1 is the smallest wavenumber in the box.

The effect of compressibility on the onset of nonhelical turbulent dynamo action is investigated using both direct simulations as well as simulations with shock-capturing viscosities, keeping however the regular magnetic diffusivity. It is found that the critical magnetic Reynolds number increases from about 35 in the subsonic regime to about 70 in the supersonic regime. In the high resolution direct simulations the shock structures are much sharper than in the low resolution shock-capturing simulations. Nevertheless, the magnetic field looks roughly similar in both cases and does not show shock structures. Similarly, the onset of dynamo action is not significantly affected by the shock-capturing viscosity.

Large scale dynamos produce small scale current helicity as a waste product that quenches the large scale dynamo process (alpha effect). This quenching can be catastrophic (i.e.intensify with magnetic Reynolds number) unless one has fluxes of small scale magnetic (or current) helicity out of the system. We derive the form of helicity fluxes in turbulent dynamos, taking also into account the nonlinear effects of Lorentz forces due to fluctuating fields. We confirm the form of an earlier derived magnetic helicity flux term, and also show that it is not renormalized by the small scale magnetic field, just like turbulent diffusion. Additional nonlinear fluxes are identified, which are driven by the anisotropic and antisymmetric parts of the magnetic correlations. These could provide further ways for turbulent dynamos to transport out small scale magnetic helicity, so as to avoid catastrophic quenching.

The excitation of internal gravity waves by an entropy bubble oscillating in an isothermal atmosphere is investigated using direct two-dimensional numerical simulations. The oscillation field is measured by a projection of the simulated velocity field onto the anelastic solutions of the linear eigenvalue problem for the perturbations. This facilitates a quantitative study of both the spectrum and the amplitudes of excited g-modes.

Nonhelical hydromagnetic forced turbulence is investigated using large scale simulations on up to 256 processors and $1024^3$ meshpoints. The magnetic Prandtl number is varied between 1/8 and 30, although in most cases it is unity. When the magnetic Reynolds number is based on the inverse forcing wavenumber, the critical value for dynamo action is shown to be around 35 for magnetic Prandtl number of unity. For small magnetic Prandtl numbers we find the critical magnetic Reynolds number to increase with decreasing magnetic Prandtl number. The Kazantsev $k^{3/2}$ spectrum for magnetic energy is confirmed for the kinematic regime, i.e. when nonlinear effects are still unimportant and when the magnetic Prandtl number is unity. In the nonlinear regime, the energy budget converges for large Reynolds numbers (around 1000) such that for our parameters about 70% is in kinetic energy and about 30% is in magnetic energy. The energy dissipation rates are converged to 30% viscous dissipation and 70% resistive dissipation. Second order structure functions of the Elsasser variables give evidence for a $k^{-5/3}$ spectrum. Nevertheless, the three-dimensional spectrum is close to $k^{-3/2}$, but we argue that this is due to the bottleneck effect. The bottleneck effect is shown to be equally strong both for magnetic and nonmagnetic turbulence, but it is far weaker in one-dimensional spectra that are normally studied in laboratory turbulence. Structure function exponents for other orders are well described by the She-Leveque formula, but the velocity field is significantly less intermittent and the magnetic field is more intermittent than the Elsasser variables.

The standard theory of the solar cycle in terms of an alpha-Omega dynamo hinges on a proper understanding of the nonlinear alpha effect. Boundary conditions play a surprisingly important role in determining the magnitude of alpha. For closed boundaries, the total magnetic helicity is conserved, and since the alpha effect produces magnetic helicity of one sign in the large scale field, it must simultaneously produce magnetic helicity of the opposite sign. It is this secondary magnetic helicity that suppresses the dynamo in a potentially catastrophic fashion. Open boundaries allow magnetic helicity to be lost. Simulations are presented that allow an estimate of alpha in the presence of open or closed boundaries, either with or without solar-like differential rotation. In all cases the sign of the magnetic helicity agrees with that observed at the solar surface (negative in the north, positive in the south), where significant amounts of magnetic helicity can be ejected via coronal mass ejections. It is shown that open boundaries tend to alleviate catastrophic alpha quenching. The importance of looking at current helicity instead of magnetic helicity is emphasized and the conceptual advantages are discussed.

A new simulation set-up is proposed for studying mean field dynamo action. The model combines the computational advantages of local cartesian geometry with the ability to include a shear profile that resembles the sun's differential rotation at low latitudes. It is shown that in a two-dimensional mean field model this geometry produces cyclic solutions with dynamo waves traveling away from the equator -- as expected for a positive alpha effect in the northern hemisphere. In three dimensions with turbulence driven by a helical forcing function, an alpha effect is self-consistently generated in the presence of a finite imposed toroidal magnetic field. The results suggest that, due to a finite flux of current helicity out of the domain, alpha quenching appears to be non-catastrophic -- at least for intermediate values of the magnetic Reynolds number. For larger values of the magnetic Reynolds number, however, there is evidence for a reversal of the trend and that $\alpha$ may decrease with increasing magnetic Reynolds number. Control experiments with closed boundaries confirm that in the absence of a current helicity flux, but with shear as before, alpha quenching is always catastrophic and alpha decreases inversely proportional to the magnetic Reynolds number. For solar parameters, our results suggest a current helicity flux of about 0.001 G^2/s. This corresponds to a magnetic helicity flux, integrated over the northern hemisphere and over the 11 year solar cycle, of about 10^{46}Mx^2.

We compute the thermal axino production rate in supersymmetric QCD to leading order in the gauge coupling. Using hard thermal loop resummation and the Braaten-Yuan prescription, we obtain a finite result in a gauge-invariant way, which takes into account Debye screening in the hot quark-gluon-squark-gluino plasma. The relic axino density from thermal reactions in the early Universe is evaluated assuming the axino is the lightest supersymmetric particle and stable due to R-parity conservation. From the comparison with the WMAP results, we find that axinos could provide the dominant part of cold dark matter, for example, for an axino mass of 100 keV and a reheating temperature of 10^6 GeV.

Reaction-diffusion equations based on a polymerization model are solved to simulate the spreading of hypothetic left and right handed life forms on the Earth's surface. The equations exhibit front-like behavior as is familiar from the theory of the spreading of epidemics. It is shown that the relevant time scale for achieving global homochirality is not, however, the time scale of front propagation, but the much longer global diffusion time. The process can be sped up by turbulence and large scale flows. It is speculated that, if the deep layers of the early ocean were sufficiently quiescent, there may have been the possibility of competing early life forms with opposite handedness.

The axino is the fermionic superpartner of the axion. Assuming the axino is the lightest supersymmetric particle and stable due to R-parity conservation, we compute the relic axino density from thermal reactions in the early Universe. From the comparison with the WMAP results, we find that thermally produced axinos could provide the dominant part of cold dark matter, for example, for an axino mass of 100 keV and a reheating temperature of 10^6 GeV.

In helical hydromagnetic turbulence with an imposed magnetic field (which is constant in space and time) the magnetic helicity of the field within a periodic domain is no longer an invariant of the ideal equations. Alternatively, there is a generalized magnetic helicity that is an invariant of the ideal equations. It is shown that this quantity is not gauge invariant and that it can therefore not be used in practice. Instead, the evolution equation of the magnetic helicity of the field describing the deviation from the imposed field is shown to be a useful tool. It is demonstrated that this tool can determine steady state quenching of the alpha-effect. A simple three-scale model is derived to describe the evolution of the magnetic helicity and to predict its sign as a function of the imposed field strength. The results of the model agree favorably with simulations.

Decaying turbulence is studied numerically using as initial condition a random flow whose shell-integrated energy spectrum increases with wavenumber k like k^q. Alternatively, initial conditions are generated from a driven turbulence simulation by simply stopping the driving. It is known that the dependence of the decaying energy spectrum on wavenumber, time, and viscosity can be collapsed onto a unique scaling function that depends only on two parameters. This is confirmed using three-dimensional simulations and the dependence of the scaling function on its two arguments is determined.

The stability and conservation properties of a recently proposed polymerization model are studied. The achiral (racemic) solution is linearly unstable once the relevant control parameter (here the fidelity of the catalyst) exceeds a critical value. The growth rate is calculated for different fidelity parameters and cross-inhibition rates. A chirality parameter is defined and shown to be conserved by the nonlinear terms of the model. Finally, a truncated version of the model is used to derive a set of two ordinary differential equations and it is argued that these equations are more realistic than those used in earlier models of that form.

The dimensionless kinetic energy dissipation rate C_epsilon is estimated from numerical simulations of statistically stationary isotropic box turbulence that is slightly compressible. The Taylor microscale Reynolds number Re_lambda range is 20 < Re_lambda < 220 and the statistical stationarity is achieved with a random phase forcing method. The strong Re_lambda dependence of C_epsilon abates when Re_lambda approx. 100 after which C_epsilon slowly approaches approx 0.5 a value slightly different to previously reported simulations but in good agreement with experimental results. If C_epsilon is estimated at a specific time step from the time series of the quantities involved it is necessary to account for the time lag between energy injection and energy dissipation. Also, the resulting value can differ from the ensemble averaged value by up to +-30%. This may explain the spread in results from previously published estimates of C_epsilon.

The interaction between a protostellar magnetosphere and a surrounding dynamo-active accretion disc is investigated using an axisymmetric mean-field model. In all models investigated, the dynamo-generated magnetic field in the disc arranges itself such that in the corona, the field threading the disc is anti-aligned with the central dipole so that no X-point forms. When the magnetospheric field is strong enough (stellar surface field strength around 2 kG or larger), accretion happens in a recurrent fashion with periods of around 15 to 30 days, which is somewhat longer than the stellar rotation period of around 10 days. In the case of a stellar surface field strength of at least a few 100 G, the star is being spun up by the magnetic torque exerted on the star. The stellar accretion rates are always reduced by the presence of a magnetosphere which tends to divert a much larger fraction of the disc material into the wind. Both, a pressure-driven stellar wind and a disc wind form. In all our models with disc dynamo, the disc wind is structured and driven by magneto-centrifugal as well as pressure forces.

Evidence for non-Fickian diffusion of a passive scalar is presented using direct simulations of homogeneous isotropic turbulence. The results compare favorably with an explicitly time-dependent closure model based on the tau approximation. In the numerical experiments three different cases are considered: (i) zero mean concentration with finite initial concentration flux, (ii) an initial top hat profile for the concentration, and (iii) an imposed background concentration gradient. All cases agree in the resulting relaxation time in the tau approximation relating the triple correlation to the concentration flux. The first order smoothing approximation is shown to be inapplicable.

The excitation of gravity waves by penetrative convective plumes is investigated using 2D direct simulations of compressible convection. The oscillation field is measured by a new technique based on the projection of our simulation data onto the theoretical g-modes solutions of the associated linear eigenvalue problem. This allows us to determine both the excited modes and their corresponding amplitudes accurately.

Forced turbulence simulations are used to determine the turbulent kinematic viscosity, nu_t, from the decay rate of a large scale velocity field. Likewise, the turbulent magnetic diffusivity, eta_t, is determined from the decay of a large scale magnetic field. In the kinematic regime, when the field is weak, the turbulent magnetic Prandtl number, nu_t/eta_t, is about unity. When the field is nonhelical, eta_t is quenched when magnetic and kinetic energies become comparable. For helical fields the quenching is stronger and can be described by a dynamical quenching formula.

Local three-dimensional shearing box simulations of the compressible coupled dust-gas equations are used in the fluid approximation to study the evolution of different initial vortex configurations in a protoplanetary disc and their dust-trapping capabilities. The initial conditions for the gas are derived from an analytic solution to the compressible Euler equation and the continuity equation. The solution is valid if there is a vacuum outside the vortex. In the simulations the vortex is either embedded in a hot corona, or it is extended in a cylindrical fashion in the vertical direction. Both configurations are found to survive for at least one orbit and lead to accumulation of dust inside the vortex. This confirms earlier findings that dust accumulates in anticyclonic vortices, indicating that this is a viable mechanism for planetesimal formation.

It is argued that much of the observed magnetic helicity losses at the solar surface may represent a reduction of an otherwise more dominant nonlinearity of solar and stellar dynamos. This nonlinearity is proportional to the internal twist (as opposed to writhe) of helical and sigmoidal surface structures.

Nonhelical hydromagnetic turbulence without an imposed magnetic field is considered in the case where the magnetic Prandtl number is unity. The magnetic field is entirely due to dynamo action. The magnetic energy spectrum peaks at a wavenumber of about 5 times the minimum wavenumber in the domain, and not at the resistive scale, as has previously been argued. Throughout the inertial range the spectral magnetic energy exceeds the kinetic energy by a factor of about 2.5, and both spectra are approximately parallel. At first glance, the total energy spectrum seems to be close to k^{-3/2}, but there is a strong bottleneck effect and it is suggested that the asymptotic spectrum is k^{-5/3}. This is supported by the value of the second order structure function exponent that is found to be \zeta_2=0.70, suggesting a k^{-1.70} spectrum.

Numerical simulations of the multi-phase interstellar medium have been carried out, using a 3D, nonlinear, magnetohydrodynamic, shearing-box model, with random motions driven by supernova explosions. These calculations incorporate the effects of magnetic fields and rotation in 3D; these play important dynamical roles in the galaxy, but are neglected in many other simulations. The supernovae driving the motions are not arbitrarily imposed, but occur where gas accumulates into cold, dense clouds; their implementation uses a physically motivated model for the evolution of such clouds. The process is self-regulating, and produces mean supernova rates as part of the solution. Simulations with differing mean density show a power law relation between the supernova rate and density, with exponent 1.7; this value is within the range suggested from observations (taking star formation rate as a proxy for supernova rate). The global structure of the supernova driven medium is strongly affected by the presence of magnetic fields; e.g. for one solution the filling factor of hot gas is found to vary from 0.19 (with no field) to 0.12 (with initial mid-plane field B_0=6 muG).

Considerable progress has been made in understanding the hydrodynamics, but only to a certain extent the magnetohydrodynamics, of shaping bipolar outflows forming protoplanetary nebulae (PPNs) and planetary nebulae (PNs). In particular, Blackman et al. (2001b,2001a) point out two problems related to the formation of PNs and PPNs, regarding the formation of multipolar structures and the origin of the nebulae. They propose a solution by giving a semi-quantitative physical model which should be investigated by numerical simulations.

At numerical resolutions around $512^3$ and above, three-dimensional energy spectra from turbulence simulations begin to show noticeably shallower spectra than $k^{-5/3}$ near the dissipation wavenumber (`bottleneck effect'). This effect is shown to be significantly weaker in one-dimensional spectra such as those obtained in wind tunnel turbulence. The difference can be understood in terms of the transformation between one-dimensional and three-dimensional energy spectra under the assumption that the turbulent velocity field is isotropic. Transversal and longitudinal energy spectra are similar and can both accurately be computed from the full three-dimensional spectra. Second-order structure functions are less susceptible to the bottleneck effect and may be better suited for inferring the scaling exponent from numerical simulation data.

In the past few years suggestions have emerged that the solar magnetic field might have a bi-helical component with oppositely polarized magnetic fields at large and small scales, and that the shedding of such fields may be crucial for the operation of the dynamo. It is shown that, if a bi-helical field is shed into the solar wind, positive and negative contributions of the magnetic helicity spectrum tend to mix and decay. Even in the absence of turbulence, mixing and decay can occur on a time scale faster than the resistive one provided the two signs of magnetic helicity originate from a single tube. In the presence of turbulence, positively and negatively polarized contributions mix rapidly in such a way that the ratio of magnetic helicity to magnetic energy is largest both at the largest scale and in the dissipation range. In absolute units the small scale excess of helical fields is however negligible.

An axisymmetric model of a cool, dynamo-active accretion disc is applied to protostellar discs. Thermally and magnetically driven outflows develop that are not collimated within 0.1 AU. In the presence of a central magnetic field from the protostar, accretion onto the protostar is highly episodic, which is in agreement with earlier work.

According to the kinematic theory of nonhelical dynamo action the magnetic energy spectrum increases with wavenumber and peaks at the resistive cutoff wavenumber. It has previously been argued that even in the dynamical case the magnetic energy peaks at the resistive scale. Using high resolution simulations (up to 1024^3 meshpoints) with no large scale imposed field we show that the magnetic energy peaks at a wavenumber that is independent of the magnetic Reynolds number and about 5 times larger than the forcing wavenumber. Throughout the inertial range the spectral magnetic energy exceeds the kinetic energy by a factor of 2 to 3. Both spectra are approximately parallel. The total energy spectrum seems to be close to k^{-3/2}, but there is a strong bottleneck effect and it is suggested that the asymptotic spectrum is instead k^{-5/3}. This is supported by the value of the second order structure function exponent that is found to be zeta_2=0.70, suggesting a k^{-1.70} spectrum. The third order structure function scaling exponent is very close to unity, in agreement with Goldreich-Sridhar theory. Adding an imposed field tends to suppress the small scale magnetic field. We find that at large scales the magnetic power spectrum follows then a k^{-1} slope. When the strength of the imposed field is of the same order as the dynamo generated field, we find almost equipartition between the magnetic and kinetic power spectra.

There is mounting evidence that the ejection of magnetic helicity from the solar surface is important for the solar dynamo. Observations suggest that in the northern hemisphere the magnetic helicity flux is negative. We propose that this magnetic helicity flux is mostly due to small scale magnetic fields; in contrast to the more systematic large scale field of the 11 year cycle, whose helicity flux may be of opposite sign, and may be excluded from the observational interpretation. Using idealized simulations of MHD turbulence as well as a simple two-scale model, we show that shedding small scale (helical) field has two important effects. (i) The strength of the large scale field reaches the observed levels. (ii) The evolution of the large scale field proceeds on time scales shorter than the resistive time scale, as would otherwise be enforced by magnetic helicity conservation. In other words, the losses ensure that the solar dynamo is always in the near-kinematic regime. This requires, however, that the ratio of small scale to large scale losses cannot be too small, for otherwise the large scale field in the near-kinematic regime will not reach the observed values.

Isotropic homogeneous hydromagnetic turbulence is studied using numerical simulations at resolutions of up to 1024^3 meshpoints. It is argued that, in contrast to the kinematic regime, the nonlinear regime is characterized by a spectral magnetic power that is decreasing with increasing wavenumber, regardless of whether or not the turbulence has helicity. This means that the root-mean-square field strength converges to a limiting value at large magnetic Reynolds numbers. The total (magnetic and kinetic) energy spectrum tends to be somewhat shallower than k^{-5/3}, in agreement with the findings of other groups. In the presence of helicity, an inverse cascade develops, provided the scale separation between the size of the computational box and the scale of the energy carrying eddies exceeds a ratio of at least two. Finally, the constraints imposed by magnetic helicity conservation on mean-field theory are reviewed and new results of simulations are presented.

We study the evolution of a stochastic helical magnetic field generated in the early Universe after the electroweak phase transition, using standard magnetohydrodynamics (MHD). We find how the coherence length xi, magnetic energy E_M and magnetic helicity H evolve with time. We show that the self-similarity of the magnetic power spectrum alone implies that xi ~ t^{1/2}. This in turn implies that magnetic helicity decays as H ~ t^{-2s}, and that the magnetic energy decays as E_M ~ t^{-0.5-2s}, where s is inversely proportional to the magnetic Reynolds number Re_M. These laws improve on several previous estimates.

Statistical parameters of the ISM driven by thermal energy injections from supernova explosions have been obtained from 3D, nonlinear, magnetohydrodynamic, shearing-box simulations for spiral arm and interarm regions. The density scale height obtained for the interarm regions is 50% larger than within the spiral arms because of the higher gas temperature. The filling factor of the hot gas is also significantly larger between the arms and depends sensitively on magnetic field strength.

Over the past few years there has been growing interest in helical magnetic field structures seen at the solar surface, in coronal mass ejections, as well as in the solar wind. Although there is a great deal of randomness in the data, on average the extended structures are mostly left-handed on the northern hemisphere and right-handed on the southern. Surface field structures are also classified as dextral (= right bearing) and sinistral (= left bearing) occurring preferentially in the northern and southern hemispheres respectively. Of particular interest here is a quantitative measurement of the associated emergence rates of helical structures, which translate to magnetic helicity fluxes. In this review, we give a brief survey of what has been found so far and what is expected based on models. Particular emphasis is put on the scale dependence of the associated fields and an attempt is made to estimate the helicity flux of the mean field vs. fluctuating field.

Two questions about the solar magnetic field might be answered together once their connection is identified. The first is important for large scale dynamo theory: what prevents the magnetic backreaction forces from shutting down the dynamo cycle? The second question is: what determines the handedness of twist and writhe in magnetized coronal ejecta? Magnetic helicity conservation is important for answering both questions. Conservation implies that dynamo generation of large scale writhed structures is accompanied by the oppositely signed twist along these structures. The latter is associated with the backreaction force. We suggest that coronal mass ejections (CME's) simultaneously liberate small scale twist and large scale writhe of opposite sign, helping to prevent the cycle from quenching and enabling a net magnetic flux change in each hemisphere. Observations and helicity spectrum measurements from a simulation of a rising flux ribbon support this idea. We show a new pictorial of dynamo flux generation that includes the backreaction and magnetic helicity conservation by characterizing the field as a 2-D ribbon rather than a 1-D line.

The principle of the Karlsruhe dynamo experiment is closely related to that of the Roberts dynamo working with a simple fluid flow which is, with respect to proper Cartesian co-ordinates x, y and z, periodic in x and y and independent of z. A modified Roberts dynamo problem is considered with a flow more similar to that in the experimental device. Solutions are calculated numerically, and on this basis an estimate of the excitation condition of the experimental dynamo is given. The modified Roberts dynamo problem is also considered in the framework of the mean-field dynamo theory, in which the crucial induction effect of the fluid motion is an anisotropic alpha-effect. Numerical results are given for the dependence of the mean-field coefficients on the fluid flow rates. The excitation condition of the dynamo is also discussed within this framework. The behavior of the dynamo in the nonlinear regime, i.e. with backreaction of the magnetic field on the fluid flow, depends on the effect of the Lorentz force on the flow rates. The quantities determining this effect are calculated numerically. The results for the mean-field coefficients and the quantities describing the backreaction provide corrections to earlier results, which were obtained under simplifying assumptions.

We present an axisymmetric numerical model of a dynamo active accretion disc. If the dynamo-generated magnetic field in the disc is sufficiently strong (close to equipartition with thermal energy), a fast magneto-centrifugally driven outflow develops within a conical shell near the rotation axis, together with a slower pressure driven outflow from the outer parts of the disc as well as around the axis. Our results show that a dynamo active accretion disc can contribute to driving an outflow even without any external magnetic field. The fast outflow in the conical shell is confined by the azimuthal field which is produced by the dynamo in the disc and advected to the disc corona. This part of the outflow has high angular momentum and is cooler and less dense than its surroundings. The conical shell's half-opening angle is typically about 30 degrees near the disc and decreases slightly with height. The slow outflow is hotter and denser.

Various approaches to estimate turbulent transport coefficients from numerical simulations of hydromagnetic turbulence are discussed. A quantitative comparison between the averaged magnetic field obtained from a specific three-dimensional simulation of a rotating turbulent shear flow in a slab and a simple one-dimensional alpha-omega dynamo model is given. A direct determination of transport coefficients is attempted by calculating the correlation matrix of different components of the field and its derivatives. The alpha-effect operating on the toroidal field is found to be negative and of similar magnitude as the value obtained previously by correlating the electromotive force with the mean magnetic field. The turbulent diffusion of the toroidal field is comparable to the kinematic viscosity that was determined earlier by comparing the stress with the shear. However, the turbulent diffusion of the radial field component is smaller and can even be formally negative. The method is then modified to obtain the spectral dependence of the turbulent transport coefficients on the wavenumber. There is evidence for nonlocal behavior in that most of the response comes from the smallest wavenumbers corresponding to the largest scale possible in the simulation.

We explore the dependence of the amplitude of stellar dynamo cycle variability (as seen in the Mount Wilson Ca II HK timeseries data) on other stellar parameters. We find that the fractional cycle amplitude A_cyc (i.e. the ratio of the peak-to-peak variation to the average) decreases somewhat with mean activity, increases with decreasing effective temperature, but is not correlated with inverse Rossby number Ro^{-1}. We find that A_cyc increases with the ratio of cycle and rotational frequencies omega_cyc/Omega along two, nearly parallel branches.

The connection between helically isotropic MHD turbulence and mean-field dynamo theory is reviewed. The nonlinearity in the mean-field theory is not yet well established, but detailed comparison with simulations begin to help select viable forms of the nonlinearity. The crucial discriminant is the magnetic helicity, which is known to evolve only on a slow resistive time scale in the limit of large magnetic Reynolds number. Particular emphasis is put on the possibility of memory effects, which means that an additional explicitly time-dependent equation for the nonlinearity is solved simultaneously with the mean-field equations. This approach leads to better agreement with the simulations, while it would also produce more favorable agreement between models and stellar dynamos.

Recent progress in the theory of solar and stellar dynamos is reviewed. Particular emphasis is placed on the mean-field theory which tries to describe the collective behavior of the magnetic field. In order to understand solar and stellar activity, a quantitatively reliable theory is necessary. Much of the new developments center around magnetic helicity conservation which is seen to be important in numerical simulations. Only a dynamical, explicitly time dependent theory of alpha-quenching is able to describe this behavior correctly.

We supplement the mean field dynamo growth equation with the total magnetic helicity evolution equation. This provides an explicitly time dependent model for alpha quenching in dynamo theory. For dynamos without shear, this approach accounts for the observed large scale field growth and saturation in numerical simulations. After a significant kinematic phase, the dynamo is resistively quenched, i.e. the saturation time depends on the microscopic resistivity. This is independent of whether or not the turbulent diffusivity is resistively quenched. We find that the approach is also successful for dynamos that include shear and exhibit migratory waves (cycles). In this case however, whether or not the cycle period remains of the order of the dynamical time scale at large magnetic Reynolds numbers does depend how on how the turbulent magnetic diffusivity quenches. Since this is unconstrained by magnetic helicity conservation, the diffusivity is presently an input parameter. Comparison to current numerical experiments suggests a turbulent diffusivity that depends only weakly on the magnetic Reynolds number, but higher resolution simulations are needed.

We study the pumping of magnetic flux in three-dimensional compressible magnetoconvection in the context of stellar dynamos. The simulation domain represents a rectangular section from the lower part of a stellar convection zone plus the underlying stably stratified layer, with a total depth of up to five pressure scale heights. Once convection has attained a statistically stationary state, a magnetic field is introduced. The magnetic field is subsequently modified by the convective motions, and the resulting pumping effects are isolated by calculating various coefficients of the expansion of the electromotive force, uxb, in terms of components of the mean magnetic field. The dependence of the pumping effects on rotation, latitude and other parameters is studied. First numerical evidence is found for the existence of pumping effects in the horizontal directions, unless the rotation axis coincides with the vertical axis, as is the case on the poles. Evidence is found that the pumping effects act differently on different components of the mean magnetic field. Latitudinal pumping is mainly equatorward for toroidal field, and can be poleward for poloidal field. Longitudinal pumping is mainly retrograde for the radial field but prograde for the latitudinal field. The pumping effect in the vertical direction is found to be dominated by the diamagnetic effect, equivalent to a predominating downward advection with a maximum speed in the turbulent case of about 10 percent of the rms convective velocity. Where possible, an attempt is made to identify the physical origin of the effect. Finally, some consequences of the results for stellar dynamos are discussed.

The theory of large scale dynamos is reviewed with particular emphasis on the magnetic helicity constraint in the presence of closed and open boundaries. In the presence of closed or periodic boundaries, helical dynamos respond to the helicity constraint by developing small scale separation in the kinematic regime, and by showing long time scales in the nonlinear regime where the scale separation has grown to the maximum possible value. A resistively limited evolution towards saturation is also found at intermediate scales before the largest scale of the system is reached. Larger aspect ratios can give rise to different structures of the mean field which are obtained at early times, but the final saturation field strength is still decreasing with decreasing resistivity. In the presence of shear, cyclic magnetic fields are found whose period is increasing with decreasing resistivity, but the saturation energy of the mean field is in strong super-equipartition with the turbulent energy. It is shown that artificially induced losses of small scale field of opposite sign of magnetic helicity as the large scale field can, at least in principle, accelerate the production of large scale (poloidal) field. Based on mean field models with an outer potential field boundary condition in spherical geometry, we verify that the sign of the magnetic helicity flux from the large scale field agrees with the sign of alpha. For solar parameters, typical magnetic helicity fluxes lie around 10^{47} Mx^2 per cycle.

The self-excitation of magnetic field by a spiral Couette flow between two coaxial cylinders is considered. We solve numerically the fully nonlinear, three-dimensional MHD equations for magnetic Prandtl numbers Pm (ratio of kinematic viscosity to magnetic diffusivity) between 0.14 and 10 and kinematic and magnetic Reynolds numbers up to about 2000. In the initial stage of exponential field growth (kinematic dynamo regime), we find that the dynamo switches from one distinct regime to another as the radial width \delta r_B$ of the magnetic field distribution becomes smaller than the separation of the field maximum from the flow boundary. The saturation of magnetic field growth is due to a reduction in the velocity shear resulting mainly from the azimuthally and axially averaged part of the Lorentz force, which agrees with an asymptotic result for the limit of Pm<<1. In the parameter regime considered, the magnetic energy decreases with kinematic Reynolds number as Re^{-0.84}, which is approximately as predicted by the nonlinear asymptotic theory (\sim Re^{-1}). However, when the velocity field is maintained by a volume force (rather than by viscous stress) the dependence of magnetic energy on the kinematic Reynolds number is much weaker.

We perform direct numerical simulations of three-dimensional freely decaying magnetohydrodynamic (MHD) turbulence. For helical magnetic fields an inverse cascade effect is observed in which power is transfered from smaller scales to larger scales. The magnetic field reaches a scaling regime with self-similar evolution, and power law behavior at high wavenumbers. We also find power law decay in the magnetic and kinematic energies, and power law growth in the characteristic length scale of the magnetic field.

We have carried out numerical simulations of freely decaying magnetohydrodynamic (MHD) turbulence in three dimensions, which can be applied to the evolution of stochastic magnetic fields in the early Universe. For helical magnetic fields an inverse cascade effect is observed in which magnetic helicity and energy is transfered from smaller scales to larger scales, accompanied by power law growth in the characteristic length scale of the magnetic field. The magnetic field quickly reaches a scaling regime with self-similar evolution, and power law behaviour at high wavenumbers. We also find power law decay in the magnetic and kinematic energies.

In ordinary turbulence research it has been a long standing tradition to solve the equations in spectral space giving the best possible accuracy. This is indeed a natural choice for incompressible problems with periodic boundaries, but it is no longer optimal in many astrophysical circumstances. It is argued that lower order spatial derivatives schemes are unacceptable in view of their low overall accuracy, even when mass, momentum, and energy are conserved to machine accuracy. High order finite difference schemes are therefore found to be quite efficient and physically appropriate. They are also easily and efficiently implemented on massively parallel computers. High order schemes also yield sufficient overall accuracy. Our code uses centered finite differences which make the adaptation to other problems simple. Since the code is not written in conservative form, conservation of mass, energy and momentum can be used to monitor to quality of the solution. A third order Runge-Kutta scheme with 2N-storage is used for calculating the time advance.

In numerical studies of turbulence, hyperviscosity is often used as a tool to extend the inertial subrange and to reduce the dissipative subrange. By analogy, hyperdiffusivity (or hyperresistivity) is sometimes used in magnetohydrodynamics. The underlying assumption is that only the small scales are affected by this manipulation. In the present paper, possible side effects on the evolution of the large scale magnetic field are investigated. It is found that for turbulent flows with helicity, hyperdiffusivity causes the dynamo-generated magnetic field to saturate at a higher level than normal diffusivity. This result is successfully interpreted in terms of magnetic helicity conservation, which also predicts that full saturation is only reached after a time comparable to the large scale magnetic (hyper)diffusion time.

The advantages of high-order finite difference scheme for astrophysical MHD and turbulence simulations are highlighted. A number of one-dimensional test cases are presented ranging from various shock tests to Parker-type wind solutions. Applications to magnetized accretion discs and their associated outflows are discussed. Particular emphasis is placed on the possibility of dynamo action in three-dimensional turbulent convection and shear flows, which is relevant to stars and astrophysical discs. The generation of large scale fields is discussed in terms of an inverse magnetic cascade and the consequences imposed by magnetic helicity conservation are reviewed with particular emphasis on the issue of alpha-quenching.

We present numerical simulations of three-dimensional compressible magnetoconvection in a rotating rectangular box that represents a section of the solar convection zone. The box contains a convectively unstable layer, surrounded by stably stratified layers with overshooting convection. The magnetic Reynolds number, Rm, is chosen subcritical, thus excluding spontaneous growth of the magnetic field through dynamo action, and the magnetic energy is maintained by introducing a constant magnetic field into the box, once convection has attained a statistically stationary state. Under the influence of the Coriolis force, the advection of the magnetic field results in a non-vanishing contribution to the mean electric field, given by uxb. From this electric field, we calculate the alpha-effect, separately for the stably and the unstably stratified layers, by averaging over time and over suitably defined volumes. From the variation of alpha we derive an error estimate, and the dependence of alpha on rotation and magnetic field strength is studied. Evidence is found for rotational quenching of the vertical alpha-effect, and for a monotonic increase of the horizontal alpha-effect with increasing rotation. For Rm~30, our results for both vertical and horizontal alpha-effect are consistent with magnetic quenching by a factor 1/[1+Rm(B_0/B_eq)^2]. The signs of the small-scale current helicity and of the vertical component of alpha are found to be opposite to those for isotropic turbulence.

The possibility of non-helical large scale dynamo action is investigated using three-dimensional simulations of global accretion discs as well as idealized local simulations without rotation and only shear. Particular emphasis is placed on a certain correlation between vorticity and azimuthal velocity gradient which has been predicted to drive large scale dynamo action, independent of the presence or absence of kinetic helicity. In the global disc simulations two types of behaviours are found: those which do show this type of velocity correlation and those which do not. The former ones are typically also the cases where the resistivity is larger. The latter ones show signs typical of dynamo action based on the usual helicity effect. In the idealized simulations without rotation and just shear the above correlation is found to be particularly strong. In both cases there is, as expected, a systematic flux of magnetic helicity through the midplane. However, very little helicity flux leaves the domain through the top and bottom boundaries. The idealized simulations reveal that much of this systematic flux comes from the rotational component of the helicity flux and does not contribute to its divergence.

The generation of large scale flows by the anisotropic kinetic alpha (AKA) effect is investigated in simulations with a suitable time-dependent space- and time-periodic anisotropic forcing lacking parity invariance. The forcing pattern moves relative to the fluid, which leads to a breaking of the Galilean invariance as required for the AKA effect to exist. The AKA effect is found to produce a clear large scale flow pattern when the Reynolds number, R, is small as only a few modes are excited in linear theory. In this case the non-vanishing components of the AKA tensor are dynamically independent of the Reynolds number. For larger values of R, many more modes are excited and the components of the AKA tensor are found to decrease rapidly with increasing value of R. However, once there is a magnetic field (imposed and of sufficient strength, or dynamo-generated and saturated) the field begins to suppress the AKA effect, regardless of the value of R. It is argued that the AKA effect is unlikely to be astrophysically significant unless the magnetic field is weak and R is small.

We further investigate relationships between activity cycle periods in cool stars and rotation to include new cycle data, and explore different parameterizations of the problem. We find that relations between cycle and rotational frequencies (omega_cyc vs. Omega) and between their ratio and the inverse Rossby number (omega_cyc/Omega vs. Ro^{-1}) show many similarities, including three branches and similar rms scatter. We briefly discuss some implications for dynamo models.

Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turbulent stress acting on the epicyclic motion is comparable in size to the stress that drives the accretion, however an important ingredient in the damping of the epicyclic motion is its parametric decay into inertial waves.

We consider the effect of vertical outflows on the mean-field dynamo in a thin disk. These outflows could be due to winds or magnetic buoyancy. We analyse both two-dimensional finite-difference numerical solutions of the axisymmetric dynamo equations and a free-decay mode expansion using the thin-disk approximation. Contrary to expectations, a vertical velocity can enhance dynamo action, provided it is not too strong. In the nonlinear regime this can lead to super-exponential growth of the magnetic field.

The buoyant rise of thermals (i.e. bubbles of enhanced entropy, but initially in pressure equilibrium) is investigated numerically in three dimensions for the case of an adiabatically stratified layer covering 6-9 pressure scale heights. It is found that these bubbles can travel to large heights before being braked by the excess pressure that builds up in order to drive the gas sideways in the head of the bubble. Until this happens, the momentum of the bubble grows as described by the time-integrated buoyancy force. This validates the simple theory of bubble dynamics whereby the mass entrainment of the bubble provides an effective braking force well before the bubble stops ascending. This is quantified by an entrainment parameter alpha which is calculated from the simulations and is found to be in good agreement with the experimental measurements. This work is discussed in the context of contact binaries whose secondaries could be subject to dissipative heating in the outermost layers.

Dynamo action is investigated in simulations of locally isotropic and homogeneous turbulence in a slab between open boundaries. It is found that a `pseudo-vacuum' boundary condition (where the field is vertical) leads to strong helicity fluxes which significantly reduce the amplitude of the resulting large scale field. On the other hand, if there is a conducting halo outside the dynamo-active region the large scale field amplitude can reach larger values, but the time scale after which this field is reached increases linearly with the magnetic Reynolds number. In both cases most of the helicity flux is found to occur on large scales. From the variety of models considered we conclude that open boundaries tend to lower the saturation field strength compared to the case with periodic boundaries. The rate at which this lower saturation field strength is attained is roughly independent of the strength of the turbulence and of the boundary conditions. For dynamos with less helicity, however, significant field strengths could be reached in a shorter time.

The evolution of magnetic fields is studied using simulations of forced helical turbulence with strong imposed shear. After some initial exponential growth, the magnetic field develops a large scale travelling wave pattern. The resulting field structure possesses magnetic helicity, which is conserved in a periodic box by the ideal MHD equations and can hence only change on a resistive time scale. This constrains strongly the growth time of the large scale magnetic field, but less strongly the length of the cycle period. Comparing with the case without shear, the time scale for large scale field amplification is shortened by a factor Q, which depends on the relative importance of shear and helical turbulence, and which controls also the ratio of toroidal to poloidal field. The results of the simulations can be reproduced qualitatively and quantitatively with a mean-field alpha-Omega dynamo model with alpha-effect and the turbulent magnetic diffusivity coefficients that are less strongly quenched than in the corresponding alpha^2-dynamo.

The effect of shear on the growth of large scale magnetic fields in helical turbulence is investigated. The resulting large-scale magnetic field is also helical and continues to evolve, after saturation of the small scale field, on a slow resistive time scale. This is a consequence of magnetic helicity conservation. Because of shear, the time scale needed to reach an equipartition-strength large scale field is shortened proportionally to the ratio of the resulting toroidal to poloidal large scale fields.

Helicity, a measure of the linkage of flux lines, has subtle and largely unknown effects upon dynamics. Both magnetic and hydrodynamic helicity are conserved for ideal systems and could suppress nonlinear dynamics. What actually happens is not clear because in a fully three-dimensional system there are additional channels whereby intense, small-scale dynamics can occur. This contribution shows one magnetic and one hydrodynamic case where for each the presence of helicity does not suppress small-scale intense dynamics of the type that might lead to reconnection.

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.

The emergence of a large scale magnetic field from randomly forced isotropic strongly helical flows is discussed in terms of the inverse cascade of magnetic helicity and the alpha-effect. In simulations of such flows the maximum field strength exceeds the equipartition field strength for large scale separation. However, helicity conservation controls the speed at which this final state is reached. In the presence of open boundaries magnetic helicity fluxes out of the domain are possible. This reduces the timescales of the field growth, but it also tends to reduce the maximum attainable field strength.

We study the effects of vertical outflows on mean-field dynamos in disks. These outflows could be due to thermal winds or magnetic buoyancy. We analyse numerical solutions of the nonlinear mean-field dynamo equations using a two-dimensional finite-difference model. Contrary to expectations, a modest vertical velocity can enhance dynamo action. This can lead to super-exponential growth of the magnetic field and to higher magnetic energies at saturation in the nonlinear regime.

A number of problems of solar and stellar dynamo theory are briefly reviewed and the current status of possible solutions is discussed. Results of direct numerical simulations are described in view of mean-field dynamo theory and the relation between the alpha-effect and the inverse cascade of magnetic helicity is highlighted. The possibility of `catastrophic' quenching of the alpha-effect is explained in terms of the constraint placed by the conservation of magnetic helicity.

Using a periodic box calculation it is shown that, owing to helicity conservation, a large scale field can only develop on a resistive timescale. This behaviour can be reproduced by a mean-field dynamo with alpha and eta_t quenchings that are equally strong and `catastrophic'.

The dynamical friction experienced by a body moving in a gaseous medium is different from the friction in the case of a collisionless stellar system. Here we consider the orbital evolution of a gravitational perturber inside a gaseous sphere using three-dimensional simulations, ignoring however self-gravity. The results are analysed in terms of a `local' formula with the associated Coulomb logarithm taken as a free parameter. For forced circular orbits, the asymptotic value of the component of the drag force in the direction of the velocity is a slowly varying function of the Mach number in the range 1.0-1.6. The dynamical friction timescale for free decay orbits is typically only half as long as in the case of a collisionless background, which is in agreement with E.C. Ostriker's recent analytic result. The orbital decay rate is rather insensitive to the past history of the perturber. It is shown that, similar to the case of stellar systems, orbits are not subject to any significant circularization. However, the dynamical friction timescales are found to increase with increasing orbital eccentricity for the Plummer model, whilst no strong dependence on the initial eccentricity is found for the isothermal sphere.

It is shown that ambipolar diffusion as a toy nonlinearity leads to very similar behaviour of large scale turbulent dynamos as full MHD. This is demonstrated using both direct simulations in a periodic box and a closure model for the magnetic correlation functions applicable to infinite space. Large scale fields develop via a nonlocal inverse cascade as described by the alpha-effect. However, because magnetic helicity can only change on a resistive timescale, the time it takes to organize the field into large scales increases with magnetic Reynolds number.

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, $\alpha_{\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, $\alpha_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $\alpha_{\rm v}$ is approximately equal to $\alpha_{\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 $\alpha_{\rm v}$ and $\alpha_{\rm h}$ we conclude that for $\beta = 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.

Numerical evidence for a finite-time singularity in ideal 3D magnetohydrodynamics (MHD) is presented. The simulations start from two interlocking magnetic flux rings with no initial velocity. The magnetic curvature force causes the flux rings to shrink until they come into contact. This produces a current sheet between them. In the ideal compressible calculations, the evidence for a singularity in a finite time $t_c$ is that the peak current density behaves like $|J|_\infty \sim 1/(t_c-t)$ for a range of sound speeds (or plasma betas). For the incompressible calculations consistency with the compressible calculations is noted and evidence is presented that there is convergence to a self-similar state. In the resistive reconnection calculations the magnetic helicity is nearly conserved and energy is dissipated.

We study the changes in the dynamical behaviour of axisymmetric spherical mean-field dynamo models produced by changes in their geometry and topology, by considering a two parameter family of models, ranging from a full sphere to spherical shell, torus and disc-like configurations, within a unified framework. We find that the two parameter space of the family of models considered here separates into at least three different regions with distinct characteristics for the onset of dynamo action. In two of these regions, the most easily excited fields are oscillatory, in one case with dipolar symmetry, and in the other with quadrupolar, whereas in the third region the most easily excited field is steady and quadrupolar. In the nonlinear regime, we find that topological changes can alter significantly the dynamical behaviour, whilst modest changes in geometry can produce qualitative changes, particularly for thin disc-like configurations. This is of potential importance, since the exact shapes of astrophysical bodies, especially accretion discs and galaxies, are usually not precisely known.

Calculations for mean field dynamo models (in both full spheres and spherical shells), with both algebraic and dynamic $\alpha$--quenchings, show qualitative as well as quantitative differences and similarities in the dynamical behaviour of these models. We summarise and enhance recent results with extra examples. Overall, the effect of using a dynamic $\alpha$ appears to be complicated and is affected by the region of parameter space examined.

We study axisymmetric mean field spherical and spherical shell dynamo models, with both dynamic and algebraic $\alpha$--quenchings. Our results show that there are qualitative as well as quantitative differences and similarities between these models. Regarding similarities, both groups of models exhibit symmetric, antisymmetric and mixed modes of behaviour. As regards differences, the important feature in the full sphere models is the occurrence of chaotic behaviour in the algebraic $\alpha$--quenching models. For the spherical shell models with dynamic $\alpha$ the main features include the possibility of multi-attractor regimes with final state sensitivity with respect to small changes in the magnitude of $\alpha$ and the initial parity. We find the effect of introducing a dynamic $\alpha$ is likely to be complicated and depend on the region of the parameter space considered, rather than a uniform change towards simplicity or complexity.

We investigate the behaviour of $\alpha\Omega$ dynamos with a dynamic $\alpha$, whose evolution is governed by the imbalance between a driving and a damping term. We focus on truncated versions of such dynamo models which are often studied in connection with solar and stellar variability. Given the approximate nature of such models, it is important to study how robust they are with respect to reasonable changes in the formulation of the driving and damping terms. For each case, we also study the effects of changes of the dynamo number and its sign, the truncation order and initial conditions. Our results show that changes in the formulation of the driving term have important consequences for the dynamical behaviour of such systems, with the detailed nature of these effects depending crucially on the form of the driving term assumed, the value and the sign of the dynamo number and the initial conditions. On the other hand, the change in the damping term considered here seems to produce little qualitative effect.

In a recent work (Covas et al. 1997), the behaviour and the robustness of truncated $\alpha\Omega$ dynamos with a dynamic $\alpha$ were studied with respect to a number of changes in the driving term of the dynamic $\alpha$ equation, which was considered previously by Schmalz and Stix (1991) to be of the form $\sim A_{\phi}B_{\phi}$. Here we review and extend our previous work and consider the effect of adding a quadratic quenching term of the form $\alpha|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.

We study the effects of plasma viscosity on the dynamics of primordial magnetic fields by simulating magnetohydrodynamics in the early universe by appropriate non-linear cascade models. We find numerically that even in the presence of large kinetic viscosity, magnetic energy is transferred to large length scales. There are indications, however, that the inverse cascade stops at a given time which depends on the magnitude of viscosity. For realistic viscosities we do not find equipartition between magnetic and kinetic energies.

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.

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.