-
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.
-
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.
-
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.
-
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.
-
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.
-
Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets
Authors:
José Ortuño-Macías,
Krzysztof Nalewajko,
Dmitri A. Uzdensky,
Mitchell C. Begelman,
Gregory R. Werner,
Alexander Y. Chen,
Bhupendra Mishra
Abstract:
Relativistic magnetized jets, such as those from AGN, GRBs and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and non-thermal radiation. Here we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron-positron pair plasma.…
▽ More
Relativistic magnetized jets, such as those from AGN, GRBs and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and non-thermal radiation. Here we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron-positron pair plasma. Generalizing recent treatments by Alves et al. (2018) and Davelaar et al. (2020), we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. (2018) is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.
△ Less
Submitted 6 May, 2022; v1 submitted 4 May, 2022;
originally announced May 2022.
-
Spontaneous magnetization of collisionless plasma through the action of a shear flow
Authors:
Muni Zhou,
Vladimir Zhdankin,
Matthew W. Kunz,
Nuno F. Loureiro,
Dmitri A. Uzdensky
Abstract:
We study in a fully kinetic framework the generation of seed magnetic fields through the Weibel instability driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and its saturation when the se…
▽ More
We study in a fully kinetic framework the generation of seed magnetic fields through the Weibel instability driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and its saturation when the seed magnetic fields become strong enough to instigate gyromotion of particles and thereby inhibit their free-streaming. The predicted scaling dependencies of the saturated seed fields on key parameters (e.g., ratio of system scale to electron skin depth, the forcing amplitude) are confirmed by 3D and 2D particle-in-cell simulations using an electron-positron plasma. This work demonstrates the spontaneous magnetization of a collisionless plasma through large-scale motions as simple as a shear flow, and therefore has important implications for magnetogenesis in dilute astrophysical systems.
△ Less
Submitted 3 October, 2021;
originally announced October 2021.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Pulsar Radio Emission Mechanism: Radio Nanoshots as a Low Frequency Afterglow of Relativistic Magnetic Reconnection
Authors:
Alexander Philippov,
Dmitri A. Uzdensky,
Anatoly Spitkovsky,
Benoît Cerutti
Abstract:
In this Letter we propose that coherent radio emission of Crab, other young energetic pulsars, and millisecond pulsars is produced in the magnetospheric current sheet beyond the light cylinder. We carry out global and local two-dimensional kinetic plasma simulations of reconnection to illustrate the coherent emission mechanism. Reconnection in the current sheet beyond the light cylinder proceeds i…
▽ More
In this Letter we propose that coherent radio emission of Crab, other young energetic pulsars, and millisecond pulsars is produced in the magnetospheric current sheet beyond the light cylinder. We carry out global and local two-dimensional kinetic plasma simulations of reconnection to illustrate the coherent emission mechanism. Reconnection in the current sheet beyond the light cylinder proceeds in the very efficient plasmoid-dominated regime, and current layer gets fragmented into a dynamic chain of plasmoids which undergo successive coalescence. Mergers of sufficiently large plasmoids produce secondary perpendicular current sheets, which are also plasmoid-unstable. Collisions of plasmoids with each other and with the upstream magnetic field eject fast-magnetosonic waves, which propagate upstream across the background field and successfully escape from the plasma as electromagnetic waves that fall in the radio band. This model successfully explains many important features of the observed radio emission from Crab and other pulsars with high magnetic field at the light cylinder: phase coincidence with the high-energy emission, nano-second duration (nanoshots), and extreme instantaneous brightness of individual pulses.
△ Less
Submitted 20 February, 2019;
originally announced February 2019.
-
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.
-
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.
-
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.
-
GRB and blazar jets shining through their stripes
Authors:
Dimitrios Giannios,
Dmitri A. Uzdensky
Abstract:
Black-hole driven relativistic astrophysical jets, such as blazars and gamma-ray bursts (GRBs), are powerful sources of electromagnetic radiation. Their emission is powered by some energy dissipation and particle acceleration mechanism operating in-situ at large distances from the black hole. We propose that the formation of the dissipative structures in the jet is controlled by the time variabili…
▽ More
Black-hole driven relativistic astrophysical jets, such as blazars and gamma-ray bursts (GRBs), are powerful sources of electromagnetic radiation. Their emission is powered by some energy dissipation and particle acceleration mechanism operating in-situ at large distances from the black hole. We propose that the formation of the dissipative structures in the jet is controlled by the time variability of the accretion disc. We argue that the open magnetic field lines through the black hole, which drive a strongly magnetized jet, may have their polarity reversing over time-scales related to the growth of the magneto-rotational dynamo in the disc. Consequently, the jet is characterized by an alternating toroidal field polarity along its propagation axis, i.e., it is a "striped jet". Magnetic reconnection in the current sheets that form between the stripes dissipates the alternating-field energy and thus powers further jet acceleration. Here, we consider a jet with a broad distribution of stripe widths $l>l_{\rm min}$, above a dominant scale $l_{\rm min}$. We find that the bulk acceleration of the jet, driven by the annihilation of the stripes, is very gradual. The dissipation rate peaks at a distance $z_{\rm peak}\sim 10^6 R_g\, (Γ_\infty/30)^2\, (l_{\rm min}/1000R_g)$, where $R_{g}$ is the black-hole's gravititional radius and $Γ_{\infty}$ the jet's asymptotic Lorentz factor, and exhibits a very broad plateau extending by $\sim 4-5$ orders of magnitude in distance. The prolonged energy dissipation accounts for the flat-to-inverted long-wavelength spectra commonly observed in jets. The model can also account for the broad range of flaring timescales of blazars and the fact that their bulk acceleration appears to continue out to ~100 pc scales. In GRB jets, the model predicts comparable power for the photopsheric and Thomson-thin emission components.
△ Less
Submitted 23 May, 2018;
originally announced May 2018.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Non-thermal particle acceleration in collisionless relativistic electron-proton reconnection
Authors:
G. R. Werner,
D. A. Uzdensky,
M. C. Begelman,
B. Cerutti,
K. Nalewajko
Abstract:
Magnetic reconnection in relativistic collisionless plasmas can accelerate particles and power high-energy emission in various astrophysical systems. Whereas most previous studies focused on relativistic reconnection in pair plasmas, less attention has been paid to electron-ion plasma reconnection, expected in black hole accretion flows and relativistic jets. We report a comprehensive particle-in-…
▽ More
Magnetic reconnection in relativistic collisionless plasmas can accelerate particles and power high-energy emission in various astrophysical systems. Whereas most previous studies focused on relativistic reconnection in pair plasmas, less attention has been paid to electron-ion plasma reconnection, expected in black hole accretion flows and relativistic jets. We report a comprehensive particle-in-cell numerical investigation of reconnection in an electron-ion plasma, spanning a wide range of ambient ion magnetizations $σ_i$, from the semirelativistic regime (ultrarelativistic electrons but nonrelativistic ions, 0.001<<$σ_i$<<1) to the fully relativistic regime (both species are ultrarelativistic, $σ_i$>>1). We investigate how the reconnection rate, electron and ion plasma flows, electric and magnetic field structures, electron/ion energy partitioning, and nonthermal particle acceleration depend on $σ_i$. Our key findings are: (1) the reconnection rate is about 0.1 of the Alfvenic rate across all regimes; (2) electrons can form concentrated moderately relativistic outflows even in the semirelativistic, small-$σ_i$ regime; (3) while the released magnetic energy is partitioned equally between electrons and ions in the ultrarelativistic limit, the electron energy fraction declines gradually with decreased $σ_i$ and asymptotes to about 0.25 in the semirelativistic regime; (4) reconnection leads to efficient nonthermal electron acceleration with a $σ_i$-dependent power-law index, $p(σ_i) \simeq $const$+0.7 {σ_i}^{-1/2}$. These findings are important for understanding black hole systems and lend support to semirelativistic reconnection models for powering nonthermal emission in blazar jets, offering a natural explanation for the spectral indices observed in these systems.
△ Less
Submitted 22 December, 2017; v1 submitted 14 December, 2016;
originally announced December 2016.
-
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.
-
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.
-
Radiative Magnetic Reconnection in Astrophysics
Authors:
Dmitri A. Uzdensky
Abstract:
I review a new rapidly growing area of high-energy plasma astrophysics --- radiative magnetic reconnection, i.e., a reconnection regime where radiation reaction influences reconnection dynamics, energetics, and nonthermal particle acceleration. This influence be may be manifested via a number of astrophysically important radiative effects, such as radiation-reaction limits on particle acceleration…
▽ More
I review a new rapidly growing area of high-energy plasma astrophysics --- radiative magnetic reconnection, i.e., a reconnection regime where radiation reaction influences reconnection dynamics, energetics, and nonthermal particle acceleration. This influence be may be manifested via a number of astrophysically important radiative effects, such as radiation-reaction limits on particle acceleration, radiative cooling, radiative resistivity, braking of reconnection outflows by radiation drag, radiation pressure, viscosity, and even pair creation at highest energy densities. Self-consistent inclusion of these effects in magnetic reconnection theory and modeling calls for serious modifications to our overall theoretical approach to the problem. In addition, prompt reconnection-powered radiation often represents our only observational diagnostic tool for studying remote astrophysical systems; this underscores the importance of developing predictive modeling capabilities to connect the underlying physical conditions in a reconnecting system to observable radiative signatures. This Chapter gives an overview of recent theoretical progress in developing basic physical understanding of radiative reconnection, with a special emphasis on astrophysically important radiation mechanisms like synchrotron, curvature, and inverse-Compton. It also offers a broad review of key high-energy astrophysical applications of radiative reconnection, such as: pulsar wind nebulae and magnetospheres, accreting black-hole coronae and jets in XRBs and AGN, magnetospheres of magnetars, and Gamma-Ray Bursts. Finally, this Chapter discusses the most critical open questions and outlines the directions for future research in this exciting new frontier of plasma astrophysics.
△ Less
Submitted 19 October, 2015;
originally announced October 2015.
-
On the distribution of particle acceleration sites in plasmoid-dominated relativistic magnetic reconnection
Authors:
Krzysztof Nalewajko,
Dmitri A. Uzdensky,
Benoît Cerutti,
Gregory R. Werner,
Mitchell C. Begelman
Abstract:
We investigate the distribution of particle acceleration sites, independently of the actual acceleration mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no in…
▽ More
We investigate the distribution of particle acceleration sites, independently of the actual acceleration mechanism, during plasmoid-dominated, relativistic collisionless magnetic reconnection by analyzing the results of a particle-in-cell numerical simulation. The simulation is initiated with Harris-type current layers in pair plasma with no guide magnetic field, negligible radiative losses, no initial perturbation, and using periodic boundary conditions. We find that the plasmoids develop a robust internal structure, with colder dense cores and hotter outer shells, that is recovered after each plasmoid merger on a dynamical time scale. We use spacetime diagrams of the reconnection layers to probe the evolution of plasmoids, and in this context we investigate the individual particle histories for a representative sample of energetic electrons. We distinguish three classes of particle acceleration sites associated with (1) magnetic X-points, (2) regions between merging plasmoids, and (3) the trailing edges of accelerating plasmoids. We evaluate the contribution of each class of acceleration sites to the final energy distribution of energetic electrons -- magnetic X-points dominate at moderate energies, and the regions between merging plasmoids dominate at higher energies. We also identify the dominant acceleration scenarios, in order of decreasing importance -- (1) single acceleration between merging plasmoids, (2) single acceleration at a magnetic X-point, and (3) acceleration at a magnetic X-point followed by acceleration in a plasmoid. Particle acceleration is absent only in the vicinity of stationary plasmoids. The effect of magnetic mirrors due to plasmoid contraction does not appear to be significant in relativistic reconnection.
△ Less
Submitted 11 November, 2015; v1 submitted 10 August, 2015;
originally announced August 2015.
-
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.
-
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.
-
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.
-
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.
-
The extent of power-law energy spectra in collisionless relativistic magnetic reconnection in pair plasmas
Authors:
G. R. Werner,
D. A. Uzdensky,
B. Cerutti,
K. Nalewajko,
M. C. Begelman
Abstract:
Using two-dimensional particle-in-cell simulations, we characterize the energy spectra of particles accelerated by relativistic magnetic reconnection (without guide field) in collisionless electron-positron plasmas, for a wide range of upstream magnetizations $σ$ and system sizes $L$. The particle spectra are well-represented by a power law $γ^{-α}$, with a combination of exponential and super-exp…
▽ More
Using two-dimensional particle-in-cell simulations, we characterize the energy spectra of particles accelerated by relativistic magnetic reconnection (without guide field) in collisionless electron-positron plasmas, for a wide range of upstream magnetizations $σ$ and system sizes $L$. The particle spectra are well-represented by a power law $γ^{-α}$, with a combination of exponential and super-exponential high-energy cutoffs, proportional to $σ$ and $L$, respectively. For large $L$ and $σ$, the power-law index $α$ approaches about 1.2.
△ Less
Submitted 21 December, 2015; v1 submitted 29 September, 2014;
originally announced September 2014.
-
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.
-
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.
-
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.
-
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.
-
Statistical Analysis of Current Sheets in Three-Dimensional Magnetohydrodynamic Turbulence
Authors:
Vladimir Zhdankin,
Dmitri A. Uzdensky,
Jean C. Perez,
Stanislav Boldyrev
Abstract:
We develop a framework for studying the statistical properties of current sheets in numerical simulations of 3D magnetohydrodynamic (MHD) turbulence. We describe an algorithm that identifies current sheets in a simulation snapshot and then determines their geometrical properties (including length, width, and thickness) and intensities (peak current density and total energy dissipation rate). We th…
▽ More
We develop a framework for studying the statistical properties of current sheets in numerical simulations of 3D magnetohydrodynamic (MHD) turbulence. We describe an algorithm that identifies current sheets in a simulation snapshot and then determines their geometrical properties (including length, width, and thickness) and intensities (peak current density and total energy dissipation rate). We then apply this procedure to simulations of reduced MHD turbulence and perform a statistical analysis on the obtained population of current sheets. We evaluate the role of reconnection by separately studying the populations of current sheets which contain magnetic X-points and those which do not. We find that the statistical properties of the two populations are different in general. We compare the scaling of these properties to phenomenological predictions obtained for the inertial range of MHD turbulence. Finally, we test whether the reconnecting current sheets are consistent with the Sweet-Parker model.
△ Less
Submitted 6 February, 2013;
originally announced February 2013.
-
Physical Conditions in the Reconnection Layer in Pulsar Magnetospheres
Authors:
Dmitri A. Uzdensky,
Anatoly Spitkovsky
Abstract:
The magnetosphere of a rotating pulsar naturally develops a current sheet beyond the light cylinder (LC). Magnetic reconnection in this current sheet inevitably dissipates a nontrivial fraction of the pulsar spin-down power within a few LC radii. We develop a basic physical picture of reconnection in this environment and discuss its implications for the observed pulsed gamma-ray emission. We argue…
▽ More
The magnetosphere of a rotating pulsar naturally develops a current sheet beyond the light cylinder (LC). Magnetic reconnection in this current sheet inevitably dissipates a nontrivial fraction of the pulsar spin-down power within a few LC radii. We develop a basic physical picture of reconnection in this environment and discuss its implications for the observed pulsed gamma-ray emission. We argue that reconnection proceeds in the plasmoid-dominated regime, via an hierarchical chain of multiple secondary islands/flux ropes. The inter-plasmoid reconnection layers are subject to strong synchrotron cooling, leading to significant plasma compression. Using the conditions of pressure balance across these current layers, the balance between the heating by magnetic energy dissipation and synchrotron cooling, and Ampere's law, we obtain simple estimates for key parameters of the layers --- temperature, density, and layer thickness. In the comoving frame of the relativistic pulsar wind just outside of the equatorial current sheet, these basic parameters are uniquely determined by the strength of the reconnecting upstream magnetic field. For the case of the Crab pulsar, we find them to be of order 10 GeV, $10^{13} cm^{-3}$, and 10 cm, respectively. After accounting for the bulk Doppler boosting due to the pulsar wind, the synchrotron and inverse-Compton emission from the reconnecting current sheet can explain the observed pulsed high-energy (GeV) and VHE (~100 GeV) radiation, respectively. Also, we suggest that the rapid relative motions of the secondary plasmoids in the hierarchical chain may contribute to the production of the pulsar radio emission.
△ Less
Submitted 11 October, 2012;
originally announced October 2012.
-
Plasmoid and Kelvin-Helmholtz instabilities in Sweet-Parker current sheets
Authors:
N. F. Loureiro,
A. A. Schekochihin,
D. A. Uzdensky
Abstract:
A 2D linear theory of the instability of Sweet-Parker (SP) current sheets is developed in the framework of Reduced MHD. A local analysis is performed taking into account the dependence of a generic equilibrium profile on the outflow coordinate. The plasmoid instability [Loureiro et al, Phys. Plasmas {\bf 14}, 100703 (2007)] is recovered, i.e., current sheets are unstable to the formation of a larg…
▽ More
A 2D linear theory of the instability of Sweet-Parker (SP) current sheets is developed in the framework of Reduced MHD. A local analysis is performed taking into account the dependence of a generic equilibrium profile on the outflow coordinate. The plasmoid instability [Loureiro et al, Phys. Plasmas {\bf 14}, 100703 (2007)] is recovered, i.e., current sheets are unstable to the formation of a large-wave-number chain of plasmoids ($k_{\rm max}\Lsheet \sim S^{3/8}$, where $k_{\rm max}$ is the wave-number of fastest growing mode, $S=\Lsheet V_A/η$ is the Lundquist number, $\Lsheet$ is the length of the sheet, $V_A$ is the Alfvén speed and $η$ is the plasma resistivity), which grows super-Alfvénically fast ($\gmaxτ_A\sim S^{1/4}$, where $\gmax$ is the maximum growth rate, and $τ_A=\Lsheet/V_A$). For typical background profiles, the growth rate and the wave-number are found to {\it increase} in the outflow direction. This is due to the presence of another mode, the Kelvin-Helmholtz (KH) instability, which is triggered at the periphery of the layer, where the outflow velocity exceeds the Alfvén speed associated with the upstream magnetic field. The KH instability grows even faster than the plasmoid instability, $\gmax τ_A \sim k_{\rm max} \Lsheet\sim S^{1/2}$. The effect of viscosity ($ν$) on the plasmoid instability is also addressed. In the limit of large magnetic Prandtl numbers, $Pm=ν/η$, it is found that $\gmax\sim S^{1/4}Pm^{-5/8}$ and $k_{\rm max} \Lsheet\sim S^{3/8}Pm^{-3/16}$, leading to the prediction that the critical Lundquist number for plasmoid instability in the $Pm\gg1$ regime is $\Scrit\sim 10^4Pm^{1/2}$. These results are verified via direct numerical simulation of the linearized equations, using a new, analytical 2D SP equilibrium solution.
△ Less
Submitted 4 August, 2012;
originally announced August 2012.
-
Vertical Structure and Coronal Power of Accretion Disks Powered by MRI Turbulence
Authors:
Dmitri A. Uzdensky
Abstract:
In this paper we consider two outstanding intertwined problems in modern high-energy astrophysics: (1) the vertical thermal structure of an optically thick accretion disk heated by the dissipation of magnetohydrodynamic (MHD) turbulence driven by the magneto-rotational instability (MRI), and (2) determining the fraction of the accretion power released in the corona above the disk. For simplicity,…
▽ More
In this paper we consider two outstanding intertwined problems in modern high-energy astrophysics: (1) the vertical thermal structure of an optically thick accretion disk heated by the dissipation of magnetohydrodynamic (MHD) turbulence driven by the magneto-rotational instability (MRI), and (2) determining the fraction of the accretion power released in the corona above the disk. For simplicity, we consider a gas-pressure-dominated disk and assume a constant opacity. We argue that the local turbulent dissipation rate due to the disruption of MRI channel flows by secondary parasitic instabilities should be uniform across most of the disk, almost up to the disk photosphere. We then obtain a self-consistent analytical solution for the vertical thermal structure of the disk, governed by the balance between the heating by MRI turbulence and the cooling by radiative diffusion. Next, we argue that the coronal power fraction is determined by the competition between the Parker instability, viewed as a parasitic instability feeding off of MRI channel flows, and other parasitic instabilities. We show that the Parker instability inevitably becomes important near the disk surface, leading to a certain lower limit on the coronal power. While most of the analysis in this paper focuses on the case of a disk threaded by an externally imposed vertical magnetic field, we also discuss the zero-net-flux case, in which the magnetic field is produced by the MRI dynamo itself, and show that most of our arguments and conclusions should be valid in this case as well.
△ Less
Submitted 29 May, 2012;
originally announced May 2012.
-
Beaming and rapid variability of high-energy radiation from relativistic pair plasma reconnection
Authors:
Benoit Cerutti,
Gregory R. Werner,
Dmitri A. Uzdensky,
Mitchell C. Begelman
Abstract:
We report on the first study of the angular distribution of energetic particles and radiation generated in relativistic collisionless electron-positron pair plasma reconnection, using two-dimensional particle-in-cell simulations. We discover a strong anisotropy of the particles accelerated by reconnection and the associated strong beaming of their radiation. The focusing of particles and radiation…
▽ More
We report on the first study of the angular distribution of energetic particles and radiation generated in relativistic collisionless electron-positron pair plasma reconnection, using two-dimensional particle-in-cell simulations. We discover a strong anisotropy of the particles accelerated by reconnection and the associated strong beaming of their radiation. The focusing of particles and radiation increases with their energy; in this sense, this "kinetic beaming" effect differs fundamentally from the relativistic Doppler beaming usually invoked in high-energy astrophysics, in which all photons are focused and boosted achromatically. We also present, for the first time, the modeling of the synchrotron emission as seen by an external observer during the reconnection process. The expected lightcurves comprise several bright symmetric sub-flares emitted by the energetic beam of particles sweeping across the line of sight intermittently, and exhibit super-fast time variability as short as about one tenth of the system light-crossing time. The concentration of the energetic particles into compact regions inside magnetic islands and particle anisotropy explain the rapid variability. This radiative signature of reconnection can account for the brightness and variability of the gamma-ray flares in the Crab Nebula and in blazars.
△ Less
Submitted 3 July, 2012; v1 submitted 14 May, 2012;
originally announced May 2012.
-
Energetic Constraints on a Rapid Gamma-Ray Flare in PKS 1222+216
Authors:
Krzysztof Nalewajko,
Mitchell C. Begelman,
Benoit Cerutti,
Dmitri A. Uzdensky,
Marek Sikora
Abstract:
We study theoretical implications of a rapid Very-High-Energy (VHE) flare detected by MAGIC in the Flat-Spectrum Radio Quasar PKS 1222+216. The minimum distance from the jet origin at which this flare could be produced is 0.5 pc. A moderate Doppler factor of the VHE source, D_{VHE} ~ 20, is allowed by all opacity constraints. The concurrent High-Energy (HE) emission observed by Fermi provides esti…
▽ More
We study theoretical implications of a rapid Very-High-Energy (VHE) flare detected by MAGIC in the Flat-Spectrum Radio Quasar PKS 1222+216. The minimum distance from the jet origin at which this flare could be produced is 0.5 pc. A moderate Doppler factor of the VHE source, D_{VHE} ~ 20, is allowed by all opacity constraints. The concurrent High-Energy (HE) emission observed by Fermi provides estimates of the total jet power and the jet magnetic field strength. Energetic constraints for the VHE flare are extremely tight: for an isotropic particle distribution they require a huge co-moving energy density in the emitting region and a very efficient radiative process. We disfavor hadronic processes due to their low radiative efficiency, as well as the synchrotron scenario recently proposed for the case of HE flares in the Crab Nebula, since the parameters needed to overcome the radiative losses are quite extreme. The VHE emission can be explained by the Synchrotron Self-Compton (SSC) mechanism for D_{VHE} ~ 20 or by the External Radiation Compton (ERC) mechanism involving the infrared radiation of the dusty torus for D_{VHE} ~ 50. After discussing several alternative scenarios, we propose that the extreme energy density constraint can be satisfied when the emission comes from highly anisotropic short-lived bunches of particles formed by the kinetic beaming mechanism in magnetic reconnection sites. By focusing the emitting particles into very narrow beams, this mechanism allows one to relax the causality constraint on the source size, decreasing the required energy density by 4 orders of magnitude.
△ Less
Submitted 13 July, 2012; v1 submitted 9 February, 2012;
originally announced February 2012.
-
Extreme particle acceleration in magnetic reconnection layers. Application to the gamma-ray flares in the Crab Nebula
Authors:
Benoit Cerutti,
Dmitri A. Uzdensky,
Mitchell C. Begelman
Abstract:
The gamma-ray space telescopes AGILE and Fermi detected short and bright synchrotron gamma-ray flares at photon energies above 100 MeV in the Crab Nebula. This discovery suggests that electron-positron pairs in the nebula are accelerated to PeV energies in a milliGauss magnetic field, which is difficult to explain with classical models of particle acceleration and pulsar wind nebulae. We investiga…
▽ More
The gamma-ray space telescopes AGILE and Fermi detected short and bright synchrotron gamma-ray flares at photon energies above 100 MeV in the Crab Nebula. This discovery suggests that electron-positron pairs in the nebula are accelerated to PeV energies in a milliGauss magnetic field, which is difficult to explain with classical models of particle acceleration and pulsar wind nebulae. We investigate whether particle acceleration in a magnetic reconnection layer can account for the puzzling properties of the flares. We numerically integrate relativistic test-particle orbits in the vicinity of the layer, including the radiation reaction force, and using analytical expressions for the large-scale electromagnetic fields. As they get accelerated by the reconnection electric field, the particles are focused deep inside the current layer where the magnetic field is small. The electrons suffer less from synchrotron losses and are accelerated to extremely high energies. Population studies show that, at the end of the layer, the particle distribution piles up at the maximum energy given by the electric potential drop and is focused into a thin fan beam. Applying this model to the Crab Nebula, we find that the emerging synchrotron emission spectrum peaks above 100 MeV and is close to the spectral shape of a single electron. The flare inverse Compton emission is negligible and no detectable emission is expected at other wavelengths. This mechanism provides a plausible explanation for the gamma-ray flares in the Crab Nebula and could be at work in other astrophysical objects such as relativistic jets in active galactic nuclei.
△ Less
Submitted 5 January, 2012; v1 submitted 3 October, 2011;
originally announced October 2011.
-
Magnetic reconnection and stochastic plasmoid chains in high-Lundquist-number plasmas
Authors:
N. F. Loureiro,
R. Samtaney,
A. A. Schekochihin,
D. A. Uzdensky
Abstract:
A numerical study of magnetic reconnection in the large-Lundquist-number ($S$), plasmoid-dominated regime is carried out for $S$ up to $10^7$. The theoretical model of Uzdensky {\it et al.} [Phys. Rev. Lett. {\bf 105}, 235002 (2010)] is confirmed and partially amended. The normalized reconnection rate is $\normEeff\sim 0.02$ independently of $S$ for $S\gg10^4$. The plasmoid flux ($Ψ$) and half-wid…
▽ More
A numerical study of magnetic reconnection in the large-Lundquist-number ($S$), plasmoid-dominated regime is carried out for $S$ up to $10^7$. The theoretical model of Uzdensky {\it et al.} [Phys. Rev. Lett. {\bf 105}, 235002 (2010)] is confirmed and partially amended. The normalized reconnection rate is $\normEeff\sim 0.02$ independently of $S$ for $S\gg10^4$. The plasmoid flux ($Ψ$) and half-width ($w_x$) distribution functions scale as $f(Ψ)\sim Ψ^{-2}$ and $f(w_x)\sim w_x^{-2}$. The joint distribution of $Ψ$ and $w_x$ shows that plasmoids populate a triangular region $w_x\gtrsimΨ/B_0$, where $B_0$ is the reconnecting field. It is argued that this feature is due to plasmoid coalescence. Macroscopic "monster" plasmoids with $w_x\sim 10$% of the system size are shown to emerge in just a few Alfvén times, independently of $S$, suggesting that large disruptive events are an inevitable feature of large-$S$ reconnection.
△ Less
Submitted 19 August, 2011;
originally announced August 2011.
-
Reconnection-Powered Linear Accelerator and Gamma-Ray Flares in the Crab Nebula
Authors:
Dmitri A. Uzdensky,
Benoit Cerutti,
Mitchell C. Begelman
Abstract:
The recent discovery of day-long gamma-ray flares in the Crab Nebula, presumed to be synchrotron emission by PeV (10^{15} eV) electrons in milligauss magnetic fields, presents a strong challenge to particle acceleration models. The observed photon energies exceed the upper limit (~100 MeV) obtained by balancing the acceleration rate and synchrotron radiation losses under standard conditions where…
▽ More
The recent discovery of day-long gamma-ray flares in the Crab Nebula, presumed to be synchrotron emission by PeV (10^{15} eV) electrons in milligauss magnetic fields, presents a strong challenge to particle acceleration models. The observed photon energies exceed the upper limit (~100 MeV) obtained by balancing the acceleration rate and synchrotron radiation losses under standard conditions where the electric field is smaller than the magnetic field. We argue that a linear electric accelerator, operating at magnetic reconnection sites, is able to circumvent this difficulty. Sufficiently energetic electrons have gyroradii so large that their motion is insensitive to small-scale turbulent structures in the reconnection layer and is controlled only by large-scale fields. We show that such particles are guided into the reconnection layer by the reversing magnetic field as they are accelerated by the reconnection electric field. As these electrons become confined within the current sheet, they experience a decreasing perpendicular magnetic field that may drop below the accelerating electric field. This enables them to reach higher energies before suffering radiation losses and hence to emit synchrotron radiation in excess of the 100 MeV limit, providing a natural resolution to the Crab gamma-ray flare paradox.
△ Less
Submitted 4 May, 2011;
originally announced May 2011.
-
Magnetic Reconnection in Extreme Astrophysical Environments
Authors:
Dmitri A. Uzdensky
Abstract:
Magnetic reconnection is a basic plasma process of dramatic rearrangement of magnetic topology, often leading to a violent release of magnetic energy. It is important in magnetic fusion and in space and solar physics --- areas that have so far provided the context for most of reconnection research. Importantly, these environments consist just of electrons and ions and the dissipated energy always…
▽ More
Magnetic reconnection is a basic plasma process of dramatic rearrangement of magnetic topology, often leading to a violent release of magnetic energy. It is important in magnetic fusion and in space and solar physics --- areas that have so far provided the context for most of reconnection research. Importantly, these environments consist just of electrons and ions and the dissipated energy always stays with the plasma. In contrast, in this paper I introduce a new direction of research, motivated by several important problems in high-energy astrophysics --- reconnection in high energy density (HED) radiative plasmas, where radiation pressure and radiative cooling become dominant factors in the pressure and energy balance. I identify the key processes distinguishing HED reconnection: special-relativistic effects; radiative effects (radiative cooling, radiation pressure, and Compton resistivity); and, at the most extreme end, QED effects, including pair creation. I then discuss the main astrophysical applications --- situations with magnetar-strength fields (exceeding the quantum critical field of about 4 x 10^13 G): giant SGR flares and magnetically-powered central engines and jets of GRBs. Here, magnetic energy density is so high that its dissipation heats the plasma to MeV temperatures. Electron-positron pairs are then copiously produced, making the reconnection layer highly collisional and dressing it in a thick pair coat that traps radiation. The pressure is dominated by radiation and pairs. Yet, radiation diffusion across the layer may be faster than the global Alfvén transit time; then, radiative cooling governs the thermodynamics and reconnection becomes a radiative transfer problem, greatly affected by the ultra-strong magnetic field. This overall picture is very different from our traditional picture of reconnection and thus represents a new frontier in reconnection research.
△ Less
Submitted 12 January, 2011;
originally announced January 2011.