Black Hole Magnetosphere Research Articles

Enhanced gamma radiation towards the rotation axis from the immediate vicinity of extremely rotating black holes

Abstract We investigate the acceleration of electrons and positrons by magnetic-field-aligned electric fields in the polar funnel of an accreting black hole (BH). Applying the pulsar outer-gap theory to BH magnetospheres, we find that such a lepton accelerator arises in the immediate vicinity of the event horizon due to frame-dragging, and that their gamma-ray luminosity increases with decreasing accretion rate. Furthermore, we demonstrate that the gamma-ray flux is enhanced along the rotation axis by more than an order of magnitude if the BH spin increases from a = 0.90M to a = 0.9999M. As a result, if a ten-solar-mass, almost-maximally rotating BH is located within 3 kpc, when its accretion rate is between 0.005 and 0.01 per cent of the Eddington rate, its high-energy flare becomes detectable with the Fermi/Large Area Telescope, provided that the flare lasts longer than 1.2 months and that we view the source nearly along the rotation axis. In addition, its very high energy flux is marginally detectable with the Cherenkov Telescope Array, provided that the flare lasts longer than a night and that our viewing angle is about 45 deg with respect to the rotation axis.

Strong constraints on hadronic models of blazar activity fromFermiand IceCube stacking analysis

High-energy emission from blazars is produced by electrons which are either accelerated directly (the assumption of leptonic models of blazar activity) or produced in interactions of accelerated protons with matter and radiation fields (the assumption of hadronic models). The hadronic models predict that gamma-ray emission is accompanied by neutrino emission with comparable energy flux but with a different spectrum. We derive constraints on the hadronic models of activity of blazars imposed by non-detection of neutrino flux from a population of gamma-ray emitting blazars. We stack the gamma-ray and muon neutrino flux from 749 blazars situated in the declination strip above -5 degrees. Non-detection of neutrino flux from the stacked blazar sample rules out the proton induced cacade models in which the high-energy emission is powered by interactions of shock-accelerated proton beam in the AGN jet with the ambient matter or with the radiation field of the black hole accretion disk. The result remains valid also for the case of interactions in the scattered radiation field in the broad line region. IceCube constraint could be avoided if the spectrum of accelerated protons is sharply peaking in the ultra-high-energy cosmic ray range, as in the models of acceleration in the magnetic reconnection regions or in the vacuum gaps of black hole magnetospheres. Models based on these acceleration mechanisms are consistent with the data only if characteristic energies of accelerated protons are higher than 1e19 eV.

Scenarios for Ultrafast Gamma-Ray Variability in AGN

Abstract We analyze three scenarios to address the challenge of ultrafast gamma-ray variability reported from active galactic nuclei. We focus on the energy requirements imposed by these scenarios: (i) external cloud in the jet, (ii) relativistic blob propagating through the jet material, and (iii) production of high-energy gamma-rays in the magnetosphere gaps. We show that while the first two scenarios are not constrained by the flare luminosity, there is a robust upper limit on the luminosity of flares generated in the black hole magnetosphere. This limit depends weakly on the mass of the central black hole and is determined by the accretion disk magnetization, viewing angle, and the pair multiplicity. For the most favorable values of these parameters, the luminosity for 5-minute flares is limited by , which excludes a black hole magnetosphere origin of the flare detected from IC 310. In the scopes of scenarios (i) and (ii), the jet power, which is required to explain the IC 310 flare, exceeds the jet power estimated based on the radio data. To resolve this discrepancy in the framework of scenario (ii), it is sufficient to assume that the relativistic blobs are not distributed isotropically in the jet reference frame. A realization of scenario (i) demands that the jet power during the flare exceeds by a factor 102 the power of the radio jet relevant to a timescale of 108 years.

Effects of inner Alfvén surface location on black hole energy extraction in the limit of a force-free magnetosphere

An energy extracting black hole magnetosphere can be defined by the location of its inner Alfv\'{e}n surface, which determines the rate of black hole energy extraction along a given magnetic field line. We study how the location of the inner Alfv\'{e}n surface can modify the structure of energy extracting black hole magnetospheres in the force-free limit. Hundreds of magnetospheres are numerically computed via a general relativistic extension of the Newtonian magnetofrictional method for a full range of black hole spins and flow parameters. We find that jet-like structures naturally form very close to the horizon for Alfv\'{e}n surfaces near the boundary of the ergosphere and that energy is extracted towards the equatorial plane for Alfv\'{e}n surfaces close to the horizon. This suggests that two broad classes of energy extracting black hole magnetospheres might exist; those that transmit extracted energy directly to distant observers, and those that transmit extracted energy to nearby accreting matter. Applied to transient high energy phenomena, we find that they might also differ in timescale by a factor of 20 or more.

New mechanism of acceleration of particles by stellar black holes

In this paper, we study the efficiency of particle acceleration in the magnetospheres of stellar mass black holes. For this purpose, we consider the linearized set of the Euler equation, continuity equation and Poisson equation, respectively. After introducing the varying relativistic centrifugal force, we show that the charge separation undergoes the parametric instability, leading to generation of centrifugally-excited Langmuir waves. It is shown that these waves, via the Langmuir collapse damp by means of the Landau damping, as a result, energy transfers to particles accelerating them to energies of the order of [Formula: see text] eV.

LEPTON ACCELERATION IN THE VICINITY OF THE EVENT HORIZON: HIGH-ENERGY AND VERY-HIGH-ENERGY EMISSIONS FROM ROTATING BLACK HOLES WITH VARIOUS MASSES

ABSTRACT We investigate the electrostatic acceleration of electrons and positrons in the vicinity of the event horizon, applying the pulsar outer-gap model to black hole (BH) magnetospheres. During a low accretion phase, the radiatively inefficient accretion flow (RIAF) cannot emit enough MeV photons that are needed to sustain the force-free magnetosphere via two-photon collisions. In such a charge-starved region (or a gap), an electric field arises along the magnetic field lines to accelerate electrons and positrons into ultra-relativistic energies. These relativistic leptons emit copious gamma rays via curvature and inverse-Compton (IC) processes. Some of such gamma rays collide with the submillimeter-IR photons emitted from the RIAF to materialize as pairs, which polarize to partially screen the original acceleration electric field. It is found that the gap gamma-ray luminosity increases with decreasing accretion rate. However, if the accretion rate decreases too much, the diminished RIAF soft photon field can no longer sustain a stationary pair production within the gap. As long as a stationary gap is formed, the magnetosphere becomes force-free outside the gap by the cascaded pairs, irrespective of the BH mass. If a nearby stellar-mass BH is in quiescence, or if a galactic intermediate-mass BH is in a very low accretion state, its curvature and IC emissions are found to be detectable with Fermi/LAT and imaging atmospheric Cherenkov telescopes (IACT). If a low-luminosity active galactic nucleus is located within about 30 Mpc, the IC emission from its supermassive BH is marginally detectable with IACT.

Riemann solvers and Alfven waves in black hole magnetospheres.

In the magnetosphere of a rotating black hole, an inner Alfven critical surface (IACS) must be crossed by inflowing plasma. Inside the IACS, Alfven waves are inward directed toward the black hole. The majority of the proper volume of the active region of spacetime (the ergosphere) is inside of the IACS. The charge and the totally transverse momentum flux (the momentum flux transverse to both the wave normal and the unperturbed magnetic field) are both determined exclusively by the Alfven polarization. Thus, it is important for numerical simulations of black hole magnetospheres to minimize the dissipation of Alfven waves. Elements of the dissipated wave emerge in adjacent cells regardless of the IACS, there is no mechanism to prevent Alfvenic information from crossing outward. Thus, numerical dissipation can affect how simulated magnetospheres attain the substantial Goldreich-Julian charge density associated with the rotating magnetic field. In order to help minimize dissipation of Alfven waves in relativistic numerical simulations we have formulated a one-dimensional Riemann solver, called HLLI, which incorporates the Alfven discontinuity and the contact discontinuity. We have also formulated a multidimensional Riemann solver, called MuSIC, that enables low dissipation propagation of Alfven waves in multiple dimensions. The importance of higher order schemes in lowering the numerical dissipation of Alfven waves is also catalogued.

Ultrafast VHE Gamma-Ray Flares of IC 310

AbstractIn 2012 November MAGIC detected a bright flare from IC 310. The flare consisted of two sharp peaks with a typical duration of ~ 5 min. The energy released during that event has been estimated to be at the level of 2 × 1044 erg s−1.In this work we derive an upper limit on the possible luminosity of flares generated in black hole (BH) magnetosphere, which depends very weakly on the mass of BH and is determined by disk magnetisation, viewing angle, and pair multiplicity. Since all these parameters are smaller than a unit, the luminosity 2 × 1043 erg s−1 can be taken as a strict upper limit for flare luminosity for several minutes variability time. This upper limit appears to be approximately an order of magnitude below the value measured with MAGIC. Thus, we conclude that it seems very unfeasible that the magnetospheric processes can be indeed behind the bright flaring activity recorded from IC 310.

Particle acceleration in the vacuum gaps in black hole magnetospheres

We consider particle acceleration in vacuum gaps in magnetospheres of black holes powered through Blandford-Znajek mechanism and embedded into radiatively-inefficient accretion flow (RIAF) environment. In such situation the gap height is limited by the onset of gamma-gamma pair production on the infrared photons originating from the RIAF. We numerically calculate acceleration and propagation of charged particles taking into account the detailed structure of electric and magnetic field in the gap and in the entire black hole magnetosphere, radiative energy losses and interactions of gamma rays produced by the propagated charged particles with the background radiation field of RIAF. We show that the presence of the vacuum gap has clear observational signatures. The spectra of emission from gaps embedded into a relatively high luminosity RIAF are dominated by the inverse Compton emission with a sharp, super-exponential cut-off in the very-high-energy gamma-ray band. The cut-off energy is determined by the properties of the RIAF and is largely independent of the structure of magnetosphere and geometry of the gap.The spectra of the gap residing in low-luminosity RIAFs are dominated by synchrotron / curvature emission with the spectra extending into 1-100 GeV energy range. We also consider the effect of possible acceleration of protons in the gap and find that proton energies could reach the ultra-high-energy cosmic ray (UHECR) range only in extremely low luminosity RIAFs.

ANALYTIC PROPERTIES OF FORCE-FREE JETS IN THE KERR SPACETIME—II

ABSTRACT We reinvestigate the structure of a steady axisymmetric force-free magnetosphere around a Kerr black hole (BH). The BH magnetosphere structure is governed by a second-order differential equation of A ϕ depending on two “free” functions Ω and I, where A ϕ is the ϕ component of the vector potential of the electromagnetic field, Ω is the angular velocity of the magnetic field lines, and I is the poloidal electric current. While the two functions Ω and I are not arbitrarily given, they do need to be self-consistently determined along with the differential equation. Based on the perturbation approach we proposed in paper I, in this paper, we self-consistently sort out two boundary conditions governing Ω and I, and interpret these conditions mathematically and physically. Making use of the boundary conditions, we prove that all magnetic field lines crossing the infinite-redshift surface also penetrate the event horizon. Furthermore, we argue that the BH Meissner effect does not work in a force-free magnetosphere due to the perfect conductivity.

ANALYTIC PROPERTIES OF FORCE-FREE JETS IN THE KERR SPACETIME—I

The Blandford–Znajek (BZ) mechanism describes a process for extracting the rotation energy from a spinning black hole (BH) via magnetic field lines penetrating the event horizon of the central BH. In this paper, we present a perturbation approach to study force-free jets launched by the BZ mechanism, and then discuss its two immediate applications: (1) we present a high-order split-monopole perturbation solution to the BZ mechanism which accurately pins down the energy extraction rate and well describes the structure of the BH magnetosphere for the entire range of BH spins (0 ≤ a ≤ 1); (2) the approach yields an exact constraint for the monopole field configuration in Kerr spacetime, where Aϕ is the ϕ component of the vector potential of electromagnetic field, Ω is the angular velocity of the magnetic field lines, and I is the poloidal electric current. The constraint is of particular importance to benchmark the accuracy of numerical simulations.

Black hole jet power from impedance matching

Black hole jet power depends on the angular velocity of magnetic field lines, $\Omega_F$. Force-free black hole magnetospheres typically have $\Omega_F/\Omega_H \approx 0.5$, where $\Omega_H$ is the angular velocity of the horizon. We give a streamlined proof of this result using an extension of the classical black hole membrane paradigm. The proof is based on an impedance-matching argument between membranes at the horizon and infinity. Then we consider a general relativistic magnetohydrodynamic simulation of an accreting, spinning black hole and jet. We find that the theory correctly describes the simulation in the jet region. However, the field lines threading the horizon near the equator have much smaller $\Omega_F/\Omega_H$ because the force-free approximation breaks down in the accretion flow.

AXISYMMETRIC, NONSTATIONARY BLACK HOLE MAGNETOSPHERES: REVISITED

An axisymmetric, stationary, general-relativistic, electrodynamic engine model of an active galactic nucleus was formulated by Macdonald and Thorne that consisted of a supermassive black hole surrounded by a plasma magnetosphere and a magnetized accretion disk. Based on this initial formulation, a nonstationary, force-free version of their model was constructed by Park & Vishniac (PV), with the simplifying assumption that the poloidal component of the magnetic field line velocity be confined along the radial direction in cylindrical polar coordinates. In this paper, we derive the new, nonstationary Transfield Equation, which was not specified in PV. If we can solve this Transfield Equation numerically, then we will understand the axisymmetric, nonstationary black hole magnetosphere in more rigorous ways.

Publisher’s Note: Note on bunching of field lines in black hole magnetospheres [Phys. Rev. D92, 044053 (2015)

Received 14 September 2015DOI:https://doi.org/10.1103/PhysRevD.92.069902© 2015 American Physical Society

Note on bunching of field lines in black hole magnetospheres

Numerical simulations of Blandford-Znajek energy extraction at high spin have revealed that field lines tend to bunch near the poles of the event horizon. We show that this behavior can be derived analytically from the assumption of fixed functional dependence of current and field line rotation on magnetic flux. The argument relies crucially on the existence of the Znajek condition, which offers non-trivial information about the fields on the horizon without requiring a full force-free solution. We also provide some new analytic expressions for the parabolic field configuration.

Spectrum for the electric dipole which nonradially falling into a black hole

In this paper, the electromagnetic bremsstrahlung spectrum of the dipole which falls along a spiral orbit into the Schwarzschild black hole was found. The characteristic features of this electromagnetic spectrum can be used for determination of the black hole mass. This new way (if implemented) provides higher accuracy of the determination of the black hole mass. Also these features in the spectrum can be used for the determination of the some physical processes in the black hole magnetosphere and in the accretion disk. It is also shown that the asymptotic behavior of this spectrum (at high frequencies) is practically independent from the impact parameter of the falling dipole.

Multi-epoch study of the gamma-ray emission within the M87 magnetosphere model

M87 is a nearby radio galaxy that has been detected at energies ranging from radio to very high energy (VHE) gamma-rays. Its proximity and its jet, misaligned from the line of sight allow detailed morphological studies. The imaging atmospheric Cherenkov technique (from 100 GeV to 10 TeV) provides insufficient angular resolution (few arc-minutes) to resolve the M87 emission region. However, the short time scale variability observed by MAGIC, HESS and VERITAS suggests the TeV emission is coming from a very small region, most likely close to the core. We propose that the variable TeV emission may be produced in a pair-starved region of the central black hole (BH) magnetosphere, i.e. a region where the density of the electron-positron plasma is not sufficient to completely screen the accelerating electric field. The funnel, a low density and magnetically dominated region around the poles, appears as a favourable site of low-density where a Blandford-Znajek process may explain the main properties of the TeV γ-ray emission from M87. We produce a broadband spectral energy distribution (SED) of the resulting radiation and compare the model with the observed fluxes from the nucleus of M87, for both low and high γ-ray activities. We finish with a brief discussion on the connection between the accretion rate and the intermittence in the formation of gaps in the magnetosphere.

Fourth-order split monopole perturbation solutions to the Blandford-Znajek mechanism

The Blandford-Znajek (BZ) mechanism describes a physical process for the energy extraction from a spinning black hole (BH), which is believed to power a great variety of astrophysical sources, such as active galactic nuclei and gamma ray bursts. The only known analytic solution to the BZ mechanism is a split monopole perturbation solution up to $O({a}^{2})$, where $a$ is the spin parameter of a Kerr black hole. In this paper, we extend the monopole solution to higher order $\ensuremath{\sim}O({a}^{4})$. We carefully investigate the structure of the BH magnetosphere, including the angular velocity of magnetic field lines $\mathrm{\ensuremath{\Omega}}$, the toroidal magnetic field ${B}^{\ensuremath{\phi}}$, as well as the poloidal electric current $I$. In addition, the relevant energy extraction rate $\stackrel{\ifmmode \dot{}\else \textperiodcentered \fi{}}{E}$ and the stability of this high-order monopole perturbation solution are also examined.

Regular and chaotic motion in general relativity: The case of an inclined black hole magnetosphere

Dynamics of charged particles in the vicinity of a rotating black hole embedded in the external large-scale magnetic field is numerically investigated. In particular, we consider a non- axisymmetric model in which the asymptotically uniform magnetic field is inclined with respect to the axis of rotation. We study the effect of inclination onto the prevailing dynamic regime of particle motion, i.e. we ask whether the inclined field allows regular trajectories or if instead, the deterministic chaos dominates the motion. In this contribution we further discuss the role of initial condition, particularly, the initial azimuthal angle. To characterize the measure of chaoticness we compute maximal Lyapunov exponents and employ the method of Recurrence Quantification Analysis.

Black hole jets without large-scale net magnetic flux

Abstract We propose a scenario for launching relativistic jets from rotating black holes, in which small-scale magnetic flux loops, sustained by disc turbulence, are forced to inflate and open by differential rotation between the black hole and the accretion flow. This mechanism does not require a large-scale net magnetic flux in the accreting plasma. Estimates suggest that the process could operate effectively in many systems, and particularly naturally and efficiently when the accretion flow is retrograde. We present the results of general-relativistic force-free electrodynamic simulations demonstrating the time evolution of the black hole's magnetosphere, the cyclic formation of jets, and the effect of magnetic reconnection. The jets are highly variable on time-scales ∼10–103rg/c, where rg is the black hole's gravitational radius. The reconnecting current sheets observed in the simulations may be responsible for the hard X-ray emission from accreting black holes.