Using hybrid-kinetic particle-in-cell simulation, we study the evolution of an expanding, collisionless, magnetized plasma in which strong Alfvenic turbulence is persistently driven. Temperature anisotropy generated adiabatically by the plasma expans
ion (and consequent decrease in the mean magnetic-field strength) gradually reduces the effective elasticity of the field lines, causing reductions in the linear frequency and residual energy of the Alfv{e}nic fluctuations. In response, these fluctuations modify their interactions and spatial anisotropy to maintain a scale-by-scale critical balance between their characteristic linear and nonlinear frequencies. Once the temperature anisotropy is sufficiently negative, the plasma becomes unstable to kinetic firehose instabilities, which excite rapidly growing magnetic fluctuations at ion-Larmor scales. The consequent pitch-angle scattering of particles maintains the temperature anisotropy near marginal stability, even as the turbulent plasma continues to expand. The resulting evolution of parallel and perpendicular temperatures does not satisfy double-adiabatic conservation laws, but is described accurately by a simple model that includes anomalous scattering. Our results have implications for understanding the complex interplay between macro- and micro-scale physics in various hot, dilute, astrophysical plasmas, and offer predictions concerning power spectra, residual energy, ion-Larmor-scale spectral breaks, and non-Maxwellian features in ion distribution functions that may be tested by measurements taken in high-beta regions of the solar wind.
Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $delta n/n_0 gtrsim 2/beta$ (where $betagg{1}$ is the ratio of thermal to magn
etic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more `fluid-like than kinetic. A nonlinear fluctuation-dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels, and weakly collisional thermodynamics when IAWs with $delta n/n_0 gtrsim 2/beta$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.
The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle--particle colli
sions allows departures from local thermodynamic equilibrium. These departures exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present a numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart. If instead the parallel viscous stress is left unabated -- a situation relevant to recent kinetic simulations of the fluctuation dynamo and to the early stages of the dynamo in a magnetized ICM -- the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics of the saturated state of the large-Prandtl-number (${rm Pm}gtrsim{1}$) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches magnetic-energy spectra. A prediction of this model, confirmed by our simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. In all cases that do exhibit a dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.
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.
To explain energy dissipation via turbulence in collisionless, magnetized plasmas, the existence of a dual real- and velocity-space cascade of ion-entropy fluctuations below the ion gyroradius has been proposed. Such a dual cascade, predicted by the
gyrokinetic theory, has previously been observed in gyrokinetic simulations of two-dimensional, electrostatic turbulence. For the first time we show evidence for a dual phase-space cascade of ion-entropy fluctuations in a three-dimensional simulation of hybrid-kinetic, electromagnetic turbulence. Some of the scalings observed in the energy spectra are consistent with a generalized theory for the cascade that accounts for the spectral anisotropy of critically balanced, intermittent, sub-ion-Larmor-scale fluctuations. The observed velocity-space cascade is also anisotropic with respect to the magnetic-field direction, with linear phase mixing along magnetic-field lines proceeding mainly at spatial scales above the ion gyroradius and nonlinear phase mixing across magnetic-field lines proceeding at perpendicular scales below the ion gyroradius. Such phase-space anisotropy could be sought in heliospheric and magnetospheric data of solar-wind turbulence and has far-reaching implications for the dissipation of turbulence in weakly collisional astrophysical plasmas.
The nonlinear state of a high-beta collisionless plasma is investigated when an imposed linear shear amplifies or diminishes a uniform magnetic field, driving pressure anisotropies and hence firehose/mirror instabilities. The evolution of the resulti
ng microscale turbulence is considered when the shear is switched off or reversed after one shear time (mimicking local behaviour of a macroscopic flow), so a new macroscale configuration is superimposed on the microscale state left behind by the previous one. There is a threshold value of plasma beta: when $betallOmega/S$ (ion cyclotron frequency/shear rate), the emergence of firehose/mirror fluctuations driven unstable by shear and their disappearance when the shear is removed/reversed are quasi-instantaneous compared to the shear time, viz., the decay time of these fluctuations is $simbeta/Omega ll 1/S$ (this result follows from the free decay of the fluctuations being constrained by the same marginal-stability thresholds as their growth). In contrast, when $betagtrsimOmega/S$ (ultra-high beta), the old microscale state can only be removed on the shear timescale. In this regime, driven firehose fluctuations grow secularly to order-unity amplitudes, compensating for the decay of the mean field and so pinning pressure anisotropy at marginal stability with no appreciable scattering of particles---which is unlike what happens at moderate $beta$. When the shear reverses, the shearing away of this firehose turbulence compensates for the increase in the mean field and thus prevents growth of the pressure anisotropy, stopping the system from going mirror-unstable. Therefore, at ultra-high beta, the system stays close to the firehose threshold, the mirror instability is largely suppressed, while the mean magnetic energy barely changes at all. Implications for plasma dynamo and thus the origin of cosmic magnetism are discussed.
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.
Cold fronts -- contact discontinuities in the intracluster medium (ICM) of galaxy clusters -- should be disrupted by Kelvin-Helmholtz (K-H) instabilities due to the associated shear velocity. However, many observed cold fronts appear stable. This ope
ns the possibility to place constraints on microphysical mechanisms that stabilize them, such as the ICM viscosity and/or magnetic fields. We performed exploratory high-resolution simulations of cold fronts arising from subsonic gas sloshing in cluster cores using the grid-based Athena MHD code, comparing the effects of isotropic Spitzer and anisotropic Braginskii viscosity (expected in a magnetized plasma). Magnetized simulations with full Braginskii viscosity or isotropic Spitzer viscosity reduced by a factor f ~ 0.1 are both in qualitative agreement with observations in terms of suppressing K-H instabilities. The RMS velocity of turbulence within the sloshing region is only modestly reduced by Braginskii viscosity. We also performed unmagnetized simulations with and without viscosity and find that magnetic fields have a substantial effect on the appearance of the cold fronts, even if the initial field is weak and the viscosity is the same. This suggests that determining the dominant suppression mechanism of a given cold front from X-ray observations (e.g. viscosity or magnetic fields) by comparison with simulations is not straightforward. Finally, we performed simulations including anisotropic thermal conduction, and find that including Braginskii viscosity in these simulations does not significant affect the evolution of cold fronts; they are rapidly smeared out by thermal conduction, as in the inviscid case.
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.
