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

We present a new algorithm for the discretization of the Vlasov-Maxwell system of equations for the study of plasmas in the kinetic regime. Using the discontinuous Galerkin finite element method for the spatial discretization, we obtain a high order accurate solution for the plasma's distribution function. Time stepping for the distribution function is done explicitly with a third order strong-stability preserving Runge-Kutta method. Since the Vlasov equation in the Vlasov-Maxwell system is a high dimensional transport equation, up to six dimensions plus time, we take special care to note various features we have implemented to reduce the cost while maintaining the integrity of the solution, including the use of a reduced high-order basis set. A series of benchmarks, from simple wave and shock calculations, to a five dimensional turbulence simulation, are presented to verify the efficacy of our set of numerical methods, as well as demonstrate the power of the implemented features.

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

We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihin et al., Astrophys. J. Suppl. 182:310 (2009)]. Two main applications of these equations are magnetised (Alfvenic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge-Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor-corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag-Tang-type decaying turbulence, both in fluid and kinetic regimes.

Results of the first validation of large guide field, $B_g / \delta B_0 \gg 1$, gyrokinetic simulations of magnetic reconnection at a fusion and solar corona relevant $\beta_i = 0.01$ and solar wind relevant $\beta_i = 1$ are presented, where $\delta B_0$ is the reconnecting field. Particle-in-cell (PIC) simulations scan a wide range of guide magnetic field strength to test for convergence to the gyrokinetic limit. The gyrokinetic simulations display a high degree of morphological symmetry, to which the PIC simulations converge when $\beta_i B_g / \delta B_0 \gtrsim 1$ and $B_g / \delta B_0 \gg 1$. In the regime of convergence, the reconnection rate, relative energy conversion, and overall magnitudes are found to match well between the PIC and gyrokinetic simulations, implying that gyrokinetics is capable of making accurate predictions well outside its regime of formal applicability. These results imply that in the large guide field limit many quantities resulting from the nonlinear evolution of reconnection scale linearly with the guide field.

A unique method of driving Alfvenic turbulence via an oscillating Langevin antenna is presented. This method of driving is motivated by a desire to inject energy into a finite domain numerical simulation in a manner that models the nonlinear transfer of energy from fluctuations in the turbulent cascade at scales larger than the simulation domain.. The oscillating Langevin antenna is shown to capture the essential features of the larger scale turbulence and efficiently couple to the plasma, generating steady-state turbulence within one characteristic turnaround time. The antenna is also sufficiently flexible to explore both strong and weak regimes of Alfvenic plasma turbulence.

This paper presents the numerical verification of an asymptotic analytical solution for the nonlinear interaction between counterpropagating Alfven waves, the fundamental building block of astrophysical plasma turbulence. The analytical solution, derived in the weak turbulence limit using the equations of incompressible MHD, is compared to a nonlinear gyrokinetic simulation of an Alfven wave collision. The agreement between these methods signifies that the incompressible solution satisfactorily describes the essential dynamics of the nonlinear energy transfer, even under the weakly collisional plasma conditions relevant to many astrophysical environments.

A refined cascade model for kinetic turbulence in weakly collisional astrophysical plasmas is presented that includes both the transition between weak and strong turbulence and the effect of nonlocal interactions on the nonlinear transfer of energy. The model describes the transition between weak and strong MHD turbulence and the complementary transition from strong kinetic Alfven wave (KAW) turbulence to weak dissipating KAW turbulence, a new regime of weak turbulence in which the effects of shearing by large scale motions and kinetic dissipation play an important role. The inclusion of the effect of nonlocal motions on the nonlinear energy cascade rate in the dissipation range, specifically the shearing by large-scale motions, is proposed to explain the nearly power-law energy spectra observed in the dissipation range of both kinetic numerical simulations and solar wind observations.

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.

The gyrokinetic simulation code AstroGK is developed to study fundamental aspects of kinetic plasmas and for applications mainly to astrophysical problems. AstroGK is an Eulerian slab code that solves the electromagnetic Gyrokinetic-Maxwell equations in five-dimensional phase space, and is derived from the existing gyrokinetics code GS2 by removing magnetic geometry effects. Algorithms used in the code are described. The code is benchmarked using linear and nonlinear problems. Serial and parallel performance scalings are also presented.

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.

Different results for the cascade power in strong, incompressible MHD turbulence with nonzero cross helicity appear in the literature. In this paper, we discuss the conditions under which these different results are valid. We define z+ to be the rms amplitude of Alfven waves propagating parallel to the background magnetic field, and z- to be the rms amplitude of Alfven waves propagating anti-parallel to the background magnetic field. Nonzero cross helicity implies that z+ and z- differ, and we take z- to be less than z+. We find that the mechanism that generates the z- fluctuations strongly affects the cascade power, because it controls the coherence time for interactions between oppositely directed wave packets at the outer scale. In particular, for fixed values of z+ and z-, the cascade power is in many cases larger when the z- fluctuations are generated by the reflection of z+ fluctuations than when the z- fluctuations are generated by forcing that is only weakly correlated with the z+ fluctuations.

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.

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.

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

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.

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.

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.

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.

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

Differential rotation occurs in conducting flows in accretion disks and planetary cores. In such systems, the magnetorotational instability can arise from coupling Lorentz and centrifugal forces to cause large radial angular momentum fluxes. We present the first experimental observation of the magnetorotational instability. Our system consists of liquid sodium between differentially rotating spheres, with an imposed coaxial magnetic field. We characterize the observed patterns, dynamics and torque increases, and establish that this instability can occur from a hydrodynamic turbulent background.

We consider the linear axisymmetric stability of a differentially rotating collisionless plasma in the presence of a weak magnetic field; we restrict our analysis to wavelengths much larger than the proton Larmor radius. This is the kinetic version of the magnetorotational instability explored extensively as mechanism for magnetic field amplification and angular momentum transport in accretion disks. The kinetic calculation is appropriate for hot accretion flows onto compact objects and for the growth of very weak magnetic fields, where the collisional mean free path is larger than the wavelength of the unstable modes. We show that the kinetic instability criterion is the same as in MHD, namely that the angular velocity decrease outwards. However, nearly every mode has a linear kinetic growth rate that differs from its MHD counterpart. The kinetic growth rates also depend explicitly on beta, i.e., on the ratio of the gas pressure to the pressure of the seed magnetic field. For beta ~ 1 the kinetic growth rates are similar to the MHD growth rates while for beta >> 1 they differ significantly. For beta >> 1, the fastest growing mode has a growth rate of sqrt{3} Omega for a Keplerian disk, larger than its MHD counterpart; there are also many modes whose growth rates are negligible, < beta^{-1/2} Omega << Omega. We provide a detailed physical interpretation of these results and show that gas pressure forces, rather than just magnetic forces, are central to the behavior of the magnetorotational instability in a collisionless plasma. We also discuss the astrophysical implications of our analysis.