-
Relaxation to universal non-Maxwellian equilibria in a collisionless plasma
Authors:
Robert J. Ewart,
Michael L. Nastac,
Pablo J. Bilbao,
Thales Silva,
Luís O. Silva,
Alexander A. Schekochihin
Abstract:
Generic equilibria are derived for turbulent relaxing plasmas via an entropy-maximization procedure that accounts for the short-time conservation of certain collisionless invariants. The conservation of these collisionless invariants endows the system with a partial `memory' of its prior conditions, but is imperfect on long time scales due to the development of a turbulent cascade to small scales,…
▽ More
Generic equilibria are derived for turbulent relaxing plasmas via an entropy-maximization procedure that accounts for the short-time conservation of certain collisionless invariants. The conservation of these collisionless invariants endows the system with a partial `memory' of its prior conditions, but is imperfect on long time scales due to the development of a turbulent cascade to small scales, which breaks the precise conservation of phase volume, making this memory imprecise. The equilibria are still determined by the short-time collisionless invariants, but the invariants themselves are driven to a universal form by the nature of the turbulence. This is numerically confirmed for the case of beam instabilities in one-dimensional electrostatic plasmas, where sufficiently strong turbulence appears to cause the distribution function of particle energies to develop a universal power-law tail, with exponent -2.
△ Less
Submitted 3 September, 2024;
originally announced September 2024.
-
Measurement of zero-frequency fluctuations generated by coupling between Alfvén modes in the JET tokamak
Authors:
Juan Ruiz Ruiz,
Jeronimo Garcia,
Michael Barnes,
Mykola Dreval,
Carine Giroud,
Valerian H. Hall-Chen,
Michael R. Hardman,
Jon C. Hillesheim,
Yevgen Kazakov,
Samuele Mazzi,
Felix I. Parra,
Bhavin S. Patel,
Alexander A. Schekochihin,
Ziga Stancar,
the JET Contributors,
the EUROfusion Tokamak Exploitation Team
Abstract:
We report the first experimental detection of a zero-frequency fluctuation that is pumped by an Alfvén mode in a magnetically confined plasma. Core-localized bidirectional Alfvén modes of frequency inside the toroidicity-induced gap (and its harmonics) exhibit three-wave coupling interactions with a zero-frequency fluctuation. The observation of the zero-frequency fluctuation is consistent with th…
▽ More
We report the first experimental detection of a zero-frequency fluctuation that is pumped by an Alfvén mode in a magnetically confined plasma. Core-localized bidirectional Alfvén modes of frequency inside the toroidicity-induced gap (and its harmonics) exhibit three-wave coupling interactions with a zero-frequency fluctuation. The observation of the zero-frequency fluctuation is consistent with theoretical and numerical predictions of zonal modes pumped by Alfvén modes, and is correlated with an increase in the deep core ion temperature, temperature gradient, and confinement factor $H_{89,P}$. Despite the energetic particle transport induced by the Alfvén eigenmodes, the generation of a zero-frequency fluctuation that can suppress the turbulence leads to an overall improvement of confinement.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
Collisional whistler instability and electron temperature staircase in inhomogeneous plasma
Authors:
N. A. Lopez,
A. F. A. Bott,
A. A. Schekochihin
Abstract:
High-beta magnetized plasmas often exhibit anomalously structured temperature profiles, as seen from galaxy cluster observations and recent experiments. It is well known that when such plasmas are collisionless, temperature gradients along the magnetic field can excite whistler waves that efficiently scatter electrons to limit their heat transport. Only recently has it been shown that parallel tem…
▽ More
High-beta magnetized plasmas often exhibit anomalously structured temperature profiles, as seen from galaxy cluster observations and recent experiments. It is well known that when such plasmas are collisionless, temperature gradients along the magnetic field can excite whistler waves that efficiently scatter electrons to limit their heat transport. Only recently has it been shown that parallel temperature gradients can excite whistler waves also in collisional plasmas. Here we develop a Wigner--Moyal theory for the collisional whistler instability starting from Braginskii-like fluid equations in a slab geometry. This formalism is necessary because, for a large region in parameter space, the fastest-growing whistler waves have wavelengths comparable to the background temperature gradients. We find additional damping terms in the expression for the instability growth rate involving inhomogeneous Nernst advection and resistivity. They (i) enable whistler waves to re-arrange the electron temperature profile via growth, propagation, and subsequent dissipation, and (ii) allow non-constant temperature profiles to exist stably. For high-beta plasmas, the marginally stable solutions take the form of a temperature staircase along the magnetic field lines. The electron heat flux can also be suppressed by the Ettingshausen effect when the whistler intensity profile is sufficiently peaked and oriented opposite the background temperature gradient. This mechanism allows cold fronts without magnetic draping, might reduce parallel heat losses in inertial fusion experiments, and generally demonstrates that whistler waves can regulate transport even in the collisional limit.
△ Less
Submitted 16 June, 2024;
originally announced June 2024.
-
Suppression of temperature-gradient-driven turbulence by sheared flows in fusion plasmas
Authors:
P. G. Ivanov,
T. Adkins,
D. Kennedy,
M. Giacomin,
M. Barnes,
A. A. Schekochihin
Abstract:
Starting from the assumption that saturation of plasma turbulence driven by temperature-gradient instabilities in fusion plasmas is achieved by a local energy cascade between a long-wavelength outer scale, where energy is injected into the fluctuations, and a small-wavelength dissipation scale, where fluctuation energy is thermalized by particle collisions, we formulate a detailed phenomenological…
▽ More
Starting from the assumption that saturation of plasma turbulence driven by temperature-gradient instabilities in fusion plasmas is achieved by a local energy cascade between a long-wavelength outer scale, where energy is injected into the fluctuations, and a small-wavelength dissipation scale, where fluctuation energy is thermalized by particle collisions, we formulate a detailed phenomenological theory for the influence of perpendicular flow shear on magnetized-plasma turbulence. Our theory introduces two distinct regimes, called the weak-shear and strong-shear regimes, each with its own set of scaling laws for the scale and amplitude of the fluctuations and for the level of turbulent heat transport. We discover that the ratio of the typical radial and poloidal wavenumbers of the fluctuations (i.e., their aspect ratio) at the outer scale plays a central role in determining the dependence of the turbulent transport on the imposed flow shear. Our theoretical predictions are found to be in excellent agreement with numerical simulations of two paradigmatic models of fusion-relevant plasma turbulence: (i) an electrostatic fluid model of slab electron-scale turbulence, and (ii) Cyclone-base-case gyrokinetic ion-scale turbulence. Additionally, our theory envisions a potential mechanism for the suppression of electron-scale turbulence by perpendicular ion-scale flows based on the role of the aforementioned aspect ratio of the electron-scale fluctuations.
△ Less
Submitted 1 May, 2024;
originally announced May 2024.
-
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.
-
Laboratory realization of relativistic pair-plasma beams
Authors:
C. D. Arrowsmith,
P. Simon,
P. Bilbao,
A. F. A. Bott,
S. Burger,
H. Chen,
F. D. Cruz,
T. Davenne,
I. Efthymiopoulos,
D. H. Froula,
A. M. Goillot,
J. T. Gudmundsson,
D. Haberberger,
J. Halliday,
T. Hodge,
B. T. Huffman,
S. Iaquinta,
F. Miniati,
B. Reville,
S. Sarkar,
A. A. Schekochihin,
L. O. Silva,
R. Simpson,
V. Stergiou,
R. M. G. M. Trines
, et al. (4 additional authors not shown)
Abstract:
Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black holes and neutron star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with such pair plasmas. Their behaviour is quite different from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components and their role in the dy…
▽ More
Relativistic electron-positron plasmas are ubiquitous in extreme astrophysical environments such as black holes and neutron star magnetospheres, where accretion-powered jets and pulsar winds are expected to be enriched with such pair plasmas. Their behaviour is quite different from typical electron-ion plasmas due to the matter-antimatter symmetry of the charged components and their role in the dynamics of such compact objects is believed to be fundamental. So far, our experimental inability to produce large yields of positrons in quasi-neutral beams has restricted the understanding of electron-positron pair plasmas to simple numerical and analytical studies which are rather limited. We present first experimental results confirming the generation of high-density, quasi-neutral, relativistic electron-positron pair beams using the 440 GeV/c beam at CERN's Super Proton Synchrotron (SPS) accelerator. The produced pair beams have a volume that fills multiple Debye spheres and are thus able to sustain collective plasma oscillations. Our work opens up the possibility of directly probing the microphysics of pair plasmas beyond quasi-linear evolution into regimes that are challenging to simulate or measure via astronomical observations.
△ Less
Submitted 8 December, 2023;
originally announced December 2023.
-
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.
-
Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence
Authors:
Michael L. Nastac,
Robert J. Ewart,
Wrick Sengupta,
Alexander A. Schekochihin,
Michael Barnes,
William D. Dorland
Abstract:
We consider a nearly collisionless plasma consisting of a species of `test particles' in 1D-1V, stirred by an externally imposed stochastic electric field. The mean effect on the particle distribution function is stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant…
▽ More
We consider a nearly collisionless plasma consisting of a species of `test particles' in 1D-1V, stirred by an externally imposed stochastic electric field. The mean effect on the particle distribution function is stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant $C_2 \propto \iint dx dv \, \langle f^2 \rangle$. We find that $C_2$ is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of $C_2$ in Fourier space, with $k$ and $s$ denoting spatial and velocity wavenumbers, respectively. Whereas even the linear phase mixing alone would lead to a constant flux of $C_2$ to high $s$ (towards the collisional dissipation range) at every $k$, the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of $C_2$ is to both high $k$ and high $s$ simultaneously. Integrating over velocity (spatial) wavenumbers, the $k$-space ($s$-space) flux of $C_2$ is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the $(k,s)$ plane, with power-law asymptotics at large $k$ and $s$. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one.
△ Less
Submitted 14 June, 2024; v1 submitted 27 October, 2023;
originally announced October 2023.
-
Kinetic stability of Chapman-Enskog plasmas
Authors:
Archie F. A. Bott,
Steven C. Cowley,
Alexander A. Schekochihin
Abstract:
In this paper, we investigate the kinetic stability of classical, collisional plasma - that is, plasma in which the mean-free-path $λ$ of constituent particles is short compared to the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared to the evolution time. Fluid equations are typically used to describe such plasmas, si…
▽ More
In this paper, we investigate the kinetic stability of classical, collisional plasma - that is, plasma in which the mean-free-path $λ$ of constituent particles is short compared to the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared to the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman-Enskog (CE) expansion in $λ/L \ll 1$, and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $β$, the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in magnetised plasma, we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the 'whisper instability') whose growth rate in some parameter regimes is large compared to other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $λ$, the electron inertial scale $d_e$ and $β$. This work is of general consequence in emphasising the fact that high-$β$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $β> L/λ$. In this situation, the determination of transport coefficients via the standard CE approach is not valid.
△ Less
Submitted 26 October, 2023;
originally announced October 2023.
-
Intermittency of density fluctuations and zonal-flow generation in MAST edge plasmas
Authors:
A. Sladkomedova,
I. Cziegler,
A. R. Field,
A. A. Schekochihin,
D. Dunai,
P. G. Ivanov,
the MAST-U Team,
the EUROfusion MST1 Team
Abstract:
The properties of the edge ion-scale turbulence are studied using the beam emission spectroscopy (BES) diagnostic on MAST. Evidence of the formation of large-scale high-amplitude coherent structures, filamentary density blobs and holes, 2$-$4 cm inside the plasma separatrix is presented. Measurements of radial velocity and skewness of the density fluctuations indicate that density holes propagate…
▽ More
The properties of the edge ion-scale turbulence are studied using the beam emission spectroscopy (BES) diagnostic on MAST. Evidence of the formation of large-scale high-amplitude coherent structures, filamentary density blobs and holes, 2$-$4 cm inside the plasma separatrix is presented. Measurements of radial velocity and skewness of the density fluctuations indicate that density holes propagate radially inwards, with the skewness profile peaking at 7$-$10 cm inside the separatrix. Poloidal velocities of the density fluctuations measured using cross-correlation time delay estimation (CCTDE) are found to exhibit an intermittent behaviour. Zonal-flow analysis reveals the presence of poloidally symmetric coherent oscillations $-$ low-frequency (LF) zonal flows and geodesic acoustic modes (GAM). Shearing rates of the observed zonal flows are found to be comparable to the turbulence decorrelation rate. The observed bursts in density-fluctuation power are followed by quiescent periods with a transient increase in the power of sheared flows. Three-wave interactions between broadband turbulence and a GAM are illustrated using the autobispectral technique. It is shown that the zonal flows and the density-fluctuation field are nonlinearly coupled and LF zonal flows mediate the energy transfer from high- to low-frequency density fluctuations.
△ Less
Submitted 21 December, 2023; v1 submitted 11 June, 2023;
originally announced June 2023.
-
Non-thermal particle acceleration and power-law tails via relaxation to universal Lynden-Bell equilibria
Authors:
Robert J. Ewart,
Michael L. Nastac,
Alexander A. Schekochihin
Abstract:
Collisionless and weakly collisional plasmas often exhibit non-thermal quasi-equilibria. Among these quasi-equilibria, distributions with power-law tails are ubiquitous. It is shown that the statistical-mechanical approach originally suggested by Lynden-Bell (1967) can easily recover such power-law tails. Moreover, we show that, despite the apparent diversity of Lynden-Bell equilibria, a generic f…
▽ More
Collisionless and weakly collisional plasmas often exhibit non-thermal quasi-equilibria. Among these quasi-equilibria, distributions with power-law tails are ubiquitous. It is shown that the statistical-mechanical approach originally suggested by Lynden-Bell (1967) can easily recover such power-law tails. Moreover, we show that, despite the apparent diversity of Lynden-Bell equilibria, a generic form of the equilibrium distribution at high energies is a `hard' power-law tail $\propto \varepsilon^{-2}$, where $\varepsilon$ is the particle energy. The shape of the `core' of the distribution, located at low energies, retains some dependence on the initial condition but it is the tail (or `halo') that contains most of the energy. Thus, a degree of universality exists in collisionless plasmas.
△ Less
Submitted 3 September, 2023; v1 submitted 7 April, 2023;
originally announced April 2023.
-
Scale invariance and critical balance in electrostatic drift-kinetic turbulence
Authors:
T. Adkins,
P. G. Ivanov,
A. A. Schekochihin
Abstract:
The equations of electrostatic drift kinetics are observed to possess a symmetry associated with their intrinsic scale invariance. Under the assumptions of spatial periodicity, stationarity, and locality, this symmetry implies a particular scaling of the turbulent heat flux with the system's parallel size, from which its scaling with the equilibrium temperature gradient can be deduced under some a…
▽ More
The equations of electrostatic drift kinetics are observed to possess a symmetry associated with their intrinsic scale invariance. Under the assumptions of spatial periodicity, stationarity, and locality, this symmetry implies a particular scaling of the turbulent heat flux with the system's parallel size, from which its scaling with the equilibrium temperature gradient can be deduced under some additional assumptions. This macroscopic transport prediction is then confirmed numerically for a reduced model of electron-temperature-gradient-driven turbulence in slab geometry. The system realises this scaling through a turbulent cascade from large to small perpendicular spatial scales. The route of this cascade through wavenumber space (i.e. the relationship between parallel and perpendicular scales in the inertial range) is shown to be determined by a balance between nonlinear-decorrelation and parallel-dissipation timescales. This type of `critically balanced' cascade, which maintains a constant energy flux despite the presence of parallel dissipation throughout the inertial range (as well as order-unity dissipative losses at the outer scale) is expected to be a generic feature of plasma turbulence. The outer scale of the turbulence, on which the turbulent heat flux depends, is determined by the breaking of drift-kinetic scale invariance due to the existence of large-scale parallel inhomogeneity (the parallel system size).
△ Less
Submitted 21 August, 2023; v1 submitted 26 March, 2023;
originally announced March 2023.
-
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.
-
The gyrokinetic dispersion relation of microtearing modes in collisionless toroidal plasmas
Authors:
Benjamin D. G. Chandran,
Alexander A. Schekochihin
Abstract:
We solve the linearized gyrokinetic equation, quasineutrality condition, and Ampere's law to obtain the dispersion relation of microtearing modes (MTMs) in collisionless low-beta toroidal plasmas. Consistent with past studies, we find that MTMs are driven unstable by the electron temperature gradient and that this instability drive is mediated by magnetic drifts. The dispersion relation that we de…
▽ More
We solve the linearized gyrokinetic equation, quasineutrality condition, and Ampere's law to obtain the dispersion relation of microtearing modes (MTMs) in collisionless low-beta toroidal plasmas. Consistent with past studies, we find that MTMs are driven unstable by the electron temperature gradient and that this instability drive is mediated by magnetic drifts. The dispersion relation that we derive can be evaluated numerically very quickly and may prove useful for devising strategies to mitigate MTM instability in fusion devices.
△ Less
Submitted 26 February, 2024; v1 submitted 3 November, 2022;
originally announced November 2022.
-
Indirect Measurements of Gas Velocities in Galaxy Clusters: Effects of Ellipticity and Cluster Dynamic State
Authors:
Irina Zhuravleva,
Mandy C. Chen,
Eugene Churazov,
Alexander A. Schekochihin,
Congyao Zhang,
Daisuke Nagai
Abstract:
While awaiting direct velocity measurement of gas motions in the hot intracluster medium, we rely on indirect probes, including gas perturbations in galaxy clusters. Using a sample of $\sim 80$ clusters in different dynamic states from Omega500 cosmological simulations, we examine scaling relations between the fluctuation amplitudes of gas density, $δρ/ρ$, pressure, $δP/P$, X-ray surface brightnes…
▽ More
While awaiting direct velocity measurement of gas motions in the hot intracluster medium, we rely on indirect probes, including gas perturbations in galaxy clusters. Using a sample of $\sim 80$ clusters in different dynamic states from Omega500 cosmological simulations, we examine scaling relations between the fluctuation amplitudes of gas density, $δρ/ρ$, pressure, $δP/P$, X-ray surface brightness, Sunyaev-Zeldovich (SZ) y-parameter, and the characteristic Mach number of gas motions, $M_{\rm 1d}$. In relaxed clusters, accounting for halo ellipticities reduces $δρ/ρ$ or $δP/P$ by a factor of up to 2 within $r_{500c}$. We confirm a strong linear correlation between $δρ/ρ$ (or $δP/P$) and $M_{\rm 1d}$ in relaxed clusters, with the proportionality coefficient $η\approx 1$. For unrelaxed clusters, the correlation is less strong and has a larger $η\approx 1.3\pm 0.5$ ($1.5\pm0.5$) for $δρ/ρ$ ($δP/P$). Examination of the power-law scaling of $M_{\rm 1d}$ with $δρ/ρ$ shows that it is almost linear for relaxed clusters, while for the unrelaxed ones, it is closer to $δρ/ρ\propto M_{\rm 1d}^2$, supporting an increasing role of non-linear terms and compressive modes. In agreement with previous studies, we observe a strong correlation of $M_{\rm 1d}$ with radius. Correcting for these correlations leaves a residual scatter in $M_{\rm 1d}$ of $\sim 4 (7)$ per cent for relaxed (perturbed) clusters. Hydrostatic mass bias correlates with $M_{\rm 1d}$ as strongly as with $δρ/ρ$ in relaxed clusters. The residual scatters after correcting for derived trends is $\sim 6-7$ per cent. These predictions can be verified with existing X-ray and SZ observations of galaxy clusters combined with forthcoming velocity measurements with X-ray microcalorimeters.
△ Less
Submitted 20 October, 2022;
originally announced October 2022.
-
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.
-
Bubble-driven Gas Uplift in Galaxy Clusters and its Velocity Features
Authors:
Congyao Zhang,
Irina Zhuravleva,
Marie-Lou Gendron-Marsolais,
Eugene Churazov,
Alexander A. Schekochihin,
William R. Forman
Abstract:
Buoyant bubbles of relativistic plasma are essential for active galactic nucleus feedback in galaxy clusters, stirring and heating the intracluster medium (ICM). Observations suggest that these rising bubbles maintain their integrity and sharp edges much longer than predicted by hydrodynamic simulations. In this study, we assume that bubbles can be modeled as rigid bodies and demonstrate that inta…
▽ More
Buoyant bubbles of relativistic plasma are essential for active galactic nucleus feedback in galaxy clusters, stirring and heating the intracluster medium (ICM). Observations suggest that these rising bubbles maintain their integrity and sharp edges much longer than predicted by hydrodynamic simulations. In this study, we assume that bubbles can be modeled as rigid bodies and demonstrate that intact bubbles and their long-term interactions with the ambient ICM play an important role in shaping gas kinematics, forming thin gaseous structures (e.g., H$α$ filaments), and generating internal waves in cluster cores. We find that well-developed eddies are formed in the wake of a buoyantly rising bubble, and it is these eddies, rather than the Darwin drift, that are responsible for most of the gas mass uplift. The eddies gradually elongate along the bubble's direction of motion due to the strong density stratification of the atmosphere and eventually detach from the bubble, quickly evolving into a high-speed jet-like stream propagating towards the cluster center. This picture naturally explains the presence of long straight and horseshoe-shaped H$α$ filaments in the Perseus cluster, inward and outward motions of the gas, and the X-ray-weighted gas velocity distributions near the northwestern bubble observed by Hitomi. Our model reproduces the observed H$α$ velocity structure function of filaments, providing a simple interpretation for its steep scaling and normalization: laminar gas flows and large eddies within filaments driven by the intact bubbles, rather than spatially homogeneous small-scale turbulence, are sufficient to produce a structure function consistent with observations.
△ Less
Submitted 27 September, 2022; v1 submitted 8 March, 2022;
originally announced March 2022.
-
Cosmic-void observations reconciled with primordial magnetogenesis
Authors:
David N. Hosking,
Alexander A. Schekochihin
Abstract:
It has been suggested that the weak magnetic field hosted by the intergalactic medium in cosmic voids could be a relic from the early Universe. However, accepted models of turbulent magnetohydrodynamic decay predict that the present-day strength of fields originally generated at the electroweak phase transition (EWPT) without parity violation would be too low to explain the observed scattering of…
▽ More
It has been suggested that the weak magnetic field hosted by the intergalactic medium in cosmic voids could be a relic from the early Universe. However, accepted models of turbulent magnetohydrodynamic decay predict that the present-day strength of fields originally generated at the electroweak phase transition (EWPT) without parity violation would be too low to explain the observed scattering of $γ$-rays from TeV blazars. Here, we propose that the decay is mediated by magnetic reconnection and conserves the mean square fluctuation level of magnetic helicity. We find that the relic fields would be stronger by several orders of magnitude under this theory than was indicated by previous treatments, which restores the consistency of the EWPT-relic hypothesis with the observational constraints. Moreover, efficient EWPT magnetogenesis would produce relics at the strength required to resolve the Hubble tension via magnetic effects at recombination and seed galaxy-cluster fields close to their present-day strength.
△ Less
Submitted 12 January, 2023; v1 submitted 7 March, 2022;
originally announced March 2022.
-
Dimits transition in three-dimensional ion-temperature-gradient turbulence
Authors:
Plamen G. Ivanov,
Alexander A. Schekochihin,
William Dorland
Abstract:
We extend our previous work on the 2D Dimits transition in ion-scale turbulence (Ivanov et al. 2020) to include variations along the magnetic field. We consider a three-field fluid model for the perturbations of electrostatic potential, ion temperature, and ion parallel flow in a constant-magnetic-curvature geometry without magnetic shear. It is derived in the cold-ion, long-wavelength asymptotic…
▽ More
We extend our previous work on the 2D Dimits transition in ion-scale turbulence (Ivanov et al. 2020) to include variations along the magnetic field. We consider a three-field fluid model for the perturbations of electrostatic potential, ion temperature, and ion parallel flow in a constant-magnetic-curvature geometry without magnetic shear. It is derived in the cold-ion, long-wavelength asymptotic limit of the gyrokinetic theory. Just as in the 2D model, a low-transport (Dimits) regime exists and is found to be dominated by a quasi-static staircase-like arrangement of strong zonal flows and zonal temperature. This zonal staircase is formed and maintained by a negative turbulent viscosity for the zonal flows. Unlike the 2D model, the 3D one does not suffer from an unphysical blow up beyond the Dimits threshold where the staircase becomes nonlinearly unstable. Instead, a well-defined finite-amplitude saturated state is established. This qualitative difference between 2D and 3D is due to the appearance of small-scale `parasitic' modes that exist only if we allow perturbations to vary along the magnetic field lines. These modes extract energy from the large-scale perturbations and provide an effective enhancement of large-scale thermal diffusion, thus aiding the energy transfer from large injection scales to small dissipative ones. We show that in our model, the parasitic modes always favour a zonal-flow-dominated state. In fact, a Dimits state with a zonal staircase is achieved regardless of the strength of the linear drive provided the system is sufficiently extended along the magnetic field and sufficient parallel resolution is provided.
△ Less
Submitted 16 July, 2022; v1 submitted 5 March, 2022;
originally announced March 2022.
-
Three-Dimensional Inhomogeneity of Electron-Temperature-Gradient Turbulence in the Edge of Tokamak Plasmas
Authors:
J. F. Parisi,
F. I. Parra,
C. M. Roach,
M. R. Hardman,
A. A. Schekochihin,
I. G. Abel,
N. Aiba,
J. Ball,
M. Barnes,
B. Chapman-Oplopoiou,
D. Dickinson,
W. Dorland,
C. Giroud,
D. R. Hatch,
J. C. Hillesheim,
J. Ruiz Ruiz,
S. Saarelma,
D. St-Onge
Abstract:
Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a…
▽ More
Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport -- in our case by $\sim$40 \% -- through multiscale interactions.
△ Less
Submitted 2 July, 2022; v1 submitted 1 March, 2022;
originally announced March 2022.
-
Emergence of long-range correlations and thermal spectra in forced turbulence
Authors:
David N. Hosking,
Alexander A. Schekochihin
Abstract:
Recent numerical studies have shown that forced, statistically isotropic turbulence develops a `thermal equilibrium' spectrum, $\mathcal{E}(k) \propto k^2$, at large scales. This behaviour presents a puzzle, as it appears to imply the growth of a non-zero Saffman integral, which would require the longitudinal velocity correlation function, $χ(r)$, to satisfy $χ(r\to \infty)\propto r^{-3}$. As is w…
▽ More
Recent numerical studies have shown that forced, statistically isotropic turbulence develops a `thermal equilibrium' spectrum, $\mathcal{E}(k) \propto k^2$, at large scales. This behaviour presents a puzzle, as it appears to imply the growth of a non-zero Saffman integral, which would require the longitudinal velocity correlation function, $χ(r)$, to satisfy $χ(r\to \infty)\propto r^{-3}$. As is well known, the Saffman integral is an invariant of decaying turbulence, precisely because non-local interactions (i.e., interactions via exchange of pressure waves) are too weak to generate such correlations. Subject to certain restrictions on the nature of the forcing, we argue that the same should be true for forced turbulence. We show that long-range correlations and a $k^2$ spectrum arise as a result of the turbulent diffusion of linear momentum, and extend only up to a maximum scale that grows slowly with time. This picture has a number of interesting consequences. First, if the forcing generates eddies with significant linear momentum (as in so-called Saffman turbulence), a thermal spectrum is not reached - instead, a shallower spectrum develops. Secondly, the energy of turbulence that is forced for a while and then allowed to decay obeys Saffman's decay laws for a period that is much longer than the duration of the forcing stage.
△ Less
Submitted 1 August, 2023; v1 submitted 1 February, 2022;
originally announced February 2022.
-
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.
-
Electromagnetic instabilities and plasma turbulence driven by electron-temperature gradient
Authors:
T. Adkins,
A. A. Schekochihin,
P. G. Ivanov,
C. M. Roach
Abstract:
Electromagnetic (EM) instabilities and turbulence driven by the electron-temperature gradient are considered in a local slab model of a tokamak-like plasma. The model describes perturbations at scales both larger and smaller than the flux-freezing scale $d_e$, and so captures both electrostatic and EM regimes of turbulence. The well-known electrostatic instabilities -- slab and curvature-mediated…
▽ More
Electromagnetic (EM) instabilities and turbulence driven by the electron-temperature gradient are considered in a local slab model of a tokamak-like plasma. The model describes perturbations at scales both larger and smaller than the flux-freezing scale $d_e$, and so captures both electrostatic and EM regimes of turbulence. The well-known electrostatic instabilities -- slab and curvature-mediated ETG -- are recovered, and a new instability is found in the EM regime, called the Thermo-Alfvénic instability (TAI). It exists in both a slab version (sTAI, destabilising kinetic Alfvén waves) and a curvature-mediated version (cTAI), which is a cousin of the (electron-scale) kinetic ballooning mode (KBM). The cTAI turns out to be dominant at the largest scales covered by the model (greater than $d_e$ but smaller than $ρ_i$), its physical mechanism hinging on the fast equalisation of the total temperature along perturbed magnetic field lines (in contrast to KBM, which is pressure balanced). A turbulent cascade theory is then constructed, with two energy-injection scales: $d_e$, where the drivers are slab ETG and sTAI, and a larger (parallel system size dependent) scale, where the driver is cTAI. The latter dominates the turbulent transport if the temperature gradient is greater than a certain critical value, which scales inversely with the electron beta. The resulting heat flux scales more steeply with the temperature gradient than that due to electrostatic ETG turbulence, giving rise to stiffer transport. This can be viewed as a physical argument in favour of near-marginal steady-state in electron-transport-controlled plasmas (e.g., the pedestal). While the model is simplistic, the new physics that is revealed by it should be of interest to those attempting to model the effect of EM turbulence in tokamak-relevant configurations with high beta and large electron temperature gradients.
△ Less
Submitted 1 July, 2022; v1 submitted 14 January, 2022;
originally announced January 2022.
-
Collisionless relaxation of a Lynden-Bell plasma
Authors:
Robert J. Ewart,
Andrew Brown,
Toby Adkins,
Alexander A. Schekochihin
Abstract:
Plasmas whose Coulomb-collision rates are very small may relax on shorter time scales to non-Maxwellian quasi-equilibria, which, nevertheless, have a universal form, with dependence on initial conditions retained only via an infinite set of Casimir invariants enforcing phase-volume conservation. These are distributions derived by Lynden-Bell (1967) via a statistical-mechanical entropy-maximisation…
▽ More
Plasmas whose Coulomb-collision rates are very small may relax on shorter time scales to non-Maxwellian quasi-equilibria, which, nevertheless, have a universal form, with dependence on initial conditions retained only via an infinite set of Casimir invariants enforcing phase-volume conservation. These are distributions derived by Lynden-Bell (1967) via a statistical-mechanical entropy-maximisation procedure, assuming perfect mixing of phase-space elements. To show that these equilibria are reached dynamically, one must derive an effective `collisionless collision integral' for which they are fixed points -- unique and inevitable provided the integral has an appropriate H-theorem. We describe how such collision integrals are derived and what assumptions are required for them to have a closed form, how to prove the H-theorems for them, and why, for a system carrying sufficiently large electric-fluctuation energy, collisionless relaxation should be fast. It is suggested that collisionless dynamics may favour maximising entropy locally in phase space before converging to global maximum-entropy states. Relaxation due to interspecies interaction is examined, leading, inter alia, to spontaneous transient generation of electron currents. The formalism also allows efficient recovery of `true' collision integrals for both classical and quantum plasmas.
△ Less
Submitted 8 July, 2022; v1 submitted 10 January, 2022;
originally announced January 2022.
-
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.
-
Energy partition between Alfvénic and compressive fluctuations in magnetorotational turbulence with near-azimuthal mean magnetic field
Authors:
Y. Kawazura,
A. A. Schekochihin,
M. Barnes,
W. Dorland,
S. A. Balbus
Abstract:
The theory of magnetohydrodynamic (MHD) turbulence predicts that Alfvénic and slow-mode-like compressive fluctuations are energetically decoupled at small scales in the inertial range. The partition of energy between these fluctuations determines the nature of dissipation, which, in many astrophysical systems, happens on scales where plasma is collisionless. However, when the magnetorotational ins…
▽ More
The theory of magnetohydrodynamic (MHD) turbulence predicts that Alfvénic and slow-mode-like compressive fluctuations are energetically decoupled at small scales in the inertial range. The partition of energy between these fluctuations determines the nature of dissipation, which, in many astrophysical systems, happens on scales where plasma is collisionless. However, when the magnetorotational instability (MRI) drives the turbulence, it is difficult to resolve numerically the scale at which both types of fluctuations start to be decoupled because the MRI energy injection occurs in a broad range of wavenumbers, and both types of fluctuations are usually expected to be coupled even at relatively small scales. In this study, we focus on collisional MRI turbulence threaded by a near-azimuthal mean magnetic field, which is naturally produced by the differential rotation of a disc. We show that, in such a case, the decoupling scales are reachable using a reduced MHD model that includes differential-rotation effects. In our reduced MHD model, the Alfvénic and compressive fluctuations are coupled only through the linear terms that are proportional to the angular velocity of the accretion disc. We numerically solve for the turbulence in this model and show that the Alfvénic and compressive fluctuations are decoupled at the small scales of our simulations as the nonlinear energy transfer dominates the linear coupling below the MRI-injection scale. In the decoupling scales, the energy flux of compressive fluctuations contained in the small scales is almost double that of Alfvénic fluctuations. Finally, we discuss the application of this result to prescriptions of ion-to-electron heating ratio in hot accretion flows.
△ Less
Submitted 6 May, 2022; v1 submitted 24 October, 2021;
originally announced October 2021.
-
Bistable turbulence in strongly magnetised plasmas with a sheared mean flow
Authors:
Nicolas Christen,
Michael Barnes,
Michael R. Hardman,
Alexander A. Schekochihin
Abstract:
The prevailing paradigm for plasma turbulence associates a unique stationary state to given equilibrium parameters. We report the discovery of bistable turbulence in a strongly magnetised plasma. Two distinct states, obtained with identical equilibrium parameters in first-principle gyrokinetic simulations, have turbulent fluxes of particles, momentum and energy that differ by an order of magnitude…
▽ More
The prevailing paradigm for plasma turbulence associates a unique stationary state to given equilibrium parameters. We report the discovery of bistable turbulence in a strongly magnetised plasma. Two distinct states, obtained with identical equilibrium parameters in first-principle gyrokinetic simulations, have turbulent fluxes of particles, momentum and energy that differ by an order of magnitude, with the low-transport state agreeing with experimental observations. Occurrences of the two states are regulated by the competition between an externally imposed mean flow shear and "zonal" flows generated by the plasma. With small turbulent amplitudes, zonal flows have little impact, and the mean shear causes turbulence to saturate in a low-transport state. With larger amplitudes, the zonal shear can (partially) oppose the effect of the mean shear, allowing the system to sustain a high-transport state. This poses a new challenge for research that has so far assumed a uniquely defined turbulent state.
△ Less
Submitted 8 July, 2022; v1 submitted 16 September, 2021;
originally announced September 2021.
-
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.
-
Strong suppression of heat conduction in a laboratory replica of galaxy-cluster turbulent plasmas
Authors:
J. Meinecke,
P. Tzeferacos,
J. S. Ross,
A. F. A. Bott,
S. Feister,
H. -S. Park,
A. R. Bell,
R. Blandford,
R. L. Berger,
R. Bingham,
A. Casner,
L. E. Chen,
J. Foster,
D. H. Froula,
C. Goyon,
D. Kalantar,
M. Koenig,
B. Lahmann,
C. -K. Li,
Y. Lu,
C. A. J. Palmer,
R. Petrasso,
H. Poole,
B. Remington,
B. Reville
, et al. (10 additional authors not shown)
Abstract:
Galaxy clusters are filled with hot, diffuse X-ray emitting plasma, with a stochastically tangled magnetic field whose energy is close to equipartition with the energy of the turbulent motions \cite{zweibel1997, Vacca}. In the cluster cores, the temperatures remain anomalously high compared to what might be expected considering that the radiative cooling time is short relative to the Hubble time \…
▽ More
Galaxy clusters are filled with hot, diffuse X-ray emitting plasma, with a stochastically tangled magnetic field whose energy is close to equipartition with the energy of the turbulent motions \cite{zweibel1997, Vacca}. In the cluster cores, the temperatures remain anomalously high compared to what might be expected considering that the radiative cooling time is short relative to the Hubble time \cite{cowie1977,fabian1994}. While feedback from the central active galactic nuclei (AGN) \cite{fabian2012,birzan2012,churazov2000} is believed to provide most of the heating, there has been a long debate as to whether conduction of heat from the bulk to the core can help the core to reach the observed temperatures \cite{narayan2001,ruszkowski2002,kunz2011}, given the presence of tangled magnetic fields. Interestingly, evidence of very sharp temperature gradients in structures like cold fronts implies a high degree of suppression of thermal conduction \cite{markevitch2007}. To address the problem of thermal conduction in a magnetized and turbulent plasma, we have created a replica of such a system in a laser laboratory experiment. Our data show a reduction of local heat transport by two orders of magnitude or more, leading to strong temperature variations on small spatial scales, as is seen in cluster plasmas \cite{markevitch2003}.
△ Less
Submitted 18 May, 2021;
originally announced May 2021.
-
Reconnection-controlled decay of magnetohydrodynamic turbulence and the role of invariants
Authors:
David N. Hosking,
Alexander A. Schekochihin
Abstract:
We present a new theoretical picture of magnetically dominated, decaying turbulence in the absence of a mean magnetic field. We demonstrate that such turbulence is governed by the reconnection of magnetic structures, and not by ideal dynamics, as has previously been assumed. We obtain predictions for the magnetic-energy decay laws by proposing that turbulence decays on reconnection timescales, whi…
▽ More
We present a new theoretical picture of magnetically dominated, decaying turbulence in the absence of a mean magnetic field. We demonstrate that such turbulence is governed by the reconnection of magnetic structures, and not by ideal dynamics, as has previously been assumed. We obtain predictions for the magnetic-energy decay laws by proposing that turbulence decays on reconnection timescales, while respecting the conservation of certain integral invariants representing topological constraints satisfied by the reconnecting magnetic field. As is well known, the magnetic helicity is such an invariant for initially helical field configurations, but does not constrain non-helical decay, where the volume-averaged magnetic-helicity density vanishes. For such a decay, we propose a new integral invariant, analogous to the Loitsyansky and Saffman invariants of hydrodynamic turbulence, that expresses the conservation of the random (scaling as $\mathrm{volume}^{1/2}$) magnetic helicity contained in any sufficiently large volume. Our treatment leads to novel predictions for the magnetic-energy decay laws: in particular, while we expect the canonical $t^{-2/3}$ power law for helical turbulence when reconnection is fast (i.e., plasmoid-dominated or stochastic), we find a shallower $t^{-4/7}$ decay in the slow `Sweet-Parker' reconnection regime, in better agreement with existing numerical simulations. For non-helical fields, for which there currently exists no definitive theory, we predict power laws of $t^{-10/9}$ and $t^{-20/17}$ in the fast- and slow-reconnection regimes, respectively. We formulate a general principle of decay of turbulent systems subject to conservation of Saffman-like invariants, and propose how it may be applied to MHD turbulence with a strong mean magnetic field and to isotropic MHD turbulence with initial equipartition between the magnetic and kinetic energies.
△ Less
Submitted 2 September, 2021; v1 submitted 2 December, 2020;
originally announced December 2020.
-
MHD Turbulence: A Biased Review
Authors:
Alexander A. Schekochihin
Abstract:
This review puts the developments of the last few years in the context of the canonical time line (Kolmogorov to Iroshnikov-Kraichnan to Goldreich-Sridhar to Boldyrev). It is argued that Beresnyak's objection that Boldyrev's alignment theory violates the RMHD rescaling symmetry can be reconciled with alignment if the latter is understood as an intermittency effect. Boldyrev's scalings, recovered i…
▽ More
This review puts the developments of the last few years in the context of the canonical time line (Kolmogorov to Iroshnikov-Kraichnan to Goldreich-Sridhar to Boldyrev). It is argued that Beresnyak's objection that Boldyrev's alignment theory violates the RMHD rescaling symmetry can be reconciled with alignment if the latter is understood as an intermittency effect. Boldyrev's scalings, recovered in this interpretation, are thus an example of a physical theory of intermittency in a turbulent system. Emergence of aligned structures brings in reconnection physics, so the theory of MHD turbulence intertwines with the physics of tearing and current-sheet disruption. Recent work on this by Loureiro, Mallet et al. is reviewed and it is argued that we finally have a reasonably complete picture of MHD cascade all the way to the dissipation scale. This picture appears to reconcile Beresnyak's Kolmogorov scaling of the dissipation cutoff with Boldyrev's aligned cascade. These ideas also enable some progress in understanding saturated MHD dynamo, argued to be controlled by reconnection and to contain, at small scales, a tearing-mediated cascade similar to its strong-mean-field counterpart. On the margins of this core narrative, standard weak-MHD-turbulence theory is argued to require adjustment - and a scheme for it is proposed - to take account of the part that a spontaneously emergent 2D condensate plays in mediating the Alfven-wave cascade. This completes the picture of the MHD cascade at large scales. A number of outstanding issues are surveyed, concerning imbalanced MHD turbulence (for which a new theory is proposed), residual energy, subviscous and decaying regimes of MHD turbulence (where reconnection again features prominently). Finally, it is argued that the natural direction of research is now away from MHD and into kinetic territory.
△ Less
Submitted 19 July, 2022; v1 submitted 1 October, 2020;
originally announced October 2020.
-
Violation of the zeroth law of turbulence in space plasmas
Authors:
Romain Meyrand,
Jonathan Squire,
Alexander A. Schekochihin,
William Dorland
Abstract:
The zeroth law of turbulence states that, for fixed energy input into large-scale motions, the statistical steady state of a turbulent system is independent of microphysical dissipation properties. The behavior, which is fundamental to nearly all fluid-like systems from industrial processes to galaxies, occurs because nonlinear processes generate smaller and smaller scales in the flow, until the d…
▽ More
The zeroth law of turbulence states that, for fixed energy input into large-scale motions, the statistical steady state of a turbulent system is independent of microphysical dissipation properties. The behavior, which is fundamental to nearly all fluid-like systems from industrial processes to galaxies, occurs because nonlinear processes generate smaller and smaller scales in the flow, until the dissipation -- no matter how small -- can thermalize the energy input. Using direct numerical simulations and theoretical arguments, we show that in strongly magnetized plasma turbulence such as that recently observed by the Parker Solar Probe (PSP) spacecraft, the zeroth law is routinely violated. Namely, when such turbulence is "imbalanced" -- when the large-scale energy input is dominated by Alfvén waves propagating in one direction (the most common situation in space plasmas) -- nonlinear conservation laws imply the existence of a "barrier" at scales near the ion gyroradius. This causes energy to build up over time at large scales. The resulting magnetic-energy spectra bear a strong similarity to those observed in situ, exhibiting a sharp, steep kinetic transition range above and around the ion-Larmor scale, with flattening at yet smaller scales, thus resolving the decade-long puzzle of the position and variability of ion-kinetic spectral breaks in plasma turbulence. The "barrier" effect also suggests that how a plasma is forced at large scales (the imbalance) may have a crucial influence on thermodynamic properties such as the ion-to-electron heating ratio.
△ Less
Submitted 10 June, 2021; v1 submitted 6 September, 2020;
originally announced September 2020.
-
Inefficient magnetic-field amplification in supersonic laser-plasma turbulence
Authors:
A. F. A. Bott,
L. Chen,
G. Boutoux,
T. Caillaud,
A. Duval,
M. Koenig,
B. Khiar,
I. Lantuéjoul,
L. Le-Deroff,
B. Reville,
R. Rosch,
D. Ryu,
C. Spindloe,
B. Vauzour,
B. Villette,
A. A. Schekochihin,
D. Q. Lamb,
P. Tzeferacos,
G. Gregori,
A. Casner
Abstract:
We report a laser-plasma experiment that was carried out at the LMJ-PETAL facility and realized the first magnetized, turbulent, supersonic plasma with a large magnetic Reynolds number ($\mathrm{Rm} \approx 45$) in the laboratory. Initial seed magnetic fields were amplified, but only moderately so, and did not become dynamically significant. A notable absence of magnetic energy at scales smaller t…
▽ More
We report a laser-plasma experiment that was carried out at the LMJ-PETAL facility and realized the first magnetized, turbulent, supersonic plasma with a large magnetic Reynolds number ($\mathrm{Rm} \approx 45$) in the laboratory. Initial seed magnetic fields were amplified, but only moderately so, and did not become dynamically significant. A notable absence of magnetic energy at scales smaller than the outer scale of the turbulent cascade was also observed. Our results support the notion that moderately supersonic, low-magnetic-Prandtl-number plasma turbulence is inefficient at amplifying magnetic fields.
△ Less
Submitted 14 August, 2020;
originally announced August 2020.
-
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.
-
Elasticity of tangled magnetic fields
Authors:
D. N. Hosking,
A. A. Schekochihin,
S. A. Balbus
Abstract:
The fundamental difference between incompressible ideal magnetohydrodynamics and the dynamics of a non-conducting fluid is that magnetic fields exert a tension force that opposes their bending; magnetic fields behave like elastic strings threading the fluid. It is natural, therefore, to expect that a magnetic field tangled at small length scales should resist a large-scale shear in an elastic way,…
▽ More
The fundamental difference between incompressible ideal magnetohydrodynamics and the dynamics of a non-conducting fluid is that magnetic fields exert a tension force that opposes their bending; magnetic fields behave like elastic strings threading the fluid. It is natural, therefore, to expect that a magnetic field tangled at small length scales should resist a large-scale shear in an elastic way, much as a ball of tangled elastic strings responds elastically to an impulse. Furthermore, a tangled field should support the propagation of `magnetoelastic waves', the isotropic analogue of Alfvén waves on a straight magnetic field. Here, we study magnetoelasticity in the idealised context of an equilibrium tangled field configuration. In contrast to previous treatments, we explicitly account for intermittency of the Maxwell stress, and show that this intermittency necessarily decreases the frequency of magnetoelastic waves in a stable field configuration. We develop a mean-field formalism to describe magnetoelastic behaviour, retaining leading-order corrections due to the coupling of large- and small-scale motions, and solve the initial-value problem for viscous fluids subjected to a large-scale shear, showing that the development of small-scale motions results in anomalous viscous damping of large-scale waves. Finally, we test these analytic predictions using numerical simulations of standing waves on tangled, linear force-free magnetic-field equilibria.
△ Less
Submitted 26 June, 2020;
originally announced June 2020.
-
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.
-
Ion versus electron heating in compressively driven astrophysical gyrokinetic turbulence
Authors:
Y. Kawazura,
A. A. Schekochihin,
M. Barnes,
J. M. TenBarge,
Y. Tong,
K. G. Klein,
W. Dorland
Abstract:
The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio $Q_\rmi/Q_\rme$ as a function of the compressive-to-Alfvénic driving power ratio $P_\compr/P_\AW$, of the ratio of ion thermal pressure to…
▽ More
The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio $Q_\rmi/Q_\rme$ as a function of the compressive-to-Alfvénic driving power ratio $P_\compr/P_\AW$, of the ratio of ion thermal pressure to magnetic pressure $β_\rmi$, and of the ratio of ion-to-electron background temperatures $T_\rmi/T_\rme$. It is shown that $Q_\rmi/Q_\rme$ is an increasing function of $P_\compr/P_\AW$. When the compressive driving is sufficiently large, $Q_\rmi/Q_\rme$ approaches $\simeq P_\compr/P_\AW$. This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfvén waves is absent for both low and high $β_\rmi$, meaning that the compressive driving is directly connected to the ion entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low $β_\rmi$ and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.
△ Less
Submitted 23 October, 2020; v1 submitted 10 April, 2020;
originally announced April 2020.
-
Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence
Authors:
P. G. Ivanov,
A. A. Schekochihin,
W. Dorland,
A. R. Field,
F. I. Parra
Abstract:
The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$-pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a…
▽ More
The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$-pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal-temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations. A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear. The break up of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system -- the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol{E}\times\boldsymbol{B}$ flow. By analysing the linear ITG modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.
△ Less
Submitted 16 July, 2020; v1 submitted 8 April, 2020;
originally announced April 2020.
-
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.
-
Suppressed effective viscosity in the bulk intergalactic plasma
Authors:
I. Zhuravleva,
E. Churazov,
A. A. Schekochihin,
S. W. Allen,
A. Vikhlinin,
N. Werner
Abstract:
Transport properties, such as viscosity and thermal conduction, of the hot intergalactic plasma in clusters of galaxies, are largely unknown. While for laboratory plasmas these characteristics are derived from the gas density and temperature, such recipes can be fundamentally different for the intergalactic plasma due to a low rate of particle collisions and a weak magnetic field. In numerical sim…
▽ More
Transport properties, such as viscosity and thermal conduction, of the hot intergalactic plasma in clusters of galaxies, are largely unknown. While for laboratory plasmas these characteristics are derived from the gas density and temperature, such recipes can be fundamentally different for the intergalactic plasma due to a low rate of particle collisions and a weak magnetic field. In numerical simulations, one often cuts through these unknowns by modeling these plasmas as hydrodynamic fluids, even though local, non-hydrodynamic features observed in clusters contradict this assumption. Using deep Chandra observations of the Coma Cluster, we probe gas fluctuations in intergalactic medium down to spatial scales where the transport processes should prominently manifest themselves - at least if hydrodynamic models with pure Coulomb collision rates were indeed adequate. We find that they do not, implying that the effective isotropic viscosity is orders of magnitude smaller than naively expected. This indicates an enhanced collision rate in the plasma due to particle scattering off microfluctuations caused by plasma instabilities, or that the transport processes are anisotropic with respect to local magnetic field. For that reason, numerical models with high Reynolds number appear more consistent with observations. Our results also demonstrate that observations of turbulence in clusters are becoming a branch of astrophysics that can sharpen theoretical views on such plasmas.
△ Less
Submitted 14 June, 2019;
originally announced June 2019.
-
[Plasma 2020 Decadal] Disentangling the Spatiotemporal Structure of Turbulence Using Multi-Spacecraft Data
Authors:
J. M. TenBarge,
O. Alexandrova,
S. Boldyrev,
F. Califano,
S. S. Cerri,
C. H. K. Chen,
G. G. Howes,
T. Horbury,
P. A. Isenberg,
H. Ji,
K. G. Klein,
C. Krafft,
M. Kunz,
N. F. Loureiro,
A. Mallet,
B. A. Maruca,
W. H. Matthaeus,
R. Meyrand,
E. Quataert,
J. C. Perez,
O. W. Roberts,
F. Sahraoui,
C. S. Salem,
A. A. Schekochihin,
H. Spence
, et al. (4 additional authors not shown)
Abstract:
This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and…
▽ More
This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.
△ Less
Submitted 13 March, 2019;
originally announced March 2019.
-
[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.
-
Constraints on ion vs. electron heating by plasma turbulence at low beta
Authors:
A. A. Schekochihin,
Y. Kawazura,
M. A. Barnes
Abstract:
It is shown that in low-beta, weakly collisional plasmas, such as the solar corona, some instances of the solar wind, the aurora, inner regions of accretion discs, their coronae, and some laboratory plasmas, Alfvénic fluctuations produce no ion heating within the gyrokinetic approximation, i.e., as long as their amplitudes (at the Larmor scale) are small and their frequencies stay below the ion La…
▽ More
It is shown that in low-beta, weakly collisional plasmas, such as the solar corona, some instances of the solar wind, the aurora, inner regions of accretion discs, their coronae, and some laboratory plasmas, Alfvénic fluctuations produce no ion heating within the gyrokinetic approximation, i.e., as long as their amplitudes (at the Larmor scale) are small and their frequencies stay below the ion Larmor frequency (even as their spatial scales can be above or below the ion Larmor scale). Thus, all low-frequency ion heating in such plasmas is due to compressive fluctuations ("slow modes"). Because these fluctuations energetically decouple from the Alfvénic ones already in the inertial range, the above conclusion means that the energy partition between ions and electrons in low-beta plasmas is decided at the outer scale, where turbulence is launched, and can be determined from magnetohydrodynamic (MHD) models of the relevant astrophysical systems. Any additional ion heating must come from non-gyrokinetic mechanisms such as cyclotron heating or the stochastic heating owing to distortions of ions' Larmor orbits. An exception to these conclusions occurs in the Hall limit, i.e., when the ratio of the ion to electron temperatures is as low as the ion beta (equivalently, the electron beta is order unity). In this regime, slow modes couple to Alfvénic ones well above the Larmor scale (viz., at the ion inertial or ion sound scale), so the Alfvénic and compressive cascades join and then separate again into two cascades of fluctuations that linearly resemble kinetic Alfvén and ion cyclotron waves, with the former heating electrons and the latter ions. The two cascades are shown to decouple, scalings for them are derived, and it is argued physically that the two species will be heated by them at approximately equal rates.
△ Less
Submitted 10 April, 2019; v1 submitted 23 December, 2018;
originally announced December 2018.
-
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.
-
Transport of high-energy charged particles through spatially-intermittent turbulent magnetic fields
Authors:
L. E. Chen,
A. F. A. Bott,
P. Tzeferacos,
A. Rigby,
A. Bell,
R. Bingham,
C. Graziani,
J. Katz,
M. Koenig,
C. K. Li,
R. Petrasso,
H. -S. Park,
J. S. Ross,
D. Ryu,
T. G. White,
B. Reville,
J. Matthews,
J. Meinecke,
F. Miniati,
E. G. Zweibel,
S. Sarkar,
A. A. Schekochihin,
D. Q. Lamb,
D. H. Froula,
G. Gregori
Abstract:
Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results sho…
▽ More
Identifying the sources of the highest energy cosmic rays requires understanding how they are deflected by the stochastic, spatially intermittent intergalactic magnetic field. Here we report measurements of energetic charged-particle propagation through a laser-produced magnetized plasma with these properties. We characterize the diffusive transport of the particles experimentally. The results show that the transport is diffusive and that, for the regime of interest for the highest-energy cosmic rays, the diffusion coefficient is unaffected by the spatial intermittency of the magnetic field.
△ Less
Submitted 30 March, 2020; v1 submitted 13 August, 2018;
originally announced August 2018.
-
Fluidization of collisionless plasma turbulence
Authors:
Romain Meyrand,
Anjor Kanekar,
William Dorland,
Alexander A. Schekochihin
Abstract:
In a collisionless, magnetized plasma, particles may stream freely along magnetic-field lines, leading to phase "mixing" of their distribution function and consequently to smoothing out of any "compressive" fluctuations (of density, pressure, etc.,). This rapid mixing underlies Landau damping of these fluctuations in a quiescent plasma-one of the most fundamental physical phenomena that make plasm…
▽ More
In a collisionless, magnetized plasma, particles may stream freely along magnetic-field lines, leading to phase "mixing" of their distribution function and consequently to smoothing out of any "compressive" fluctuations (of density, pressure, etc.,). This rapid mixing underlies Landau damping of these fluctuations in a quiescent plasma-one of the most fundamental physical phenomena that make plasma different from a conventional fluid. Nevertheless, broad power-law spectra of compressive fluctuations are observed in turbulent astrophysical plasmas (most vividly, in the solar wind) under conditions conducive to strong Landau damping. Elsewhere in nature, such spectra are normally associated with fluid turbulence, where energy cannot be dissipated in the inertial scale range and is therefore cascaded from large scales to small. By direct numerical simulations and theoretical arguments, it is shown here that turbulence of compressive fluctuations in collisionless plasmas strongly resembles one in a collisional fluid and does have broad power-law spectra. This "fluidization" of collisionless plasmas occurs because phase mixing is strongly suppressed on average by "stochastic echoes", arising due to nonlinear advection of the particle distribution by turbulent motions. Besides resolving the long-standing puzzle of observed compressive fluctuations in the solar wind, our results suggest a conceptual shift for understanding kinetic plasma turbulence generally: rather than being a system where Landau damping plays the role of dissipation, a collisionless plasma is effectively dissipationless except at very small scales. The universality of "fluid" turbulence physics is thus reaffirmed even for a kinetic, collisionless system.
△ Less
Submitted 13 August, 2018;
originally announced August 2018.
-
Thermal disequilibration of ions and electrons by collisionless plasma turbulence
Authors:
Yohei Kawazura,
Michael Barnes,
Alexander A. Schekochihin
Abstract:
Does overall thermal equilibrium exist between ions and electrons in a weakly collisional, magnetised, turbulent plasma---and, if not, how is thermal energy partitioned between ions and electrons? This is a fundamental question in plasma physics, the answer to which is also crucial for predicting the properties of far-distant astronomical objects such as accretion discs around black holes. In the…
▽ More
Does overall thermal equilibrium exist between ions and electrons in a weakly collisional, magnetised, turbulent plasma---and, if not, how is thermal energy partitioned between ions and electrons? This is a fundamental question in plasma physics, the answer to which is also crucial for predicting the properties of far-distant astronomical objects such as accretion discs around black holes. In the context of discs, this question was posed nearly two decades ago and has since generated a sizeable literature. Here we provide the answer for the case in which energy is injected into the plasma via Alfvénic turbulence: collisionless turbulent heating typically acts to disequilibrate the ion and electron temperatures. Numerical simulations using a hybrid fluid-gyrokinetic model indicate that the ion-electron heating-rate ratio is an increasing function of the thermal-to-magnetic energy ratio, $β_\mathrm{i}$: it ranges from $\sim0.05$ at $β_\mathrm{i}=0.1$ to at least $30$ for $β_\mathrm{i} \gtrsim 10$. This energy partition is approximately insensitive to the ion-to-electron temperature ratio $T_\mathrm{i}/T_\mathrm{e}$. Thus, in the absence of other equilibrating mechanisms, a collisionless plasma system heated via Alfvénic turbulence will tend towards a nonequilibrium state in which one of the species is significantly hotter than the other, viz., hotter ions at high $β_\mathrm{i}$, hotter electrons at low $β_\mathrm{i}$. Spectra of electromagnetic fields and the ion distribution function in 5D phase space exhibit an interesting new magnetically dominated regime at high $β_i$ and a tendency for the ion heating to be mediated by nonlinear phase mixing ("entropy cascade") when $β_\mathrm{i}\lesssim1$ and by linear phase mixing (Landau damping) when $β_\mathrm{i}\gg1$
△ Less
Submitted 15 January, 2019; v1 submitted 19 July, 2018;
originally announced July 2018.
-
Generation of Internal Waves by Buoyant Bubbles in Galaxy Clusters and Heating of Intracluster Medium
Authors:
Congyao Zhang,
Eugene Churazov,
Alexander A. Schekochihin
Abstract:
Buoyant bubbles of relativistic plasma in cluster cores plausibly play a key role in conveying the energy from a supermassive black hole to the intracluster medium (ICM) - the process known as radio-mode AGN feedback. Energy conservation guarantees that a bubble loses most of its energy to the ICM after crossing several pressure scale heights. However, actual processes responsible for transferring…
▽ More
Buoyant bubbles of relativistic plasma in cluster cores plausibly play a key role in conveying the energy from a supermassive black hole to the intracluster medium (ICM) - the process known as radio-mode AGN feedback. Energy conservation guarantees that a bubble loses most of its energy to the ICM after crossing several pressure scale heights. However, actual processes responsible for transferring the energy to the ICM are still being debated. One attractive possibility is the excitation of internal waves, which are trapped in the cluster's core and eventually dissipate. Here we show that a sufficient condition for efficient excitation of these waves in stratified cluster atmospheres is flattening of the bubbles in the radial direction. In our numerical simulations, we model the bubbles phenomenologically as rigid bodies buoyantly rising in the stratified cluster atmosphere. We find that the terminal velocities of the flattened bubbles are small enough so that the Froude number ${\rm Fr}\lesssim 1$. The effects of stratification make the dominant contribution to the total drag force balancing the buoyancy force. In particular, clear signs of internal waves are seen in the simulations. These waves propagate horizontally and downwards from the rising bubble, spreading their energy over large volumes of the ICM. If our findings are scaled to the conditions of the Perseus cluster, the expected terminal velocity is $\sim100-200{\,\rm km\,s^{-1}}$ near the cluster cores, which is in broad agreement with direct measurements by the Hitomi satellite.
△ Less
Submitted 25 February, 2018;
originally announced February 2018.
-
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.
-
A solvable model of Vlasov-kinetic plasma turbulence in Fourier-Hermite phase space
Authors:
T. Adkins,
A. A. Schekochihin
Abstract:
A class of simple kinetic systems is considered, described by the 1D Vlasov-Landau equation with Poisson or Boltzmann electrostatic response and an energy source. Assuming a stochastic electric field, a solvable model is constructed for the phase-space turbulence of the particle distribution. The model is a kinetic analog of the Kraichnan-Batchelor model of chaotic advection. The solution of the m…
▽ More
A class of simple kinetic systems is considered, described by the 1D Vlasov-Landau equation with Poisson or Boltzmann electrostatic response and an energy source. Assuming a stochastic electric field, a solvable model is constructed for the phase-space turbulence of the particle distribution. The model is a kinetic analog of the Kraichnan-Batchelor model of chaotic advection. The solution of the model is found in Fourier-Hermite space and shows that the free-energy flux from low to high Hermite moments is suppressed, with phase mixing cancelled on average by anti-phase-mixing (stochastic plasma echo). This implies that Landau damping is an ineffective route to dissipation (i.e., to thermalisation of electric energy via velocity space). The full Fourier-Hermite spectrum is derived. Its asymptotics are $m^{-3/2}$ at low wave numbers and high Hermite moments ($m$) and $m^{-1/2}k^{-2}$ at low Hermite moments and high wave numbers ($k$). These conclusions hold at wave numbers below a certain cut off (analog of Kolmogorov scale), which increases with the amplitude of the stochastic electric field and scales as inverse square of the collision rate. The energy distribution and flows in phase space are a simple and, therefore, useful example of competition between phase mixing and nonlinear dynamics in kinetic turbulence, reminiscent of more realistic but more complicated multi-dimensional systems that have not so far been amenable to complete analytical solution.
△ Less
Submitted 7 January, 2018; v1 submitted 10 September, 2017;
originally announced September 2017.