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 pre
ssure-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.
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 describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for the study of astrophysical plasma dynamics. The code incorporates an energy-conserving particle integrator into a stable, second-order--accurate, three-stage predictor-predi
ctor-corrector integration algorithm. The constrained transport method is used to enforce the divergence-free constraint on the magnetic field. A delta-f scheme is included to facilitate a reduced-noise study of systems in which only small departures from an initial distribution function are anticipated. The effects of rotation and shear are implemented through the shearing-sheet formalism with orbital advection. These algorithms are embedded within an architecture similar to that used in the popular astrophysical magnetohydrodynamics code Athena, one that is modular, well-documented, easy to use, and efficiently parallelized for use on thousands of processors. We present a series of tests in one, two, and three spatial dimensions that demonstrate the fidelity and versatility of the code.
The intracluster medium of galaxy clusters is a weakly collisional, high-beta plasma in which the transport of heat and momentum occurs primarily along magnetic-field lines. Anisotropic heat conduction allows convective instabilities to be driven by
temperature gradients of either sign, the magnetothermal instability (MTI) in the outskirts of non-isothermal clusters and the heat-flux buoyancy-driven instability (HBI) in their cooling cores. We employ the Athena MHD code to investigate the nonlinear evolution of these instabilities, self-consistently including the effects of anisotropic viscosity (i.e. Braginskii pressure anisotropy), anisotropic conduction, and radiative cooling. We highlight the importance of the microscale instabilities that inevitably accompany and regulate the pressure anisotropies generated by the HBI and MTI. We find that, in all but the innermost regions of cool-core clusters, anisotropic viscosity significantly impairs the ability of the HBI to reorient magnetic-field lines orthogonal to the temperature gradient. Thus, while radio-mode feedback appears necessary in the central few tens of kpc, conduction may be capable of offsetting radiative losses throughout most of a cool core over a significant fraction of the Hubble time. Magnetically-aligned cold filaments are then able to form by local thermal instability. Viscous dissipation during the formation of a cold filament produces accompanying hot filaments, which can be searched for in deep Chandra observations of nearby cool-core clusters. In the case of the MTI, anisotropic viscosity maintains the coherence of magnetic-field lines over larger distances than in the inviscid case, providing a natural lower limit for the scale on which the field can fluctuate freely. In the nonlinear state, the magnetic field exhibits a folded structure in which the field-line curvature and field strength are anti-correlated.
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 r
espect 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.
