Black Hole Magnetosphere Research Articles

Scaling up global kinetic models of pulsar magnetospheres using a hybrid force-free-PIC numerical approach

Context. The particle-in-cell approach has proven effective in modeling neutron-star and black-hole magnetospheres from first principles, but global simulations are plagued with an unrealistically small separation between the scales where microphysics operates and the system-size scales due to limited numerical resources. A legitimate concern is whether the scale separation achieved to date is large enough for results to be safely extrapolated to realistic scales. Aims. In this work, our aim is to explore the effect of scaling up physical parameters and to check whether salient features uncovered by pure kinetic models at smaller scales are still valid, with a special emphasis on particle acceleration and high-energy radiation emitted beyond the light cylinder. Methods. To reach this objective, we developed a new hybrid numerical scheme coupling the ideal force-free and the particle-in-cell methods to optimize the numerical cost of global models. We propose a domain decomposition of the magnetosphere based on the magnetic-field topology using the flux function. The force-free model is enforced along open field lines while the particle-in-cell model is restricted to the reconnecting field line region. Results. As a proof of concept, this new hybrid model is applied to simulate a weak millisecond pulsar magnetosphere with realistic scales using high-resolution axisymmetric simulations. Magnetospheric features reported by previous kinetic models are recovered, and strong synchrotron radiation above 100MeV consistent with the Fermi-LAT gamma-ray pulsar population is successfully reproduced. Conclusions.This work further consolidates the shining-reconnecting current sheet scenario as the origin of the gamma-ray emission in pulsars, as well as firmly establishing pulsar magnetospheres as at least teraelectronvolt particle accelerators.

Mapping the Distribution of the Magnetic Field Strength along the NGC 315 Jet

We study magnetic field strengths along the jet in NGC 315. First, we estimated the angular velocity of rotation in the jet magnetosphere by comparing the measured velocity profile of NGC 315 with the magnetohydrodynamic jet model proposed by Tomimatsu and Takahashi. Similar to the case of M87, we find that the model can reproduce the logarithmic feature of the velocity profile and suggest a slowly rotating black hole magnetosphere for NGC 315. By substituting the estimated Ω F into the jet power predicted by the Blandford–Znajek mechanism, we estimate the magnetic field strength near the event horizon of the central black hole as 5 × 103 G ≲ B H ≲ 2 × 104 G. We then estimate magnetic field strengths along the jet by comparing the spectral index distribution obtained from very long baseline interferometry (VLBI) observations with a synchrotron-emitting jet model. Then we constrain the magnetic field strength at a deprojected distance z from the black hole to be in the range 0.06 G ≲ B(z) ≲ 0.9 G for 5.2 × 103 r g ≲ z ≲ 4.9 × 104 r g , where r g represents the gravitational radius. By combining the obtained field strengths at the event horizon and the downstream section of the jet, we find that the accretion flow at the jet base is consistent with a magnetically arrested disk. We discuss a comparison of the jet power and the magnetic flux anchored to the event horizon in NGC 315 and M87.

Current Sheet Alignment in Oblique Black Hole Magnetospheres: A Black Hole Pulsar?

We study the magnetospheric evolution of a nonaccreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of 3D general relativistic magnetohydrodynamic simulations. This serves as a model for a neutron star (NS) collapse or a BH–NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped wind. We show that the initially inclined split monopolar current sheet aligns over time with the BH equatorial plane. The inclination angle evolves exponentially toward alignment, with an alignment timescale that is inversely proportional to the square of the BH angular velocity, where higher spin results in faster alignment. Furthermore, magnetic reconnection in the current sheet leads to exponential decay of event-horizon-penetrating magnetic flux with nearly the same timescale for all considered BH spins. In addition, we present relations for the BH mass and spin in terms of the period and alignment timescale of the striped wind. The explored scenario of a rotating, aligning, and reconnecting current sheet can potentially lead to multimessenger electromagnetic counterparts to a gravitational-wave event due to the acceleration of particles powering high-energy radiation, plasmoid mergers resulting in coherent radio signals, and pulsating emission due to the initial misalignment of the BH magnetosphere.

One-dimensional General Relativistic Particle-in-cell Simulations of Stellar-mass Black Hole Magnetospheres: A Semianalytic Model of Gamma-Rays from Gaps

In the absence of a sufficient amount of plasma injection into the black hole (BH) magnetosphere, the force-free state of the magnetosphere cannot be maintained, leading to the emergence of strong, time-dependent, longitudinal electric fields (i.e., spark gaps). Recent studies of supermassive BH magnetospheres using analytical methods and particle-in-cell (PIC) simulations propose the possibility of efficient particle acceleration and consequent gamma-ray emission in the spark gap. In this work, we perform 1D general relativistic PIC simulations to examine the gamma-ray emission from stellar-mass BH magnetospheres. We find that intermittent spark gaps emerge and particles are efficiently accelerated in a similar manner to the supermassive BH case. We build a semianalytic model of the plasma dynamics and radiative processes, which reproduces the maximum electron energies and peak gamma-ray luminosities of the simulation results. Based on this model, we show that the gamma-ray signals from stellar-mass BHs wandering through the interstellar medium could be detected by gamma-ray telescopes such as the Fermi Large Area Telescope or the Cherenkov Telescope Array.

Jump-starting Relativistic Flows and the M87 Jet

We point out the dominant importance of plasma injection effects of relativistic winds from pulsars and black holes. We demonstrate that outside the light cylinder, the magnetically dominated outflows sliding along the helical magnetic field move nearly radially with very large Lorentz factors, γ 0 ≫ 1, imprinted into the flow during pair production within the gaps. Only at larger distances, r ≥ γ 0(c/Ω), does MHD acceleration Γ ∝ r take over. As a result, Blandford–Znajek (BZ)-driven outflows produce spine-brightened images. The best-resolved case of the jet in M87 shows both edge-brightened features, as well as weaker spine-brightened features. Only the spine-brightened component can be BZ driven/originate from the black hole's magnetosphere.

Pair Plasma Cascade in Rotating Black Hole Magnetospheres with a Split Monopole Flux Model

An electron-positron cascade in the magnetospheres of Kerr black holes (BHs) is a fundamental ingredient to fueling the relativistic γ-ray jets seen at the polar regions of galactic supermassive BHs (SMBHs). This leptonic cascade occurs in the spark gap region of a BH magnetosphere where the unscreened electric field parallel to the magnetic field is present; hence, it is affected by the magnetic field structure. A previous study explored the case of a thin accretion disk, representative of active galactic nuclei. Here we explore the case of a quasi-spherical gas distribution, as is expected to be present around the SMBH Sgr A* in the center of our Milky Way galaxy, for example. The properties and efficiency of the leptonic cascade are studied. The findings of our study and the implications for SMBH systems in various spectral and accretion states are discussed. The relationships and scalings derived from varying the mass of the BH and background photon spectra are further used to analyze the leptonic cascade process to power jets seen in astronomical observations. In particular, one finds the efficiency of the cascade in a quasi-spherical gas distribution peaks at the jet axis. Observationally, this should lead to a more prominent jet core, in contrast to the thin disk accretion case, where it peaks around the jet–disk interface. One also finds the spectrum of the background photons plays a key role. The cascade efficiency is maximum for a spectral index of 2, while harder and softer spectra lead to a less efficient cascade.

Dynamics of M87 jet

AbstractFlows originating from black hole magnetospheres via Blandford-Znajek (BZ) process start highly relativistically, with very large Lorentz factors γ0 1, imprinted into the flow during pair production within the gaps. As a result, BZ-driven outflows would produce spine-brightened images, contrary to observations of the edge-brightened jet in M87. We conclude that M87 jet is not BZ-driven.

Reconnection-driven flares in 3D black hole magnetospheres

Context. Low-luminosity supermassive and stellar-mass black holes (BHs) may be embedded in a collisionless and highly magnetized plasma. They show nonthermal flares indicative of particles being accelerated up to relativistic speeds by dissipative processes in the vicinity of the BH. During near-infrared flares from the supermassive BH Sagittarius A* (Sgr A*), the GRAVITY Collaboration detected circular motion and polarization evolution, which suggest the presence of transient synchrotron-emitting hot spots moving around the BH. Aims. We study 3D reconnecting current layers in the magnetosphere of spinning BHs to determine whether plasma-loaded flux ropes which are formed near the event horizon could reproduce the hot spot observations and help constrain the BH spin. Methods. We performed global 3D particle-in-cell simulations in Kerr spacetime of a pair plasma embedded in a strong and large-scale magnetic field originating in a perfectly conducting disk in prograde Keplerian rotation. Results. A cone-shaped current layer develops which surrounds the twisted open magnetic field lines threading the event horizon. Spinning magnetic field lines coupling the disk to the BH inflate and reconnect a few gravitational radii above the disk. This quasi-periodic cycle accelerates particles, which accumulate in a few macroscopic flux ropes rotating with the outermost coupling magnetic field line. Once flux ropes detach, they propagate in the current layer following what appears as a rapidly opening spiral when seen face-on. A single flux rope carries enough relativistic electrons and positrons to emit synchrotron radiation at levels suitable to reproduce the peak-luminosity of the flares of Sgr A* but it quickly fades away as it flows away. Conclusions. Our kinematic analysis of the flux ropes motion favors a BH spin of 0.65 to 0.8 for Sgr A*. The duration of the flares of Sgr A* can only be explained provided the underlying magnetic loop seeded in the disk mid-plane has a finite lifetime and azimuthal extension. In this scenario, the hot spot corresponds to a spinning arc along which multiple reconnection sites power the net emission as flux ropes episodically detach.

Linear analysis of the Kelvin–Helmholtz instability in relativistic magnetized symmetric flows

ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$, where $v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.

Magnetic reconnection plasmoid model for Sagittarius A* flares

Context. Sagittarius A*, the supermassive black hole at the center of our Galaxy, exhibits episodic near-infrared flares. The recent monitoring of three such events with the GRAVITY instrument has shown that some flares are associated with orbital motions in the close environment of the black hole. The GRAVITY data analysis indicates a super-Keplerian azimuthal velocity, while (sub-) Keplerian velocity is expected for the hot flow surrounding the black hole. Aims. We develop a semi-analytic model of the Sagittarius A* flares based on an ejected large plasmoid, inspired by recent particle-in-cell global simulations of black hole magnetospheres. We model the infrared astrometric and photometric signatures associated with this model. Methods. We considered a spherical macroscopic hot plasma region that we call a large plasmoid. This structure was ejected along a conical orbit in the vicinity of the black hole. This plasmoid was assumed to be formed by successive mergers of smaller plasmoids produced through magnetic reconnection that we did not model. Nonthermal electrons were injected into the plasmoid. We computed the evolution of the electron-distribution function under the influence of synchrotron cooling. We solved the radiative transfer problem associated with this scenario and transported the radiation along null geodesics of the Schwarzschild space time. We also took the quiescent radiation of the accretion flow into account, on top of which the flare evolves. Results. For the first time, we successfully account for the astrometric and flux variations of the GRAVITY data with a flare model that incorporates an explicit modeling of the emission mechanism. The prediction of our model and recent data agree well. In particular, the azimuthal velocity of the plasmoid is set by the magnetic field line to which it belongs, which is anchored in the inner parts of the accretion flow, hence the super-Keplerian motion. The astrometric track is also shifted with respect to the center of mass due to the quiescent radiation, in agreement with the difference measured with the GRAVITY data. Conclusions. These results support the hypothesis that magnetic reconnection in a black hole magnetosphere is a viable model for the infrared flares of Sagittarius A*.

Energy-dependent flavor ratios, cascade/track spectrum tension and high-energy neutrinos from magnetospheres of supermassive black holes

The IceCube neutrino observatory measures the diffuse flux of high-energy astrophysical neutrinos by means of various techniques, and there exists a mild tension between spectra obtained in different analyses. The spectrum derived from reconstruction of muon tracks is harder than that from cascades, dominated by electron and tau neutrinos. If confirmed, this tension may provide a clue to the origin of these neutrinos, which remains uncertain. Here we investigate the possibility that this tension may be caused by the change of the flavor content of astrophysical neutrinos with energy. We assume that at higher energies, the flux contains more muon neutrinos than expected in the usually assumed flavor equipartition. This may happen if the neutrinos are produced in regions of the magnetic field so strong that muons, born in pi-meson decays, cool by synchrotron radiation faster than decay. The magnetic field of ∼104 G is required for this mechanism to be relevant for the IceCube results. We note that these field values are reachable in the immediate vicinity of supermassive black holes in active galactic nuclei and present a working toy model of the population of these potential neutrino sources. While this model predicts the required flavor ratios and describes the high-energy spectrum, it needs an additional component to explain the observed neutrino flux at lower energies.

Radiative Reconnection-powered TeV Flares from the Black Hole Magnetosphere in M87

Active galactic nuclei in general, and the supermassive black hole in M87 in particular, show bright and rapid gamma-ray flares up to energies of 100 GeV and above. For M87, the flares show multiwavelength components, and the variability timescale is comparable to the dynamical time of the event horizon, suggesting that the emission may come from a compact region near the nucleus. However, the emission mechanism for these flares is not well understood. Recent high-resolution general-relativistic magnetohydrodynamic simulations show the occurrence of episodic magnetic reconnection events that can power flares near the black hole event horizon. In this work, we analyze the radiative properties of the reconnecting current layer under the extreme plasma conditions applicable to the black hole in M87 from first principles. We show that abundant pair production is expected in the vicinity of the reconnection layer, to the extent that the produced secondary pair plasma dominates the reconnection dynamics. Using analytic estimates backed by two-dimensional particle-in-cell simulations, we demonstrate that in the presence of strong synchrotron cooling, reconnection can produce a hard power-law distribution of pair plasma imprinted in the outgoing synchrotron (up to a few tens of MeV) and the inverse-Compton signal (up to TeV). We produce synthetic radiation spectra from our simulations, which can be directly compared with the results of future multiwavelength observations of M87* flares.

Synthetic Images of Magnetospheric Reconnection-Powered Radiation around Supermassive Black Holes.

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

Implications from the Velocity Profile of the M87 Jet: A Possibility of a Slowly Rotating Black Hole Magnetosphere

Motivated by the measured velocity profile of the M87 jet using the Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) by Park et al. indicating that the starting position of the jet acceleration is farther from the central engine of the jet than predicted in general relativistic magnetohydrodynamic simulations, we explore how to mitigate the apparent discrepancy between the simulations and the KaVA observation. We use a semi-analytic jet model proposed by Tomimatsu & Takahashi consistently solving the transmagnetic field structure but neglecting any dissipation effects. By comparing the jet model with the observed M87 jet velocity profile, we find that the model can reproduce the logarithmic feature of the velocity profile and can fit the observed data when choosing c/(100r g ) ≲ Ω F ≲ c/(70r g ), where r g is the gravitational radius. While a total specific energy () of the jet changes the terminal bulk Lorentz factor of the jet, a slower angular velocity of the black hole magnetosphere (Ω F ) makes a light-cylinder radius (r lc) larger, and it consequently pushes out a location of a starting point of the jet acceleration. Using the estimated Ω F we further estimate the magnetic field strength on the event horizon scale in M87 by assuming Blandford–Znajek process is in action. The corresponding magnetic flux threading the event horizon of M87 is in good agreement with a magnetically arrested disk regime.

Magnetic Reconnection in Black Hole Magnetospheres: Lepton Loading into Jets, Superluminal Radio Blobs, and Multiwavelength Flares

Supermassive black holes in active galactic nuclei launch relativistic jets, as indicated by observed superluminal radio blobs. The energy source of these jets is widely discussed in the theoretical framework of the Blandford–Znajek process, the electromagnetic energy extraction from rotating black holes (BHs), while the formation mechanism of the radio blobs in the electromagnetically dominated jets has been a long-standing problem. Recent high-resolution magnetohydrodynamic simulations of magnetically arrested disks exhibited magnetic reconnection in a transient magnetically dominated part of the equatorial disk near the BH horizon, which led to a promising scenario of efficient MeV gamma-ray production and subsequent electron–positron pair loading into the BH magnetosphere. We develop this scenario to build a theoretical framework on energetics, timescales, and particle number density of the superluminal radio blobs and discuss observable signatures in other wave bands. We analytically show that the nonthermal electrons emit broadband photons from optical to multi-MeV bands. The electron–positron pairs produced in the magnetosphere are optically thick for synchrotron self-absorption, so that the injected energy is stored in the plasma. The stored energy is enough to power the superluminal radio blobs observed in M87. This scenario predicts rather dim radio blobs around Sgr A*, which are consistent with no clear detection by current facilities. In addition, this scenario inevitably produces strong X-ray flares in a short timescale, which will be detectable by future X-ray satellites.

Chiral anomalous processes in magnetospheres of pulsars and black holes

We propose that chirally asymmetric plasma can be produced in the gap regions of the magnetospheres of pulsars and black holes. We show that, in the case of supermassive black holes situated in active galactic nuclei, the chiral charge density and the chiral chemical potential are very small and unlikely to have any observable effects. In contrast, the chiral asymmetry produced in the magnetospheres of magnetars can be substantial. It can trigger the chiral plasma instability that, in turn, can lead to observable phenomena in magnetars. In particular, the instability should trigger circularly polarized electromagnetic radiation in a wide window of frequencies, spanning from radio to near-infrared. As such, the produced chiral charge has the potential to affect some features of fast radio bursts.

Acceleration and cooling of the corona during X-ray flares from the Seyfert galaxy I Zw 1

ABSTRACTWe report on X-ray flares that were observed from the active galactic nucleus I Zwicky 1 (I Zw 1) in 2020 January by the NuSTAR and XMM-Newton observatories. The X-ray spectrum is well-described by a model comprised of the continuum emission from the corona and its reflection from the accretion disc around a rapidly spinning (a > 0.94) black hole. In order to model the broadband spectrum, it is necessary to account for the variation in ionization across the disc. Analysis of the X-ray spectrum in time periods before, during and after the flares reveals the underlying changes to the corona associated with the flaring. During the flares, the reflection fraction drops significantly, consistent with the acceleration of the corona away from the accretion disc. We find the first evidence that during the X-ray flares, the temperature drops from $140_{-20}^{+100}$ keV before to $45_{-9}^{+40}$ keV during the flares. The profile of the iron K line reveals the emissivity profile of the accretion disc, showing it to be illuminated by a compact corona extending no more than $7_{-2}^{+4}$ rg over the disc before the flares, but with tentative evidence that the corona expands as it is accelerated during the flares. Once the flares subsided, the corona had collapsed to a radius of $6_{-2}^{+2}$ rg. The rapid time-scale of the flares suggests that they arise within the black hole magnetosphere rather than in the accretion disc, and the variation of the corona is consistent with the continuum arising from the Comptonization of seed photons from the disc.

The Response of Black Hole Spark Gaps to External Changes: A Production Mechanism of Rapid TeV Flares?

Abstract We study the response of a starved Kerr black hole magnetosphere to abrupt changes in the intensity of disk emission and in the global magnetospheric current, by means of one-dimensional general relativistic particle-in-cell simulations. Such changes likely arise from the intermittency of the accretion process. We find that in cases where the pair-production opacity contributed by the soft disk photons is modest, as in, e.g., M87, such changes can give rise to delayed, strong teraelectronvolt (TeV) flares, dominated by curvature emission of particles accelerated in the gap. The flare rise time, and the delay between the external variation and the onset of the flare emitted from the outer gap boundary, are of the order of the light-crossing time of the gap. The rapid, large-amplitude TeV flares observed in M87, and perhaps, other active galactic nuclei may be produced by such a mechanism.

Electrically charged black holes and the Blandford–Znajek mechanism

ABSTRACT Recently, it was claimed by King & Pringle that accretion of electric charge by a black hole (BH) rotating in an aligned external magnetic field results in a ‘dead’ vacuum magnetosphere, where the electric field is totally screened, no vacuum breakdown is possible, and the Blandford–Znajek mechanism cannot operate. Here, we study in details the properties of the Wald solution for electrically charged BHs discussed in their paper. Our results show that the claim is erroneous as in the solution with the critical charge q0 = 2aB0, there exists a drop of electrostatic potential along all magnetic field lines except the one coinciding with the symmetry axis. It is also found that while uncharged rotating BHs expel external vacuum magnetic field from their event horizon (the Meissner effect), electric charging of BHs pulls the magnetic field lines back on it, resembling what has been observed in some previous force-free, Relativistic Magnetohydrodynamics (RMHD), and Particle-in-Cell (PIC) simulations of BH magnetospheres. This suggests that accretion of electric charge may indeed be a feature of the BH electrodynamics. However, our analysis shows that the value q0 of the BH charge given by Wald is likely to be only an upper limit, and that the actual value depends of the details of the magnetospheric physics.

Magnetic Hair and Reconnection in Black Hole Magnetospheres.

The no-hair theorem of general relativity states that isolated black holes are characterized by three parameters: mass, spin, and charge. In this Letter we consider Kerr black holes endowed with highly magnetized plasma-filled magnetospheres. Using general relativistic kinetic plasma and resistive magnetohydrodynamics simulations, we show that a dipole magnetic field on the event horizon opens into a split monopole and reconnects in a plasmoid-unstable current sheet. The no-hair theorem is satisfied, in the sense that all components of the stress-energy tensor decay exponentially in time. We measure the decay time of magnetic flux on the event horizon for plasmoid-dominated reconnection in collisionless and collisional plasma. The reconnecting magnetosphere should be a powerful source of hard x-ray emission when the magnetic field is strong.