The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle--particle collisions allows departures from local thermodynamic equilibrium. These departures exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present a numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart. If instead the parallel viscous stress is left unabated -- a situation relevant to recent kinetic simulations of the fluctuation dynamo and to the early stages of the dynamo in a magnetized ICM -- the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics of the saturated state of the large-Prandtl-number (${\rm Pm}\gtrsim{1}$) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches magnetic-energy spectra. A prediction of this model, confirmed by our simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. In all cases that do exhibit a dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

The fundamental difference between incompressible ideal magnetohydrodynamics and the dynamics of a non-conducting fluid is that magnetic fields exert a tension force that opposes their bending; magnetic fields behave like elastic strings threading the fluid. It is natural, therefore, to expect that a magnetic field tangled at small length scales should resist a large-scale shear in an elastic way, much as a ball of tangled elastic strings responds elastically to an impulse. Furthermore, a tangled field should support the propagation of `magnetoelastic waves', the isotropic analogue of Alfv\'en waves on a straight magnetic field. Here, we study magnetoelasticity in the idealised context of an equilibrium tangled field configuration. In contrast to previous treatments, we explicitly account for intermittency of the Maxwell stress, and show that this intermittency necessarily decreases the frequency of magnetoelastic waves in a stable field configuration. We develop a mean-field formalism to describe magnetoelastic behaviour, retaining leading-order corrections due to the coupling of large- and small-scale motions, and solve the initial-value problem for viscous fluids subjected to a large-scale shear, showing that the development of small-scale motions results in anomalous viscous damping of large-scale waves. Finally, we test these analytic predictions using numerical simulations of standing waves on tangled, linear force-free magnetic-field equilibria.

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $\delta n/n_0 \gtrsim 2/\beta$ (where $\beta\gg{1}$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves' otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more `fluid-like' than kinetic. A nonlinear fluctuation-dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels, and weakly collisional thermodynamics when IAWs with $\delta n/n_0 \gtrsim 2/\beta$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio $Q_{\text{i}}/Q_{\text{e}}$ as a function of the compressive-to-Alfv\'enic driving power ratio $P_{\text{compr}}/P_{\text{AW}}$, of the ratio of ion thermal pressure to magnetic pressure $\beta_{\text{i}}$, and of the ratio of ion-to-electron background temperatures $T_{\text{i}}/T_{\text{e}}$. It is shown that $Q_{\text{i}}/Q_{\text{e}}$ is an increasing function of $P_{\text{compr}}/P_{\text{AW}}$. When the compressive driving is sufficiently large, $Q_{\text{i}}/Q_{\text{e}}$ approaches $\simeq P_{\text{compr}}/P_{\text{AW}}$. This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfv\'en waves is absent for both low and high $\beta_{\text{i}}$, meaning that the compressive driving is directly connected to the ion entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low $\beta_{\text{i}}$ and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.

A multi-institutional, multi-national science team will soon submit a NASA proposal to build a constellation of spacecraft to fly into the near-Earth solar wind in a swarm spanning a multitude of scales in order to obtain critically needed measurements that will reveal the underlying dynamics of magnetized turbulence. This white paper, submitted to the Plasma 2020 Decadal Survey Committee, provides a brief overview of turbulent systems that constitute an area of compelling plasma physics research, including why this mission is needed, and how this mission will achieve the goal of revealing how energy is transferred across scales and boundaries in plasmas throughout the universe.

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.

Transport properties, such as viscosity and thermal conduction, of the hot intergalactic plasma in clusters of galaxies, are largely unknown. While for laboratory plasmas these characteristics are derived from the gas density and temperature, such recipes can be fundamentally different for the intergalactic plasma due to a low rate of particle collisions and a weak magnetic field. In numerical simulations, one often cuts through these unknowns by modeling these plasmas as hydrodynamic fluids, even though local, non-hydrodynamic features observed in clusters contradict this assumption. Using deep Chandra observations of the Coma Cluster, we probe gas fluctuations in intergalactic medium down to spatial scales where the transport processes should prominently manifest themselves - at least if hydrodynamic models with pure Coulomb collision rates were indeed adequate. We find that they do not, implying that the effective isotropic viscosity is orders of magnitude smaller than naively expected. This indicates an enhanced collision rate in the plasma due to particle scattering off microfluctuations caused by plasma instabilities, or that the transport processes are anisotropic with respect to local magnetic field. For that reason, numerical models with high Reynolds number appear more consistent with observations. Our results also demonstrate that observations of turbulence in clusters are becoming a branch of astrophysics that can sharpen theoretical views on such plasmas.

It is shown that in low-beta, weakly collisional plasmas, such as the solar corona, some instances of the solar wind, the aurora, inner regions of accretion discs, their coronae, and some laboratory plasmas, Alfv\'enic fluctuations produce no ion heating within the gyrokinetic approximation, i.e., as long as their amplitudes (at the Larmor scale) are small and their frequencies stay below the ion Larmor frequency (even as their spatial scales can be above or below the ion Larmor scale). Thus, all low-frequency ion heating in such plasmas is due to compressive fluctuations ("slow modes"). Because these fluctuations energetically decouple from the Alfv\'enic ones already in the inertial range, the above conclusion means that the energy partition between ions and electrons in low-beta plasmas is decided at the outer scale, where turbulence is launched, and can be determined from magnetohydrodynamic (MHD) models of the relevant astrophysical systems. Any additional ion heating must come from non-gyrokinetic mechanisms such as cyclotron heating or the stochastic heating owing to distortions of ions' Larmor orbits. An exception to these conclusions occurs in the Hall limit, i.e., when the ratio of the ion to electron temperatures is as low as the ion beta (equivalently, the electron beta is order unity). In this regime, slow modes couple to Alfv\'enic ones well above the Larmor scale (viz., at the ion inertial or ion sound scale), so the Alfv\'enic and compressive cascades join and then separate again into two cascades of fluctuations that linearly resemble kinetic Alfv\'en and ion cyclotron waves, with the former heating electrons and the latter ions. The two cascades are shown to decouple, scalings for them are derived, and it is argued physically that the two species will be heated by them at approximately equal rates.

Galaxy clusters are massive dark matter-dominated systems filled with X-ray emitting, optically thin plasma. Their large size and relative simplicity (at least as astrophysical objects go) make them a unique laboratory to measure some of the interesting plasma properties that are inaccessible by other means but fundamentally important for understanding and modeling many astrophysical phenomena -- from solar flares to black hole accretion to galaxy formation and the emergence of the cosmological Large Scale Structure. While every cluster astrophysicist is eagerly anticipating the direct gas velocity measurements from the forthcoming microcalorimeters onboard XRISM, Athena and future missions such as Lynx, a number of those plasma properties can best be probed by high-resolution X-ray imaging of galaxy clusters. Chandra has obtained some trailblazing results, but only grazed the surface of such studies. In this white paper, we discuss why we need arcsecond-resolution, high collecting area, low relative background X-ray imagers (with modest spectral resolution), such as the proposed AXIS and the imaging detector of Lynx.

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.

This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed.

Does overall thermal equilibrium exist between ions and electrons in a weakly collisional, magnetised, turbulent plasma---and, if not, how is thermal energy partitioned between ions and electrons? This is a fundamental question in plasma physics, the answer to which is also crucial for predicting the properties of far-distant astronomical objects such as accretion discs around black holes. In the context of discs, this question was posed nearly two decades ago and has since generated a sizeable literature. Here we provide the answer for the case in which energy is injected into the plasma via Alfv\'enic turbulence: collisionless turbulent heating typically acts to disequilibrate the ion and electron temperatures. Numerical simulations using a hybrid fluid-gyrokinetic model indicate that the ion-electron heating-rate ratio is an increasing function of the thermal-to-magnetic energy ratio, $\beta_\mathrm{i}$: it ranges from $\sim0.05$ at $\beta_\mathrm{i}=0.1$ to at least $30$ for $\beta_\mathrm{i} \gtrsim 10$. This energy partition is approximately insensitive to the ion-to-electron temperature ratio $T_\mathrm{i}/T_\mathrm{e}$. Thus, in the absence of other equilibrating mechanisms, a collisionless plasma system heated via Alfv\'enic turbulence will tend towards a nonequilibrium state in which one of the species is significantly hotter than the other, viz., hotter ions at high $\beta_\mathrm{i}$, hotter electrons at low $\beta_\mathrm{i}$. Spectra of electromagnetic fields and the ion distribution function in 5D phase space exhibit an interesting new magnetically dominated regime at high $\beta_i$ and a tendency for the ion heating to be mediated by nonlinear phase mixing ("entropy cascade") when $\beta_\mathrm{i}\lesssim1$ and by linear phase mixing (Landau damping) when $\beta_\mathrm{i}\gg1$

We propose that pressure anisotropy causes weakly collisional turbulent plasmas to self-organize so as to resist changes in magnetic-field strength. We term this effect "magneto-immutability" by analogy with incompressibility (resistance to changes in pressure). The effect is important when the pressure anisotropy becomes comparable to the magnetic pressure, suggesting that in collisionless, weakly magnetized (high-$\beta$) plasmas its dynamical relevance is similar to that of incompressibility. Simulations of magnetized turbulence using the weakly collisional Braginskii model show that magneto-immutable turbulence is surprisingly similar, in most statistical measures, to critically balanced MHD turbulence. However, in order to minimize magnetic-field variation, the flow direction becomes more constrained than in MHD, and the turbulence is more strongly dominated by magnetic energy (a nonzero "residual energy"). These effects represent key differences between pressure-anisotropic and fluid turbulence, and should be observable in the $\beta\gtrsim1$ turbulent solar wind.

A heat flux in a high-$\beta$ plasma with low collisionality triggers the whistler instability. Quasilinear theory predicts saturation of the instability in a marginal state characterized by a heat flux that is fully controlled by electron scattering off magnetic perturbations. This marginal heat flux does not depend on the temperature gradient and scales as $1/\beta$. We confirm this theoretical prediction by performing numerical particle-in-cell simulations of the instability. We further calculate the saturation level of magnetic perturbations and the electron scattering rate as functions of $\beta$ and the temperature gradient to identify the saturation mechanism as quasilinear. Suppression of the heat flux is caused by oblique whistlers with magnetic-energy density distributed over a wide range of propagation angles. This result can be applied to high-$\beta$ astrophysical plasmas, such as the intracluster medium, where thermal conduction at sharp temperature gradients along magnetic-field lines can be significantly suppressed. We provide a convenient expression for the amount of suppression of the heat flux relative to the classical Spitzer value as a function of the temperature gradient and $\beta$. For a turbulent plasma, the additional independent suppression by the mirror instability is capable of producing large total suppression factors (several tens in galaxy clusters) in regions with strong temperature gradients.

We present a theoretical framework for describing electromagnetic kinetic turbulence in a multi-species, magnetized, pressure-anisotropic plasma. Turbulent fluctuations are assumed to be small compared to the mean field, to be spatially anisotropic with respect to it, and to have frequencies small compared to the ion cyclotron frequency. At scales above the ion Larmor radius, the theory reduces to the pressure-anisotropic generalization of kinetic reduced magnetohydrodynamics (KRMHD) formulated by Kunz et al. (2015). At scales at and below the ion Larmor radius, three main objectives are achieved. First, we analyse the linear response of the pressure-anisotropic gyrokinetic system, and show it to be a generalisation of previously explored limits. The effects of pressure anisotropy on the stability and collisionless damping of Alfvenic and compressive fluctuations are highlighted, with attention paid to the spectral location and width of the frequency jump that occurs as Alfven waves transition into kinetic Alfven waves. Secondly, we derive and discuss a general free-energy conservation law, which captures both the KRMHD free-energy conservation at long wavelengths and dual cascades of kinetic Alfven waves and ion entropy at sub-ion-Larmor scales. We show that non-Maxwellian features in the distribution function change the amount of phase mixing and the efficiency of magnetic stresses, and thus influence the partitioning of free energy amongst the cascade channels. Thirdly, a simple model is used to show that pressure anisotropy can cause large variations in the ion-to-electron heating ratio due to the dissipation of Alfvenic turbulence. Our theory provides a foundation for determining how pressure anisotropy affects the turbulent fluctuation spectra, the differential heating of particle species, and the ratio of parallel and perpendicular phase mixing in space and astrophysical plasmas.

Buoyant bubbles of relativistic plasma in cluster cores plausibly play a key role in conveying the energy from a supermassive black hole to the intracluster medium (ICM) - the process known as radio-mode AGN feedback. Energy conservation guarantees that a bubble loses most of its energy to the ICM after crossing several pressure scale heights. However, actual processes responsible for transferring the energy to the ICM are still being debated. One attractive possibility is the excitation of internal waves, which are trapped in the cluster's core and eventually dissipate. Here we show that a sufficient condition for efficient excitation of these waves in stratified cluster atmospheres is flattening of the bubbles in the radial direction. In our numerical simulations, we model the bubbles phenomenologically as rigid bodies buoyantly rising in the stratified cluster atmosphere. We find that the terminal velocities of the flattened bubbles are small enough so that the Froude number ${\rm Fr}\lesssim 1$. The effects of stratification make the dominant contribution to the total drag force balancing the buoyancy force. In particular, clear signs of internal waves are seen in the simulations. These waves propagate horizontally and downwards from the rising bubble, spreading their energy over large volumes of the ICM. If our findings are scaled to the conditions of the Perseus cluster, the expected terminal velocity is $\sim100-200{\,\rm km\,s^{-1}}$ near the cluster cores, which is in broad agreement with direct measurements by the Hitomi satellite.

A class of simple kinetic systems is considered, described by the 1D Vlasov-Landau equation with Poisson or Boltzmann electrostatic response and an energy source. Assuming a stochastic electric field, a solvable model is constructed for the phase-space turbulence of the particle distribution. The model is a kinetic analog of the Kraichnan-Batchelor model of chaotic advection. The solution of the model is found in Fourier-Hermite space and shows that the free-energy flux from low to high Hermite moments is suppressed, with phase mixing cancelled on average by anti-phase-mixing (stochastic plasma echo). This implies that Landau damping is an ineffective route to dissipation (i.e., to thermalisation of electric energy via velocity space). The full Fourier-Hermite spectrum is derived. Its asymptotics are $m^{-3/2}$ at low wave numbers and high Hermite moments ($m$) and $m^{-1/2}k^{-2}$ at low Hermite moments and high wave numbers ($k$). These conclusions hold at wave numbers below a certain cut off (analog of Kolmogorov scale), which increases with the amplitude of the stochastic electric field and scales as inverse square of the collision rate. The energy distribution and flows in phase space are a simple and, therefore, useful example of competition between phase mixing and nonlinear dynamics in kinetic turbulence, reminiscent of more realistic but more complicated multi-dimensional systems that have not so far been amenable to complete analytical solution.

Using two-dimensional hybrid-kinetic simulations, we explore the nonlinear "interruption" of standing and traveling shear-Alfv\'en waves in collisionless plasmas. Interruption involves a self-generated pressure anisotropy removing the restoring force of a linearly polarized Alfv\'enic perturbation, and occurs for wave amplitudes $\delta B_{\perp}/B_{0}\gtrsim \beta^{\,-1/2}$ (where $\beta$ is the ratio of thermal to magnetic pressure). We use highly elongated domains to obtain maximal scale separation between the wave and the ion gyroscale. For standing waves above the amplitude limit, we find that the large-scale magnetic field of the wave decays rapidly. The dynamics are strongly affected by the excitation of oblique firehose modes, which transition into long-lived parallel fluctuations at the ion gyroscale and cause significant particle scattering. Traveling waves are damped more slowly, but are also influenced by small-scale parallel fluctuations created by the decay of firehose modes. Our results demonstrate that collisionless plasmas cannot support linearly polarized Alfv\'en waves above $\delta B_{\perp}/B_{0}\sim \beta^{\,-1/2}$. They also provide a vivid illustration of two key aspects of low-collisionality plasma dynamics: (i) the importance of velocity-space instabilities in regulating plasma dynamics at high $\beta$, and (ii) how nonlinear collisionless processes can transfer mechanical energy directly from the largest scales into thermal energy and microscale fluctuations, without the need for a scale-by-scale turbulent cascade.

We describe polarization of the Sunyaev-Zel'dovich (SZ) effect associated with electron pressure anisotropy likely present in the intracluster medium (ICM). The ICM is an astrophysical example of a weakly collisional plasma where the Larmor frequencies of charged particles greatly exceed their collision frequencies. This permits formation of pressure anisotropies, driven by evolving magnetic fields via adiabatic invariance, or by heat fluxes. SZ polarization arises in the process of Compton scattering of the cosmic microwave background (CMB) photons off the thermal ICM electrons due to the difference in the characteristic thermal velocities of the electrons along two mutually orthogonal directions in the sky plane. The signal scales linearly with the optical depth of the region containing large-scale correlated anisotropy, and with the degree of anisotropy itself. It has the same spectral dependence as the polarization induced by cluster motion with respect to the CMB frame (kinematic SZ effect polarization), but can be distinguished by its spatial pattern. { For the illustrative case of a galaxy cluster with a cold front, where electron transport is mediated by Coulomb collisions, we estimate the CMB polarization degree at the level of 10$^{-8}$ ($\sim 10$ nK). An increase of the effective electron collisionality due to plasma instabilities will reduce the effect. Such polarization, therefore, may be an independent probe of the electron collisionality in the ICM, which is one of the key properties of a high-$\beta$ weakly collisional plasma from the point of view of both astrophysics and plasma theory.

We calculate the disruption scale $\lambda_{\rm D}$ at which sheet-like structures in dynamically aligned Alfv\'enic turbulence are destroyed by the onset of magnetic reconnection in a low-$\beta$ collisionless plasma. The scaling of $\lambda_{\rm D}$ depends on the order of the statistics being considered, with more intense structures being disrupted at larger scales. The disruption scale for the structures that dominate the energy spectrum is $\lambda_{\rm D}\sim L_\perp^{1/9}(d_e\rho_s)^{4/9}$, where $d_e$ is the electron inertial scale, $\rho_s$ is the ion sound scale, and $L_\perp$ is the outer scale of the turbulence. When $\beta_e$ and $\rho_s/L_\perp$ are sufficiently small, the scale $\lambda_{\rm D}$ is larger than $\rho_s$ and there is a break in the energy spectrum at $\lambda_{\rm D}$, rather than at $\rho_s$. We propose that the fluctuations produced by the disruption are circularised flux ropes, which may have already been observed in the solar wind. We predict the relationship between the amplitude and radius of these structures and quantify the importance of the disruption process to the cascade in terms of the filling fraction of undisrupted structures and the fractional reduction of the energy contained in them at the ion sound scale $\rho_s$. Both of these fractions depend strongly on $\beta_e$, with the disrupted structures becoming more important at lower $\beta_e$. Finally, we predict that the energy spectrum between $\lambda_{\rm D}$ and $\rho_s$ is steeper than $k_\perp^{-3}$, when this range exists. Such a steep "transition range" is sometimes observed in short intervals of solar-wind turbulence. The onset of collisionless magnetic reconnection may therefore significantly affect the nature of plasma turbulence around the ion gyroscale.

We propose a mechanism whereby the intense, sheet-like structures naturally formed by dynamically aligning Alfv\'enic turbulence are destroyed by magnetic reconnection at a scale $\hat{\lambda}_{\rm D}$, larger than the dissipation scale predicted by models of intermittent, dynamically aligning turbulence. The reconnection process proceeds in several stages: first, a linear tearing mode with $N$ magnetic islands grows and saturates, and then the $X$-points between these islands collapse into secondary current sheets, which then reconnect until the original structure is destroyed. This effectively imposes an upper limit on the anisotropy of the structures within the perpendicular plane, which means that at scale $\hat{\lambda}_{\rm D}$ the turbulent dynamics change: at scales larger than $\hat{\lambda}_{\rm D}$, the turbulence exhibits scale-dependent dynamic alignment and a spectral index approximately equal to $-3/2$, while at scales smaller than $\hat{\lambda}_{\rm D}$, the turbulent structures undergo a succession of disruptions due to reconnection, limiting dynamic alignment, steepening the effective spectral index and changing the final dissipation scale. The scaling of $\hat{\lambda}_{\rm D}$ with the Lundquist (magnetic Reynolds) number $S_{L_\perp}$ depends on the order of the statistics being considered, and on the specific model of intermittency; the transition between the two regimes in the energy spectrum is predicted at approximately $\hat{\lambda}_{\rm D} \sim S_{L_\perp}^{-0.6}$. The spectral index below $\hat{\lambda}_{\rm D}$ is bounded between $-5/3$ and $-2.3$. The final dissipation scale is at $\hat{\lambda}_{\eta,\infty}\sim S_{L_\perp}^{-3/4}$, the same as the Kolmogorov scale arising in theories of turbulence that do not involve scale-dependent dynamic alignment.

The Sun provides us with the only spatially well-resolved astrophysical example of turbulent thermal convection. While various aspects of solar photospheric turbulence, such as granulation (one-Megameter horizontal scale), are well understood, the questions of the physical origin and dynamical organization of larger-scale flows, such as the 30-Megameters supergranulation and flows deep in the solar convection zone, remain largely open in spite of their importance for solar dynamics and magnetism. Here, we present a new critical global observational characterization of multiscale photospheric flows and subsequently formulate an anisotropic extension of the Bolgiano-Obukhov theory of hydrodynamic stratified turbulence that may explain several of their distinctive dynamical properties. Our combined analysis suggests that photospheric flows in the horizontal range of scales between supergranulation and granulation have a typical vertical correlation scale of 2.5 to 4 Megameters and operate in a strongly anisotropic, self-similar, nonlinear, buoyant dynamical regime. While the theory remains speculative at this stage, it lends itself to quantitative comparisons with future high-resolution acoustic tomography of subsurface layers and advanced numerical models. Such a validation exercise may also lead to new insights into the asymptotic dynamical regimes in which other, unresolved turbulent anisotropic astrophysical fluid systems supporting waves or instabilities operate.

Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $\mu$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. Here we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization of the Universe.

The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and study dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics.

This work, which extends Squire et al. [ApJL, 830 L25 (2016)], explores the effect of self-generated pressure anisotropy on linearly polarized shear-Alfv\'en fluctuations in low-collisionality plasmas. Such anisotropies lead to stringent limits on the amplitude of magnetic perturbations in high-beta plasmas, above which a fluctuation can destabilize itself through the parallel firehose instability. This causes the wave frequency to approach zero, "interrupting" the wave and stopping its oscillation. These effects are explored in detail in the collisionless and weakly collisional "Braginskii" regime, for both standing and traveling waves. The focus is on simplified models in one dimension, on scales much larger than the ion gyroradius. The effect has interesting implications for the physics of magnetized turbulence in the high-beta conditions that are prevalent in many astrophysical plasmas.

We propose a simple statistical model of three-dimensionally anisotropic, intermittent, strong Alfv\'enic turbulence, incorporating both critical balance and dynamic alignment. Our model is based on log-Poisson statistics for Elsasser-field increments along the magnetic field. We predict the scalings of Elsasser-field conditional two-point structure functions with point separations in all three directions in a coordinate system locally aligned with the direction of the magnetic field and of the fluctuating fields and obtain good agreement with numerical simulations. We also derive a scaling of the parallel coherence scale of the fluctuations, $l_\parallel \propto \lambda^{1/2}$, where $\lambda$ is the perpendicular scale. This is indeed observed for the bulk of the fluctuations in numerical simulations.

It is shown that low-collisionality plasmas cannot support linearly polarized shear-Alfv\'en fluctuations above a critical amplitude $\delta B_{\perp}/B_{0} \sim \beta^{\,-1/2}$, where $\beta$ is the ratio of thermal to magnetic pressure. Above this cutoff, a developing fluctuation will generate a pressure anisotropy that is sufficient to destabilize itself through the parallel firehose instability. This causes the wave frequency to approach zero, interrupting the fluctuation before any oscillation. The magnetic field lines rapidly relax into a sequence of angular zig-zag structures. Such a restrictive bound on shear-Alfv\'en-wave amplitudes has far-reaching implications for the physics of magnetized turbulence in the high-$\beta$ conditions prevalent in many astrophysical plasmas, as well as for the solar wind at $\sim 1 \mathrm{AU}$ where $\beta \gtrsim 1$.

The firehose and mirror instabilities are thought to arise in a variety of space and astrophysical plasmas, constraining the pressure anisotropies and drifts between particle species. The plasma stability depends on all species simultaneously, meaning that a combined analysis is required. Here, we present the first such analysis in the solar wind, using the long-wavelength stability parameters to combine the anisotropies and drifts of all major species (core and beam protons, alphas, and electrons). At the threshold, the firehose parameter was found to be dominated by protons (67%), but also to have significant contributions from electrons (18%) and alphas (15%). Drifts were also found to be important, contributing 57% in the presence of a proton beam. A similar situation was found for the mirror, with contributions of 61%, 28%, and 11% for protons, electrons, and alphas, respectively. The parallel electric field contribution, however, was found to be small at 9%. Overall, the long-wavelength thresholds constrain the data well (<1% unstable), and the implications of this are discussed.

Astrophysical plasmas are typically magnetized, with the Larmor radii of the charged particles many orders of magnitude smaller than their collisional mean free paths. The fundamental properties of such plasmas, e.g., conduction and viscosity, may depend on the instabilities driven by the anisotropy of the particle distribution functions and operating at scales comparable to the Larmor scales. We discuss a possibility that the pressure anisotropy of thermal electrons could produce polarization of thermal bremsstrahlung emission. In particular, we consider coherent large-scale motions in galaxy clusters to estimate the level of anisotropy driven by stretching of the magnetic-field lines by plasma flow and by heat fluxes associated with thermal gradients. Our estimate of the degree of polarization is $\sim 0.1 \%$ at energies $\gtrsim kT$. While this value is too low for the forthcoming generation of X-ray polarimeters, it is potentially an important proxy for the processes taking place at extremely small scales, which are impossible to resolve spatially. The absence of the effect at the predicted level may set a lower limit on the electron collisionality in the ICM. At the same time, the small value of the effect implies that it does not preclude the use of clusters as (unpolarized) calibration sources for X-ray polarimeters at this level of accuracy.

The ICM plasma is subject to firehose and mirror instabilities at scales of order the ion Larmor radius. The mirror instability generates fluctuations of magnetic-field strength $\delta B / B \sim 1$. These fluctuations act as magnetic traps for the heat-conducting electrons, suppressing their transport. We calculate the effective parallel thermal conductivity in the ICM in the presence of the mirror fluctuations for different stages of the evolution of the instability. The mirror fluctuations are limited in amplitude by the maximum and minimum values of the field strength, with no large deviations from the mean value. This key property leads to a finite suppression of thermal conduction at large scales. We find suppression down to $\approx 0.2$ of the Spitzer value for the secular phase of the perturbations' growth, and $\approx 0.3$ for their saturated phase. The effect operates in addition to other suppression mechanisms and independently of them. Globally, fluctuations $\delta B / B \sim 1$ can be present on much larger scales, of the order of the scale of turbulent motions. However, we do not expect large suppression of thermal conduction by these, because their scale is considerably larger than the collisional mean free path of the ICM electrons. The obtained suppression of thermal conduction by a factor of $\sim 5$ appears to be characteristic and potentially universal for a weakly collisional mirror-unstable plasma.

We discuss a novel technique of manipulating X-ray images of galaxy clusters to reveal the nature of small-scale density/temperature perturbations in the intra cluster medium (ICM). As we show, this technique can be used to differentiate between sound waves and isobaric perturbations in Chandra images of the Perseus and M87/Virgo clusters. The comparison of the manipulated images with the radio data and with the results of detailed spectral analysis shows that this approach successfully classifies the types of perturbations and helps to reveal their nature. For the central regions (5-100 kpc) of the M87 and Perseus clusters this analysis suggests that observed images are dominated by isobaric perturbations, followed by perturbations caused by bubbles of relativistic plasma and weak shocks. Such a hierarchy is best explained in a "slow" AGN feedback scenario, when much of the mechanical energy output of a central black hole is captured by the bubble enthalpy that is gradually released during buoyant rise of the bubbles. The "image arithmetic" works best for prominent structure and for datasets with excellent statistics, visualizing the perturbations with a given effective equation of state. The same approach can be extended to faint perturbations via cross-spectrum analysis of surface brightness fluctuations in X-ray images in different energy bands.

A scaling theory of long-wavelength electrostatic turbulence in a magnetised, weakly collisional plasma (e.g., ITG turbulence) is proposed, with account taken both of the nonlinear advection of the perturbed particle distribution by fluctuating ExB flows and of its phase mixing, which is caused by the streaming of the particles along the mean magnetic field and, in a linear problem, would lead to Landau damping. It is found that it is possible to construct a consistent theory in which very little free energy leaks into high velocity moments of the distribution function, rendering the turbulent cascade in the energetically relevant part of the wave-number space essentially fluid-like. The velocity-space spectra of free energy expressed in terms of Hermite-moment orders are steep power laws and so the free-energy content of the phase space does not diverge at infinitesimal collisionality (while it does for a linear problem); collisional heating due to long-wavelength perturbations vanishes in this limit (also in contrast with the linear problem, in which it occurs at the finite rate equal to the Landau-damping rate). The ability of the free energy to stay in the low velocity moments of the distribution function is facilitated by the "anti-phase-mixing" effect, whose presence in the nonlinear system is due to the stochastic version of the plasma echo (the advecting velocity couples the phase-mixing and anti-phase-mixing perturbations). The partitioning of the wave-number space between the (energetically dominant) region where this is the case and the region where linear phase mixing wins its competition with nonlinear advection is governed by the "critical balance" between linear and nonlinear timescales (which for high Hermite moments splits into two thresholds, one demarcating the wave-number region where phase mixing predominates, the other where plasma echo does).

Magnetic fields pervade the entire Universe and affect the formation and evolution of astrophysical systems from cosmological to planetary scales. The generation and dynamical amplification of extragalactic magnetic fields through cosmic times, up to $\mu$Gauss levels reported in nearby galaxy clusters, near equipartition with kinetic energy of plasma motions and on scales of at least tens of kiloparsecs, is a major puzzle largely unconstrained by observations. A dynamo effect converting kinetic flow energy into magnetic energy is often invoked in that context, however extragalactic plasmas are weakly collisional (as opposed to magnetohydrodynamic fluids), and whether magnetic-field growth and sustainment through an efficient turbulent dynamo instability is possible in such plasmas is not established. Fully kinetic numerical simulations of the Vlasov equation in a six-dimensional phase space necessary to answer this question have until recently remained beyond computational capabilities. Here, we show by means of such simulations that magnetic-field amplification via a dynamo instability does occur in a stochastically-driven, non-relativistic subsonic flow of initially unmagnetized collisionless plasma. We also find that the dynamo self-accelerates and becomes entangled with kinetic instabilities as magnetization increases. The results suggest that such a plasma dynamo may be realizable in laboratory experiments, support the idea that intracluster medium (ICM) turbulence may have significantly contributed to the amplification of cluster magnetic fields up to near-equipartition levels on a timescale shorter than the Hubble time, and emphasize the crucial role of multiscale kinetic physics in high-energy astrophysical plasmas.

Cores of relaxed galaxy clusters are often disturbed by AGN. Their Chandra observations revealed a wealth of structures induced by shocks, subsonic gas motions, bubbles of relativistic plasma, etc. In this paper, we determine the nature and energy content of gas fluctuations in the Perseus core by probing statistical properties of emissivity fluctuations imprinted in the soft- and hard-band X-ray images. About 80 per cent of the total variance of perturbations on ~ 8-70 kpc scales in the inner region have an isobaric nature, i.e., are consistent with slow displacements of the gas in pressure equilibrium with ambient medium. Observed variance translates to the ratio of non-thermal to thermal energy of ~13 per cent. In the region dominated by weak "ripples", about half of the total variance is also associated with isobaric perturbations on scales ~ a few tens of kpc. If these isobaric perturbations are induced by buoyantly rising bubbles, then these results suggest that most of the AGN-injected energy should first go into bubbles rather than into shocks. Using simulations of a shock propagating through the Perseus atmosphere, we found that models reproducing the observed features of a central shock have more than 50 per cent of the AGN-injected energy associated with the bubble enthalpy and only about 20 per cent is carried away with the shock. Such energy partition is consistent with the AGN-feedback model, mediated by bubbles of relativistic plasma, and supports the importance of turbulence in the balance between gas heating and radiative cooling.

The nonlinear state of a high-beta collisionless plasma is investigated when an imposed linear shear amplifies or diminishes a uniform magnetic field, driving pressure anisotropies and hence firehose/mirror instabilities. The evolution of the resulting microscale turbulence is considered when the shear is switched off or reversed after one shear time (mimicking local behaviour of a macroscopic flow), so a new macroscale configuration is superimposed on the microscale state left behind by the previous one. There is a threshold value of plasma beta: when $\beta\ll\Omega/S$ (ion cyclotron frequency/shear rate), the emergence of firehose/mirror fluctuations driven unstable by shear and their disappearance when the shear is removed/reversed are quasi-instantaneous compared to the shear time, viz., the decay time of these fluctuations is $\sim\beta/\Omega \ll 1/S$ (this result follows from the free decay of the fluctuations being constrained by the same marginal-stability thresholds as their growth). In contrast, when $\beta\gtrsim\Omega/S$ ("ultra-high" beta), the old microscale state can only be removed on the shear timescale. In this regime, driven firehose fluctuations grow secularly to order-unity amplitudes, compensating for the decay of the mean field and so pinning pressure anisotropy at marginal stability with no appreciable scattering of particles---which is unlike what happens at moderate $\beta$. When the shear reverses, the shearing away of this firehose turbulence compensates for the increase in the mean field and thus prevents growth of the pressure anisotropy, stopping the system from going mirror-unstable. Therefore, at ultra-high beta, the system stays close to the firehose threshold, the mirror instability is largely suppressed, while the mean magnetic energy barely changes at all. Implications for plasma dynamo and thus the origin of cosmic magnetism are discussed.

We measure the local anisotropy of numerically simulated strong Alfv\'enic turbulence with respect to two local, physically relevant directions: along the local mean magnetic field and along the local direction of one of the fluctuating Elsasser fields. We find significant scaling anisotropy with respect to both these directions: the fluctuations are "ribbon-like" --- statistically, they are elongated along both the mean magnetic field and the fluctuating field. The latter form of anisotropy is due to scale-dependent alignment of the fluctuating fields. The intermittent scalings of the $n$th-order conditional structure functions in the direction perpendicular to both the local mean field and the fluctuations agree well with the theory of Chandran et al. 2015, while the parallel scalings are consistent with those implied by the critical-balance conjecture. We quantify the relationship between the perpendicular scalings and those in the fluctuation and parallel directions, and find that the scaling exponent of the perpendicular anisotropy (i.e., of the aspect ratio of the Alfv\'enic structures in the plane perpendicular to the mean magnetic field) depends on the amplitude of the fluctuations. This is shown to be equivalent to the anticorrelation of fluctuation amplitude and alignment at each scale. The dependence of the anisotropy on amplitude is shown to be more significant for the anisotropy between the perpendicular and fluctuation-direction scales than it is between the perpendicular and parallel scales.

Central jetted active galactic nuclei (AGN) appear to heat the core regions of the intracluster medium (ICM) in cooling-core galaxy clusters and groups, thereby preventing a cooling catastrophe. However, the physical mechanism(s) by which the directed flow of kinetic energy is thermalized throughout the ICM core remains unclear. We examine one widely discussed mechanism whereby the AGN induces subsonic turbulence in the ambient medium, the dissipation of which provides the ICM heat source. Through controlled inviscid 3-d hydrodynamic simulations, we verify that explosive AGN-like events can launch gravity waves (g-modes) into the ambient ICM which in turn decay to volume-filling turbulence. In our model, however, this process is found to be inefficient, with less than 1% of the energy injected by the AGN activity actually ending up in the turbulence of the ambient ICM. This efficiency is an order of magnitude or more too small to explain the observations of AGN-feedback in galaxy clusters and groups with short central cooling times. Atmospheres in which the g-modes are strongly trapped/confined have an even lower efficiency since, in these models, excitation of turbulence relies on the g-modes' ability to escape from the center of the cluster into the bulk ICM. Our results suggest that, if AGN-induced turbulence is indeed the mechanism by which the AGN heats the ICM core, its driving may rely on physics beyond that captured in our ideal hydrodynamic model.

We present numerical evidence that in strong Alfvenic turbulence, the critical balance principle---equality of the nonlinear decorrelation and linear propagation times---is scale invariant, in the sense that the probability distribution of the ratio of these times is independent of scale. This result only holds if the local alignment of the Elsasser fields is taken into account in calculating the nonlinear time. At any given scale, the degree of alignment is found to increase with fluctuation amplitude, supporting the idea that the cause of alignment is mutual dynamical shearing of Elsasser fields. The scale-invariance of critical balance (while all other quantities of interest are strongly intermittent, i.e., have scale-dependent distributions) suggests that it is the most robust of the scaling principles used to describe Alfvenic turbulence. The quality afforded by situ fluctuation measurements in the solar wind allows for direct verification of this fundamental principle.

The Wisconsin Plasma Astrophysics Laboratory (WiPAL) is a flexible user facility designed to study a range of astrophysically relevant plasma processes as well as novel geometries that mimic astrophysical systems. A multi-cusp magnetic bucket constructed from strong samarium cobalt permanent magnets now confines a 10 m$^3$, fully ionized, magnetic-field free plasma in a spherical geometry. Plasma parameters of $ T_{e}\approx5$ to $20$ eV and $n_{e}\approx10^{11}$ to $5\times10^{12}$ cm$^{-3}$ provide an ideal testbed for a range of astrophysical experiments including self-exciting dynamos, collisionless magnetic reconnection, jet stability, stellar winds, and more. This article describes the capabilities of WiPAL along with several experiments, in both operating and planning stages, that illustrate the range of possibilities for future users.

A theoretical framework for low-frequency electromagnetic (drift-)kinetic turbulence in a collisionless, multi-species plasma is presented. The result generalises reduced magnetohydrodynamics (RMHD) and kinetic RMHD (Schekochihin et al. 2009) for pressure-anisotropic plasmas, allowing for species drifts---a situation routinely encountered in the solar wind and presumably ubiquitous in hot dilute astrophysical plasmas (e.g. intracluster medium). Two main objectives are achieved. First, in a non-Maxwellian plasma, the relationships between fluctuating fields (e.g., the Alfven ratio) are order-unity modified compared to the more commonly considered Maxwellian case, and so a quantitative theory is developed to support quantitative measurements now possible in the solar wind. The main physical feature of low-frequency plasma turbulence survives the generalisation to non-Maxwellian distributions: Alfvenic and compressive fluctuations are energetically decoupled, with the latter passively advected by the former; the Alfvenic cascade is fluid, satisfying RMHD equations (with the Alfven speed modified by pressure anisotropy and species drifts), whereas the compressive cascade is kinetic and subject to collisionless damping. Secondly, the organising principle of this turbulence is elucidated in the form of a generalised kinetic free-energy invariant. It is shown that non-Maxwellian features in the distribution function reduce the rate of phase mixing and the efficacy of magnetic stresses; these changes influence the partitioning of free energy amongst the various cascade channels. As the firehose or mirror instability thresholds are approached, the dynamics of the plasma are modified so as to reduce the energetic cost of bending magnetic-field lines or of compressing/rarefying them. Finally, it is shown that this theory can be derived as a long-wavelength limit of non-Maxwellian slab gyrokinetics.

X-ray surface brightness fluctuations in the core of the Perseus Cluster are analyzed, using deep observations with the Chandra observatory. The amplitude of gas density fluctuations on different scales is measured in a set of radial annuli. It varies from 8 to 12 per cent on scales of ~10-30 kpc within radii of 30-160 kpc from the cluster center and from 9 to 7 per cent on scales of ~20-30 kpc in an outer, 60-220 kpc annulus. Using a statistical linear relation between the observed amplitude of density fluctuations and predicted velocity, the characteristic velocity of gas motions on each scale is calculated. The typical amplitudes of the velocity outside the central 30 kpc region are 90-140 km/s on ~20-30 kpc scales and 70-100 km/s on smaller scales ~7-10 kpc. The velocity power spectrum is consistent with cascade of turbulence and its slope is in a broad agreement with the slope for canonical Kolmogorov turbulence. The gas clumping factor estimated from the power spectrum of the density fluctuations is lower than 7-8 per cent for radii ~30-220 kpc from the center, leading to a density bias of less than 3-4 per cent in the cluster core. Uncertainties of the analysis are examined and discussed. Future measurements of the gas velocities with the Astro-H, Athena and Smart-X observatories will directly measure the gas density-velocity perturbation relation and further reduce systematic uncertainties in these quantities.

The hot, X-ray-emitting intracluster medium (ICM) is the dominant baryonic constituent of clusters of galaxies. In the cores of many clusters, radiative energy losses from the ICM occur on timescales significantly shorter than the age of the system. Unchecked, this cooling would lead to massive accumulations of cold gas and vigorous star formation, in contradiction to observations. Various sources of energy capable of compensating these cooling losses have been proposed, the most promising being heating by the supermassive black holes in the central galaxies through inflation of bubbles of relativistic plasma. Regardless of the original source of energy, the question of how this energy is transferred to the ICM has remained open. Here we present a plausible solution to this question based on deep Chandra X-ray observatory data and a new data-analysis method that enables us to evaluate directly the ICM heating rate due to the dissipation of turbulence. We find that turbulent heating is sufficient to offset radiative cooling and indeed appears to balance it locally at each radius - it might therefore be the key element in resolving the gas cooling problem in cluster cores and, more universally, in atmospheres of X-ray gas-rich systems.

Slow dynamical changes in magnetic-field strength and invariance of the particles' magnetic moments generate ubiquitous pressure anisotropies in weakly collisional, magnetized astrophysical plasmas. This renders them unstable to fast, small-scale mirror and firehose instabilities, which are capable of exerting feedback on the macroscale dynamics of the system. By way of a new asymptotic theory of the early nonlinear evolution of the mirror instability in a plasma subject to slow shearing or compression, we show that the instability does not saturate quasilinearly at a steady, low-amplitude level. Instead, the trapping of particles in small-scale mirrors leads to nonlinear secular growth of magnetic perturbations, $\delta B/B \propto t^{2/3}$. Our theory explains recent collisionless simulation results, provides a prediction of the mirror evolution in weakly collisional plasmas and establishes a foundation for a theory of nonlinear mirror dynamics with trapping, valid up to $\delta B/B =O(1)$.

Hybrid-kinetic numerical simulations of firehose and mirror instabilities in a collisionless plasma are performed in which pressure anisotropy is driven as the magnetic field is changed by a persistent linear shear $S$. For a decreasing field, it is found that mostly oblique firehose fluctuations grow at ion Larmor scales and saturate with energies $\sim$$S^{1/2}$; the pressure anisotropy is pinned at the stability threshold by particle scattering off microscale fluctuations. In contrast, nonlinear mirror fluctuations are large compared to the ion Larmor scale and grow secularly in time; marginality is maintained by an increasing population of resonant particles trapped in magnetic mirrors. After one shear time, saturated order-unity magnetic mirrors are formed and particles scatter off their sharp edges. Both instabilities drive sub-ion-Larmor--scale fluctuations, which appear to be kinetic-Alfv\'{e}n-wave turbulence. Our results impact theories of momentum and heat transport in astrophysical and space plasmas, in which the stretching of a magnetic field by shear is a generic process.

We address the problem of evaluating the power spectrum of the velocity field of the ICM using only information on the plasma density fluctuations, which can be measured today by Chandra and XMM-Newton observatories. We argue that for relaxed clusters there is a linear relation between the rms density and velocity fluctuations across a range of scales, from the largest ones, where motions are dominated by buoyancy, down to small, turbulent scales: $(\delta\rho_k/\rho)^2 = \eta_1^2 (V_{1,k}/c_s)^2$, where $\delta\rho_k/\rho$ is the spectral amplitude of the density perturbations at wave number $k$, $V_{1,k}^2=V_k^2/3$ is the mean square component of the velocity field, $c_s$ is the sound speed, and $\eta_1$ is a dimensionless constant of order unity. Using cosmological simulations of relaxed galaxy clusters, we calibrate this relation and find $\eta_1\approx 1 \pm 0.3$. We argue that this value is set at large scales by buoyancy physics, while at small scales the density and velocity power spectra are proportional because the former are a passive scalar advected by the latter. This opens an interesting possibility to use gas density power spectra as a proxy for the velocity power spectra in relaxed clusters, across a wide range of scales.

X-ray observations of hot gas in galaxy clusters often show steeper temperature gradients across cold fronts -- contact discontinuities, driven by the differential gas motions. These sharp (a few kpc wide) surface brightness/temperature discontinuities would be quickly smeared out by the electron thermal conduction in unmagnetized plasma, suggesting significant suppression of the heat flow across the discontinuities. In fact, the character of the gas flow near cold fronts is favorable for suppression of conduction by aligning magnetic field lines along the discontinuities. We argue that a similar mechanism is operating in the bulk of the gas. Generic 3D random isotropic and incompressible motions increase the temperature gradients (in some places) and at the same time suppress the local conduction by aligning the magnetic field lines perpendicular to the local temperature gradient. We show that the suppression of the effective conductivity in the bulk of the gas can be linked to the increase of the frozen magnetic field energy density. On average the rate of decay of the temperature fluctuations $\mathrm{d}\langle \delta T^2 \rangle /\mathrm{d}t$ decreases as $\langle B^2 \rangle ^{-1/5}$.

We develop an analytic model of intermittent, three-dimensional, strong, reduced magnetohydrodynamic turbulence with zero cross helicity. We take the fluctuation amplitudes to have a log-Poisson distribution and incorporate into the model a new phenomenology of scale-dependent dynamic alignment. The log-Poisson distribution in our model is characterized by two parameters. To calculate these parameters, we make use of three assumptions: that the energy cascade rate is independent of scale within the inertial range, that the most intense coherent structures at scale $\lambda$ are sheet-like with a volume filling factor proportional to $\lambda$, and that most of the cascade power arises from interactions between exceptionally intense fluctuations and much weaker fluctuations. We then compute the scalings of the power spectrum, the kurtosis, higher-order structure functions, and three different average alignment angles. These scalings appear to be consistent with numerous results from direct numerical simulations.

In weakly collisional extragalactic plasmas such as the intracluster medium, viscous stress and the rate of change of the magnetic-field strength are proportional to the local pressure anisotropy, so subject to constraints imposed by the pressure-anisotropy-driven mirror and firehose instabilities and controlled by the local instantaneous plasma beta. The dynamics of such plasmas is dramatically different from a conventional MHD fluid. The plasma is expected to stay locally in a marginal state with respect to the instabilities, but how it does this is an open question. Two models of magnetic-field evolution are investigated. In the first, marginality is achieved via suppression of the rate of change of the field. In the second, the instabilities give rise to anomalous collisionality, reducing pressure anisotropy to marginal - at the same time decreasing viscosity and so increasing the turbulent rate of strain. Implications of these models are studied in a simplified 0D setting. In the first model, the field grows explosively but on a time scale that scales with initial beta, while in the second, dynamical field strength can be reached in one large-scale turbulence turn-over time regardless of the initial seed. Both models produce very intermittent fields. Both also suffer from strong constraints on their applicability: for typical cluster-core conditions, scale separation between the fluid motions and the microscale fluctuations breaks down at beta~10^5-10^4. At larger beta (weaker fields), a fully collisionless plasma dynamo theory is needed in order to justify the growth of the field from a tiny primordial seed. However, the models discussed here are appropriate for studying the structure of the currently observed field as well as large-scale dynamics and thermodynamics of the magnetized ICM or similarly dilute astrophysical plasmas.

The hot X-ray emitting plasma in galaxy clusters is predicted to have turbulent motions which can contribute around ten percent of the cluster's central energy density. We report deep Chandra X-ray Observatory observations of the Coma cluster core, showing the presence of quasi-linear high-density arms spanning 150 kpc, consisting of low-entropy material likely stripped from merging subclusters. Two appear to be connected with a subgroup of galaxies at 650 kpc radius that is merging into the cluster, implying coherence over several hundred Myr. Such long lifetime implies that strong isotropic turbulence and conduction are suppressed in the core, despite the unrelaxed state of the cluster. Magnetic fields are presumably responsible. The structures seen in Coma present insight into the past Gyr of subcluster merger activity.

Cool-core clusters (e.g., Perseus or M87) often possess a network of bright gaseous filaments, observed in radio, infrared, optical and X-ray bands. We propose that these filaments are powered by the reconnection of the magnetic field in the wakes of buoyant bubbles. AGN-inflated bubbles of relativistic plasma rise buoyantly in the cluster atmosphere, stretching and amplifying the field in the wake to values of $\displaystyle \beta =8\pi P_{gas}/B^2\sim 1$. The field lines in the wake have opposite directions and are forced together as the bubble motion stretches the filament. This setup bears strong similarity to the coronal loops on the Sun or to the Earth's magneto-tail. The reconnection process naturally explains both the required level of local dissipation rate in filaments and the overall luminosity of filaments. The original source of power for the filaments is the potential energy of buoyant bubbles, inflated by the central AGN.

Magnetic reconnection in strongly magnetized (low-beta), weakly collisional plasmas is investigated using a novel fluid-kinetic model [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] which retains non-isothermal electron kinetics. It is shown that electron heating via Landau damping (linear phase mixing) is the dominant dissipation mechanism. In time, electron heating occurs after the peak of the reconnection rate; in space, it is concentrated along the separatrices of the magnetic island. For sufficiently large systems, the peak reconnection rate is $cE_{max}\approx 0.2v_AB_{y,0}$, where $v_A$ is the Alfv\'en speed based on the reconnecting field $B_{y,0}$. The island saturation width is the same as in MHD models except for small systems, when it becomes comparable to the kinetic scales.

We present a measurement of the scale-dependent, three-dimensional structure of the magnetic field fluctuations in inertial range solar wind turbulence with respect to a local, physically motivated coordinate system. The Alfvenic fluctuations are three-dimensionally anisotropic, with the sense of this anisotropy varying from large to small scales. At the outer scale, the magnetic field correlations are longest in the local fluctuation direction, consistent with Alfven waves. At the proton gyroscale, they are longest along the local mean field direction and shortest in the direction perpendicular to the local mean field and the local field fluctuation. The compressive fluctuations are highly elongated along the local mean field direction, although axially symmetric perpendicular to it. Their large anisotropy may explain why they are not heavily damped in the solar wind.

A 2D linear theory of the instability of Sweet-Parker (SP) current sheets is developed in the framework of Reduced MHD. A local analysis is performed taking into account the dependence of a generic equilibrium profile on the outflow coordinate. The plasmoid instability [Loureiro et al, Phys. Plasmas {\bf 14}, 100703 (2007)] is recovered, i.e., current sheets are unstable to the formation of a large-wave-number chain of plasmoids ($k_{\rm max}\Lsheet \sim S^{3/8}$, where $k_{\rm max}$ is the wave-number of fastest growing mode, $S=\Lsheet V_A/\eta$ is the Lundquist number, $\Lsheet$ is the length of the sheet, $V_A$ is the Alfv\'en speed and $\eta$ is the plasma resistivity), which grows super-Alfv\'enically fast ($\gmax\tau_A\sim S^{1/4}$, where $\gmax$ is the maximum growth rate, and $\tau_A=\Lsheet/V_A$). For typical background profiles, the growth rate and the wave-number are found to {\it increase} in the outflow direction. This is due to the presence of another mode, the Kelvin-Helmholtz (KH) instability, which is triggered at the periphery of the layer, where the outflow velocity exceeds the Alfv\'en speed associated with the upstream magnetic field. The KH instability grows even faster than the plasmoid instability, $\gmax \tau_A \sim k_{\rm max} \Lsheet\sim S^{1/2}$. The effect of viscosity ($\nu$) on the plasmoid instability is also addressed. In the limit of large magnetic Prandtl numbers, $Pm=\nu/\eta$, it is found that $\gmax\sim S^{1/4}Pm^{-5/8}$ and $k_{\rm max} \Lsheet\sim S^{3/8}Pm^{-3/16}$, leading to the prediction that the critical Lundquist number for plasmoid instability in the $Pm\gg1$ regime is $\Scrit\sim 10^4Pm^{1/2}$. These results are verified via direct numerical simulation of the linearized equations, using a new, analytical 2D SP equilibrium solution.

X-ray surface brightness fluctuations in the core ($650 \times 650$ kpc) region of the Coma cluster observed with XMM-Newton and Chandra are analyzed using a 2D power spectrum approach. The resulting 2D spectra are converted to 3D power spectra of gas density fluctuations. Our independent analyses of the XMM-Newton and Chandra observations are in excellent agreement and provide the most sensitive measurements of surface brightness and density fluctuations for a hot cluster. We find that the characteristic amplitude of the volume filling density fluctuations relative to the smooth underlying density distribution varies from 7-10% on scales of $\sim$500 kpc down to $\sim$5% at scales $\sim$ 30 kpc. On smaller spatial scales, projection effects smear the density fluctuations by a large factor, precluding strong limits on the fluctuations in 3D. On the largest scales probed (hundreds of kpc), the dominant contributions to the observed fluctuations most likely arise from perturbations of the gravitational potential by the two most massive galaxies in Coma, NGC4874 and NGC4889, and the low entropy gas brought to the cluster by an infalling group. Other plausible sources of X-ray surface brightness fluctuations are discussed, including turbulence, metal abundance variations, and unresolved sources. Despite a variety of possible origins for density fluctuations, the gas in the Coma cluster core is remarkably homogeneous on scales from $\sim$ 500 to $\sim$30 kpc.

Weak Alfvenic turbulence in a periodic domain is considered as a mixed state of Alfven waves interacting with the two-dimensional (2D) condensate. Unlike in standard treatments, no spectral continuity between the two is assumed and indeed none is found. If the 2D modes are not directly forced, k^{-2} and k^{-1} spectra are found for the Alfven waves and the 2D modes, respectively, with the latter less energetic than the former. The wave number at which their energies become comparable marks the transition to strong turbulence. For imbalanced energy injection, the spectra are similar and the Elsasser ratio scales as the ratio of the energy fluxes in the counterpropagting Alfven waves. If the 2D modes are forced, a 2D inverse cascade dominates the dynamics at the largest scales, but at small enough scales, the same weak and then strong regimes as described above are achieved.

A minimal model for magnetic reconnection and, generally, low-frequency dynamics in low-beta plasmas is proposed. The model combines analytical and computational simplicity with physical realizability: it is a rigorous limit of gyrokinetics for plasma beta of order the electron-ion mass ratio. The model contains collisions and can be used both in the collisional and collisionless reconnection regimes. It includes gyrokinetic ions (not assumed cold) and allows for the topological rearrangement of the magnetic field lines by either resistivity or electron inertia, whichever predominates. The two-fluid dynamics are coupled to electron kinetics --- electrons are not assumed isothermal and are described by a reduced drift-kinetic equation. The model therefore allows for irreversibility and conversion of magnetic energy into electron heat via parallel phase mixing in velocity space. An analysis of the exchanges between various forms of free energy and its conversion into electron heat is provided. It is shown how all relevant linear waves and regimes of the tearing instability (collisionless, semicollisional and fully resistive) are recovered in various limits of our model. An efficient way to simulate our equations numerically is proposed, via the Hermite representation of the velocity space. It is shown that small scales in velocity space will form, giving rise to a shallow Hermite-space spectrum, whence it is inferred that, for steady-state or sufficiently slow dynamics, the electron heating rate will remain finite in the limit of vanishing collisionality.

A rigorous theory for the generation of a large-scale magnetic field by random non-helically forced motions of a conducting fluid combined with a linear shear is presented in the analytically tractable limit of low Rm and weak shear. The dynamo is kinematic and due to fluctuations in the net (volume-averaged) electromotive force. This is a minimal proof-of-concept quasilinear calculation aiming to put the shear dynamo, a new effect recently found in numerical experiments, on a firm theoretical footing. Numerically observed scalings of the wavenumber and growth rate of the fastest growing mode, previously not understood, are derived analytically. The simplicity of the model suggests that shear dynamo action may be a generic property of sheared magnetohydrodynamic turbulence.

A numerical study of magnetic reconnection in the large-Lundquist-number ($S$), plasmoid-dominated regime is carried out for $S$ up to $10^7$. The theoretical model of Uzdensky {\it et al.} [Phys. Rev. Lett. {\bf 105}, 235002 (2010)] is confirmed and partially amended. The normalized reconnection rate is $\normEeff\sim 0.02$ independently of $S$ for $S\gg10^4$. The plasmoid flux ($\Psi$) and half-width ($w_x$) distribution functions scale as $f(\Psi)\sim \Psi^{-2}$ and $f(w_x)\sim w_x^{-2}$. The joint distribution of $\Psi$ and $w_x$ shows that plasmoids populate a triangular region $w_x\gtrsim\Psi/B_0$, where $B_0$ is the reconnecting field. It is argued that this feature is due to plasmoid coalescence. Macroscopic "monster" plasmoids with $w_x\sim 10$% of the system size are shown to emerge in just a few Alfv\'en times, independently of $S$, suggesting that large disruptive events are an inevitable feature of large-$S$ reconnection.

We investigate the anisotropy of Alfv\'enic turbulence in the inertial range of slow solar wind and in both driven and decaying reduced magnetohydrodynamic simulations. A direct comparison is made by measuring the anisotropic second-order structure functions in both data sets. In the solar wind, the perpendicular spectral index of the magnetic field is close to -5/3. In the forced simulation, it is close to -5/3 for the velocity and -3/2 for the magnetic field. In the decaying simulation, it is -5/3 for both fields. The spectral index becomes steeper at small angles to the local magnetic field direction in all cases. We also show that when using the global rather than local mean field, the anisotropic scaling of the simulations cannot always be properly measured.

It is proposed that critical balance - a scale-by-scale balance between the linear propagation and nonlinear interaction time scales - can be used as a universal scaling conjecture for determining the spectra of strong turbulence in anisotropic wave systems. Magnetohydrodynamic (MHD), rotating and stratified turbulence are considered under this assumption and, in particular, a novel and experimentally testable energy cascade scenario and a set of scalings of the spectra are proposed for low-Rossby-number rotating turbulence. It is argued that in neutral fluids, the critically balanced anisotropic cascade provides a natural path from strong anisotropy at large scales to isotropic Kolmogorov turbulence at very small scales. It is also argued that the kperp^{-2} spectra seen in recent numerical simulations of low-Rossby-number rotating turbulence may be analogous to the kperp^{-3/2} spectra of the numerical MHD turbulence in the sense that they could be explained by assuming that fluctuations are polarised (aligned) approximately as inertial waves (Alfven waves for MHD).

Weakly collisional plasmas dynamically develop pressure anisotropies with respect to the magnetic field. These anisotropies trigger plasma instabilities at scales just above the ion Larmor radius \rho_i and much below the mean free path \lambda_{mfp}. They have growth rates of a fraction of the ion cyclotron frequency - much faster than either the global dynamics or local turbulence. The instabilities dramatically modify the transport properties and, therefore, the macroscopic dynamics of the plasma. Their nonlinear evolution drives pressure anisotropies towards marginal stability, controlled by the plasma beta \beta_i. Here this nonlinear evolution is worked out for the simplest analytically tractable example - the parallel firehose instability. In the nonlinear regime, both analytical theory and the numerical solution predict secular growth of magnetic fluctuations. They develop a k^{-3} spectrum, extending from scales somewhat larger than \rho_i to the maximum scale that grows secularly with time (~t^{1/2}); the relative pressure anisotropy (\pperp-\ppar)/\ppar tends to the marginal value -2/\beta_i. The marginal state is achieved via changes in the magnetic field, not particle scattering. When a parallel ion heat flux is present, the firehose mutates into the new gyrothermal instability (GTI), which continues to exist up to firehose-stable values of pressure anisotropy, which can be positive and are limited by the heat flux. The nonlinear evolution of the GTI also features secular growth of magnetic fluctuations, but the spectrum is eventually dominated by modes around the scale ~\rho_i l_T/\lambda_{mfp}, where l_T is the scale of the parallel temperature variation. Implications for momentum and heat transport are speculated about. This study is motivated by the dynamics of galaxy cluster plasmas.

The first three-dimensional, nonlinear gyrokinetic simulation of plasma turbulence resolving scales from the ion to electron gyroradius with a realistic mass ratio is presented, where all damping is provided by resolved physical mechanisms. The resulting energy spectra are quantitatively consistent with a magnetic power spectrum scaling of $k^{-2.8}$ as observed in \emph{in situ} spacecraft measurements of the "dissipation range" of solar wind turbulence. Despite the strongly nonlinear nature of the turbulence, the linear kinetic \Alfven wave mode quantitatively describes the polarization of the turbulent fluctuations. The collisional ion heating is measured at sub-ion-Larmor radius scales, which provides the first evidence of the ion entropy cascade in an electromagnetic turbulence simulation.

We present the first measurement of the scale-dependent power anisotropy of Elsasser variables in imbalanced fast solar wind turbulence. The dominant Elsasser mode is isotropic at lower spacecraft frequencies but becomes increasingly anisotropic at higher frequencies. The sub-dominant mode is anisotropic throughout, but in a scale-independent way (at higher frequencies). There are two distinct subranges exhibiting different scalings within what is normally considered the inertial range. The low Alfven ratio and shallow scaling of the sub-dominant Elsasser mode suggest an interpretation of the observed discrepancy between the velocity and magnetic field scalings. The total energy is dominated by the latter. These results do not appear to be fully explained by any of the current theories of incompressible imbalanced MHD turbulence.

We consider the problem of self-regulated heating and cooling in galaxy clusters and the implications for cluster magnetic fields and turbulence. Viscous heating of a weakly collisional magnetised plasma is regulated by the pressure anisotropy with respect to the local direction of the magnetic field. The intracluster medium is a high-beta plasma, where pressure anisotropies caused by the turbulent stresses and the consequent local changes in the magnetic field will trigger very fast microscale instabilities. We argue that the net effect of these instabilities will be to pin the pressure anisotropies at a marginal level, controlled by the plasma beta parameter. This gives rise to local heating rates that turn out to be comparable to the radiative cooling rates. Furthermore, we show that a balance between this heating and Bremsstrahlung cooling is thermally stable, unlike the often conjectured balance between cooling and thermal conduction. Given a sufficient (and probably self-regulating) supply of turbulent power, this provides a physical mechanism for mitigating cooling flows and preventing cluster core collapse. For observed density and temperature profiles, the assumed balance of viscous heating and radiative cooling allows us to predict magnetic-field strengths, turbulent velocities and turbulence scales as functions of distance from the centre. Specific predictions and comparisons with observations are given for several different clusters. Our predictions can be further tested by future observations of cluster magnetic fields and turbulent velocities.

A conceptual model of resistive magnetic reconnection via a stochastic plasmoid chain is proposed. The global reconnection rate is shown to be independent of the Lundquist number. The distribution of fluxes in the plasmoids is shown to be an inverse square law. It is argued that there is a finite probability of emergence of abnormally large plasmoids, which can disrupt the chain (and may be responsible for observable large abrupt events in solar flares and sawtooth crashes). A criterion for the transition from magnetohydrodynamic to collisionless regime is provided.

Two-dimensional electrostatic turbulence in magnetized weakly-collisional plasmas exhibits a cascade of entropy in phase space [Phys. Rev. Lett. 103, 015003 (2009)]. At scales smaller than the gyroradius, this cascade is characterized by the dimensionless ratio D of the collision time to the eddy turnover time measured at the scale of the thermal Larmor radius. When D >> 1, a broad spectrum of fluctuations at sub-Larmor scales is found in both position and velocity space. The distribution function develops structure as a function of v_{perp}, the velocity coordinate perpendicular to the local magnetic field. The cascade shows a local-scale nonlinear interaction in both position and velocity spaces, and Kolmogorov's scaling theory can be extended into phase space.

The anisotropy of turbulence in the fast solar wind, between the ion and electron gyroscales, is directly observed using a multispacecraft analysis technique. Second order structure functions are calculated at different angles to the local magnetic field, for magnetic fluctuations both perpendicular and parallel to the mean field. In both components, the structure function value at large angles to the field S_perp is greater than at small angles S_par: in the perpendicular component S_perp/S_par = 5 +- 1 and in the parallel component S_perp/S_par > 3, implying spatially anisotropic fluctuations, k_perp > k_par. The spectral index of the perpendicular component is -2.6 at large angles and -3 at small angles, in broad agreement with critically balanced whistler and kinetic Alfven wave predictions. For the parallel component, however, it is shallower than -1.9, which is considerably less steep than predicted for a kinetic Alfven wave cascade.

Both global dynamics and turbulence in magnetized weakly collisional cosmic plasmas are described by general magnetofluid equations that contain pressure anisotropies and heat fluxes that must be calculated from microscopic plasma kinetic theory. It is shown that even without a detailed calculation of the pressure anisotropy or the heat fluxes, one finds the macroscale dynamics to be generically unstable to microscale Alfvenically polarized fluctuations. Two instabilities are considered in detail: the parallel firehose instability (including the finite-Larmor-radius effects that determine the fastest growing mode) and the gyrothermal instability (GTI). The latter is a new result - it is shown that a parallel ion heat flux destabilizes Alfvenically polarized fluctuations even in the absence of the negative pressure anisotropy required for the firehose. The main conclusion is that both pressure anisotropies and heat fluxes trigger plasma microinstabilities and, therefore, their values will likely be set by the nonlinear evolution of these instabilities. Ideas for understanding this nonlinear evolution are discussed. It is argued that cosmic plasmas will generically be "three-scale systems," comprising global dynamics, mesoscale turbulence and microscale plasma fluctuations. The astrophysical example of cool cores of galaxy clusters is considered and it is noted that observations point to turbulence in clusters being in a marginal state with respect to plasma microinstabilities and so it is the plasma microphysics that is likely to set the heating and conduction properties of the intracluster medium. In particular, a lower bound on the scale of temperature fluctuations implied by the GTI is derived.

We describe a technique for probing the statistical properties of cosmic magnetic fields based on radio polarimetry data. Second-order magnetic field statistics like the power spectrum cannot always distinguish between magnetic fields with essentially different spatial structure. Synchrotron polarimetry naturally allows certain 4th-order magnetic field statistics to be inferred from observational data, which lifts this degeneracy and can thereby help us gain a better picture of the structure of the cosmic fields and test theoretical scenarios describing magnetic turbulence. In this work we show that a 4th-order correlator of physical interest, the tension-force spectrum, can be recovered from the polarized synchrotron emission data. We develop an estimator for this quantity based on polarized-emission observations in the Faraday-rotation-free frequency regime. We consider two cases: a statistically isotropic field distribution, and a statistically isotropic field superimposed on a weak mean field. In both cases the tension force power spectrum is measurable; in the latter case, the magnetic power spectrum may also be obtainable. The method is exact in the idealized case of a homogeneous relativistic-electron distribution that has a power-law energy spectrum with a spectral index p=3, and assumes statistical isotropy of the turbulent field. We carry out tests of our method using synthetic data generated from numerically simulated magnetic fields. We show that the method is valid, that it is not prohibitively sensitive to the value of the electron spectral index, and that the observed tension-force spectrum allows one to distinguish between, e.g., a randomly tangled magnetic field (a default assumption in many studies) and a field organized in folded flux sheets or filaments.

A relationship is derived between power anisotropy and wavevector anisotropy in turbulent fluctuations. This can be used to interpret plasma turbulence measurements, for example in the solar wind. If fluctuations are anisotropic in shape then the ion gyroscale break point in spectra in the directions parallel and perpendicular to the magnetic field would not occur at the same frequency, and similarly for the electron gyroscale break point. This is an important consideration when interpreting solar wind observations in terms of anisotropic turbulence theories. Model magnetic field power spectra are presented assuming a cascade of critically balanced Alfven waves in the inertial range and kinetic Alfven waves in the dissipation range. The variation of power anisotropy with scale is compared to existing solar wind measurements and the similarities and differences are discussed.

We measure the power and spectral index anisotropy of high speed solar wind turbulence from scales larger than the outer scale down to the ion gyroscale, thus covering the entire inertial range. We show that the power and spectral indices at the outer scale of turbulence are approximately isotropic. The turbulent cascade causes the power anisotropy at smaller scales manifested by anisotropic scalings of the spectrum: close to k^{-5/3} across and k^{-2} along the local magnetic field, consistent with a critically balanced Alfvenic turbulence. By using data at different radial distances from the Sun, we show that the width of the inertial range does not change with heliocentric distance and explain this by calculating the radial dependence of the ratio of the outer scale to the ion gyroscale. At the smallest scales of the inertial range, close to the ion gyroscale, we find an enhancement of power parallel to the magnetic field direction coincident with a decrease in the perpendicular power. This is most likely related to energy injection by ion kinetic modes such as the firehose instability and also marks the beginning of the dissipation range of solar wind turbulence.

Two-dimensional numerical simulations of the effect of background turbulence on 2D resistive magnetic reconnection are presented. For sufficiently small values of the resistivity ($\eta$) and moderate values of the turbulent power ($\epsilon$), the reconnection rate is found to have a much weaker dependence on $\eta$ than the Sweet-Parker scaling of $\eta^{1/2}$ and is even consistent with an $\eta-$independent value. For a given value of $\eta$, the dependence of the reconnection rate on the turbulent power exhibits a critical threshold in $\epsilon$ above which the reconnection rate is significantly enhanced.

We present a theoretical framework for plasma turbulence in astrophysical plasmas (solar wind, interstellar medium, galaxy clusters, accretion disks). The key assumptions are that the turbulence is anisotropic with respect to the mean magnetic field and frequencies are low compared to the ion cyclotron frequency. The energy injected at the outer scale scale has to be converted into heat, which ultimately cannot be done without collisions. A KINETIC CASCADE develops that brings the energy to collisional scales both in space and velocity. Its nature depends on the physics of plasma fluctuations. In each of the physically distinct scale ranges, the kinetic problem is systematically reduced to a more tractable set of equations. In the "inertial range" above the ion gyroscale, the kinetic cascade splits into a cascade of Alfvenic fluctuations, which are governed by the RMHD equations at both the collisional and collisionless scales, and a passive cascade of compressive fluctuations, which obey a linear kinetic equation along the moving field lines associated with the Alfvenic component. In the "dissipation range" between the ion and electron gyroscales, there are again two cascades: the kinetic-Alfven-wave (KAW) cascade governed by two fluid-like Electron RMHD equations and a passive phase-space cascade of ion entropy fluctuations. The latter cascade brings the energy of the inertial-range fluctuations that was damped by collisionless wave-particle interaction at the ion gyroscale to collisional scales in the phase space and leads to ion heating. The KAW energy is similarly damped at the electron gyroscale and converted into electron heat. Kolmogorov-style scaling relations are derived for these cascades. Astrophysical and space-physical applications are discussed in detail.

A detailed numerical study of magnetic reconnection in resistive MHD for very large, previously inaccessible, Lundquist numbers ($10^4\le S\le 10^8$) is reported. Large-aspect-ratio Sweet-Parker current sheets are shown to be unstable to super-Alfv\'enically fast formation of plasmoid (magnetic-island) chains. The plasmoid number scales as $S^{3/8}$ and the instability growth rate in the linear stage as $S^{1/4}$, in agreement with the theory by Loureiro et al. [Phys. Plasmas {\bf 14}, 100703 (2007)]. In the nonlinear regime, plasmoids continue to grow faster than they are ejected and completely disrupt the reconnection layer. These results suggest that high-Lundquist-number reconnection is inherently time-dependent and hence call for a substantial revision of the standard Sweet-Parker quasi-stationary picture for $S>10^4$.

Howes et al. Reply to Comment on "Kinetic Simulations of Magnetized Turbulence in Astrophysical Plasmas" arXiv:0711.4355

This paper studies the turbulent cascade of magnetic energy in weakly collisional magnetized plasmas. A cascade model is presented, based on the assumptions of local nonlinear energy transfer in wavenumber space, critical balance between linear propagation and nonlinear interaction times, and the applicability of linear dissipation rates for the nonlinearly turbulent plasma. The model follows the nonlinear cascade of energy from the driving scale in the MHD regime, through the transition at the ion Larmor radius into the kinetic Alfven wave regime, in which the turbulence is dissipated by kinetic processes. The turbulent fluctuations remain at frequencies below the ion cyclotron frequency due to the strong anisotropy of the turbulent fluctuations, k_parallel << k_perp (implied by critical balance). In this limit, the turbulence is optimally described by gyrokinetics; it is shown that the gyrokinetic approximation is well satisfied for typical slow solar wind parameters. Wave phase velocity measurements are consistent with a kinetic Alfven wave cascade and not the onset of ion cyclotron damping. The conditions under which the gyrokinetic cascade reaches the ion cyclotron frequency are established. Cascade model solutions imply that collisionless damping provides a natural explanation for the observed range of spectral indices in the dissipation range of the solar wind. The dissipation range spectrum is predicted to be an exponential fall off; the power-law behavior apparent in observations may be an artifact of limited instrumental sensitivity. The cascade model is motivated by a programme of gyrokinetic simulations of turbulence and particle heating in the solar wind.

This paper describes a conceptual framework for understanding kinetic plasma turbulence as a generalized form of energy cascade in phase space. It is emphasized that conversion of turbulent energy into thermodynamic heat is only achievable in the presence of some (however small) degree of collisionality. The smallness of the collision rate is compensated by the emergence of small-scale structure in the velocity space. For gyrokinetic turbulence, a nonlinear perpendicular phase mixing mechanism is identified and described as a turbulent cascade of entropy fluctuations simultaneously occurring at spatial scales smaller than the ion gyroscale and in velocity space. Scaling relations for the resulting fluctuation spectra are derived. An estimate for the collisional cutoff is provided. The importance of adequately modeling and resolving collisions in gyrokinetic simulations is biefly discussed, as well as the relevance of these results to understanding the dissipation-range turbulence in the solar wind and the electrostatic microturbulence in fusion plasmas.

Numerical simulations of forced turbulence in elongated shearing boxes are carried out to demonstrate that a nonhelical turbulence in conjunction with a linear shear can give rise to a mean-field dynamo. Exponential growth of magnetic field at scales larger than the outer (forcing) scale of the turbulence is found. Over a range of values of the shearing rate S spanning approximately two orders of magnitude, the growth rate of the magnetic field is proportional to the imposed shear, gamma ~ S, while the characteristic spatial scale of the field is l_b ~ S^(-1/2). The effect is quite general: earlier results for the nonrotating case by Yousef et al. 2008 (PRL 100, 184501) are extended to shearing boxes with Keplerian rotation; it is also shown that the shear dynamo mechanism operates both below and above the threshold for the fluctuation dynamo. The apparently generic nature of the shear dynamo effect makes it an attractive object of study for the purpose of understanding the generation of magnetic fields in astrophysical systems.

The feasibility of a mean-field dynamo in nonhelical turbulence with superimposed linear shear is studied numerically in elongated shearing boxes. Exponential growth of magnetic field at scales much larger than the outer scale of the turbulence is found. The charateristic scale of the field is l_B ~ S^{-1/2} and growth rate is gamma ~ S, where S is the shearing rate. This newly discovered shear dynamo effect potentially represents a very generic mechanism for generating large-scale magnetic fields in a broad class of astrophysical systems with spatially coherent mean flows.

In turbulent high-beta astrophysical plasmas (exemplified by the galaxy cluster plasmas), pressure-anisotropy-driven firehose and mirror fluctuations grow nonlinearly to large amplitudes, dB/B ~ 1, on a timescale comparable to the turnover time of the turbulent motions. The principle of their nonlinear evolution is to generate secularly growing small-scale magnetic fluctuations that on average cancel the temporal change in the large-scale magnetic field responsible for the pressure anisotropies. The presence of small-scale magnetic fluctuations may dramatically affect the transport properties and, thereby, the large-scale dynamics of the high-beta astrophysical plasmas.

This letter presents the first ab initio, fully electromagnetic, kinetic simulations of magnetized turbulence in a homogeneous, weakly collisional plasma at the scale of the ion Larmor radius (ion gyroscale). Magnetic and electric-field energy spectra show a break at the ion gyroscale; the spectral slopes are consistent with scaling predictions for critically balanced turbulence of Alfven waves above the ion gyroscale (spectral index -5/3) and of kinetic Alfven waves below the ion gyroscale (spectral indices of -7/3 for magnetic and -1/3 for electric fluctuations). This behavior is also qualitatively consistent with in situ measurements of turbulence in the solar wind. Our findings support the hypothesis that the frequencies of turbulent fluctuations in the solar wind remain well below the ion cyclotron frequency both above and below the ion gyroscale.

Current sheets formed in magnetic reconnection events are found to be unstable to high-wavenumber perturbations. The instability is very fast: its maximum growth rate scales as S^{1/4} v_A/L, where L is the length of the sheet, v_A the Alfven speed and S the Lundquist number. As a result, a chain of plasmoids (secondary islands) is formed, whose number scales as S^{3/8}.

We present a cascade model for turbulence in weakly collisional plasmas that follows the nonlinear cascade of energy from the large scales of driving in the MHD regime to the small scales of the kinetic Alfven wave regime where the turbulence is dissipated by kinetic processes. Steady-state solutions of the model for the slow solar wind yield three conclusions: (1) beyond the observed break in the magnetic energy spectrum, one expects an exponential cut-off; (2) the widely held interpretation that this dissipation range obeys power-law behavior is an artifact of instrumental sensitivity limitations; and, (3) over the range of parameters relevant to the solar wind, the observed variation of dissipation range spectral indices from -2 to -4 is naturally explained by the varying effectiveness of Landau damping, from an undamped prediction of -7/3 to a strongly damped index around -4.

Direct numerical simulations of incompressible nonhelical randomly forced MHD turbulence are used to demonstrate for the first time that the fluctuation dynamo exists in the limit of large magnetic Reynolds number Rm>>1 and small magnetic Prandtl number Pm<<1. The dependence of the critical Rm_c for dynamo on the hydrodynamic Reynolds number Re is obtained for 1<Re<6700. In the limit Pm<<1, Rm_c is about three times larger than for the previously well established dynamo at large and moderate Prandtl numbers: Rm_c<200 for Re>6000 compared to Rm_c~60 for Pm>=1. Is is not as yet possible to determine numerically whether the growth rate of the magnetic energy is ~Rm^{1/2} in the limit Rm->infinity, as it should if the dynamo is driven by the inertial-range motions at the resistive scale.

We discuss the current understanding of the most basic properties of astrophysical MHD turbulence and trace the origins of the modern views and theoretical uncertainties to the ideas set forth in 1950s and 1960s by Iroshnikov, Kraichnan, Batchelor, Schlueter and Biermann. Universal aspects of the theory are emphasised. Two main astrophysical applications are touched upon: turbulence in the solar wind and in clusters of galaxies. These are, in a certain (very approximate) sense, two ``pure'' cases of small-scale turbulence, where theoretical models of the two main regimes of MHD turbulence -- with and without a strong mean magnetic field -- can be put to the test. They are also good examples of a complication that is more or less generic in astrophysical plasmas: the MHD description is, in fact, insufficient for astrophysical turbulence and plasma physics must make an entrance.

This paper is a detailed report on a programme of simulations used to settle a long-standing issue in the dynamo theory and demonstrate that the fluctuation dynamo exists in the limit of large magnetic Reynolds number Rm>>1 and small magnetic Prandtl number Pm<<1. The dependence of the critical Rm_c vs. the hydrodynamic Reynolds number Re is obtained for 1<Re<6700. In the limit Pm<<1, Rm_c is ~3 times larger than for Pm>1. The stability curve Rm_c(Re) (and, it is argued, the nature of the dynamo) is substantially different from the case of the simulations and liquid-metal experiments with a mean flow. It is not as yet possible to determine numerically whether the growth rate is ~Rm^{1/2} in the limit Re>>Rm>>1, as should be the case if the dynamo is driven by the inertial-range motions. The magnetic-energy spectrum in the low-Pm regime is qualitatively different from the Pm>1 case and appears to develop a negative spectral slope, although current resolutions are insufficient to determine its asymptotic form. At 1<Rm<Rm_c, the magnetic fluctuations induced via the tangling by turbulence of a weak mean field are investigated and the possibility of a k^{-1} spectrum above the resistive scale is examined. At low Rm<1, the induced fluctuations are well described by the quasistatic approximation; the k^{-11/3} spectrum is confirmed for the first time in direct numerical simulations.

This paper examines the consistency of the exact scaling laws for isotropic MHD turbulence in numerical simulations with large magnetic Prandtl numbers Pm and with Pm=1. The exact laws are used to elucidate the structure of the magnetic and velocity fields. Despite the linear scaling of certain third-order correlation functions, the situation is not analogous to the case of Kolmogorov turbulence. The magnetic field is adequately described by a model of stripy (folded) field with direction reversals at the resistive scale. At currently available resolutions, the cascade of kinetic energy is short-circuited by the direct exchange of energy between the forcing-scale motions and the stripy magnetic fields. This nonlocal interaction is the defining feature of isotropic MHD turbulence.

Observations of galaxy clusters show that the intracluster medium (ICM) is likely to be turbulent and is certainly magnetized. The properties of this magnetized turbulence are determined both by fundamental nonlinear magnetohydrodynamic interactions and by the plasma physics of the ICM, which has very low collisionality. Cluster plasma threaded by weak magnetic fields is subject to firehose and mirror instabilities. These saturate and produce fluctuations at the ion gyroscale, which can scatter particles, increasing the effective collision rate and, therefore, the effective Reynolds number of the ICM. A simple way to model this effect is proposed. The model yields a self-accelerating fluctuation dynamo whereby the field grows explosively fast, reaching the observed, dynamically important, field strength in a fraction of the cluster lifetime independent of the exact strength of the seed field. It is suggested that the saturated state of the cluster turbulence is a combination of the conventional isotropic magnetohydrodynamic turbulence, characterized by folded, direction-reversing magnetic fields and an Alfv\'en-wave cascade at collisionless scales. An argument is proposed to constrain the reversal scale of the folded field. The picture that emerges appears to be in qualitative agreement with observations of magnetic fields in clusters.

Theoretical approaches to low-frequency magnetized turbulence in collisionless and weakly collisional astrophysical plasmas are reviewed. The proper starting point for an analytical description of these plasmas is kinetic theory, not fluid equations. The anisotropy of the turbulence is used to systematically derive a series of reduced analytical models. Above the ion gyroscale, it is shown rigourously that the Alfven waves decouple from the electron-density and magnetic-field-strength fluctuations and satisfy the Reduced MHD equations. The density and field-strength fluctuations (slow waves and the entropy mode in the fluid limit), determined kinetically, are passively mixed by the Alfven waves. The resulting hybrid fluid-kinetic description of the low-frequency turbulence is valid independently of collisionality. Below the ion gyroscale, the turbulent cascade is partially converted into a cascade of kinetic Alfven waves, damped at the electron gyroscale. This cascade is described by a pair of fluid-like equations, which are a reduced version of the Electron MHD. The development of these theoretical models is motivated by observations of the turbulence in the solar wind and interstellar medium. In the latter case, the turbulence is spatially inhomogeneous and the anisotropic Alfvenic turbulence in the presence of a strong mean field may coexist with isotropic MHD turbulence that has no mean field.

Magnetohydrodynamic (MHD) turbulence is encountered in a wide variety of astrophysical plasmas, including accretion disks, the solar wind, and the interstellar and intracluster medium. On small scales, this turbulence is often expected to consist of highly anisotropic fluctuations with frequencies small compared to the ion cyclotron frequency. For a number of applications, the small scales are also collisionless, so a kinetic treatment of the turbulence is necessary. We show that this anisotropic turbulence is well described by a low frequency expansion of the kinetic theory called gyrokinetics. This paper is the first in a series to examine turbulent astrophysical plasmas in the gyrokinetic limit. We derive and explain the nonlinear gyrokinetic equations and explore the linear properties of gyrokinetics as a prelude to nonlinear simulations. The linear dispersion relation for gyrokinetics is obtained and its solutions are compared to those of hot-plasma kinetic theory. These results are used to validate the performance of the gyrokinetic simulation code {\tt GS2} in the parameter regimes relevant for astrophysical plasmas. New results on global energy conservation in gyrokinetics are also derived. We briefly outline several of the problems to be addressed by future nonlinear simulations, including particle heating by turbulence in hot accretion flows and in the solar wind, the magnetic and electric field power spectra in the solar wind, and the origin of small-scale density fluctuations in the interstellar medium.

The Planck mission will permit measurements of the polarization of the cosmic microwave background and of polarized foregrounds such as our own Galaxy with an unprecedented combination of accuracy and completeness. This will provide information on cosmological and galactic magnetic fields. The latter can be studied in detail via nearly Faraday-rotation free synchrotron and polarized dust emission. Methods are discussed to extract physically relevant information on the magnetic turbulence from Planck data and other measurements.

We propose a model of magnetic-field growth in galaxy clusters whereby the field is amplified by a factor of about 10^8 over a cosmologically short time of ~10^8 yr. Our model is based on the idea that the viscosity of the intracluster medium during the field-amplification epoch is determined not by particle collisions but by plasma microinstabilities: these give rise to small-scale fluctuations, which scatter particles, increasing their effective collision rate and, therefore, the effective Reynolds number. This gives rise to a bootstrap effect as the growth of the field triggers the instabilities which increase the Reynolds number which, in turn, accelerates the growth of the field. The growth is explosive and the result is that the observed field strength is reached over a fraction of the cluster lifetime independent of the exact strength of the seed field (which only needs to be above ~10^{-15} G to trigger the explosive growth).

We show that under very general conditions, cluster plasmas threaded by weak magnetic fields are subject to very fast growing plasma instabilities driven by the anisotropy of the plasma pressure (viscous stress) with respect to the local direction of the magnetic field. Such an anisotropy will naturally arise in any weakly magnetized plasma that has low collisionality and is subject to stirring. The magnetic field must be sufficiently weak for the instabilities to occur, viz., beta>Re^{1/2}. The instabilities are captured by the extended MHD model with Braginskii viscosity. However, their growth rates are proportional to the wavenumber down to the ion gyroscale, so MHD equations with Braginskii viscosity are not well posed and a fully kinetic treatment is necessary. The instabilities can lead to magnetic fields in clusters being amplified from seed strength of ~10^{-18} G to dynamically important strengths of ~10 microG on cosmologically trivial time scales (~10^8 yr). The fields produced during the amplification stage are at scales much smaller than observed. Predicting the saturated field scale and structure will require a kinetic theory of magnetized cluster turbulence.

We study the nonlinear evolution of the resistive tearing mode in slab geometry in two dimensions. We show that, in the strongly driven regime (large Delta'), a collapse of the X-point occurs once the island width exceeds a certain critical value ~1/Delta'. A current sheet is formed and the reconnection is exponential in time with a growth rate ~eta^1/2, where eta is the resistivity. If the aspect ratio of the current sheet is sufficiently large, the sheet can itself become tearing-mode unstable, giving rise to secondary islands, which then coalesce with the original island. The saturated state depends on the value of Delta'. For small Delta', the saturation amplitude is ~Delta' and quantitatively agrees with the theoretical prediction. If Delta' is large enough for the X-point collapse to have occured, the saturation amplitude increases noticeably and becomes independent of Delta'.

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

We report an extensive numerical study of the small-scale turbulent dynamo at large magnetic Prandtl numbers Pm. A Pm scan is given for the model case of low-Reynolds-number turbulence. We concentrate on three topics: magnetic-energy spectra and saturation levels, the structure of the field lines, and the field-strength distribution. The main results are (1) the folded structure (direction reversals at the resistive scale, field lines curved at the scale of the flow) persists from the kinematic to the nonlinear regime; (2) the field distribution is self-similar and appears to be lognormal during the kinematic regime and exponential in the saturated state; and (3) the bulk of the magnetic energy is at the resistive scale in the kinematic regime and remains there after saturation, although the spectrum becomes much shallower. We propose an analytical model of saturation based on the idea of partial two-dimensionalization of the velocity gradients with respect to the local direction of the magnetic folds. The model-predicted spectra are in excellent agreement with numerical results. Comparisons with large-Re, moderate-Pm runs are carried out to confirm the relevance of these results. New features at large Re are elongation of the folds in the nonlinear regime from the viscous scale to the box scale and the presence of an intermediate nonlinear stage of slower-than-exponential magnetic-energy growth accompanied by an increase of the resistive scale and partial suppression of the kinetic-energy spectrum in the inertial range. Numerical results for the saturated state do not support scale-by-scale equipartition between magnetic and kinetic energies, with a definite excess of magnetic energy at small scales. A physical picture of the saturated state is proposed.

We consider a solvable model of the decay of scalar variance in a single-scale random velocity field. We show that if there is a separation between the flow scale k_flow^{-1} and the box size k_box^{-1}, the decay rate lambda ~ (k_box/k_flow)^2 is determined by the turbulent diffusion of the box-scale mode. Exponential decay at the rate lambda is preceded by a transient powerlike decay (the total scalar variance ~ t^{-5/2} if the Corrsin invariant is zero, t^{-3/2} otherwise) that lasts a time t~1/\lambda. Spectra are sharply peaked at k=k_box. The box-scale peak acts as a slowly decaying source to a secondary peak at the flow scale. The variance spectrum at scales intermediate between the two peaks (k_box<<k<<k_flow) is ~ k + a k^2 + ... (a>0). The mixing of the flow-scale modes by the random flow produces, for the case of large Peclet number, a k^{-1+delta} spectrum at k>>k_flow, where delta ~ lambda is a small correction. Our solution thus elucidates the spectral make up of the ``strange mode,'' combining small-scale structure and a decay law set by the largest scales.

We consider the problem of incompressible, forced, nonhelical, homogeneous, isotropic MHD turbulence with no mean magnetic field. This problem is essentially different from the case with externally imposed uniform mean field. There is no scale-by-scale equipartition between magnetic and kinetic energies as would be the case for the Alfven-wave turbulence. The isotropic MHD turbulence is the end state of the turbulent dynamo which generates folded fields with small-scale direction reversals. We propose that the statistics seen in numerical simulations of isotropic MHD turbulence could be explained as a superposition of these folded fields and Alfven-like waves that propagate along the folds.

Cluster plasmas are magnetised already at very low magnetic field strength. Low collisionality implies that conservation of the first adiabatic invariant results in an anisotropic viscous stress (Braginskii viscosity) or, equivalently, anisotropic plasma pressure. This triggers firehose and mirror instabilities, which have growth rates proportional to the wavenumber down to scales of the order of ion Larmor radius. This means that MHD equations with Braginskii viscosity are not well posed and fully kinetic description is necessary. In this paper, we review the basic picture of small-scale dynamo in the cluster plasma and attempt to reconcile it with the existence of plasma instabilities at collisionless scales.

We consider the problem of incompressible, forced, nonhelical, homogeneous and isotropic MHD turbulence with no mean magnetic field and large magnetic Prandtl number. This type of MHD turbulence is the end state of the turbulent dynamo, which generates folded fields with small-scale direction reversals. We propose a model in which saturation is achieved as a result of the velocity statistics becoming anisotropic with respect to the local direction of the magnetic folds. The model combines the effects of weakened stretching and quasi-two-dimensional mixing and produces magnetic-energy spectra in remarkable agreement with numerical results at least in the case of a one-scale flow. We conjecture that the statistics seen in numerical simulations could be explained as a superposition of these folded fields and Alfv\'en-like waves that propagate along the folds.

The amplification of magnetic fields in a highly conducting fluid is studied numerically. During growth, the magnetic field is spatially intermittent: it does not uniformly fill the volume, but is concentrated in long thin folded structures. Contrary to a commonly held view, intermittency of the folded field does not increase indefinitely throughout the growth stage if diffusion is present. Instead, as we show, the probability-density function (PDF) of the field strength becomes self-similar. The normalized moments increase with magnetic Prandtl number in a powerlike fashion. We argue that the self-similarity is to be expected with a finite flow scale and system size. In the nonlinear saturated state, intermittency is reduced and the PDF is exponential. Parallels are noted with self-similar behavior recently observed for passive-scalar mixing and for map dynamos.

We report a series of numerical simulations showing that the critical magnetic Reynolds number Rm_c for the nonhelical small-scale dynamo depends on the Reynolds number Re. Namely, the dynamo is shut down if the magnetic Prandtl number Pr=Rm/Re is less than some critical value Pr_c<1 even for Rm for which dynamo exists at Pr>=1. We argue that, in the limit of Re->infinity, a finite Pr_c may exist. The second possibility is that Pr_c->0 as Re->infinity, while Rm_c tends to a very large constant value inaccessible at current resolutions. If there is a finite Pr_c, the dynamo is sustainable only if magnetic fields can exist at scales smaller than the flow scale, i.e., it is always effectively a large-Pr dynamo. If there is a finite Rm_c, our results provide a lower bound: Rm_c<220 for Pr<=1/8. This is larger than Rm in many planets and in all liquid-metal experiments.

The growth and saturation of magnetic field in conducting turbulent media with large magnetic Prandtl numbers are investigated. This regime is very common in low-density hot astrophysical plasmas. During the early (kinematic) stage, weak magnetic fluctuations grow exponentially and concentrate at the resistive scale, which lies far below the hydrodynamic viscous scale. The evolution becomes nonlinear when the magnetic energy is comparable to the kinetic energy of the viscous-scale eddies. A physical picture of the ensuing nonlinear evolution of the MHD dynamo is proposed. Phenomenological considerations are supplemented with a simple Fokker--Planck model of the nonlinear evolution of the magnetic-energy spectrum. It is found that, while the shift of the bulk of the magnetic energy from the subviscous scales to the velocity scales may be possible, it occurs very slowly -- at the resistive, rather than dynamical, time scale (for galaxies, this means that generation of large-scale magnetic fields cannot be explained by this mechanism). The role of Alfvenic motions and the implications for the fully developed isotropic MHD turbulence are discussed.

We study the intermittency and field-line structure of the MHD turbulence in plasmas with very large magnetic Prandtl numbers. In this regime, which is realized in the interstellar medium, some accretion disks, protogalaxies, galaxy-cluster gas, early Universe, etc., magnetic fluctuations can be excited at scales below the viscous cutoff. The salient feature of the resulting small-scale magnetic turbulence is the folded structure of the fields. It is characterized by very rapid transverse spatial oscillation of the field direction, while the field lines remain largely unbent up to the scale of the flow. Quantitatively, the fluctuation level and the field-line geometry can be studied in terms of the statistics of the field strength and of the field-line curvature. In the kinematic limit, the distribution of the field strength is an expanding lognormal, while that of the field-line curvature K is stationary and has a power tail K^{-13/7}. The field strength and curvature are anticorrelated, i.e. the growing fields are mostly flat, while the sharply curved fields remain relatively weak. The field, therefore, settles into a reduced-tension state. Numerical simulations demonstrate three essential features of the nonlinear regime. First, the total magnetic energy is equal to the total kinetic energy. Second, the intermittency is partially suppressed compared to the kinematic case, as the fields become more volume-filling and their distribution develops an exponential tail. Third, the folding structure of the field is unchanged from the kinematic case: the anticorrelation between the field strength and the curvature persists and the distribution of the latter retains the same power tail. We propose a model of back reaction based on the folding picture that reproduces all of the above numerical results.

This is a brief review of the main results of our recent studies of the nonlinear evolution of the small-scale MHD dynamo in the high-Prandtl-number regime and of the structure of the resulting saturated state of the isotropic homogeneous MHD turbulence. It is emphasized that the MHD regime without a uniform mean field (as is the case in our studies) is fundamentally different from the one in which such a field is externally imposed. The ability of the turbulence to bend and fold the magnetic-field lines leads to the emergence of a distinctive small-scale structure. The fields are organized in folds of characteristic length comparable to the size of the largest turbulent eddies with spatial-direction reversals at the resistive scale. These folds are very hard to destroy. In the nonlinear regime, the folding structure coexists with Alfven waves propagating along the folds. The turbulent energy injected by the forcing is dissipated in part resistively via the small-scale magnetic fields, and in part viscously via the Alfven waves.

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.

A weak fluctuating magnetic field embedded into a turbulent conducting medium grows exponentially while its characteristic scale decays. In the ISM and protogalactic plasmas, the magnetic Pr is very large, so a broad spectrum of growing magnetic fluctuations is excited at subviscous scales. We study the statistical correlations that are set up in the field pattern and show that the magnetic-field lines possess a folding structure, where most of the scale decrease is due to rapid transverse field direction reversals, while the scale of the field variation along itself stays approximately constant. Specifically, we find that the field strength and the field-line curvature are anticorrelated, and the curvature possesses a stationary limiting distribution with the bulk located at the values of curvature comparable to the characteristic wave number of the velocity field and a power tail extending to large values of curvature. The regions of large curvature, therefore, occupy only a small fraction of the total volume of the system. Our theoretical results are corroborated by direct numerical simulations. The implication of the folding effect is that the advent of the Lorentz back reaction occurs when the magnetic energy approaches that of the smallest turbulent eddies. Our results also directly apply to the problem of statistical geometry of the material lines in a random flow.

Most of the theoretical results on the kinematic amplification of small-scale magnetic fluctuations by turbulence have been confined to the model of white-noise-like advecting turbulent velocity field. In this work, the statistics of the passive magnetic field in the diffusion-free regime are considered for the case when the advecting flow is finite-time correlated. A new method is developed that allows one to systematically construct the correlation-time expansion for statistical characteristics of the field. The expansion is valid provided the velocity correlation time is smaller than the characteristic growth time of the magnetic fluctuations. This expansion is carried out up to first order in the general case of a d-dimensional arbitrarily compressible advecting flow. The growth rates for all moments of the magnetic field are derived. The effect of the first-order corrections is to reduce these growth rates. It is shown that introducing a finite correlation time leads to the loss of the small-scale statistical universality, which was present in the limit of the delta-correlated velocity field. Namely, the shape of the velocity time-correlation profile and the large-scale spatial structure of the flow become important. The latter is a new effect, that implies, in particular, that the approximation of a locally-linear shear flow does not fully capture the effect of nonvanishing correlation time. Physical applications of this theory include the small-scale kinematic dynamo in the interstellar medium and protogalactic plasmas.

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