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.

Solar wind provides an example of a weakly collisional plasma expanding from a thermal source in the presence of spatially diverging magnetic field lines. Observations show that in the inner heliosphere, the electron temperature declines with the distance approximately as $T_{e}(r)\sim r^{-0.3} \dots r^{-0.7}$, which is significantly slower than the adiabatic expansion law $ \sim r^{-4/3}$. Motivated by such observations, we propose a kinetic theory that addresses the non-adiabatic evolution of a nearly collisionless plasma expanding from a central thermal source. We concentrate on the dynamics of energetic electrons propagating along a radially diverging magnetic flux tube. Due to the conservation of their magnetic moments, the electrons form a beam collimated along the magnetic field lines. Due to weak energy exchange with the background plasma, the beam population slowly loses its energy and heats the background plasma. We propose that no matter how weak the collisions are, at large enough distances from the source a universal regime of expansion is established where the electron temperature declines as $T_e(r)\propto r^{-2/5}$. This is close to the observed scaling of the solar wind temperature in the inner heliosphere. Our first-principle kinetic derivation may thus provide an explanation for the slower-than-adiabatic temperature decline in the solar wind. More broadly, it may be useful for describing magnetized winds from G-type stars.

Recently, it has been realized that magnetic plasma turbulence and magnetic field reconnection are inherently related phenomena. Turbulent fluctuations generate regions of a sheared magnetic field that become unstable to the tearing instability and reconnection, thus modifying turbulence at the corresponding scales. In this contribution, we give a brief discussion of some recent results on tearing-mediated magnetic turbulence. We illustrate the main ideas of this rapidly developing field of study by concentrating on two important examples -- magnetohydrodynamic Alfv\'en turbulence and small-scale kinetic-Alfv\'en turbulence. Due to various potential applications of these phenomena in space physics and astrophysics, we specifically try not to overload the text by heavy analytical derivations, but rather present a qualitative discussion accessible to a non-expert in the theories of turbulence and reconnection.

Recent analytical works on strong magnetized plasma turbulence have hypothesized the existence of a range of scales where the tearing instability may govern the energy cascade. In this paper, we estimate the conditions under which such tearing may give rise to full nonlinear magnetic reconnection in the turbulent eddies, thereby enabling significant energy conversion and dissipation. When those conditions are met, a new turbulence regime is accessed where reconnection-driven energy dissipation becomes common, rather than the rare feature that it must be when they are not.

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.

Recent studies have suggested that the tearing instability may play a significant role in magnetic turbulence. In this work, we review the theory of the magnetohydrodynamic tearing instability in the general case of an arbitrary tearing parameter, which is relevant for applications in turbulence. We discuss a detailed derivation of the results for the standard Harris profile and accompany it by the derivation of the results for a lesser known sine-shaped profile. We devote special attention to the exact solution of the inner equation, which is the central result in the theory of tearing instability. We also briefly discuss the influence of shear flows on tearing instability in magnetic structures. Our presentation is self-contained; we expect it to be accessible to researchers in plasma turbulence who are not experts in magnetic reconnection.

In this work, we analyze the kinetic stability of a solar wind electron distribution composed of core and strahl subpopulations. The core is modeled by a drifting Maxwellian distribution, while the strahl is modeled by an analytic function recently derived in (Horaites et al. 2018) from the collisional kinetic equation. We perform a numerical linear stability analysis using the LEOPARD solver (Astfalk & Jenko 2017), which allows for arbitrary gyrotropic distribution functions in a magnetized plasma. We find that for typical solar wind conditions, the core-strahl distribution is unstable to the kinetic Alfv\'en and magnetosonic modes. The maximum growth rates for these instabilities occur at wavenumbers $k d_i \lesssim 1$, at moderately oblique angles of propagation, thus providing a potential source of kinetic-scale turbulence. In contrast with previous reports, we however do not find evidence for a whistler instability directly associated with the electron strahl. This may be related to the more realistic shape of the electron strahl distribution function adopted in our work. We therefore suggest that the whistler modes often invoked to explain anomalous scattering of strahl particles could appear as a result of nonlinear mode coupling and turbulent cascade originating at scales $k d_i \lesssim 1$.

Alfv\'enic-type turbulence in strongly magnetized, low-beta pair plasmas is investigated. A coupled set of equations for the evolution of the magnetic and flow potentials are derived, covering both fluid and kinetic scales. In the fluid (MHD) range those equations are the same as for electron-ion plasmas, so turbulence at those scales is expected to be of the Alfv\'enic nature, exhibiting critical balance, dynamic alignment, and transition to a tearing mediated regime at small scales. The critical scale at which a transition to a tearing-mediated range occurs is derived, and the spectral slope in that range is predicted to be $k_\perp^{-8/3}$ (or $k_\perp^{-3}$, depending on details of the reconnecting configuration assumed). At scales below the electron (and positron) skin depth, it is argued that turbulence is dictated by a cascade of the inertial Alfv\'en wave, which we show to result in the magnetic energy spectrum $\propto k_\perp^{-11/3}$.

It has recently been proposed (Loureiro & Boldyrev 2017; Mallet et al. 2017) that the inertial interval in magnetohydrodynamic (MHD) turbulence is terminated at small scales not by a Kolmogorov-like dissipation region, but rather by a new sub-inertial interval mediated by tearing instability. However, many astrophysical plasmas are nearly collisionless so that the MHD approximation is not applicable to turbulence at small scales. In this Letter, we propose the extension of the theory of reconnection-mediated turbulence to plasmas which are so weakly collisional that the reconnection occurring in the turbulent eddies is caused by electron inertia rather than by resistivity. We find that the transition scale to reconnection-mediated turbulence depends on the plasma beta and on the assumptions of the plasma turbulence model. However, in all cases analyzed, the energy spectra in the reconnection-mediated interval range from $E(k)dk\propto k^{-8/3}dk$ to $E(k)dk\propto k^{-3}dk$.

We present a combined observational and theoretical analysis to investigate the nature of plasma turbulence at kinetic scales in the Earth's magnetosheath. In the first decade of the kinetic range, just below the ion gyroscale, the turbulence was found to be similar to that in the upstream solar wind: predominantly anisotropic, low-frequency and kinetic Alfv\'en in nature. A key difference, however, is that the magnetosheath ions are typically much hotter than the electrons, $T_\mathrm{i}\gg T_\mathrm{e}$, which, together with $\beta_\mathrm{i}\sim 1$, leads to a change in behaviour in the second decade, close to electron scales. The turbulence here is characterised by an increased magnetic compressibility, following a mode we term the inertial kinetic Alfv\'en wave, and a steeper spectrum of magnetic fluctuations, consistent with the prediction $E_B(k_\perp)\propto k_\perp^{-11/3}$ that we obtain from a set of nonlinear equations. This regime of plasma turbulence may also be relevant for other astrophysical environments with $T_\mathrm{i}\gg T_\mathrm{e}$, such as the solar corona, hot accretion flows, and regions downstream of collisionless shocks.

Magnetic dynamo action caused by the magnetorotational instability is studied in the shearing-box approximation with no imposed net magnetic flux. Consistent with recent studies, the dynamo action is found to be sensitive to the aspect ratio of the box: it is much easier to obtain in tall boxes (stretched in the direction normal to the disk plane) than in long boxes (stretched in the radial direction). Our direct numerical simulations indicate that the dynamo is possible in both cases, given a large enough magnetic Reynolds number. To explain the relatively larger effort required to obtain the dynamo action in a long box, we propose that the turbulent eddies caused by the instability most efficiently fold and mix the magnetic field lines in the radial direction. As a result, in the long box the scale of the generated strong azimuthal (stream-wise directed) magnetic field is always comparable to the scale of the turbulent eddies. In contrast, in the tall box the azimuthal magnetic flux spreads in the vertical direction over a distance exceeding the scale of the turbulent eddies. As a result, different vertical sections of the tall box are permeated by large-scale nonzero azimuthal magnetic fluxes, facilitating the instability. In agreement with this picture, the cases when the dynamo is efficient are characterized by a strong intermittency of the local azimuthal magnetic fluxes.

In magnetohydrodynamic (MHD) turbulence, the large-scale magnetic field sets a preferred local direction for the small-scale dynamics, altering the statistics of turbulence from the isotropic case. This happens even in the absence of a total magnetic flux, since MHD turbulence forms randomly oriented large-scale domains of strong magnetic field. It is therefore customary to study small-scale magnetic plasma turbulence by assuming a strong background magnetic field relative to the turbulent fluctuations. This is done, for example, in reduced models of plasmas, such as reduced MHD, reduced-dimension kinetic models, gyrokinetics, etc., which make theoretical calculations easier and numerical computations cheaper. Recently, however, it has become clear that the turbulent energy dissipation is concentrated in the regions of strong magnetic field variations. A significant fraction of the energy dissipation may be localized in very small volumes corresponding to the boundaries between strongly magnetized domains. In these regions the reduced models are not applicable. This has important implications for studies of particle heating and acceleration in magnetic plasma turbulence. The goal of this work is to systematically investigate the relationship between local magnetic field variations and magnetic energy dissipation, and to understand its implications for modeling energy dissipation in realistic turbulent plasmas.

The intermittent small-scale structure of turbulence governs energy dissipation in many astrophysical plasmas and is often believed to have universal properties for sufficiently large systems. In this work, we argue that small-scale turbulence in accretion disks is universal in the sense that it is insensitive to the magnetorotational instability (MRI) and background shear, and therefore indistinguishable from standard homogeneous magnetohydrodynamic (MHD) turbulence at small scales. We investigate the intermittency of current density, vorticity, and energy dissipation in numerical simulations of incompressible MHD turbulence driven by the MRI in a shearing box. We find that the simulations exhibit a similar degree of intermittency as in standard MHD turbulence. We perform a statistical analysis of intermittent dissipative structures and find that energy dissipation is concentrated in thin sheet-like structures that span a wide range of scales up to the box size. We show that these structures exhibit strikingly similar statistical properties to those in standard MHD turbulence. Additionally, the structures are oriented in the toroidal direction with a characteristic tilt of approximately 17.5 degrees, implying an effective guide field in that direction.

The current understanding of MHD turbulence envisions turbulent eddies which are anisotropic in all three directions. In the plane perpendicular to the local mean magnetic field, this implies that such eddies become current-sheet-like structures at small scales. We analyze the role of magnetic reconnection in these structures and conclude that reconnection becomes important at a scale $\lambda\sim L S_L^{-4/7}$, where $S_L$ is the outer-scale ($L$) Lundquist number and $\lambda$ is the smallest of the field-perpendicular eddy dimensions. This scale is larger than the scale set by the resistive diffusion of eddies, therefore implying a fundamentally different route to energy dissipation than that predicted by the Kolmogorov-like phenomenology. In particular, our analysis predicts the existence of the sub-inertial, reconnection interval of MHD turbulence, with the Fourier energy spectrum $E(k_\perp)\propto k_\perp^{-5/2}$, where $k_\perp$ is the wave number perpendicular to the local mean magnetic field. The same calculation is also performed for high (perpendicular) magnetic Prandtl number plasmas ($Pm$), where the reconnection scale is found to be $\lambda/L\sim S_L^{-4/7}Pm^{-2/7}$.

Turbulence is ubiquitous in plasmas, leading to rich dynamics characterized by irregularity, irreversibility, energy fluctuations across many scales, and energy transfer across many scales. Another fundamental and generic feature of turbulence, although sometimes overlooked, is the inhomogeneous dissipation of energy in space and in time. This is a consequence of intermittency, the scale-dependent inhomogeneity of dynamics caused by fluctuations in the turbulent cascade. Intermittency causes turbulent plasmas to self-organize into coherent dissipative structures, which may govern heating, diffusion, particle acceleration, and radiation emissions. In this paper, we present recent progress on understanding intermittency in incompressible magnetohydrodynamic turbulence with a strong guide field. We focus on the statistical analysis of intermittent dissipative structures, which occupy a small fraction of the volume but arguably account for the majority of energy dissipation. We show that, in our numerical simulations, intermittent structures in the current density, vorticity, and Els\"{a}sser vorticities all have nearly identical statistical properties. We propose phenomenological explanations for the scalings based on general considerations of Els\"{a}sser vorticity structures. Finally, we examine the broader implications of intermittency for astrophysical systems.

The local properties of turbulence driven by the magnetorotational instability (MRI) in rotating, shearing flows are studied in the framework of a shearing-box model. Based on numerical simulations, we propose that the MRI-driven turbulence comprises two components: the large-scale shear-aligned strong magnetic field and the small-scale fluctuations resembling magnetohydrodynamic (MHD) turbulence. The energy spectrum of the large-scale component is close to $k^{-2}$, whereas the spectrum of the small-scale component agrees with the spectrum of strong MHD turbulence $k^{-3/2}$. While the spectrum of the fluctuations is universal, the outer-scale characteristics of the turbulence are not; they depend on the parameters of the system, such as the net magnetic flux. However, there is remarkable universality among the allowed turbulent states -- their intensity $v_0$ and their outer scale $\lambda_0$ satisfy the balance condition $v_0/\lambda_0\sim \mathrm d\Omega/\mathrm d\ln r$, where $\mathrm d\Omega/\mathrm d\ln r$ is the local orbital shearing rate of the flow. Finally, we find no sustained dynamo action in the $\mathrm{Pm}=1$ zero net-flux case for Reynolds numbers as high as $45\,000$, casting doubts on the existence of an MRI dynamo in the $\mathrm{Pm}\leq 1$ regime.

In this work we study the process of energy dissipation triggered by a slow large scale motion of a magnetized conducting fluid. Our consideration is motivated by the problem of heating the solar corona, which is believed to be governed by fast reconnection events set off by the slow motion of magnetic field lines anchored in the photospheric plasma. To elucidate the physics governing the disruption of the imposed laminar motion and the energy transfer to small scales, we propose a simplified model where the large-scale motion of magnetic field lines is prescribed not at the footpoints but rather imposed volumetrically. As a result, the problem can be treated numerically with an efficient, highly-accurate spectral method, allowing us to use a resolution and statistical ensemble exceeding those of the previous work. We find that, even though the large-scale deformations are slow, they eventually lead to reconnection events that drive a turbulent state at smaller scales. The small-scale turbulence displays many of the universal features of field-guided MHD turbulence like a well developed inertial range spectrum. Based on these observations, we construct a phenomenological model that gives the scalings of the amplitude of the fluctuations and the energy dissipation rate as functions of the input parameters. We find a good agreement between the numerical results and the predictions of the model.

Recent observations reveal that magnetic turbulence in the nearly colisionless solar wind plasma extends to scales smaller than the plasma microscales, such as ion gyroradius and ion inertial length. Measured breaks in the spectra of magnetic and density fluctuations at high frequencies are thought to be related to the transition from large-scale hydromagnetic to small-scale kinetic turbulence. The scales of such transitions and the responsible physical mechanisms are not well understood however. In the present work we emphasize the crucial role of the plasma parameters in the transition to kinetic turbulence, such as the ion and electron plasma beta, the electron to ion temperature ratio, the degree of obliquity of turbulent fluctuations. We then propose an explanation for the spectral breaks reported in recent observations.

Energy dissipation is highly intermittent in turbulent plasmas, being localized in coherent structures such as current sheets. The statistical analysis of spatial dissipative structures is an effective approach to studying turbulence. In this paper, we generalize this methodology to investigate four-dimensional spatiotemporal structures, i.e., dissipative processes representing sets of interacting coherent structures, which correspond to flares in astrophysical systems. We develop methods for identifying and characterizing these processes, and then perform a statistical analysis of dissipative processes in numerical simulations of driven magnetohydrodynamic turbulence. We find that processes are often highly complex, long-lived, and weakly asymmetric in time. They exhibit robust power-law probability distributions and scaling relations, including a distribution of dissipated energy with power-law index near -1.75, indicating that intense dissipative events dominate the overall energy dissipation. We compare our results with the previously observed statistical properties of solar flares.

Energy dissipation in magnetohydrodynamic (MHD) turbulence is known to be highly intermittent in space, being concentrated in sheet-like coherent structures. Much less is known about intermittency in time, another fundamental aspect of turbulence which has great importance for observations of solar flares and other space/astrophysical phenomena. In this Letter, we investigate the temporal intermittency of energy dissipation in numerical simulations of MHD turbulence. We consider four-dimensional spatiotemporal structures, "flare events", responsible for a large fraction of the energy dissipation. We find that although the flare events are often highly complex, they exhibit robust power-law distributions and scaling relations. We find that the probability distribution of dissipated energy has a power law index close to -1.75, similar to observations of solar flares, indicating that intense dissipative events dominate the heating of the system. We also discuss the temporal asymmetry of flare events as a signature of the turbulent cascade.

The inertial-interval energy spectrum of strong magnetohydrodynamic (MHD) turbulence with a uniform background magnetic field was observed numerically to be close to $k^{-3/2}$ by a number of independent groups. A dissenting opinion has been voiced by Beresnyak, A. 2011, PRL, 106, 075001-. 2012, MNRAS, 422, 3495-. 2014, ApJ, 784, L20 that the spectral scaling is close to $k^{-5/3}$. The conclusions of these papers are however incorrect as they are based on numerical simulations that are drastically unresolved, so that the discrete numerical scheme does not approximate the physical solution at the scales where the measurements are performed. These results have been rebutted in our more detailed papers Perez, J. C., Mason, J., Boldyrev, S., & Cattaneo, F. 2012, PRX, 2, 041005-. 2014, ApJL, 793, L13; here, by popular demand, we present a brief and simple explanation of our major criticism of Beresnyak's work.

A variety of studies of magnetised plasma turbulence invoke theories for the advection of a passive scalar by turbulent fluctuations. Examples include modelling the electron density fluctuations in the interstellar medium, understanding the chemical composition of galaxy clusters and the intergalactic medium, and testing the prevailing phenomenological theories of magnetohydrodynamic turbulence. While passive scalar turbulence has been extensively studied in the hydrodynamic case, its counterpart in MHD turbulence is significantly less well understood. Herein we conduct a series of high-resolution direct numerical simulations of incompressible, field-guided, MHD turbulence in order to establish the fundamental properties of passive scalar evolution. We study the scalar anisotropy, establish the scaling relation analogous to Yaglom's law, and measure the intermittency of the passive scalar statistics. We also assess to what extent the pseudo Alfven fluctuations in strong MHD turbulence can be modelled as a passive scalar. The results suggest that the dynamics of a passive scalar in MHD turbulence is considerably more complicated than in the hydrodynamic case.

We investigate the intermittency of energy dissipation in magnetohydrodynamic (MHD) turbulence by identifying dissipative structures and measuring their characteristic scales. We find that the probability distribution of energy dissipation rates exhibits a power law tail with index very close to the critical value of -2.0, which indicates that structures of all intensities contribute equally to energy dissipation. We find that energy dissipation is uniformly spread among coherent structures with lengths and widths in the inertial range. At the same time, these structures have thicknesses deep within the dissipative regime. As the Reynolds number is increased, structures become thinner and more numerous, while the energy dissipation continues to occur mainly in large-scale coherent structures. This implies that in the limit of high Reynolds number, energy dissipation occurs in thin, tightly packed current sheets which nevertheless span a continuum of scales up to the system size, exhibiting features of both coherent structures and nanoflares previously conjectured as a coronal heating mechanism.

Magnetohydrodynamic (MHD) turbulence in the majority of natural systems, including the interstellar medium, the solar corona, and the solar wind, has Reynolds numbers far exceeding the Reynolds numbers achievable in numerical experiments. Much attention is therefore drawn to the universal scaling properties of small-scale fluctuations, which can be reliably measured in the simulations and then extrapolated to astrophysical scales. However, in contrast with hydrodynamic turbulence, where the universal structure of the inertial and dissipation intervals is described by the Kolmogorov self-similarity, the scaling for MHD turbulence cannot be established based solely on dimensional arguments due to the presence of an intrinsic velocity scale -- the Alfven velocity. In this Letter, we demonstrate that the Kolmogorov first self-similarity hypothesis cannot be formulated for MHD turbulence in the same way it is formulated for the hydrodynamic case. Besides profound consequences for the analytical consideration, this also imposes stringent conditions on numerical studies of MHD turbulence. In contrast with the hydrodynamic case, the discretization scale in numerical simulations of MHD turbulence should decrease faster than the dissipation scale, in order for the simulations to remain resolved as the Reynolds number increases.

The nature of subproton scale fluctuations in the solar wind is an open question, partly because two similar types of electromagnetic turbulence can occur: kinetic Alfven turbulence and whistler turbulence. These two possibilities, however, have one key qualitative difference: whistler turbulence, unlike kinetic Alfven turbulence, has negligible power in density fluctuations. In this Letter, we present new observational data, as well as analytical and numerical results, to investigate this difference. The results show, for the first time, that the fluctuations well below the proton scale are predominantly kinetic Alfven turbulence, and, if present at all, the whistler fluctuations make up only a small fraction of the total energy.

Observations reveal nearly power-law spectra of magnetic and density plasma fluctuations at subproton scales in the solar wind, which indicates the presence of a turbulent cascade. We discuss the three-field and two-field models for micro-scale plasma fluctuations, and then present the results of numerical simulations of a two-field model of kinetic-Alfven turbulence, which models plasma motion at sub-proton scales.

We develop a framework for studying the statistical properties of current sheets in numerical simulations of 3D magnetohydrodynamic (MHD) turbulence. We describe an algorithm that identifies current sheets in a simulation snapshot and then determines their geometrical properties (including length, width, and thickness) and intensities (peak current density and total energy dissipation rate). We then apply this procedure to simulations of reduced MHD turbulence and perform a statistical analysis on the obtained population of current sheets. We evaluate the role of reconnection by separately studying the populations of current sheets which contain magnetic X-points and those which do not. We find that the statistical properties of the two populations are different in general. We compare the scaling of these properties to phenomenological predictions obtained for the inertial range of MHD turbulence. Finally, we test whether the reconnecting current sheets are consistent with the Sweet-Parker model.

The statistical properties of magnetic discontinuities in the solar wind are investigated by measuring fluctuations in the magnetic field direction, given by the rotation Delta theta that the magnetic field vector undergoes during time interval Delta t. We show that the probability density function for rotations, P(Delta theta), can be described by a simple model in which the magnetic field vector rotates with a relative increment (Delta B)/B that is lognormally distributed. We find that the probability density function of increments, P((Delta B)/B), has a remarkable scaling property: the normalized variable x=[(Delta B)/B]*[(Delta t)/(Delta t_0)]^-a has a universal lognormal distribution for all time intervals Delta t. We then compare measurements from the solar wind with those from direct numerical simulations of magnetohydrodynamic (MHD) turbulence. We find good agreement for P(Delta theta) obtained in the two cases when the magnetic guide-field to fluctuations ratio B_0/b_rms is chosen accordingly. However, the scale invariance of P((Delta B)/B) is broken in the MHD simulations with relatively limited inertial interval, which causes P(Delta theta) to scale with measurement interval differently than in the solar wind.

The energy spectrum of magnetohydrodynamic turbulence attracts interest due to its fundamental importance and its relevance for interpreting astrophysical data. Here we present measurements of the energy spectra from a series of high-resolution direct numerical simulations of MHD turbulence with a strong guide field and for increasing Reynolds number. The presented simulations, with numerical resolutions up to 2048^3 mesh points and statistics accumulated over 30 to 150 eddy turnover times, constitute, to the best of our knowledge, the largest statistical sample of steady state MHD turbulence to date. We study both the balanced case, where the energies associated with Alfv\'en modes propagating in opposite directions along the guide field, E^+ and $E^-, are equal, and the imbalanced case where the energies are different. In the balanced case, we find that the energy spectrum converges to a power law with exponent -3/2 as the Reynolds number is increased, consistent with phenomenological models that include scale-dependent dynamic alignment. For the imbalanced case, with E^+>E^-, the simulations show that E^- ~ k_{\perp}^{-3/2} for all Reynolds numbers considered, while E^+ has a slightly steeper spectrum at small Re. As the Reynolds number increases, E^+ flattens. Since both E^+ and E^- are pinned at the dissipation scale and anchored at the driving scales, we postulate that at sufficiently high Re the spectra will become parallel in the inertial range and scale as E^+ ~ E^- ~ k_{\perp}^{-3/2}. Questions regarding the universality of the spectrum and the value of the "Kolmogorov constant" are discussed.

A model for strong kinetic Alfv\'en plasma turbulence at scales smaller than the ion gyroscale is proposed. It is argued that magnetic and density fluctuations are concentrated mostly at two-dimensional structures, which leads to their Fourier energy spectra $E(k_\perp)\propto k_{\perp}^{-8/3}$, where $k_\perp$ is the wave-vector component normal to the strong background magnetic field. The results are shown to be in good agreement with numerical simulations, and they can explain recent observations of magnetic and density fluctuations in the solar wind at sub-proton scales.

Recent measurements of solar wind turbulence report the presence of intermittent, exponentially distributed angular discontinuities in the magnetic field. In this Letter, we study whether such discontinuities can be produced by magnetohydrodynamic (MHD) turbulence. We detect the discontinuities by measuring the fluctuations of the magnetic field direction, Delta theta, across fixed spatial increments Delta x in direct numerical simulations of MHD turbulence with an imposed uniform guide field B_0. A large region of the probability density function (pdf) for Delta theta is found to follow an exponential decay, proportional to exp(-Delta theta/theta_*), with characteristic angle theta_* ~ (14 deg) (b_rms/B_0)^0.65 for a broad range of guide-field strengths. We find that discontinuities observed in the solar wind can be reproduced by MHD turbulence with reasonable ratios of b_rms/B_0. We also observe an excess of small angular discontinuities when Delta x becomes small, possibly indicating an increasing statistical significance of dissipation-scale structures. The structure of the pdf in this case closely resembles the two-population pdf seen in the solar wind. We thus propose that strong discontinuities are associated with inertial-range MHD turbulence, while weak discontinuities emerge from near-dissipation-range turbulence. In addition, we find that the structure functions of the magnetic field direction exhibit anomalous scaling exponents, which indicates the existence of intermittent structures.

Magnetised plasma turbulence pervades the universe and is likely to play an important role in a variety of astrophysical settings. Magnetohydrodynamics (MHD) provides the simplest theoretical framework in which phenomenological models for the turbulent dynamics can be built. Numerical simulations of MHD turbulence are widely used to guide and test the theoretical predictions; however, simulating MHD turbulence and accurately measuring its scaling properties is far from straightforward. Computational power limits the calculations to moderate Reynolds numbers and often simplifying assumptions are made in order that a wider range of scales can be accessed. After describing the theoretical predictions and the numerical approaches that are often employed in studying strong incompressible MHD turbulence, we present the findings of a series of high-resolution direct numerical simulations. We discuss the effects that insufficiencies in the computational approach can have on the solution and its physical interpretation.

Recent numerical and observational studies revealed that spectra of magnetic and velocity fluctuations in MHD turbulence have different scaling indexes. This intriguing feature has been recently explained in the case of weak MHD turbulence, that is, turbulence consisting of weakly interacting Alfven waves. However, astrophysical turbulence is strong in majority of cases. In the present work, we propose a unifying picture that allows one to address weak and strong MHD turbulence on the same footing. We argue that magnetic and kinetic energies are different in both weak and strong MHD turbulence. Their difference, the so-called residual energy, is spontaneously generated by turbulence, it has the Fourier spectrum E_r(k)=E_v(k)-E_b(k) \propto -f_w(k_||/k_perp) k_perp^{-2} in weak turbulence, and E_r(k) \propto -f_s(k_||/k_perp) k_perp^{-3} in strong turbulence. Here f_w,s(x) are functions declining fast for x>C_w,s and not significantly varying for $x<C_w,s$ with some constants C_w,s, and k_|| and k_perp the field-parallel and field-perpendicular wave vectors with respect to the applied strong uniform magnetic field.

Recent observations indicate that kinetic and magnetic energies are not in equipartition in the solar wind turbulence. Rather, magnetic fluctuations are more energetic and have somewhat steeper energy spectrum compared to the velocity fluctuations. This leads to the presence of the so-called residual energy E_r=E_v-E_b in the inertial interval of turbulence. This puzzling effect is addressed in the present paper in the framework of weak turbulence theory. Using a simple model of weakly colliding Alfv\'en waves, we demonstrate that the kinetic-magnetic equipartition indeed gets broken as a result of nonlinear interaction of Alfv\'en waves. We establish that magnetic energy is indeed generated more efficiently as a result of these interactions, which proposes an explanation for the solar wind observations.

There is mounting evidence in solar wind observations and in numerical simulations that kinetic and magnetic energies are not in equipartition in magnetohydrodynamic turbulence. The origin of their mismatch, the residual energy E_r=E_v-E_b, is not well understood. In the present work this effect is studied analytically in the regime of weak magnetohydrodynamic turbulence. We find that residual energy is spontaneously generated by turbulent dynamics, and it has a negative sign, in good agreement with the observations. We obtain that the residual energy condenses around k_||=0 with its k_||-spectrum broadening linearly with k_perp, where k_|| and k_perp are the wavenumbers parallel and perpendicular to the background magnetic field, and the field-perpendicular spectrum of the residual energy has the scaling E_r(k_perp)\propto k_perp^{-1} in the inertial interval. These results are found to be in agreement with numerical simulations. We propose that residual energy plays a fundamental role in Alfvenic turbulence and it should be taken into account for correct interpretation of observational and numerical data.

The question is addressed to what extent incompressible magnetohydrodynamics (MHD) can describe random magnetic and velocity fluctuations measured in the solar wind. It is demonstrated that distributions of spectral indices for the velocity, magnetic field, and total energy obtained from high resolution numerical simulations are qualitatively and quantitatively similar to solar wind observations at 1 AU. Both simulations and observations show that in the inertial range the magnetic field spectrum E_b is steeper than the velocity spectrum E_v with E_b >~ E_v and that the residual energy E_R = E_b-E_v decreases nearly following a k_perp^-2 scaling.

Magnetised turbulence is ubiquitous in astrophysical systems, where it notoriously spans a broad range of spatial scales. Phenomenological theories of MHD turbulence describe the self-similar dynamics of turbulent fluctuations in the inertial range of scales. Numerical simulations serve to guide and test these theories. However, the computational power that is currently available restricts the simulations to Reynolds numbers that are significantly smaller than those in astrophysical settings. In order to increase computational efficiency and, therefore, probe a larger range of scales, one often takes into account the fundamental anisotropy of field-guided MHD turbulence, with gradients being much slower in the field-parallel direction. The simulations are then optimised by employing the reduced MHD equations and relaxing the field-parallel numerical resolution. In this work we explore a different possibility. We propose that there exist certain quantities that are remarkably stable with respect to the Reynolds number. As an illustration, we study the alignment angle between the magnetic and velocity fluctuations in MHD turbulence, measured as the ratio of two specially constructed structure functions. We find that the scaling of this ratio can be extended surprisingly well into the regime of relatively low Reynolds number. However, the extended scaling becomes easily spoiled when the dissipation range in the simulations is under-resolved. Thus, taking the numerical optimisation methods too far can lead to spurious numerical effects and erroneous representation of the physics of MHD turbulence, which in turn can affect our ability to correctly identify the physical mechanisms that are operating astrophysical systems.

A review of MHD dynamos and turbulence.

Hydromagnetic dynamo theory provides the prevailing theoretical description for the origin of magnetic fields in the universe. Here we consider the problem of kinematic, small-scale dynamo action driven by a random, incompressible, non-helical, homogeneous and isotropic flow. In the Kazantsev dynamo model the statistics of the driving flow are assumed to be instantaneously correlated in time. Here we compare the results of the model with the dynamo properties of a simulated flow that has equivalent spatial characteristics as the Kazantsev flow but different temporal statistics. In particular, the simulated flow is a solution of the forced Navier-Stokes equations and hence has a finite correlation time. We find that the Kazantsev model typically predicts a larger magnetic growth rate and a magnetic spectrum that peaks at smaller scales. However, we show that by filtering the diffusivity spectrum at small scales it is possible to bring the growth rates into agreement and simultaneously align the magnetic spectra.

Amplification of magnetic field due to kinematic turbulent dynamo action is studied in the regime of small magnetic Prandtl numbers. Such a regime is relevant for planets and stars interiors, as well as for liquid metal laboratory experiments. A comprehensive analysis based on the Kazantsev-Kraichnan model is reported, which establishes the dynamo threshold and the dynamo growth rates for varying kinetic helicity of turbulent fluctuations. It is proposed that in contrast with the case of large magnetic Prandtl numbers, the kinematic dynamo action at small magnetic Prandtl numbers is significantly affected by kinetic helicity, and it can be made quite efficient with an appropriate choice of the helicity spectrum.

Magnetohydrodynamics (MHD) provides the simplest description of magnetic plasma turbulence in a variety of astrophysical and laboratory systems. MHD turbulence with nonzero cross helicity is often called imbalanced, as it implies that the energies of Alfv\'en fluctuations propagating parallel and anti-parallel the background field are not equal. Recent analytical and numerical studies have revealed that at every scale, MHD turbulence consists of regions of positive and negative cross helicity, indicating that such turbulence is inherently locally imbalanced. In this paper, results from high resolution numerical simulations of steady-state incompressible MHD turbulence, with and without cross helicity are presented. It is argued that the inertial range scaling of the energy spectra (E^+ and E^-) of fluctuations moving in opposite directions is independent of the amount of cross-helicity. When cross helicity is nonzero, E^+ and E^- maintain the same scaling, but have differing amplitudes depending on the amount of cross-helicity.

Magnetohydrodynamics (MHD) is invoked to address turbulent fluctuations in a variety of astrophysical systems. MHD turbulence in nature is often anisotropic and imbalanced, in that Alfvenic fluctuations moving in opposite directions along the background magnetic field carry unequal energies. This work formulates specific requirements for effective numerical simulations of strong imbalanced MHD turbulence with a guide field B0 High-resolution simulations are then performed and they suggest that the spectra of the counter-propagating Alfven modes do not differ from the balanced case, while their amplitudes and the corresponding rates of energy cascades are significantly affected by the imbalance. It is further proposed that the stronger the imbalance the larger the magnetic Reynolds number that is required in numerical simulations in order to correctly reproduce the turbulence spectrum. This may explain current discrepancies among numerical simulations and observations of imbalanced MHD turbulence.

Turbulence of magnetohydrodynamic waves in nature and in the laboratory is generally cross-helical or non-balanced, in that the energies of Afv\'en waves moving in opposite directions along the guide magnetic field are unequal. Based on high-resolution numerical simulations it is proposed that such turbulence spontaneously generates a condensate of the residual energy $E_v-E_b$ at small field-parallel wave numbers. As a result, the energy spectra of counter-propagating Alfv\'en waves are not scale-invariant. In the limit of infinite Reynolds number, the universality is asymptotically restored at large wave numbers, and both spectra attain the scaling $E(k)\propto k_{\perp}^{-2}$. The generation of condensate is apparently related to the breakdown of mirror symmetry in non-balanced turbulence.

Magnetohydrodynamic (MHD) turbulence is pervasive in astrophysical systems. Recent high-resolution numerical simulations suggest that the energy spectrum of strong incompressible MHD turbulence is $E(k_{\perp})\propto k_{\perp}^{-3/2}$. So far, there has been no phenomenological theory that simultaneously explains this spectrum and satisfies the exact analytic relations for MHD turbulence due to Politano & Pouquet. Indeed, the Politano-Pouquet relations are often invoked to suggest that the spectrum of MHD turbulence instead has the Kolmogorov scaling -5/3. Using geometrical arguments and numerical tests, here we analyze this seeming contradiction and demonstrate that the -3/2 scaling and the Politano-Pouquet relations are reconciled by the phenomenon of scale-dependent dynamic alignment that was recently discovered in MHD turbulence.

We investigate the structure of magnetic field amplified by turbulent velocity fluctuations, in the framework of the kinematic Kazantsev-Kraichnan model. We consider Kolmogorov distribution of velocity fluctuations, and assume that both Reynolds number and magnetic Reynolds number are very large. We present the full numerical solution of the model for the spectra and the growth rates of magnetic fluctuations. We consider astrophysically relevant limits of large and small magnetic Prandtl numbers, and address both helical and nonhelical cases.

A combination of observation, theory, modeling, and laboratory plasma experiments provides a multifaceted approach to develop a much greater understanding of how magnetic fields arise in galactic settings and how these magnetic fields mediate important processes that affect the dynamics, distribution, and composition of galactic plasmas. An important emphasis below is the opportunity to connect laboratory experiments to astrophysics. This approach is especially compelling for the galactic neighborhood, where the distribution and character of magnetic fields can be observed with greater detail than what is possible elsewhere in the universe. The ability to produce laboratory plasmas with unparalleled accessibility permits an even greater level of detail to be assessed and exposed. Theory and modeling provide fundamental ways to understand important processes, and they act as the bridge to connect experimental validation to astronomical observations. In many cases the studies that utilize this approach can make use of existing laboratory facilities, resulting in a cost that is quite small compared to the cost of measurements in dedicated space missions.

Strong incompressible three-dimensional magnetohydrodynamic turbulence is investigated by means of high resolution direct numerical simulations. The simulations show that the configuration space is characterized by regions of positive and negative cross-helicity, corresponding to highly aligned or anti-aligned velocity and magnetic field fluctuations, even when the average cross-helicity is zero. To elucidate the role of cross-helicity, the spectra and structure of turbulence are obtained in imbalanced regions where cross-helicity is non-zero. When averaged over regions of positive and negative cross-helicity, the result is consistent with the simulations of balanced turbulence. An analytical explanation for the obtained results is proposed.

We report experimental observation of an instability in a Couette-Taylor flow of a polymer fluid in a thin gap between two coaxially rotating cylinders in a regime where their angular velocity decreases with the radius while the specific angular momentum increases with the radius. In the considered regime, neither the inertial Rayleigh instability nor the purely elastic instability are possible. We propose that the observed "elasto-rotational" instability is an analog of the magnetorotational instability which plays a fundamental role in astrophysical Keplerian accretion disks.

We investigate the growth and structure of magnetic fields amplified by kinematic dynamo action in turbulence with non-zero kinetic helicity. We assume a simple Gaussian velocity correlation tensor, which allows us to consider very large magnetic Reynolds numbers, up to one trillion. We use the kinematic Kazantsev-Kraichnan model of dynamo and find a complete numerical solution for the correlation functions of growing magnetic fields.

We investigate magnetic field amplification in a turbulent velocity field with nonzero helicity, in the framework of the kinematic Kazantsev-Kraichnan model. We present the numerical solution of the model for the practically important case of Kolmogorov distribution of velocity fluctuations, with a large magnetic Reynolds number. We find that in contrast to the nonhelical case where growing magnetic fields are described by a few bound eigenmodes concentrated inside the inertial interval of the velocity field, in the helical case the number of bound eigenmodes considerably increases; moreover, new unbound eigenmodes appear. Both bound and unbound eigenmodes contribute to the large-scale magnetic field. This indicates a limited applicability of the conventional alpha model of a large-scale dynamo action, which captures only unbound modes.

Recent numerical and observational studies contain conflicting reports on the spectrum of magnetohydrodynamic turbulence. In an attempt to clarify the issue we investigate anisotropic incompressible magnetohydrodynamic turbulence with a strong guide field $B_0$. We perform numerical simulations of the reduced MHD equations in a special setting that allows us to elucidate the transition between weak and strong turbulent regimes. Denote $k_{\|}$, $k_\perp$ characteristic field-parallel and field-perpendicular wavenumbers of the fluctuations, and $b_{\lambda}$ the fluctuating field at the scale $\lambda\sim 1/k_{\perp}$. We find that when the critical balance condition, $k_{\|}B_0\sim k_{\perp} b_{\lambda}$, is satisfied, the turbulence is strong, and the energy spectrum is $E(k_{\perp})\propto k^{-3/2}_{\perp}$. As the $k_{\|}$ width of the spectrum increases, the turbulence rapidly becomes weaker, and in the limit $k_{\|}B_0\gg k_{\perp} b_{\lambda}$, the spectrum approaches $E(k_{\perp})\propto k_{\perp}^{-2}$. The observed sensitivity of the spectrum to the balance of linear and nonlinear interactions may explain the conflicting numerical and observational findings where this balance condition is not well controlled.

We report the results of an extensive set of direct numerical simulations of forced, incompressible, magnetohydrodynamic turbulence with a strong guide field. The aim is to resolve the controversy regarding the power law exponent (\alpha, say) of the field perpendicular energy spectrum E(k_\perp) \propto k_\perp ^ {\alpha}. The two main theoretical predictions, \alpha=-3/2 and \alpha=-5/3, have both received some support from numerical simulations carried out by different groups. Our simulations have a resolution of 512^3 mesh points, a strong guide field, an anisotropic simulation domain, and implement a broad range of large-scale forcing routines, including those previously reported in the literature. Our findings indicate that the spectrum of well developed, strong incompressible MHD turbulence with a strong guide field is E(k_{\perp})\propto k_{\perp}^{-3/2}.

Motivated by recent analytic predictions, we report numerical evidence showing that in driven incompressible magnetohydrodynamic turbulence the magnetic- and velocity-field fluctuations locally tend to align the directions of their polarizations. This dynamic alignment is stronger at smaller scales with the angular mismatch between the polarizations decreasing with the scale \lambda approximately as \theta_\lambda ~ \lambda^{1/4}. This can naturally lead to a weakening of the nonlinear interactions and provide an explanation for the energy spectrum E(k) ~ k^{-3/2} that is observed in numerical experiments of strongly magnetized turbulence.

We sought to determine the scale of scintillation in the interstellar plasma of PSR J0437-4715. We used the Very Long Baseline Array to obtain scintillation amplitude and phase data, from dynamic spectra at 327 MHz. We observe two scales of scintillation of pulsar PSR J0437-4715, differing by more than an order of magnitude in scintillation bandwidth. The wider-bandwidth scale of scintillation that we observe indicates less scattering for this pulsar than for other nearby pulsars, except for PSR B0950+08.

We propose a phenomenological theory of strong incompressible magnetohydrodynamic turbulence in the presence of a strong large-scale external magnetic field. We argue that in the inertial range of scales, magnetic-field and velocity-field fluctuations tend to align the directions of their polarizations. However, the perfect alignment cannot be reached, it is precluded by the presence of a constant energy flux over scales. As a consequence, the directions of fluid and magnetic-field fluctuations at each scale $\lambda$ become effectively aligned within the angle $\phi_{\lambda}\propto \lambda^{1/4}$, which leads to scale-dependent depletion of nonlinear interaction and to the field-perpendicular energy spectrum $E(k_{\perp})\propto k_{\perp}^{-3/2}$. Our results may be universal, i.e., independent of the external magnetic field, since small-scale fluctuations locally experience a strong field produced by large-scale eddies.

It was recently suggested by Boldyrev & Gwinn that the characteristics of radio scintillations from distant pulsars are best understood if the interstellar electron-density fluctuations that cause the time broadening of the radio pulses obey non-Gaussian statistics. In this picture the density fluctuations are inferred to be strong on very small scales ($\sim 10^8-10^{10} {cm}$). We argue that such density structures could correspond to the ionized boundaries of molecular regions (clouds) and demonstrate that the power-law distribution of scattering angles that is required to match the observations arises naturally from the expected intersections of our line of sight with randomly distributed, thin, approximately spherical ionized shells of this type. We show that the observed change in the time-broadening behavior for pulsar dispersion measures $\lesssim 30 {\rm pc} {\rm cm}^{-3}$ is consistent with the expected effect of the general ISM turbulence, which should dominate the scattering for nearby pulsars. We also point out that if the clouds are ionized by nearby stars, then their boundaries may become turbulent on account of an ionization front instability. This turbulence could be an alternative cause of the inferred density structures. An additional effect that might contribute to the strength of the small-scale fluctuations in this case is the expected flattening of the turbulent density spectrum when the eddy sizes approach the proton gyroscale.

A great deal of study has been carried out over the last twenty years on the origin of the magnetic activity in the Galactic center. One of the most popular hypotheses assumes milli-Gauss magnetic field with poloidal geometry, pervading the inner few hundred parsecs of the Galactic-center region. However, there is a growing observational evidence for the large-scale distribution of a much weaker field of B \lesssim 10 micro G in this region. Here, we propose that the Galactic-center magnetic field originates from turbulent activity that is known to be extreme in the central hundred parsecs. In this picture the spatial distribution of the magnetic field energy is highly intermittent, and the regions of strong field have filamentary structures. We propose that the observed nonthermal radio filaments appear in (or, possibly, may be identified with) such strongly magnetized regions. At the same time, the large-scale diffuse magnetic field is weak. Both results of our model can explain the magnetic field measurements of the the Galactic-center region. In addition, we discuss the role of ionized outflow from stellar clusters in producing the long magnetized filaments perpendicular to the Galactic plane.

We investigate analytically the amplification of a weak magnetic field in a homogeneous and isotropic turbulent flow lacking reflectional symmetry (helical turbulence). We propose that the spectral distributions of magnetic energy and magnetic helicity can be found as eigenmodes of a self-adjoint, Schr\"odinger-type system of evolution equations. We argue that large-scale and small-scale magnetic fluctuations cannot be effectively separated, and that the conventional alpha-model is, in general, not an adequate description of the large-scale dynamo mechanism. As a consequence, the correct numerical modeling of such processes should resolve magnetic fluctuations down to the very small, resistive scales.

Radio waves propagating from distant pulsars in the interstellar medium (ISM), are refracted by electron density inhomogeneities, so that the intensity of observed pulses fluctuates with time. The theory relating the observed pulse time-shapes to the electron-density correlation function has developed for 30 years, however, two puzzles have remained. First, observational scaling of pulse broadening with the pulsar distance is anomalously strong; it is consistent with the standard model only when non-uniform statistics of electron fluctuations along the line of sight are assumed. Second, the observed pulse shapes are consistent with the standard model only when the scattering material is concentrated in a narrow slab between the pulsar and the Earth. We propose that both paradoxes are resolved at once if one assumes stationary and uniform, but non-Gaussian statistics of the electron-density distribution. Such statistics must be of Levy type, and the propagating ray should exhibit a Levy flight. We propose that a natural realization of such statistics may be provided by the interstellar medium with random electron-density discontinuities. We develop a theory of wave propagation in such a non-Gaussian random medium, and demonstrate its good agreement with observations. The qualitative introduction of the approach and the resolution of the anomalous-scaling paradox was presented earlier in [PRL 91, 131101 (2003); ApJ 584, 791 (2003)].

We propose a phenomenological model for incompressible magnetohydrodynamic turbulence. We argue that nonlinear-wave interaction weakens as the energy cascade proceeds to small scales, however, the anisotropy of fluctuations along the large-scale magnetic field increases, which makes turbulence strong at all scales. To explain the weakening of the interaction, we propose that small-scale fluctuations of the velocity and magnetic fields become increasingly dynamically aligned as their scale decreases, so that turbulent eddies become locally anisotropic in the plane perpendicular to the large-scale magnetic field. In the limit of weak anisotropy, that is, weak large-scale magnetic field, our model reproduces the Goldreich-Sridhar spectrum, while the limit of strong anisotropy, that is, strong large-scale magnetic field, corresponds to the Iroshnikov-Kraichnan scaling of the spectrum. This is in good agreement with recent numerical results.

We consider the one-dimensional Burgers equation randomly stirred at large scales by a Gaussian short-time correlated force. Using the method of dissipative anomalies, we obtain velocity and velocity-difference probability density functions and confirm the results with high-resolution numerical simulations.

Supersonic turbulent flows of magnetized gas are believed to play an important role in the dynamics of star-forming clouds in galaxies. Understanding statistical properties of such flows is crucial for developing a theory of star formation. In this letter we propose a unified approach for obtaining the velocity scaling in compressible and super--Alfv\'{e}nic turbulence, valid for arbitrary sonic Mach number, \ms. We demonstrate with numerical simulations that the scaling can be described with the She--L\'{e}v\^{e}que formalism, where only one parameter, interpreted as the Hausdorff dimension of the most intense dissipative structures, needs to be varied as a function of \ms. Our results thus provide a method for obtaining the velocity scaling in interstellar clouds once their Mach numbers have been inferred from observations.

We analyze the initial, kinematic stage of magnetic field evolution in an isotropic and homogeneous turbulent conducting fluid with a rough velocity field, v(l) ~ l^alpha, alpha<1. We propose that in the limit of small magnetic Prandtl number, i.e. when ohmic resistivity is much larger than viscosity, the smaller the roughness exponent, alpha, the larger the magnetic Reynolds number that is needed to excite magnetic fluctuations. This implies that numerical or experimental investigations of magnetohydrodynamic turbulence with small Prandtl numbers need to achieve extremely high resolution in order to describe magnetic phenomena adequately.

Observations of radio signals from distant pulsars provide a valuable tool for investigation of interstellar turbulence. The time-shapes of the signals are the result of pulse broadening by the fluctuating electron density in the interstellar medium. While the scaling of the shapes with the signal frequency is well understood, the observed anomalous scaling with respect to the pulsar distance has remained a puzzle for more than 30 years. We propose a new model for interstellar electron density fluctuations, which explains the observed scaling relations. We suggest that these fluctuations obey Levy statistics rather than Gaussian statistics, as assumed in previous treatments of interstellar scintillations.

We develop a model of thin turbulent accretion discs supported by magnetic pressure of turbulent magnetic fields. This applies when the turbulent kinetic and magnetic energy densities are greater than the thermal energy density in the disc. Whether such discs survive in nature or not remains to be determined, but here we simply demonstrate that self-consistent solutions exist when the alpha-prescription for the viscous stress, similar to that of the original Shakura-Sunyaev model, is used. We show that \alpha \sim 1 for the strongly magnetized case and we calculate the radial structure and emission spectra from the disc in the regime when it is optically thick. Strongly magnetized optically thick discs can apply to the full range of disc radii for objects < 10^{-2} of the Eddington luminosity or for the outer parts of discs in higher luminosity sources. In the limit that the magnetic pressure is equal to the thermal or radiation pressure, our strongly magnetized disc model transforms into the Shakura-Sunyaev model with \alpha=1. Our model produces spectra quite similar to those of standard Shakura-Sunyaev models. In our comparative study, we also discovered a small discrepancy in the spectral calculations of Shakura and Sunyaev (1973).

We compute the structure function scaling of the integrated intensity images of two J=1-0 13CO maps of Taurus and Perseus. The scaling exponents of the structure functions follow the velocity scaling of supersonic turbulence, suggesting that turbulence plays an important role in the fragmentation of cold interstellar clouds. The data also allows to verify the validity of the two basic assumptions of the hierarchical symmetry model, originally proposed for the derivation of the velocity structure function scaling. This shows that the same hierarchical symmetry holds for the projected density field of cold interstellar clouds.

Temporal broadening of pulsar signals results from electron density fluctuations in the interstellar medium that cause the radiation to travel along paths of different lengths. The Gaussian theory of fluctuations predicts that the pulse temporal broadening should scale with the wavelength as lambda^4, and with the dispersion measure (corresponding to distance to the pulsar) as DM^2. For large dispersion measure, DM > 20 pc/cm^3, the observed scaling is lambda^4 DM^4, contradicting the conventional theory. Although the problem has existed for 30 years, there has been no resolution to this paradox. We suggest that scintillations for distant pulsars are caused by non-Gaussian, spatially intermittent density fluctuations with a power-like probability distribution. This probability distribution does not have a second moment in a large range of density fluctuations, and therefore the previously applied conventional Fokker-Planck theory does not hold. Instead, we propose to apply the theory of Levy distributions (so-called Levy flights). Using the scaling analysis (confirmed by numerical simulations of ray propagation) we show that the observed scaling is recovered for large DM, if the density differences, delta N, have Levy distribution decaying as |delta N|^{-5/3}.

The interstellar medium (ISM) provides a unique laboratory for highly supersonic, driven hydrodynamics turbulence. We present a theory of such turbulence, confirm it by numerical simulations, and use the results to explain observational properties of interstellar molecular clouds, the regions where stars are born.

The existence of a weak galactic magnetic field has been repeatedly confirmed by observational data. The origin of this field has not as yet been explained in a fully satisfactory way and represents one of the main challenges of the astrophysical dynamo theory. In both the galactic dynamo theory and the primordial-origin theory, a major influence is exerted by the small-scale magnetic fluctuations. This article is devoted to constructing a systematic second-order statistical theory of such small-scale fields. The statistics of these fields are studied in the kinematic approximation and for the case of large Prandtl numbers, which is relevant for the galactic and protogalactic plasma. The advecting velocity field is assumed to be Gaussian and short-time correlated. Theoretical understanding of this kinematic dynamo model is a necessary prerequisite for any prospective nonlinear dynamo theory. The theory is developed for an arbitrary degree of compressibility and formally in d dimensions, which generalizes the previously known results, elicits the structure of the solutions, and uncovers a number of new effects. The magnetic energy spectra are studied as they grow and spread over scales during the initial stage of the field amplification. Exact Green's-function solutions are obtained. The spectral theory is supplemented by the study of magnetic-field correlation functions in the configuration space, where the dynamo problem can be mapped onto a particular one-dimensional quantum-mechanical problem. The latter approach is most suitable for the description of the kinematic dynamo in the long-time limit, i.e. when the magnetic excitation has spread over all scales present in the system. A simple way of calculating the growth rates of the magnetic fields in this long-time limit is proposed.

We present a direct numerical and analytical study of driven supersonic MHD turbulence that is believed to govern the dynamics of star-forming molecular clouds. We describe statistical properties of the turbulence by measuring the velocity difference structure functions up to the fifth order. In particular, the velocity power spectrum in the inertial range is found to be close to E(k) \~ k^{-1.74}, and the velocity difference scales as <|\Delta u|>~ L^{0.42}. The results agree well with the Kolmogorov--Burgers analytical model suggested for supersonic turbulence in [astro-ph/0108300]. We then generalize the model to more realistic, fractal structure of molecular clouds, and show that depending on the fractal dimension of a given molecular cloud, the theoretical value for the velocity spectrum spans the interval [-1.74 ... -1.89], while the corresponding window for the velocity difference scaling exponent is [0.42 ... 0.78].

The process of star formation in interstellar molecular clouds is believed to be controlled by driven supersonic magnetohydrodynamic turbulence. We suggest that in the inertial range such turbulence obeys the Kolmogorov law, while in the dissipative range it behaves as Burgers turbulence developing shock singularities. On the base of the She-Leveque analytical model we then predict the velocity power spectrum in the inertial range to be E_k ~ k^{-1.74}. This result reproduces the observational Larson law, <u^2_l>~ l^{0.74...0.76}, [Larson, MNRAS 194 (1981) 809] and agrees well with recent numerical findings by Padoan and Nordlund [astro-ph/0011465]. The application of the model to more general dissipative structures, with higher fractal dimensionality, leads to better agreement with recent observational results.

We formulate a solvable model that describes generation and saturation of mean magnetic field in a dynamo with kinetic helicity, in the limit of large magnetic Prandtl number. This model is based on the assumption that the stochastic part of the velocity field is Gaussian and white in time (the Kazantsev-Kraichnan ensemble), while the regular part describing the back reaction of the magnetic field is chosen from balancing the viscous and Lorentz stresses in the MHD Navier-Stokes equation. The model provides an analytical explanation for previously obtained numerical results.

The Batchelor passive advection is an advection by a smooth velocity field. If the velocity field is a delta-correlated in time random Gaussian process, then the problem is reduced to quantum mechanics of fluctuating velocity gradient u'(t). For the finite-time correlated velocity field, such a reduction does not exist. To illustrate this point, the second moment of a passively advected magnetic field is considered, and the stochastic calculus is used to find finite-time corrections to its growth rate. The growth rate depends on large scale properties of the velocity field. Moreover, the problem is not universal with respect to the short-time regularization: different regularizations give different answers for the growth rate.

We study geometric properties of a random Gaussian short-time correlated velocity field by considering statistics of a passively advected metric tensor. That describes universal properties of fluctuations of tensor objects frozen into the fluid and passively advected by it. The problem of one-point statistics of co- and contravariant tensors is solved exactly, provided the advected fields do not reach dissipative scales, which would break the symmetry of the problem. Asymptotic in time duality of the problem is established, which in the three-dimensional case relates the probabilities of the volume deformations into "tubes" and into "sheets".

The randomly driven Navier-Stokes equation without pressure in d-dimensional space is considered as a model of strong turbulence in a compressible fluid. We derive a closed equation for the velocity-gradient probability density function. We find the asymptotics of this function for the case of the gradient velocity field (Burgers turbulence), and provide a numerical solution for the two-dimensional case. Application of these results to the velocity-difference probability density function is discussed.