Search for: All records

Abstract The black hole at the center of M87 is observed to flare regularly in the very high-energy (VHE) band, with photon energies ≳100 GeV. The rapid variability, which can be as short as 2 days in the VHE lightcurve, constrains some of the flares to originate close to the black hole. Magnetic reconnection is a promising candidate for explaining the flares, where the VHE emission comes from background soft photons that inverse Compton scatter off of high-energy electron–positron pairs in the reconnecting current sheet. In this work, we ray trace photons from a current sheet near the black hole event horizon during a flux eruption in a magnetically arrested state in a general relativistic magnetohydrodynamics simulation. We incorporate beaming of the Compton up-scattered photons, based on results from radiative kinetic simulations of relativistic reconnection. We then construct VHE lightcurves that account for the dynamics of the current sheet and lensing from general-relativistic effects. We find that most of the flux originates in the inner 5 gravitational radii, and beaming is essential to explain the observed flux from the strongest VHE flares. The ray traced lightcurves show features resulting from the changing volume of the reconnecting current sheet on timescales that can be consistent with observations. Furthermore, we find that the amount of beaming depends strongly on two effects: (i) the current sheet inclination with respect to the observer and (ii) the anisotropy in the direction of motion of the accelerated particles. 

Abstract The merger of a black hole (BH) and a neutron star (NS) in most cases is expected to leave no material around the remnant BH; therefore, such events are often considered as sources of gravitational waves without electromagnetic counterparts. However, a bright counterpart can emerge if the NS is strongly magnetized, as its external magnetosphere can experience radiative shocks and magnetic reconnection during/after the merger. We use magnetohydrodynamic simulations in the dynamical spacetime of a merging BH–NS binary to investigate its magnetospheric dynamics. We find that compressive waves excited in the magnetosphere develop into monster shocks as they propagate outward. After swallowing the NS, the BH acquires a magnetosphere that quickly evolves into a split-monopole configuration and then undergoes an exponential decay (balding), enabled by magnetic reconnection and also assisted by the ringdown of the remnant BH. This spinning BH drags the split monopole into rotation, forming a transient pulsar-like state. It emits a striped wind if the swallowed magnetic-dipole moment is inclined to the spin axis. We predict two types of transients from this scenario: (1) a fast radio burst emitted by the shocks as they expand to large radii; and (2) an X-ray/γ-ray burst emitted by thee^±outflow heated by magnetic dissipation. 

Abstract Supermassive binary black holes in galactic centers are potential multimessenger sources in gravitational waves and electromagnetic radiation. To find such objects, isolating unique electromagnetic signatures of their accretion flow is key. With the aid of three-dimensional general-relativistic magnetohydrodynamic simulations that utilize an approximate, semianalytic, superimposed spacetime metric, we identify two such signatures for merging binaries. Both involve magnetic reconnection and are analogous to plasma processes observed in the solar corona. The first, like colliding flux tubes that can cause solar flares, involves colliding jets that form an extended reconnection layer, dissipating magnetic energy and causing the two jets to merge. The second, akin to coronal mass ejection events, involves the accretion of magnetic field lines onto both black holes; these magnetic fields then twist, inflate, and form a trailing current sheet, ultimately reconnecting and driving a hot outflow. We provide estimates for the associated electromagnetic emission for both processes, showing that they likely accelerate electrons to high energies and are promising candidates for continuous, stochastic, and/or quasi-periodic higher-energy electromagnetic emission. We also show that the accretion flows around each black hole can display features associated with the magnetically arrested state. However, simulations with black hole spins misaligned with the orbital plane and simulations with larger Bondi radii saturate at lower values of horizon-penetrating magnetic flux than standard magnetically arrested disks, leading to weaker, intermittent jets owing to feedback from the weak jets or equatorial flux tubes ejected by reconnecting field lines near the horizon. 

Abstract Using two-dimensional general relativistic resistive magnetohydrodynamic simulations, we investigate the properties of the sheath separating the black hole jet from the surrounding medium. We find that the electromagnetic power flowing through the jet sheath is comparable to the overall accretion power of the black hole. The sheath is an important site of energy dissipation as revealed by the copious appearance of reconnection layers and plasmoid chains. About 20% of the sheath power is dissipated between 2 and 10 gravitational radii. The plasma in the dissipative sheath moves along a nearly paraboloidal surface with transrelativistic bulk motions dominated by the radial component, whose dimensionless 4-velocity is ∼1.2 ± 0.5. In the frame moving with the mean (radially dependent) velocity, the distribution of stochastic bulk motions resembles a Maxwellian with an “effective bulk temperature” of ∼100 keV. Scaling the global simulation to Cygnus X-1 parameters gives a rough estimate of the Thomson optical depth across the jet sheath, ∼0.01–0.1, and it may increase in future magnetohydrodynamic simulations with self-consistent radiative losses. These properties suggest that the dissipative jet sheath may be a viable “coronal” region, capable of upscattering seed soft photons into a hard, nonthermal tail, as seen during the hard states of X-ray binaries and active galactic nuclei. 

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. 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. 

Abstract A variety of high-energy astrophysical phenomena are powered by the release—via magnetic reconnection—of the energy stored in oppositely directed fields. Single-fluid resistive magnetohydrodynamic (MHD) simulations with uniform resistivity yield dissipation rates that are much lower (by nearly 1 order of magnitude) than equivalent kinetic calculations. Reconnection-driven phenomena could be accordingly modeled in resistive MHD employing a nonuniform, “effective” resistivity informed by kinetic calculations. In this work, we analyze a suite of fully kinetic particle-in-cell (PIC) simulations of relativistic pair-plasma reconnection—where the magnetic energy is greater than the rest mass energy—for different strengths of the guide field orthogonal to the alternating component. We extract an empirical prescription for the effective resistivity, ${\eta }_{\mathrm{eff}}=\alpha B_0\mid \mathit{J}{\mid }^p/\left(\mid \mathit{J}{\mid }^{p+1}+{\left(e{n}_{t}c\right)}^{p+1}\right)$, whereB₀is the reconnecting magnetic field strength,Jis the current density,n_tis the lab-frame total number density,eis the elementary charge, andcis the speed of light. The guide field dependence is encoded inαandp, which we fit to PIC data. This resistivity formulation—which relies only on single-fluid MHD quantities—successfully reproduces the spatial structure and strength of nonideal electric fields and thus provides a promising strategy for enhancing the reconnection rate in resistive MHD simulations. 

Abstract The time-variable emission from the accretion flow of Sgr A*, the supermassive black hole at the Galactic center, has long been examined in the radio-to-millimeter, near-infrared (NIR), and X-ray regimes of the electromagnetic spectrum. However, until now, sensitivity and angular resolution have been insufficient in the crucial mid-infrared (MIR) regime. The MIRI instrument on JWST has changed that, and we report the first MIR detection of Sgr A*. The detection was during a flare that lasted about 40 minutes, a duration similar to NIR and X-ray flares, and the source's spectral index steepened as the flare ended. The steepening suggests that synchrotron cooling is an important process for Sgr A*'s variability and implies magnetic fields strengths ~ 40–70 G in the emission zone. Observations at 1.3 mm with the Submillimeter Array revealed a counterpart flare lagging the MIR flare by ≈10 minutes. The observations can be self-consistently explained as synchrotron radiation from a single population of gradually cooling high-energy electrons accelerated through (a combination of) magnetic reconnection and/or magnetized turbulence. 

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

ABSTRACT The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current microphysical CR transport models in magnetohydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying the alternative regime of large-amplitude turbulence with δB/B0 ≫ 1, in which intermittent small-scale magnetic field reversals are ubiquitous. We study particle transport in such turbulence by integrating trajectories in stationary snapshots. To quantify spatial diffusion, we use a set-up with continuous particle injection and escape, which we term the turbulent leaky box. We find that particle transport is very different from the strong guide-field case. Low-energy particles are better confined than high-energy particles, despite less efficient pitch-angle isotropization at small energies. In the limit of weak guide field, energy-dependent confinement is driven by the energy-dependent (in)ability to follow reversing magnetic field lines exactly and by the scattering in regions of ‘resonant curvature’, where the field line bends on a scale that is of the order of the local particle gyro-radius. We derive a heuristic model of particle transport in magnetic folds that approximately reproduces the energy dependence of transport found numerically. We speculate that CR propagation in the Galaxy is regulated by the intermittent field reversals highlighted here and discuss the implications of our findings for CR transport in the Milky Way. 

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