Showing 1–50 of 74 results for author: Uzdensky, D

X-ray imaging and electron temperature evolution in laser-driven magnetic reconnection experiments at the National Ignition Facility

Authors: V. Valenzuela-Villaseca, J. M. Molina, D. B. Schaeffer, S. Malko, J. Griff-McMahon, K. Lezhnin, M. J. Rosenberg, S. X. Hu, D. Kalantar, C. Trosseille, H. -S. Park, B. A. Remington, G. Fiksel, D. Uzdensky, A. Bhattacharjee, W. Fox

Abstract: We present results from X-ray imaging of high-aspect-ratio magnetic reconnection experiments driven at the National Ignition Facility. Two parallel, self-magnetized, elongated laser-driven plumes are produced by tiling 40 laser beams. A magnetic reconnection layer is formed by the collision of the plumes. A gated X-ray framing pinhole camera with micro-channel plate (MCP) detector produces multipl… ▽ More We present results from X-ray imaging of high-aspect-ratio magnetic reconnection experiments driven at the National Ignition Facility. Two parallel, self-magnetized, elongated laser-driven plumes are produced by tiling 40 laser beams. A magnetic reconnection layer is formed by the collision of the plumes. A gated X-ray framing pinhole camera with micro-channel plate (MCP) detector produces multiple images through various filters of the formation and evolution of both the plumes and current sheet. As the diagnostic integrates plasma self-emission along the line of sight, 2-dimensional electron temperature maps $\langle T_e \rangle_Y$ are constructed by taking the ratio of intensity of these images obtained with different filters. The plumes have a characteristic temperature $\langle T_e \rangle_Y = 240 \pm 20$ eV at 2 ns after the initial laser irradiation and exhibit a slow cooling up to 4 ns. The reconnection layer forms at 3 ns with a temperature $\langle T_e \rangle_Y = 280 \pm 50$ eV as the result of the collision of the plumes. The error bars of the plumes and current sheet temperatures separate at $4$ ns, showing the heating of the current sheet from colder inflows. Using a semi-analytical model, we find that the observed heating of the current sheet is consistent with being produced by electron-ion drag, rather than the conversion of magnetic to kinetic energy. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: Submitted to Physics of Plasmas. 19 pages (total), 14 figures, 2 tables

Radiatively Cooled Magnetic Reconnection Experiments Driven by Pulsed Power

Authors: R Datta, K Chandler, C E Myers, J P Chittenden, A J Crilly, C Aragon, D J Ampleford, J T Banasek, A Edens, W R Fox, S B Hansen, E C Harding, C A Jennings, H Ji, C C Kuranz, S V Lebedev, Q Looker, S G Patel, A J Porwitzky, G A Shipley, D A Uzdensky, D A Yager-Elorriaga, J D Hare

Abstract: We present evidence for strong radiative cooling in a pulsed-power-driven magnetic reconnection experiment. Two aluminum exploding wire arrays, driven by a 20 MA peak current, 300 ns rise time pulse from the Z machine (Sandia National Laboratories), generate strongly-driven plasma flows ($M_A \approx 7$) with anti-parallel magnetic fields, which form a reconnection layer ($S_L \approx 120$) at the… ▽ More We present evidence for strong radiative cooling in a pulsed-power-driven magnetic reconnection experiment. Two aluminum exploding wire arrays, driven by a 20 MA peak current, 300 ns rise time pulse from the Z machine (Sandia National Laboratories), generate strongly-driven plasma flows ($M_A \approx 7$) with anti-parallel magnetic fields, which form a reconnection layer ($S_L \approx 120$) at the mid-plane. The net cooling rate far exceeds the Alfvénic transit rate ($τ_{\text{cool}}^{-1}/τ_{\text{A}}^{-1} > 100$), leading to strong cooling of the reconnection layer. We determine the advected magnetic field and flow velocity using inductive probes positioned in the inflow to the layer, and inflow ion density and temperature from analysis of visible emission spectroscopy. A sharp decrease in X-ray emission from the reconnection layer, measured using filtered diodes and time-gated X-ray imaging, provides evidence for strong cooling of the reconnection layer after its initial formation. X-ray images also show localized hotspots, regions of strong X-ray emission, with velocities comparable to the expected outflow velocity from the reconnection layer. These hotspots are consistent with plasmoids observed in 3D radiative resistive magnetohydrodynamic simulations of the experiment. X-ray spectroscopy further indicates that the hotspots have a temperature (170 eV) much higher than the bulk layer ($\leq$ 75 eV) and inflow temperatures (about 2 eV), and that these hotspots generate the majority of the high-energy (> 1 keV) emission. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Plasmoid formation and strong radiative cooling in a driven magnetic reconnection experiment

Authors: R. Datta, K. Chandler, C. E. Myers, J. P. Chittenden, A. J. Crilly, C. Aragon, D. J. Ampleford, J. T. Banasek, A. Edens, W. R. Fox, S. B. Hansen, E. C. Harding, C. A. Jennings, H. Ji, C. C. Kuranz, S. V. Lebedev, Q. Looker, S. G. Patel, A. Porwitzky, G. A. Shipley, D. A. Uzdensky, D. A. Yager-Elorriaga, J. D. Hare

Abstract: We present results from the first experimental study of strongly radiatively-cooled magnetic reconnection. Two exploding aluminum wire arrays, driven simultaneously by the Z machine ($I_{max} = 20 \, \text{MA}$, $t_{\text{rise}} = 300 \, \text{ns}$), generate a radiatively-cooled reconnection layer ($S_L \approx 120$) in which the total cooling rate exceeds the hydrodynamic transit rate (… ▽ More We present results from the first experimental study of strongly radiatively-cooled magnetic reconnection. Two exploding aluminum wire arrays, driven simultaneously by the Z machine ($I_{max} = 20 \, \text{MA}$, $t_{\text{rise}} = 300 \, \text{ns}$), generate a radiatively-cooled reconnection layer ($S_L \approx 120$) in which the total cooling rate exceeds the hydrodynamic transit rate ($τ_{\text{hydro}}/τ_{\text{cool}} > 100$). Measurements of X-ray emission from the reconnection layer using a filtered diode ($>1$ keV) show a narrow (50 ns FWHM) burst of emission at 220 ns after current start, consistent with the formation and subsequent rapid cooling of the reconnection layer. Time-gated X-ray images of the reconnection layer show fast-moving (up to 50 km/s) hotspots inside the layer, consistent with the presence of plasmoids observed in 3D resistive magnetohydrodynamic simulations. X-ray spectroscopy shows that these hotspots generate the majority of Al K-shell emission (at around 1.6 keV) prior to the onset of cooling, and exhibit temperatures of 170 eV, much greater than the temperature of the plasma inflows and the rest of the reconnection layer. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Simulations of Radiatively Cooled Magnetic Reconnection Driven by Pulsed Power

Authors: Rishabh Datta, Aidan J. Crilly, Jeremy P. Chittenden, Simran Chowdhry, Katherine Chandler, Nikita Chaturvedi, Clayton E. Myers, William R. Fox, Stephanie B. Hansen, Christopher A. Jennings, Hantao Ji, Carolyn C. Kuranz, Sergey V. Lebedev, Dmitri A. Uzdensky, Jack D. Hare

Abstract: Magnetic reconnection is an important process in astrophysical environments, as it re-configures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of th… ▽ More Magnetic reconnection is an important process in astrophysical environments, as it re-configures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the layer. In this paper, we perform 2D & 3D simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively-cooled regime. These simulations are performed in GORGON, an Eulerian resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P$_{1/3}$ multi-group radiation transport. The load produces highly collisional, super-Alfvénic $(M_{A} \approx 1.5)$, supersonic $(M_S \approx 4-5)$ plasma flows which generate a reconnection layer ($L/δ \approx 100, S_L \approx 400$). The reconnection layer undergoes radiative collapse when the radiative losses exceed Ohmic and compressional heating $τ_{cool}^{-1}/τ_A^{-1} \approx 100$; this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Finally, the current sheet is unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer. △ Less

Submitted 3 January, 2024; originally announced January 2024.

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

Electron and proton energization in 3D reconnecting current sheets in semirelativistic plasma with guide magnetic field

Authors: Gregory R. Werner, Dmitri A. Uzdensky

Abstract: Using 3D particle-in-cell simulation, we characterize energy conversion, as a function of guide magnetic field, in a thin current sheet in semirelativistic plasma, with relativistic electrons and subrelativistic protons. There, magnetic reconnection, the drift-kink instability (DKI), and the flux-rope kink instability all compete and interact in their nonlinear stages to convert magnetic energy to… ▽ More Using 3D particle-in-cell simulation, we characterize energy conversion, as a function of guide magnetic field, in a thin current sheet in semirelativistic plasma, with relativistic electrons and subrelativistic protons. There, magnetic reconnection, the drift-kink instability (DKI), and the flux-rope kink instability all compete and interact in their nonlinear stages to convert magnetic energy to plasma energy. We compare fully 3D simulations with 2D in two different planes to isolate reconnection and DKI effects. In zero guide field, these processes yield distinct energy conversion signatures: ions gain more energy than electrons in 2Dxy (reconnection), while the opposite is true in 2Dyz (DKI), and the 3D result falls in between. The flux-rope instability, which occurs only in 3D, allows more magnetic energy to be released than in 2D, but the rate of energy conversion in 3D tends to be lower. Increasing the guide magnetic field strongly suppresses DKI, and in all cases slows and reduces the overall amount of energy conversion; it also favors electron energization through a process by which energy is first stored in the motional electric field of flux ropes before energizing particles. Understanding the evolution of the energy partition thus provides insight into the role of various plasma processes, and is important for modeling radiation from astrophysical sources such as accreting black holes and their jets. △ Less

Submitted 5 April, 2024; v1 submitted 24 November, 2023; originally announced November 2023.

Comments: 10 pages, 3 figures

Journal ref: Astrophys. J. Lett. 964(2), L21 (2024)

Magnetic field generation in laser-solid interactions at strong-field QED relevant intensities

Authors: Brandon K. Russell, Marija Vranic, Paul T. Campbell, Alexander G. R. Thomas, Kevin M. Schoeffler, Dmitri A. Uzdensky, Louise Willingale

Abstract: Magnetic field generation in ultra-intense laser-solid interactions is studied over a range of laser intensities relevant to next-generation laser facilities ($a_0 = 50-500$) using 2D particle-in-cell simulations. It is found that fields on the order of 0.1 MT (1 GigaGauss) may be generated by relativistic electrons traveling along the surface of the target. However a significant fraction of the e… ▽ More Magnetic field generation in ultra-intense laser-solid interactions is studied over a range of laser intensities relevant to next-generation laser facilities ($a_0 = 50-500$) using 2D particle-in-cell simulations. It is found that fields on the order of 0.1 MT (1 GigaGauss) may be generated by relativistic electrons traveling along the surface of the target. However a significant fraction of the energy budget is converted to high-energy photons, ~38% at $a_0=500$, greatly reducing the available energy for field generation. A model for the evolution of the target-surface fields and their scaling with $a_0$ is developed using laser parameters and assumed values for the average radial electron velocity and reflectivity. The model and empirical scaling allow for the estimation of field strengths on the next generation of laser facilities, a necessary component to the proposal of any future magnetized experiment. △ Less

Submitted 10 September, 2023; originally announced September 2023.

Comments: 8 pages, 8 figures

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

Kinetic simulations and gamma-ray signatures of Klein-Nishina relativistic magnetic reconnection

Authors: J. Mehlhaff, G. Werner, B. Cerutti, D. Uzdensky, M. Begelman

Abstract: Black hole and neutron star environments often comprise collisionless plasmas immersed in strong magnetic fields and intense baths of low-frequency radiation. In such conditions, relativistic magnetic reconnection can tap the magnetic field energy, accelerating high-energy particles that rapidly cool by inverse Compton (IC) scattering the dense photon background. At the highest particle energies r… ▽ More Black hole and neutron star environments often comprise collisionless plasmas immersed in strong magnetic fields and intense baths of low-frequency radiation. In such conditions, relativistic magnetic reconnection can tap the magnetic field energy, accelerating high-energy particles that rapidly cool by inverse Compton (IC) scattering the dense photon background. At the highest particle energies reached in bright gamma-ray sources, IC scattering can stray into the Klein-Nishina regime. Here, the Comptonized photons exceed pair-production threshold with the radiation background and may thus return their energy to the reconnecting plasma as fresh electron-positron pairs. To reliably characterize observable signatures of such Klein-Nishina reconnection, in this work, we present first-principles particle-in-cell simulations of pair-plasma relativistic reconnection coupled to Klein-Nishina and pair-production physics. The simulations show substantial differences between the observable signatures of Klein-Nishina reconnection and reconnection coupled only to low-energy Thomson IC cooling (without pair production). The latter regime exhibits strong harder-when-brighter behaviour; the former involves a stable spectral shape independent of overall brightness. This spectral stability is reminiscent of flat-spectrum radio quasar (FSRQ) GeV high states, furnishing evidence that Klein-Nishina radiative physics operates in FSRQs. The simulated Klein-Nishina reconnection pair yield spans from low to order-unity and follows an exponential scaling law in a single governing parameter. Pushing this parameter beyond its range studied here might give way to a copious pair-creation regime. Besides FSRQs, we discuss potential applications to accreting black hole X-ray binaries, the M87$^*$ magnetosphere, and gamma-ray binaries. △ Less

Submitted 25 June, 2023; originally announced June 2023.

Comments: 38 pages, 19 figures, submitted to MNRAS

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

Authors: K. M. Schoeffler, T. Grismayer, D. Uzdensky, L. O. Silva

Abstract: The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with two- and three-dimensional (2D and 3D) particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiat… ▽ More The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with two- and three-dimensional (2D and 3D) particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field $B$ and plasma density~$n$, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the $n\mbox{-}B$ space are developed and used to visualize and quantify the effects of compression. The $n\mbox{-}B$ histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic (QED) effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres. △ Less

Submitted 29 March, 2023; originally announced March 2023.

Comments: 31 pages, 23 figures

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

Particle Injection and Nonthermal Particle Acceleration in Relativistic Magnetic Reconnection

Authors: Omar French, Fan Guo, Qile Zhang, Dmitri Uzdensky

Abstract: Magnetic reconnection in the relativistic regime has been proposed as an important process for the efficient production of nonthermal particles and high-energy emissions. Using fully kinetic particle-in-cell simulations, we investigate how guide-field strength and domain size affect characteristic spectral features and acceleration processes. We study two stages of acceleration: energization up un… ▽ More Magnetic reconnection in the relativistic regime has been proposed as an important process for the efficient production of nonthermal particles and high-energy emissions. Using fully kinetic particle-in-cell simulations, we investigate how guide-field strength and domain size affect characteristic spectral features and acceleration processes. We study two stages of acceleration: energization up until the injection energy $γ_{\rm inj}$ and further acceleration that generates a power-law spectrum. Stronger guide fields increase the power-law index and $γ_{\rm inj}$, which suppresses acceleration efficiency. These quantities seemingly converge with increasing domain size, suggesting that our findings can be extended to large-scale systems. We find that three distinct mechanisms contribute to acceleration during injection: particle streaming along the parallel electric field, Fermi reflection, and the pickup process. Fermi and pickup processes, related to the electric field perpendicular to the magnetic field, govern the injection for weak guide fields and larger domains. Meanwhile, parallel electric fields are important for injection in the strong guide field regime. In the post-injection stage, we find that perpendicular electric fields dominate particle acceleration in the weak guide field regime, whereas parallel electric fields control acceleration for strong guide fields. These findings will help explain the nonthermal acceleration and emissions in high-energy astrophysics, including black hole jets and pulsar wind nebulae. △ Less

Submitted 1 March, 2023; v1 submitted 15 October, 2022; originally announced October 2022.

Comments: 27 pages, 15 figures, accepted for publication to The Astrophysical Journal

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

A novel kinetic mechanism for the onset of fast magnetic reconnection and plasmoid instabilities in collisionless plasmas

Authors: W. Fox, G. Fiksel, D. B. Schaeffer, D. Haberberger, J. Matteucci, K. Lezhnin, A. Bhattacharjee, M. J. Rosenberg, S. X. Hu, A. Howard, D. Uzdensky, K. Germaschewski

Abstract: Magnetic reconnection can explosively release magnetic energy when opposing magnetic fields merge and annihilate through a current sheet, driving plasma jets and accelerating non-thermal particle populations to high energy, in plasmas ranging from space and astrophysical to laboratory scales. Through laboratory experiments and spacecraft observations, significant experimental progress has been mad… ▽ More Magnetic reconnection can explosively release magnetic energy when opposing magnetic fields merge and annihilate through a current sheet, driving plasma jets and accelerating non-thermal particle populations to high energy, in plasmas ranging from space and astrophysical to laboratory scales. Through laboratory experiments and spacecraft observations, significant experimental progress has been made in demonstrating how fast dissipation and reconnection occurs in narrow, kinetic-scale current sheets. However, a challenge has been to demonstrate what triggers reconnection and how it proceeds rapidly and efficiently as part of a global system much larger than these kinetic scales. Here we show experimentally the full development of a process where the current sheet forms and then breaks up into multiple current-carrying structures at the ion kinetic scale. The results are consistent with tearing of the current sheet, however modified by collisionless kinetic ion effects, which leads to a larger growth rate and number of plasmoids than observed in previous experiments or compared to predictions from standard tearing instability theory and previous non-linear kinetic reconnection simulations. This effect will increase the role of plasmoid instabilities in many natural reconnection systems and should be considered in triggering rapid reconnection in a broad range of natural plasmas with collisionless, compressible flows, including at the Earth's magnetosheath and magnetotail and at the heliopause, in accretion disks, and in turbulent high-Mach-number collisionless shocks. △ Less

Submitted 6 December, 2021; originally announced December 2021.

Comments: 18 pages, 4 figures, with Supplementary materials

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

Pair-Regulated Klein-Nishina Relativistic Magnetic Reconnection with Applications to Blazars and Accreting Black Holes

Authors: J. M. Mehlhaff, G. R. Werner, D. A. Uzdensky, M. C. Begelman

Abstract: Relativistic magnetic reconnection is a powerful agent through which magnetic energy can be tapped in astrophysics, energizing particles that then produce observed radiation. In some systems, the highest energy photons come from particles Comptonizing an ambient radiation bath supplied by an external source. If the emitting particle energies are high enough, this inverse Compton (IC) scattering en… ▽ More Relativistic magnetic reconnection is a powerful agent through which magnetic energy can be tapped in astrophysics, energizing particles that then produce observed radiation. In some systems, the highest energy photons come from particles Comptonizing an ambient radiation bath supplied by an external source. If the emitting particle energies are high enough, this inverse Compton (IC) scattering enters the Klein-Nishina regime, which differs from the low-energy Thomson IC limit in two significant ways. First, radiative losses become inherently discrete, with particles delivering an order-unity fraction of their energies to single photons. Second, Comptonized photons may pair-produce with the ambient radiation, opening up another channel for radiative feedback on magnetic reconnection. We analytically study externally illuminated highly magnetized reconnecting systems for which both of these effects are important. We identify a universal (initial magnetization-independent) quasi-steady state in which gamma-rays emitted from the reconnection layer are absorbed in the upstream region, and the resulting hot pairs dominate the energy density of the inflow plasma. However, a true pair cascade is unlikely, and the number density of created pairs remains subdominant to that of the original plasma for a wide parameter range. Future particle-in-cell simulation studies may test various aspects. Pair-regulated Klein-Nishina reconnection may explain steep spectra (quiescent and flaring) from flat-spectrum radio quasars and black hole accretion disc coronae. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 43 pages, 22 figures, submitted to MNRAS

Reconnection and particle acceleration in 3D current sheet evolution in moderately-magnetized astrophysical pair plasma

Authors: Gregory R. Werner, Dmitri A. Uzdensky

Abstract: Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially 2D) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how thi… ▽ More Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially 2D) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in 3D, we conduct a comprehensive particle-in-cell simulation study comparing 2D and 3D evolution of long, thin current sheets in moderately-magnetized, collisionless, relativistically-hot electron-positron plasma, and find dramatic differences. We first systematically characterize this process in 2D, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths, and system sizes. We then show that 3D simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3D current-sheet evolution is not always fundamentally classical reconnection with perturbing 3D effects, but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details, and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in 2D. Intriguingly, nonthermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in 2D and 3D. Though the variety of underlying current-sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable. △ Less

Submitted 15 December, 2021; v1 submitted 5 June, 2021; originally announced June 2021.

Comments: 92 pages, 42 figures, 4 tables

Journal ref: Journal of Plasma Physics 87(6), 905870613 (2021)

Statistical description of coalescing magnetic islands via magnetic reconnection

Authors: Muni Zhou, David H. Wu, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath, and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral… ▽ More The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath, and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in 2D. We solve our island kinetic equation numerically for three different types of initial distribution: delta-distribution, Gaussian and power-law distribution. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde t^{-1}$, the number of islands decreases as $\tilde t^{-1}$, and the averaged area of islands grows as $\tilde t$, where $\tilde t$ is the time normalized to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. △ Less

Submitted 30 October, 2021; v1 submitted 28 April, 2021; originally announced April 2021.

Comments: 38 pages; 10 figures; Accepted by Journal of Plasma Physics

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Authors: H. Ji, J. Karpen, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, A. Bhattacharjee, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, B. Chen, L. -J. Chen, Y. Chen, A. Chien, L. Comisso, D. Craig, J. Dahlin, W. Daughton, E. DeLuca , et al. (83 additional authors not shown)

Abstract: Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagen… ▽ More Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

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

Relativistic Nonthermal Particle Acceleration in Two-Dimensional Collisionless Magnetic Reconnection

Authors: Dmitri A. Uzdensky

Abstract: Magnetic reconnection, especially in the relativistic regime, provides an efficient mechanism for accelerating relativistic particles and thus offers an attractive physical explanation for nonthermal high-energy emission from various astrophysical sources. I present a simple analytical model that elucidates key physical processes responsible for reconnection-driven relativistic nonthermal particle… ▽ More Magnetic reconnection, especially in the relativistic regime, provides an efficient mechanism for accelerating relativistic particles and thus offers an attractive physical explanation for nonthermal high-energy emission from various astrophysical sources. I present a simple analytical model that elucidates key physical processes responsible for reconnection-driven relativistic nonthermal particle acceleration (NTPA) in the large-system, plasmoid-dominated regime in two dimensions. The model aims to explain the numerically-observed dependencies of the power-law index $p$ and high-energy cutoff $γ_c$ of the resulting nonthermal particle energy spectrum $f(γ)$ on the ambient plasma magnetization $σ$, and (for $γ_c$) on the system size $L$. In this self-similar model, energetic particles are continuously accelerated by the out-of-plane reconnection electric field $E_{\rm rec}$ until they become magnetized by the reconnected magnetic field and eventually trapped in plasmoids large enough to confine them. The model also includes diffusive Fermi acceleration by particle bouncing off rapidly moving plasmoids. I argue that the balance between electric acceleration and magnetization controls the power-law index, while trapping in plasmoids governs the cutoff, thus tying the particle energy spectrum to the plasmoid distribution. △ Less

Submitted 18 July, 2020; originally announced July 2020.

Comments: Submitted for publication; 49 pages, 4 figures

Journal ref: Journal of Plasma Physics, vol. 88, Issue 1, article id.905880114 (2022)

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena throughout the Universe

Authors: H. Ji, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, L. -J. Chen, Y. Chen, A. Chien, D. Craig, J. Dahlin, W. Daughton, E. DeLuca, C. F. Dong, S. Dorfman, J. Drake, F. Ebrahimi , et al. (75 additional authors not shown)

Abstract: This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 6 pages, 1 figure, white paper submitted to both Plasma 2020 and Astro 2020 Decadal Surveys

Fast magnetic reconnection in highly-extended current sheets at the National Ignition Facility

Authors: W. Fox, D. B. Schaeffer, M. J. Rosenberg, G. Fiksel, J. Matteucci, H. -S. Park, A. F. A. Bott, K. Lezhnin, A. Bhattacharjee, D. Kalantar, B. A. Remington, D. Uzdensky, C. K. Li, F. H. Séguin, S. X. Hu

Abstract: Fast magnetic reconnection was observed between magnetized laser-produced plasmas at the National Ignition Facility. Two highly-elongated plasma plumes were produced by tiling two rows of lasers, with magnetic fields generated in each plume by the Biermann battery effect. Detailed magnetic field observations, obtained from proton radiography using a D$^3$He capsule implosion, reveal reconnection o… ▽ More Fast magnetic reconnection was observed between magnetized laser-produced plasmas at the National Ignition Facility. Two highly-elongated plasma plumes were produced by tiling two rows of lasers, with magnetic fields generated in each plume by the Biermann battery effect. Detailed magnetic field observations, obtained from proton radiography using a D$^3$He capsule implosion, reveal reconnection occurring in an extended, quasi-1D current sheet with large aspect ratio $\sim 100$. The 1-D geometry allowed a rigorous and unique reconstruction of the magnetic field, which showed a reconnection current sheet that thinned down to a half-width close to the electron gyro-scale. Despite the large aspect ratio, a large fraction of the magnetic flux reconnected, suggesting fast reconnection supported by the non-gyrotropic electron pressure tensor. △ Less

Submitted 13 March, 2020; originally announced March 2020.

Comments: 5 pages, 5 Figures, Supplementary Info

Kinetic Beaming in Radiative Relativistic Magnetic Reconnection: A Mechanism for Rapid Gamma-Ray Flares in Jets

Authors: J. M. Mehlhaff, G. R. Werner, D. A. Uzdensky, M. C. Begelman

Abstract: Rapid gamma-ray flares pose an astrophysical puzzle, requiring mechanisms both to accelerate energetic particles and to produce fast observed variability. These dual requirements may be satisfied by collisionless relativistic magnetic reconnection. On the one hand, relativistic reconnection can energize gamma-ray emitting electrons. On the other, as previous kinetic simulations have shown, the rec… ▽ More Rapid gamma-ray flares pose an astrophysical puzzle, requiring mechanisms both to accelerate energetic particles and to produce fast observed variability. These dual requirements may be satisfied by collisionless relativistic magnetic reconnection. On the one hand, relativistic reconnection can energize gamma-ray emitting electrons. On the other, as previous kinetic simulations have shown, the reconnection acceleration mechanism preferentially focuses high-energy particles -- and their emitted photons -- into beams, which may create rapid blips in flux as they cross a telescope's line of sight. Using a series of 2D pair-plasma particle-in-cell simulations, we explicitly demonstrate the critical role played by radiative cooling in mediating the observable signatures of this `kinetic beaming' effect. Only in our efficiently cooled simulations do we measure kinetic beaming beyond one light crossing time of the reconnection layer. We find a correlation between the cooling strength and the photon energy range across which persistent kinetic beaming occurs: stronger cooling coincides with a wider range of beamed photon energies. We also apply our results to rapid gamma-ray flares in flat-spectrum radio quasars, suggesting that a paradigm of radiatively efficient kinetic beaming constrains relevant emission models. In particular, beaming-produced variability may be more easily realized in two-zone (e.g. spine-sheath) set-ups, with Compton seed photons originating in the jet itself, rather than in one-zone external Compton scenarios. △ Less

Submitted 25 August, 2020; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: 25 pages, 23 figures, accepted by MNRAS

Multi-scale dynamics of magnetic flux tubes and inverse magnetic energy transfer

Authors: Muni Zhou, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: We report on an analytical and numerical study of the dynamics of a three-dimensional array of identical magnetic flux tubes in the reduced-magnetohydrodynamic description of the plasma. We propose that the long-time evolution of this system is dictated by flux-tube mergers and that such mergers are dynamically constrained by the conservation of the pertinent (ideal) invariants, {\it viz.} the mag… ▽ More We report on an analytical and numerical study of the dynamics of a three-dimensional array of identical magnetic flux tubes in the reduced-magnetohydrodynamic description of the plasma. We propose that the long-time evolution of this system is dictated by flux-tube mergers and that such mergers are dynamically constrained by the conservation of the pertinent (ideal) invariants, {\it viz.} the magnetic potential and axial fluxes of each tube. We also propose that in the direction perpendicular to the merging plane, flux tubes evolve in critically-balanced fashion. These notions allow us to construct an analytical model for how quantities such as the magnetic energy and the energy-containing scale evolve as functions of time. Of particular importance is the conclusion that, like its two-dimensional counterpart, this system exhibits an inverse transfer of magnetic energy that terminates only at the system scale. We perform direct numerical simulations that confirm these predictions and reveal other interesting aspects of the evolution of the system. We find, for example, that the early time evolution is characterized by a sharp decay of the initial magnetic energy, which we attribute to the ubiquitous formation of current sheets. We also show that a quantitatively similar inverse transfer of magnetic energy is observed when the initial condition is a random, small-scale magnetic seed field. △ Less

Submitted 13 June, 2020; v1 submitted 20 January, 2020; originally announced January 2020.

Comments: 33 pages, 19 figures, accepted for publication in Journal of Plasma Physics

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

Magnetic island merger as a mechanism for inverse magnetic energy transfer

Authors: Muni Zhou, Pallavi Bhat, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: Magnetic energy transfer from small to large scales due to successive magnetic island coalescence is investigated. A solvable analytical model is introduced and shown to correctly capture the evolution of the main quantities of interest, as borne out by numerical simulations. Magnetic reconnection is identified as the key mechanism enabling the inverse transfer, and setting its properties: magneti… ▽ More Magnetic energy transfer from small to large scales due to successive magnetic island coalescence is investigated. A solvable analytical model is introduced and shown to correctly capture the evolution of the main quantities of interest, as borne out by numerical simulations. Magnetic reconnection is identified as the key mechanism enabling the inverse transfer, and setting its properties: magnetic energy decays as $\tilde t^{-1}$, where $\tilde t$ is time normalized to the (appropriately defined) reconnection timescale; and the correlation length of the field grows as $\tilde t^{1/2}$. The magnetic energy spectrum is self-similar, and evolves as $\propto \tilde t ^{-3/2}k^{-2}$, where the $k$-dependence is imparted by the formation of thin current sheets. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: 6 pages, 5 figures, submitted for publication

Journal ref: Phys. Rev. Research 1, 012004 (2019)

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)

Analytical Properties of Linear Electrostatic Waves in Two-Component Quantum and Classical Plasmas

Authors: Shane Rightley, Dmitri Uzdensky

Abstract: We examine the properties of linear electrostatic waves in unmagnetized quantum and classical plasmas consisting of one or two populations of electrons with analytically tractable distribution functions in the presence of a stationary neutralising ion background. Beginning with the kinetic quantum plasma longitudinal susceptibility, we assess the effects due to increasing complexity of the backgro… ▽ More We examine the properties of linear electrostatic waves in unmagnetized quantum and classical plasmas consisting of one or two populations of electrons with analytically tractable distribution functions in the presence of a stationary neutralising ion background. Beginning with the kinetic quantum plasma longitudinal susceptibility, we assess the effects due to increasing complexity of the background distribution function. Firstly, we analyse dispersion and Landau damping in one-component plasmas and consider distribution functions with a variety of analytical properties: the Dirac delta function, the Cauchy profile with two complex first-order poles, the squared Cauchy profile with two second order poles, and the inverse-quartic profile with four first-order poles; we also briefly discuss the non-meromorphic totally- and arbitrarily-degenerate Fermi-Dirac distribution. In order to study electrostatic instabilities, we then turn to plasmas with two populations of electrons streaming relative to each other in two cases: a symmetric case of two counter-streaming identical populations and a bump-on-tail case with a primary population and a delta-function beam. We obtain the corresponding linear kinetic dispersion relations and evaluate the properties of instabilities when the electron distribution functions are of the delta function, Cauchy, squared-Cauchy, or inverse-quartic types. In agreement with other studies, we find that in general quantum effects reduce the range of wavelengths for unstable modes at long wavelengths. We also find a second window of instability at shorter wavelengths and elucidate its nature as being due to quantum recoil. We note the possible implications for studies of laboratory and astrophysical quantum plasmas. △ Less

Submitted 4 July, 2018; originally announced July 2018.

Comments: 44 pages, 22 figures, 2 tables

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)

Particle acceleration in relativistic magnetic reconnection with strong inverse-Compton cooling in pair plasmas

Authors: Gregory R. Werner, Alexander A. Philippov, Dmitri A. Uzdensky

Abstract: Particle-in-cell (PIC) simulations have shown that relativistic collisionless magnetic reconnection drives nonthermal particle acceleration (NTPA), potentially explaining high-energy (X-ray/$γ$-ray) synchrotron and/or inverse Compton (IC) radiation observed from various astrophysical sources. The radiation back-reaction force on radiating particles has been neglected in most of these simulations,… ▽ More Particle-in-cell (PIC) simulations have shown that relativistic collisionless magnetic reconnection drives nonthermal particle acceleration (NTPA), potentially explaining high-energy (X-ray/$γ$-ray) synchrotron and/or inverse Compton (IC) radiation observed from various astrophysical sources. The radiation back-reaction force on radiating particles has been neglected in most of these simulations, even though radiative cooling considerably alters particle dynamics in many astrophysical environments where reconnection may be important. We present a radiative PIC study examining the effects of external IC cooling on the basic dynamics, NTPA, and radiative signatures of relativistic reconnection in pair plasmas. We find that, while the reconnection rate and overall dynamics are basically unchanged, IC cooling significantly influences NTPA: the particle spectra still show a hard power law (index $\geq -2$) as in nonradiative reconnection, but transition to a steeper power law that extends to a cooling-dependent cutoff. The steep power-law index fluctuates in time between roughly $-$3 and $-$5. The time-integrated photon spectra display corresponding power laws with indices $\approx -0.5$ and $\approx -1.1$, similar to those observed in hard X-ray spectra of accreting black holes. △ Less

Submitted 12 November, 2018; v1 submitted 4 May, 2018; originally announced May 2018.

Comments: 5 pages, 5 figures

Journal ref: published in MNRAS Lett. 482(1), L60-L64 (2019)

Development of Tearing Instability in a Current Sheet Forming by Sheared Incompressible Flow

Authors: Elizabeth A. Tolman, Nuno F. Loureiro, Dmitri A. Uzdensky

Abstract: Sweet-Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvénic incompressible flows which narrow the sheet. Standard tearing theory (Furth et al. 1963; Coppi et al. 1976; R… ▽ More Sweet-Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvénic incompressible flows which narrow the sheet. Standard tearing theory (Furth et al. 1963; Coppi et al. 1976; Rutherford 1973) is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal MHD, a coordinate transformation reveals that the flow compresses and stretches perturbations. A multiple scale formulation describes how linear tearing mode theory (Furth et al. 1963; Coppi et al. 1976) can be applied to an equilibrium changing under flow, showing that the flow affects the separable exponential growth only implicitly, by making the standard scalings time-dependent. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet. △ Less

Submitted 28 February, 2018; v1 submitted 24 November, 2017; originally announced November 2017.

Journal ref: Journal of Plasma Physics 84.1 (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)

Nonthermal Particle Acceleration in 3D Relativistic Magnetic Reconnection in Pair Plasma

Authors: Gregory R. Werner, Dmitri A. Uzdensky

Abstract: As a fundamental process converting magnetic to plasma energy in high-energy astrophysical plasmas, relativistic magnetic reconnection is a leading explanation for the acceleration of particles to the ultrarelativistic energies necessary to power nonthermal emission (especially X-rays and gamma-rays) in pulsar magnetospheres and pulsar wind nebulae, coronae and jets of accreting black holes, and g… ▽ More As a fundamental process converting magnetic to plasma energy in high-energy astrophysical plasmas, relativistic magnetic reconnection is a leading explanation for the acceleration of particles to the ultrarelativistic energies necessary to power nonthermal emission (especially X-rays and gamma-rays) in pulsar magnetospheres and pulsar wind nebulae, coronae and jets of accreting black holes, and gamma-ray bursts. An important objective of plasma astrophysics is therefore the characterization of nonthermal particle acceleration (NTPA) effected by reconnection. Reconnection-powered NTPA has been demonstrated over a wide range of physical conditions using large two-dimensional (2D) kinetic simulations. However, its robustness in realistic 3D reconnection -- in particular, whether the 3D relativistic drift-kink instability (RDKI) disrupts NTPA -- has not been systematically investigated, although pioneering 3D simulations have observed NTPA in isolated cases. Here we present the first comprehensive study of NTPA in 3D relativistic reconnection in collisionless electron-positron plasmas, characterizing NTPA as the strength of 3D effects is varied systematically via the length in the third dimension and the strength of the guide magnetic field. We find that, while the RDKI prominently perturbs 3D reconnecting current sheets, it does not suppress particle acceleration, even for zero guide field; fully 3D reconnection robustly and efficiently produces nonthermal power-law particle spectra closely resembling those obtained in 2D. This finding provides strong support for reconnection as the key mechanism powering high-energy flares in various astrophysical systems. We also show that strong guide fields significantly inhibit NTPA, slowing reconnection and limiting the energy available for plasma energization, yielding steeper and shorter power-law spectra. △ Less

Submitted 11 July, 2017; v1 submitted 15 May, 2017; originally announced May 2017.

Comments: 10 pages, 4 figures, published in ApJL 843:L27 (2017)

Journal ref: ApJL 843:L27 (2017)

Relativistic Turbulence with Strong Synchrotron and Synchrotron-Self-Compton Cooling

Authors: Dmitri A. Uzdensky

Abstract: Many relativistic plasma environments in high-energy astrophysics, including pulsar wind nebulae, hot accretion flows onto black holes, relativistic jets in active galactic nuclei and gamma-ray bursts, and giant radio lobes, are naturally turbulent. The plasma in these environments is often so hot that synchrotron and inverse-Compton (IC) radiative cooling becomes important. In this paper we inves… ▽ More Many relativistic plasma environments in high-energy astrophysics, including pulsar wind nebulae, hot accretion flows onto black holes, relativistic jets in active galactic nuclei and gamma-ray bursts, and giant radio lobes, are naturally turbulent. The plasma in these environments is often so hot that synchrotron and inverse-Compton (IC) radiative cooling becomes important. In this paper we investigate the general thermodynamic and radiative properties (and hence the observational appearance) of an optically thin relativistically hot plasma stirred by driven magnetohydrodynamic (MHD) turbulence and cooled by radiation. We find that if the system reaches a statistical equilibrium where turbulent heating is balanced by radiative cooling, the effective electron temperature tends to attain a universal value $θ= kT_e/m_e c^2 \sim 1/\sqrt{τ_T}$, where $τ_T=n_eσ_T L \ll 1$ is the system's Thomson optical depth, essentially independent of the strength of turbulent driving or magnetic field. This is because both MHD turbulent dissipation and synchrotron cooling are proportional to the magnetic energy density. We also find that synchrotron self-Compton (SSC) cooling and perhaps a few higher-order IC components are automatically comparable to synchrotron in this regime. The overall broadband radiation spectrum then consists of several distinct components (synchrotron, SSC, etc.), well separated in photon energy (by a factor $\sim τ_T^{-1}$) and roughly equal in power. The number of IC peaks is checked by Klein-Nishina effects and depends logarithmically on $τ_T$ and the magnetic field. We also examine the limitations due to synchrotron self-absorption, explore applications to Crab PWN and blazar jets, and discuss links to radiative magnetic reconnection. △ Less

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

Comments: 12 pages, 1 figure; replaced with a slightly modified revised version matching the final published (in MNRAS) article

Journal ref: MNRAS, vol. 477, issue 3, pp. 2849-2857 (2018)

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)

Magnetic reconnection: from the Sweet-Parker model to stochastic plasmoid chains

Authors: N. F. Loureiro, D. A. Uzdensky

Abstract: (abridged) Magnetic reconnection is the topological reconfiguration of the magnetic field in a plasma, accompanied by the violent release of energy and particle acceleration. Reconnection is as ubiquitous as plasmas themselves, with solar flares perhaps the most popular example. Over the last few years, the theoretical understanding of magnetic reconnection in large-scale fluid systems has undergo… ▽ More (abridged) Magnetic reconnection is the topological reconfiguration of the magnetic field in a plasma, accompanied by the violent release of energy and particle acceleration. Reconnection is as ubiquitous as plasmas themselves, with solar flares perhaps the most popular example. Over the last few years, the theoretical understanding of magnetic reconnection in large-scale fluid systems has undergone a major paradigm shift. The steady-state model of reconnection described by the famous Sweet-Parker (SP) theory, which dominated the field for ~50 years, has been replaced with an essentially time-dependent, bursty picture of the reconnection layer, dominated by the continuous formation and ejection of multiple secondary islands (plasmoids). Whereas in the SP model reconnection was predicted to be slow, a major implication of this new paradigm is that reconnection in fluid systems is fast (i.e., independent of the Lundquist number), provided that the system is large enough. This conceptual shift hinges on the realization that SP-like current layers are violently unstable to the plasmoid instability - implying, therefore, that such current sheets are super-critically unstable and thus can never form in the first place. This suggests that the formation of a current sheet and the subsequent reconnection process cannot be decoupled, as is commonly assumed. This paper provides an introductory-level overview of the recent developments in reconnection theory and simulations that led to this essentially new framework. We briefly discuss the role played by the plasmoid instability in selected applications, and describe some of the outstanding challenges that remain at the frontier of this subject. Amongst these are the analytical and numerical extension of the plasmoid instability to (i) 3D and (ii) non-MHD regimes. New results are reported in both cases. △ Less

Submitted 28 July, 2015; originally announced July 2015.

Comments: Submitted to Plasma Phys. Control. Fusion, EPS 2015 Special Issue

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)

Landau Damping of Electrostatic Waves in Arbitrarily Degenerate Quantum Plasmas

Authors: Shane Rightley, Dmitri Uzdensky

Abstract: We carry out a systematic study of the dispersion relation for linear electrostatic waves in an arbitrarily degenerate quantum electron plasma. We solve for the complex frequency spectrum for arbitrary values of wavenumber $k$ and level of degeneracy $μ$. Our finding is that for large $k$ and high $μ$ the real part of the frequency $ω_{r}$ grows linearly with $k$ and scales with $μ$ only because o… ▽ More We carry out a systematic study of the dispersion relation for linear electrostatic waves in an arbitrarily degenerate quantum electron plasma. We solve for the complex frequency spectrum for arbitrary values of wavenumber $k$ and level of degeneracy $μ$. Our finding is that for large $k$ and high $μ$ the real part of the frequency $ω_{r}$ grows linearly with $k$ and scales with $μ$ only because of the scaling of the Fermi energy. In this regime the relative Landau damping rate $γ/ω_{r}$ becomes independent of $k$ and varies inversly with $μ$. Thus, damping is weak but finite at moderate levels of degeneracy for short wavelengths. △ Less

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

Comments: 5 pages, 2 figures, replaced with the final version accepted for publication in Physics of Plasmas Letters

Large-Eddy Simulations of Magnetohydrodynamic Turbulence in Heliophysics and Astrophysics

Authors: Mark S. Miesch, William H. Matthaeus, Axel Brandenburg, Arakel Petrosyan, Annick Pouquet, Claude Cambon, Frank Jenko, Dmitri Uzdensky, James Stone, Steve Tobias, Juri Toomre, Marco Velli

Abstract: We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as i… ▽ More We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions. △ Less

Submitted 18 July, 2015; v1 submitted 7 May, 2015; originally announced May 2015.

Comments: 51 pages, 6 figures (Figs 2, 3, and 4 color), accepted to Space Science Reviews (in press). The paper is a product of a workshop on "LES Modeling in MHD Turbulence" held in Boulder, CO in May, 2013, sponsored by the Geophysical Turbulence Program at the National Center for Atmospheric Research

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)

Magnetic Reconnection Onset via Disruption of a Forming Current Sheet by the Tearing Instability

Authors: Dmitri A. Uzdensky, Nuno F. Loureiro

Abstract: The recent realization that Sweet-Parker current sheets are violently unstable to the secondary tearing (plasmoid) instability implies that such current sheets cannot occur in real systems. This suggests that, in order to understand the onset of magnetic reconnection, one needs to consider the growth of the tearing instability in a current layer as it is being formed. Such an analysis is performed… ▽ More The recent realization that Sweet-Parker current sheets are violently unstable to the secondary tearing (plasmoid) instability implies that such current sheets cannot occur in real systems. This suggests that, in order to understand the onset of magnetic reconnection, one needs to consider the growth of the tearing instability in a current layer as it is being formed. Such an analysis is performed here in the context of nonlinear resistive MHD for a generic time-dependent equilibrium representing a gradually forming current sheet. It is shown that two onset regimes, single-island and multi-island, are possible, depending on the rate of current sheet formation. A simple model is used to compute the criterion for transition between these two regimes, as well as the reconnection onset time and the current sheet parameters at that moment. For typical solar corona parameters this model yields results consistent with observations. △ Less

Submitted 19 February, 2016; v1 submitted 16 November, 2014; originally announced November 2014.

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

Journal ref: Phys. Rev. Lett. 116, 105003 (2016)

Plasma Physics of Extreme Astrophysical Environments

Authors: Dmitri A. Uzdensky, Shane Rightley

Abstract: Certain classes of astrophysical objects, namely magnetars and central engines of supernovae and gamma-ray bursts (GRBs), are characterized by extreme physical conditions not encountered elsewhere in the Universe. In particular, they possess magnetic fields that exceed the critical quantum field of 44 teragauss. Figuring out how these complex ultra-magnetized systems work requires understanding va… ▽ More Certain classes of astrophysical objects, namely magnetars and central engines of supernovae and gamma-ray bursts (GRBs), are characterized by extreme physical conditions not encountered elsewhere in the Universe. In particular, they possess magnetic fields that exceed the critical quantum field of 44 teragauss. Figuring out how these complex ultra-magnetized systems work requires understanding various plasma processes, both small-scale kinetic and large-scale magnetohydrodynamic (MHD). However, an ultra-strong magnetic field modifies the underlying physics to such an extent that many relevant plasma-physical problems call for building QED-based relativistic quantum plasma physics. In this review, after describing the extreme astrophysical systems of interest and identifying the key relevant plasma-physical problems, we survey the recent progress in the development of such a theory. We discuss how a super-critical field modifies the properties of vacuum and matter and outline the basic theoretical framework for describing both non-relativistic and relativistic quantum plasmas. We then turn to astrophysical applications of relativistic QED plasma physics relevant to magnetar magnetospheres and central engines of supernovae and long GRBs. Specifically, we discuss propagation of light through a magnetar magnetosphere; large-scale MHD processes driving magnetar activity and GRB jet launching and propagation; energy-transport processes governing the thermodynamics of extreme plasma environments; micro-scale kinetic plasma processes important in the interaction of intense magnetospheric electric currents with a magnetar's surface; and magnetic reconnection of ultra-strong magnetic fields. Finally, we point out that future progress will require the development of numerical modeling capabilities. △ Less

Submitted 20 January, 2014; originally announced January 2014.

Comments: 56 pages, 1 figure. Invited review accepted for publication in Reports on Progress in Physics

Gamma-ray flares in the Crab Nebula: A case of relativistic reconnection?

Authors: Benoit Cerutti, Gregory R. Werner, Dmitri A. Uzdensky, Mitchell C. Begelman

Abstract: The Crab Nebula was formed after the collapse of a massive star about a thousand years ago, leaving behind a pulsar that inflates a bubble of ultra-relativistic electron-positron pairs permeated with magnetic field. The observation of brief but bright flares of energetic gamma rays suggests that pairs are accelerated to PeV energies within a few days; such rapid acceleration cannot be driven by sh… ▽ More The Crab Nebula was formed after the collapse of a massive star about a thousand years ago, leaving behind a pulsar that inflates a bubble of ultra-relativistic electron-positron pairs permeated with magnetic field. The observation of brief but bright flares of energetic gamma rays suggests that pairs are accelerated to PeV energies within a few days; such rapid acceleration cannot be driven by shocks. Here, it is argued that the flares may be the smoking gun of magnetic dissipation in the Nebula. Using 2D and 3D particle-in-cell simulations, it is shown that the observations are consistent with relativistic magnetic reconnection, where pairs are subject to strong radiative cooling. The Crab flares may highlight the importance of relativistic magnetic reconnection in astrophysical sources. △ Less

Submitted 13 January, 2014; originally announced January 2014.

Comments: 8 pages, 5 figures, Invited paper accepted for publication in the annual Special Issue of Physics of Plasmas, based on an invited talk at the 55th APS DPP meeting, November 11-15, 2013, held in Denver, CO

Three-dimensional relativistic pair plasma reconnection with radiative feedback in the Crab Nebula

Authors: Benoit Cerutti, Gregory R. Werner, Dmitri A. Uzdensky, Mitchell C. Begelman

Abstract: The discovery of rapid synchrotron gamma-ray flares above 100 MeV from the Crab Nebula has attracted new interest in alternative particle acceleration mechanisms in pulsar wind nebulae. Diffuse shock-acceleration fails to explain the flares because particle acceleration and emission occur during a single or even sub-Larmor timescale. In this regime, the synchrotron energy losses induce a drag forc… ▽ More The discovery of rapid synchrotron gamma-ray flares above 100 MeV from the Crab Nebula has attracted new interest in alternative particle acceleration mechanisms in pulsar wind nebulae. Diffuse shock-acceleration fails to explain the flares because particle acceleration and emission occur during a single or even sub-Larmor timescale. In this regime, the synchrotron energy losses induce a drag force on the particle motion that balances the electric acceleration and prevents the emission of synchrotron radiation above 160 MeV. Previous analytical studies and 2D particle-in-cell (PIC) simulations indicate that relativistic reconnection is a viable mechanism to circumvent the above difficulties. The reconnection electric field localized at X-points linearly accelerates particles with little radiative energy losses. In this paper, we check whether this mechanism survives in 3D, using a set of large PIC simulations with radiation reaction force and with a guide field. In agreement with earlier works, we find that the relativistic drift kink instability deforms and then disrupts the layer, resulting in significant plasma heating but few non-thermal particles. A moderate guide field stabilizes the layer and enables particle acceleration. We report that 3D magnetic reconnection can accelerate particles above the standard radiation reaction limit, although the effect is less pronounced than in 2D with no guide field. We confirm that the highest energy particles form compact bunches within magnetic flux ropes, and a beam tightly confined within the reconnection layer, which could result in the observed Crab flares when, by chance, the beam crosses our line of sight. △ Less

Submitted 13 January, 2014; v1 submitted 11 November, 2013; originally announced November 2013.

Comments: 16 pages, 16 figures, Accepted for publication in The Astrophysical Journal

Simulations of particle acceleration beyond the classical synchrotron burnoff limit in magnetic reconnection: An explanation of the Crab flares

Authors: Benoit Cerutti, Gregory R. Werner, Dmitri A. Uzdensky, Mitchell C. Begelman

Abstract: It is generally accepted that astrophysical sources cannot emit synchrotron radiation above 160 MeV in their rest frame. This limit is given by the balance between the accelerating electric force and the radiation reaction force acting on the electrons. The discovery of synchrotron gamma-ray flares in the Crab Nebula, well above this limit, challenges this classical picture of particle acceleratio… ▽ More It is generally accepted that astrophysical sources cannot emit synchrotron radiation above 160 MeV in their rest frame. This limit is given by the balance between the accelerating electric force and the radiation reaction force acting on the electrons. The discovery of synchrotron gamma-ray flares in the Crab Nebula, well above this limit, challenges this classical picture of particle acceleration. To overcome this limit, particles must accelerate in a region of high electric field and low magnetic field. This is possible only with a non-ideal magnetohydrodynamic process, like magnetic reconnection. We present the first numerical evidence of particle acceleration beyond the synchrotron burnoff limit, using a set of 2D particle-in-cell simulations of ultra-relativistic pair plasma reconnection. We use a new code, Zeltron, that includes self-consistently the radiation reaction force in the equation of motion of the particles. We demonstrate that the most energetic particles move back and forth across the reconnection layer, following relativistic Speiser orbits. These particles then radiate >160 MeV synchrotron radiation rapidly, within a fraction of a full gyration, after they exit the layer. Our analysis shows that the high-energy synchrotron flux is highly variable in time because of the strong anisotropy and inhomogeneity of the energetic particles. We discover a robust positive correlation between the flux and the cut-off energy of the emitted radiation, mimicking the effect of relativistic Doppler amplification. A strong guide field quenches the emission of >160 MeV synchrotron radiation. Our results are consistent with the observed properties of the Crab flares, supporting the reconnection scenario. △ Less

Submitted 7 May, 2013; v1 submitted 25 February, 2013; originally announced February 2013.

Comments: 15 pages, 16 figures, Accepted for publication in The Astrophysical Journal