Showing 1–34 of 34 results for author: Zhdankin, V

First-principles measurement of ion and electron energization in collisionless accretion flow

Authors: Evgeny A. Gorbunov, Fabio Bacchini, Vladimir Zhdankin, Gregory R. Werner, Mitchell C. Begelman, Dmitri A. Uzdensky

Abstract: We present the largest 3D Particle-in-Cell shearing-box simulations of turbulence driven by the magnetorotational instability, for the first time employing the realistic proton-to-electron mass ratio. We investigate the energy partition between relativistically hot electrons and subrelativistic ions in turbulent accretion flows, a regime relevant to collisionless, radiatively inefficient accretion… ▽ More We present the largest 3D Particle-in-Cell shearing-box simulations of turbulence driven by the magnetorotational instability, for the first time employing the realistic proton-to-electron mass ratio. We investigate the energy partition between relativistically hot electrons and subrelativistic ions in turbulent accretion flows, a regime relevant to collisionless, radiatively inefficient accretion flows of supermassive black holes, such as those targeted by the Event Horizon Telescope. We provide a simple empirical formula to describe the measured heating ratio between ions and electrons, which can be used for more accurate global modeling of accretion flows with standard fluid approaches such as general-relativistic magnetohydrodynamics. △ Less

Submitted 3 October, 2024; originally announced October 2024.

Scale Statistics of Current Sheets in Relativistic Collisionless Plasma Turbulence

Authors: Roberto F. Serrano, Joonas Nättilä, Vladimir Zhdankin

Abstract: We analyze distributions of the spatial scales of coherent intermittent structures -- current sheets -- obtained from fully kinetic, two-dimensional simulations of relativistic plasma turbulence using unsupervised machine-learning data dissection. We find that the distribution functions of sheet length $\ell$ (longest scale of the analyzed structure in the direction perpendicular to the dominant g… ▽ More We analyze distributions of the spatial scales of coherent intermittent structures -- current sheets -- obtained from fully kinetic, two-dimensional simulations of relativistic plasma turbulence using unsupervised machine-learning data dissection. We find that the distribution functions of sheet length $\ell$ (longest scale of the analyzed structure in the direction perpendicular to the dominant guide field) and curvature $r_c$ (radius of a circle fitted to the structures) can be well-approximated by power-law distributions, indicating self-similarity of the structures. The distribution for the sheet width $w$ (shortest scale of the structure) peaks at the kinetic scales and decays exponentially at larger values. The data shows little or no correlation between $w$ and $\ell$, as expected from theoretical considerations. The typical $r_c$ depends linearly on $\ell$, which indicates that the sheets all have a similar curvature relative to their sizes. We find a weak correlation between $r_c$ and $w$. Our results can be used to inform realistic magnetohydrodynamic sub-grid models for plasma turbulence in high-energy astrophysics. △ Less

Submitted 22 August, 2024; originally announced August 2024.

Comments: 19 pages, 10 figures

Collisionless Magnetorotational Turbulence in Pair Plasmas: Steady-state Dynamics, Particle Acceleration, and Radiative Cooling

Authors: Fabio Bacchini, Vladimir Zhdankin, Evgeny A. Gorbunov, Gregory R. Werner, Lev Arzamasskiy, Mitchell C. Begelman, Dmitri A. Uzdensky

Abstract: We present 3D fully kinetic shearing-box simulations of pair-plasma magnetorotational turbulence with unprecedented macro-to-microscopic scale separation. While retrieving the expected fluid behavior of the plasma at large scales, we observe a steepening of turbulent spectra at kinetic scales and substantial angular-momentum transport linked with kinetic processes. For the first time, we provide a… ▽ More We present 3D fully kinetic shearing-box simulations of pair-plasma magnetorotational turbulence with unprecedented macro-to-microscopic scale separation. While retrieving the expected fluid behavior of the plasma at large scales, we observe a steepening of turbulent spectra at kinetic scales and substantial angular-momentum transport linked with kinetic processes. For the first time, we provide a definitive demonstration of nonthermal particle acceleration in kinetic magnetorotational turbulence agnostically of initial conditions, by means of a novel strategy exploiting synchrotron cooling. △ Less

Submitted 14 June, 2024; v1 submitted 2 January, 2024; originally announced January 2024.

Comments: Accepted in PRL

Dimensional measures of generalized entropy

Authors: Vladimir Zhdankin

Abstract: Entropy is useful in statistical problems as a measure of irreversibility, randomness, mixing, dispersion, and number of microstates. However, there remains ambiguity over the precise mathematical formulation of entropy, generalized beyond the additive definition pioneered by Boltzmann, Gibbs, and Shannon (applicable to thermodynamic equilibria). For generalized entropies to be applied rigorously… ▽ More Entropy is useful in statistical problems as a measure of irreversibility, randomness, mixing, dispersion, and number of microstates. However, there remains ambiguity over the precise mathematical formulation of entropy, generalized beyond the additive definition pioneered by Boltzmann, Gibbs, and Shannon (applicable to thermodynamic equilibria). For generalized entropies to be applied rigorously to nonequilibrium statistical mechanics, we suggest that there is a need for a physically interpretable (dimensional) framework that can be connected to dynamical processes operating in phase space. In this work, we introduce dimensional measures of entropy that admit arbitrary invertible weight functions (subject to curvature and convergence requirements). These "dimensional entropies" have physical dimensions of phase-space volume and represent the extent of level sets of the distribution function. Dimensional entropies with power-law weight functions (related to Rényi and Tsallis entropies) are particularly robust, as they do not require any internal dimensional parameters due to their scale invariance. We also point out the existence of composite entropy measures that can be constructed from functionals of dimensional entropies. We calculate the response of the dimensional entropies to perturbations, showing that for a structured distribution, perturbations have the largest impact on entropies weighted at a similar phase-space scale. This elucidates the link between dynamics (perturbations) and statistics (entropies). Finally, we derive corresponding generalized maximum-entropy distributions. Dimensional entropies may be useful as a diagnostic (for irreversibility) and for theoretical modeling (if the underlying irreversible processes in phase space are understood) in chaotic and complex systems, such as collisionless systems of particles with long-range interactions. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: 16 pages, accepted for publication in Journal of Physics A: Math. Theor 2023

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

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

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

Submitted 28 July, 2023; originally announced August 2023.

Comments: 16 pages, 10 figures

Particle acceleration by magnetic Rayleigh-Taylor instability: mechanism for flares in black-hole accretion flows

Authors: Vladimir Zhdankin, Bart Ripperda, Alexander A. Philippov

Abstract: We study the magnetic Rayleigh-Taylor instability in relativistic collisionless plasma, as an astrophysical process for nonthermal particle acceleration. We consider dense plasma on top of a highly magnetized cavity with sheared magnetic field. Using particle-in-cell simulations, we show that small plumes grow and merge progressively to form a large-scale plume, which broadens to drive rapid magne… ▽ More We study the magnetic Rayleigh-Taylor instability in relativistic collisionless plasma, as an astrophysical process for nonthermal particle acceleration. We consider dense plasma on top of a highly magnetized cavity with sheared magnetic field. Using particle-in-cell simulations, we show that small plumes grow and merge progressively to form a large-scale plume, which broadens to drive rapid magnetic reconnection in the cavity. We find that this leads to efficient particle acceleration capable of explaining flares from the inner accretion flow onto the black hole Sgr A*. △ Less

Submitted 10 October, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 10 pages, 7 figures; accepted for publication in Physical Review Research

Journal ref: Phys. Rev. Research 5, 043023 (2023)

Synchrotron Firehose Instability

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

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

Submitted 30 October, 2022; originally announced October 2022.

Comments: 17 pages, 8 figures, submitted for publication

Fully kinetic shearing-box simulations of magnetorotational turbulence in 2D and 3D. I. Pair plasmas

Authors: Fabio Bacchini, Lev Arzamasskiy, Vladimir Zhdankin, Gregory R. Werner, Mitchell C. Begelman, Dmitri A. Uzdensky

Abstract: The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e. kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temper… ▽ More The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e. kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic Particle-in-Cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e. MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron-ion simulations, on which we will focus in future work. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: Submitted to ApJ

Nonthermal particle acceleration from maximum entropy in collisionless plasmas

Authors: Vladimir Zhdankin

Abstract: Dissipative processes cause collisionless plasmas in many systems to develop nonthermal particle distributions with broad power-law tails. The prevalence of power-law energy distributions in space/astrophysical observations and kinetic simulations of systems with a variety of acceleration and trapping (or escape) mechanisms poses a deep mystery. We consider the possibility that such distributions… ▽ More Dissipative processes cause collisionless plasmas in many systems to develop nonthermal particle distributions with broad power-law tails. The prevalence of power-law energy distributions in space/astrophysical observations and kinetic simulations of systems with a variety of acceleration and trapping (or escape) mechanisms poses a deep mystery. We consider the possibility that such distributions can be modeled from maximum-entropy principles, when accounting for generalizations beyond the Boltzmann-Gibbs entropy. Using a dimensional representation of entropy (related to the Renyi and Tsallis entropies), we derive generalized maximum-entropy distributions with a power-law tail determined by the characteristic energy scale at which irreversible dissipation occurs. By assuming that particles are typically energized by an amount comparable to the free energy (per particle) before equilibrating, we derive a formula for the power-law index as a function of plasma parameters for magnetic dissipation in systems with sufficiently complex topologies. The model reproduces several results from kinetic simulations of relativistic turbulence and magnetic reconnection. △ Less

Submitted 29 June, 2022; v1 submitted 24 March, 2022; originally announced March 2022.

Comments: 13 pages, 4 figures, accepted for publication in Journal of Plasma Physics

Journal ref: Journal of Plasma Physics, 88 (3), 175880303 (2022)

Generalized entropy production in collisionless plasma flows and turbulence

Authors: Vladimir Zhdankin

Abstract: Collisionless plasmas exhibit nonthermal particle distributions after being energized; as a consequence, they enter a state of low Boltzmann-Gibbs (BG) entropy relative to the thermal state. The Vlasov equations predict that in a collisionless plasma with closed boundaries, BG entropy is formally conserved, along with an infinite set of Casimir invariants; this provides a seemingly strong constrai… ▽ More Collisionless plasmas exhibit nonthermal particle distributions after being energized; as a consequence, they enter a state of low Boltzmann-Gibbs (BG) entropy relative to the thermal state. The Vlasov equations predict that in a collisionless plasma with closed boundaries, BG entropy is formally conserved, along with an infinite set of Casimir invariants; this provides a seemingly strong constraint that may explain how plasmas maintain low entropy. Nevertheless, it is commonly believed that entropy production is enabled by phase mixing or nonlinear entropy cascades. The question of whether such anomalous entropy production occurs, and of how to characterize it quantitatively, is a fundamental problem in plasma physics. We construct a new theoretical framework for characterizing entropy production (in a generalized sense) based on ideally conserved dimensional quantities derived from the Casimir invariants; these are referred to as ``Casimir momenta'' and they generalize BG entropy. The growth of Casimir momenta relative to the average particle momentum indicates entropy production. We apply this framework to quantify entropy production in particle-in-cell simulations of laminar flows and turbulent flows driven in relativistic plasma, where efficient nonthermal particle acceleration is enabled. We demonstrate that a large amount of anomalous entropy is produced by turbulence despite nonthermal features. These results have implications for reduced modeling of nonthermal particle acceleration and for diagnosing irreversible dissipation in collisionless plasmas such as the solar wind and Earth's magnetosphere. Dimensional representations of generalized entropy analogous to Casimir momenta may be useful for other problems in statistical physics. △ Less

Submitted 7 July, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: 21 pages, 9 figures, accepted for publication in Physical Review X; abstract shortened

Journal ref: Phys. Rev. X 12, 031011 (2022)

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

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

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

Submitted 3 October, 2021; originally announced October 2021.

Comments: 20 pages, 16 figures

Kinetic Simulations of Imbalanced Turbulence in a Relativistic Plasma: Net Flow and Particle Acceleration

Authors: Amelia Hankla, Vladimir Zhdankin, Gregory Werner, Dmitri Uzdensky, Mitchell Begelman

Abstract: Turbulent high-energy astrophysical systems often feature asymmetric energy injection: for instance, Alfven waves propagating from an accretion disk into its corona. Such systems are "imbalanced": the energy fluxes parallel and anti-parallel to the large-scale magnetic field are unequal. In the past, numerical studies of imbalanced turbulence have focused on the magnetohydrodynamic regime. In the… ▽ More Turbulent high-energy astrophysical systems often feature asymmetric energy injection: for instance, Alfven waves propagating from an accretion disk into its corona. Such systems are "imbalanced": the energy fluxes parallel and anti-parallel to the large-scale magnetic field are unequal. In the past, numerical studies of imbalanced turbulence have focused on the magnetohydrodynamic regime. In the present study, we investigate externally-driven imbalanced turbulence in a collisionless, ultrarelativistically hot, magnetized pair plasma using three-dimensional particle-in-cell (PIC) simulations. We find that the injected electromagnetic momentum efficiently converts into plasma momentum, resulting in net motion along the background magnetic field with speeds up to a significant fraction of lightspeed. This discovery has important implications for the launching of accretion disk winds. We also find that although particle acceleration in imbalanced turbulence operates on a slower timescale than in balanced turbulence, it ultimately produces a power-law energy distribution similar to balanced turbulence. Our results have ramifications for black hole accretion disk coronae, winds, and jets. △ Less

Submitted 12 December, 2021; v1 submitted 31 August, 2021; originally announced August 2021.

Comments: 16 pages, accepted to MNRAS

Particle energization in relativistic plasma turbulence: solenoidal versus compressive driving

Authors: Vladimir Zhdankin

Abstract: Many high-energy astrophysical systems contain magnetized collisionless plasmas with relativistic particles, in which turbulence can be driven by an arbitrary mixture of solenoidal and compressive motions. For example, turbulence in hot accretion flows may be driven solenoidally by the magnetorotational instability or compressively by spiral shock waves. It is important to understand the role of t… ▽ More Many high-energy astrophysical systems contain magnetized collisionless plasmas with relativistic particles, in which turbulence can be driven by an arbitrary mixture of solenoidal and compressive motions. For example, turbulence in hot accretion flows may be driven solenoidally by the magnetorotational instability or compressively by spiral shock waves. It is important to understand the role of the driving mechanism on kinetic turbulence and the associated particle energization. In this work, we compare particle-in-cell simulations of solenoidally driven turbulence with similar simulations of compressively driven turbulence. We focus on plasma that has an initial beta of unity, relativistically hot electrons, and varying ion temperature. Apart from strong large-scale density fluctuations in the compressive case, the turbulence statistics are similar for both drives, and the bulk plasma is described reasonably well by an isothermal equation of state. We find that nonthermal particle acceleration is more efficient when turbulence is driven compressively. In the case of relativistically hot ions, both driving mechanisms ultimately lead to similar power-law particle energy distributions, but over a different duration. In the case of non-relativistic ions, there is significant nonthermal particle acceleration only for compressive driving. Additionally, we find that the electron-to-ion heating ratio is less than unity for both drives, but takes a smaller value for compressive driving. We demonstrate that this additional ion energization is associated with the collisionless damping of large-scale compressive modes via perpendicular electric fields. △ Less

Submitted 4 December, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: 29 pages, 28 figures, accepted for publication in ApJ

Journal ref: ApJ 922, 172 (2021)

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

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

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

Submitted 23 July, 2020; originally announced July 2020.

Comments: 8 pages, 5 figures, submitted for publication

Kinetic turbulence in shining pair plasma: intermittent beaming and thermalization by radiative cooling

Authors: Vladimir Zhdankin, Dmitri A. Uzdensky, Gregory R. Werner, Mitchell C. Begelman

Abstract: High-energy astrophysical systems frequently contain collisionless relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an… ▽ More High-energy astrophysical systems frequently contain collisionless relativistic plasmas that are heated by turbulent cascades and cooled by emission of radiation. Understanding the nature of this radiative turbulence is a frontier of extreme plasma astrophysics. In this paper, we use particle-in-cell simulations to study the effects of external inverse Compton radiation on turbulence driven in an optically thin, relativistic pair plasma. We focus on the statistical steady state (where injected energy is balanced by radiated energy) and perform a parameter scan spanning from low magnetization to high magnetization ($0.04 \lesssim σ\lesssim 11$). We demonstrate that the global particle energy distributions are quasi-thermal in all simulations, with only a modest population of nonthermal energetic particles (extending the tail by a factor of $\sim 2$). This indicates that nonthermal particle acceleration (observed in similar non-radiative simulations) is quenched by strong radiative cooling. The quasi-thermal energy distributions are well fit by analytic models in which stochastic particle acceleration (due to, e.g., second-order Fermi mechanism or gyroresonant interactions) is balanced by the radiation reaction force. Despite the efficient thermalization of the plasma, nonthermal energetic particles do make a conspicuous appearance in the anisotropy of the global momentum distribution as highly variable, intermittent beams (for high magnetization cases). The beamed high-energy particles are spatially coincident with intermittent current sheets, suggesting that localized magnetic reconnection may be a mechanism for kinetic beaming. This beaming phenomenon may explain rapid flares observed in various astrophysical systems (such as blazar jets, the Crab nebula, and Sagittarius A*). △ Less

Submitted 27 January, 2020; v1 submitted 21 August, 2019; originally announced August 2019.

Comments: 26 pages, 30 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Extreme Plasma Astrophysics

Authors: D. Uzdensky, M. Begelman, A. Beloborodov, R. Blandford, S. Boldyrev, B. Cerutti, F. Fiuza, D. Giannios, T. Grismayer, M. Kunz, N. Loureiro, M. Lyutikov, M. Medvedev, M. Petropoulou, A. Philippov, E. Quataert, A. Schekochihin, K. Schoeffler, L. Silva, L. Sironi, A. Spitkovsky, G. Werner, V. Zhdankin, J. Zrake, E. Zweibel

Abstract: This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects. This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects. △ Less

Submitted 20 July, 2019; v1 submitted 13 March, 2019; originally announced March 2019.

Comments: A science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. 7 pages including cover page and references. Paper updated in late March 2019 to include a several additional co-authors and references, and a few small changes

First-principles demonstration of diffusive-advective particle acceleration in kinetic simulations of relativistic plasma turbulence

Authors: Kai Wong, Vladimir Zhdankin, Dmitri A. Uzdensky, Gregory R. Werner, Mitchell C. Begelman

Abstract: Nonthermal relativistic plasmas are ubiquitous in astrophysical systems like pulsar wind nebulae and active galactic nuclei, as inferred from their emission spectra. The underlying nonthermal particle acceleration (NTPA) processes have traditionally been modeled with a Fokker-Planck (FP) diffusion-advection equation in momentum space. In this paper, we directly test the FP framework in ab-initio k… ▽ More Nonthermal relativistic plasmas are ubiquitous in astrophysical systems like pulsar wind nebulae and active galactic nuclei, as inferred from their emission spectra. The underlying nonthermal particle acceleration (NTPA) processes have traditionally been modeled with a Fokker-Planck (FP) diffusion-advection equation in momentum space. In this paper, we directly test the FP framework in ab-initio kinetic simulations of driven magnetized turbulence in relativistic pair plasma. By statistically analyzing the motion of tracked particles, we demonstrate the diffusive nature of NTPA and measure the FP energy diffusion ($D$) and advection ($A$) coefficients as functions of particle energy $γm_e c^2$. We find that $D(γ)$ scales as $γ^2$ in the high-energy nonthermal tail, in line with 2nd-order Fermi acceleration theory, but has a much weaker scaling at lower energies. We also find that $A$ is not negligible and reduces NTPA by tending to pull particles towards the peak of the particle energy distribution. This study provides strong support for the FP picture of turbulent NTPA, thereby enhancing our understanding of space and astrophysical plasmas. △ Less

Submitted 22 March, 2020; v1 submitted 10 January, 2019; originally announced January 2019.

Comments: 7 pages, 5 figures, Accepted to Astrophysical Journal Letters

Electron and ion energization in relativistic plasma turbulence

Authors: Vladimir Zhdankin, Dmitri A. Uzdensky, Gregory R. Werner, Mitchell C. Begelman

Abstract: Electron and ion energization (i.e., heating and nonthermal acceleration) is a fundamental, but poorly understood, outcome of plasma turbulence. In this work, we present new results on this topic from particle-in-cell simulations of driven turbulence in collisionless, relativistic electron-ion plasma. We focus on temperatures such that ions (protons) are sub-relativistic and electrons are ultra-re… ▽ More Electron and ion energization (i.e., heating and nonthermal acceleration) is a fundamental, but poorly understood, outcome of plasma turbulence. In this work, we present new results on this topic from particle-in-cell simulations of driven turbulence in collisionless, relativistic electron-ion plasma. We focus on temperatures such that ions (protons) are sub-relativistic and electrons are ultra-relativistic, a regime relevant for high-energy astrophysical systems such as hot accretion flows onto black holes. We find that ions tend to be preferentially heated, gaining up to an order of magnitude more energy than electrons, and propose a simple empirical formula to describe the electron-ion energy partition as a function of the ratio of electron-to-ion gyroradii (which in turn is a function of initial temperatures and plasma beta). We also find that while efficient nonthermal particle acceleration occurs for both species in the ultra-relativistic regime, nonthermal electron populations are diminished with decreasing temperature whereas nonthermal ion populations are essentially unchanged. These results have implications for modeling and interpreting observations of hot accretion flows. △ Less

Submitted 17 January, 2019; v1 submitted 6 September, 2018; originally announced September 2018.

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

Journal ref: Phys. Rev. Lett. 122, 055101 (2019)

System-size convergence of nonthermal particle acceleration in relativistic plasma turbulence

Authors: Vladimir Zhdankin, Dmitri A. Uzdensky, Gregory R. Werner, Mitchell C. Begelman

Abstract: We apply collisionless particle-in-cell simulations of relativistic pair plasmas to explore whether driven turbulence is a viable high-energy astrophysical particle accelerator. We characterize nonthermal particle distributions for varying system sizes up to $L/2πρ_{e0} = 163$, where $L/2π$ is the driving scale and $ρ_{e0}$ is the initial characteristic Larmor radius. We show that turbulent partic… ▽ More We apply collisionless particle-in-cell simulations of relativistic pair plasmas to explore whether driven turbulence is a viable high-energy astrophysical particle accelerator. We characterize nonthermal particle distributions for varying system sizes up to $L/2πρ_{e0} = 163$, where $L/2π$ is the driving scale and $ρ_{e0}$ is the initial characteristic Larmor radius. We show that turbulent particle acceleration produces power-law energy distributions that, when compared at a fixed number of large-scale dynamical times, slowly steepen with increasing system size. We demonstrate, however, that convergence is obtained by comparing the distributions at different times that increase with system size (approximately logarithmically). We suggest that the system-size dependence arises from the time required for particles to reach the highest accessible energies via Fermi acceleration. The converged power-law index of the energy distribution, $α\approx 3.0$ for magnetization $σ= 3/8$, makes turbulence a possible explanation for nonthermal spectra observed in systems such as the Crab nebula. △ Less

Submitted 22 May, 2018; originally announced May 2018.

Comments: 7 pages, 4 figures, submitted to Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, 867, L18 (2018)

Numerical investigation of kinetic turbulence in relativistic pair plasmas I: Turbulence statistics

Authors: Vladimir Zhdankin, Dmitri A. Uzdensky, Gregory R. Werner, Mitchell C. Begelman

Abstract: We describe results from particle-in-cell simulations of driven turbulence in collisionless, magnetized, relativistic pair plasma. This physical regime provides a simple setting for investigating the basic properties of kinetic turbulence and is relevant for high-energy astrophysical systems such as pulsar wind nebulae and astrophysical jets. In this paper, we investigate the statistics of turbule… ▽ More We describe results from particle-in-cell simulations of driven turbulence in collisionless, magnetized, relativistic pair plasma. This physical regime provides a simple setting for investigating the basic properties of kinetic turbulence and is relevant for high-energy astrophysical systems such as pulsar wind nebulae and astrophysical jets. In this paper, we investigate the statistics of turbulent fluctuations in simulations on lattices of up to $1024^3$ cells and containing up to $2 \times 10^{11}$ particles. Due to the absence of a cooling mechanism in our simulations, turbulent energy dissipation reduces the magnetization parameter to order unity within a few dynamical times, causing turbulent motions to become sub-relativistic. In the developed stage, our results agree with predictions from magnetohydrodynamic turbulence phenomenology at inertial-range scales, including a power-law magnetic energy spectrum with index near $-5/3$, scale-dependent anisotropy of fluctuations described by critical balance, log-normal distributions for particle density and internal energy density (related by a $4/3$ adiabatic index, as predicted for an ultra-relativistic ideal gas), and the presence of intermittency. We also present possible signatures of a kinetic cascade by measuring power-law spectra for the magnetic, electric, and density fluctuations at sub-Larmor scales. △ Less

Submitted 14 September, 2017; originally announced September 2017.

Comments: 24 pages, 26 figures, submitted for publication

Journal ref: MNRAS, 474, 2, 2514-2535 (2018)

Influence of a large-scale field on energy dissipation in magnetohydrodynamic turbulence

Authors: Vladimir Zhdankin, Stanislav Boldyrev, Joanne Mason

Abstract: In magnetohydrodynamic (MHD) turbulence, the large-scale magnetic field sets a preferred local direction for the small-scale dynamics, altering the statistics of turbulence from the isotropic case. This happens even in the absence of a total magnetic flux, since MHD turbulence forms randomly oriented large-scale domains of strong magnetic field. It is therefore customary to study small-scale magne… ▽ More In magnetohydrodynamic (MHD) turbulence, the large-scale magnetic field sets a preferred local direction for the small-scale dynamics, altering the statistics of turbulence from the isotropic case. This happens even in the absence of a total magnetic flux, since MHD turbulence forms randomly oriented large-scale domains of strong magnetic field. It is therefore customary to study small-scale magnetic plasma turbulence by assuming a strong background magnetic field relative to the turbulent fluctuations. This is done, for example, in reduced models of plasmas, such as reduced MHD, reduced-dimension kinetic models, gyrokinetics, etc., which make theoretical calculations easier and numerical computations cheaper. Recently, however, it has become clear that the turbulent energy dissipation is concentrated in the regions of strong magnetic field variations. A significant fraction of the energy dissipation may be localized in very small volumes corresponding to the boundaries between strongly magnetized domains. In these regions the reduced models are not applicable. This has important implications for studies of particle heating and acceleration in magnetic plasma turbulence. The goal of this work is to systematically investigate the relationship between local magnetic field variations and magnetic energy dissipation, and to understand its implications for modeling energy dissipation in realistic turbulent plasmas. △ Less

Submitted 9 March, 2017; originally announced March 2017.

Comments: 6 pages, 5 figures, to appear in Monthly Notices of the Royal Astronomical Society

Universal small-scale structure in turbulence driven by magnetorotational instability

Authors: Vladimir Zhdankin, Justin Walker, Stanislav Boldyrev, Geoffroy Lesur

Abstract: The intermittent small-scale structure of turbulence governs energy dissipation in many astrophysical plasmas and is often believed to have universal properties for sufficiently large systems. In this work, we argue that small-scale turbulence in accretion disks is universal in the sense that it is insensitive to the magnetorotational instability (MRI) and background shear, and therefore indisting… ▽ More The intermittent small-scale structure of turbulence governs energy dissipation in many astrophysical plasmas and is often believed to have universal properties for sufficiently large systems. In this work, we argue that small-scale turbulence in accretion disks is universal in the sense that it is insensitive to the magnetorotational instability (MRI) and background shear, and therefore indistinguishable from standard homogeneous magnetohydrodynamic (MHD) turbulence at small scales. We investigate the intermittency of current density, vorticity, and energy dissipation in numerical simulations of incompressible MHD turbulence driven by the MRI in a shearing box. We find that the simulations exhibit a similar degree of intermittency as in standard MHD turbulence. We perform a statistical analysis of intermittent dissipative structures and find that energy dissipation is concentrated in thin sheet-like structures that span a wide range of scales up to the box size. We show that these structures exhibit strikingly similar statistical properties to those in standard MHD turbulence. Additionally, the structures are oriented in the toroidal direction with a characteristic tilt of approximately 17.5 degrees, implying an effective guide field in that direction. △ Less

Submitted 9 February, 2017; originally announced February 2017.

Comments: 9 pages, 11 figures, to appear in Monthly Notices of the Royal Astronomical Society

Kinetic turbulence in relativistic plasma: from thermal bath to non-thermal continuum

Authors: Vladimir Zhdankin, Gregory R. Werner, Dmitri A. Uzdensky, Mitchell C. Begelman

Abstract: We present results from particle-in-cell simulations of driven turbulence in magnetized, collisionless, and relativistic pair plasma. We find that fluctuations are consistent with the classical $k_\perp^{-5/3}$ magnetic energy spectrum at fluid scales and a steeper $k_\perp^{-4}$ spectrum at sub-Larmor scales, where $k_\perp$ is the wavevector perpendicular to the mean field. We demonstrate the de… ▽ More We present results from particle-in-cell simulations of driven turbulence in magnetized, collisionless, and relativistic pair plasma. We find that fluctuations are consistent with the classical $k_\perp^{-5/3}$ magnetic energy spectrum at fluid scales and a steeper $k_\perp^{-4}$ spectrum at sub-Larmor scales, where $k_\perp$ is the wavevector perpendicular to the mean field. We demonstrate the development of a non-thermal, power-law particle energy distribution, $f(E) \sim E^{-α}$, with index $α$ that decreases with increasing magnetization and increases with increasing system size (relative to the characteristic Larmor radius). Our simulations indicate that turbulence can be a viable source of energetic particles in high-energy astrophysical systems, such as pulsar wind nebulae, if scalings asymptotically become insensitive to the system size. △ Less

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

Comments: 6 pages, 4 figures, to appear in PRL

Journal ref: Phys. Rev. Lett. 118, 055103 (2017)

Scalings of intermittent structures in magnetohydrodynamic turbulence

Authors: Vladimir Zhdankin, Stanislav Boldyrev, Dmitri A. Uzdensky

Abstract: Turbulence is ubiquitous in plasmas, leading to rich dynamics characterized by irregularity, irreversibility, energy fluctuations across many scales, and energy transfer across many scales. Another fundamental and generic feature of turbulence, although sometimes overlooked, is the inhomogeneous dissipation of energy in space and in time. This is a consequence of intermittency, the scale-dependent… ▽ More Turbulence is ubiquitous in plasmas, leading to rich dynamics characterized by irregularity, irreversibility, energy fluctuations across many scales, and energy transfer across many scales. Another fundamental and generic feature of turbulence, although sometimes overlooked, is the inhomogeneous dissipation of energy in space and in time. This is a consequence of intermittency, the scale-dependent inhomogeneity of dynamics caused by fluctuations in the turbulent cascade. Intermittency causes turbulent plasmas to self-organize into coherent dissipative structures, which may govern heating, diffusion, particle acceleration, and radiation emissions. In this paper, we present recent progress on understanding intermittency in incompressible magnetohydrodynamic turbulence with a strong guide field. We focus on the statistical analysis of intermittent dissipative structures, which occupy a small fraction of the volume but arguably account for the majority of energy dissipation. We show that, in our numerical simulations, intermittent structures in the current density, vorticity, and Elsässer vorticities all have nearly identical statistical properties. We propose phenomenological explanations for the scalings based on general considerations of Elsässer vorticity structures. Finally, we examine the broader implications of intermittency for astrophysical systems. △ Less

Submitted 16 February, 2016; originally announced February 2016.

Comments: 24 pages, 6 figures, to appear in Physics of Plasmas

Journal ref: Physics of Plasmas 23, 055705 (2016)

Intermittency of Energy Dissipation in Alfvenic Turbulence

Authors: Vladimir Zhdankin, Stanislav Boldyrev, Christopher H. K. Chen

Abstract: We investigate the intermittency of energy dissipation in Alfvenic turbulence by considering the statistics of the coarse-grained energy dissipation rate, using direct measurements from numerical simulations of magnetohydrodynamic turbulence and surrogate measurements from the solar wind. We compare the results to the predictions of the log-normal and log-Poisson random cascade models. We find tha… ▽ More We investigate the intermittency of energy dissipation in Alfvenic turbulence by considering the statistics of the coarse-grained energy dissipation rate, using direct measurements from numerical simulations of magnetohydrodynamic turbulence and surrogate measurements from the solar wind. We compare the results to the predictions of the log-normal and log-Poisson random cascade models. We find that, to a very good approximation, the log-normal model describes the probability density function for the energy dissipation over a broad range of scales, but does not accurately describe the scaling exponents of the moments. The log-Poisson model better describes the scaling exponents of the moments, while the comparison with the probability density function is not straightforward. △ Less

Submitted 22 December, 2015; originally announced December 2015.

Comments: To appear in MNRAS Letters. 5 pages, 4 figures

Journal ref: MNRAS 457, L69-L73 (2016)

Spectral breaks of Alfvenic turbulence in a collisionless plasma

Authors: Stanislav Boldyrev, Christopher H. K. Chen, Qian Xia, Vladimir Zhdankin

Abstract: Recent observations reveal that magnetic turbulence in the nearly colisionless solar wind plasma extends to scales smaller than the plasma microscales, such as ion gyroradius and ion inertial length. Measured breaks in the spectra of magnetic and density fluctuations at high frequencies are thought to be related to the transition from large-scale hydromagnetic to small-scale kinetic turbulence. Th… ▽ More Recent observations reveal that magnetic turbulence in the nearly colisionless solar wind plasma extends to scales smaller than the plasma microscales, such as ion gyroradius and ion inertial length. Measured breaks in the spectra of magnetic and density fluctuations at high frequencies are thought to be related to the transition from large-scale hydromagnetic to small-scale kinetic turbulence. The scales of such transitions and the responsible physical mechanisms are not well understood however. In the present work we emphasize the crucial role of the plasma parameters in the transition to kinetic turbulence, such as the ion and electron plasma beta, the electron to ion temperature ratio, the degree of obliquity of turbulent fluctuations. We then propose an explanation for the spectral breaks reported in recent observations. △ Less

Submitted 1 July, 2015; originally announced July 2015.

Comments: 9 pages, 7 figures. A few typos found in the published version are corrected

Journal ref: The Astrophysical Journal, 806, 238, 2015

Temporal Analysis of Dissipative Structures in Magnetohydrodynamic Turbulence

Authors: Vladimir Zhdankin, Dmitri A. Uzdensky, Stanislav Boldyrev

Abstract: Energy dissipation is highly intermittent in turbulent plasmas, being localized in coherent structures such as current sheets. The statistical analysis of spatial dissipative structures is an effective approach to studying turbulence. In this paper, we generalize this methodology to investigate four-dimensional spatiotemporal structures, i.e., dissipative processes representing sets of interacting… ▽ More Energy dissipation is highly intermittent in turbulent plasmas, being localized in coherent structures such as current sheets. The statistical analysis of spatial dissipative structures is an effective approach to studying turbulence. In this paper, we generalize this methodology to investigate four-dimensional spatiotemporal structures, i.e., dissipative processes representing sets of interacting coherent structures, which correspond to flares in astrophysical systems. We develop methods for identifying and characterizing these processes, and then perform a statistical analysis of dissipative processes in numerical simulations of driven magnetohydrodynamic turbulence. We find that processes are often highly complex, long-lived, and weakly asymmetric in time. They exhibit robust power-law probability distributions and scaling relations, including a distribution of dissipated energy with power-law index near -1.75, indicating that intense dissipative events dominate the overall energy dissipation. We compare our results with the previously observed statistical properties of solar flares. △ Less

Submitted 28 June, 2015; originally announced June 2015.

Comments: Submitted to the Astrophysical Journal. 48 pages, 18 figures

Journal ref: The Astrophysical Journal, 811, 6 (2015)

Temporal Intermittency of Energy Dissipation in Magnetohydrodynamic Turbulence

Authors: Vladimir Zhdankin, Dmitri A. Uzdensky, Stanislav Boldyrev

Abstract: Energy dissipation in magnetohydrodynamic (MHD) turbulence is known to be highly intermittent in space, being concentrated in sheet-like coherent structures. Much less is known about intermittency in time, another fundamental aspect of turbulence which has great importance for observations of solar flares and other space/astrophysical phenomena. In this Letter, we investigate the temporal intermit… ▽ More Energy dissipation in magnetohydrodynamic (MHD) turbulence is known to be highly intermittent in space, being concentrated in sheet-like coherent structures. Much less is known about intermittency in time, another fundamental aspect of turbulence which has great importance for observations of solar flares and other space/astrophysical phenomena. In this Letter, we investigate the temporal intermittency of energy dissipation in numerical simulations of MHD turbulence. We consider four-dimensional spatiotemporal structures, "flare events", responsible for a large fraction of the energy dissipation. We find that although the flare events are often highly complex, they exhibit robust power-law distributions and scaling relations. We find that the probability distribution of dissipated energy has a power law index close to -1.75, similar to observations of solar flares, indicating that intense dissipative events dominate the heating of the system. We also discuss the temporal asymmetry of flare events as a signature of the turbulent cascade. △ Less

Submitted 7 January, 2015; originally announced January 2015.

Comments: To appear in Physical Review Letters. 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 114, 065002 (2015)

Energy dynamics and current sheet structure in fluid and kinetic simulations of decaying magnetohydrodynamic turbulence

Authors: K. D. Makwana, V. Zhdankin, H. Li, W. Daughton, F. Cattaneo

Abstract: Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfvén waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numeric… ▽ More Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfvén waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numerical dissipation whereas the kinetic code has kinetic dissipation. The inertial range power spectrum index is similar in both the codes. The fluid code shows a perpendicular wavenumber spectral slope of $k_{\perp}^{-1.3}$. The kinetic code shows a spectral slope of $k_{\perp}^{-1.5}$ for smaller simulation domain, and $k_{\perp}^{-1.3}$ for larger domain. We estimate that collisionless damping mechanisms in the kinetic code can account for the dissipation of the observed nonlinear energy cascade. Current sheets are geometrically characterized. Their lengths and widths are in good agreement between the two codes. The length scales linearly with the driving scale of the turbulence. In the fluid code, their thickness is determined by the grid resolution as there is no explicit diffusivity. In the kinetic code, their thickness is very close to the skin-depth, irrespective of the grid resolution. This work shows that kinetic codes can reproduce the MHD inertial range dynamics at large scales, while at the same time capturing important kinetic physics at small scales. △ Less

Submitted 1 June, 2015; v1 submitted 15 December, 2014; originally announced December 2014.

Journal ref: Physics of Plasmas, vol. 22, 042902 (2015)

Energy Dissipation in Magnetohydrodynamic Turbulence: Coherent Structures or "Nanoflares"?

Authors: Vladimir Zhdankin, Stanislav Boldyrev, Jean Carlos Perez, Steven M. Tobias

Abstract: We investigate the intermittency of energy dissipation in magnetohydrodynamic (MHD) turbulence by identifying dissipative structures and measuring their characteristic scales. We find that the probability distribution of energy dissipation rates exhibits a power law tail with index very close to the critical value of -2.0, which indicates that structures of all intensities contribute equally to en… ▽ More We investigate the intermittency of energy dissipation in magnetohydrodynamic (MHD) turbulence by identifying dissipative structures and measuring their characteristic scales. We find that the probability distribution of energy dissipation rates exhibits a power law tail with index very close to the critical value of -2.0, which indicates that structures of all intensities contribute equally to energy dissipation. We find that energy dissipation is uniformly spread among coherent structures with lengths and widths in the inertial range. At the same time, these structures have thicknesses deep within the dissipative regime. As the Reynolds number is increased, structures become thinner and more numerous, while the energy dissipation continues to occur mainly in large-scale coherent structures. This implies that in the limit of high Reynolds number, energy dissipation occurs in thin, tightly packed current sheets which nevertheless span a continuum of scales up to the system size, exhibiting features of both coherent structures and nanoflares previously conjectured as a coronal heating mechanism. △ Less

Submitted 10 September, 2014; originally announced September 2014.

Comments: To appear in The Astrophysical Journal. 12 pages, 10 figures

Journal ref: The Astrophysical Journal, 795, 127 (2014)

Statistical Analysis of Current Sheets in Three-Dimensional Magnetohydrodynamic Turbulence

Authors: Vladimir Zhdankin, Dmitri A. Uzdensky, Jean C. Perez, Stanislav Boldyrev

Abstract: We develop a framework for studying the statistical properties of current sheets in numerical simulations of 3D magnetohydrodynamic (MHD) turbulence. We describe an algorithm that identifies current sheets in a simulation snapshot and then determines their geometrical properties (including length, width, and thickness) and intensities (peak current density and total energy dissipation rate). We th… ▽ More We develop a framework for studying the statistical properties of current sheets in numerical simulations of 3D magnetohydrodynamic (MHD) turbulence. We describe an algorithm that identifies current sheets in a simulation snapshot and then determines their geometrical properties (including length, width, and thickness) and intensities (peak current density and total energy dissipation rate). We then apply this procedure to simulations of reduced MHD turbulence and perform a statistical analysis on the obtained population of current sheets. We evaluate the role of reconnection by separately studying the populations of current sheets which contain magnetic X-points and those which do not. We find that the statistical properties of the two populations are different in general. We compare the scaling of these properties to phenomenological predictions obtained for the inertial range of MHD turbulence. Finally, we test whether the reconnecting current sheets are consistent with the Sweet-Parker model. △ Less

Submitted 6 February, 2013; originally announced February 2013.

Comments: 19 pages, 19 figures

Journal ref: The Astrophysical Journal, 771, 124 (2013)

Distribution of Magnetic Discontinuities in the Solar Wind and in MHD Turbulence

Authors: Vladimir Zhdankin, Stanislav Boldyrev, Joanne Mason

Abstract: The statistical properties of magnetic discontinuities in the solar wind are investigated by measuring fluctuations in the magnetic field direction, given by the rotation Delta theta that the magnetic field vector undergoes during time interval Delta t. We show that the probability density function for rotations, P(Delta theta), can be described by a simple model in which the magnetic field vector… ▽ More The statistical properties of magnetic discontinuities in the solar wind are investigated by measuring fluctuations in the magnetic field direction, given by the rotation Delta theta that the magnetic field vector undergoes during time interval Delta t. We show that the probability density function for rotations, P(Delta theta), can be described by a simple model in which the magnetic field vector rotates with a relative increment (Delta B)/B that is lognormally distributed. We find that the probability density function of increments, P((Delta B)/B), has a remarkable scaling property: the normalized variable x=[(Delta B)/B]*[(Delta t)/(Delta t_0)]^-a has a universal lognormal distribution for all time intervals Delta t. We then compare measurements from the solar wind with those from direct numerical simulations of magnetohydrodynamic (MHD) turbulence. We find good agreement for P(Delta theta) obtained in the two cases when the magnetic guide-field to fluctuations ratio B_0/b_rms is chosen accordingly. However, the scale invariance of P((Delta B)/B) is broken in the MHD simulations with relatively limited inertial interval, which causes P(Delta theta) to scale with measurement interval differently than in the solar wind. △ Less

Submitted 12 October, 2012; originally announced October 2012.

Comments: To appear in Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, 760, L22 (2012)

Magnetic Discontinuities in Magnetohydrodynamic Turbulence and in the Solar Wind

Authors: Vladimir Zhdankin, Stanislav Boldyrev, Joanne Mason, Jean Carlos Perez

Abstract: Recent measurements of solar wind turbulence report the presence of intermittent, exponentially distributed angular discontinuities in the magnetic field. In this Letter, we study whether such discontinuities can be produced by magnetohydrodynamic (MHD) turbulence. We detect the discontinuities by measuring the fluctuations of the magnetic field direction, Delta theta, across fixed spatial increme… ▽ More Recent measurements of solar wind turbulence report the presence of intermittent, exponentially distributed angular discontinuities in the magnetic field. In this Letter, we study whether such discontinuities can be produced by magnetohydrodynamic (MHD) turbulence. We detect the discontinuities by measuring the fluctuations of the magnetic field direction, Delta theta, across fixed spatial increments Delta x in direct numerical simulations of MHD turbulence with an imposed uniform guide field B_0. A large region of the probability density function (pdf) for Delta theta is found to follow an exponential decay, proportional to exp(-Delta theta/theta_*), with characteristic angle theta_* ~ (14 deg) (b_rms/B_0)^0.65 for a broad range of guide-field strengths. We find that discontinuities observed in the solar wind can be reproduced by MHD turbulence with reasonable ratios of b_rms/B_0. We also observe an excess of small angular discontinuities when Delta x becomes small, possibly indicating an increasing statistical significance of dissipation-scale structures. The structure of the pdf in this case closely resembles the two-population pdf seen in the solar wind. We thus propose that strong discontinuities are associated with inertial-range MHD turbulence, while weak discontinuities emerge from near-dissipation-range turbulence. In addition, we find that the structure functions of the magnetic field direction exhibit anomalous scaling exponents, which indicates the existence of intermittent structures. △ Less

Submitted 19 April, 2012; originally announced April 2012.

Comments: To appear in Physical Review Letters

Journal ref: Phys. Rev. Lett. 108, 175004 (2012)

Residual energy in magnetohydrodynamic turbulence and in the solar wind

Authors: Stanislav Boldyrev, Jean Carlos Perez, Vladimir Zhdankin

Abstract: Recent observations indicate that kinetic and magnetic energies are not in equipartition in the solar wind turbulence. Rather, magnetic fluctuations are more energetic and have somewhat steeper energy spectrum compared to the velocity fluctuations. This leads to the presence of the so-called residual energy E_r=E_v-E_b in the inertial interval of turbulence. This puzzling effect is addressed in th… ▽ More Recent observations indicate that kinetic and magnetic energies are not in equipartition in the solar wind turbulence. Rather, magnetic fluctuations are more energetic and have somewhat steeper energy spectrum compared to the velocity fluctuations. This leads to the presence of the so-called residual energy E_r=E_v-E_b in the inertial interval of turbulence. This puzzling effect is addressed in the present paper in the framework of weak turbulence theory. Using a simple model of weakly colliding Alfvén waves, we demonstrate that the kinetic-magnetic equipartition indeed gets broken as a result of nonlinear interaction of Alfvén waves. We establish that magnetic energy is indeed generated more efficiently as a result of these interactions, which proposes an explanation for the solar wind observations. △ Less

Submitted 30 August, 2011; originally announced August 2011.

Comments: To appear in "Physics of the Heliosphere: a 10-year Retrospective", Jacob Heerikhuisen, Gang Li, and Gary P. Zank, Eds