Showing 1–50 of 159 results for author: Beloborodov, A

Giant Hall Waves Launched by Superconducting Phase Transition in Pulsars

Authors: Ashley Bransgrove, Yuri Levin, Andrei M. Beloborodov

Abstract: The cores of pulsars are expected to become superconducting soon after birth. The transition to type-II superconductivity is associated with the bunching of magnetic field lines into discrete superconducting flux tubes which possess enormous tension. The coupling of the crust to the flux tubes implies the existence of huge tangential magnetic fields at the crust-core interface. We show that the tr… ▽ More The cores of pulsars are expected to become superconducting soon after birth. The transition to type-II superconductivity is associated with the bunching of magnetic field lines into discrete superconducting flux tubes which possess enormous tension. The coupling of the crust to the flux tubes implies the existence of huge tangential magnetic fields at the crust-core interface. We show that the transition to superconductivity triggers a highly non-linear response in the Hall drift of the crustal magnetic field, an effect which was neglected in previous numerical modelling. We argue that at the time of the phase transition giant Hall waves are launched from the crust-core interface toward the surface. Our models show that if the crust contains a multipolar magnetic field $\sim 10^{13}$ G, the amplitude of the Hall waves is $\sim 10^{15}$ G. The elastic deformation of the lattice is included in our models, which allows us to track the time-dependent shear stresses everywhere in the crust. The simulations indicate that the Hall waves may be strong enough to break the crust, and could cause star quakes which trigger rotation glitches and changes in the radio pulse profile. The Hall waves also couple to slow magnetospheric changes which cause anomalous braking indices. The emission of the giant Hall waves from the crust-core interface facilitates fast flux expulsion from the superconducting core, provided that the flux tubes in the core are themselves sufficiently mobile. For all of the flux tube mobility prescriptions implemented in this work, the core approaches the Meissner state with B=0 at late times. △ Less

Submitted 20 August, 2024; originally announced August 2024.

Comments: 13 pages, 5 figures, submitted for publication

Tayler Instability Revisited

Authors: Valentin A. Skoutnev, Andrei M. Beloborodov

Abstract: Tayler instability of toroidal magnetic fields $B_φ$ is broadly invoked as a trigger for turbulence and angular momentum transport in stars. This paper presents a systematic revision of the linear stability analysis for a rotating, magnetized, and stably stratified star. For plausible configurations of $B_φ$, instability requires diffusive processes: viscosity, magnetic diffusivity, or thermal/com… ▽ More Tayler instability of toroidal magnetic fields $B_φ$ is broadly invoked as a trigger for turbulence and angular momentum transport in stars. This paper presents a systematic revision of the linear stability analysis for a rotating, magnetized, and stably stratified star. For plausible configurations of $B_φ$, instability requires diffusive processes: viscosity, magnetic diffusivity, or thermal/compositional diffusion. Our results reveal a new physical picture, demonstrating how different diffusive effects independently trigger instability of two types of waves in the rotating star: magnetostrophic waves and inertial waves. It develops via overstability of the waves, whose growth rate sharply peaks at some characteristic wavenumbers. We determine instability conditions for each wave branch and find the characteristic wavenumbers. The results are qualitatively different for stars with magnetic Prandtl number $Pm\ll 1$ (e.g. the Sun) and $Pm\gg 1$ (e.g. protoneutron stars). The parameter dependence of unstable modes suggests a non-universal scaling of the possible Tayler-Spruit dynamo. △ Less

Submitted 5 September, 2024; v1 submitted 29 April, 2024; originally announced April 2024.

Comments: Accepted in ApJ

Monster Shocks, Gamma-Ray Bursts, and Black Hole Quasi-normal Modes from Neutron-star Collapse

Authors: Elias R. Most, Andrei M. Beloborodov, Bart Ripperda

Abstract: We perform the first magnetohydrodynamic simulation tracking the magnetosphere of a collapsing magnetar. The collapse is expected for massive rotating magnetars formed in merger events, and may occur many hours after the merger. Our simulation suggests a novel mechanism for a gamma-ray burst (GRB) which is uncollimated and forms a delayed high-energy counterpart of the merger gravitational waves.… ▽ More We perform the first magnetohydrodynamic simulation tracking the magnetosphere of a collapsing magnetar. The collapse is expected for massive rotating magnetars formed in merger events, and may occur many hours after the merger. Our simulation suggests a novel mechanism for a gamma-ray burst (GRB) which is uncollimated and forms a delayed high-energy counterpart of the merger gravitational waves. The simulation shows that the collapse launches an outgoing magnetospheric shock, and a hot magnetized outflow forms behind the shock. The outflow is baryon-free and uncollimated, and its power peaks on a millisecond timescale. Then, the outflow becomes modulated by the ring-down of the nascent black hole, imprinting its kilohertz quasi-normal modes on the GRB tail. △ Less

Submitted 9 October, 2024; v1 submitted 1 April, 2024; originally announced April 2024.

Comments: 9 pages, 4 figures; accepted version; to appear in ApJL

Journal ref: Astrophysical Journal Letters 974 L12 (2024)

The High Energy X-ray Probe (HEX-P): Magnetars and Other Isolated Neutron Stars

Authors: J. A. J. Alford, G. A. Younes, Z. Wadiasingh, M. Abdelmaguid, H. An, M. Bachetti, M. Baring, A. Beloborodov, A. Y. Chen, T. Enoto, J. A. García, J. D. Gelfand, E. V. Gotthelf, A. Harding, C. -P. Hu, A. D. Jaodand, V. Kaspi, C. Kim, C. Kouveliotou, L. Kuiper, K. Mori, M. Nynka, J. Park, D. Stern, J. Valverde , et al. (1 additional authors not shown)

Abstract: The hard X-ray emission from magnetars and other isolated neutron stars remains under-explored. An instrument with higher sensitivity to hard X-rays is critical to understanding the physics of neutron star magnetospheres and also the relationship between magnetars and Fast Radio Bursts (FRBs). High sensitivity to hard X-rays is required to determine the number of magnetars with hard X-ray tails, a… ▽ More The hard X-ray emission from magnetars and other isolated neutron stars remains under-explored. An instrument with higher sensitivity to hard X-rays is critical to understanding the physics of neutron star magnetospheres and also the relationship between magnetars and Fast Radio Bursts (FRBs). High sensitivity to hard X-rays is required to determine the number of magnetars with hard X-ray tails, and to track transient non-thermal emission from these sources for years post-outburst. This sensitivity would also enable previously impossible studies of the faint non-thermal emission from middle-aged rotation-powered pulsars (RPPs), and detailed phase-resolved spectroscopic studies of younger, bright RPPs. The High Energy X-ray Probe (HEX-P) is a probe-class mission concept that will combine high spatial resolution X-ray imaging ($<5$ arcsec half-power diameter (HPD) at 0.2--25 keV) and broad spectral coverage (0.2--80 keV) with a sensitivity superior to current facilities (including XMM-Newton and NuSTAR). HEX-P has the required timing resolution to perform follow-up observations of sources identified by other facilities and positively identify candidate pulsating neutron stars. Here we discuss how HEX-P is ideally suited to address important questions about the physics of magnetars and other isolated neutron stars. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 22 pages, 8 figures, submitted to Frontiers

Damping of strong GHz waves near magnetars and the origin of fast radio bursts

Authors: Andrei M. Beloborodov

Abstract: We investigate how a fast radio burst (FRB) emitted near a magnetar would propagate through its surrounding dipole magnetosphere at radii $r=10^7$-$10^9$ cm. First, we show that a GHz burst emitted in the O-mode with luminosity $L\gg 10^{40}$ erg/s is immediately damped for all propagation directions except a narrow cone along the magnetic axis. Then we examine bursts in the X-mode. GHz waves prop… ▽ More We investigate how a fast radio burst (FRB) emitted near a magnetar would propagate through its surrounding dipole magnetosphere at radii $r=10^7$-$10^9$ cm. First, we show that a GHz burst emitted in the O-mode with luminosity $L\gg 10^{40}$ erg/s is immediately damped for all propagation directions except a narrow cone along the magnetic axis. Then we examine bursts in the X-mode. GHz waves propagating near the magnetic equator behave as magnetohydrodynamic (MHD) waves if they have $L\gg 10^{40}$ erg/s. The waves develop plasma shocks in each oscillation and dissipate at $r\sim 3 \times 10^8 L_{42}^{-1/4}$ cm. Waves with lower $L$ or propagation directions closer to the magnetic axis do not obey MHD. Instead, they interact with individual particles and require a kinetic description. The kinetic interaction quickly accelerates particles to Lorentz factors $10^4$-$10^5$ at the expense of the wave energy, which again results in strong damping of the wave. In either propagation regime, MHD or kinetic, the dipole magnetosphere surrounding the FRB source acts as a pillow absorbing the radio burst and re-radiating the absorbed energy in X-rays. These results constrain the origin of observed FRBs. We argue that the observed FRBs avoid damping because they are emitted by relativistic outflows from magnetospheric explosions, so that the GHz waves do not need to propagate through the outer equilibrium magnetosphere surrounding the magnetar. △ Less

Submitted 31 July, 2024; v1 submitted 22 July, 2023; originally announced July 2023.

Comments: 31 pages, 11 figures, accepted to ApJ

Polarized X-rays from windy accretion in Cygnus X-1

Authors: Juri Poutanen, Alexandra Veledina, Andrei M. Beloborodov

Abstract: Recent X-ray polarimetric data on the prototypical black hole X-ray binary Cyg~X-1 from the Imaging X-ray Polarimetry Explorer present tight constraints on accretion geometry in the hard spectral state. Contrary to general expectations of a low, <1% polarization degree (PD), the observed average PD was found to be a factor of 4 higher. Aligned with the jet position angle on the sky, the observed p… ▽ More Recent X-ray polarimetric data on the prototypical black hole X-ray binary Cyg~X-1 from the Imaging X-ray Polarimetry Explorer present tight constraints on accretion geometry in the hard spectral state. Contrary to general expectations of a low, <1% polarization degree (PD), the observed average PD was found to be a factor of 4 higher. Aligned with the jet position angle on the sky, the observed polarization favors geometry of the X-ray emission region stretched normally to the jet in the accretion disk plane. The high PD is, however, difficult to reconcile with the low orbital inclination of the binary $i\approx30$ deg. We suggest that this puzzle can be explained if the emitting plasma is outflowing with a mildly relativistic velocity $\gtrsim0.4\,c$. Our radiative transfer simulations show that Comptonization in the outflowing medium elongated in the plane of the disk radiates X-rays with the degree and direction of polarization consistent with observations at $i\approx30$ deg. △ Less

Submitted 11 May, 2023; v1 submitted 22 February, 2023; originally announced February 2023.

Comments: 11 pages, 6 figures, ApJL, in press

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)

Monster radiative shocks in the perturbed magnetospheres of neutron stars

Authors: Andrei M. Beloborodov

Abstract: Magnetospheres of neutron stars can be perturbed by star quakes, interaction in a binary system, or sudden collapse of the star. The perturbations are typically in the kHz band and excite magnetohydrodynamic waves. We show that compressive magnetospheric waves steepen into monster shocks, possibly the strongest shocks in the universe. The shocks are radiative, i.e. the plasma energy is radiated be… ▽ More Magnetospheres of neutron stars can be perturbed by star quakes, interaction in a binary system, or sudden collapse of the star. The perturbations are typically in the kHz band and excite magnetohydrodynamic waves. We show that compressive magnetospheric waves steepen into monster shocks, possibly the strongest shocks in the universe. The shocks are radiative, i.e. the plasma energy is radiated before it crosses the shock. As the kHz wave with the radiative shock expands through the magnetosphere, it produces a bright X-ray burst. Then, it launches an approximately adiabatic blast wave, which will expand far from the neutron star. These results suggest a new mechanism for X-ray bursts from magnetars and support the connection of magnetar X-ray activity with fast radio bursts. Similar shocks may persistently occur in magnetized neutron-star binaries, generating an X-ray precursor of the binary merger. Powerful radiative shocks are also predicted in the magnetosphere of a neutron star when it collapses into a black hole, producing a bright X-ray transient. △ Less

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

Comments: 27 pages, ApJ, in press

Saturation of the filamentation instability and dispersion measure of Fast Radio Bursts

Authors: Emanuele Sobacchi, Yuri Lyubarsky, Andrei M. Beloborodov, Lorenzo Sironi, Masanori Iwamoto

Abstract: Nonlinear effects are crucial for the propagation of Fast Radio Bursts (FRBs) near the source. We study the filamentation of FRBs in the relativistic winds of magnetars, which are commonly invoked as the most natural FRB progenitors. As a result of filamentation, the particle number density and the radiation intensity develop strong gradients along the direction of the wind magnetic field. A stead… ▽ More Nonlinear effects are crucial for the propagation of Fast Radio Bursts (FRBs) near the source. We study the filamentation of FRBs in the relativistic winds of magnetars, which are commonly invoked as the most natural FRB progenitors. As a result of filamentation, the particle number density and the radiation intensity develop strong gradients along the direction of the wind magnetic field. A steady state is reached when the plasma pressure balances the ponderomotive force. In such a steady state, particles are confined into periodically spaced thin sheets, and electromagnetic waves propagate between them as in a waveguide. We show that: (i) The dispersion relation resembles that in the initial homogeneous plasma, but the effective plasma frequency is determined by the separation of the sheets, not directly by the mean particle density. (ii) The contribution of relativistic magnetar winds to the dispersion measure of FRBs could be several orders of magnitude larger than previously thought. The dispersion measure of the wind depends on the properties of individual bursts (e.g. the luminosity), and therefore can change significantly among different bursts from repeating FRBs. (iii) Induced Compton scattering is suppressed because most of the radiation propagates in near vacuum regions. △ Less

Submitted 12 January, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 9 pages, accepted for publication in ApJL

Radio Emission and Electric Gaps in Pulsar Magnetospheres

Authors: Ashley Bransgrove, Andrei M. Beloborodov, Yuri Levin

Abstract: The origin of pulsar radio emission is one of the old puzzles in theoretical astrophysics. In this Letter we present a global kinetic plasma simulation which shows from first-principles how and where radio emission can be produced in pulsar magnetospheres. We observe the self-consistent formation of electric gaps which periodically ignite electron-positron discharge. The gaps form above the polar-… ▽ More The origin of pulsar radio emission is one of the old puzzles in theoretical astrophysics. In this Letter we present a global kinetic plasma simulation which shows from first-principles how and where radio emission can be produced in pulsar magnetospheres. We observe the self-consistent formation of electric gaps which periodically ignite electron-positron discharge. The gaps form above the polar-cap, and in the bulk return-current. Discharge of the gaps excites electromagnetic modes which share several features with the radio emission of real pulsars. We also observe the excitation of plasma waves and charge bunches by streaming instabilities in the outer magnetosphere. Our numerical experiment demonstrates that global kinetic models can provide deep insight into the emission physics of pulsars, and may help interpret their multi-wavelength observations. △ Less

Submitted 8 December, 2023; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 9 pages, accepted ApJL

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

Comptonization by Reconnection Plasmoids in Black Hole Coronae II: Electron-Ion Plasma

Authors: Navin Sridhar, Lorenzo Sironi, Andrei M. Beloborodov

Abstract: We perform two-dimensional particle-in-cell simulations of magnetic reconnection in electron-ion plasmas subject to strong Compton cooling and calculate the X-ray spectra produced by this process. The simulations are performed for trans-relativistic reconnection with magnetization $1\leq σ\leq 3$ (defined as the ratio of magnetic tension to plasma rest-mass density), which is expected in the coron… ▽ More We perform two-dimensional particle-in-cell simulations of magnetic reconnection in electron-ion plasmas subject to strong Compton cooling and calculate the X-ray spectra produced by this process. The simulations are performed for trans-relativistic reconnection with magnetization $1\leq σ\leq 3$ (defined as the ratio of magnetic tension to plasma rest-mass density), which is expected in the coronae of accretion disks around black holes. We find that magnetic dissipation proceeds with inefficient energy exchange between the heated ions and the Compton-cooled electrons. As a result, most electrons are kept at a low temperature in Compton equilibrium with radiation, and so thermal Comptonization cannot reach photon energies $\sim 100$ keV observed from accreting black holes. Nevertheless, magnetic reconnection efficiently generates $\sim 100$ keV photons because of mildly relativistic bulk motions of the plasmoid chain formed in the reconnection layer. Comptonization by the plasmoid motions dominates the radiative output and controls the peak of the radiation spectrum $E_{\rm pk}$. We find $E_{\rm pk}\sim 40$ keV for $σ=1$ and $E_{\rm pk}\sim100$ keV for $σ=3$. In addition to the X-ray peak around 100 keV, the simulations show a non-thermal MeV tail emitted by a non-thermal electron population generated near X-points of the reconnection layer. The results are consistent with the typical hard state of accreting black holes. In particular, we find that the spectrum of Cygnus~X-1 is well explained by electron-ion reconnection with $σ\sim 3$. △ Less

Submitted 26 September, 2022; v1 submitted 5 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS: 16 pages, 15 figures, 1 table, 8 appendices. Part-I can be accessed at arXiv:2107.00263

Microphysics of Relativistic Collisionless Electron-ion-positron Shocks

Authors: Daniel Groselj, Lorenzo Sironi, Andrei M. Beloborodov

Abstract: We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizations $σ\lesssim 10^{-4}$ and pair-loading factors $Z_\pm \lesssim 10$ are studied, where $Z_\pm$ is the number of loaded electrons and positrons per ion. We find the following. (1) The shock becomes mediated by the ion Larmor… ▽ More We perform particle-in-cell simulations to elucidate the microphysics of relativistic weakly magnetized shocks loaded with electron-positron pairs. Various external magnetizations $σ\lesssim 10^{-4}$ and pair-loading factors $Z_\pm \lesssim 10$ are studied, where $Z_\pm$ is the number of loaded electrons and positrons per ion. We find the following. (1) The shock becomes mediated by the ion Larmor gyration in the mean field when $σ$ exceeds a critical value $σ_{\rm L}$ that decreases with $Z_\pm$. At $σ\lesssimσ_{\rm L}$ the shock is mediated by particle scattering in the self-generated microturbulent fields, the strength and scale of which decrease with $Z_\pm$, leading to lower $σ_{\rm L}$. (2) The energy fraction carried by the post-shock pairs is robustly in the range between 20% and 50% of the upstream ion energy. The mean energy per post-shock electron scales as $\overline{E}_{\rm e}\propto (Z_\pm+1)^{-1}$. (3) Pair loading suppresses nonthermal ion acceleration at magnetizations as low as $σ\approx 5\times 10^{-6}$. The ions then become essentially thermal with mean energy $\overline{E}_{\rm i}$, while electrons form a nonthermal tail, extending from $E\sim (Z_\pm + 1)^{-1}\overline{E}_{\rm i}$ to $\overline{E}_{\rm i}$. When $σ= 0$, particle acceleration is enhanced by the formation of intense magnetic cavities that populate the precursor during the late stages of shock evolution. Here, the maximum energy of the nonthermal ions and electrons keeps growing over the duration of the simulation. Alongside the simulations, we develop theoretical estimates consistent with the numerical results. Our findings have important implications for models of early gamma-ray burst afterglows. △ Less

Submitted 4 February, 2022; originally announced February 2022.

Comments: 20 pages, 16 figures, submitted to ApJ

Journal ref: Astrophys. J. 933, 74 (2022)

Formation of Magnetically Truncated Accretion Disks in 3D Radiation-Transport Two-Temperature GRMHD Simulations

Authors: M. T. P. Liska, G. Musoke, A. Tchekhovskoy, O. Porth, A. M. Beloborodov

Abstract: Multi-wavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nuclei (AGN) transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole. However, the formation, structure and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magn… ▽ More Multi-wavelength observations suggest that the accretion disk in the hard and intermediate states of X-ray binaries (XRBs) and active galactic nuclei (AGN) transitions from a cold, thin disk at large distances into a hot, thick flow close to the black hole. However, the formation, structure and dynamics of such truncated disks are poorly constrained due to the complexity of the thermodynamic, magnetic, and radiative processes involved. We present the first radiation-transport two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of truncated disks radiating at ~35% of the Eddington luminosity with and without large-scale poloidal magnetic flux. We demonstrate that when a geometrically-thin accretion disk is threaded by large-scale net poloidal magnetic flux, it self-consistently transitions at small radii into a two-phase medium of cold gas clumps floating through a hot, magnetically dominated corona. This transition occurs at a well-defined truncation radius determined by the distance out to which the disk is saturated with magnetic flux. The average ion and electron temperatures in the semi-opaque corona reach, respectively, T_i ~ 10^10K and T_e ~ 5 10^8K. The system produces radiation, powerful collimated jets and broader winds at the total energy efficiency exceeding ~90%, the highest ever energy extraction efficiency from a spinning black hole by a radiatively efficient flow in a GRMHD simulation. This is consistent with jetted ejections observed during XRB outbursts. The two-phase medium may naturally lead to broadened iron line emission observed in the hard state. △ Less

Submitted 9 August, 2022; v1 submitted 10 January, 2022; originally announced January 2022.

Comments: 8 pages, 5 figures, YouTube playlist: https://www.youtube.com/playlist?list=PLDO1oeU33GwlkbTgY9gybClXqqv0ky7Ha, Published in Astrophysical Journal Letters

Filamentation of Fast Radio Bursts in magnetar winds

Authors: Emanuele Sobacchi, Yuri Lyubarsky, Andrei Beloborodov, Lorenzo Sironi

Abstract: Magnetars are the most promising progenitors of Fast Radio Bursts (FRBs). Strong radio waves propagating through the magnetar wind are subject to non-linear effects, including modulation/filamentation instabilities. We derive the dispersion relation for modulations of strong waves propagating in magnetically-dominated pair plasmas focusing on dimensionless strength parameters $a_0\lesssim 1$, and… ▽ More Magnetars are the most promising progenitors of Fast Radio Bursts (FRBs). Strong radio waves propagating through the magnetar wind are subject to non-linear effects, including modulation/filamentation instabilities. We derive the dispersion relation for modulations of strong waves propagating in magnetically-dominated pair plasmas focusing on dimensionless strength parameters $a_0\lesssim 1$, and discuss implications for FRBs. As an effect of the instability, the FRB radiation intensity develops sheets perpendicular to the direction of the wind magnetic field. When the FRB front expands outside the radius where the instability ends, the radiation sheets are scattered due to diffraction. The FRB scattering timescale depends on the properties of the magnetar wind. In a cold wind, the typical scattering timescale is $τ_{\rm sc}\sim{\rm\; μs-ms}$ at the frequency $ν\sim 1{\rm\; GHz}$. The scattering timescale increases at low frequencies, with the scaling $τ_{\rm sc}\proptoν^{-2}$. The frequency-dependent broadening of the brightest pulse of FRB 181112 is consistent with this scaling. From the scattering timescale of the pulse, one can estimate that the wind Lorentz factor is larger than a few tens. In a warm wind, the scattering timescale can approach $τ_{\rm sc}\sim{\rm\; ns}$. Then scattering produces a frequency modulation of the observed intensity with a large bandwidth, $Δν\sim 1/τ_{\rm sc}\gtrsim 100{\rm\; MHz}$. Broadband frequency modulations observed in FRBs could be due to scattering in a warm magnetar wind. △ Less

Submitted 25 January, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: Section 3 is modified. Accepted for publication in MNRAS

Heating of Magnetically Dominated Plasma by Alfvén-Wave Turbulence

Authors: J. Nättilä, A. M. Beloborodov

Abstract: Magnetic energy around astrophysical compact objects can strongly dominate over plasma rest mass. Emission observed from these systems may be fed by dissipation of Alfvén wave turbulence, which cascades to small damping scales, energizing the plasma. We use 3D kinetic simulations to investigate this process. When the cascade is excited naturally, by colliding large-scale Alfvén waves, we observe q… ▽ More Magnetic energy around astrophysical compact objects can strongly dominate over plasma rest mass. Emission observed from these systems may be fed by dissipation of Alfvén wave turbulence, which cascades to small damping scales, energizing the plasma. We use 3D kinetic simulations to investigate this process. When the cascade is excited naturally, by colliding large-scale Alfvén waves, we observe quasithermal heating with no nonthermal particle acceleration. We also find that the particles are energized along the magnetic field lines and so are poor producers of synchrotron radiation. At low plasma densities, our simulations show the transition to "charge-starved" cascades, with a distinct damping mechanism. △ Less

Submitted 13 May, 2022; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 6 pages, 3 figures, 1 appendix. Accepted to PRL

Axisymmetric pulsar magnetosphere revisited

Authors: Rui Hu, Andrei M. Beloborodov

Abstract: We present a global kinetic plasma simulation of an axisymmetric pulsar magnetosphere with self-consistent $e^\pm$ pair production. We use the particle-in-cell method and log-spherical coordinates with a grid size $4096\times 4096$. This allows us to achieve a high voltage induced by the pulsar rotation and investigate pair creation in a young pulsar far from the death line. We find the following.… ▽ More We present a global kinetic plasma simulation of an axisymmetric pulsar magnetosphere with self-consistent $e^\pm$ pair production. We use the particle-in-cell method and log-spherical coordinates with a grid size $4096\times 4096$. This allows us to achieve a high voltage induced by the pulsar rotation and investigate pair creation in a young pulsar far from the death line. We find the following. (1) The energy release and $e^\pm$ creation are strongly concentrated in the thin, Y-shaped current sheet, with a peak localized in a small volume at the Y-point. (2) The Y-point is shifted inward from the light cylinder by $\sim 15\%$, and "breathes" with a small amplitude. (3) The dense $e^\pm$ cloud at the Y-point is in ultra-relativistic rotation, which we call super-rotation, because it exceeds co-rotation with the star. The cloud receives angular momentum flowing from the star along the poloidal magnetic lines. (4) Gamma-ray emission peaks at the Y-point and is collimated in the azimuthal direction, tangent to the Y-point circle. (5) The separatrix current sheet between the closed magnetosphere and the open magnetic field lines is sustained by the electron backflow from the Y-point cloud. Its thickness is self-regulated to marginal charge starvation. (6) Only a small fraction of dissipation occurs in the separatrix inward of the Y-point. A much higher power is released in the equatorial plane, especially at the Y-point where the created dense $e^\pm$ plasma is spun up and intermittently ejected through the nozzle between the two open magnetic fluxes. △ Less

Submitted 17 October, 2022; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 16 pages, ApJ, in press

Relativistic wind from a magnetic monopole rotator

Authors: Rui Hu, Andrei M. Beloborodov, Alexander Y. Chen

Abstract: A rotating star with a monopole (or split monopole) magnetic field gives the simplest, prototype model of a rotationally driven stellar wind. Winds from compact objects, in particular neutron stars, carry strong magnetic fields with modest plasma loading, and develop ultra-relativistic speeds. We investigate the relativistic wind launched from a dense, gravitationally bound, atmosphere on the stel… ▽ More A rotating star with a monopole (or split monopole) magnetic field gives the simplest, prototype model of a rotationally driven stellar wind. Winds from compact objects, in particular neutron stars, carry strong magnetic fields with modest plasma loading, and develop ultra-relativistic speeds. We investigate the relativistic wind launched from a dense, gravitationally bound, atmosphere on the stellar surface. We first examine the problem analytically and then perform global kinetic plasma simulations. Our results show how the wind acceleration mechanism changes from centrifugal (magnetohydrodynamic) to electrostatic (charge-separated) depending on the parameters of the problem. The two regimes give winds with different angular distributions and different scalings with the magnetization parameter. △ Less

Submitted 8 September, 2021; originally announced September 2021.

Comments: 14 pages, submitted to ApJ

Can a strong radio burst escape the magnetosphere of a magnetar?

Authors: Andrei M. Beloborodov

Abstract: We examine the possibility that fast radio bursts (FRBs) are emitted inside the magnetosphere of a magnetar. On its way out, the radio wave must interact with a low-density $e^\pm$ plasma in the outer magnetosphere at radii $R=10^9$-$10^{10}\,$cm. In this region, the magnetospheric particles have a huge cross section for scattering the wave. As a result, the wave strongly interacts with the magnet… ▽ More We examine the possibility that fast radio bursts (FRBs) are emitted inside the magnetosphere of a magnetar. On its way out, the radio wave must interact with a low-density $e^\pm$ plasma in the outer magnetosphere at radii $R=10^9$-$10^{10}\,$cm. In this region, the magnetospheric particles have a huge cross section for scattering the wave. As a result, the wave strongly interacts with the magnetosphere and compresses it, depositing the FRB energy into the compressed field and the scattered radiation. The scattered spectrum extends to the $γ$-ray band and triggers $e^\pm$ avalanche, further boosting the opacity. These processes choke FRBs, disfavoring scenarios with a radio source confined at $R\ll 10^{10}\,$cm. Observed FRBs can be emitted by magnetospheric flare ejecta transporting energy to large radii. △ Less

Submitted 18 October, 2021; v1 submitted 17 August, 2021; originally announced August 2021.

Comments: 6 pages, 1 figure, accepted to ApJL

Scattering of ultrastrong electromagnetic waves by magnetized particles

Authors: Andrei M. Beloborodov

Abstract: Observations of powerful radio waves from neutron star magnetospheres raise the question of how strong waves interact with particles in a strong background magnetic field $B_{bg}$. This problem is examined by solving the particle motion in the wave. Remarkably, waves with amplitudes $E_0>B_{bg}$ pump particle energy via repeating resonance events, quickly reaching the radiation reaction limit. As… ▽ More Observations of powerful radio waves from neutron star magnetospheres raise the question of how strong waves interact with particles in a strong background magnetic field $B_{bg}$. This problem is examined by solving the particle motion in the wave. Remarkably, waves with amplitudes $E_0>B_{bg}$ pump particle energy via repeating resonance events, quickly reaching the radiation reaction limit. As a result, the wave is scattered with a huge cross section. This fact has implications for models of fast radio bursts and magnetars. Particles accelerated in the wave emit gamma-rays, which can trigger an $e^\pm$ avalanche and, instead of silent escape, the wave will produce X-ray fireworks. △ Less

Submitted 7 April, 2022; v1 submitted 11 August, 2021; originally announced August 2021.

Comments: 5 pages, PRL, in press

Coherent Electromagnetic Emission from Relativistic Magnetized Shocks

Authors: Lorenzo Sironi, Illya Plotnikov, Joonas Nättilä, Andrei M. Beloborodov

Abstract: Relativistic magnetized shocks are a natural source of coherent emission, offering a plausible radiative mechanism for Fast Radio Bursts (FRBs). We present first-principles 3D simulations that provide essential information for the FRB models based on shocks: the emission efficiency, spectrum, and polarization. The simulated shock propagates in an $e^\pm$ plasma with magnetization $σ>1$. The measur… ▽ More Relativistic magnetized shocks are a natural source of coherent emission, offering a plausible radiative mechanism for Fast Radio Bursts (FRBs). We present first-principles 3D simulations that provide essential information for the FRB models based on shocks: the emission efficiency, spectrum, and polarization. The simulated shock propagates in an $e^\pm$ plasma with magnetization $σ>1$. The measured fraction of shock energy converted to coherent radiation is $\simeq 10^{-3} \, σ^{-1}$, and the energy-carrying wavenumber of the wave spectrum is $\simeq 4 \,ω_{\rm c}/c$, where $ω_{\rm c}$ is the upstream gyrofrequency. The ratio of the O-mode and X-mode energy fluxes emitted by the shock is $\simeq 0.4\,σ^{-1}$. The dominance of the X-mode at $σ\gg 1$ is particularly strong, approaching 100% in the spectral band around $2\,ω_{\rm c}$. We also provide a detailed description of the emission mechanism for both X- and O-modes. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 6 pages, 5 figures, accepted to PRL

Comptonization by Reconnection Plasmoids in Black Hole Coronae I: Magnetically Dominated Pair Plasma

Authors: Navin Sridhar, Lorenzo Sironi, Andrei M. Beloborodov

Abstract: We perform two-dimensional particle-in-cell simulations of reconnection in magnetically dominated electron-positron plasmas subject to strong Compton cooling. We vary the magnetization $σ\gg1$, defined as the ratio of magnetic tension to plasma inertia, and the strength of cooling losses. Magnetic reconnection under such conditions can operate in magnetically dominated coronae around accreting bla… ▽ More We perform two-dimensional particle-in-cell simulations of reconnection in magnetically dominated electron-positron plasmas subject to strong Compton cooling. We vary the magnetization $σ\gg1$, defined as the ratio of magnetic tension to plasma inertia, and the strength of cooling losses. Magnetic reconnection under such conditions can operate in magnetically dominated coronae around accreting black holes, which produce hard X-rays through Comptonization of seed soft photons. We find that the particle energy spectrum is dominated by a peak at mildly relativistic energies, which results from bulk motions of cooled plasmoids. The peak has a quasi-Maxwellian shape with an effective temperature of $\sim 100$ keV, which depends only weakly on the flow magnetization and the strength of radiative cooling. The mean bulk energy of the reconnected plasma is roughly independent of $σ$, whereas the variance is larger for higher magnetizations. The spectra also display a high-energy tail, which receives $\sim 25$% of the dissipated reconnection power for $σ=10$ and $\sim 40$% for $σ=40$. We complement our particle-in-cell studies with a Monte-Carlo simulation of the transfer of seed soft photons through the reconnection layer, and find the escaping X-ray spectrum. The simulation demonstrates that Comptonization is dominated by the bulk motions in the chain of Compton-cooled plasmoids and, for $σ\sim 10$, yields a spectrum consistent with the typical hard state of accreting black holes. △ Less

Submitted 14 September, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: Accepted for publication in MNRAS. 17 pages, 11 figures, 4 appendices, 1 table

Synchrotron-Self-Compton radiation from magnetically-dominated turbulent plasmas in relativistic jets

Authors: Emanuele Sobacchi, Lorenzo Sironi, Andrei M. Beloborodov

Abstract: Relativistic jets launched by rotating black holes are powerful emitters of non-thermal radiation. Extraction of the rotational energy via electromagnetic stresses produces magnetically-dominated jets, which may become turbulent. Studies of magnetically-dominated plasma turbulence from first principles show that most of the accelerated particles have small pitch angles, i.e. the particle velocity… ▽ More Relativistic jets launched by rotating black holes are powerful emitters of non-thermal radiation. Extraction of the rotational energy via electromagnetic stresses produces magnetically-dominated jets, which may become turbulent. Studies of magnetically-dominated plasma turbulence from first principles show that most of the accelerated particles have small pitch angles, i.e. the particle velocity is nearly aligned with the local magnetic field. We examine synchrotron-self-Compton radiation from anisotropic particles in the fast cooling regime. The small pitch angles reduce the synchrotron cooling rate and promote the role of inverse Compton (IC) cooling, which can occur in two different regimes. In the Thomson regime, both synchrotron and IC components have soft spectra, $νF_ν\proptoν^{1/2}$. In the Klein-Nishina regime, synchrotron radiation has a hard spectrum, typically $νF_ν\proptoν$, over a broad range of frequencies. Our results have implications for the modelling of BL Lacs and Gamma-Ray Bursts (GRBs). BL Lacs produce soft synchrotron and IC spectra, as expected when Klein-Nishina effects are minor. The observed synchrotron and IC luminosities are typically comparable, which indicates a moderate anisotropy with pitch angles $θ\gtrsim0.1$. Rare orphan gamma-ray flares may be produced when $θ\ll0.1$. The hard spectra of GRBs may be consistent with synchrotron radiation when the emitting particles are IC cooling in the Klein-Nishina regime, as expected for pitch angles $θ\sim0.1$. Blazar and GRB spectra can be explained by turbulent jets with a similar electron plasma magnetisation parameter, $σ_{\rm e}\sim10^4$, which for electron-proton plasmas corresponds to an overall magnetisation $σ=(m_{\rm e}/m_{\rm p})σ_{\rm e}\sim10$. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: accepted for publication by MNRAS

Fast dissipation of Colliding Alfvén Waves in a Magnetically Dominated Plasma

Authors: Xinyu Li, Andrei M. Beloborodov, Lorenzo Sironi

Abstract: Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object luminosity. We describe a dissipation mechanism which works faster than magnetic reconnection. The mechanism involves two strong Alfvén waves with anti-aligned magnetic fields $\boldsymbol{B}_1$ and $\boldsymbol{B}_2$ that propagate in opposite directions along background magnetic… ▽ More Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object luminosity. We describe a dissipation mechanism which works faster than magnetic reconnection. The mechanism involves two strong Alfvén waves with anti-aligned magnetic fields $\boldsymbol{B}_1$ and $\boldsymbol{B}_2$ that propagate in opposite directions along background magnetic field $\boldsymbol{B}_0$ and collide. The collision forms a thin current sheet perpendicular to $\boldsymbol{B}_0$, which absorbs the incoming waves. The current sheet is sustained by electric field $\boldsymbol{E}$ breaking the magnetohydrodynamic condition $E<B$ and accelerating particles to high energies. We demonstrate this mechanism with kinetic plasma simulations using a simple setup of two symmetric plane waves with amplitude $A=B_1/B_0=B_2/B_0$ propagating in a uniform $\boldsymbol{B}_0$. The mechanism is activated when $A>1/2$. It dissipates a large fraction of the wave energy, $f=(2A-1)/A^2$, reaching $100\%$ when $A=1$. The plane geometry allows one to see the dissipation process in a one-dimensional simulation. We also perform two-dimensional simulations, enabling spontaneous breaking of the plane symmetry by the tearing instability of the current sheet. At moderate $A$ of main interest the tearing instability is suppressed. Dissipation transitions to normal, slower, magnetic reconnection at $A\gg 1$. The fast dissipation described in this paper may occur in various objects with perturbed magnetic fields, including magnetars, jets from accreting black holes, and pulsar wind nebulae. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Comments: 12 pages, 10 figures, submitted to ApJ

Conceptual design of the Spin Physics Detector

Authors: V. M. Abazov, V. Abramov, L. G. Afanasyev, R. R. Akhunzyanov, A. V. Akindinov, N. Akopov, I. G. Alekseev, A. M. Aleshko, V. Yu. Alexakhin, G. D. Alexeev, M. Alexeev, A. Amoroso, I. V. Anikin, V. F. Andreev, V. A. Anosov, A. B. Arbuzov, N. I. Azorskiy, A. A. Baldin, V. V. Balandina, E. G. Baldina, M. Yu. Barabanov, S. G. Barsov, V. A. Baskov, A. N. Beloborodov, I. N. Belov , et al. (270 additional authors not shown)

Abstract: The Spin Physics Detector, a universal facility for studying the nucleon spin structure and other spin-related phenomena with polarized proton and deuteron beams, is proposed to be placed in one of the two interaction points of the NICA collider that is under construction at the Joint Institute for Nuclear Research (Dubna, Russia). At the heart of the project there is huge experience with polarize… ▽ More The Spin Physics Detector, a universal facility for studying the nucleon spin structure and other spin-related phenomena with polarized proton and deuteron beams, is proposed to be placed in one of the two interaction points of the NICA collider that is under construction at the Joint Institute for Nuclear Research (Dubna, Russia). At the heart of the project there is huge experience with polarized beams at JINR. The main objective of the proposed experiment is the comprehensive study of the unpolarized and polarized gluon content of the nucleon. Spin measurements at the Spin Physics Detector at the NICA collider have bright perspectives to make a unique contribution and challenge our understanding of the spin structure of the nucleon. In this document the Conceptual Design of the Spin Physics Detector is presented. △ Less

Submitted 2 February, 2022; v1 submitted 31 January, 2021; originally announced February 2021.

Radiative turbulent flares in magnetically-dominated plasmas

Authors: J. Nättilä, A. M. Beloborodov

Abstract: We perform 2D and 3D kinetic simulations of reconnection-mediated turbulent flares in a magnetized electron-positron plasma, with weak and strong radiative cooling. Such flares can be generated around neutron stars and accreting black holes. We focus on the magnetically-dominated regime where tension of the background magnetic field lines exceeds the plasma rest mass by a factor $σ_0 > 1$. In the… ▽ More We perform 2D and 3D kinetic simulations of reconnection-mediated turbulent flares in a magnetized electron-positron plasma, with weak and strong radiative cooling. Such flares can be generated around neutron stars and accreting black holes. We focus on the magnetically-dominated regime where tension of the background magnetic field lines exceeds the plasma rest mass by a factor $σ_0 > 1$. In the simulations, turbulence is excited on a macroscopic scale $l_0$, and we observe that it develops by forming thin, dynamic current sheets on various scales. The deposited macroscopic energy dissipates by energizing thermal and nonthermal particles. The particle energy distribution is shaped by impulsive acceleration in reconnecting current sheets, gradual stochastic acceleration, and radiative losses. We parameterize radiative cooling by the ratio $A$ of light-crossing time $l_0/c$ to a cooling timescale, and study the effect of increasing $A$ on the flare. When radiative losses are sufficiently weak, $A<σ_0^{-1}$, the produced emission is dominated by stochastically accelerated particles, and the radiative power depends logarithmically on $A$. The resulting flare radiation spectrum is broad and anisotropic. In the strong-cooling regime, $A>σ_0^{-1}$, stochastic acceleration is suppressed while impulsive acceleration in the current sheets continues to operate. As $A$ increases further, the emission becomes dominated by thermal particles. Our simulations offer a new tool to study particle acceleration by turbulence, especially at high energies, where cooling competes with acceleration. We find that the particle distribution is influenced by strong intermittency of dissipation, and stochastic acceleration cannot be described by a universal diffusion coefficient. △ Less

Submitted 29 September, 2021; v1 submitted 5 December, 2020; originally announced December 2020.

Comments: 22 pages, 17 figures. Accepted to ApJ

A first-principle simulation of blast wave emergence at the photosphere of a neutron star merger

Authors: Christoffer Lundman, Andrei Beloborodov

Abstract: We present the first ab initio simulation of a radiation-mediated shock emerging at the photosphere of a relativistic outflow. The simulation is performed using our code radshock that follows fluid dynamics coupled to time-dependent radiative transfer, calculated with the Monte-Carlo method. We use the code to examine the radiative blast wave emerging from neutron star merger GW~170817. It was pre… ▽ More We present the first ab initio simulation of a radiation-mediated shock emerging at the photosphere of a relativistic outflow. The simulation is performed using our code radshock that follows fluid dynamics coupled to time-dependent radiative transfer, calculated with the Monte-Carlo method. We use the code to examine the radiative blast wave emerging from neutron star merger GW~170817. It was previously proposed that the merger ejected a dark, relativistically expanding, homologous envelope, and then an explosion inside the envelope produced the observed gamma-ray burst GRB~170817A. Our simulation demonstrates how the shock wave generates radiation as it propagates through the envelope, approaches its photosphere, releases the radiation, and collapses, splitting into two collisionless shocks of a microscopic thickness. We find the light curve and the spectral evolution of the produced gamma-ray burst; both are similar to the observed GRB~170817A. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: 5 pages, submitted to ApJL

Emission of magnetar bursts and precursors of neutron star mergers

Authors: Andrei M. Beloborodov

Abstract: Magnetar bursts can be emitted by Alfvén waves growing in the outer magnetosphere to nonlinear amplitudes, $δB/B\sim 1$, and triggering magnetic reconnection. Similar magnetic flares should occur quasi-periodically in a magnetized neutron star binary nearing merger. In both cases, fast dissipation in the magnetic flare creates optically thick $e^\pm$ plasma, whose heat capacity is negligible compa… ▽ More Magnetar bursts can be emitted by Alfvén waves growing in the outer magnetosphere to nonlinear amplitudes, $δB/B\sim 1$, and triggering magnetic reconnection. Similar magnetic flares should occur quasi-periodically in a magnetized neutron star binary nearing merger. In both cases, fast dissipation in the magnetic flare creates optically thick $e^\pm$ plasma, whose heat capacity is negligible compared with the generated radiation energy. Magnetic dissipation then involves photon viscosity and acts through Compton drag on the plasma bulk motions in the reconnection region. The effective temperature of the resulting Comptonization process is self-regulated to tens of keV. The generated X-ray emission is calculated using time-dependent radiative transfer simulations, which follow the creation of $e^\pm$ pairs and the production, Comptonization, and escape of photons. The simulations show how the dissipation region becomes dressed in an $e^\pm$ coat, and how the escaping spectrum is shaped by radiative transfer through the coat. The results are compared with observed magnetar bursts, including the recent activity of SGR 1935+2154 accompanied by a fast radio burst. Predictions are made for X-ray precursors of magnetized neutron star mergers. △ Less

Submitted 24 July, 2021; v1 submitted 14 November, 2020; originally announced November 2020.

Comments: 14 pages, accepted to ApJ

Relativistic Alfvén Waves Entering Charge Starvation in the Magnetospheres of Neutron Stars

Authors: Alexander Y. Chen, Yajie Yuan, Andrei M. Beloborodov, Xinyu Li

Abstract: Instabilities in a neutron star can generate Alfvén waves in its magnetosphere. Propagation along the curved magnetic field lines strongly shears the wave, boosting its electric current $j_{\rm A}$. We derive an analytic expression for the evolution of the wave vector $\boldsymbol{k}$ and the growth of $j_{\rm A}$. In the strongly sheared regime, $j_{\rm A}$ may exceed the maximum current $j_{0}$… ▽ More Instabilities in a neutron star can generate Alfvén waves in its magnetosphere. Propagation along the curved magnetic field lines strongly shears the wave, boosting its electric current $j_{\rm A}$. We derive an analytic expression for the evolution of the wave vector $\boldsymbol{k}$ and the growth of $j_{\rm A}$. In the strongly sheared regime, $j_{\rm A}$ may exceed the maximum current $j_{0}$ that can be supported by the background $e^{\pm}$ plasma. We investigate these "charge-starved" waves, first using a simplified two-fluid analytic model, then with first-principles kinetic simulations. We find that the Alfvén wave continues to propagate successfully even when $κ\equiv j_{\rm A}/j_{0} \gg 1$. It sustains $j_{\rm A}$ by compressing and advecting the plasma along the magnetic field lines with particle Lorentz factors $\sim κ^{1/2}$. The simulations show how plasma instabilities lead to gradual dissipation of the wave energy, giving a dissipation power $L_{\rm diss}\sim 10^{35}(κ/100)^{1/2} (B_w/10^{11}\,{\rm G})\,\mathrm{erg/s}$, where $B_w$ is the wave amplitude. Our results imply that dissipation due to charge starvation is not sufficient to power observed fast radio bursts (FRBs), in contrast to recent proposals. △ Less

Submitted 29 October, 2020; originally announced October 2020.

Comments: 9 pages, 5 figures, submitted to ApJ

Self-modulation of Fast Radio Bursts

Authors: Emanuele Sobacchi, Yuri Lyubarsky, Andrei M. Beloborodov, Lorenzo Sironi

Abstract: Fast Radio Bursts (FRBs) are extreme astrophysical phenomena entering the realm of non-linear optics, a field developed in laser physics. A classical non-linear effect is self-modulation. We examine the propagation of FRBs through the circumburst environment using the idealised setup of a monochromatic linearly-polarised GHz wave propagating through a uniform plasma slab of density $N$ at distance… ▽ More Fast Radio Bursts (FRBs) are extreme astrophysical phenomena entering the realm of non-linear optics, a field developed in laser physics. A classical non-linear effect is self-modulation. We examine the propagation of FRBs through the circumburst environment using the idealised setup of a monochromatic linearly-polarised GHz wave propagating through a uniform plasma slab of density $N$ at distance $R$ from the source. We find that self-modulation occurs if the slab is located within a critical radius $R_{\rm crit}\sim 10^{17}(N/10^2{\rm\; cm}^{-3})(L/10^{42}{\rm\; erg\; s}^{-1}){\rm\; cm}$, where $L$ is the isotropic equivalent of the FRB luminosity. Self-modulation breaks the burst into pancakes transverse to the radial direction. When $R\lesssim R_{\rm crit}$, the transverse size of the pancakes is smaller than the Fresnel scale. The pancakes are strongly diffracted as the burst exits the slab, and interference between the pancakes produces a frequency modulation of the observed intensity with a sub-GHz bandwidth. When $R\sim R_{\rm crit}$, the transverse size of the pancakes becomes comparable with the Fresnel scale, and the effect of diffraction is weaker. The observed intensity is modulated on a timescale of ten microseconds, which corresponds to the radial width of the pancakes. Our results suggest that self-modulation may cause the temporal and frequency structure observed in FRBs. △ Less

Submitted 16 October, 2020; originally announced October 2020.

Comments: accepted for publication in MNRAS

Plasmoid ejection by Alfvén waves and the fast radio bursts from SGR 1935+2154

Authors: Yajie Yuan, Andrei M. Beloborodov, Alexander Y. Chen, Yuri Levin

Abstract: Using numerical simulations we show that low-amplitude Alfvén waves from a magnetar quake propagate to the outer magnetosphere and convert to "plasmoids" (closed magnetic loops) which accelerate from the star, driving blast waves into the magnetar wind. Quickly after its formation, the plasmoid becomes a thin relativistic pancake. It pushes out the magnetospheric field lines, and they gradually re… ▽ More Using numerical simulations we show that low-amplitude Alfvén waves from a magnetar quake propagate to the outer magnetosphere and convert to "plasmoids" (closed magnetic loops) which accelerate from the star, driving blast waves into the magnetar wind. Quickly after its formation, the plasmoid becomes a thin relativistic pancake. It pushes out the magnetospheric field lines, and they gradually reconnect behind the pancake, generating a variable wind far stronger than the normal spindown wind of the magnetar. Repeating ejections drive blast waves in the amplified wind. We suggest that these ejections generate the simultaneous X-ray and radio bursts detected from SGR 1935+2154. A modest energy budget of the magnetospheric perturbation $\sim 10^{40}$ erg is sufficient to produce the observed bursts. Our simulation predicts a narrow (a few ms) X-ray spike from the magnetosphere, arriving almost simultaneously with the radio burst emitted far outside the magnetosphere. This timing is caused by the extreme relativistic motion of the ejecta. △ Less

Submitted 26 August, 2020; v1 submitted 5 June, 2020; originally announced June 2020.

Comments: 8 pages, 5 figures, ApJL in press

Journal ref: The Astrophysical Journal Letters 900 L21 (2020)

Precessing flaring magnetar as a source of repeating FRB 180916.J0158+65

Authors: Yuri Levin, Andrei M. Beloborodov, Ashley Bransgrove

Abstract: Recently, CHIME detected periodicity in the bursting rate of the repeating FRB 180916.J0158+65. In a popular class of models, the fast radio bursts (FRBs) are created by giant magnetic flares of a hyper-active magnetar driven by fast ambipolar diffusion in the core. We point out that in this scenario the magnetar is expected to precess freely with a period of hours to weeks. The internal magnetic… ▽ More Recently, CHIME detected periodicity in the bursting rate of the repeating FRB 180916.J0158+65. In a popular class of models, the fast radio bursts (FRBs) are created by giant magnetic flares of a hyper-active magnetar driven by fast ambipolar diffusion in the core. We point out that in this scenario the magnetar is expected to precess freely with a period of hours to weeks. The internal magnetic field $B\sim 10^{16}$G deforms the star, and magnetic flares induce sudden changes in magnetic stresses. The resulting torques and displacements of the principal axes of inertia are capable of pumping a significant amplitude of precession. The anisotropy of the flaring FRB activity, combined with precession, implies a strong periodic modulation of the visible bursting rate. The ultra-strong field invoked in the magnetar model provides: (1) energy for the frequent giant flares, (2) the high rate of ambipolar diffusion, releasing the magnetic energy on the timescale $\sim 10^9$s, (3) the core temperature $T\approx 10^9$K, likely above the critical temperature for neutron superfluidity, (4) strong magnetospheric torques, which efficiently spin down the star, and (5) deformation with ellipticity $ε> 10^{-6}$, much greater than the rotational deformation. These conditions result in a precession with negligible viscous damping, and can explain the observed 16 day period in FRB 180916.J0158+65. △ Less

Submitted 23 April, 2020; v1 submitted 11 February, 2020; originally announced February 2020.

Comments: 5 pages, improved discussion of excitation of free precession, accepted to ApJ Letters

A Quake Quenching the Vela Pulsar

Authors: Ashley Bransgrove, Andrei M. Beloborodov, Yuri Levin

Abstract: The remarkable null pulse coincident with the 2016 glitch in Vela rotation indicates a dynamical event involving the crust and the magnetosphere of the neutron star. We propose that a crustal quake associated with the glitch strongly disturbed the Vela magnetosphere and thus interrupted its radio emission. We present the first global numerical simulations of a neutron starquake. Our code resolves… ▽ More The remarkable null pulse coincident with the 2016 glitch in Vela rotation indicates a dynamical event involving the crust and the magnetosphere of the neutron star. We propose that a crustal quake associated with the glitch strongly disturbed the Vela magnetosphere and thus interrupted its radio emission. We present the first global numerical simulations of a neutron starquake. Our code resolves the elasto-dynamics of the entire crust and follows the evolution of Alfvén waves excited in the magnetosphere. We observe Rayleigh surface waves propagating away from the epicenter of the quake, around the circumference of the crust - an instance of the so-called whispering gallery modes. The Rayleigh waves set the initial spatial scale of the magnetospheric disturbance. Once launched, the Alfvén waves bounce in the closed magnetosphere, become de-phased, and generate strong electric currents, capable of igniting electric discharge. Most likely, the discharge floods the magnetosphere with electron-positron plasma, quenching the radio emission. We find that the observed $\sim 0.2$ s disturbance is consistent with the damping time of the crustal waves if the crust is magnetically coupled to the superconducting core of the neutron star. The quake is expected to produce a weak X-ray burst of short duration. △ Less

Submitted 30 September, 2020; v1 submitted 23 January, 2020; originally announced January 2020.

Comments: 21 pages, 10 figures, accepted to ApJ

Kinetic Simulations of Radiative Magnetic Reconnection in the Coronae of Accreting Black Holes

Authors: Lorenzo Sironi, Andrei M. Beloborodov

Abstract: We perform two- and three-dimensional particle-in-cell simulations of reconnection in magnetically-dominated pair plasmas subject to strong Compton cooling. Reconnection under such conditions operates in accretion disk coronae around black holes, which produce hard X-rays through Comptonization. Our simulations show that most of the plasma in the reconnection layer is kept cold by Compton losses a… ▽ More We perform two- and three-dimensional particle-in-cell simulations of reconnection in magnetically-dominated pair plasmas subject to strong Compton cooling. Reconnection under such conditions operates in accretion disk coronae around black holes, which produce hard X-rays through Comptonization. Our simulations show that most of the plasma in the reconnection layer is kept cold by Compton losses and locked in magnetically-dominated plasmoids with small thermal pressure. Compton drag clears cavities inside plasmoids and affects their bulk motions. These effects, however, weakly change the reconnection rate and the plasmoid size distribution from those in non-radiative reconnection. This demonstrates that the reconnection dynamics is governed by similar magnetic stresses in both cases and weakly affected by thermal pressure. We examine the energy distribution of particles energized by radiative reconnection and observe two components. (1) A mildly-relativistic peak, which results from bulk motions of cooled plasmoids. This component receives most of the dissipated reconnection power and dominates the output X-ray emission. The peak has a quasi-Maxwellian shape with an effective temperature of 100 keV. Thus, it mimics thermal Comptonization used previously to fit hard-state spectra of accreting black holes. (2) A high-energy tail, which receives 20% of dissipated reconnection power. It is populated by particles accelerated impulsively at X-points or "picked up" by fast outflows from X-points. The high-energy particles immediately cool, and their inverse Compton emission explains the MeV spectral tail detected in the hard state of Cyg X-1. Our first-principle simulations support reconnection as a mechanism powering hard X-ray emission from accreting black holes. △ Less

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

Comments: 18 pages, 15 figures, ApJ accepted; added discussion of 3D simulations

Blast Waves from Magnetar Flares and Fast Radio Bursts

Authors: Andrei M. Beloborodov

Abstract: Magnetars younger than one century are expected to be hyper active. Besides winds powered by rotation they generate frequent magnetic flares, which launch powerful blast waves into the wind. These internal shocks act as masers producing fast (millisecond) radio bursts (FRBs) with the following properties. (1) GHz radio emission occurs at radii $r\sim 10^{14}$ cm and lasts $\lesssim 1$ ms in observ… ▽ More Magnetars younger than one century are expected to be hyper active. Besides winds powered by rotation they generate frequent magnetic flares, which launch powerful blast waves into the wind. These internal shocks act as masers producing fast (millisecond) radio bursts (FRBs) with the following properties. (1) GHz radio emission occurs at radii $r\sim 10^{14}$ cm and lasts $\lesssim 1$ ms in observer's time. (2) Induced scattering in the surrounding wind does not suppress the radio burst. (3) The emission has linear polarization set by the magnetar rotation axis. (4) The emission drifts to lower frequencies during the burst, and its duration broadens at lower frequencies. (5) Blast waves in inhomogeneous winds may emit variable bursts; periodicity might appear on sub-ms timescales if the magnetar rotates with $\sim 1$ s period. However, the observed FRB structure is likely changed by lensing effects during propagation through the host galaxy. (6) The FRBs from magnetars are expected to repeat, with rare strong bursts (up to $\sim 10^{43}$ erg) or more frequent weak bursts. (7) When a repeating flare strikes the wind bubble in the tail of a previous flare, the FRB turns into a bright optical flash. Its luminosity may approach that of a supernova Ia and last seconds. The rate of these optical flashes in the universe is much lower than the FRB rate, however it may exceed the supernova rate. Locations of hyper-active magnetars in their host galaxies depend on how they form: magnetars created in supernovae explosions will trace star formation regions, and magnetars formed in mergers of compact objects will be offset. The merger magnetars are expected to be most energetic and particularly hyper-active. △ Less

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

Comments: 20 pages, 6 figures, ApJ, in press

X-ray Monitoring of the Magnetar CXOU J171405.7-381031 in SNR CTB 37B

Authors: E. V. Gotthelf, J. P. Halpern, K. Mori, A. M. Beloborodov

Abstract: We present the results of our 8 year X-ray monitoring campaign on CXOU J171405.7-381031, the magnetar associated with the faint supernova remnant (SNR) CTB 37B. It is among the youngest by inferred spin-down age, and most energetic in spin-down power of magnetars, and may contribute, at least partially, to the GeV and TeV emission coincident with the SNR. We use a series of Chandra, XMM-Newton, an… ▽ More We present the results of our 8 year X-ray monitoring campaign on CXOU J171405.7-381031, the magnetar associated with the faint supernova remnant (SNR) CTB 37B. It is among the youngest by inferred spin-down age, and most energetic in spin-down power of magnetars, and may contribute, at least partially, to the GeV and TeV emission coincident with the SNR. We use a series of Chandra, XMM-Newton, and NuSTAR observations to characterize the timing and spectral properties of the magnetar. The spin-down rate of the pulsar almost doubled in <1 year and then decreased slowly to a more stable value. Its X-ray flux varied by approx, 50%, possibly correlated with the spin down rate. The 1-79 keV spectrum is well-characterized by an absorbed blackbody plus power-law model with an average temperature of kT=0.62+/-0.04 keV and photon index Gamma=0.92+/-0.16, or by a Comptonized blackbody with kT=0.55+/-0.04 keV and an additional hard power law with Gamma=0.70+/-0.20, In contrast with most magnetars, the pulsed signal is found to decrease with energy up to 6 keV, which is apparently caused by mixing with the hard spectral component that is pulse-phase shifted by approx. 0.43 cycles from the soft X-rays. We also analyze the spectrum of the nearby, diffuse nonthermal source XMMU J171410.8-381442, whose relation to the SNR is uncertain. △ Less

Submitted 30 July, 2019; originally announced July 2019.

Comments: 10 pages, 10 figures, 2 Tables; Latex, aastex63 style. To appear in the Astrophysical Journal

Pulsar wind-heated accretion disk and the origin of modes in transitional millisecond pulsar PSR J1023+0038

Authors: Alexandra Veledina, Joonas Nättilä, Andrei M. Beloborodov

Abstract: Transitional millisecond pulsars provide a unique set of observational data for understanding accretion at low rates onto magnetized neutron stars. In particular, PSR~J1023+0038 exhibits a remarkable bimodality of the X-ray luminosity (low and high modes), pulsations extending from the X-ray to the optical band, GeV emission, and occasional X-ray flares. We discuss a scenario for the pulsar intera… ▽ More Transitional millisecond pulsars provide a unique set of observational data for understanding accretion at low rates onto magnetized neutron stars. In particular, PSR~J1023+0038 exhibits a remarkable bimodality of the X-ray luminosity (low and high modes), pulsations extending from the X-ray to the optical band, GeV emission, and occasional X-ray flares. We discuss a scenario for the pulsar interaction with the accretion disk capable of explaining the observed behavior. We suggest that during the high mode the disk is truncated outside the light cylinder, allowing the pulsar wind to develop near the equatorial plane and strike the disk. The dissipative wind-disk collision energizes the disk particles and generates synchrotron emission, which peaks in the X-ray band and extends down to the optical band. The emission is modulated by the pulsar wind rotation, resulting in a pulse profile with two peaks 180\degr\ apart. This picture explains the high-mode luminosity, spectrum, and pulse profile (X-ray and optical) of PSR J1023+0038. It may also explain the X-ray flares as events of sudden increase in the effective disk cross section intercepting the wind. In contrast to previously proposed models, we suggest that the disk penetrates the light cylinder only during the low X-ray mode. This penetration suppresses the dissipation caused by the pulsar wind-disk collision, and the system enters the propeller regime. The small duty cycle of the propeller explains the low spindown rate of the pulsar. △ Less

Submitted 6 June, 2019; originally announced June 2019.

Comments: 14 pages, 4 figures, submitted to ApJ

Report number: NORDITA 2019-055

Magnetars as Astrophysical Laboratories of Extreme Quantum Electrodynamics: The Case for a Compton Telescope

Authors: Zorawar Wadiasingh, George Younes, Matthew G. Baring, Alice K. Harding, Peter L. Gonthier, Kun Hu, Alexander van der Horst, Silvia Zane, Chryssa Kouveliotou, Andrei M. Beloborodov, Chanda Prescod-Weinstein, Tanmoy Chattopadhyay, Sunil Chandra, Constantinos Kalapotharakos, Kyle Parfrey, Harsha Blumer, Demos Kazanas

Abstract: A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This p… ▽ More A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: 11 pages, 4 figures, Astro2020 Science White Paper submitted to the National Academies of Sciences

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

Radio Time-Domain Signatures of Magnetar Birth

Authors: Casey J. Law, Ben Margalit, Nipuni T. Palliyaguru, Brian D. Metzger, Lorenzo Sironi, Yong Zheng, Edo Berger, Raffaella Margutti, Andrei Beloborodov, Matt Nicholl, Tarraneh Eftekhari, Indrek Vurm, Peter K. G. Williams

Abstract: The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from millise… ▽ More The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from milliseconds to gigaseconds will be used to address diverse physical questions related to the stellar evolution and death, the circum- and intergalactic medium, and fundamental physics on scales ranging from the quantum to cosmological. We argue that high-sensitivity radio interferometer concepts in development will be critical to accessing this science. There is also a growing need for scientists with expertise in new statistical tools and software techniques. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Science white paper for Astro2020

Gravity and Light: Combining Gravitational Wave and Electromagnetic Observations in the 2020s

Authors: R. J. Foley, K. D. Alexander, I. Andreoni, I. Arcavi, K. Auchettl, J. Barnes, G. Baym, E. C. Bellm, A. M. Beloborodov, N. Blagorodnova, J. P. Blakeslee, P. R. Brady, M. Branchesi, J. S. Brown, N. Butler, M. Cantiello, R. Chornock, D. O. Cook, J. Cooke, D. L. Coppejans, A. Corsi, S. M. Couch, M. W. Coughlin, D. A. Coulter, P. S. Cowperthwaite , et al. (88 additional authors not shown)

Abstract: As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse s… ▽ More As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: submitted to Astro2020

Relativistic envelopes and gamma-rays from neutron star mergers

Authors: Andrei M. Beloborodov, Christoffer Lundman, Yuri Levin

Abstract: We suggest that neutron star mergers eject an ultra-relativistic envelope of mass $m\sim 10^{-7}M_\odot$, which helps explain the gamma-ray burst from GW170817. One ejection mechanism is the ablation of the neutron star surface by the burst of neutrinos in the first $30\,μ$s of the merger. Another, more efficient, mechanism for inflating the ultra-relativistic envelope is an internal shock in the… ▽ More We suggest that neutron star mergers eject an ultra-relativistic envelope of mass $m\sim 10^{-7}M_\odot$, which helps explain the gamma-ray burst from GW170817. One ejection mechanism is the ablation of the neutron star surface by the burst of neutrinos in the first $30\,μ$s of the merger. Another, more efficient, mechanism for inflating the ultra-relativistic envelope is an internal shock in the massive ejecta from the merger. A strong shock is expected if the merger product is a magnetar, which emits a centrifugally accelerated wind. The shock propagates outward through the ejecta and accelerates in its outer layers at radii $r\sim 10^9-10^{10}$cm, launching an ultra-relativistic opaque envelope filled with $\sim 10^4$ photons per nucleon. The Lorentz factor profile of the envelope rises outward and determines its homologous expansion, which adiabatically cools the trapped photons. Once the magnetar loses its differential rotation and collapses into a black hole, a powerful jet forms. It drives a blast wave into the envelope, chasing its outer layers and eventually catching up with the envelope photosphere at $r\sim 10^{12}$cm. The ultra-relativistic photospheric breakout of the delayed blast wave emits a gamma-ray burst in a broad solid angle around the merger axis. This model explains the gamma-ray pulse from merger GW170817 with luminosity $L_γ\sim 10^{47}$erg/s, duration $Δt_{\rm obs}\sim 0.5$s, and characteristic photon energy $\sim 100$keV. The blast wave Lorentz factor at the envelope photosphere is consistent with $Γ\geq 5$ that we derive from the observed light curve of the burst. We suggest future tests of the model. △ Less

Submitted 24 July, 2021; v1 submitted 28 December, 2018; originally announced December 2018.

Comments: 18 pages, 8 figures, published in ApJ 2020

Journal ref: The Astrophysical Journal, Volume 897, Issue 2, id.141 (2020)

Dissipation of Alfvén Waves in Relativistic Magnetospheres of Magnetars

Authors: Xinyu Li, Jonathan Zrake, Andrei M. Beloborodov

Abstract: Magnetar flares excite strong Alfvén waves in the magnetosphere of the neutron star. The wave energy can (1) dissipate in the magnetosphere, (2) convert to "fast modes" and possibly escape, and (3) penetrate the neutron star crust and dissipate there. We examine and compare the three options. Particularly challenging are nonlinear interactions between strong waves, which develop a cascade to small… ▽ More Magnetar flares excite strong Alfvén waves in the magnetosphere of the neutron star. The wave energy can (1) dissipate in the magnetosphere, (2) convert to "fast modes" and possibly escape, and (3) penetrate the neutron star crust and dissipate there. We examine and compare the three options. Particularly challenging are nonlinear interactions between strong waves, which develop a cascade to small dissipative scales. This process can be studied in the framework of force-free electrodynamics (FFE). We perform three-dimensional FFE simulations to investigate Alfvén wave dissipation, how long it takes, and how it depends on the initial wave amplitude on the driving scale. In the simulations, we launch two large Alfvén wave packets that keep bouncing on closed magnetic field lines and collide repeatedly until the full turbulence spectrum develops. Besides dissipation due to the turbulent cascade, we find that in some simulations spurious energy losses occur immediately in the first collisions. This effect occurs in special cases where the FFE description breaks. It is explained with a simple one-dimensional model, which we examine in both FFE and full magnetohydrodynamic settings. We find that magnetospheric dissipation through nonlinear wave interactions is relatively slow and more energy is drained into the neutron star. The wave energy deposited into the star is promptly dissipated through plastic crustal flows induced at the bottom of the liquid ocean, and a fraction of the generated heat is radiated from the stellar surface. △ Less

Submitted 24 October, 2018; originally announced October 2018.

Sub-photospheric turbulence as a heating mechanism in gamma-ray bursts

Authors: Jonathan Zrake, Andrei M. Beloborodov, Christoffer Lundman

Abstract: We examine the possible role of turbulence in feeding the emission of gamma-ray bursts (GRBs). Turbulence may develop in a GRB jet as the result of hydrodynamic or current-driven instabilities. The jet carries dense radiation and the turbulence cascade can be damped by Compton drag, passing kinetic fluid energy to photons through scattering. We identify two regimes of turbulence dissipation: (1) "… ▽ More We examine the possible role of turbulence in feeding the emission of gamma-ray bursts (GRBs). Turbulence may develop in a GRB jet as the result of hydrodynamic or current-driven instabilities. The jet carries dense radiation and the turbulence cascade can be damped by Compton drag, passing kinetic fluid energy to photons through scattering. We identify two regimes of turbulence dissipation: (1) "Viscous" - the turbulence cascade is Compton damped on a scale $\ell_{\rm damp}$ greater than the photon mean free path $\ell_\star$. Then turbulence energy is passed to photons via bulk Comptonization by smooth shear flows on scale $\ell_\star<\ell_{\rm damp}$. (2) "Collisionless" - the cascade avoids Compton damping and extends to microscopic plasma scales much smaller than $\ell_\star$. The collisionless dissipation energizes plasma particles, which radiate the received energy; how the dissipated power is partitioned between particles needs further investigation with kinetic simulations. We show that the dissipation regime switches from viscous to collisionless during the jet expansion, at a critical value of the jet optical depth which depends on the amplitude of turbulence. Turbulent GRB jets are expected to emit nonthermal photospheric radiation. Our analysis also suggests revisions of turbulent Comptonization in black hole accretion disks discussed in previous works. △ Less

Submitted 4 October, 2018; originally announced October 2018.

Radiation mediated shocks in gamma-ray bursts: Subshock photon production

Authors: Christoffer Lundman, Andrei Beloborodov

Abstract: Internal shocks provide a plausible heating mechanism in the jets of gamma-ray bursts (GRBs). Shocks occurring below the jet photosphere are mediated by radiation. It was previously found that radiation mediated shocks (RMSs) inside GRB jets are inefficient photon producers, and the photons that mediate the RMS must originate from an earlier stage of the explosion. We show that this conclusion is… ▽ More Internal shocks provide a plausible heating mechanism in the jets of gamma-ray bursts (GRBs). Shocks occurring below the jet photosphere are mediated by radiation. It was previously found that radiation mediated shocks (RMSs) inside GRB jets are inefficient photon producers, and the photons that mediate the RMS must originate from an earlier stage of the explosion. We show that this conclusion is valid only for non-magnetized jets. RMSs that propagate in moderately magnetized plasma develop a collisionless subshock which locally heats the plasma to a relativistic temperature, and the hot electrons emit copious synchrotron photons inside the RMS. We find that this mechanism is generally effective for mildly relativistic shocks and may be the main source of photons observed in GRBs. We derive a simple analytical formula for the generated photon number per proton, $Z$, which gives $Z=10^5$-$10^6$, consistent with observations. The photons are initially injected with low energies, well below the observed GRB peak. Their number is controlled by two main factors: (1) the abundance of electron-positron pairs created in the shock, which is self-consistently calculated, and (2) the upper limit on the brightness temperature of soft radiation set by induced downscattering. The injected soft photons that survive induced downscattering gain energy in the RMS via bulk Comptonization and shape its nonthermal spectrum. △ Less

Submitted 9 April, 2018; originally announced April 2018.

Comments: 10 pages, submitted to ApJ

A strong limit on the very-high-energy emission from GRB 150323A

Authors: A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, V. Bugaev, M. P. Connolly, W. Cui, M. Errando, A. Falcone, Q. Feng, J. P. Finley, A. Flinders, L. Fortson, A. Furniss, G. H. Gillanders, M. Hütten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson, P. Kaaret , et al. (41 additional authors not shown)

Abstract: On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit i… ▽ More On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant since the very-high-energy (VHE) observation started only $\sim$2 minutes after the prompt emission peaked, and $Fermi$-LAT observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB~150323A ($z=0.593$) limits the attenuation by the extragalactic background light to $\sim 50$ % at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below $\sim100$ GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be $A\gtrsim 3\times 10^{11}$ g cm$^{-1}$, consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the ISM, which therefore cannot be ruled out as the environment of GRB 150323A. △ Less

Submitted 3 March, 2018; originally announced March 2018.

Comments: 6 pages, 1 figure. Accepted for publication in The Astrophysical Journal

Magnetic Field Evolution of Neutron Stars I: Basic formalism, numerical techniques, and first results

Authors: Ashley Bransgrove, Yuri Levin, Andrei Beloborodov

Abstract: In this work we explore the evolution of magnetic fields inside strongly magnetized neutron stars in axisymmetry. We model numerically the coupled field evolution in the core and the crust. Our code models the Hall drift and Ohmic effects in the crust, the back-reaction on the field from magnetically-induced elastic deformation of the crust, the magnetic twist exchange between the crust and the co… ▽ More In this work we explore the evolution of magnetic fields inside strongly magnetized neutron stars in axisymmetry. We model numerically the coupled field evolution in the core and the crust. Our code models the Hall drift and Ohmic effects in the crust, the back-reaction on the field from magnetically-induced elastic deformation of the crust, the magnetic twist exchange between the crust and the core, and the drift of superconducting flux tubes inside the core. The correct hydromagnetic equilibrium is enforced in the core. We find that i) The Hall attractor found by Gourgouliatos and Cumming in the crust exists also for configurations when the B-field penetrates into the core. However, the evolution timescale for the core-penetrating fields is dramatically different than that of the fields confined to the crust. ii) The combination of Jones' flux tube drift and Ohmic dissipation in the crust can deplete the pulsar magnetic fields on the timescale of $150$ Myr if the crust is hot ($T\sim2\times10^8$ K), but acts on much slower timescales for cold neutron stars, such as recycled pulsars ($\sim 1.8 $ Gyr, depending on impurity levels). iii) The outward motion of superfluid vortices during the rapid spin-down of a young highly magnetized pulsar, can result in a partial expulsion of flux from the core when $B\lesssim 10^{13}$ G. However for $B\gtrsim2\times 10^{13}$ G, the combination of a stronger magnetic field and a longer spin period implies that the core field cannot be expelled. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: 22 pages, 10 figures, accepted for publication in MNRAS

A Nova Outburst Powered by Shocks

Authors: Kwan-Lok Li, Brian D. Metzger, Laura Chomiuk, Indrek Vurm, Jay Strader, Thomas Finzell, Andrei M. Beloborodov, Thomas Nelson, Benjamin J. Shappee, Christopher S. Kochanek, Jose L. Prieto, Stella Kafka, Thomas W. -S. Holoien, Todd A. Thompson, Paul J. Luckas, Hiroshi Itoh

Abstract: Classical novae are runaway thermonuclear burning events on the surfaces of accreting white dwarfs in close binary star systems, sometimes appearing as new naked-eye sources in the night sky. The standard model of novae predicts that their optical luminosity derives from energy released near the hot white dwarf which is reprocessed through the ejected material. Recent studies with the Fermi Large… ▽ More Classical novae are runaway thermonuclear burning events on the surfaces of accreting white dwarfs in close binary star systems, sometimes appearing as new naked-eye sources in the night sky. The standard model of novae predicts that their optical luminosity derives from energy released near the hot white dwarf which is reprocessed through the ejected material. Recent studies with the Fermi Large Area Telescope have shown that many classical novae are accompanied by gigaelectronvolt gamma-ray emission. This emission likely originates from strong shocks, providing new insights into the properties of nova outflows and allowing them to be used as laboratories to study the unknown efficiency of particle acceleration in shocks. Here we report gamma-ray and optical observations of the Milky Way nova ASASSN-16ma, which is among the brightest novae ever detected in gamma-rays. The gamma-ray and optical light curves show a remarkable correlation, implying that the majority of the optical light comes from reprocessed emission from shocks rather than the white dwarf. The ratio of gamma-ray to optical flux in ASASSN-16ma directly constrains the acceleration efficiency of non-thermal particles to be ~0.005, favouring hadronic models for the gamma-ray emission. The need to accelerate particles up to energies exceeding 100 gigaelectronvolts provides compelling evidence for magnetic field amplification in the shocks. △ Less

Submitted 3 September, 2017; originally announced September 2017.

Comments: Published in Nature Astronomy. This is the authors' version with 55 pages, 8 figures, and 3 tables

Journal ref: Nature Astronomy, 1, 0222 (2017)

Radiation mediated shocks in gamma-ray bursts: Pair creation

Authors: Christoffer Lundman, Andrei Beloborodov, Indrek Vurm

Abstract: Sub-photospheric shock dissipation is one of the main proposed mechanisms for producing the prompt gamma-ray burst (GRB) emission. Such shocks are mediated by scattering of radiation. We introduce a time dependent, special relativistic code which dynamically couples Monte Carlo radiative transfer to the flow hydrodynamics. The code also self-consistently implements electron-positron pair productio… ▽ More Sub-photospheric shock dissipation is one of the main proposed mechanisms for producing the prompt gamma-ray burst (GRB) emission. Such shocks are mediated by scattering of radiation. We introduce a time dependent, special relativistic code which dynamically couples Monte Carlo radiative transfer to the flow hydrodynamics. The code also self-consistently implements electron-positron pair production and annihilation. We simulate shocks with properties relevant for GRBs and study the steady-state solutions, which are accurate deep below the jet photosphere. The shock generates a power-law photon spectrum through the first-order Fermi mechanism, extending upwards from the typical upstream photon energy. Strong shocks (for which the downstream pressure is much larger than the upstream pressure) have rising $νF_ν$ shock spectra. The spectrum extends up to $ε_{max} \equiv E_{max}/m_e c^2 \sim v^2$ for non-relativistic shocks, where $m_e$ is the electron rest mass and $v$ is the relative speed between the upstream and downstream in units of the speed of light $c$. For mildly relativistic shocks the power law softens at $ε\gtrsim 10^{-1}$ due to Klein-Nishina effects, and shocks with $vγ\gtrsim 1$, where $γ\equiv (1-v^2)^{-1/2}$, produce electron-positron pairs. As an example, a strong shock with $vγ= 3$ and a photon-to-proton ratio of $n_γ/n_p = 2 \times 10^5$ has a peak pair-to-proton ratio of $Z_\pm \approx 225$. The main effect of pairs in a steady-state shock is to decrease its spatial width by a factor of $\sim Z_\pm$. The post-shock spectrum thermalizes in the downstream. In absence of emission and absorption processes, kinetic equilibrium at temperature $θ_d \equiv kT_d/m_e c^2 \approx ε_d/3$ is reached at an optical depth of $τ\gg θ_d^{-1}$ behind the shock, where $ε_d$ is the average downstream photon energy. △ Less

Submitted 8 August, 2017; originally announced August 2017.

Comments: 14 pages, 8 figures, submitted to ApJ

Magnetars

Authors: Victoria M. Kaspi, Andrei Beloborodov

Abstract: Magnetars are young and highly magnetized neutron stars which display a wide array of X-ray activity including short bursts, large outbursts, giant flares and quasi-periodic oscillations, often coupled with interesting timing behavior including enhanced spin-down, glitches and anti-glitches. The bulk of this activity is explained by the evolution and decay of an ultrastrong magnetic field, stressi… ▽ More Magnetars are young and highly magnetized neutron stars which display a wide array of X-ray activity including short bursts, large outbursts, giant flares and quasi-periodic oscillations, often coupled with interesting timing behavior including enhanced spin-down, glitches and anti-glitches. The bulk of this activity is explained by the evolution and decay of an ultrastrong magnetic field, stressing and breaking the neutron star crust, which in turn drives twists of the external magnetosphere and powerful magnetospheric currents. The population of detected magnetars has grown to about 30 objects and shows unambiguous phenomenological connection with very highly magnetized radio pulsars. Recent progress in magnetar theory includes explanation of the hard X-ray component in the magnetar spectrum and development of surface heating models, explaining the sources' remarkable radiative output. △ Less

Submitted 28 February, 2017; originally announced March 2017.

Comments: 40 pages, Annual Review of Astronomy and Astrophysics, in press