Showing 1–50 of 72 results for author: Philippov, A

Emergence of a radio jet in the changing-look AGN 1ES 1927+654

Authors: Eileen T. Meyer, Sibasish Laha, Onic I. Shuvo, Agniva Roychowdhury, David A. Green, Lauren Rhodes, Amelia M. Hankla, Alexander Philippov, Rostom Mbarek, Ari laor, Mitchell C. Begelman, Dev R. Sadaula, Ritesh Ghosh, Gabriele Bruni, Francesca Panessa, Matteo Guainazzi, Ehud Behar, Megan Masterson, Haocheng Zhang, Xiaolong Yang, Mark A. Gurwell, Garrett K. Keating, David Williams-Baldwin, Justin D. Bray, Emmanuel K. Bempong-Manful , et al. (10 additional authors not shown)

Abstract: We present multi-frequency (5-345 GHz) and multi-resolution radio observations of 1ES 1927+654, widely considered one of the most unusual and extreme changing-look active galactic nuclei (CL-AGN). The source was first designated a CL-AGN after an optical outburst in late 2017 and has since displayed considerable changes in X-ray emission, including the destruction and rebuilding of the X-ray coron… ▽ More We present multi-frequency (5-345 GHz) and multi-resolution radio observations of 1ES 1927+654, widely considered one of the most unusual and extreme changing-look active galactic nuclei (CL-AGN). The source was first designated a CL-AGN after an optical outburst in late 2017 and has since displayed considerable changes in X-ray emission, including the destruction and rebuilding of the X-ray corona in 2019-2020. Radio observations prior to 2023 show a faint and compact radio source not unusual for radio-quiet AGN. Starting in February 2023, 1ES 1927+654 began exhibiting a radio flare with a steep exponential rise, reaching a peak 60 times previous flux levels, and has maintained this higher level of radio emission for nearly a year. The 5-23 GHz spectrum is broadly similar to gigahertz-peaked radio sources, which are understood to be young radio jets less than ~1000 years old. Recent high-resolution VLBA observations at 23.5 GHz now show resolved extensions on either side of the core, with a separation of 0.14 pc, consistent with a new and mildly relativistic bipolar outflow. A steady increase in the soft X-ray band (0.3-2 keV) concurrent with the radio may be consistent with jet-driven shocked gas, though further observations are needed to test alternate scenarios. This source joins a growing number of CL-AGN and tidal disruption events which show late-time radio activity, years after the initial outburst. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: Submitted to ApJ Letters 24 June 2023. 15 pages, 4 figures, 1 table

Radio polarization of millisecond pulsars with multipolar magnetic fields

Authors: Ankan Sur, Yajie Yuan, Alexander Philippov

Abstract: NICER has observed a few millisecond pulsars where the geometry of the X-ray emitting hotspots on the neutron star is analyzed in order to constrain the mass and radius from X-ray light curve modeling. One example, PSR J0030+0451, is shown to possibly have significant multipolar magnetic fields at the stellar surface. Using force-free simulations of the magnetosphere structure, it has been shown t… ▽ More NICER has observed a few millisecond pulsars where the geometry of the X-ray emitting hotspots on the neutron star is analyzed in order to constrain the mass and radius from X-ray light curve modeling. One example, PSR J0030+0451, is shown to possibly have significant multipolar magnetic fields at the stellar surface. Using force-free simulations of the magnetosphere structure, it has been shown that the radio, X-ray, and gamma-ray light curves can be modeled simultaneously with appropriate field configuration. An even more stringent test is to compare predictions of the force-free magnetosphere model with observations of the radio polarization. This paper attempts to reproduce the radio polarization of PSR J0030+0451 using a force-free magnetospheric solution. As a result of our modeling, we can reproduce certain features of the polarization well. △ Less

Submitted 9 April, 2024; v1 submitted 18 February, 2024; originally announced February 2024.

Comments: 12 pages, 8 figures, matched version of ApJ

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.

Comments: Accepted in MNRAS

Pulsar Nulling and Vacuum Radio Emission from Axion Clouds

Authors: Andrea Caputo, Samuel J. Witte, Alexander A. Philippov, Ted Jacobson

Abstract: Non-relativistic axions can be efficiently produced in in the polar caps of pulsars, resulting in the formation of a dense cloud of gravitationally bound axions. Here, we investigate the interplay between such an axion cloud and the electrodynamics in the pulsar magnetosphere, focusing specifically on the dynamics in the polar caps, where the impact of the axion cloud is expected to be most pronou… ▽ More Non-relativistic axions can be efficiently produced in in the polar caps of pulsars, resulting in the formation of a dense cloud of gravitationally bound axions. Here, we investigate the interplay between such an axion cloud and the electrodynamics in the pulsar magnetosphere, focusing specifically on the dynamics in the polar caps, where the impact of the axion cloud is expected to be most pronounced. For sufficiently light axions $m_a \lesssim 10^{-7}$ eV, we show that the axion cloud can occasionally screen the local electric field responsible for particle acceleration and pair production, inducing a periodic nulling of the pulsar's intrinsic radio emission. At larger axion masses, the small-scale fluctuations in the axion field tend to suppress the back-reaction of the axion on the electrodynamics; however, we point out that the incoherent oscillations of the axion in short-lived regions of vacuum near the neutron star surface can produce a narrow radio line, which provides a complementary source of radio emission to the plasma-resonant emission processes identified in previous work. While this work focuses on the leading order correction to pair production in the magnetosphere, we speculate that there can exist dramatic deviations in the electrodynamics of these systems when the axion back-reaction becomes non-linear. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: v1: 6 pages, 1 figure, Supplementary Material (7 pages)

The Interplay between accelerated Protons, X-rays and Neutrinos in the Corona of NGC 1068: Constraints from Kinetic Plasma Simulations

Authors: Rostom Mbarek, Alexander Philippov, Alexander Chernoglazov, Amir Levinson, Richard Mushotzky

Abstract: We examine properties of accelerated protons potentially responsible for the neutrino excess observed in the direction of NGC 1068, using constraints from kinetic Particle-in-Cell (PIC) simulations. We find that i) coronal X-rays and Optical/Ultra-Violet light in the inner disk lead to efficient absorption of hadronic $γ$-rays within 100 Schwarzschild radii from the black hole; ii) protons acceler… ▽ More We examine properties of accelerated protons potentially responsible for the neutrino excess observed in the direction of NGC 1068, using constraints from kinetic Particle-in-Cell (PIC) simulations. We find that i) coronal X-rays and Optical/Ultra-Violet light in the inner disk lead to efficient absorption of hadronic $γ$-rays within 100 Schwarzschild radii from the black hole; ii) protons accelerated from the coronal thermal pool cannot account for the observed neutrinos; and iii) explaining the observed signal requires an injection of protons with a hard spectrum, peaking at $γ_p\sim 10^3-10^4$, into the turbulent magnetically-dominated corona, where they are confined and re-accelerated. The resulting neutrino signal can be consistent with IceCube observations. In our most favorable scenario, the injected protons are pre-accelerated in intermittent current sheets in the vicinity of the black hole, occurring either at the boundary between the disk and the outflow or during magnetic flux eruption events. △ Less

Submitted 24 April, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 11 pages with 8 figures, along with Supplementary Material; Accepted

Journal ref: Phys. Rev. D 109, L101306 Published 24 May 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.

Comments: 31 pages, 17 figures, 1 table. Accepted for publication in ApJ

Synchrotron polarization signatures of surface waves in supermassive black hole jets

Authors: Jordy Davelaar, Bart Ripperda, Lorenzo Sironi, Alexander A. Philippov, Hector Olivares, Oliver Porth, Bram van den Berg, Thomas Bronzwaer, Koushik Chatterjee, Matthew Liska

Abstract: Supermassive black holes in active galactic nuclei (AGN) are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and are thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of non-thermal electrons. In this Letter, we present a global general relat… ▽ More Supermassive black holes in active galactic nuclei (AGN) are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and are thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of non-thermal electrons. In this Letter, we present a global general relativistic magnetohydrodynamical (GRMHD) simulation of a magnetically arrested disk (MAD). After the simulation reaches the MAD state, we show that waves are continuously launched from the vicinity of the black hole and propagate along the interface between the jet and the wind. At this interface, a steep gradient in velocity is present between the mildly relativistic wind and the highly relativistic jet. The interface is, therefore, a shear layer, and due to the shear, the waves generate roll-ups that alter the magnetic field configuration and the shear layer geometry. We then perform polarized radiation transfer calculations of our GRMHD simulation and find signatures of the waves in both total intensity and linear polarization, effectively lowering the fully resolved polarization fraction. The tell-tale polarization signatures of the waves could be observable by future Very Long Baseline Interferometric observations, e.g., by the next-generation Event Horizon Telescope. △ Less

Submitted 24 October, 2023; v1 submitted 14 September, 2023; originally announced September 2023.

Comments: 20 pages, 17 figures, accepted for publication in ApJL

Electromagnetic precursors to black hole - neutron star gravitational wave events: Flares and reconnection-powered fast-radio transients from the late inspiral

Authors: Elias R. Most, Alexander A. Philippov

Abstract: The presence of magnetic fields in the late inspiral of black hole -- neutron star binaries could lead to potentially detectable electromagnetic precursor transients. Using general-relativistic force-free electrodynamics simulations, we investigate pre-merger interactions of the common magnetosphere of black hole -- neutron star systems. We demonstrate that these systems can feature copious electr… ▽ More The presence of magnetic fields in the late inspiral of black hole -- neutron star binaries could lead to potentially detectable electromagnetic precursor transients. Using general-relativistic force-free electrodynamics simulations, we investigate pre-merger interactions of the common magnetosphere of black hole -- neutron star systems. We demonstrate that these systems can feature copious electromagnetic flaring activity, which we find depends on the magnetic field orientation but not on black hole spin. Due to interactions with the surrounding magnetosphere, these flares could lead to Fast Radio Burst-like transients and X-ray emission, with $\mathcal{L}_{\rm EM} \lesssim 10^{41} \left( B_\ast/ 10^{12}\, \rm G \right)^2\, \rm erg/ s$ as an upper bound for the luminosity, where $B_\ast$ is the magnetic field strength on the surface of the neutron star. △ Less

Submitted 28 September, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 13 pages, 6 figures; accepted version, to appear in ApJL

Journal ref: ApJL 956 L33 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.

Comments: 25 pages, 18 figures, 1 table; Accepted for publication in ApJ

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.

Comments: Published in MNRAS, 16 pages, 9 Figures

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.

Comments: 21 pages, 11 figures, Accepted for publication by ApJL

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

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

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

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

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

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

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.

Comments: Accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 132, 085202 (2024)

Millimeter Observational Signatures of Flares in Magnetically Arrested Black Hole Accretion Models

Authors: He Jia, Bart Ripperda, Eliot Quataert, Christopher J. White, Koushik Chatterjee, Alexander Philippov, Matthew Liska

Abstract: In general relativistic magneto-hydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes. The change in 230 GHz flux during the expe… ▽ More In general relativistic magneto-hydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating $γ\gtrsim 100$ ($T_e \gtrsim 10^{11}$ K) electrons. For models in which the electrons are heated to $T_e \sim 10^{11}$ K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multi-wavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization. △ Less

Submitted 23 September, 2023; v1 submitted 21 January, 2023; originally announced January 2023.

Comments: 11+7 pages, 9+7 figures, 1 table, accepted by MNRAS

Collisionless accretion onto black holes: dynamics and flares

Authors: Alisa Galishnikova, Alexander Philippov, Eliot Quataert, Fabio Bacchini, Kyle Parfrey, Bart Ripperda

Abstract: We study the accretion of collisionless plasma onto a rotating black hole from first principles using axisymmetric general-relativistic particle-in-cell simulations. We carry out a side-by-side comparison of these results to analogous general-relativistic magnetohydrodynamic simulations. Although there are many similarities in the overall flow dynamics, three key differences between the kinetic an… ▽ More We study the accretion of collisionless plasma onto a rotating black hole from first principles using axisymmetric general-relativistic particle-in-cell simulations. We carry out a side-by-side comparison of these results to analogous general-relativistic magnetohydrodynamic simulations. Although there are many similarities in the overall flow dynamics, three key differences between the kinetic and fluid simulations are identified. Magnetic reconnection is more efficient, and rapidly accelerates a nonthermal particle population, in our kinetic approach. In addition, the plasma in the kinetic simulations develops significant departures from thermal equilibrium, including pressure anisotropy that excites kinetic-scale instabilities, and a large field-aligned heat flux near the horizon that approaches the free-streaming value. We discuss the implications of our results for modeling event-horizon scale observations of Sgr A* and M87 by GRAVITY and the Event Horizon Telescope. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 10 pages, 8 figures, submitted

Magnetically driven coupling in relativistic radiation-mediated shocks

Authors: J. F. Mahlmann, A. Vanthieghem, A. A. Philippov, A. Levinson, E. Nakar, F. Fiuza

Abstract: The radiation drag in photon-rich environments of cosmic explosions can seed kinetic instabilities by inducing velocity spreads between relativistically streaming plasma components. Such microturbulence is likely imprinted on the breakout signals of radiation-mediated shocks. However, large-scale, transverse magnetic fields in the deceleration region of the shock transition can suppress the domina… ▽ More The radiation drag in photon-rich environments of cosmic explosions can seed kinetic instabilities by inducing velocity spreads between relativistically streaming plasma components. Such microturbulence is likely imprinted on the breakout signals of radiation-mediated shocks. However, large-scale, transverse magnetic fields in the deceleration region of the shock transition can suppress the dominant kinetic instabilities by preventing the development of velocity separations between electron-positron pairs and a heavy ion species. We use a one-dimensional (1D) five-fluid radiative transfer code to generate self-consistent profiles of the radiation drag force and plasma composition in the deceleration region. For increasing magnetization, our models predict rapidly growing pair multiplicities and a substantial radiative drag developing self-similarly throughout the deceleration region. We extract the critical magnetization parameter $σ_{c}$, determining the limiting magnetic field strength at which a three-species plasma can develop kinetic instabilities before reaching the isotropized downstream. For a relativistic, single ion plasma drifting with $γ_{u} = 10$ in the upstream of a relativistic radiation-mediated shock, we find the threshold $σ_{c}\approx 10^{-7}$ for the onset of microturbulence. Suppression of plasma instabilities in the case of multi-ion composition would likely require much higher values of $σ_{c}$. Identifying high-energy signatures of microturbulence in shock-breakout signals and combining them with the magnetization limits provided in this work will allow a deeper understanding of the magnetic environment of cosmic explosions like supernovae, gamma-ray bursts, and neutron star binary mergers. △ Less

Submitted 24 January, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 12 pages, 11 figures, Accepted for publication by MNRAS

Plasmoid Instability in the Multiphase Interstellar Medium

Authors: Drummond B. Fielding, Bart Ripperda, Alexander A. Philippov

Abstract: The processes controlling the complex clump structure, phase distribution, and magnetic field geometry that develops across a broad range of scales in the turbulent interstellar medium remains unclear. Using unprecedentedly high resolution three-dimensional magnetohydrodynamic simulations of thermally unstable turbulent systems, we show that large current sheets unstable to plasmoid-mediated recon… ▽ More The processes controlling the complex clump structure, phase distribution, and magnetic field geometry that develops across a broad range of scales in the turbulent interstellar medium remains unclear. Using unprecedentedly high resolution three-dimensional magnetohydrodynamic simulations of thermally unstable turbulent systems, we show that large current sheets unstable to plasmoid-mediated reconnection form regularly throughout the volume. The plasmoids form in three distinct environments: (i) within cold clumps, (ii) at the asymmetric interface of the cold and warm phases, and (iii) within the warm, volume-filling phase. We then show that the complex magneto-thermal phase structure is characterized by a predominantly highly magnetized cold phase, but that regions of high magnetic curvature, which are the sites of reconnection, span a broad range in temperature. Furthermore, we show that thermal instabilities change the scale dependent anisotropy of the turbulent magnetic field, reducing the increase in eddy elongation at smaller scales. Finally, we show that most of the mass is contained in one contiguous cold structure surrounded by smaller clumps that follow a scale free mass distribution. These clumps tend to be highly elongated and exhibit a size versus internal velocity relation consistent with supersonic turbulence, and a relative clump distance-velocity scaling consistent with subsonic motion. We discuss the striking similarity of cold plasmoids to observed tiny scale atomic and ionized structures and HI fibers, and consider how the prevalence of plasmoids will modify the motion of charged particles thereby impacting cosmic ray transport and thermal conduction in the ISM and other similar systems. △ Less

Submitted 29 April, 2024; v1 submitted 11 November, 2022; originally announced November 2022.

Comments: 21 pages, 13 figures. For associated movies, see https://dfielding14.github.io/movies/

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.

Comments: 22 pages, 12 figures

Journal ref: ApJ 943 (2023) 105

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.

Comments: 16 pages, 5 figures, 1 table

Journal ref: ApJL 943 (2023) L29

Reconnection-powered fast radio transients from coalescing neutron star binaries

Authors: Elias R. Most, Alexander A. Philippov

Abstract: It is an open question whether and how gravitational wave events involving neutron stars can be preceded by electromagnetic counterparts. This work shows that the collision of two neutron stars with magnetic fields well below magnetar-level strengths can produce millisecond Fast-Radio-Burst-like transients. Using global force-free electrodynamics simulations, we demonstrate that electromagnetic fl… ▽ More It is an open question whether and how gravitational wave events involving neutron stars can be preceded by electromagnetic counterparts. This work shows that the collision of two neutron stars with magnetic fields well below magnetar-level strengths can produce millisecond Fast-Radio-Burst-like transients. Using global force-free electrodynamics simulations, we demonstrate that electromagnetic flares, produced by overtwisted common flux tubes in the binary magnetosphere, collide with the orbital current sheet and compress it, resulting in enhanced magnetic reconnection. As a result, the current sheet fragments into a sequence of plasmoids, which collide with each other leading to the emission of coherent electromagnetic waves. The resulting millisecond-long burst of radiation should have frequencies in the range of $10-20\,\rm GHz$ for magnetic fields of $B^{\ast}=10^{11}\, \rm G$ at the stellar surfaces. △ Less

Submitted 17 April, 2023; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 10 pages, 4 figures, version accepted by PRL

Electromagnetic precursor flares from the late inspiral of neutron star binaries

Authors: Elias R. Most, Alexander A. Philippov

Abstract: The coalescence of two neutron stars is accompanied by the emission of gravitational waves, and can also feature electromagnetic counterparts powered by mass ejecta and the formation of a relativistic jet after the merger. Since neutron stars can feature strong magnetic fields, the non-trivial interaction of the neutron star magnetospheres might fuel potentially powerful electromagnetic transients… ▽ More The coalescence of two neutron stars is accompanied by the emission of gravitational waves, and can also feature electromagnetic counterparts powered by mass ejecta and the formation of a relativistic jet after the merger. Since neutron stars can feature strong magnetic fields, the non-trivial interaction of the neutron star magnetospheres might fuel potentially powerful electromagnetic transients prior to merger. A key process powering those precursor transients is relativistic reconnection in strong current sheets formed between the two stars. In this work, we provide a detailed analysis of how the twisting of the common magnetosphere of the binary leads to an emission of electromagnetic flares, akin to those produced in the solar corona. By means of relativistic force-free electrodynamics simulations, we clarify the role of different magnetic field topologies in the process. We conclude that flaring will always occur for suitable magnetic field alignments, unless one of the neutron stars has a magnetic field significantly weaker than the other. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Comments: 19 pages, 10 figures

Magnetar bursts due to Alfvén wave nonlinear breakout

Authors: Yajie Yuan, Andrei M. Beloborodov, Alexander Y. Chen, Yuri Levin, Elias R. Most, Alexander A. Philippov

Abstract: The most common form of magnetar activity is short X-ray bursts, with durations from milliseconds to seconds, and luminosities ranging from $10^{36}$ to $10^{43}\ {\rm erg}\,{\rm s}^{-1}$. Recently, an X-ray burst from the galactic magnetar SGR 1935+2154 was detected to be coincident with two fast radio burst (FRB) like events from the same source, providing evidence that FRBs may be linked to mag… ▽ More The most common form of magnetar activity is short X-ray bursts, with durations from milliseconds to seconds, and luminosities ranging from $10^{36}$ to $10^{43}\ {\rm erg}\,{\rm s}^{-1}$. Recently, an X-ray burst from the galactic magnetar SGR 1935+2154 was detected to be coincident with two fast radio burst (FRB) like events from the same source, providing evidence that FRBs may be linked to magnetar bursts. Using fully 3D force-free electrodynamics simulations, we show that such magnetar bursts may be produced by Alfvén waves launched from localized magnetar quakes: a wave packet propagates to the outer magnetosphere, becomes nonlinear, and escapes the magnetosphere, forming an ultra-relativistic ejecta. The ejecta pushes open the magnetospheric field lines, creating current sheets behind it. Magnetic reconnection can happen at these current sheets, leading to plasma energization and X-ray emission. The angular size of the ejecta can be compact, $\lesssim 1$ sr if the quake launching region is small, $\lesssim 0.01$ sr at the stellar surface. We discuss implications for the FRBs and the coincident X-ray burst from SGR 1935+2154. △ Less

Submitted 1 June, 2022; v1 submitted 18 April, 2022; originally announced April 2022.

Comments: 14 pages, 11 figures, accepted for publication in ApJ

Journal ref: ApJ (2022) 933, 174

Electromagnetic fireworks: Fast radio bursts from rapid reconnection in the compressed magnetar wind

Authors: J. F. Mahlmann, A. A. Philippov, A. Levinson, A. Spitkovsky, H. Hakobyan

Abstract: One scenario for the generation of fast radio bursts (FRBs) is magnetic reconnection in a current sheet of the magnetar wind. Compressed by a strong magnetic pulse induced by a magnetar flare, the current sheet fragments into a self-similar chain of magnetic islands. Time-dependent plasma currents at their interfaces produce coherent radiation during their hierarchical coalescence. We investigate… ▽ More One scenario for the generation of fast radio bursts (FRBs) is magnetic reconnection in a current sheet of the magnetar wind. Compressed by a strong magnetic pulse induced by a magnetar flare, the current sheet fragments into a self-similar chain of magnetic islands. Time-dependent plasma currents at their interfaces produce coherent radiation during their hierarchical coalescence. We investigate this scenario using 2D radiative relativistic particle-in-cell simulations to compute the efficiency of the coherent emission and to obtain frequency scalings. Consistent with expectations, a fraction of the reconnected magnetic field energy, $f\sim 0.002$, is converted to packets of high-frequency fast magnetosonic waves which can escape from the magnetar wind as radio emission. In agreement with analytical estimates, we find that magnetic pulses of $10^{47}\text{erg}\;\text{s}^{-1}$ can trigger relatively narrowband GHz emission with luminosities of approximately $10^{42}\text{erg}\;\text{s}^{-1}$, sufficient to explain bright extragalactic FRBs. The mechanism provides a natural explanation for a downward frequency drift of burst signals, as well as the $\sim 100\;\text{ns}$ substructure recently detected in FRB 20200120E. △ Less

Submitted 17 June, 2022; v1 submitted 8 March, 2022; originally announced March 2022.

Comments: 19 pages, 9 figures, Accepted for publication by ApJL

Journal ref: ApJL 932 L20 (2022)

Images of magnetospheric reconnection-powered radiation around supermassive black holes

Authors: Benjamin Crinquand, Benoît Cerutti, Guillaume Dubus, Kyle Parfrey, Alexander A. Philippov

Abstract: Accreting supermassive black holes can now be observed at the event-horizon scale at mm wavelengths. Current predictions for the image rely on hypotheses (fluid modeling, thermal electrons) which might not always hold in the vicinity of the black hole, so that a full kinetic treatment is in order. In this letter, we describe the first 3D global general-relativistic particle-in-cell simulation of a… ▽ More Accreting supermassive black holes can now be observed at the event-horizon scale at mm wavelengths. Current predictions for the image rely on hypotheses (fluid modeling, thermal electrons) which might not always hold in the vicinity of the black hole, so that a full kinetic treatment is in order. In this letter, we describe the first 3D global general-relativistic particle-in-cell simulation of a black-hole magnetosphere. The system displays a persistent equatorial current sheet. Synthetic images are computed by ray-tracing synchrotron emission from nonthermal particles accelerated in this current sheet by magnetic reconnection. We identify several time-dependent features of the image at moderate viewing angles: a variable radius of the ring, and hot spots moving along it. In this regime, our model predicts that most of the flux of the image lies inside the critical curve. These results could help understand future observations of black-hole magnetospheres at improved temporal and spatial resolution. △ Less

Submitted 28 October, 2022; v1 submitted 9 February, 2022; originally announced February 2022.

Comments: 6 pages + Supplemental Material, 3 figures, accepted for publication in PRL

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.

Comments: Updated to accepted version

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.

Comments: 12 pages, 9 figures

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.

Comments: accepted for publication in the ApJ Letters

Modeling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

Authors: Elias R. Most, Jorge Noronha, Alexander A. Philippov

Abstract: Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modeling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially… ▽ More Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modeling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially true when the mean free path in the plasma is large compared to the system size, and kinetic effects start to become important. Going beyond common approaches used in the literature, we describe a fully relativistic covariant 14-moment based two-fluid system appropriate for the study of electron-ion or electron-positron plasmas. This generalized Israel-Stewart-like system of equations of motion is obtained directly from the relativistic Boltzmann-Vlasov equation. Crucially, this new formulation can account for non-ideal effects, such as anisotropic pressures and heat fluxes. We show that a relativistic two-fluid plasma can be recast as a single fluid coupled to electromagnetic fields with (potentially large) out-of-equilibrium corrections. In particular, we keep all electron degrees of freedom, which provide self-consistent evolution equations for electron temperature and momentum. The equations outlined in this paper are able to capture the full two-fluid character of collisionless plasmas found in black hole accretion and flaring processes around compact objects, as well Braginskii-like two-fluid magnetohydrodynamics applicable to weakly collisional plasmas inside accretion disks. This new formulation will be instrumental in the construction of a large class of next-generation simulations of relativistic transient phenomena produced around black holes and neutron stars. △ Less

Submitted 2 July, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: 15 pages, version accepted by MNRAS

Boundary Layers of Accretion Disks: Wave-Driven Transport and Disk Evolution

Authors: Matthew S. B. Coleman, Roman R. Rafikov, Alexander A. Philippov

Abstract: Astrophysical objects possessing a material surface (white dwarfs, young stars, etc.) may accrete gas from the disc through the so-called surface boundary layer (BL), in which the angular velocity of the accreting gas experiences a sharp drop. Acoustic waves excited by the supersonic shear in the BL play an important role in mediating the angular momentum and mass transport through that region. He… ▽ More Astrophysical objects possessing a material surface (white dwarfs, young stars, etc.) may accrete gas from the disc through the so-called surface boundary layer (BL), in which the angular velocity of the accreting gas experiences a sharp drop. Acoustic waves excited by the supersonic shear in the BL play an important role in mediating the angular momentum and mass transport through that region. Here we examine the characteristics of the angular momentum transport produced by the different types of wave modes emerging in the inner disc, using the results of a large suite of hydrodynamic simulations of the BLs. We provide a comparative analysis of the transport properties of different modes across the range of relevant disc parameters. In particular, we identify the types of modes which are responsible for the mass accretion onto the central object. We find the correlated perturbations of surface density and radial velocity to provide an important contribution to the mass accretion rate. Although the wave-driven transport is intrinsically non-local, we do observe a clear correlation between the angular momentum flux injected into the disc by the waves and the mass accretion rate through the BL. We find the efficiency of angular momentum transport (normalized by thermal pressure) to be a weak function of the flow Mach number. We also quantify the wave-driven evolution of the inner disc, in particular the modification of the angular frequency profile in the disc. Our results pave the way for understanding wave-mediated transport in future three-dimensional, magnetohydrodynamic studies of the BLs. △ Less

Submitted 10 March, 2022; v1 submitted 4 November, 2021; originally announced November 2021.

Comments: 16 pages, 9 figures, submitted to MNRAS

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.

Comments: ApJ Letters (2022), 924, L32

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.

Comments: 5 pages, 4 figures, accepted PRL

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.

Journal ref: The Astrophysical Journal (2021), 921(1), 75

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.

Comments: Accepted

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.

Comments: Accepted in JPP, accompanied by Paper I: Weak Alfvénic turbulence in relativistic plasmas I: dynamical equations and basic dynamics of interacting resonant triads

Boundary Layers of Accretion Disks: Discovery of Vortex-Driven Modes and Other Waves

Authors: Matthew S. B. Coleman, Roman R. Rafikov, Alexander A. Philippov

Abstract: Disk accretion onto weakly magnetized objects possessing a material surface must proceed via the so-called boundary layer (BL) - a region at the inner edge of the disk, in which the velocity of accreting material abruptly decreases from its Keplerian value. Supersonic shear arising in the BL is known to be conducive to excitation of acoustic waves that propagate into both the accretor and the disk… ▽ More Disk accretion onto weakly magnetized objects possessing a material surface must proceed via the so-called boundary layer (BL) - a region at the inner edge of the disk, in which the velocity of accreting material abruptly decreases from its Keplerian value. Supersonic shear arising in the BL is known to be conducive to excitation of acoustic waves that propagate into both the accretor and the disk, enabling angular momentum and mass transport across the BL. We carry out a numerical exploration of different wave modes that operate near the BL, focusing on their morphological characteristics in the innermost parts of accretion disk. Using a large suite of simulations covering a broad range of Mach numbers (of the supersonic shear flow in the BL), we provide accurate characterization of the different types of modes, verifying their properties against analytical results, when available. We discover new types of modes, in particular, global spiral density waves launched by vortices forming in the disk near the BL as a result of the Rossby wave instability; this instability is triggered by the vortensity production in that region caused by the nonlinear damping of acoustic waves. Azimuthal wavenumbers of the dominant modes that we observe appear to increase monotonically with the Mach number of the runs, but a particular mix of modes found in a simulation is mildly stochastic. Our results provide a basis for better understanding of the angular momentum and mass transport across the BL as well as the emission variability in accreting objects. △ Less

Submitted 11 October, 2021; v1 submitted 22 March, 2021; originally announced March 2021.

Comments: 25 pages, 16 figures, changed title and reorganized structure, replaced to match the version accepted by MNRAS

Synthetic gamma-ray lightcurves of Kerr black-hole magnetospheric activity from particle-in-cell simulations

Authors: Benjamin Crinquand, Benoît Cerutti, Guillaume Dubus, Kyle Parfrey, Alexander Philippov

Abstract: Context: The origin of ultra-rapid flares of very high-energy radiation from active galactic nuclei remains elusive. Magnetospheric processes, occurring in the close vicinity of the central black hole, could account for these flares. Aims: We aim to bridge the gap between simulations and observations by synthesizing gamma-ray lightcurves in order to characterize the activity of a black-hole magn… ▽ More Context: The origin of ultra-rapid flares of very high-energy radiation from active galactic nuclei remains elusive. Magnetospheric processes, occurring in the close vicinity of the central black hole, could account for these flares. Aims: We aim to bridge the gap between simulations and observations by synthesizing gamma-ray lightcurves in order to characterize the activity of a black-hole magnetosphere, using kinetic simulations. Methods: We perform global axisymmetric two-dimensional general-relativistic particle-in-cell simulations of a Kerr black-hole magnetosphere. We include a self-consistent treatment of radiative processes and plasma supply, as well as a realistic magnetic configuration, with a large-scale equatorial current sheet. We couple our particle-in-cell code with a ray-tracing algorithm, in order to produce synthetic lightcurves. Results: These simulations show a highly dynamic magnetosphere, as well as very efficient dissipation of the magnetic energy. An external supply of magnetic flux is found to maintain the magnetosphere in a dynamic state, otherwise the magnetosphere settles in a quasi-steady Wald-like configuration. The dissipated energy is mostly converted to gamma-ray photons. The lightcurves at low viewing angle (face-on) mainly trace the spark gap activity and exhibit high variability. On the other hand, no significant variability is found at high viewing angle (edge-on), where the main contribution comes from the reconnecting current sheet. Conclusions: We observe that black-hole magnetospheres with a current sheet are characterized by a very high radiative efficiency. The typical amplitude of the flares in our simulations is lower than what is detected in active galactic nuclei. Such flares could result from the variation of parameters external to the black hole △ Less

Submitted 20 May, 2021; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 12 pages, 10 figures. Accepted in Astronomy and Astrophysics

Journal ref: A&A 650, A163 (2021)

Simulations of the radio polarization of a precessing pulsar PSR J1906+0746

Authors: A. K. Galishnikova, A. A. Philippov, V. S. Beskin

Abstract: The recently constructed theory of radio wave propagation in the pulsar magnetosphere outlines the general aspects of the radio light curve and polarization formation. It allows us to describe general properties of mean profiles, such as the position angle ($PA$) of the linear polarization, and the circular polarization for the realistic structure of the pair creation region in the pulsar magnetos… ▽ More The recently constructed theory of radio wave propagation in the pulsar magnetosphere outlines the general aspects of the radio light curve and polarization formation. It allows us to describe general properties of mean profiles, such as the position angle ($PA$) of the linear polarization, and the circular polarization for the realistic structure of the pair creation region in the pulsar magnetosphere. In this work, we present an application of the radio wave propagation theory to the radio observations of pulsar PSR J1906+0746. This pulsar is particularly interesting because observations of relativistic spin-precession in a binary system allows us to put strong constraints on its geometry. Because it is an almost orthogonal rotator, the pulsar allows us to observe both magnetic poles; as we show, this is crucial for testing the theory of radio wave propagation and obtaining constraints on the parameters of magnetospheric plasma. Our results show that plasma parameters are qualitatively consistent with theories of pair plasma production in polar cap discharges. Specifically, for PSR J1906+0746, we constrain the plasma multiplicity $λ\sim 10^3$ and the Lorentz-factor of secondary plasma $γ\sim $ a few hundred. △ Less

Submitted 1 September, 2020; v1 submitted 31 August, 2020; originally announced August 2020.

Comments: 9 pages, 4 figures. Published in MNRAS

Journal ref: July 2020 Monthly Notices of the Royal Astronomical Society, Volume 497, Issue 3, pp.2831-2838

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.

Comments: 11 pages, 10 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 642, A204 (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.

Alfvén wave mode conversion in pulsar magnetospheres

Authors: Yajie Yuan, Yuri Levin, Ashley Bransgrove, Alexander Philippov

Abstract: The radio emission anomaly coincident with the 2016 glitch of the Vela pulsar may be caused by a star quake that launches Alfvén waves into the magnetosphere, disturbing the original radio emitting region. To quantify the lifetime of the Alfvén waves, we investigate a possible energy loss mechanism, the conversion of Alfvénwaves into fast magnetosonic waves. Using axisymmetric force-free simulatio… ▽ More The radio emission anomaly coincident with the 2016 glitch of the Vela pulsar may be caused by a star quake that launches Alfvén waves into the magnetosphere, disturbing the original radio emitting region. To quantify the lifetime of the Alfvén waves, we investigate a possible energy loss mechanism, the conversion of Alfvénwaves into fast magnetosonic waves. Using axisymmetric force-free simulations, we follow the propagation of Alfvén waves launched from the stellar surface with small amplitude into the closed zone of a force-free dipolar pulsar magnetosphere. We observe mode conversion happening in the ideal force-free regime. The conversion efficiency during the first passage of the Alfvén wave through the equator can be large, for waves that reach large amplitudes as they travel away from the star, or propagate on the field lines passing close to the Y-point. However, the conversion efficiency is reduced due to dephasing on subsequent passages and considerable Alfvén power on the closed field lines remains. Thus while some leakage into the fast mode happens, we need detailed understanding of the original quenching in order to say whether mode conversion alone can lead to reactivation of the pulsar on a short timescale. △ Less

Submitted 14 December, 2020; v1 submitted 22 July, 2020; originally announced July 2020.

Comments: 14 pages, 14 figures, accepted for publication in ApJ

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.

Comments: Added Sgr A* flare analysis, flare energetics and statistics, retrograde and Schwarzschild cases, new figure on convergence of adaptive mesh refinement

Journal ref: Astrophysical Journal, 900, 2 (2020)

Multi-dimensional simulations of ergospheric pair discharges around black holes

Authors: Benjamin Crinquand, Benoît Cerutti, Alexander Philippov, Kyle Parfrey, Guillaume Dubus

Abstract: Black holes are known to launch powerful relativistic jets and emit highly variable gamma radiation. How these jets are loaded with plasma remains poorly understood. Spark gaps are thought to drive particle acceleration and pair creation in the black-hole magnetosphere. In this paper, we perform 2D axisymmetric general-relativistic particle-in-cell simulations of a monopole black-hole magnetospher… ▽ More Black holes are known to launch powerful relativistic jets and emit highly variable gamma radiation. How these jets are loaded with plasma remains poorly understood. Spark gaps are thought to drive particle acceleration and pair creation in the black-hole magnetosphere. In this paper, we perform 2D axisymmetric general-relativistic particle-in-cell simulations of a monopole black-hole magnetosphere with a realistic treatment of inverse Compton scattering and pair production. We find that the magnetosphere can self-consistently fill itself with plasma and activate the Blandford-Znajek mechanism. A highly time-dependent spark gap opens near the inner light surface which injects pair plasma into the magnetosphere. These results may account for the high-energy activity observed in active galactic nuclei and explain the origin of plasma at the base of the jet. △ Less

Submitted 7 March, 2020; originally announced March 2020.

Comments: 5 pages, 3 figures + 5 pages, 2 figures in the Supplemental Material. Accepted in Physical Review Letters

Electromagnetic precursors to gravitational wave events: Numerical simulations of flaring in pre-merger binary neutron star magnetospheres

Authors: Elias R. Most, Alexander A. Philippov

Abstract: The detection of gravitational waves from neutron star merger events has opened up a new field of multi-messenger astronomy linking gravitational waves events to short-gamma ray bursts and kilonova afterglows. A further - yet to be discovered - electromagnetic counterpart is precursor emission produced by the non-trivial interaction of the magnetospheres of the two neutron stars prior to merger. B… ▽ More The detection of gravitational waves from neutron star merger events has opened up a new field of multi-messenger astronomy linking gravitational waves events to short-gamma ray bursts and kilonova afterglows. A further - yet to be discovered - electromagnetic counterpart is precursor emission produced by the non-trivial interaction of the magnetospheres of the two neutron stars prior to merger. By performing special-relativistic force-free simulations of orbiting neutron stars we discuss the effect of different magnetic field orientations and show how the emission can be significantly enhanced by differential motion present in the binary, either due to stellar spins or misaligned stellar magnetospheres. We find that the built-up of twist in the magnetic flux tube connecting the two stars can lead to the repeated emission of powerful flares for a variety of orbital configurations. We also discuss potential coherent radio emission mechanisms in the flaring process. △ Less

Submitted 14 April, 2020; v1 submitted 16 January, 2020; originally announced January 2020.

Comments: 9 pages, 5 figures, matches version accepted by ApJL, corrected typo in appendix

Journal ref: ApJL 893 L6, 2020

Origin of Pulsar Radio Emission

Authors: Alexander Philippov, Andrey Timokhin, Anatoly Spitkovsky

Abstract: Since pulsars were discovered as emitters of bright coherent radio emission more than half a century ago, the cause of the emission has remained a mystery. In this Letter we demonstrate that coherent radiation can be directly generated in non-stationary pair plasma discharges which are responsible for filling the pulsar magnetosphere with plasma. By means of large-scale two-dimensional kinetic pla… ▽ More Since pulsars were discovered as emitters of bright coherent radio emission more than half a century ago, the cause of the emission has remained a mystery. In this Letter we demonstrate that coherent radiation can be directly generated in non-stationary pair plasma discharges which are responsible for filling the pulsar magnetosphere with plasma. By means of large-scale two-dimensional kinetic plasma simulations, we show that if pair creation is non-uniform across magnetic field lines, the screening of electric field by freshly produced pair plasma is accompanied by the emission of waves which are electromagnetic in nature. Using localized simulations of the screening process, we identify these waves as superluminal ordinary (O) modes, which should freely escape from the magnetosphere as the plasma density drops along the wave path. The spectrum of the waves is broadband and the frequency range is comparable to that of observed pulsar radio emission. △ Less

Submitted 25 June, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: 9 pages, 4 figures

Journal ref: Phys. Rev. Lett. 124, 245101 (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.

Comments: 9 pages, 5 figures, accepted for publication in ApJL

Kink instability: evolution and energy dissipation in Relativistic Force-Free Non-Rotating Jets

Authors: Omer Bromberg, Chandra B. Singh, Jordy Davelaar, Alexander A. Philippov

Abstract: We study the evolution of kink instability in a force-free, non-rotating plasma column of high magnetization. The main dissipation mechanism is identified as reconnection of magnetic field-lines with various intersection angles, driven by the compression of the growing kink lobes. We measure dissipation rates ${\rm d} U_{Bφ}/{{\rm d}t} \approx -0.1 U_{Bφ}/τ$, where $τ$ is the linear growth time of… ▽ More We study the evolution of kink instability in a force-free, non-rotating plasma column of high magnetization. The main dissipation mechanism is identified as reconnection of magnetic field-lines with various intersection angles, driven by the compression of the growing kink lobes. We measure dissipation rates ${\rm d} U_{Bφ}/{{\rm d}t} \approx -0.1 U_{Bφ}/τ$, where $τ$ is the linear growth time of the kink instability. This value is consistent with the expansion velocity of the kink mode, which drives the reconnection. The relaxed state is close to a force-free Taylor state. We constraint the energy of that state using considerations from linear stability analysis. Our results are important for understanding magnetic field dissipation in various extreme astrophysical objects, most notably in relativistic jets. We outline the evolution of the kink instability in such jets and derive constrains on the conditions that allow for the kink instability to grow in these systems. △ Less

Submitted 27 August, 2019; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: 18 pages, 13 figures

Extreme Plasma Astrophysics

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

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

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

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

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.

Comments: 7 pages, 4 figures

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.

Comments: 6 pages, 4 figures. Submitted to PRL

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

Effects of synchrotron cooling and pair production on collisionless relativistic reconnection

Authors: Hayk Hakobyan, Alexander Philippov, Anatoly Spitkovsky

Abstract: High energy radiation from nonthermal particles accelerated in relativistic magnetic reconnection is thought to be important in many astrophysical systems, ranging from blazar jets and black hole accretion disk coronae to pulsars and magnetar flares. The presence of a substantial density of high energy photons ($>$MeV) in these systems can make two-photon pair production ($γγ\to e^-e^+$) an additi… ▽ More High energy radiation from nonthermal particles accelerated in relativistic magnetic reconnection is thought to be important in many astrophysical systems, ranging from blazar jets and black hole accretion disk coronae to pulsars and magnetar flares. The presence of a substantial density of high energy photons ($>$MeV) in these systems can make two-photon pair production ($γγ\to e^-e^+$) an additional source of plasma particles and can affect the radiative properties of these objects. We present the results of novel particle-in-cell simulations that track both the radiated synchrotron photons and the created pairs, with which we study the evolution of a two-dimensional reconnecting current sheet in pair plasma. Synchrotron radiation from accelerated particles in the current sheet produces hot secondary pairs in the upstream which are later advected into the current sheet where they are reaccelerated and produce more photons. In the optically thin regime, when most of the radiation is leaving the upstream unaffected, this process is self-regulating and depends only on the background magnetic field and the optical depth of photons to pair production. The extra plasma loading also affects the properties of reconnection. We study how the inflow of the secondary plasma, with multiplicities up to several hundred, reduces the effective magnetization of the plasma, suppressing the acceleration and thus decreasing the high energy photon spectrum cutoff. This offers an explanation for the weak dependence of the observed gamma-ray cutoff in pulsars on the magnetic field at the light cylinder. △ Less

Submitted 27 September, 2018; originally announced September 2018.

Comments: 22 pages, 15 figures, submitted to ApJ