The Cassini mission entered the Grand Finale phase in April 2017 and executed 22.5 highly inclined, close-in orbits around Saturn before diving into the planet on September 15th 2017. Here we present our analysis of the Cassini Grand Finale magnetometer (MAG) dataset, focusing on Saturn's internal magnetic field. These measurements demonstrate that Saturn's internal magnetic field is exceptionally axisymmetric, with a dipole tilt less than 0.007 degrees (25.2 arcsecs). Saturn's magnetic equator was directly measured to be shifted northward by ~ 0.0468 +/- 0.00043 (1-sigma) $R_S$, 2820 +/- 26 km, at cylindrical radial distances between 1.034 and 1.069 $R_S$ from the spin-axis. Although almost perfectly axisymmetric, Saturn's internal magnetic field exhibits features on many characteristic length scales in the latitudinal direction. Examining Br at the a=0.75 $R_S$, c=0.6993 $R_S$ isobaric surface, the degrees 4 to 11 contributions correspond to latitudinally banded magnetic perturbations with characteristic width similar to that of the off-equatorial zonal jets observed in the atmosphere of Saturn. Saturn's internal magnetic field beyond 60 degrees latitude, in particular the small-scale features, are less well constrained by the available measurements, mainly due to incomplete spatial coverage in the polar region. A stably stratified layer thicker than 2500 km likely exists above Saturn's deep dynamo to filter out the non-axisymmetric internal magnetic field. A heat transport mechanism other than pure conduction, e.g. double diffusive convection, must be operating within this layer to be compatible with Saturn's observed luminosity. The latitudinally banded magnetic perturbations likely arise from a shallow secondary dynamo action with latitudinally banded differential rotation in the semi-conducting layer.

Saturn's magnetospheric magnetic field, planetary radio emissions, plasma populations and magnetospheric structure are all known to be modulated at periods close to the assumed rotation period of the planetary interior. These oscillations are readily apparent despite the high degree of axi-symmetry in the internally produced magnetic field of the planet, and have different rotation periods in the northern and southern hemispheres. In this paper we study the spatial structure of (near-) planetary period magnetic field oscillations in Saturn's equatorial magnetosphere. Extending previous analyses of these phenomena, we include all suitable data from the entire Cassini mission during its orbital tour of the planet, so as to be able to quantify both the amplitude and phase of these field oscillations throughout Saturn's equatorial plane, to distances of 30 planetary radii. We study the structure of these field oscillations in view of both independently rotating northern and southern systems, finding spatial variations in both magnetic fields and inferred currents flowing north-south that are common to both systems. With the greatly expanded coverage of the equatorial plane achieved during the latter years of the mission, we are able to present a complete survey of dawn-dusk and day-night asymmetries in the structure of the oscillating fields and currents. We show that the general structure of the rotating currents is simpler than previously reported, and that the relatively enhanced nightside equatorial fields and currents are due in part to related periodic vertical motion of Saturn's magnetotail current sheet.

We identify multiple periodic dusty structures in Saturn's Roche Division, a faint region spanning the $\sim3000$ km between the A and F rings. The locations and extent of these features vary over Cassini's tour of the Saturn system, being visible in 2006 and 2016-2017, but not in 2012-2014. These changes can be correlated with variations in Saturn's magnetospheric periods. In 2006 and 2016-2017, one of the drifting magnetospheric periods would produce a 3:4 resonance within the Roche Division, but in 2012-2014 these resonances would move into the A ring as the magnetospheric periods converged. A simple model of magnetic perturbations indicates that the magnetic field oscillations responsible for these structures have amplitudes of a few nanotesla, comparable to the magnetic field oscillation amplitudes of planetary period oscillations measured by the magnetometer onboard Cassini. However, some previously unnoticed features at higher radii have expected pattern speeds that are much slower than the magnetospheric periodicities. These structures may reflect an unexpectedly long-range propagation of resonant perturbations within dusty rings.

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.

The aurorae of Uranus were recently detected in the far ultraviolet with the Hubble Space Telescope (HST) providing a new, so far unique, means to remotely study the asymmetric Uranian magnetosphere from Earth. We analyze here two new HST Uranus campaigns executed in Sept. 2012 and Nov. 2014 with different temporal coverage and under variable solar wind conditions numerically predicted by three different MHD codes. Overall, HST images taken with the Space Telescope Imaging Spectrograph reveal auroral emissions in three pairs of successive images (one pair acquired in 2012 and two in 2014), hence six additional auroral detections in total, including the most intense Uranian aurorae ever seen with HST. The detected emissions occur close the expected arrival of interplanetary shocks. They appear as extended spots at southern latitudes, rotating with the planet. They radiate 5-24 kR and 1.3-8.8 GW of ultraviolet emission from H$_2$, last for tens of minutes and vary on timescales down to a few seconds. Fitting the 2014 observations with model auroral ovals constrains the longitude of the southern (northern) magnetic pole to 104+\-26{\deg} (284-\-26{\deg}) in the Uranian Longitude System. We suggest that the Uranian near-equinoctial aurorae are pulsed cusp emissions possibly triggered by large-scale magnetospheric compressions.

The magnetospheric cusps are important sites of the coupling of a magnetosphere with the solar wind. The combination of both ground- and space-based observations at Earth have enabled considerable progress to be made in understanding the terrestrial cusp and its role in the coupling of the magnetosphere to the solar wind via the polar magnetosphere. Voyager 2 fully explored Neptune's cusp in 1989 but highly inclined orbits of the Cassini spacecraft at Saturn present the most recent opportunity to repeatedly studying the polar magnetosphere of a rapidly rotating planet. In this paper we discuss observations made by Cassini during two passes through Saturn's southern polar magnetosphere. Our main findings are that i) Cassini directly encounters the southern polar cusp with evidence for the entry of magnetosheath plasma into the cusp via magnetopause reconnection, ii) magnetopause reconnection and entry of plasma into the cusp can occur over a range of solar wind conditions, and iii) double cusp morphologies are consistent with the position of the cusp oscillating in phase with Saturn's global magnetospheric periodicities.

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.

We present new observations of four closely-spaced NUV transits of the hot Jupiter-like exoplanet WASP-12b using HST/COS, significantly increasing the phase resolution of the observed NUV light curve relative to previous observations, while minimising the temporal variation of the system. We observe significant excess NUV absorption during the transit, with mean normalised in-transit fluxes of $F_\mathrm{norm}\simeq0.97$, i.e. $\simeq$2-5 $\sigma$ deeper than the optical transit level of $\simeq0.986$ for a uniform stellar disk (the exact confidence level depending on the normalisation method used). We further observe an asymmetric transit shape, such that the post-conjunction fluxes are overall $\simeq$2-3 $\sigma$ higher than pre-conjunction values, and characterised by rapid variations in count rate between the pre-conjunction and out of transit levels. We do not find evidence for an early ingress to the NUV transit as suggested by earlier HST observations. However, we show that the NUV count rate observed prior to the optical transit is highly variable, but overall $\simeq$2.2-3.0 $\sigma$ below the post-transit values and comparable in depth to the optical transit, possibly forming a variable region of NUV absorption from at least phase $\phi\simeq$0.83, limited by the data coverage.

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

We present an analysis of recent high spatial and spectral resolution ground-based infrared observations of H3+ spectra obtained with the 10-metre Keck II telescope in April 2011. We observed H3+ emission from Saturn's northern and southern auroral regions, simultaneously, over the course of more than two hours, obtaining spectral images along the central meridian as Saturn rotates. Previous ground-based work has derived only an average temperature of an individual polar region, summing an entire night of observations. Here we analyse 20 H3+ spectra, 10 for each hemisphere, providing H3+ temperature, column density and total emission in both the northern and southern polar regions simultaneously, improving on past results in temporal cadence and simultaneity. We find that: 1) the average thermospheric temperatures are 527+/-18 K in northern Spring and 583+/-13 K in southern Autumn, respectively; 2) this asymmetry in temperature is likely to be the result of an inversely proportional relationship between the total thermospheric heating rate (Joule heating and ion drag) and magnetic field strength - i.e. the larger northern field strength leads to reduced total heating rate and a reduced temperature, irrespective of season, and 3) this implies that thermospheric heating and temperatures are relatively insensitive to seasonal effects.

Saturn's ionosphere is produced when the otherwise neutral atmosphere is exposed to a flow of energetic charged particles or solar radiation. At low latitudes the latter should result in a weak planet-wide glow in infrared (IR), corresponding to the planet's uniform illumination by the Sun. The observed low-latitude ionospheric electron density is lower and the temperature higher than predicted by models. A planet-ring magnetic connection has been previously suggested in which an influx of water from the rings could explain the lower than expected electron densities in Saturn's atmosphere. Here we report the detection of a pattern of features, extending across a broad latitude band from ~25 to 60 degrees, that is superposed on the lower latitude background glow, with peaks in emission that map along the planet's magnetic field lines to gaps in Saturn's rings. This pattern implies the transfer of charged water products from the ring-plane to the ionosphere, revealing the influx on a global scale, flooding between 30 to 43% of the planet's upper-atmospheric surface. This ring `rain' plays a fundamental role in modulating ionospheric emissions and suppressing electron densities.

A number of ultra-cool dwarfs emit circularly polarised radio waves generated by the electron cyclotron maser instability. In the solar system such radio is emitted from regions of strong auroral magnetic field-aligned currents. We thus apply ideas developed for Jupiter's magnetosphere, being a well-studied rotationally-dominated analogue in our solar system, to the case of fast-rotating UCDs. We explain the properties of the radio emission from UCDs by showing that it would arise from the electric currents resulting from an angular velocity shear in the fast-rotating magnetic field and plasma, i.e. by an extremely powerful analogue of the process which causes Jupiter's auroras. Such a velocity gradient indicates that these bodies interact significantly with their space environment, resulting in intense auroral emissions. These results strongly suggest that auroras occur on bodies outside our solar system.

Planetary plasma and magnetic field environments can be studied by in situ measurements or by remote sensing. While the former provide precise information about plasma behaviour, instabilities and dynamics on local scales, the latter offers the global view necessary to understand the overall interaction of the magnetospheric plasma with the solar wind. Here we propose a novel and more elegant approach employing remote X-ray imaging techniques, which are now possible thanks to the relatively recent discovery of solar wind charge exchange X-ray emissions in the vicinity of the Earth's magnetosphere. We describe how an appropriately designed and located X-ray telescope, supported by simultaneous in situ measurements of the solar wind, can be used to image the dayside magnetosphere, magnetosheath and bow shock, with a temporal and spatial resolution sufficient to address several key outstanding questions concerning how the solar wind interacts with the Earth's magnetosphere on a global level. Our studies have led us to propose 'AXIOM: Advanced X-ray Imaging Of the Magnetosphere', a concept mission using a Vega launcher with a LISA Pathfinder-type Propulsion Module to place the spacecraft in a Lissajous orbit around the Earth - Moon L1 point. The model payload consists of an X-ray Wide Field Imager and an in situ plasma and magnetic field measurement package. This package comprises sensors designed to measure the bulk properties of the solar wind and to characterise its minor ion populations which cause charge exchange emission, and a magnetometer designed to measure the strength and direction of the solar wind magnetic field. We show simulations that demonstrate how the proposed X-ray telescope design is capable of imaging the predicted emission from the dayside magnetosphere with the sensitivity and cadence required to achieve the science goals of the mission.

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.

Plasma experiments in laboratory settings offer unique opportunities to address fundamental aspects of the solar dynamo and magnetism in the solar atmosphere. We argue here that ground-based laboratory experiments have direct connections to NASA based missions and NSF programs, and that a small investment in laboratory heliophysics may have a high payoff. We advocate for broad involvement in community-scale plasma experiments.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

Recent X-ray observations by the space mission Chandra confirmed the astonishing evidence for a diffuse, hot, thermal plasma at a temperature of 9. $10^7$ K (8 keV) found by previous surveys to extend over a few hundred parsecs in the Galactic Centre region. This plasma coexists with the usual components of the interstellar medium such as cold molecular clouds and a soft (~0.8 keV) component produced by supernova remnants, and its origin remains uncertain. First, simple calculations using a mean sound speed for a hydrogen-dominated plasma have suggested that it should not be gravitationally bound, and thus requires a huge energy source to heat it in less than the escape time. Second, an astrophysical mechanism must be found to generate such a high temperature. No known source has been identified to fulfill both requirements. Here we address the energetics problem and show that the hot component could actually be a gravitationally confined helium plasma. We illustrate the new prospects this opens by discussing the origin of this gas, and by suggesting possible heating mechanisms.

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

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.

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.

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

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.

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.

We simulate the evolution of an initially weak magnetic field in forced turbulence for a range of Prandtl numbers. The field grows exponentially with the Kulsrud-Anderson $k^{3/2}$ spectrum until the magnetic energy approaches the viscous-scale kinetic energy, where the magnetic forces then backreact on the velocity. Further growth proceeds more slowly until a saturated state is reached where the magnetic and kinetic energies are equal, and where the magnetic energy exists primarily at the resistive scale. We discuss the structure of this turbulence and the extrapolation of the results to astrophysically-large Prandtl numbers.

The effect of magnetic shear and shear flow on local gravitationally induced instabilities is investigated. A simple model is constructed allowing for an arbitrary entropy gradient and a shear plasma flow in the Boussinesq approximation. A transformation to shearing magnetic coordinates achieves a model with plasma flow along the magnetic field lines where the coordinate lines are coincident with the field lines. The solution for the normal modes of the system depends on two parameters: the Alfven Mach number of the plasma flow and the entropy gradient. The behavior of the unstable normal modes of this system is summarized by a stability diagram. Important characteristics of this stability diagram are the following: magnetic shear is stabilizing and the entropy gradient must exceed a threshold value for unstable mode growth to occur; flow acts to suppress mode growth in a substantially unstable regime as expected, yet near marginal stability it can lessen the stabilizing effect of magnetic shear and enhance the growth rates of the instability; and, as the Alfven Mach number approaches one, the instability is completely stabilized. Analytical work is presented supporting the characteristics of the stability diagram and illuminating the physical mechanisms controlling the behavior of the model. The implications of this work for astrophysical and fusion applications and the potential for future research extending the results to include compressibility are discussed.

Generation of magnetic field energy, without mean field generation, is studied. Isotropic mirror-symmetric turbulence of a conducting fluid amplifies the energy of small-scale magnetic perturbations if the magnetic Reynolds number is high, and the dimensionality of space d satisfies 2.103 < d <8.765. The result does not depend on the model of turbulence, incompressibility and isotropy being the only requirements.