This issue commemorates an outstanding scientist of the twentieth century, Yakov Borisovich Zeldovich, in connection with the centenary of his birth (8 March 1914), with a collection of reviews and research articles broadly related to large-scale random phenomena in astrophysical plasmas.

A side remark in the old paper by Zeldovich and co-authors leads to the recent discovery of a universal conservation law of turbulent dispersion.

We briefly review some aspects of magnetic turbulence intermittency observed in space plasmas. Deviation of statistical characteristics of a system (e.g. its high statistical momenta) from the Gaussian can manifest itself as domination of rare large intensity peaks often associated with the intermittency in the system's dynamics. Thirty years ago, Zeldovich stressed the importance of the non-Gaussian appearance of the sharp values of vector and scalar physical parameters in random media as a factor of magnetic field amplification in cosmic structures. Magnetic turbulence is governing the behavior of collisionless plasmas in space and especially the physics of shocks and magnetic reconnections. Clear evidence of intermittent magnetic turbulence was found in recent in-situ spacecraft measurements of magnetic fields in the near-Earth and interplanetary plasma environments. We discuss the potentially promising approaches of incorporating the knowledge gained from spacecraft in-situ measurements into modern models describing plasma dynamics and radiation in various astrophysical systems. As an example, we discuss supernova remnants (SNRs) which are known to be the sources of energy, momentum, chemical elements, and high-energy cosmic rays (CRs) in galaxies. Supernova shocks accelerate charged particles to very high energies and may strongly amplify turbulent magnetic fields via instabilities driven by CRs. Relativistic electrons accelerated in SNRs radiate polarized synchrotron emission in a broad range of frequencies spanning from the radio to gamma-rays. We discuss the effects of intermittency of magnetic turbulence on the images of polarized synchrotron X-ray emission of young SNRs and emission spectra of pulsar wind nebula.

A spectacular phenomenon of intermittency, i.e. a progressive growth of higher statistical moments of a physical field excited by an instability in a random medium, attracted the attention of Zeldovich in the last years of his life. At that time, the mathematical aspects underlying the physical description of this phenomenon were still under development and relations between various findings in the field remained obscure. Contemporary results from the theory of the product of independent random matrices (the Furstenberg theory) allowed the elaboration of the phenomenon of intermittency in a systematic way. We consider applications of the Furstenberg theory to some problems in cosmology and dynamo theory.

In this paper, we investigate the evolution of braided solar coronal loops. We assume that coronal loops consist of several internal strands which twist and braid about each other. Reconnection between the strands leads to small flares and heating of the loop to x-ray temperatures. Using a method of generating and releasing braid structure similar to a forest fire model, we show that the reconnected field lines evolve to a self-organised critical state. In this state, the frequency distributions of coherent braid sequences as well as flare energies follow power law distributions. We demonstrate how the presence of net helicity in the loop alters the distribution laws.

In this paper we review a recently developed approximate method for investigation of dynamics of compressible ellipsoidal figures. Collapse and subsequent behaviour are described by a system of ordinary differential equations for time evolution of semi-axes of a uniformly rotating, three-axis, uniform-density ellipsoid. First, we apply this approach to investigate dynamic stability of non-spherical bodies. We solve the equations that describe, in a simplified way, the Newtonian dynamics of a self-gravitating non-rotating spheroidal body. We find that, after loss of stability, a contraction to a singularity occurs only in a pure spherical collapse, and deviations from spherical symmetry prevent the contraction to the singularity through a stabilizing action of nonlinear non-spherical oscillations. The development of instability leads to the formation of a regularly or chaotically oscillating body, in which dynamical motion prevents the formation of the singularity. We find regions of chaotic and regular pulsations by constructing a Poincaré diagram. A real collapse occurs after damping of the oscillations because of energy losses, shock wave formation or viscosity. We use our approach to investigate approximately the first stages of collapse during the large scale structure formation. The theory of this process started from ideas of Ya. B. Zeldovich, concerning the formation of strongly non-spherical structures during nonlinear stages of the development of gravitational instability, known as ‘Zeldovich’s pancakes’. In this paper the collapse of non-collisional dark matter and the formation of pancake structures are investigated approximately. Violent relaxation, mass and angular momentum losses are taken into account phenomenologically. We estimate an emission of very long gravitational waves during the collapse, and discuss the possibility of gravitational lensing and polarization of the cosmic microwave background by these waves.

Fluctuation dynamos are generic to turbulent astrophysical systems. The only analytical model of the fluctuation dynamo, due to Kazantsev, assumes the velocity to be delta-correlated in time. This assumption breaks down for any realistic turbulent flow. We generalize the analytic model of fluctuation dynamos to include the effects of a finite correlation time, ${\it\tau}$, using renewing flows. The generalized evolution equation for the longitudinal correlation function $M_{L}$ leads to the standard Kazantsev equation in the ${\it\tau}\rightarrow 0$ limit, and extends it to the next order in ${\it\tau}$. We find that this evolution equation also involves third and fourth spatial derivatives of $M_{L}$, indicating that the evolution for finite-${\it\tau}$ will be non-local in general. In the perturbative case of small-${\it\tau}$ (or small Strouhal number), it can be recast using the Landau–Lifschitz approach, to one with at most second derivatives of $M_{L}$. Using both a scaling solution and the WKBJ approximation, we show that the dynamo growth rate is reduced when the correlation time is finite. Interestingly, to leading order in ${\it\tau}$, we show that the magnetic power spectrum preserves the Kazantsev form, $M(k)\propto k^{3/2}$, in the large-$k$ limit, independent of ${\it\tau}$.

Zeldovich’s stretch–twist–fold (STF) dynamo provided a breakthrough in conceptual understanding of fast dynamos, including the small-scale fluctuation dynamos. We study the evolution and saturation behaviour of two types of generalized Baker’s map dynamos, which have been used to model Zeldovich’s STF dynamo process. Using such maps allows one to analyse dynamos at much higher magnetic Reynolds numbers $\mathit{Re}_{M}$ as compared to direct numerical simulations. In the two-strip map dynamo there is constant constructive folding, while the four-strip map dynamo also allows the possibility of a destructive reversal of the field. Incorporating a diffusive step parametrized by $\mathit{Re}_{M}$ into the map, we find that the magnetic field $B(x)$ is amplified only above a critical $\mathit{Re}_{M}=R_{\mathit{crit}}\sim 4$ for both types of dynamos. The growing $B(x)$ approaches a shape-invariant eigenfunction independent of initial conditions, whose fine structure increases with increasing $\mathit{Re}_{M}$. Its power spectrum $M(k)$ displays sharp peaks reflecting the fractal nature of $B(x)$ above the diffusive scale. We explore the saturation of these dynamos in three ways: via a renormalized reduced effective $\mathit{Re}_{M}$ (case I) or due to a decrease in the efficiency of the field amplification by stretching, without changing the map (case IIa), or changing the map (case IIb), and a combination of both effects (case III). For case I, we show that $B(x)$ in the saturated state, for both types of maps, approaches the marginal eigenfunction, which is obtained for $\mathit{Re}_{M}=R_{\mathit{crit}}$ independent of the initial $\mathit{Re}_{M}=R_{M0}$. On the other hand, in case II, for the two-strip map, we show that $B(x)$ saturates, preserving the structure of the kinematic eigenfunction. Thus the energy is transferred to larger scales in case I but remains at the smallest resistive scales in case II, as can be seen from both $B(x)$ and $M(k)$. For the four-strip map, $B(x)$ oscillates with time, although with a structure similar to the kinematic eigenfunction. Interestingly, the saturated state in case III shows an intermediate behaviour, with $B(x)$ similar to the kinematic eigenfunction at an intermediate $\mathit{Re}_{M}=R_{\mathit{sat}}$, with $R_{M0}>R_{\mathit{sat}}>R_{\mathit{crit}}$. The $R_{\mathit{sat}}$ value is determined by the relative importance of the increased diffusion versus the reduced stretching. These saturation properties are akin to the range of possibilities that have been discussed in the context of fluctuation dynamos.

We derive equations for the mean entropy and the mean internal energy in low-Mach-number temperature stratified turbulence (i.e. for turbulent convection or stably stratified turbulence), and show that turbulent flux of entropy is given by $\boldsymbol{F}_{s}=\overline{{\it\rho}}\,\overline{\boldsymbol{u}s}$, where $\overline{{\it\rho}}$ is the mean fluid density, $s$ is fluctuation of entropy and overbars denote averaging over an ensemble of turbulent velocity fields, $\boldsymbol{u}$. We demonstrate that the turbulent flux of entropy is different from the turbulent convective flux, $\boldsymbol{F}_{c}=\overline{T}\,\overline{{\it\rho}}\,\overline{\boldsymbol{u}s}$, of the fluid internal energy, where $\overline{T}$ is the mean fluid temperature. This turbulent convective flux is well-known in the astrophysical and geophysical literature, and it cannot be used as a turbulent flux in the equation for the mean entropy. This result is exact for low-Mach-number temperature stratified turbulence and is independent of the model used. We also derive equations for the velocity–entropy correlation, $\overline{\boldsymbol{u}s}$, in the limits of small and large Péclet numbers, using the quasi-linear approach and the spectral ${\it\tau}$ approximation, respectively. This study is important in view of different applications to astrophysical and geophysical temperature stratified turbulence.

Accretion disc theory is less developed than stellar evolution theory although a similarly mature phenomenological picture is ultimately desired. While the interplay of theory and numerical simulations has amplified community awareness of the role of magnetic fields in angular momentum transport, there remains a long term challenge to incorporate the insights gained from simulations into improving practical models for comparison with observations. What has been learned from simulations that can lead to improvements beyond SS73 in practical models? Here, we emphasize the need to incorporate the role of non-local transport more precisely. To show where large-scale transport would fit into the theoretical framework and how it is currently missing, we review why the wonderfully practical approach of Shakura & Sunyaev (Astron. Astrophys., vol. 24, 1973, pp. 337–355, SS73) is necessarily a mean field theory, and one which does not include large-scale transport. Observations of coronae and jets, combined with the interpretation of results from shearing box simulations, of the magnetorotational instability (MRI) suggest that a significant fraction of disc transport is indeed non-local. We show that the Maxwell stresses in saturation are dominated by large-scale contributions and that the physics of MRI transport is not fully captured by a viscosity. We also clarify the standard physical interpretation of the MRI as it applies to shearing boxes. Computational limitations have so far focused most attention toward local simulations, but the next generation of global simulations should help to inform improved mean field theories. Mean field accretion theory and mean field dynamo theory should in fact be unified into a single theory that predicts the time evolution of spectra and luminosity from separate disc, corona and outflow contributions. Finally, we note that any mean field theory, including that of SS73, has a finite predictive precision that needs to be quantified when comparing the predictions to observations.

Recent advances in dynamo theory have been made by examining the competition between small- and large-scale dynamos at high magnetic Reynolds number $\mathit{Rm}$. Small-scale dynamos rely on the presence of chaotic stretching whilst the generation of large-scale fields occurs in flows lacking reflectional symmetry via a systematic electromotive force (EMF). In this paper we discuss how the statistics of the EMF (at high $\mathit{Rm}$) depend on the properties of the multi-scale velocity that is generating it. In particular, we determine that different scales of flow have different contributions to the statistics of the EMF, with smaller scales contributing to the mean without increasing the variance. Moreover, we determine when scales in such a flow act independently in their contribution to the EMF. We further examine the role of large-scale shear in modifying the EMF. We conjecture that the distribution of the EMF, and not simply the mean, largely determines the dominant scale of the magnetic field generated by the flow.

Statistical characteristics of plasma fluctuations in the solar wind (SW), the Earth’s magnetosphere and fusion devices are reviewed. The turbulence in all these media has a complicated multiscale structure and exhibits a generalized self-similarity in an extended scale range. The anomalous transport of mass and momentum is intermittent and is carried by sporadic plasma flux bursts with non-Gaussian statistics, long-range correlation and multifractality. Intermittent turbulent transport is characterized by superdiffusion with power law $\langle {\it\delta}x^{2}\rangle \propto {\it\tau}^{{\it\alpha}}$, ${\it\alpha}\approx 1.2{-}1.8$. The structure functions in all these plasma environments are well fitted by the log-Poisson model of turbulence. Intermittent plasma turbulence displays universal properties and consists of quasi-1-D singular dissipative structures.

In previous galactic dynamo theories of the origin of the magnetic field in our galaxy, the subject of flux-freezing has been omitted. As a consequence, the equation of mass flow has generally also been omitted, particularly in the halo where the galactic gravitational field will operate on the mass flow. In this paper, it has been shown that this neglect could have serious consequences for the results obtained from those galactic dynamo simulations that include the halo. A modification of these dynamo theories is proposed which involves the expulsion of very small pieces of the magnetic field lines, rather than the wholesale expulsion of the complete magnetic lines encapsulated in the previous theories. This expulsion is accomplished by a spike instability that arises from superbubbles when they break out of the galactic disc and their shells fragment. This leads to a cut in the lines of force that still remain in the disc. Subsequently, normal disc turbulence rotates the cut lines and thus dissipates their mean flux, removing them from a role in the dynamo theory. This new process takes a length of time comparable to, but slightly longer than, the previous growth time of the disc dynamo, but avoids the previous difficulties associated with flux freezing and flux expulsion.