Search for: All records

Abstract We present the first 3D global kinetic simulations of the interacting magnetospheres of premerger binary neutron stars. The stars, whose magnetic moments are antialigned, twist the field lines connecting them, leading to periodic eruptions. Each eruption consists of an expanding magnetic flux tube with a reconnecting current sheet trailing behind it, topologically analogous to coronal mass ejections. We predict two novel classes of electromagnetic precursor signals powered by the efficient dissipation of magnetic energy in these periodically forming trailing current sheets. First, particles accelerated in the sheets produce nonthermal gamma-ray signals peaking at ∼16 MeV, which escape minutes to seconds before merger while the sheets are still optically thin to pair production, with modest characteristic luminosities ofL_obs ≳ 10⁴²erg s⁻¹, detectable only for nearby mergers. Second, merging plasmoids in the sheets could produce fast-radio-burst-like transients in the final seconds before merger, with characteristic luminositiesL_radio ∼ 10³⁸⁻⁴⁰erg s⁻¹. These coherent radio precursors would be detectable by upcoming instruments, either in untargeted surveys by wide-field instruments such as CHORD, or through targeted follow-up of gravitational-wave early warning alerts with instruments such as DSA or SKA-mid. 

Abstract We develop a model of particle energization and emission from strongly turbulent black hole coronae. Our local model is based on a set of 2D radiative particle-in-cell simulations with an electron–ion plasma composition, injection and diffusive escape of photons and charged particles, and self-consistent Compton scattering. We show that a radiatively compact turbulent corona generates extended nonthermal ion distributions, while producing X-ray spectra consistent with observations. As an example, we demonstrate excellent agreement with observed X-ray spectra of NGC 4151. The predicted emission spectra feature an MeV tail, which can be studied with future MeV-band instruments. The MeV tail is shaped by nonthermal electrons accelerated at turbulent current sheets. We also find that the corona regulates itself into a two-temperature state, with ions much hotter than electrons. The ions carry away roughly two-thirds of the dissipated power, and their energization is driven by a combination of shocks and reconnecting current sheets, embedded into the turbulent flow. 

Abstract Extended periods of radio pulsations have been observed for six magnetars, displaying characteristics different from those of ordinary pulsars. In this Letter, we argue that radio emission is generated in a closed, twisted magnetic flux bundle originating near the magnetic pole and extending beyond 100 km from the magnetar. The electron–positron flow in the twisted bundle has to carry electric current and, at the same time, experiences a strong drag from the radiation field of the magnetar. This combination forces the plasma into a “radiatively locked” state with a sustained two-stream instability, generating radio emission. We demonstrate this mechanism using novel first-principles simulations that follow the plasma behavior by solving the relativistic Vlasov equation with the discontinuous Galerkin method. First, using one-dimensional simulations, we demonstrate how radiative drag induces the two-stream instability, sustaining turbulent electric fields. When extended to two dimensions, the system produces electromagnetic waves, including superluminal modes capable of escaping the magnetosphere. We measure their frequency and emitted power and incorporate the local simulation results into a global magnetospheric model. The model explains key features of the observed radio emission from magnetars: its appearance after an X-ray outburst, wide pulse profiles, luminosities ∼10³⁰erg s⁻¹, and a broad range of frequencies extending up to ∼100 GHz. 

Abstract Magnetic reconnection in current layers that form intermittently in radiatively inefficient accretion flows onto black holes is a promising mechanism for particle acceleration and high-energy emission. It has been recently proposed that such layers, arising during flux eruption events, can power the rapid TeV flares observed from the core of M87. In this scenario, inverse-Compton scattering of soft radiation from the accretion flow by energetic electron–positron pairs produced near the reconnection layer was suggested as the primary emission mechanism. However, detailed calculations show that radiation from pairs alone cannot account for the GeV emission detected by the Fermi observatory. In this work, we combine analytic estimates with 3D radiative particle-in-cell simulations of pair–proton plasmas to show that the GeV emission can be naturally explained by synchrotron radiation from protons accelerated in the current sheet. Although the exact proton content of the layer is uncertain, our model remains robust across a broad range of proton-to-pair number density ratios. While protons are subdominant in number compared to pairs, our simulations demonstrate that they can be accelerated more efficiently, leading to a self-regulated steady state in which protons dominate the energy budget. Ultimately, proton synchrotron emission accounts for approximately 5%–20% of the total dissipation power. The majority is radiated as MeV photons via pair synchrotron emission, with a smaller fraction emitted as TeV photons through inverse-Compton scattering. 

Abstract We study the propagation of cosmic rays (CRs) through a simulation of magnetohydrodynamic (MHD) turbulence at an unprecedented resolution of 10,240³. We drive turbulence that is subsonic and super-Alfvénic, characterized byδB_rms/B₀ = 2. The high resolution enables an extended inertial range such that the Alfvén scalel_A, whereδB(l_A) ≈ B₀, is well resolved. This allows us to properly capture how the cascade transitions from large amplitudes on large scales to small amplitudes on small scales. We find that sharp bends in the magnetic field are key mediators of particle transport even on small scales via resonant curvature scattering. We further find that particle scattering in the turbulence shows strong hints of self-similarity: (1) the diffusion has a weak energy dependence over almost two decades in particle energy, and (2) the particles’ random walk exhibits a broad power-law distribution of collision times such that the diffusion is dominated by the rarest, long-distance excursions. Our results suggest that large-amplitude MHD turbulence can provide efficient scattering over a wide range of CR energies and may help explain many CR observations above a ∼TeV: the flattening of the B/C spectrum, the hardening of CR primary spectra, and the weak dependence of arrival anisotropy on CR energy.