Showing 1–50 of 115 results for author: Schekochihin, A A

Cosmic-ray confinement in radio bubbles by micromirrors

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

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

Submitted 7 April, 2024; originally announced April 2024.

Comments: 10 pages, 6 figures

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.

Comments: 14 pages, 8 figures

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

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

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

Submitted 2 November, 2023; originally announced November 2023.

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

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.

Comments: Minor edits, final version accepted to PRE

Journal ref: Phys. Rev. E 109, 065210 (2024)

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.

Comments: 182 pages total (99 main text, remaining appendices), 31 figures

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.

Comments: 33 pages, 5 figures

Pressure anisotropy and viscous heating in weakly collisional plasma turbulence

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

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

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

Comments: Accepted for publication in Journal of Plasma Physics

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.

Comments: 15 pages, 11 figures, submitted

Kinetic Turbulence in Collisionless High-Beta Plasmas

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

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

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

Comments: 34 pages, 19 figures, accepted to PRX

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.

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

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.

Comments: 21 pages, including supplementary information; 6 figures

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.

Comments: 51 pages, 7 figures, accepted by J. Fluid Mech

Tearing instability and current-sheet disruption in the turbulent dynamo

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

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

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

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

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.

Comments: 46 pages, 4 figures, submitted to JPP

Insensitivity of a turbulent laser-plasma dynamo to initial conditions

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

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

Submitted 5 January, 2022; originally announced January 2022.

Comments: 13 pages, 12 figures

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.

Comments: Accepted for publication in Journal of Plasma Physics

Journal ref: Journal of Plasma Physics 88, 905880311 (2022)

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.

Journal ref: Sci. Adv. 8, eabj6799 (2022)

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.

Comments: 30 pages, 15 figures, accepted by Phys. Rev. X

Journal ref: Phys. Rev. X 11, 041005 (2021)

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.

Comments: 200 pages, 51 figures; final version accepted by JPP, contains some substantial modifications and additions compared to v2, especially in secs 6.1, 7.1.2, 7.4, 8, 9.2, 13.4.2, B, D.5.2

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.

Comments: 6 pages, 4 figures; supplemental information included (14 pages, 10 figures)

Time-resolved fast turbulent dynamo in a laser plasma

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

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

Submitted 24 July, 2020; originally announced July 2020.

Comments: 12 pages, 9 figures

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.

Comments: 26 pages, 6 figures

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

Self-sustaining sound in collisionless, high-beta plasma

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

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

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

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

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.

Comments: Accepted for publication in Phys. Rev. X

Journal ref: Phys. Rev. X 10, 041050 (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.

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

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

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.

Comments: published in Nature Astronomy

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

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

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

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

Abstract: This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood:… ▽ More This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed. △ Less

Submitted 10 March, 2019; originally announced March 2019.

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

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.

Comments: latex, 39 pages, accepted by JPP; minor improvements throughout compared to v1

Journal ref: Journal of Plasma Physics 85, 905850303 (2019)

Magneto-immutable turbulence in weakly collisional plasmas

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

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

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

Comments: Accepted for publication in J. Plasma Phys

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.

Comments: Updated Author and Reviewer Information, 23 pages 17 figures

Journal ref: Astrophys. J. 892:114, 2020

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.

Journal ref: Proceeding of the National Academy of Science 116, 771 (2019)

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.

Comments: 15 pages, 13 figures, submitted to MNRAS

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

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

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

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

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

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.

Comments: 35 pages, minor edits, final version accepted by JPP

Journal ref: Journal of Plasma Physics 84, 905840107 (2018)

Disruption of Alfvénic turbulence by magnetic reconnection in a collisionless plasma

Authors: A. Mallet, A. A. Schekochihin, B. D. G. Chandran

Abstract: We calculate the disruption scale $λ_{\rm D}$ at which sheet-like structures in dynamically aligned Alfvénic turbulence are destroyed by the onset of magnetic reconnection in a low-$β$ collisionless plasma. The scaling of $λ_{\rm D}$ depends on the order of the statistics being considered, with more intense structures being disrupted at larger scales. The disruption scale for the structures that d… ▽ More We calculate the disruption scale $λ_{\rm D}$ at which sheet-like structures in dynamically aligned Alfvénic turbulence are destroyed by the onset of magnetic reconnection in a low-$β$ collisionless plasma. The scaling of $λ_{\rm D}$ depends on the order of the statistics being considered, with more intense structures being disrupted at larger scales. The disruption scale for the structures that dominate the energy spectrum is $λ_{\rm D}\sim L_\perp^{1/9}(d_eρ_s)^{4/9}$, where $d_e$ is the electron inertial scale, $ρ_s$ is the ion sound scale, and $L_\perp$ is the outer scale of the turbulence. When $β_e$ and $ρ_s/L_\perp$ are sufficiently small, the scale $λ_{\rm D}$ is larger than $ρ_s$ and there is a break in the energy spectrum at $λ_{\rm D}$, rather than at $ρ_s$. We propose that the fluctuations produced by the disruption are circularised flux ropes, which may have already been observed in the solar wind. We predict the relationship between the amplitude and radius of these structures and quantify the importance of the disruption process to the cascade in terms of the filling fraction of undisrupted structures and the fractional reduction of the energy contained in them at the ion sound scale $ρ_s$. Both of these fractions depend strongly on $β_e$, with the disrupted structures becoming more important at lower $β_e$. Finally, we predict that the energy spectrum between $λ_{\rm D}$ and $ρ_s$ is steeper than $k_\perp^{-3}$, when this range exists. Such a steep "transition range" is sometimes observed in short intervals of solar-wind turbulence. The onset of collisionless magnetic reconnection may therefore significantly affect the nature of plasma turbulence around the ion gyroscale. △ Less

Submitted 12 October, 2017; v1 submitted 18 July, 2017; originally announced July 2017.

Comments: 21 pages, 1 figure, accepted for publication in JPP

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

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

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

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

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

Laboratory evidence of dynamo amplification of magnetic fields in a turbulent plasma

Authors: P. Tzeferacos, A. Rigby, A. Bott, A. R. Bell, R. Bingham, A. Casner, F. Cattaneo, E. M. Churazov, J. Emig, F. Fiuza, C. B. Forest, J. Foster, C. Graziani, J. Katz, M. Koenig, C. -K. Li, J. Meinecke, R. Petrasso, H. -S. Park, B. A. Remington, J. S. Ross, D. Ryu, D. Ryutov, T. G. White, B. Reville , et al. (5 additional authors not shown)

Abstract: Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $μ$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are… ▽ More Magnetic fields are ubiquitous in the Universe. Extragalactic disks, halos and clusters have consistently been shown, via diffuse radio-synchrotron emission and Faraday rotation measurements, to exhibit magnetic field strengths ranging from a few nG to tens of $μ$G. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. Here we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization of the Universe. △ Less

Submitted 9 February, 2017; originally announced February 2017.

Comments: Submitted for publication, 5 pages 4 figures with supplementary information (25 pages 9 figures)

Numerical modeling of laser-driven experiments aiming to demonstrate magnetic field amplification via turbulent dynamo

Authors: P. Tzeferacos, A. Rigby, A. Bott, A. R. Bell, R. Bingham, A. Casner, F. Cattaneo, E. M. Churazov, J. Emig, N. Flocke, F. Fiuza, C. B. Forest, J. Foster, C. Graziani, J. Katz, M. Koenig, C. -K. Li, J. Meinecke, R. Petrasso, H. -S. Park, B. A. Remington, J. S. Ross, D. Ryu, D. Ryutov, K. Weide , et al. (7 additional authors not shown)

Abstract: The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and… ▽ More The universe is permeated by magnetic fields, with strengths ranging from a femtogauss in the voids between the filaments of galaxy clusters to several teragauss in black holes and neutron stars. The standard model behind cosmological magnetic fields is the nonlinear amplification of seed fields via turbulent dynamo to the values observed. We have conceived experiments that aim to demonstrate and study the turbulent dynamo mechanism in the laboratory. Here we describe the design of these experiments through simulation campaigns using FLASH, a highly capable radiation magnetohydrodynamics code that we have developed, and large-scale three-dimensional simulations on the Mira supercomputer at Argonne National Laboratory. The simulation results indicate that the experimental platform may be capable of reaching a turbulent plasma state and study dynamo amplification. We validate and compare our numerical results with a small subset of experimental data using synthetic diagnostics. △ Less

Submitted 9 February, 2017; originally announced February 2017.

Comments: Accepted for publication on Physics of Plasmas, 15 pages 12 figures

Disruption of sheetlike structures in Alfvénic turbulence by magnetic reconnection

Authors: A. Mallet, A. A. Schekochihin, B. D. G. Chandran

Abstract: We propose a mechanism whereby the intense, sheet-like structures naturally formed by dynamically aligning Alfvénic turbulence are destroyed by magnetic reconnection at a scale $\hatλ_{\rm D}$, larger than the dissipation scale predicted by models of intermittent, dynamically aligning turbulence. The reconnection process proceeds in several stages: first, a linear tearing mode with $N$ magnetic is… ▽ More We propose a mechanism whereby the intense, sheet-like structures naturally formed by dynamically aligning Alfvénic turbulence are destroyed by magnetic reconnection at a scale $\hatλ_{\rm D}$, larger than the dissipation scale predicted by models of intermittent, dynamically aligning turbulence. The reconnection process proceeds in several stages: first, a linear tearing mode with $N$ magnetic islands grows and saturates, and then the $X$-points between these islands collapse into secondary current sheets, which then reconnect until the original structure is destroyed. This effectively imposes an upper limit on the anisotropy of the structures within the perpendicular plane, which means that at scale $\hatλ_{\rm D}$ the turbulent dynamics change: at scales larger than $\hatλ_{\rm D}$, the turbulence exhibits scale-dependent dynamic alignment and a spectral index approximately equal to $-3/2$, while at scales smaller than $\hatλ_{\rm D}$, the turbulent structures undergo a succession of disruptions due to reconnection, limiting dynamic alignment, steepening the effective spectral index and changing the final dissipation scale. The scaling of $\hatλ_{\rm D}$ with the Lundquist (magnetic Reynolds) number $S_{L_\perp}$ depends on the order of the statistics being considered, and on the specific model of intermittency; the transition between the two regimes in the energy spectrum is predicted at approximately $\hatλ_{\rm D} \sim S_{L_\perp}^{-0.6}$. The spectral index below $\hatλ_{\rm D}$ is bounded between $-5/3$ and $-2.3$. The final dissipation scale is at $\hatλ_{η,\infty}\sim S_{L_\perp}^{-3/4}$, the same as the Kolmogorov scale arising in theories of turbulence that do not involve scale-dependent dynamic alignment. △ Less

Submitted 26 July, 2017; v1 submitted 22 December, 2016; originally announced December 2016.

Comments: 10 pages, 4 figures, published in MNRAS 468 p.4862, replaced with published version

Supergranulation and multiscale flows in the solar photosphere: Global observations vs. a theory of anisotropic turbulent convection

Authors: F. Rincon, T. Roudier, A. A. Schekochihin, M. Rieutord

Abstract: The Sun provides us with the only spatially well-resolved astrophysical example of turbulent thermal convection. While various aspects of solar photospheric turbulence, such as granulation (one-Megameter horizontal scale), are well understood, the questions of the physical origin and dynamical organization of larger-scale flows, such as the 30-Megameters supergranulation and flows deep in the sola… ▽ More The Sun provides us with the only spatially well-resolved astrophysical example of turbulent thermal convection. While various aspects of solar photospheric turbulence, such as granulation (one-Megameter horizontal scale), are well understood, the questions of the physical origin and dynamical organization of larger-scale flows, such as the 30-Megameters supergranulation and flows deep in the solar convection zone, remain largely open in spite of their importance for solar dynamics and magnetism. Here, we present a new critical global observational characterization of multiscale photospheric flows and subsequently formulate an anisotropic extension of the Bolgiano-Obukhov theory of hydrodynamic stratified turbulence that may explain several of their distinctive dynamical properties. Our combined analysis suggests that photospheric flows in the horizontal range of scales between supergranulation and granulation have a typical vertical correlation scale of 2.5 to 4 Megameters and operate in a strongly anisotropic, self-similar, nonlinear, buoyant dynamical regime. While the theory remains speculative at this stage, it lends itself to quantitative comparisons with future high-resolution acoustic tomography of subsurface layers and advanced numerical models. Such a validation exercise may also lead to new insights into the asymptotic dynamical regimes in which other, unresolved turbulent anisotropic astrophysical fluid systems supporting waves or instabilities operate. △ Less

Submitted 13 March, 2017; v1 submitted 19 September, 2016; originally announced September 2016.

Comments: 10 pages, 5 figures, revised version published in A&A (high-resolution vector-graphics figures available in the journal version)

Journal ref: A&A 599, A69 (2017)

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

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

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

Submitted 8 June, 2016; originally announced June 2016.

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

A statistical model of three-dimensional anisotropy and intermittency in strong Alfvénic turbulence

Authors: A. Mallet, A. A. Schekochihin

Abstract: We propose a simple statistical model of three-dimensionally anisotropic, intermittent, strong Alfvénic turbulence, incorporating both critical balance and dynamic alignment. Our model is based on log-Poisson statistics for Elsasser-field increments along the magnetic field. We predict the scalings of Elsasser-field conditional two-point structure functions with point separations in all three dire… ▽ More We propose a simple statistical model of three-dimensionally anisotropic, intermittent, strong Alfvénic turbulence, incorporating both critical balance and dynamic alignment. Our model is based on log-Poisson statistics for Elsasser-field increments along the magnetic field. We predict the scalings of Elsasser-field conditional two-point structure functions with point separations in all three directions in a coordinate system locally aligned with the direction of the magnetic field and of the fluctuating fields and obtain good agreement with numerical simulations. We also derive a scaling of the parallel coherence scale of the fluctuations, $l_\parallel \propto λ^{1/2}$, where $λ$ is the perpendicular scale. This is indeed observed for the bulk of the fluctuations in numerical simulations. △ Less

Submitted 14 December, 2016; v1 submitted 24 May, 2016; originally announced June 2016.

Comments: 9 pages, 3 figures, accepted to MNRAS, replaced to match accepted version

A stringent limit on the amplitude of Alfvénic perturbations in high-beta low-collisionality plasmas

Authors: J. Squire, E. Quataert, A. A. Schekochihin

Abstract: It is shown that low-collisionality plasmas cannot support linearly polarized shear-Alfvén fluctuations above a critical amplitude $δB_{\perp}/B_{0} \sim β^{\,-1/2}$, where $β$ is the ratio of thermal to magnetic pressure. Above this cutoff, a developing fluctuation will generate a pressure anisotropy that is sufficient to destabilize itself through the parallel firehose instability. This causes t… ▽ More It is shown that low-collisionality plasmas cannot support linearly polarized shear-Alfvén fluctuations above a critical amplitude $δB_{\perp}/B_{0} \sim β^{\,-1/2}$, where $β$ is the ratio of thermal to magnetic pressure. Above this cutoff, a developing fluctuation will generate a pressure anisotropy that is sufficient to destabilize itself through the parallel firehose instability. This causes the wave frequency to approach zero, interrupting the fluctuation before any oscillation. The magnetic field lines rapidly relax into a sequence of angular zig-zag structures. Such a restrictive bound on shear-Alfvén-wave amplitudes has far-reaching implications for the physics of magnetized turbulence in the high-$β$ conditions prevalent in many astrophysical plasmas, as well as for the solar wind at $\sim 1 \mathrm{AU}$ where $β\gtrsim 1$. △ Less

Submitted 20 September, 2016; v1 submitted 9 May, 2016; originally announced May 2016.

Comments: Accepted for publication in Astrophys. J. Lett

Suppression of thermal conduction in a mirror-unstable plasma

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

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

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

Comments: 11 pages, 11 figures

The nature and energetics of AGN-driven perturbations in the hot gas in the Perseus Cluster

Authors: I. Zhuravleva, E. Churazov, P. Arevalo, A. A. Schekochihin, W. R. Forman, S. W. Allen, A. Simionescu, R. Sunyaev, A. Vikhlinin, N. Werner

Abstract: Cores of relaxed galaxy clusters are often disturbed by AGN. Their Chandra observations revealed a wealth of structures induced by shocks, subsonic gas motions, bubbles of relativistic plasma, etc. In this paper, we determine the nature and energy content of gas fluctuations in the Perseus core by probing statistical properties of emissivity fluctuations imprinted in the soft- and hard-band X-ray… ▽ More Cores of relaxed galaxy clusters are often disturbed by AGN. Their Chandra observations revealed a wealth of structures induced by shocks, subsonic gas motions, bubbles of relativistic plasma, etc. In this paper, we determine the nature and energy content of gas fluctuations in the Perseus core by probing statistical properties of emissivity fluctuations imprinted in the soft- and hard-band X-ray images. About 80 per cent of the total variance of perturbations on ~ 8-70 kpc scales in the inner region have an isobaric nature, i.e., are consistent with slow displacements of the gas in pressure equilibrium with ambient medium. Observed variance translates to the ratio of non-thermal to thermal energy of ~13 per cent. In the region dominated by weak "ripples", about half of the total variance is also associated with isobaric perturbations on scales ~ a few tens of kpc. If these isobaric perturbations are induced by buoyantly rising bubbles, then these results suggest that most of the AGN-injected energy should first go into bubbles rather than into shocks. Using simulations of a shock propagating through the Perseus atmosphere, we found that models reproducing the observed features of a central shock have more than 50 per cent of the AGN-injected energy associated with the bubble enthalpy and only about 20 per cent is carried away with the shock. Such energy partition is consistent with the AGN-feedback model, mediated by bubbles of relativistic plasma, and supports the importance of turbulence in the balance between gas heating and radiative cooling. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: 15 pages, 10 figures, submitted to MNRAS, comments are welcome

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

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

Abstract: The nonlinear state of a high-beta collisionless plasma is investigated when an imposed linear shear amplifies or diminishes a uniform magnetic field, driving pressure anisotropies and hence firehose/mirror instabilities. The evolution of the resulting microscale turbulence is considered when the shear is switched off or reversed after one shear time (mimicking local behaviour of a macroscopic flo… ▽ More The nonlinear state of a high-beta collisionless plasma is investigated when an imposed linear shear amplifies or diminishes a uniform magnetic field, driving pressure anisotropies and hence firehose/mirror instabilities. The evolution of the resulting microscale turbulence is considered when the shear is switched off or reversed after one shear time (mimicking local behaviour of a macroscopic flow), so a new macroscale configuration is superimposed on the microscale state left behind by the previous one. There is a threshold value of plasma beta: when $β\llΩ/S$ (ion cyclotron frequency/shear rate), the emergence of firehose/mirror fluctuations driven unstable by shear and their disappearance when the shear is removed/reversed are quasi-instantaneous compared to the shear time, viz., the decay time of these fluctuations is $\simβ/Ω\ll 1/S$ (this result follows from the free decay of the fluctuations being constrained by the same marginal-stability thresholds as their growth). In contrast, when $β\gtrsimΩ/S$ ("ultra-high" beta), the old microscale state can only be removed on the shear timescale. In this regime, driven firehose fluctuations grow secularly to order-unity amplitudes, compensating for the decay of the mean field and so pinning pressure anisotropy at marginal stability with no appreciable scattering of particles---which is unlike what happens at moderate $β$. When the shear reverses, the shearing away of this firehose turbulence compensates for the increase in the mean field and thus prevents growth of the pressure anisotropy, stopping the system from going mirror-unstable. Therefore, at ultra-high beta, the system stays close to the firehose threshold, the mirror instability is largely suppressed, while the mean magnetic energy barely changes at all. Implications for plasma dynamo and thus the origin of cosmic magnetism are discussed. △ Less

Submitted 26 December, 2015; originally announced December 2015.

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

Turbulent dynamo in a collisionless plasma

Authors: F. Rincon, F. Califano, A. A. Schekochihin, F. Valentini

Abstract: Magnetic fields pervade the entire Universe and affect the formation and evolution of astrophysical systems from cosmological to planetary scales. The generation and dynamical amplification of extragalactic magnetic fields through cosmic times, up to $μ$Gauss levels reported in nearby galaxy clusters, near equipartition with kinetic energy of plasma motions and on scales of at least tens of kilopa… ▽ More Magnetic fields pervade the entire Universe and affect the formation and evolution of astrophysical systems from cosmological to planetary scales. The generation and dynamical amplification of extragalactic magnetic fields through cosmic times, up to $μ$Gauss levels reported in nearby galaxy clusters, near equipartition with kinetic energy of plasma motions and on scales of at least tens of kiloparsecs, is a major puzzle largely unconstrained by observations. A dynamo effect converting kinetic flow energy into magnetic energy is often invoked in that context, however extragalactic plasmas are weakly collisional (as opposed to magnetohydrodynamic fluids), and whether magnetic-field growth and sustainment through an efficient turbulent dynamo instability is possible in such plasmas is not established. Fully kinetic numerical simulations of the Vlasov equation in a six-dimensional phase space necessary to answer this question have until recently remained beyond computational capabilities. Here, we show by means of such simulations that magnetic-field amplification via a dynamo instability does occur in a stochastically-driven, non-relativistic subsonic flow of initially unmagnetized collisionless plasma. We also find that the dynamo self-accelerates and becomes entangled with kinetic instabilities as magnetization increases. The results suggest that such a plasma dynamo may be realizable in laboratory experiments, support the idea that intracluster medium (ICM) turbulence may have significantly contributed to the amplification of cluster magnetic fields up to near-equipartition levels on a timescale shorter than the Hubble time, and emphasize the crucial role of multiscale kinetic physics in high-energy astrophysical plasmas. △ Less

Submitted 1 April, 2016; v1 submitted 20 December, 2015; originally announced December 2015.

Comments: revised version published in PNAS (doi: 10.1073/pnas.1525194113). Author's version pre-journal final editing, 4 pages, 4 figures, 2 movies (url links in manuscript)