Showing 1–49 of 49 results for author: Kunz, M W

Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

Authors: Evan L. Yerger, Matthew W. Kunz, Archie F. A. Bott, Anatoly Spitkovsky

Abstract: The regulation of electron heat transport in high-beta, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the… ▽ More The regulation of electron heat transport in high-beta, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as the inverse of electron beta. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of electron beta and temperature-gradient length scales under two different physical setups. The saturated heat flux in all of our simulations follows the expected inverse-electron-beta scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90-degree scattering problem in the cosmic ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Self-organization in collisionless, high-$β$ turbulence

Authors: S. Majeski, M. W. Kunz, J. Squire

Abstract: The MHD equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only po… ▽ More The MHD equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, $β\gg 1$ turbulence, focusing on the role of a self-organizational effect known as `magneto-immutability'. We describe this self-organization analytically through a high-$β$, reduced ordering of the CGL-MHD equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses, and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small ($\sim 5\%$) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-versus-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 42 pages, 13 figures

Cosmic-ray confinement in radio bubbles by micromirrors

Authors: Robert J. Ewart, Patrick Reichherzer, Archie F. A. Bott, Matthew W. Kunz, Alexander A. Schekochihin

Abstract: Radio bubbles, ubiquitous features of the intracluster medium around active galactic nuclei, are known to rise buoyantly for multiple scale heights through the intracluster medium (ICM). It is an open question how the bubbles can retain their high-energy cosmic-ray content over such distances. We propose that the enhanced scattering of cosmic rays due to micromirrors generated in the ICM, as propo… ▽ More Radio bubbles, ubiquitous features of the intracluster medium around active galactic nuclei, are known to rise buoyantly for multiple scale heights through the intracluster medium (ICM). It is an open question how the bubbles can retain their high-energy cosmic-ray content over such distances. We propose that the enhanced scattering of cosmic rays due to micromirrors generated in the ICM, as proposed recently by Reichherzer et al. (2023), is a viable mechanism for confining the cosmic rays within bubbles and can qualitatively reproduce their morphology. We discuss the observational implications of such a model of cosmic-ray confinement. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Comments: 10 pages, 6 figures

Critical magnetic Reynolds number of the turbulent dynamo in collisionless plasmas

Authors: Radhika Achikanath Chirakkara, Amit Seta, Christoph Federrath, Matthew W. Kunz

Abstract: The intracluster medium of galaxy clusters is an extremely hot and diffuse, nearly collisionless plasma, which hosts dynamically important magnetic fields of $\sim μ{\rm G}$ strength. Seed magnetic fields of much weaker strength of astrophysical or primordial origin can be present in the intracluster medium. In collisional plasmas, which can be approximated in the magneto-hydrodynamical (MHD) limi… ▽ More The intracluster medium of galaxy clusters is an extremely hot and diffuse, nearly collisionless plasma, which hosts dynamically important magnetic fields of $\sim μ{\rm G}$ strength. Seed magnetic fields of much weaker strength of astrophysical or primordial origin can be present in the intracluster medium. In collisional plasmas, which can be approximated in the magneto-hydrodynamical (MHD) limit, the turbulent dynamo mechanism can amplify weak seed fields to strong dynamical levels efficiently by converting turbulent kinetic energy into magnetic energy. However, the viability of this mechanism in weakly collisional or completely collisionless plasma is much less understood. In this study, we explore the properties of the collisionless turbulent dynamo by using three-dimensional hybrid-kinetic particle-in-cell simulations. We explore the properties of the collisionless turbulent dynamo in the kinematic regime for different values of the magnetic Reynolds number, ${\rm Rm}$, initial magnetic-to-kinetic energy ratio, $(E_{\rm{mag}}/E_{\rm{kin}})_{\rm{i}}$, and initial Larmor ratio, $(r_{\rm Larmor}/L_{\rm box})_{\rm i}$, i.e., the ratio of the Larmor radius to the size of the turbulent system. We find that in the `un-magnetised' regime, $(r_{\rm Larmor}/L_{\rm box})_{\rm i} > 1$, the critical magnetic Reynolds number for the dynamo action ${\rm Rm}_{\rm crit} \approx 107 \pm 3$. In the `magnetised' regime, $(r_{\rm Larmor}/L_{\rm box})_{\rm i} \lesssim 1$, we find a marginally higher ${\rm Rm}_{\rm crit} = 124 \pm 8$. We find that the growth rate of the magnetic energy does not depend on the strength of the seed magnetic field when the initial magnetisation is fixed. We also study the distribution and evolution of the pressure anisotropy in the collisionless plasma and compare our results with the MHD turbulent dynamo. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 17 pages, 16 figures

Journal ref: Monthly Notices of the Royal Astronomical Society (2023)

Local models of two-temperature accretion disc coronae. II. Ion thermal conduction and the absence of disc evaporation

Authors: Christopher J. Bambic, Eliot Quataert, Matthew W. Kunz, Yan-Fei Jiang

Abstract: We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray bin… ▽ More We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parameterized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc's surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: 14 pages, 7 figures, 1 table, 2 appendices, submitted to MNRAS

Tearing-mediated reconnection in magnetohydrodynamic poorly ionized plasmas. I. Onset and linear evolution

Authors: Elizabeth A. Tolman, Matthew W. Kunz, James M. Stone, Lev Arzamasskiy

Abstract: In high-Lundquist-number plasmas, reconnection proceeds via onset of tearing, followed by a nonlinear phase during which plasmoids continuously form, merge, and are ejected from the current sheet (CS). This process is understood in fully ionized, magnetohydrodynamic plasmas. However, many plasma environments, such as star-forming molecular clouds and the solar chromosphere, are poorly ionized. We… ▽ More In high-Lundquist-number plasmas, reconnection proceeds via onset of tearing, followed by a nonlinear phase during which plasmoids continuously form, merge, and are ejected from the current sheet (CS). This process is understood in fully ionized, magnetohydrodynamic plasmas. However, many plasma environments, such as star-forming molecular clouds and the solar chromosphere, are poorly ionized. We use theory and computation to study tearing-mediated reconnection in such poorly ionized systems. In this paper, we focus on the onset and linear evolution of this process. In poorly ionized plasmas, magnetic nulls on scales below $v_{\rm A,n0}/ν_{\rm ni0}$, with $v_{\rm A,n0}$ the neutral Alfvén speed and $ν_{\rm ni0}$ the neutral-ion collision frequency, will self-sharpen via ambipolar diffusion. This sharpening occurs at an increasing rate, inhibiting the onset of reconnection. Once the CS becomes thin enough, however, ions decouple from neutrals and thinning of the CS slows, allowing tearing to onset in a time of order $ν_{\rm ni0}^{-1}$. We find that the wavelength and growth rate of the mode that first disrupts the forming sheet can be predicted from a poorly ionized tearing dispersion relation; as the plasma recombination rate increases and ionization fraction decreases, the growth rate becomes an increasing multiple of $ν_{ni0}$ and the wavelength becomes a decreasing fraction of $v_{\rm A,n0}/ν_{\rm ni0}$. △ Less

Submitted 18 April, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 19 pages, 13 figures

Efficient micromirror confinement of sub-TeV cosmic rays in galaxy clusters

Authors: Patrick Reichherzer, Archie F. A. Bott, Robert J. Ewart, Gianluca Gregori, Philipp Kempski, Matthew W. Kunz, Alexander A. Schekochihin

Abstract: Recent observations suggest a stronger confinement of cosmic rays (CRs) in certain astrophysical systems than predicted by current CR-transport theories. We posit that the incorporation of microscale physics into CR-transport models can account for this enhanced CR confinement. We develop a theoretical description of the effect of magnetic microscale fluctuations originating from the mirror instab… ▽ More Recent observations suggest a stronger confinement of cosmic rays (CRs) in certain astrophysical systems than predicted by current CR-transport theories. We posit that the incorporation of microscale physics into CR-transport models can account for this enhanced CR confinement. We develop a theoretical description of the effect of magnetic microscale fluctuations originating from the mirror instability on macroscopic CR diffusion. We confirm our theory with large-dynamical-range simulations of CR transport in the intracluster medium (ICM) of galaxy clusters and kinetic simulations of CR transport in micromirror fields. We conclude that sub-TeV CR confinement in the ICM is far more effective than previously anticipated on the basis of Galactic-transport extrapolations. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: Utilizes PIC and MHD simulations, complemented by deep learning for data analysis. Currently under journal review. Comments welcome!

On hydromagnetic wave interactions in collisionless, high-$β$ plasmas

Authors: Stephen Majeski, Matthew W. Kunz

Abstract: We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high-$β$ collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew-Goldberger-Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to s… ▽ More We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high-$β$ collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew-Goldberger-Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in timescales of acoustic and Alfvénic fluctuations at high-$β$ means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed. △ Less

Submitted 14 December, 2023; v1 submitted 31 October, 2023; originally announced October 2023.

Comments: 26 pages, 7 figures

Electron-ion heating partition in imbalanced solar-wind turbulence

Authors: Jonathan Squire, Romain Meyrand, Matthew W. Kunz

Abstract: A likely candidate mechanism to heat the solar corona and solar wind is low-frequency "Alfvénic" turbulence sourced by magnetic fluctuations near the solar surface. Depending on its properties, such turbulence can heat different species via different mechanisms, and the comparison of theoretical predictions to observed temperatures, wind speeds, anisotropies, and their variation with heliocentric… ▽ More A likely candidate mechanism to heat the solar corona and solar wind is low-frequency "Alfvénic" turbulence sourced by magnetic fluctuations near the solar surface. Depending on its properties, such turbulence can heat different species via different mechanisms, and the comparison of theoretical predictions to observed temperatures, wind speeds, anisotropies, and their variation with heliocentric radius provides a sensitive test of this physics. Here we explore the importance of normalized cross helicity, or imbalance, for controlling solar-wind heating, since it is a key parameter of magnetized turbulence and varies systematically with wind speed and radius. Based on a hybrid-kinetic simulation in which the forcing's imbalance decreases with time -- a crude model for a plasma parcel entrained in the outflowing wind -- we demonstrate how significant changes to the turbulence and heating result from the "helicity barrier" effect. Its dissolution at low imbalance causes its characteristic features -- strong perpendicular ion heating with a steep "transition-range" drop in electromagnetic fluctuation spectra -- to disappear, driving more energy into electrons and parallel ion heat, and halting the emission of ion-scale waves. These predictions seem to agree with a diverse array of solar-wind observations, offering to explain a variety of complex correlations and features within a single theoretical framework. △ Less

Submitted 8 October, 2023; v1 submitted 24 August, 2023; originally announced August 2023.

Magnetogenesis in a collisionless plasma: from Weibel instability to turbulent dynamo

Authors: Muni Zhou, Vladimir Zhdankin, Matthew W. Kunz, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: We report on a first-principles numerical and theoretical study of plasma dynamo in a fully kinetic framework. By applying an external mechanical force to an initially unmagnetized plasma, we develop a self-consistent treatment of the generation of ``seed'' magnetic fields, the formation of turbulence, and the inductive amplification of fields by the fluctuation dynamo. Driven large-scale motions… ▽ More We report on a first-principles numerical and theoretical study of plasma dynamo in a fully kinetic framework. By applying an external mechanical force to an initially unmagnetized plasma, we develop a self-consistent treatment of the generation of ``seed'' magnetic fields, the formation of turbulence, and the inductive amplification of fields by the fluctuation dynamo. Driven large-scale motions in an unmagnetized, weakly collisional plasma are subject to strong phase mixing, which leads to the development of thermal pressure anisotropy. This anisotropy triggers the Weibel instability, which produces filamentary ``seed'' magnetic fields on plasma-kinetic scales. The plasma is thereby magnetized, enabling efficient stretching and folding of the fields by the plasma motions and the development of Larmor-scale kinetic instabilities such as the firehose and mirror. The scattering of particles off the associated microscale magnetic fluctuations provides an effective viscosity, regulating the field morphology and turbulence. During this process, the seed field is further amplified by the fluctuation dynamo until they reach energy equipartition with the turbulent flow. By demonstrating that equipartition magnetic fields can be generated from an initially unmagnetized plasma through large-scale turbulent flows, this work has important implications for the origin and amplification of magnetic fields in the intracluster and intergalactic mediums. △ Less

Submitted 28 July, 2023; originally announced August 2023.

Comments: 16 pages, 10 figures

Cosmic ray transport in large-amplitude turbulence with small-scale field reversals

Authors: Philipp Kempski, Drummond B. Fielding, Eliot Quataert, Alisa K. Galishnikova, Matthew W. Kunz, Alexander A. Philippov, Bart Ripperda

Abstract: The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current micro-physical CR transport models in magneto-hydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying t… ▽ More The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current micro-physical CR transport models in magneto-hydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying the alternative regime of large-amplitude turbulence with $δB/B_0 \gg 1$, in which intermittent small-scale magnetic field reversals are ubiquitous. We study particle transport in such turbulence by integrating trajectories in stationary snapshots. To quantify spatial diffusion, we use a setup with continuous particle injection and escape, which we term the turbulent leaky box. We find that particle transport is very different from the strong-guide-field case. Low-energy particles are better confined than high-energy particles, despite less efficient pitch-angle isotropization at small energies. In the limit of weak guide field, energy-dependent confinement is driven by the energy-dependent (in)ability to follow reversing magnetic field lines exactly and by the scattering in regions of ``resonant curvature", where the field line bends on a scale that is of order the local particle gyro-radius. We derive a heuristic model of particle transport in magnetic folds that approximately reproduces the energy dependence of transport found numerically. We speculate that CR propagation in the Galaxy is regulated by the intermittent field reversals highlighted here and discuss the implications of our findings for CR transport in the Milky Way. △ Less

Submitted 19 December, 2023; v1 submitted 24 April, 2023; originally announced April 2023.

Comments: Published in MNRAS, 16 pages, 9 Figures

Local models of two-temperature accretion disc coronae. I. Structure, outflows, and energetics

Authors: Christopher J. Bambic, Eliot Quataert, Matthew W. Kunz

Abstract: We use local stratified shearing-box simulations to elucidate the impact of two-temperature thermodynamics on the thermal structure of coronae in radiatively efficient accretion flows. Rather than treating the coronal plasma as an isothermal fluid, we use a simple, parameterized cooling function that models the collisional transfer of energy from the ions to the rapidly cooling leptons. Two-temper… ▽ More We use local stratified shearing-box simulations to elucidate the impact of two-temperature thermodynamics on the thermal structure of coronae in radiatively efficient accretion flows. Rather than treating the coronal plasma as an isothermal fluid, we use a simple, parameterized cooling function that models the collisional transfer of energy from the ions to the rapidly cooling leptons. Two-temperature models naturally form temperature inversions, with a hot, magnetically dominated corona surrounding a cold disc. Simulations with net vertical flux (NF) magnetic fields launch powerful magnetocentrifugal winds that would enhance accretion in a global system. The outflow rates are much better converged with increasing box height than analogous isothermal simulations, suggesting that the winds into two-temperature coronae may be sufficiently strong to evaporate a thin disc and form a radiatively inefficient accretion flow under some conditions. We find evidence for multiphase structure in the corona, with broad density and temperature distributions, and we propose criteria for the formation of a multiphase corona. The fraction of cooling in the surface layers of the disc is substantially larger for NF fields compared to zero net-flux configurations, with moderate NF simulations radiating ${\gtrsim}30$ per cent of the flow's total luminosity above two midplane scale-heights. Our work shows that NF fields may efficiently power the coronae of luminous Seyfert galaxies and quasars, providing compelling motivation for future studies of the heating mechanisms available to NF fields and the interplay of radiation with two-temperature thermodynamics. △ Less

Submitted 30 October, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: 20 pages, 11 figures, 1 table, 2 appendices, accepted for publication in MNRAS

Pressure anisotropy and viscous heating in weakly collisional plasma turbulence

Authors: Jonathan Squire, Matthew W Kunz, Lev Arzamasskiy, Zade Johnston, Eliot Quataert, Alexander A Schekochihin

Abstract: Pressure anisotropy can strongly influence the dynamics of weakly collisional, high-beta plasmas, but its effects are missed by standard magnetohydrodynamics (MHD). Small changes to the magnetic-field strength generate large pressure-anisotropy forces, heating the plasma, driving instabilities, and rearranging flows, even on scales far above the particles' gyroscales where kinetic effects are trad… ▽ More Pressure anisotropy can strongly influence the dynamics of weakly collisional, high-beta plasmas, but its effects are missed by standard magnetohydrodynamics (MHD). Small changes to the magnetic-field strength generate large pressure-anisotropy forces, heating the plasma, driving instabilities, and rearranging flows, even on scales far above the particles' gyroscales where kinetic effects are traditionally considered important. Here, we study the influence of pressure anisotropy on turbulent plasmas threaded by a mean magnetic field (Alfvénic turbulence). Extending previous results that were concerned with Braginskii MHD, we consider a wide range of regimes and parameters using a simplified fluid model based on drift kinetics with heat fluxes calculated using a Landau-fluid closure. We show that viscous (pressure-anisotropy) heating dissipates between a quarter and half of the turbulent cascade power injected at large scales; this does not depend strongly on either plasma beta or the ion-to-electron temperature ratio. This will in turn influence the plasma's thermodynamics by regulating energy partition between different dissipation channels (e.g., electron and ion heat). Due to the pressure anisotropy's rapid dynamical feedback onto the flows that create it -- an effect we term `magneto-immutability' -- the viscous heating is confined to a narrow range of scales near the forcing scale, supporting a nearly conservative, MHD-like inertial-range cascade, via which the rest of the energy is transferred to small scales. Despite the simplified model, our results -- including the viscous heating rate, distributions, and turbulent spectra -- compare favourably to recent hybrid-kinetic simulations. This is promising for the more general use of extended-fluid (or even MHD) approaches to model weakly collisional plasmas such as the intracluster medium, hot accretion flows, and the solar wind. △ Less

Submitted 6 July, 2023; v1 submitted 20 February, 2023; originally announced March 2023.

Comments: Accepted for publication in Journal of Plasma Physics

Microphysically modified magnetosonic modes in collisionless, high-$β$ plasmas

Authors: Stephen Majeski, Matthew W. Kunz, Jonathan Squire

Abstract: With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-$β$ collisionless plasmas, with particular attention to their excitation of, and reaction to, kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpen… ▽ More With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-$β$ collisionless plasmas, with particular attention to their excitation of, and reaction to, kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpendicular pressure over parallel pressure in regions where the plasma $β$ is increased. For mode amplitudes $δB/B_0 \gtrsim 0.3$, this excess excites the mirror instability. Particle scattering off these micro-scale mirrors frustrates the nonlinear saturation of transit-time damping, ensuring that large-amplitude NP modes continue their decay to small amplitudes. At asymptotically large wavelengths, we predict that the mirror-induced scattering will be large enough to interrupt transit-time damping entirely, isotropizing the pressure perturbations and morphing the collisionless NP mode into the magnetohydrodynamic (MHD) entropy mode. In fast waves, a fluctuating pressure anisotropy drives both mirror and firehose instabilities when the wave amplitude satisfies $δB/B_0 \gtrsim 2β^{-1}$. The induced particle scattering leads to delayed shock formation and MHD-like wave dynamics. Taken alongside prior work on self-interrupting Alfvén waves and self-sustaining ion-acoustic waves, our results establish a foundation for new theories of electromagnetic turbulence in low-collisionality, high-$β$ plasmas such as the intracluster medium, radiatively inefficient accretion flows, and the near-Earth solar wind. △ Less

Submitted 3 April, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 45 pages, 21 figures, accepted to Journal of Plasma Physics

Synchrotron Firehose Instability

Authors: Vladimir Zhdankin, Matthew W. Kunz, Dmitri A. Uzdensky

Abstract: We demonstrate using linear theory and particle-in-cell (PIC) simulations that a synchrotron-cooling collisionless plasma acquires pressure anisotropy and, if the plasma beta is sufficiently high, becomes unstable to the firehose instability, in a process that we dub the synchrotron firehose instability (SFHI). The SFHI channels free energy from the pressure anisotropy of the radiating, relativist… ▽ More We demonstrate using linear theory and particle-in-cell (PIC) simulations that a synchrotron-cooling collisionless plasma acquires pressure anisotropy and, if the plasma beta is sufficiently high, becomes unstable to the firehose instability, in a process that we dub the synchrotron firehose instability (SFHI). The SFHI channels free energy from the pressure anisotropy of the radiating, relativistic electrons (and/or positrons) into small-amplitude, kinetic-scale magnetic-field fluctuations, which pitch-angle scatter the particles and bring the plasma to a near-thermal state of marginal instability. The PIC simulations reveal a nonlinear cyclic evolution of firehose bursts interspersed by periods of stable cooling. We compare the SFHI for electron-positron and electron-ion plasmas. As a byproduct of the growing electron-firehose magnetic field fluctuations, magnetized ions gain a pressure anisotropy opposite to that of the electrons. If these ions are relativistically hot, we find that they also experience cooling due to collisionless thermal coupling with the electrons, which we argue is mediated by a secondary ion-cyclotron instability. We suggest that the SFHI may be activated in a number of astrophysical scenarios, such as within ejecta from black-hole accretion flows and relativistic jets, where the redistribution of energetic electrons from low to high pitch angles may cause transient bursts of radiation. △ Less

Submitted 30 October, 2022; originally announced October 2022.

Comments: 17 pages, 8 figures, submitted for publication

Kinetic Turbulence in Collisionless High-Beta Plasmas

Authors: Lev Arzamasskiy, Matthew W. Kunz, Jonathan Squire, Eliot Quataert, Alexander A. Schekochihin

Abstract: We present results from three-dimensional hybrid-kinetic simulations of Alfvénic turbulence in a high-beta, collisionless plasma. The key feature of such turbulence is the interplay between local wave--wave interactions between the fluctuations in the cascade and the non-local wave-particle interactions associated with kinetic micro-instabilities driven by anisotropy in the thermal pressure (namel… ▽ More We present results from three-dimensional hybrid-kinetic simulations of Alfvénic turbulence in a high-beta, collisionless plasma. The key feature of such turbulence is the interplay between local wave--wave interactions between the fluctuations in the cascade and the non-local wave-particle interactions associated with kinetic micro-instabilities driven by anisotropy in the thermal pressure (namely, firehose, mirror, and ion-cyclotron). We present theoretical estimates for, and calculate directly from the simulations, the effective collisionality and plasma viscosity in pressure-anisotropic high-beta turbulence, demonstrating that, for strong Alfvénic turbulence, the effective parallel-viscous scale is comparable to the driving scale of the cascade. Below this scale, the kinetic-energy spectrum indicates an Alfvénic cascade with a slope steeper than $-5/3$ due to the anisotropic viscous stress. The magnetic-energy spectrum is shallower than $-5/3$ near the ion-Larmor scale due to fluctuations produced by the firehose instability. Most of the cascade energy (80-90%) is dissipated as ion heating through a combination of Landau damping and anisotropic viscous heating. Our results have implications for models of particle heating in low-luminosity accretion onto supermassive black holes, the effective viscosity of the intracluster medium, and the interpretation of near-Earth solar-wind observations. △ Less

Submitted 2 March, 2023; v1 submitted 11 July, 2022; originally announced July 2022.

Comments: 34 pages, 19 figures, accepted to PRX

Turbulent regimes in collisions of 3D Alfvén-wave packets

Authors: Silvio Sergio Cerri, Thierry Passot, Dimitri Laveder, Pierre-Louis Sulem, Matthew W. Kunz

Abstract: Using 3D gyrofluid simulations, we revisit the problem of Alfven-wave (AW) collisions as building blocks of the Alfvenic cascade and their interplay with magnetic reconnection at magnetohydrodynamic (MHD) scales. Depending on the large-scale nonlinearity parameter $χ_0$ (the ratio between AW linear propagation time and nonlinear turnover time), different regimes are observed. For strong nonlineari… ▽ More Using 3D gyrofluid simulations, we revisit the problem of Alfven-wave (AW) collisions as building blocks of the Alfvenic cascade and their interplay with magnetic reconnection at magnetohydrodynamic (MHD) scales. Depending on the large-scale nonlinearity parameter $χ_0$ (the ratio between AW linear propagation time and nonlinear turnover time), different regimes are observed. For strong nonlinearities ($χ_0\sim1$), turbulence is consistent with a dynamically aligned, critically balanced cascade--fluctuations exhibit a scale-dependent alignment $\sinθ_k\propto k_\perp^{-1/4}$, a $k_\perp^{-3/2}$ spectrum and $k_\|\propto k_\perp^{1/2}$ spectral anisotropy. At weaker nonlinearities (small $χ_0$), a spectral break marking the transition between a large-scale weak regime and a small-scale $k_\perp^{-11/5}$ tearing-mediated range emerges, implying that dynamic alignment occurs also for weak nonlinearities. At $χ_0<1$ the alignment angle $θ_{k_\perp}$ shows a stronger scale dependence than in the $χ_0\sim1$ regime, i.e. $\sinθ_k\propto k_\perp^{-1/2}$ at $χ_0\sim0.5$, and $\sinθ_k\propto k_\perp^{-1}$ at $χ_0\sim0.1$. Dynamic alignment in the weak regime also modifies the large-scale spectrum, scaling roughly as $k_\perp^{-3/2}$ for $χ_0\sim0.5$ and as $k_\perp^{-1}$ for $χ_0\sim0.1$. A phenomenological theory of dynamically aligned turbulence at weak nonlinearities that can explain these spectra and the transition to the tearing-mediated regime is provided; at small $χ_0$, the strong scale dependence of the alignment angle combines with the increased lifetime of turbulent eddies to allow tearing to onset and mediate the cascade at scales that can be larger than those predicted for a critically balanced cascade by several orders of magnitude. Such a transition to tearing-mediated turbulence may even supplant the usual weak-to-strong transition. △ Less

Submitted 20 September, 2022; v1 submitted 9 July, 2022; originally announced July 2022.

Comments: 13 pages, 5 figures, 1 table; accepted for publication in The Astrophysical Journal

Plasma physics of the intracluster medium

Authors: Matthew W. Kunz, Thomas W. Jones, Irina Zhuravleva

Abstract: This Chapter provides a brief tutorial on some aspects of plasma physics that are fundamental to understanding the dynamics and energetics of the intracluster medium (ICM). The tutorial is split into two parts: one that focuses on the thermal plasma component -- its stability, viscosity, conductivity, and ability to amplify magnetic fields to dynamical strengths via turbulence and other plasma pro… ▽ More This Chapter provides a brief tutorial on some aspects of plasma physics that are fundamental to understanding the dynamics and energetics of the intracluster medium (ICM). The tutorial is split into two parts: one that focuses on the thermal plasma component -- its stability, viscosity, conductivity, and ability to amplify magnetic fields to dynamical strengths via turbulence and other plasma processes; and one that focuses on the non-thermal population of charged particles known as cosmic rays -- their acceleration, re-acceleration, and transport throughout the cluster volume. Observational context is woven throughout the narrative, from constraints on the strength and geometry of intracluster magnetic fields and the effective viscosity of the ICM, to examples of radio halos, radio relics, and cluster shocks that can test theories of particle acceleration. The promise of future X-ray missions to probe intracluster turbulence and discover the impact of small-scale plasma physics, coupled with sensitive, high-resolution radio observations of synchrotron-emitting plasma that reveal the properties of intracluster magnetic fields and particle-acceleration mechanisms, are likely to establish galaxy clusters as the premier cosmic laboratories for deciphering the fundamental physics of hot, dilute plasmas. △ Less

Submitted 5 May, 2022; originally announced May 2022.

Comments: This Chapter will appear in the Section "Galaxy Clusters" (Section Editors: E. Pointecouteau, E. Rasia, A. Simionescu) of the "Handbook of X-ray and Gamma-ray Astrophysics" (Editors in chief: C. Bambi and A. Santangelo)

Tearing instability and current-sheet disruption in the turbulent dynamo

Authors: Alisa K. Galishnikova, Matthew W. Kunz, Alexander A. Schekochihin

Abstract: Turbulence in a conducting plasma can amplify seed magnetic fields in what is known as the turbulent, or small-scale, dynamo. The associated growth rate and emergent magnetic-field geometry depend sensitively on the material properties of the plasma, in particular on the Reynolds number ${\rm Re}$, the magnetic Reynolds number ${\rm Rm}$, and their ratio ${\rm Pm}\equiv{\rm Rm}/{\rm Re}$. For… ▽ More Turbulence in a conducting plasma can amplify seed magnetic fields in what is known as the turbulent, or small-scale, dynamo. The associated growth rate and emergent magnetic-field geometry depend sensitively on the material properties of the plasma, in particular on the Reynolds number ${\rm Re}$, the magnetic Reynolds number ${\rm Rm}$, and their ratio ${\rm Pm}\equiv{\rm Rm}/{\rm Re}$. For ${\rm Pm} > 1$, the amplified magnetic field is gradually arranged into a folded structure, with direction reversals at the resistive scale and field lines curved at the larger scale of the flow. As the mean magnetic energy grows to come into approximate equipartition with the fluid motions, this folded structure is thought to persist. Using analytical theory and high-resolution MHD simulations with the Athena++ code, we show that these magnetic folds become unstable to tearing during the nonlinear stage of the dynamo for ${\rm Rm}\gtrsim 10^4$ and ${\rm Re}\gtrsim 10^3$. An ${\rm Rm}$- and ${\rm Pm}$-dependent tearing scale, at and below which folds are disrupted, is predicted theoretically and found to match well the characteristic field-reversal scale measured in the simulations. The disruption of folds by tearing increases the ratio of viscous-to-resistive dissipation. In the saturated state, the magnetic-energy spectrum exhibits a sub-tearing-scale steepening to a slope consistent with that predicted for tearing-mediated Alfvénic turbulence. Its spectral peak appears to be independent of the resistive scale and comparable to the driving scale of the flow, while the magnetic energy resides in a broad range of scales extending down to the field-reversal scale set by tearing. Emergence of a degree of large-scale magnetic coherence in the saturated state of the turbulent dynamo may be consistent with observations of magnetic-field fluctuations in galaxy clusters and recent laboratory experiments. △ Less

Submitted 13 October, 2022; v1 submitted 19 January, 2022; originally announced January 2022.

Comments: 15 pages, 11 figures, accepted for publication in PRX

Insensitivity of a turbulent laser-plasma dynamo to initial conditions

Authors: A. F. A. Bott, L. Chen, P. Tzeferacos, C. A. J. Palmer, A. R. Bell, R. Bingham, A. Birkel, D. H. Froula, J. Katz, M. W. Kunz, C. -K. Li, H-S. Park, R. Petrasso, J. S. Ross, B. Reville, D. Ryu, F. H. Séguin, T. G. White, A. A. Schekochihin, D. Q. Lamb, G. Gregori

Abstract: It has recently been demonstrated experimentally that a turbulent plasma created by the collision of two inhomogeneous, asymmetric, weakly magnetised laser-produced plasma jets can generate strong stochastic magnetic fields via the small-scale turbulent dynamo mechanism, provided the magnetic Reynolds number of the plasma is sufficiently large. In this paper, we compare such a plasma with one aris… ▽ More It has recently been demonstrated experimentally that a turbulent plasma created by the collision of two inhomogeneous, asymmetric, weakly magnetised laser-produced plasma jets can generate strong stochastic magnetic fields via the small-scale turbulent dynamo mechanism, provided the magnetic Reynolds number of the plasma is sufficiently large. In this paper, we compare such a plasma with one arising from two pre-magnetised plasma jets whose creation is identical save for the addition of a strong external magnetic field imposed by a pulsed magnetic field generator (`MIFEDS'). We investigate the differences between the two turbulent systems using a Thomson-scattering diagnostic, X-ray self-emission imaging and proton radiography. The Thomson-scattering spectra and X-ray images suggest that the presence of the external magnetic field has a limited effect on the plasma dynamics in the experiment. While the presence of the external magnetic field induces collimation of the flows in the colliding plasma jets and the initial strengths of the magnetic fields arising from the interaction between the colliding jets are significantly larger as a result of the external field, the energy and morphology of the stochastic magnetic fields post-amplification are indistinguishable. We conclude that, for turbulent laser-plasmas with super-critical magnetic Reynolds numbers, the dynamo-amplified magnetic fields are determined by the turbulent dynamics rather than the seed fields and modest changes in the initial flow dynamics of the plasma, a finding consistent with theoretical expectations and simulations of turbulent dynamos. △ Less

Submitted 5 January, 2022; originally announced January 2022.

Comments: 13 pages, 12 figures

Triggering tearing in a forming current sheet with the mirror instability

Authors: Himawan W. Winarto, Matthew W. Kunz

Abstract: We study the time-dependent formation and evolution of a current sheet (CS) in a magnetized, collisionless, high-beta plasma using hybrid-kinetic particle-in-cell simulations. An initially tearing-stable Harris sheet is frozen into a persistently driven incompressible flow so that its characteristic thickness gradually decreases in time. As the CS thins, the strength of the reconnecting field incr… ▽ More We study the time-dependent formation and evolution of a current sheet (CS) in a magnetized, collisionless, high-beta plasma using hybrid-kinetic particle-in-cell simulations. An initially tearing-stable Harris sheet is frozen into a persistently driven incompressible flow so that its characteristic thickness gradually decreases in time. As the CS thins, the strength of the reconnecting field increases, and adiabatic invariance in the inflowing fluid elements produces a field-biased pressure anisotropy with excess perpendicular pressure. At large plasma beta, this anisotropy excites the mirror instability, which deforms the reconnecting field on ion-Larmor scales and dramatically reduces the effective thickness of the CS. Tearing modes whose wavelengths are comparable to that of the mirrors then become unstable, triggering reconnection on smaller scales and at earlier times than would have occurred if the thinning CS were to have retained its Harris profile. A novel method for identifying and tracking X-points is introduced, yielding X-point separations that are initially intermediate between the perpendicular and parallel mirror wavelengths in the upstream plasma. These mirror-stimulated tearing modes ultimately grow and merge to produce island widths comparable to the CS thickness, an outcome we verify across a range of CS formation timescales and initial CS widths. Our results may find their most immediate application in the tearing disruption of magnetic folds generated by turbulent dynamo in weakly collisional, high-beta, astrophysical plasmas. △ Less

Submitted 14 February, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Spontaneous magnetization of collisionless plasma through the action of a shear flow

Authors: Muni Zhou, Vladimir Zhdankin, Matthew W. Kunz, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: We study in a fully kinetic framework the generation of seed magnetic fields through the Weibel instability driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and its saturation when the se… ▽ More We study in a fully kinetic framework the generation of seed magnetic fields through the Weibel instability driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and its saturation when the seed magnetic fields become strong enough to instigate gyromotion of particles and thereby inhibit their free-streaming. The predicted scaling dependencies of the saturated seed fields on key parameters (e.g., ratio of system scale to electron skin depth, the forcing amplitude) are confirmed by 3D and 2D particle-in-cell simulations using an electron-positron plasma. This work demonstrates the spontaneous magnetization of a collisionless plasma through large-scale motions as simple as a shear flow, and therefore has important implications for magnetogenesis in dilute astrophysical systems. △ Less

Submitted 3 October, 2021; originally announced October 2021.

Comments: 20 pages, 16 figures

High-frequency heating of the solar wind triggered by low-frequency turbulence

Authors: Jonathan Squire, Romain Meyrand, Matthew W. Kunz, Lev Arzamasskiy, Alexander A. Schekochihin, Eliot Quataert

Abstract: The fast solar wind's high speeds and nonthermal features require that significant heating occurs well above the Sun's surface. Two leading theories seem incompatible: low-frequency "Alfvénic" turbulence, which transports energy outwards and is observed ubiquitously by spacecraft but struggles to explain the observed dominance of ion over electron heating; and high-frequency ion-cyclotron waves (I… ▽ More The fast solar wind's high speeds and nonthermal features require that significant heating occurs well above the Sun's surface. Two leading theories seem incompatible: low-frequency "Alfvénic" turbulence, which transports energy outwards and is observed ubiquitously by spacecraft but struggles to explain the observed dominance of ion over electron heating; and high-frequency ion-cyclotron waves (ICWs), which explain the nonthermal heating of ions but lack an obvious source. Here, we argue that the recently proposed "helicity barrier" effect, which limits electron heating by inhibiting the turbulent cascade of energy to the smallest scales, can unify these two paradigms. Our six-dimensional simulations show how the helicity barrier causes the large-scale energy to grow in time, generating small parallel scales and high-frequency ICW heating from low-frequency turbulence. The resulting turbulence and ion distribution function also closely match in-situ measurements from Parker Solar Probe and other spacecraft, explaining, among other features, the decades-long puzzle of the steep "transition range" observed in magnetic fluctuation spectra. The theory predicts a causal link between plasma expansion and the ion-to-electron heating ratio. Given the observational association between wind speed and expansion, we argue that the helicity barrier could play a key role in regulating the bimodal speed distribution of the solar wind. △ Less

Submitted 23 February, 2022; v1 submitted 7 September, 2021; originally announced September 2021.

Adaptive critical balance and firehose instability in an expanding, turbulent, collisionless plasma

Authors: A. F. A. Bott, L. Arzamasskiy, M. W. Kunz, E. Quataert, J. Squire

Abstract: Using hybrid-kinetic particle-in-cell simulation, we study the evolution of an expanding, collisionless, magnetized plasma in which strong Alfvénic turbulence is persistently driven. Temperature anisotropy generated adiabatically by the plasma expansion (and consequent decrease in the mean magnetic-field strength) gradually reduces the effective elasticity of the field lines, causing reductions in… ▽ More Using hybrid-kinetic particle-in-cell simulation, we study the evolution of an expanding, collisionless, magnetized plasma in which strong Alfvénic turbulence is persistently driven. Temperature anisotropy generated adiabatically by the plasma expansion (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é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. Eventually 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. △ Less

Submitted 11 November, 2021; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: ApJL, in press. 9 pages, 5 figures

On stochastic heating and its phase-space signatures in low-$β$ kinetic turbulence

Authors: Silvio Sergio Cerri, Lev Arzamasskiy, Matthew W. Kunz

Abstract: We revisit the theory of stochastic heating of ions and investigate its phase-space signatures in kinetic turbulence of relevance to low-$β$ portions of the solar wind. We retain a full scale-dependent approach in our treatment, and consider the case in which electric-field fluctuations can be described by a generalized Ohm's law that includes Hall and thermo-electric effects. These two electric-f… ▽ More We revisit the theory of stochastic heating of ions and investigate its phase-space signatures in kinetic turbulence of relevance to low-$β$ portions of the solar wind. We retain a full scale-dependent approach in our treatment, and consider the case in which electric-field fluctuations can be described by a generalized Ohm's law that includes Hall and thermo-electric effects. These two electric-field terms provide the dominant contributions to stochastic ion heating when the ion-Larmor scale is much smaller than the ion skin depth, $ρ_{\mathrm{i}}\ll d_{\mathrm{i}}$, which is the case at $β{\ll}1$. Employing well-known spectral scaling laws for Alfvén-wave and kinetic-Alfvén-wave turbulent fluctuations, we obtain scaling relations characterizing the field-perpendicular particle-energization rate and energy diffusion coefficient associated with stochastic heating in these two regimes. Phase-space signatures of ion heating are then investigated using 3D hybrid-kinetic simulations of continuously driven Alfvénic turbulence at low $β$. In these simulations, energization of ions parallel to the magnetic field is sub-dominant compared to its perpendicular counterpart ($Q_{\parallel,\mathrm{i}}\ll Q_{\perp,\mathrm{i}}$), and the fraction of turbulent energy that goes into ion heating is ${\approx}75$\% at $β_{\mathrm{i}}=0.3$ and ${\approx}40$\% at $β_{\mathrm{i}}{\simeq}0.1$. The phase-space signatures of ion energization are consistent with Landau-resonant collisionless damping and a ($β$-dependent) combination of ion-cyclotron and stochastic heating. We demonstrate good agreement between our theory and various signatures associated with the stochastic portion of the heating. We discuss the effect of intermittency on stochastic heating and the implications of our work for the interpretation of stochastic heating in solar-wind spacecraft data. △ Less

Submitted 26 April, 2021; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: 23 pages, 9 figures (and 20 footnotes); accepted for publication in The Astrophysical Journal

Time-resolved fast turbulent dynamo in a laser plasma

Authors: A. F. A. Bott, P. Tzeferacos, L. Chen, C. A. J. Palmer, A. Rigby, A. Bell, R. Bingham, A. Birkel, C. Graziani, D. H. Froula, J. Katz, M. Koenig, M. W. Kunz, C. K. Li, J. Meinecke, F. Miniati, R. Petrasso, H. -S. Park, B. A. Remington, B. Reville, J. S. Ross, D. Ryu, D. Ryutov, F. Séguin, T. G. White , et al. (3 additional authors not shown)

Abstract: Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the universe. A theoretical framework attributing the origin and sustainment of these fields to the so-called fluctuation dynamo was recently validated by experiments on laser facilities in low-magnetic-Prandtl-number plasmas ($\mathrm{Pm} < 1$). However, the… ▽ More Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the universe. A theoretical framework attributing the origin and sustainment of these fields to the so-called fluctuation dynamo was recently validated by experiments on laser facilities in low-magnetic-Prandtl-number plasmas ($\mathrm{Pm} < 1$). However, the same framework proposes that the fluctuation dynamo should operate differently when $\mathrm{Pm} \gtrsim 1$, the regime relevant to many astrophysical environments such as the intracluster medium of galaxy clusters. This paper reports a new experiment that creates a laboratory $\mathrm{Pm} \gtrsim 1$ plasma dynamo for the first time. We provide a time-resolved characterization of the plasma's evolution, measuring temperatures, densities, flow velocities and magnetic fields, which allows us to explore various stages of the fluctuation dynamo's operation. The magnetic energy in structures with characteristic scales close to the driving scale of the stochastic motions is found to increase by almost three orders of magnitude from its initial value and saturate dynamically. It is shown that the growth of these fields occurs exponentially at a rate that is much greater than the turnover rate of the driving-scale stochastic motions. Our results point to the possibility that plasma turbulence produced by strong shear can generate fields more efficiently at the driving scale than anticipated by idealized MHD simulations of the nonhelical fluctuation dynamo; this finding could help explain the large-scale fields inferred from observations of astrophysical systems. △ Less

Submitted 24 July, 2020; originally announced July 2020.

Comments: 12 pages, 9 figures

Production and persistence of extreme two-temperature plasmas in radiative relativistic turbulence

Authors: Vladimir Zhdankin, Dmitri A. Uzdensky, Matthew W. Kunz

Abstract: Turbulence is a predominant process for energizing electrons and ions in collisionless astrophysical plasmas, and thus is responsible for shaping their radiative signatures (luminosity, spectra, and variability). To better understand the kinetic properties of a collisionless radiative plasma subject to externally driven turbulence, we investigate particle-in-cell simulations of relativistic plasma… ▽ More Turbulence is a predominant process for energizing electrons and ions in collisionless astrophysical plasmas, and thus is responsible for shaping their radiative signatures (luminosity, spectra, and variability). To better understand the kinetic properties of a collisionless radiative plasma subject to externally driven turbulence, we investigate particle-in-cell simulations of relativistic plasma turbulence with external inverse Compton cooling acting on the electrons. We find that ions continuously heat up while electrons gradually cool down (due to the net effect of radiation), and hence the ion-to-electron temperature ratio $T_i/T_e$ grows in time. We show that $T_i/T_e$ is limited only by the size and duration of the simulations (reaching $T_i/T_e \sim 10^3$), indicating that there are no efficient collisionless mechanisms of electron-ion thermal coupling. This result has implications for models of radiatively inefficient accretion flows, such as observed in the Galactic Center and in M87, for which so-called two-temperature plasmas with $T_i/T_e \gg 1$ have been invoked to explain their low luminosity. Additionally, we find that electrons acquire a quasi-thermal distribution (dictated by the competition of turbulent particle energization and radiative cooling), while ions undergo efficient nonthermal acceleration (acquiring a harder distribution than in equivalent non-radiative simulations). There is a modest nonthermal population of high-energy electrons that are beamed intermittently in space, time, and direction; these beamed electrons may explain rapid flares in certain high-energy astrophysical systems (e.g., in the Galactic Center). These numerical results demonstrate that extreme two-temperature plasmas can be produced and maintained by relativistic radiative turbulence. △ Less

Submitted 23 July, 2020; originally announced July 2020.

Comments: 8 pages, 5 figures, submitted for publication

Self-sustaining sound in collisionless, high-beta plasma

Authors: M. W. Kunz, J. Squire, A. A. Schekochihin, E. Quataert

Abstract: Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0 \gtrsim 2/β$ (where $β\gg{1}$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly… ▽ More Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0 \gtrsim 2/β$ (where $β\gg{1}$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves' otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more `fluid-like' than kinetic. A nonlinear fluctuation-dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels, and weakly collisional thermodynamics when IAWs with $δn/n_0 \gtrsim 2/β$ 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. △ Less

Submitted 17 September, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: 25 pages, 10 figures, Journal of Plasma Physics, in press

Fluctuation dynamo in a weakly collisional plasma

Authors: D. A. St-Onge, M. W. Kunz, J. Squire, A. A. Schekochihin

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

Submitted 13 July, 2020; v1 submitted 21 March, 2020; originally announced March 2020.

Comments: 62 pages, 25 figures, 1 table, accepted for publication to Journal of Plasma Physics

Journal ref: J. Plasma Phys. 86 (2020) 905860503

[Plasma 2020 Decadal] The Material Properties of Weakly Collisional, High-Beta Plasmas

Authors: M. W. Kunz, J. Squire, S. A. Balbus, S. D. Bale, C. H. K. Chen, E. Churazov, S. C. Cowley, C. B. Forest, C. F. Gammie, E. Quataert, C. S. Reynolds, A. A. Schekochihin, L. Sironi, A. Spitkovsky, J. M. Stone, I. Zhuravleva, E. G. Zweibel

Abstract: 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:… ▽ More 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. △ Less

Submitted 10 March, 2019; originally announced March 2019.

Comments: 6 pages, submitted for the Plasma 2020 Decadal Survey

Hybrid-Kinetic Simulations of Ion Heating in Alfvénic Turbulence

Authors: Lev Arzamasskiy, Matthew W. Kunz, Benjamin D. G. Chandran, Eliot Quataert

Abstract: We present three-dimensional, hybrid-kinetic numerical simulations of driven Alfvén-wave turbulence of relevance to the collisionless near-Earth solar wind. Special attention is paid to the spectral transition that occurs near the ion-Larmor scale and to the origins of preferential perpendicular ion heating and of non-thermal wings in the parallel distribution function. Several novel diagnostics a… ▽ More We present three-dimensional, hybrid-kinetic numerical simulations of driven Alfvén-wave turbulence of relevance to the collisionless near-Earth solar wind. Special attention is paid to the spectral transition that occurs near the ion-Larmor scale and to the origins of preferential perpendicular ion heating and of non-thermal wings in the parallel distribution function. Several novel diagnostics are used to show that the ion heating rate increases as the kinetic-Alfvén-wave fluctuations, which comprise the majority of the sub-ion-Larmor turbulent cascade, attain near-ion-cyclotron frequencies. We find that ${\approx}75$-$80\%$ of the cascade energy goes into heating the ions, broadly consistent with the near-Earth solar wind. This heating is accompanied by clear velocity-space signatures in the particle energization rates and the distribution functions, including a flattened core in the perpendicular-velocity distribution and non-Maxwellian wings in the parallel-velocity distribution. The latter are attributed to transit-time damping and the pitch-angle scattering of perpendicularly heated particles into the parallel direction. Accompanying these features is a steepening of the spectral index of sub-ion-Larmor magnetic-field fluctuations beyond the canonical $-2.8$, as field energy is transferred to thermal energy. These predictions may be tested by measurements in the near-Earth solar wind. △ Less

Submitted 17 May, 2019; v1 submitted 30 January, 2019; originally announced January 2019.

Comments: 14 pages, 10 figures, accepted to ApJ

Magneto-immutable turbulence in weakly collisional plasmas

Authors: Jonathan Squire, Alexander A. Schekochihin, Eliot Quataert, Matthew W. Kunz

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

Submitted 14 January, 2019; v1 submitted 29 November, 2018; originally announced November 2018.

Comments: Accepted for publication in J. Plasma Phys

Onset of magnetic reconnection in a collisionless, high-beta plasma

Authors: Andrew Alt, Matthew W. Kunz

Abstract: In a magnetized, collisionless plasma, the magnetic moment of the constituent particles is an adiabatic invariant. An increase in the magnetic-field strength in such a plasma thus leads to an increase in the thermal pressure perpendicular to the field lines. Above a $β$-dependent threshold (where $β$ is the ratio of thermal to magnetic pressure), this pressure anisotropy drives the mirror instabil… ▽ More In a magnetized, collisionless plasma, the magnetic moment of the constituent particles is an adiabatic invariant. An increase in the magnetic-field strength in such a plasma thus leads to an increase in the thermal pressure perpendicular to the field lines. Above a $β$-dependent threshold (where $β$ is the ratio of thermal to magnetic pressure), this pressure anisotropy drives the mirror instability, producing strong distortions in the field lines on ion-Larmor scales. The impact of this instability on magnetic reconnection is investigated using a simple analytical model for the formation of a current sheet (CS) and the associated production of pressure anisotropy. The difficulty in maintaining an isotropic, Maxwellian particle distribution during the formation and subsequent thinning of a CS in a collisionless plasma, coupled with the low threshold for the mirror instability in a high-$β$ plasma, imply that the geometry of reconnecting magnetic fields can differ radically from the standard Harris-sheet profile often used in simulations of collisionless reconnection. As a result, depending on the rate of CS formation and the initial CS thickness, tearing modes whose growth rates and wavenumbers are boosted by this difference may disrupt the mirror-infested CS before standard tearing modes can develop. A quantitative theory is developed to illustrate this process, which may find application in the tearing-mediated disruption of kinetic magnetorotational "channel" modes. △ Less

Submitted 6 January, 2019; v1 submitted 9 October, 2018; originally announced October 2018.

Comments: accepted for publication in Journal of Plasma Physics Letters on 4 Jan 2019

Electron-only magnetic reconnection in plasma turbulence

Authors: F. Califano, S. S. Cerri, M. Faganello, D. Laveder, M. Sisti, M. W. Kunz

Abstract: Hybrid-Vlasov-Maxwell simulations of magnetized plasma turbulence including non-linear electron-inertia effects in a generalized Ohm's law are presented. When fluctuation energy is injected on scales sufficiently close to ion-kinetic scales, the ions efficiently become de-magnetized and electron-scale current sheets largely dominate the distribution of the emerging current structures, in contrast… ▽ More Hybrid-Vlasov-Maxwell simulations of magnetized plasma turbulence including non-linear electron-inertia effects in a generalized Ohm's law are presented. When fluctuation energy is injected on scales sufficiently close to ion-kinetic scales, the ions efficiently become de-magnetized and electron-scale current sheets largely dominate the distribution of the emerging current structures, in contrast to the usual picture, where a full hierarchy of structure sizes is generally observed. These current sheets are shown to be the sites of electron-only reconnection (e-rec), in which the usual electron exhausts are unaccompanied by ion outflows and which are in qualitative agreement with those recently observed by MMS in the Earth's turbulent magnetosheath, downstream of the bow shock. Some features of the e-rec phenomenology are shown to be consistent with an electron magnetohydrodynamic description. Simulations suggest that this regime of collisionless reconnection may be found in turbulent systems where plasma processes, such as micro-instabilities and/or shocks, overpower the more customary turbulent cascade by directly injecting energy close to the ion-kinetic scales. △ Less

Submitted 16 July, 2020; v1 submitted 9 October, 2018; originally announced October 2018.

Comments: 20 pages, 8 figures; Accepted for publication

Fluctuation Dynamo in a Collisionless, Weakly Magnetized Plasma

Authors: Denis A. St-Onge, Matthew W. Kunz

Abstract: Results from a numerical study of fluctuation dynamo in a collisionless, weakly magnetized plasma are presented. The key difference between this dynamo and its magnetohydrodynamic (MHD) counterpart is the adiabatic production of magnetic-field-aligned pressure anisotropy by the amplification of a weak seed field. This in turn drives kinetic instabilities on the ion-Larmor scale---namely, firehose… ▽ More Results from a numerical study of fluctuation dynamo in a collisionless, weakly magnetized plasma are presented. The key difference between this dynamo and its magnetohydrodynamic (MHD) counterpart is the adiabatic production of magnetic-field-aligned pressure anisotropy by the amplification of a weak seed field. This in turn drives kinetic instabilities on the ion-Larmor scale---namely, firehose and mirror---which sever the adiabatic link between the thermal and magnetic pressures, thereby allowing the dynamo to proceed. After an initial phase of rapid growth driven by these instabilities, the magnetic energy grows exponentially and exhibits a $k^{3/2}$ spectrum that peaks near the resistive scale, similar to the large-magnetic-Prandtl-number ($\mathrm{Pm}\gg{1}$) MHD dynamo. The magnetic field self-organizes into a folded-sheet topology, with direction reversals at the resistive scale and field lines curved at the parallel scale of the flow. The effective $\mathrm{Pm}$ is determined by whether the ion-Larmor scale is above or below the field-reversing scale: in the former case, particles undergo Bohm-like diffusion; in the latter case, particles scatter primarily off firehose fluctuations residing at the ends of the magnetic folds, and the viscosity becomes anisotropic. The magnetic field ultimately saturates at dynamical strengths, with its spectral peak migrating towards larger scales. This feature, along with an anti-correlation of magnetic-field strength and field-line curvature and a gradual thinning of magnetic sheets into ribbons, resemble the saturated state of the large-$\mathrm{Pm}$ dynamo, the primary differences manifesting in firehose/mirror-unstable regions. These results have implications for magnetic-field growth in the weakly collisional intracluster medium of galaxy clusters. △ Less

Submitted 27 July, 2018; v1 submitted 28 June, 2018; originally announced June 2018.

Comments: 7 pages, 5 figures, submitted for publication to The Astrophysical Journal Letters

Dual phase-space cascades in 3D hybrid-Vlasov-Maxwell turbulence

Authors: S. S. Cerri, M. W. Kunz, F. Califano

Abstract: 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 e… ▽ More 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. △ Less

Submitted 8 March, 2018; v1 submitted 16 February, 2018; originally announced February 2018.

Comments: version accepted in ApJL

Astrophysical gyrokinetics: Turbulence in pressure-anisotropic plasmas at ion scales and beyond

Authors: Matthew W. Kunz, Ian G. Abel, Kristopher G. Klein, Alexander A. Schekochihin

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

Submitted 20 March, 2018; v1 submitted 6 December, 2017; originally announced December 2017.

Comments: 59 pages, 6 figures, accepted for publication in Journal of Plasma Physics (original 28 Nov 2017); abstract abridged

Pressure-anisotropy-induced nonlinearities in the kinetic magnetorotational instability

Authors: Jonathan Squire, Eliot Quataert, Matthew W. Kunz

Abstract: In collisionless and weakly collisional plasmas, such as hot accretion flows onto compact objects, the magnetorotational instability (MRI) can differ significantly from the standard (collisional) MRI. In particular, pressure anisotropy with respect to the local magnetic-field direction can both change the linear MRI dispersion relation and cause nonlinear modifications to the mode structure and gr… ▽ More In collisionless and weakly collisional plasmas, such as hot accretion flows onto compact objects, the magnetorotational instability (MRI) can differ significantly from the standard (collisional) MRI. In particular, pressure anisotropy with respect to the local magnetic-field direction can both change the linear MRI dispersion relation and cause nonlinear modifications to the mode structure and growth rate, even when the field and flow perturbations are small. This work studies these pressure-anisotropy-induced nonlinearities in the weakly nonlinear, high-ion-beta regime, before the MRI saturates into strong turbulence. Our goal is to better understand how the saturation of the MRI in a low collisionality plasma might differ from that in the collisional regime. We focus on two key effects: (i) the direct impact of self-induced pressure-anisotropy nonlinearities on the evolution of an MRI mode, and (ii) the influence of pressure anisotropy on the "parasitic instabilities" that are suspected to cause the mode to break up into turbulence. Our main conclusions are: (i) The mirror instability regulates the pressure anisotropy in such a way that the linear MRI in a collisionless plasma is an approximate nonlinear solution once the mode amplitude becomes larger than the background field (just as in MHD). This implies that differences between the collisionless and collisional MRI become unimportant at large amplitudes. (ii) The break up of large amplitude MRI modes into turbulence via parasitic instabilities is similar in collisionless and collisional plasmas. Together, these conclusions suggest that the route to magnetorotational turbulence in a collisionless plasma may well be similar to that in a collisional plasma, as suggested by recent kinetic simulations. As a supplement to these findings, we offer guidance for the design of future kinetic simulations of magnetorotational turbulence. △ Less

Submitted 13 November, 2017; v1 submitted 15 September, 2017; originally announced September 2017.

Comments: Submitted to Journal of Plasma Physics

Kinetic simulations of the interruption of large-amplitude shear-Alfvén waves in a high-beta plasma

Authors: J. Squire, M. W. Kunz, E. Quataert, A. A. Schekochihin

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

Submitted 13 December, 2017; v1 submitted 4 May, 2017; originally announced May 2017.

Journal ref: Phys. Rev. Lett. 119, 155101 (2017)

Magnetorotational Turbulence and Dynamo in a Collisionless Plasma

Authors: Matthew W. Kunz, James M. Stone, Eliot Quataert

Abstract: We present results from the first 3D kinetic numerical simulation of magnetorotational turbulence and dynamo, using the local shearing-box model of a collisionless accretion disc. The kinetic magnetorotational instability grows from a subthermal magnetic field having zero net flux over the computational domain to generate self-sustained turbulence and outward angular-momentum transport. Significan… ▽ More We present results from the first 3D kinetic numerical simulation of magnetorotational turbulence and dynamo, using the local shearing-box model of a collisionless accretion disc. The kinetic magnetorotational instability grows from a subthermal magnetic field having zero net flux over the computational domain to generate self-sustained turbulence and outward angular-momentum transport. Significant Maxwell and Reynolds stresses are accompanied by comparable viscous stresses produced by field-aligned ion pressure anisotropy, which is regulated primarily by the mirror and ion-cyclotron instabilities through particle trapping and pitch-angle scattering. The latter endow the plasma with an effective viscosity that is biased with respect to the magnetic-field direction and spatio-temporally variable. Energy spectra suggest an Alfvén-wave cascade at large scales and a kinetic-Alfvén-wave cascade at small scales, with strong small-scale density fluctuations and weak non-axisymmetric density waves. Ions undergo non-thermal particle acceleration, their distribution accurately described by a kappa distribution. These results have implications for the properties of low-collisionality accretion flows, such as that near the black hole at the Galactic center. △ Less

Submitted 19 October, 2016; v1 submitted 29 August, 2016; originally announced August 2016.

Comments: 6 pages, 6 figures, accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 117, 235101 (2016)

Linear Vlasov theory of a magnetised, thermally stratified atmosphere

Authors: Rui Xu, Matthew W. Kunz

Abstract: The stability of a collisionless, magnetised plasma to local convective disturbances is examined, with a focus on kinetic and finite-Larmor-radius effects. Specific application is made to the outskirts of galaxy clusters, which contain hot and tenuous plasma whose temperature increases in the direction of gravity. At long wavelengths (the "drift-kinetic" limit), we obtain the kinetic version of th… ▽ More The stability of a collisionless, magnetised plasma to local convective disturbances is examined, with a focus on kinetic and finite-Larmor-radius effects. Specific application is made to the outskirts of galaxy clusters, which contain hot and tenuous plasma whose temperature increases in the direction of gravity. At long wavelengths (the "drift-kinetic" limit), we obtain the kinetic version of the magnetothermal instability (MTI) and its Alfvénic counterpart (Alfvénic MTI), which were previously discovered and analysed using a magnetofluid (i.e. Braginskii) description. At sub-ion-Larmor scales, we discover an overstability driven by the electron temperature gradient of kinetic-Alfvén drift waves -- the electron MTI (eMTI) -- whose growth rate is even larger than the standard MTI. At intermediate scales, we find that ion finite-Larmor-radius effects tend to stabilise the plasma. We discuss the physical interpretation of these instabilities in detail, and compare them both with previous work on magnetised convection in a collisional plasma and with temperature-gradient-driven drift-wave instabilities well-known to the magnetic-confinement-fusion community. The implications of having both fluid and kinetic scales simultaneously driven unstable by the same temperature gradient are briefly discussed. △ Less

Submitted 18 August, 2016; originally announced August 2016.

Comments: 51 pages, 9 figures; to appear in Journal of Plasma Physics

Multi-Species Measurements of the Firehose and Mirror Instability Thresholds in the Solar Wind

Authors: C. H. K. Chen, L. Matteini, A. A. Schekochihin, M. L. Stevens, C. S. Salem, B. A. Maruca, M. W. Kunz, S. D. Bale

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

Submitted 8 June, 2016; originally announced June 2016.

Journal ref: Astrophys. J. Lett. 825 L26 (2016)

Suppression of thermal conduction in a mirror-unstable plasma

Authors: S. V. Komarov, E. M. Churazov, M. W. Kunz, A. A. Schekochihin

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

Submitted 3 June, 2016; v1 submitted 1 March, 2016; originally announced March 2016.

Comments: 11 pages, 11 figures

Pressure-anisotropy-driven microturbulence and magnetic-field evolution in shearing, collisionless plasma

Authors: S. Melville, A. A. Schekochihin, M. W. Kunz

Abstract: The nonlinear state of a high-beta collisionless plasma is investigated when an imposed linear shear amplifies or diminishes a uniform magnetic field, driving pressure anisotropies and hence firehose/mirror instabilities. The evolution of the resulting microscale turbulence is considered when the shear is switched off or reversed after one shear time (mimicking local behaviour of a macroscopic flo… ▽ More The nonlinear state of a high-beta collisionless plasma is investigated when an imposed linear shear amplifies or diminishes a uniform magnetic field, driving pressure anisotropies and hence firehose/mirror instabilities. The evolution of the resulting microscale turbulence is considered when the shear is switched off or reversed after one shear time (mimicking local behaviour of a macroscopic flow), so a new macroscale configuration is superimposed on the microscale state left behind by the previous one. There is a threshold value of plasma beta: when $β\llΩ/S$ (ion cyclotron frequency/shear rate), the emergence of firehose/mirror fluctuations driven unstable by shear and their disappearance when the shear is removed/reversed are quasi-instantaneous compared to the shear time, viz., the decay time of these fluctuations is $\simβ/Ω\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 $β\gtrsimΩ/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 $β$. 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. △ Less

Submitted 26 December, 2015; originally announced December 2015.

Comments: MNRAS-style latex, 21 pages, 37 figures, submitted to MNRAS

Inertial-Range Kinetic Turbulence in Pressure-Anisotropic Astrophysical Plasmas

Authors: M. W. Kunz, A. A. Schekochihin, C. H. K. Chen, I. G. Abel, S. C. Cowley

Abstract: 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 a… ▽ More 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. △ Less

Submitted 5 June, 2015; v1 submitted 27 January, 2015; originally announced January 2015.

Comments: 61 pages, accepted to Journal of Plasma Physics; Abstract abridged

Journal ref: J. Plasma Phys. 81 (2015) 325810501

Firehose and Mirror Instabilities in a Collisionless Shearing Plasma

Authors: Matthew W. Kunz, Alexander A. Schekochihin, James M. Stone

Abstract: 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 pinne… ▽ More 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é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. △ Less

Submitted 1 May, 2014; v1 submitted 31 January, 2014; originally announced February 2014.

Comments: 5 pages, 8 figures, accepted for publication in Physical Review Letters

Pegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics

Authors: Matthew W. Kunz, James M. Stone, Xue-Ning Bai

Abstract: 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-predictor-corrector integration algorithm. The constrained transport method is used to enforce the divergence-free constraint on the magnetic field. A delt… ▽ More 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-predictor-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. △ Less

Submitted 19 November, 2013; originally announced November 2013.

Comments: 27 pages, 12 figures, accepted for publication in Journal of Computational Physics

Magnetic self-organisation in Hall-dominated magnetorotational turbulence

Authors: Matthew W. Kunz, Geoffroy Lesur

Abstract: The magnetorotational instability (MRI) is the most promising mechanism by which angular momentum is efficiently transported outwards in astrophysical discs. However, its application to protoplanetary discs remains problematic. These discs are so poorly ionised that they may not support magnetorotational turbulence in regions referred to as `dead zones'. It has recently been suggested that the Hal… ▽ More The magnetorotational instability (MRI) is the most promising mechanism by which angular momentum is efficiently transported outwards in astrophysical discs. However, its application to protoplanetary discs remains problematic. These discs are so poorly ionised that they may not support magnetorotational turbulence in regions referred to as `dead zones'. It has recently been suggested that the Hall effect, a non-ideal magnetohydrodynamic (MHD) effect, could revive these dead zones by enhancing the magnetically active column density by an order of magnitude or more. We investigate this idea by performing local, three-dimensional, resistive Hall-MHD simulations of the MRI in situations where the Hall effect dominates over Ohmic dissipation. As expected from linear stability analysis, we find an exponentially growing instability in regimes otherwise linearly stable in resistive MHD. However, instead of vigorous and sustained magnetorotational turbulence, we find that the MRI saturates by producing large-scale, long-lived, axisymmetric structures in the magnetic and velocity fields. We refer to these structures as zonal fields and zonal flows, respectively. Their emergence causes a steep reduction in turbulent transport by at least two orders of magnitude from extrapolations based upon resistive MHD, a result that calls into question contemporary models of layered accretion. We construct a rigorous mean-field theory to explain this new behaviour and to predict when it should occur. Implications for protoplanetary disc structure and evolution, as well as for theories of planet formation, are briefly discussed. △ Less

Submitted 25 June, 2013; originally announced June 2013.

Comments: 18 pages, 16 figures, accepted for publication in MNRAS

A thermally stable heating mechanism for the intracluster medium: turbulence, magnetic fields and plasma instabilities

Authors: M. W. Kunz, A. A. Schekochihin, S. C. Cowley, J. J. Binney, J. S. Sanders

Abstract: 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 turbulen… ▽ More 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. △ Less

Submitted 21 August, 2010; v1 submitted 13 March, 2010; originally announced March 2010.

Comments: 12 pages, 4 figures, 1 table; substantial revision with essential new content