-
Fundamental, Harmonic, and Third-harmonic Plasma Emission from Beam-plasma Instabilities: A First-principles Precursor for Astrophysical Radio Bursts
Authors:
Fabio Bacchini,
Alexander A. Philippov
Abstract:
Electromagnetic fundamental and harmonic emission is ubiquitously observed throughout the heliosphere, and in particular it is commonly associated with the occurrence of Type II and III solar radio bursts. Classical analytic calculations for the plasma-emission process, though useful, are limited to idealized situations; a conclusive numerical verification of this theory is still lacking, with ear…
▽ More
Electromagnetic fundamental and harmonic emission is ubiquitously observed throughout the heliosphere, and in particular it is commonly associated with the occurrence of Type II and III solar radio bursts. Classical analytic calculations for the plasma-emission process, though useful, are limited to idealized situations; a conclusive numerical verification of this theory is still lacking, with earlier studies often providing contradicting results on e.g. the precise parameter space in which fundamental and harmonic emission can be produced. To accurately capture the chain of mechanisms underlying plasma emission - from precursor plasma processes to the generation of electromagnetic waves over long times - we perform large-scale, first-principles simulations of beam-plasma instabilities. By employing a very large number of computational particles we achieve very low numerical noise, and explore (with an array of simulations) a wide parameter space determined by the beam-plasma density ratio and the ion-to-electron temperature ratio. In particular, we observe direct evidence of both fundamental and harmonic plasma emission when the beam-to-background density ratio $\le$0.005 (with beam-to-background energy ratio ~0.5), tightly constraining this threshold. We observe that, asymptotically, in this regime ~0.1% of the initial beam energy is converted into harmonic emission, and ~0.001% into fundamental emission. In contrast with previous studies, we also find that this emission is independent of the ion-to-electron temperature ratio. In addition, we report the direct detection of third-harmonic emission in all of our simulations, at power levels compatible with observations. Our findings have important consequences for understanding the viable conditions leading to plasma emission in space systems, and for the interpretation of observed electromagnetic signals throughout the heliosphere.
△ Less
Submitted 16 February, 2024;
originally announced February 2024.
-
Polarized anisotropic synchrotron emission and absorption and its application to Black Hole Imaging
Authors:
Alisa Galishnikova,
Alexander Philippov,
Eliot Quataert
Abstract:
Low-collisionality plasma in a magnetic field generically develops anisotropy in its distribution function with respect to the magnetic field direction. Motivated by the application to radiation from accretion flows and jets, we explore the effect of temperature anisotropy on synchrotron emission. We derive analytically and provide numerical fits for the polarized synchrotron emission and absorpti…
▽ More
Low-collisionality plasma in a magnetic field generically develops anisotropy in its distribution function with respect to the magnetic field direction. Motivated by the application to radiation from accretion flows and jets, we explore the effect of temperature anisotropy on synchrotron emission. We derive analytically and provide numerical fits for the polarized synchrotron emission and absorption coefficients for a relativistic bi-Maxwellian plasma (we do not consider Faraday conversion/rotation). Temperature anisotropy can significantly change how the synchrotron emission and absorption coefficients depend on observing angle with respect to the magnetic field. The emitted linear polarization fraction does not depend strongly on anisotropy, while the emitted circular polarization does. We apply our results to black hole imaging of Sgr A* and M87* by ray-tracing a GRMHD simulation and assuming that the plasma temperature anisotropy is set by the thresholds of kinetic-scale anisotropy-driven instabilities. We find that the azimuthal asymmetry of the 230 GHz images can change by up to a factor of 3, accentuating ($T_\perp > T_\parallel$) or counteracting ($T_\perp < T_\parallel$) the image asymmetry produced by Doppler beaming. This can change the physical inferences from observations relative to models with an isotropic distribution function, e.g., by allowing for larger inclination between the line of sight and spin direction in Sgr A*. The observed image diameter and the size of the black hole shadow can also vary significantly due to plasma temperature anisotropy. We describe how the anisotropy of the plasma can affect future multi-frequency and photon ring observations. In Appendices we calculate kinetic anisotropy-driven instabilities (mirror, whistler, and firehose) for relativistically hot plasmas.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
The LHCb upgrade I
Authors:
LHCb collaboration,
R. Aaij,
A. S. W. Abdelmotteleb,
C. Abellan Beteta,
F. Abudinén,
C. Achard,
T. Ackernley,
B. Adeva,
M. Adinolfi,
P. Adlarson,
H. Afsharnia,
C. Agapopoulou,
C. A. Aidala,
Z. Ajaltouni,
S. Akar,
K. Akiba,
P. Albicocco,
J. Albrecht,
F. Alessio,
M. Alexander,
A. Alfonso Albero,
Z. Aliouche,
P. Alvarez Cartelle,
R. Amalric,
S. Amato
, et al. (1298 additional authors not shown)
Abstract:
The LHCb upgrade represents a major change of the experiment. The detectors have been almost completely renewed to allow running at an instantaneous luminosity five times larger than that of the previous running periods. Readout of all detectors into an all-software trigger is central to the new design, facilitating the reconstruction of events at the maximum LHC interaction rate, and their select…
▽ More
The LHCb upgrade represents a major change of the experiment. The detectors have been almost completely renewed to allow running at an instantaneous luminosity five times larger than that of the previous running periods. Readout of all detectors into an all-software trigger is central to the new design, facilitating the reconstruction of events at the maximum LHC interaction rate, and their selection in real time. The experiment's tracking system has been completely upgraded with a new pixel vertex detector, a silicon tracker upstream of the dipole magnet and three scintillating fibre tracking stations downstream of the magnet. The whole photon detection system of the RICH detectors has been renewed and the readout electronics of the calorimeter and muon systems have been fully overhauled. The first stage of the all-software trigger is implemented on a GPU farm. The output of the trigger provides a combination of totally reconstructed physics objects, such as tracks and vertices, ready for final analysis, and of entire events which need further offline reprocessing. This scheme required a complete revision of the computing model and rewriting of the experiment's software.
△ Less
Submitted 17 May, 2023;
originally announced May 2023.
-
High-Energy Radiation and Ion Acceleration in Three-dimensional Relativistic Magnetic Reconnection with Strong Synchrotron Cooling
Authors:
Alexander Chernoglazov,
Hayk Hakobyan,
Alexander A. Philippov
Abstract:
We present the results of 3D particle-in-cell (PIC) simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of non-radiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant com…
▽ More
We present the results of 3D particle-in-cell (PIC) simulations that explore relativistic magnetic reconnection in pair plasma with strong synchrotron cooling and a small mass fraction of non-radiating ions. Our results demonstrate that the structure of the current sheet is highly sensitive to the dynamic efficiency of radiative cooling. Specifically, stronger cooling leads to more significant compression of the plasma and magnetic field within the plasmoids. We demonstrate that ions can be efficiently accelerated to energies exceeding the plasma magnetization parameter, $\ggσ$, and form a hard power-law energy distribution, $f_i\propto γ^{-1}$. This conclusion implies a highly efficient proton acceleration in the magnetospheres of young pulsars. Conversely, the energies of pairs are limited to either $σ$ in the strong cooling regime or the radiation burnoff limit, $γ_{\rm syn}$, when cooling is weak. We find that the high-energy radiation from pairs above the synchrotron burnoff limit, $\varepsilon_c \approx 16$ MeV, is only efficiently produced in the strong cooling regime, $γ_{\rm syn} < σ$. In this regime, we find that the spectral cutoff scales as $\varepsilon_{\rm cut}\approx \varepsilon_c (σ/γ_{\rm syn})$, and the highest energy photons are beamed along the direction of the upstream magnetic field, consistent with the phenomenological models of gamma-ray emission from young pulsars. Furthermore, our results place constraints on the reconnection-driven models of gamma-ray flares in the Crab Nebula.
△ Less
Submitted 2 October, 2023; v1 submitted 3 May, 2023;
originally announced May 2023.
-
Cosmic ray transport in large-amplitude turbulence with small-scale field reversals
Authors:
Philipp Kempski,
Drummond B. Fielding,
Eliot Quataert,
Alisa K. Galishnikova,
Matthew W. Kunz,
Alexander A. Philippov,
Bart Ripperda
Abstract:
The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current micro-physical CR transport models in magneto-hydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying t…
▽ More
The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current micro-physical CR transport models in magneto-hydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying the alternative regime of large-amplitude turbulence with $δB/B_0 \gg 1$, in which intermittent small-scale magnetic field reversals are ubiquitous. We study particle transport in such turbulence by integrating trajectories in stationary snapshots. To quantify spatial diffusion, we use a setup with continuous particle injection and escape, which we term the turbulent leaky box. We find that particle transport is very different from the strong-guide-field case. Low-energy particles are better confined than high-energy particles, despite less efficient pitch-angle isotropization at small energies. In the limit of weak guide field, energy-dependent confinement is driven by the energy-dependent (in)ability to follow reversing magnetic field lines exactly and by the scattering in regions of ``resonant curvature", where the field line bends on a scale that is of order the local particle gyro-radius. We derive a heuristic model of particle transport in magnetic folds that approximately reproduces the energy dependence of transport found numerically. We speculate that CR propagation in the Galaxy is regulated by the intermittent field reversals highlighted here and discuss the implications of our findings for CR transport in the Milky Way.
△ Less
Submitted 19 December, 2023; v1 submitted 24 April, 2023;
originally announced April 2023.
-
Three-dimensional dynamics of strongly twisted magnetar magnetospheres: Kinking flux tubes and global eruptions
Authors:
J. F. Mahlmann,
A. A. Philippov,
V. Mewes,
B. Ripperda,
E. R. Most,
L. Sironi
Abstract:
The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric…
▽ More
The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, as to say, catastrophic lateral instabilities triggered by a critical footpoint displacement of $ψ_{\rm crit}\gtrsimπ$. In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to $25\%$ of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, $\lesssim 10^{43}$erg. Such events occur when significant energy builds up deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind -- for example, via relativistic magnetic reconnection or shocks.
△ Less
Submitted 25 April, 2023; v1 submitted 14 February, 2023;
originally announced February 2023.
-
Particle acceleration by magnetic Rayleigh-Taylor instability: mechanism for flares in black-hole accretion flows
Authors:
Vladimir Zhdankin,
Bart Ripperda,
Alexander A. Philippov
Abstract:
We study the magnetic Rayleigh-Taylor instability in relativistic collisionless plasma, as an astrophysical process for nonthermal particle acceleration. We consider dense plasma on top of a highly magnetized cavity with sheared magnetic field. Using particle-in-cell simulations, we show that small plumes grow and merge progressively to form a large-scale plume, which broadens to drive rapid magne…
▽ More
We study the magnetic Rayleigh-Taylor instability in relativistic collisionless plasma, as an astrophysical process for nonthermal particle acceleration. We consider dense plasma on top of a highly magnetized cavity with sheared magnetic field. Using particle-in-cell simulations, we show that small plumes grow and merge progressively to form a large-scale plume, which broadens to drive rapid magnetic reconnection in the cavity. We find that this leads to efficient particle acceleration capable of explaining flares from the inner accretion flow onto the black hole Sgr A*.
△ Less
Submitted 10 October, 2023; v1 submitted 10 February, 2023;
originally announced February 2023.
-
Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae
Authors:
Daniel Groselj,
Hayk Hakobyan,
Andrei M. Beloborodov,
Lorenzo Sironi,
Alexander Philippov
Abstract:
We report results from the first radiative particle-in-cell simulations of strong Alfvenic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accr…
▽ More
We report results from the first radiative particle-in-cell simulations of strong Alfvenic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. The method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.
△ Less
Submitted 25 January, 2024; v1 submitted 26 January, 2023;
originally announced January 2023.
-
Magnetic energy dissipation and gamma-ray emission in energetic pulsars
Authors:
Hayk Hakobyan,
Alexander Philippov,
Anatoly Spitkovsky
Abstract:
Some of the most energetic pulsars exhibit rotation-modulated gamma-ray emission in the $0.1$ to $100$ GeV band. The luminosity of this emission is typically $0.1\text{-}10\%$ of the pulsar spin-down power (gamma-ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this…
▽ More
Some of the most energetic pulsars exhibit rotation-modulated gamma-ray emission in the $0.1$ to $100$ GeV band. The luminosity of this emission is typically $0.1\text{-}10\%$ of the pulsar spin-down power (gamma-ray efficiency), implying that a significant fraction of the available electromagnetic energy is dissipated in the magnetosphere and reradiated as high-energy photons. To investigate this phenomenon we model a pulsar magnetosphere using 3D particle-in-cell simulations with strong synchrotron cooling. We particularly focus on the dynamics of the equatorial current sheet where magnetic reconnection and energy dissipation take place. Our simulations demonstrate that a fraction of the spin-down power dissipated in the magnetospheric current sheet is controlled by the rate of magnetic reconnection at microphysical plasma scales and only depends on the pulsar inclination angle. We demonstrate that the maximum energy and the distribution function of accelerated pairs is controlled by the available magnetic energy per particle near the current sheet, the magnetization parameter. The shape and the extent of the plasma distribution is imprinted in the observed synchrotron emission, in particular, in the peak and the cutoff of the observed spectrum. We study how the strength of synchrotron cooling affects the observed variety of spectral shapes. Our conclusions naturally explain why pulsars with higher spin-down power have wider spectral shapes and, as a result, lower gamma-ray efficiency.
△ Less
Submitted 17 March, 2023; v1 submitted 5 September, 2022;
originally announced September 2022.
-
Radiative reconnection-powered TeV flares from the black hole magnetosphere in M87
Authors:
Hayk Hakobyan,
Bart Ripperda,
Alexander Philippov
Abstract:
Active Galactic Nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region nearby the nucleus. However, the e…
▽ More
Active Galactic Nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region nearby the nucleus. However, the emission mechanism for these flares is not well understood. Recent high-resolution general-relativistic magnetohydrodynamic simulations show the occurrence of episodic magnetic reconnection events that can power flares nearby the black hole event horizon. In this work, we analyze the radiative properties of the reconnecting current layer under the extreme plasma conditions applicable to the black hole in M87 from the first principles. We show that abundant pair production is expected in the vicinity of the reconnection layer, to the extent that the produced secondary pair plasma dominates the reconnection dynamics. Using analytic estimates backed by two-dimensional particle-in-cell simulations we demonstrate that in the presence of strong synchrotron cooling, reconnection can produce a hard power-law distribution of pair plasma imprinted in the outgoing synchrotron (up to a few tens of MeV) and the inverse-Compton signal (up to TeV). We produce synthetic radiation spectra from our simulations, which can be directly compared with the results of future multiwavelength observations of M87* flares.
△ Less
Submitted 17 March, 2023; v1 submitted 5 September, 2022;
originally announced September 2022.
-
Electric field screening in pair discharges and generation of pulsar radio emission
Authors:
Elizabeth A. Tolman,
Alexander A. Philippov,
Andrey N. Timokhin
Abstract:
Pulsar radio emission may be generated in pair discharges which fill the pulsar magnetosphere with plasma as an accelerating electric field is screened by freshly created pairs. In this Letter we develop a simplified analytic theory for the screening of the electric field in these pair discharges and use it to estimate total radio luminosity and spectrum. The discharge has three stages. First, the…
▽ More
Pulsar radio emission may be generated in pair discharges which fill the pulsar magnetosphere with plasma as an accelerating electric field is screened by freshly created pairs. In this Letter we develop a simplified analytic theory for the screening of the electric field in these pair discharges and use it to estimate total radio luminosity and spectrum. The discharge has three stages. First, the electric field is screened for the first time and starts to oscillate. Next, a nonlinear phase occurs. In this phase, the amplitude of the electric field experiences strong damping because the field dramatically changes the momenta of newly created pairs. This strong damping ceases, and the system enters a final linear phase, when the electric field can no longer dramatically change pair momenta. Applied to pulsars, this theory may explain several aspects of radio emission, including the observed luminosity, $L_{\rm{rad}} \sim 10^{28} \rm{erg} \, \rm{s}^{-1}$, and the observed spectrum, $S_ω\sim ω^{-1.4 \pm 1.0} $.
△ Less
Submitted 29 June, 2022; v1 submitted 2 February, 2022;
originally announced February 2022.
-
The role of plasma instabilities in relativistic radiation mediated shocks: stability analysis and particle-in-cell simulations
Authors:
Arno Vanthieghem,
Jens F. Mahlmann,
Amir Levinson,
Alexander A. Philippov,
Ehud Nakar,
Frederico Fiuza
Abstract:
Relativistic radiation mediated shocks (RRMS) likely form in prodigious cosmic explosions. The structure and emission of such shocks is regulated by copious production of electron-positron pairs inside the shock transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is ex…
▽ More
Relativistic radiation mediated shocks (RRMS) likely form in prodigious cosmic explosions. The structure and emission of such shocks is regulated by copious production of electron-positron pairs inside the shock transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the nonlinear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced transport equation for the particles via pitch angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to nonadiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of nonthermal distributions, arising from the microturbulence, can affect the observed shock breakout signal in ways unaccounted for by current single-fluid models.
△ Less
Submitted 14 January, 2022;
originally announced January 2022.
-
Dynamic alignment and plasmoid formation in relativistic magnetohydrodynamic turbulence
Authors:
Alexander Chernoglazov,
Bart Ripperda,
Alexander A. Philippov
Abstract:
We present high resolution 2D and 3D simulations of magnetized decaying turbulence in relativistic resistive magneto-hydrodynamics. The simulations show dynamic formation of large scale intermittent long-lived current sheets being disrupted by the tearing instability into plasmoid chains. These current sheets are locations of enhanced magnetic field dissipation and heating of the plasma. We find m…
▽ More
We present high resolution 2D and 3D simulations of magnetized decaying turbulence in relativistic resistive magneto-hydrodynamics. The simulations show dynamic formation of large scale intermittent long-lived current sheets being disrupted by the tearing instability into plasmoid chains. These current sheets are locations of enhanced magnetic field dissipation and heating of the plasma. We find magnetic energy spectra $\propto k^{-3/2}$, together with strongly pronounced dynamic alignment of Elsasser fields and of velocity and magnetic fields, for strong guide-field turbulence, whereas we retrieve spectra $\propto k^{-5/3}$ for the case of a weak guide-field.
△ Less
Submitted 15 November, 2021;
originally announced November 2021.
-
Black hole flares: ejection of accreted magnetic flux through 3D plasmoid-mediated reconnection
Authors:
Bart Ripperda,
Matthew Liska,
Koushik Chatterjee,
Gibwa Musoke,
Alexander A. Philippov,
Sera B. Markoff,
Alexander Tchekhovskoy,
Ziri Younsi
Abstract:
Magnetic reconnection can power bright, rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high resolution ($5376\times2304\times2304$ cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current she…
▽ More
Magnetic reconnection can power bright, rapid flares originating from the inner magnetosphere of accreting black holes. We conduct extremely high resolution ($5376\times2304\times2304$ cells) general-relativistic magnetohydrodynamics simulations, capturing plasmoid-mediated reconnection in a 3D magnetically arrested disk for the first time. We show that an equatorial, plasmoid-unstable current sheet forms in a transient, non-axisymmetric, low-density magnetosphere within the inner few Schwarzschild radii. Magnetic flux bundles escape from the event horizon through reconnection at the universal plasmoid-mediated rate in this current sheet. The reconnection feeds on the highly-magnetized plasma in the jets and heats the plasma that ends up trapped in flux bundles to temperatures proportional to the jet's magnetization. The escaped flux bundles can complete a full orbit as low-density hot spots, consistent with Sgr A$^{*}$ observations by the GRAVITY interferometer. Reconnection near the horizon produces sufficiently energetic plasma to explain flares from accreting black holes, such as the TeV emission observed from M87. The drop in mass accretion rate during the flare, and the resulting low-density magnetosphere make it easier for very high energy photons produced by reconnection-accelerated particles to escape. The extreme resolution results in a converged plasmoid-mediated reconnection rate that directly determines the timescales and properties of the flare.
△ Less
Submitted 6 February, 2022; v1 submitted 30 September, 2021;
originally announced September 2021.
-
Magnetic Hair and Reconnection in Black Hole Magnetospheres
Authors:
Ashley Bransgrove,
Bart Ripperda,
Alexander Philippov
Abstract:
The no-hair theorem of general relativity states that isolated black holes are characterized by three parameters: mass, spin, and charge. In this Letter we consider Kerr black holes endowed with highly magnetized plasma-filled magnetospheres. Using general relativistic kinetic plasma and resistive magnetohydrodynamics simulations, we show that a dipole magnetic field on the event horizon opens int…
▽ More
The no-hair theorem of general relativity states that isolated black holes are characterized by three parameters: mass, spin, and charge. In this Letter we consider Kerr black holes endowed with highly magnetized plasma-filled magnetospheres. Using general relativistic kinetic plasma and resistive magnetohydrodynamics simulations, we show that a dipole magnetic field on the event horizon opens into a split-monopole and reconnects in a plasmoid-unstable current-sheet. The no-hair theorem is satisfied, in the sense that all components of the stress-energy tensor decay exponentially in time. We measure the decay time of magnetic flux on the event horizon for plasmoid-dominated reconnection in collisionless and collisional plasma. The reconnecting magnetosphere should be a powerful source of hard X-ray emission when the magnetic field is strong.
△ Less
Submitted 29 September, 2021;
originally announced September 2021.
-
Scaling of Small-Scale Dynamo Properties in the Rayleigh-Taylor Instability
Authors:
Valentin Skoutnev,
Elias R. Most,
Amitava Bhattacharjee,
Alexander A. Philippov
Abstract:
We derive scaling relations based on free-fall and isotropy assumptions for the kinematic small-scale dynamo growth rate and amplification factor over the course of the mixing, saturation, and decay phases of the Rayleigh-Taylor instability (RTI) in a fully-ionized plasma. The scaling relations are tested using sets of three dimensional, visco-resistive MHD simulations of the RTI and found to hold…
▽ More
We derive scaling relations based on free-fall and isotropy assumptions for the kinematic small-scale dynamo growth rate and amplification factor over the course of the mixing, saturation, and decay phases of the Rayleigh-Taylor instability (RTI) in a fully-ionized plasma. The scaling relations are tested using sets of three dimensional, visco-resistive MHD simulations of the RTI and found to hold in the saturation phase, but exhibit discrepancies during the mixing and decays phases, suggesting a need to relax either the free-fall or isotropy assumptions. Application of the scaling relations allows for quantitative prediction of the net amplification of magnetic energy in the kinematic dynamo phase and therefore a determination of whether the magnetic energy either remains sub-equipartition at all velocity scales or reaches equipartition with at least some scales of the turbulent kinetic energy in laboratory and astrophysical scenarios. As an example, we consider the dynamo in RTI-unstable regions of the outer envelope of a binary neutron star merger and predict that the kinematic regime of the small-scale dynamo ends on the time scale of nanoseconds and then reaches saturation on a timescale of microseconds, which are both fast compared to the millisecond relaxation time of the post-merger.
△ Less
Submitted 3 November, 2021; v1 submitted 8 June, 2021;
originally announced June 2021.
-
Weak Alfvénic turbulence in relativistic plasmas. Part 1. Dynamical equations and basic dynamics of interacting resonant triads
Authors:
J. M. TenBarge,
B. Ripperda,
A. Chernoglazov,
A. Bhattacharjee,
J. F. Mahlmann,
E. R. Most,
J. Juno,
Y. Yuan,
A. A. Philippov
Abstract:
Alfvén wave collisions are the primary building blocks of the non-relativistic turbulence that permeates the heliosphere and low-to-moderate energy astrophysical systems. However, many astrophysical systems such as gamma-ray bursts, pulsar and magnetar magnetospheres, and active galactic nuclei have relativistic flows or energy densities. To better understand these high energy systems, we derive r…
▽ More
Alfvén wave collisions are the primary building blocks of the non-relativistic turbulence that permeates the heliosphere and low-to-moderate energy astrophysical systems. However, many astrophysical systems such as gamma-ray bursts, pulsar and magnetar magnetospheres, and active galactic nuclei have relativistic flows or energy densities. To better understand these high energy systems, we derive reduced relativistic MHD equations and employ them to examine weak Alfvénic turbulence, dominated by three-wave interactions, in reduced relativistic magnetohydrodynamics, including the force-free, infinitely magnetized limit. We compare both numerical and analytical solutions to demonstrate that many of the findings from non-relativistic weak turbulence are retained in the relativistic system. But, an important distinction in the relativistic limit is the inapplicability of a formally incompressible limit, i.e, there exists finite coupling to the compressible fast mode regardless of the strength of the magnetic field. Since fast modes can propagate across field lines, this mechanism provides a route for energy to escape strongly magnetized systems, e.g., magnetar magnetospheres. However, we find that the fast-Alfvén coupling is diminished in the limit of oblique propagation.
△ Less
Submitted 23 February, 2022; v1 submitted 3 May, 2021;
originally announced May 2021.
-
Weak Alfvénic turbulence in relativistic plasmas II: Current sheets and dissipation
Authors:
B. Ripperda,
J. F. Mahlmann,
A. Chernoglazov,
J. M. TenBarge,
E. R. Most,
J. Juno,
Y. Yuan,
A. A. Philippov,
A. Bhattacharjee
Abstract:
Alfvén waves as excited in black hole accretion disks and neutron star magnetospheres are the building blocks of turbulence in relativistic, magnetized plasmas. A large reservoir of magnetic energy is available in these systems, such that the plasma can be heated significantly even in the weak turbulence regime. We perform high-resolution three-dimensional simulations of counter-propagating Alfvén…
▽ More
Alfvén waves as excited in black hole accretion disks and neutron star magnetospheres are the building blocks of turbulence in relativistic, magnetized plasmas. A large reservoir of magnetic energy is available in these systems, such that the plasma can be heated significantly even in the weak turbulence regime. We perform high-resolution three-dimensional simulations of counter-propagating Alfvén waves, showing that an $E_{B_{\perp}}(k_{\perp}) \propto k_{\perp}^{-2}$ energy spectrum develops as a result of the weak turbulence cascade in relativistic magnetohydrodynamics and its infinitely magnetized (force-free) limit. The plasma turbulence ubiquitously generates current sheets, which act as locations where magnetic energy dissipates. We show that current sheets form as a natural result of nonlinear interactions between counter-propagating Alfvén waves. These current sheets form due to the compression of elongated eddies, driven by the shear induced by growing higher order modes, and undergo a thinning process until they break-up into small-scale turbulent structures. We explore the formation of {current sheets} both in overlapping waves and in localized wave packet collisions. The relativistic interaction of localized Alfvén waves induces both Alfvén waves and fast waves and efficiently mediates the conversion and dissipation of electromagnetic energy in astrophysical systems. Plasma energization through reconnection in current sheets emerging during the interaction of Alfvén waves can potentially explain X-ray emission in black hole accretion coronae and neutron star magnetospheres.
△ Less
Submitted 25 October, 2021; v1 submitted 3 May, 2021;
originally announced May 2021.
-
Dissipation of the striped pulsar wind and non-thermal particle acceleration: 3D PIC simulations
Authors:
Benoît Cerutti,
Alexander Philippov,
Guillaume Dubus
Abstract:
The formation of a large-scale current sheet is a generic feature of pulsar magnetospheres. If the magnetic axis is misaligned with the star rotation axis, the current sheet is an oscillatory structure filling an equatorial wedge determined by the inclination angle, known as the striped wind. Relativistic reconnection could lead to significant dissipation of magnetic energy and particle accelerati…
▽ More
The formation of a large-scale current sheet is a generic feature of pulsar magnetospheres. If the magnetic axis is misaligned with the star rotation axis, the current sheet is an oscillatory structure filling an equatorial wedge determined by the inclination angle, known as the striped wind. Relativistic reconnection could lead to significant dissipation of magnetic energy and particle acceleration although the efficiency of this process is debated in this context. In this study, we aim at reconciling global models of pulsar wind dynamics and reconnection in the stripes within the same numerical framework, in order to shed new light on dissipation and particle acceleration in pulsar winds. To this end, we perform large three-dimensional particle-in-cell simulations of a split-monopole magnetosphere, from the stellar surface up to fifty light-cylinder radii away from the pulsar. Plasmoid-dominated reconnection efficiently fragments the current sheet into a dynamical network of interacting flux ropes separated by secondary current sheets which consume the field efficiently at all radii, even past the fast magnetosonic point. Our results suggest there is a universal dissipation radius solely determined by the reconnection rate in the sheet, lying well upstream the termination shock radius in isolated pair producing pulsars. The wind bulk Lorentz factor is much less relativistic than previously thought. In the comoving frame, the wind is composed of hot pairs trapped within flux ropes with a hard broad power-law spectrum, whose maximum energy is limited by the magnetization of the wind at launch. We conclude that the striped wind is most likely fully dissipated when it enters the pulsar wind nebula. The predicted wind particle spectrum after dissipation is reminiscent of the Crab Nebula radio-emitting electrons.
△ Less
Submitted 26 August, 2020;
originally announced August 2020.
-
A coupled guiding center-Boris particle pusher for magnetized plasmas in compact-object magnetospheres
Authors:
Fabio Bacchini,
Bart Ripperda,
Alexander A. Philippov,
Kyle Parfrey
Abstract:
We present a novel numerical scheme for simulating the motion of relativistic charged particles in magnetospheres of compact objects, typically filled with highly magnetized collisionless plasmas. The new algorithm is based on a dynamic switch between the full system of equations of motion and a guiding center approximation. The switch between the two formulations is based on the magnetization of…
▽ More
We present a novel numerical scheme for simulating the motion of relativistic charged particles in magnetospheres of compact objects, typically filled with highly magnetized collisionless plasmas. The new algorithm is based on a dynamic switch between the full system of equations of motion and a guiding center approximation. The switch between the two formulations is based on the magnetization of the plasma particles, such that the dynamics are accurately captured by the guiding center motion even when the gyro-frequency is under-resolved by the time step. For particles with a large gyro-radius, due to acceleration in, e.g., reconnecting current sheets, the algorithm adaptively switches to solve the full equations of motion instead. The new scheme is directly compatible with standard Particle-in-Cell codes, and is readily applicable in curved spacetimes via a dedicated covariant formulation. We test the performance of the coupled algorithm by evolving charged particles in electromagnetic configurations of reconnecting current sheets in magnetized plasma, obtained from special- and general-relativistic Particle-in-Cell simulations. The new coupled pusher is capable of producing highly accurate particle trajectories even when the time step is many orders of magnitude larger than the gyro-period, substantially reducing the restrictions of the temporal resolution.
△ Less
Submitted 3 November, 2020; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Magnetic Reconnection and Hot Spot Formation in Black Hole Accretion Disks
Authors:
Bart Ripperda,
Fabio Bacchini,
Alexander Philippov
Abstract:
Hot spots, or plasmoids, forming due to magnetic reconnection in current sheets, are conjectured to power frequent X-ray and near-infrared flares from Sgr A*, the black hole in the center of our Galaxy. It is unclear how, where, and when current sheets form in black-hole accretion disks. We perform axisymmetric general-relativistic resistive magnetohydrodynamics simulations to model reconnection a…
▽ More
Hot spots, or plasmoids, forming due to magnetic reconnection in current sheets, are conjectured to power frequent X-ray and near-infrared flares from Sgr A*, the black hole in the center of our Galaxy. It is unclear how, where, and when current sheets form in black-hole accretion disks. We perform axisymmetric general-relativistic resistive magnetohydrodynamics simulations to model reconnection and plasmoid formation in a range of accretion flows. Current sheets and plasmoids are ubiquitous features which form regardless of the initial magnetic field in the disk, the magnetization in the quasi-steady-state phase of accretion, and the spin of the black hole. Within 10 Schwarzschild radii from the event horizon, we observe plasmoids forming, after which they can merge, grow to macroscopic scales of the order of a few Schwarzschild radii, and are ultimately advected along the jet's sheath or into the disk. Large plasmoids are energized to relativistic temperatures via reconnection and contribute to the jet's limb-brightening. We find that only hot spots forming in magnetically arrested disks can potentially explain the energetics of Sgr A* flares. The flare period is determined by the reconnection rate, which we find to be between 0.01c and 0.03c in all cases, consistent with studies of reconnection in isolated Harris-type current sheets. We quantify magnetic dissipation and non-ideal electric fields which can efficiently inject non-thermal particles. We show that explicit resistivity allows for converged numerical solutions, such that the electromagnetic energy evolution and dissipation become independent of the grid scale for the extreme resolutions considered here.
△ Less
Submitted 7 September, 2020; v1 submitted 9 March, 2020;
originally announced March 2020.
-
Particle acceleration in kink-unstable jets
Authors:
Jordy Davelaar,
Alexander A. Philippov,
Omer Bromberg,
Chandra B. Singh
Abstract:
Magnetized jets in GRBs and AGNs are thought to be efficient accelerators of particles, however, the process responsible for the acceleration is still a matter of active debate. In this work, we study the kink-instability in non-rotating force-free jets using first-principle particle-in-cell simulations. We obtain similar overall evolution of the instability as found in MHD simulations. The instab…
▽ More
Magnetized jets in GRBs and AGNs are thought to be efficient accelerators of particles, however, the process responsible for the acceleration is still a matter of active debate. In this work, we study the kink-instability in non-rotating force-free jets using first-principle particle-in-cell simulations. We obtain similar overall evolution of the instability as found in MHD simulations. The instability first generates large scale current sheets, which at later times break up into small-scale turbulence. Reconnection in these sheets proceeds in the strong guide field regime, which results in a formation of steep power laws in the particle spectra. Later evolution shows heating of the plasma, which is driven by small-amplitude turbulence induced by the kink instability. These two processes energize particles due to a combination of ideal and non-ideal electric fields.
△ Less
Submitted 26 May, 2020; v1 submitted 29 October, 2019;
originally announced October 2019.
-
Extreme Plasma Astrophysics
Authors:
D. Uzdensky,
M. Begelman,
A. Beloborodov,
R. Blandford,
S. Boldyrev,
B. Cerutti,
F. Fiuza,
D. Giannios,
T. Grismayer,
M. Kunz,
N. Loureiro,
M. Lyutikov,
M. Medvedev,
M. Petropoulou,
A. Philippov,
E. Quataert,
A. Schekochihin,
K. Schoeffler,
L. Silva,
L. Sironi,
A. Spitkovsky,
G. Werner,
V. Zhdankin,
J. Zrake,
E. Zweibel
Abstract:
This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.
This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.
△ Less
Submitted 20 July, 2019; v1 submitted 13 March, 2019;
originally announced March 2019.
-
First-Principles Plasma Simulations of Black-Hole Jet Launching
Authors:
Kyle Parfrey,
Alexander Philippov,
Benoit Cerutti
Abstract:
Black holes drive powerful plasma jets to relativistic velocities. This plasma should be collisionless, and self-consistently supplied by pair creation near the horizon. We present general-relativistic collisionless plasma simulations of Kerr-black-hole magnetospheres which begin from vacuum, inject electron-positron pairs based on local unscreened electric fields, and reach steady states with ele…
▽ More
Black holes drive powerful plasma jets to relativistic velocities. This plasma should be collisionless, and self-consistently supplied by pair creation near the horizon. We present general-relativistic collisionless plasma simulations of Kerr-black-hole magnetospheres which begin from vacuum, inject electron-positron pairs based on local unscreened electric fields, and reach steady states with electromagnetically powered Blandford-Znajek jets and persistent current sheets. Particles with negative energy-at-infinity are a general feature, and can contribute significantly to black-hole rotational-energy extraction in a variant of the Penrose process. The generated plasma distribution depends on the pair-creation environment, and we describe two distinct realizations of the force-free electrodynamic solution. This sensitivity suggests that plasma kinetics will be useful in interpreting future horizon-resolving submillimeter and infrared observations.
△ Less
Submitted 10 October, 2018; v1 submitted 8 October, 2018;
originally announced October 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.