Black Hole Jet Research Articles

Coport: a new public code for polarized radiative transfer in a covariant framework

General relativistic radiative transfer calculations are essential for comparing theoretical models of black hole accretion flows and jets with observational data. In this work, we introduce Coport,[The public version of Coport is available at the following URL: https://github.com/JieweiHuang/Coport.] a novel public code specifically designed for covariant polarized ray-tracing radiative transfer computations in any spacetime. Written in Julia, Coport includes an interface for visualizing numerical results obtained from HARM, a publicly available implementation of the general relativistic magnetohydrodynamics code. We validate the precision of our code by comparing its outputs with the results from a variety of established methodologies. This includes the verification against analytical solutions, the validation through thin-disk assessments, and the evaluation via thick-disk analyses. Notably, our code employs a methodology that eliminates the need for separating the computations of spacetime propagation and plasma propagation. Instead, it directly solves the coupled, covariant, polarized radiative transfer equation in curved spacetime, seamlessly integrating the effects of gravity with plasma influences. This approach sets our code apart from the existing alternatives and enhances its accuracy and efficiency.

Black hole jets on the scale of the cosmic web.

When sustained for megayears (refs. 1,2), high-power jets from supermassive black holes (SMBHs) become the largest galaxy-made structures in the Universe3. By pumping electrons, atomic nuclei and magnetic fields into the intergalactic medium (IGM), these energetic flows affect the distribution of matter and magnetism in the cosmic web4-6 and could have a sweeping cosmological influence if they reached far at early epochs. For the past 50 years, the known size range of black hole jet pairs ended at 4.6-5.0 Mpc (refs. 7-9), or 20-30% of a cosmic void radius in the Local Universe10. An observational lack of longer jets, as well as theoretical results11, thus suggested a growth limit at about 5 Mpc (ref. 12). Here we report observations of a radio structure spanning about 7 Mpc, or roughly 66% of a coeval cosmic void radius, apparently generated by a black hole between and 6.3 Gyr after the Big Bang. The structure consists of a northern lobe, a northern jet, a core, a southern jet with an inner hotspot and a southern outer hotspot with a backflow. This system demonstrates that jets can avoid destruction by magnetohydrodynamical instabilities over cosmological distances, even at epochs when the Universe was 7 to times denser than it is today. How jets can retain such long-lived coherence is unknown at present.

Multifrequency Analysis of Favored Models for the Messier 87* Accretion Flow

The polarized images of the supermassive black hole Messier 87* (M87*) produced by the Event Horizon Telescope (EHT) provide a direct view of the near-horizon emission from a black hole accretion and jet system. The EHT theoretical analysis of the polarized M87* images compared thousands of snapshots from numerical models with a variety of spins, magnetization states, viewing inclinations, and electron energy distributions, and found a small subset consistent with the observed image. In this article, we examine two models favored by EHT analyses: a magnetically arrested disk with moderate retrograde spin and a magnetically arrested disk with high prograde spin. Both have electron distribution functions that lead to strong depolarization by cold electrons. We ray trace five snapshots from each model at 22, 43, 86, 230, 345, and 690 GHz to forecast future very long baseline interferometry (VLBI) observations and examine limitations in numerical models. We find that even at low frequencies where optical and Faraday rotation depths are large, approximately rotationally symmetric polarization persists, suggesting that shallow depths dominate the polarization signal. However, morphology and spectra suggest that the assumed thermal electron distribution is not adequate to describe emission from the jet. We find 86 GHz images show a ringlike shape determined by a combination of plasma and spacetime imprints, smaller in diameter than recent results from the Global mm-VLBI Array. We find that the photon ring becomes more apparent with increasing frequency, and is more apparent in the retrograde model, leading to large differences between models in asymmetry and polarization structure.

Numerical relativity simulations of black hole and relativistic jet formation

ABSTRACT We investigate impacts of stellar rotation and magnetic fields on black hole (BH) formation and its subsequent explosive activities, by conducting axisymmetric radiation-magnetohydrodynamics simulations of gravitational collapse of a 70 $\mathrm{M}_\odot$ star with two-moment multi energy neutrino transport in full general relativity for the first time. Due to its dense stellar structure, all models cannot avoid the eventual BH formation even though a strongly magnetized model experiences the so-called magnetorotational explosion prior to the BH formation. One intriguing phenomenon observed in the strongly magnetized model is the formation of a relativistic jet in the post-BH formation. The relativistic jet is the outcome of a combination of strong magnetic fields and low-density materials above the BH. The jet further enhances the explosion energy beyond $\sim 10^{52}$ erg, which is well exceeding the gravitational overburden ahead of the shock. Our self-consistent supernova models demonstrate that rotating magnetized massive stars at the high-mass end of supernova progenitors could be a potential candidate of hypernova and long gamma-ray burst progenitors.

Imaging the event horizon of M87* from space on different timescales

Imaging the event horizon of M87* from space on different timescales

Two Repeated Transient Quasiperiodic Oscillations in the γ-Ray Emission from the Blazar 3C 279

In this work, we report for the first time two repeated transient quasiperiodic oscillations (QPOs) in the γ-ray light curve of the TeV blazar 3C 279. We search for the periodicity in the light curve and estimate its confidence level using the weighted wavelet Z-transform, the Lomb–Scargle periodogram, and the REDFIT techniques. The main results are as follows: (1) a QPO of ∼33 days (>2.5σ) is found during the flare of 117 days (MJD 55008–55125) from 2009 June to November. Interestingly, the same QPO (∼39 days) reappeared in the flaring duration from MJD 59430 to 59585, with the confidence level of >4σ. (2) Another transient QPO of ∼91 days with a significance of >3.8σ is found during a period with 455 days (MJD 58430–58985) from 2019 February to 2020 May. Under the premise of considering the QPOs reported in the literature, the QPO of ∼40 days is repeated three times and the QPO of ∼91 days is repeated twice. We discuss several physical models explaining the QPOs of this blazar. Our study may suggest that the two QPOs originate from the twin jets of the binary black holes at the center of 3C 279. The repeated occurrence of QPOs of a similar scale strongly supports the geometric scenario of a blob spiraling within the jet. Furthermore, the hypothesis of a sheath in the jet may also be a potential explanation for the repetitive γ-ray flare patterns observed in the light curve.

The Jet and Resolved Features of the Central Supermassive Black Hole of M87 Observed with EHT in 2017—Comparison with the GMVA 86 GHz Results

M87 is the best target for studying black hole accretion and jet formation. Reanalysis of the Event Horizon Telescope (EHT) public data at 230 GHz shows a core–knots structure at the center and jet features. We here compare this with the new results of GMVA at 86 GHz showing a spatially resolved central core. There are similarities and differences between the two. At 86 GHz, “two bright regions” are seen on the ring in the core. The “core–knot–west knot” triple structure in the 230 GHz image shows apparent appearance of two peaks similar to the “two bright regions” when convolved with the GMVA beam. This similarity suggests that both frequencies reveal the same objects in the core area. Protrusions are observed on both the south and north sides of the core at both frequencies, becoming prominent and winglike at 230 GHz. The 86 GHz image shows a triple ridge jet structure, while the 230 GHz image shows only a bright central ridge with two roots. Both frequencies show a shade between the core and the central ridge. To detect the faint features from the EHT2017 data, we found that the use of all baseline data is essential. Using all baseline data, including the ultrashort data, revealed the jet and faint structures. Without the ultrashort baselines, these structures were not detectable. The lack of detection of any faint structures other than the ring in the M87 data by the EHTC is presumably due to the exclusion of ultrashort baselines from their analysis.

Magnetic field properties inside the jet of Mrk 421

Aims.We aim to probe the magnetic field geometry and particle acceleration mechanism in the relativistic jets of supermassive black holes.Methods.We conducted a polarimetry campaign from radio to X-ray wavelengths of the high-synchrotron-peak (HSP) blazar Mrk 421, including Imaging X-ray Polarimetry Explorer (IXPE) measurements from 2022 December 6–8. During the IXPE observation, we also monitored Mrk 421 usingSwift-XRT and obtained a single observation withXMM-Newtonto improve the X-ray spectral analysis. The time-averaged X-ray polarization was determined consistently using the event-by-event Stokes parameter analysis, spectropolarimetric fit, and maximum likelihood methods. We examined the polarization variability over both time and energy, the former via analysis of IXPE data obtained over a time span of 7 months.Results.We detected X-ray polarization of Mrk 421 with a degree of ΠX = 14 ± 1% and an electric-vector position angleψX = 107 ± 3° in the 2–8 keV band. From the time variability analysis, we find a significant episodic variation inψX. During the 7 months from the first IXPE pointing of Mrk 421 in 2022 May,ψXvaried in the range 0° to 180°, while ΠXremained relatively constant within ∼10–15%. Furthermore, a swing inψXin 2022 June was accompanied by simultaneous spectral variations. The results of the multiwavelength polarimetry show that ΠXwas generally ∼2–3 times greater than Π at longer wavelengths, whileψfluctuated. Additionally, based on radio, infrared, and optical polarimetry, we find that the rotation ofψoccurred in the opposite direction with respect to the rotation ofψXand over longer timescales at similar epochs.Conclusions.The polarization behavior observed across multiple wavelengths is consistent with previous IXPE findings for HSP blazars. This result favors the energy-stratified shock model developed to explain variable emission in relativistic jets. We considered two versions of the model, one with linear and the other with radial stratification geometry, to explain the rotation ofψX. The accompanying spectral variation during theψXrotation can be explained by a fluctuation in the physical conditions, for example in the energy distribution of relativistic electrons. The opposite rotation direction ofψbetween the X-ray and longer wavelength polarization accentuates the conclusion that the X-ray emitting region is spatially separated from that at longer wavelengths. Moreover, we identify a highly polarized knot of radio emission moving down the parsec-scale jet during the episode ofψXrotation, although it is unclear whether there is any connection between the two events.

Three-dimensional GRMHD Simulations of Neutron Star Jets

Neutron stars and black holes in X-ray binaries are observed to host strong collimated jets in the hard spectral state. Numerical simulations can act as a valuable tool in understanding the mechanisms behind jet formation and its properties. Although there have been significant efforts in understanding black hole jets from general relativistic magnetohydrodynamic (GRMHD) simulations in recent years, neutron star jets still remain poorly explored. We present the results from three-dimensional GRMHD simulations of accreting neutron stars with oblique magnetospheres for the very first time. The jets in our simulations are produced due to the anchored magnetic field of the rotating star in analogy with the Blandford–Znajek process. We find that for accreting stars, the star–disk magnetic field interaction plays a significant role, and as a result, the jet power becomes directly proportional to , where Φjet is the open flux in the jet. The jet power decreases with increasing stellar magnetic inclination, and finally, for an orthogonal magnetosphere, it reduces by a factor of ≃2.95 compared to the aligned case. We also find that in the strong propeller regime, with a highly oblique magnetosphere, the disk-induced collimation of the open stellar flux preserves parts of the striped wind, resulting in a striped jet.

Blandford-Znajek jets in MOdified Gravity

General relativity (GR) will be imminently challenged by upcoming experiments in the strong gravity regime, including those testing the energy extraction mechanisms for black holes. Motivated by this, we explore magnetospheric models and black hole jet emissions in Modified Gravity (MOG) scenarios. Specifically, we construct new power emitting magnetospheres in a Kerr-MOG background which are found to depend non-trivially on the MOG deformation parameter. This may allow for high-precision tests of GR. In addition, a complete set of analytic solutions for vacuum magnetic field configurations around static MOG black holes are explicitly derived, and found to comprise exclusively Heun's polynomials.

Optical circular polarization of blazar S4 0954+65 during high linear polarized states

Optical circular polarization observations can directly test the particle composition in black hole jets. We report the first observations of the BL Lac-type object S4 0954+65 in high linear polarized states. While no circular polarization was detected, we were able to place upper limits of < 0.5% at 99.7% confidence. Using a simple model and our novel optical circular polarization observations, we can constrain the allowed parameter space for the magnetic field strength and composition of the emitting particles. Our results favor models that require magnetic field strengths of only a few Gauss and models in which the jet composition is dominated by electron-positron pairs. We discuss our findings in the context of typical magnetic field strength requirements for blazar emission models.

High-energy Neutrinos from the Inner Circumnuclear Region of NGC 1068

High-energy neutrinos are detected by the IceCube Observatory in the direction of NGC 1068, the archetypical type II Seyfert galaxy. The neutrino flux, surprisingly, is more than an order of magnitude higher than the γ-ray upper limits at measured TeV energy, posing tight constraints on the physical conditions of a neutrino production site. We report an analysis of the submillimeter, mid-infrared, and ultraviolet observations of the central 50 pc of NGC 1068 and suggest that the inner dusty torus and the region where the jet interacts with the surrounding interstellar medium (ISM) may be a potential neutrino production site. Based on radiation and magnetic field properties derived from observations, we calculate the electromagnetic cascade of the γ-rays accompanying the neutrinos. When injecting protons with a hard spectrum, our model may explain the observed neutrino flux above ∼10 TeV. It predicts a unique sub-TeV γ-ray component, which could be identified by a future observation. Jet–ISM interactions are commonly observed in the proximity of jets of both supermassive and stellar-mass black holes. Our results imply that such interaction regions could be γ-ray-obscured neutrino production sites, which are needed to explain the IceCube diffuse neutrino flux.

Precessing jet nozzle connecting to a spinning black hole in M87.

The nearby radio galaxy M87 offers a unique opportunity to explore the connections between the central supermassive black hole and relativistic jets. Previous studies of the inner region of M87 revealed a wide opening angle for the jet originating near the black hole1-4. The Event Horizon Telescope resolved the central radio source and found an asymmetric ring structure consistent with expectations from general relativity5. With a baseline of 17 years of observations, there was a shift in the jet's transverse position, possibly arising from an 8- to 10-year quasi-periodicity3. However, the origin of this sideways shift remains unclear. Here we report an analysis of radio observations over 22 years that suggests a period of about 11 years for the variation in the position angle of the jet. We infer that we are seeing a spinning black hole that induces the Lense-Thirring precession of a misaligned accretion disk. Similar jet precession may commonly occur in other active galactic nuclei but has been challenging to detect owing to the small magnitude and long period of the variation.

How Supermassive Black Holes Shape Galaxies

We often think of black holes as mysterious objects that swallow everything that passes by them—from dust to entire stars. In reality, black holes “eat” in a messy way, while spinning around super quickly. They do not swallow their whole “meal” (such as interstellar gas and dust) in a mouthful. Part of it rotates around the spinning black hole, like a merry-go-round, and is eventually thrown outward as very bright, energetic streams of particles called jets. Supermassive black holes the centers of galaxies can create jets that travel very far and affect the birth of baby stars! In this article, we will describe some of the things we know about black holes, including how astronomers found out about black hole jets using powerful telescopes. We will also mention some of the things we still do not understand about black holes and their jets.

Possible contribution of X-ray binary jets to the Galactic cosmic ray and neutrino flux

ABSTRACT For over a century, the identification of high-energy cosmic ray (CR) sources remains an open question. For Galactic CRs with energy up to 1015 eV, supernova remnants (SNRs) have traditionally been thought the main candidate source. However, recent TeV γ-ray observations have questioned the SNR paradigm. Propagating CRs are deflected by the Galactic magnetic field, hence, γ-rays and neutrinos produced via inelastic hadronic interactions are the only means for unveiling the CR sources. In this work, we study the γ-ray and neutrino emission produced by CRs accelerated inside Galactic jets of stellar-mass black holes in X-ray binaries (BHXBs). We calculate the intrinsic neutrino emission of two prototypical BHXBs , Cygnus X–1 and GX 339–4, for which we have high-quality, quasi-simultaneous multiwavelength spectra. Based on these prototypical sources, we discuss the likelihood of the 35 known Galactic BHXBs to be efficient CR accelerators. Moreover, we estimate the potential contribution to the CR spectrum of a viable population of BHXBs that reside in the Galactic plane. When these BHXBs go into outburst, they may accelerate particles up to hundreds of TeV that contribute to the diffuse γ-ray and neutrino spectra while propagating in the Galactic medium. Using HERMES, an open-source code that calculates the hadronic processes along the line of sight, we discuss the contribution of BHXBs to the diffuse γ-ray and neutrino fluxes, and compare these to their intrinsic γ-ray and neutrino emissions. Finally, we discuss the contribution of BHXBs to the observed spectrum of Galactic CRs.

Particle Injection and Nonthermal Particle Acceleration in Relativistic Magnetic Reconnection* *Released January 10, 2023.

Magnetic reconnection in the relativistic regime has been proposed as an important process for the efficient production of nonthermal particles and high-energy emission. Using fully kinetic particle-in-cell simulations, we investigate how the guide-field strength and domain size affect the characteristic spectral features and acceleration processes. We study two stages of acceleration: energization up until the injection energy γ inj and further acceleration that generates a power-law spectrum. Stronger guide fields increase the power-law index and γ inj, which suppresses acceleration efficiency. These quantities seemingly converge with increasing domain size, suggesting that our findings can be extended to large-scale systems. We find that three distinct mechanisms contribute to acceleration during injection: particle streaming along the parallel electric field, Fermi reflection, and the pickup process. The Fermi and pickup processes, related to the electric field perpendicular to the magnetic field, govern the injection for weak guide fields and larger domains. Meanwhile, parallel electric fields are important for injection in the strong guide-field regime. In the post-injection stage, we find that perpendicular electric fields dominate particle acceleration in the weak guide-field regime, whereas parallel electric fields control acceleration for strong guide fields. These findings will help explain the nonthermal acceleration and emission in high-energy astrophysics, including black hole jets and pulsar wind nebulae.

The Bardeen–Petterson effect as an observable support for the Blanford–Payne process black hole jet production

The Bardeen–Petterson effect as an observable support for the Blanford–Payne process black hole jet production

A ring-like accretion structure in M87 connecting its black hole and jet

The nearby radio galaxy M87 is a prime target for studying black hole accretion and jet formation1,2. Event Horizon Telescope observations of M87 in 2017, at a wavelength of 1.3 mm, revealed a ring-like structure, which was interpreted as gravitationally lensed emission around a central black hole3. Here we report images of M87 obtained in 2018, at a wavelength of 3.5 mm, showing that the compact radio core is spatially resolved. High-resolution imaging shows a ring-like structure of {8.4}_{-1.1}^{+0.5} Schwarzschild radii in diameter, approximately 50% larger than that seen at 1.3 mm. The outer edge at 3.5 mm is also larger than that at 1.3 mm. This larger and thicker ring indicates a substantial contribution from the accretion flow with absorption effects, in addition to the gravitationally lensed ring-like emission. The images show that the edge-brightened jet connects to the accretion flow of the black hole. Close to the black hole, the emission profile of the jet-launching region is wider than the expected profile of a black-hole-driven jet, suggesting the possible presence of a wind associated with the accretion flow.

HINOTORI and Its Perspectives in the Black-Hole Jet Study

Simultaneous multi-band very long baseline interferometry (VLBI) observations at millimeter wavelengths have huge potential for various science cases. However, there exist difficulties in expanding the scientific targets, as the sensitivity of radio telescopes at millimeter wavelengths is typically lower compared to that at centimeter wavelengths. In order to realize high-sensitivity mm-VLBI observations in the East Asia region, we are promoting the HINOTORI (Hybrid Installation project in NObeyama, Triple-band ORIented) project, which aims to launch the wide-band and simultaneous triple-band (22/43/86 GHz) VLBI system with the Nobeyama 45 m Radio Telescope (NRO45). The simultaneous 22/43 GHz observation mode has already been operated for the open-use program. We have recently completed the performance evaluation of the receiver and observing system at 86 GHz. In addition, a new wide-band VLBI back-end system has been installed on the NRO45 and the performance of this receiving system has been found to be sufficient to meet scientific requirements. Currently, we are performing commissioning observations to establish regular VLBI operation with simultaneous triple-band mode together with the Korean VLBI Network. The participation of the NRO45 is expected to strengthen the mm-VLBI observation network in the East Asia region and to be a very powerful addition with respect to the science of of black hole jets.

Enabling Transformational ngEHT Science via the Inclusion of 86 GHz Capabilities

We present a case for significantly enhancing the utility and efficiency of the ngEHT by incorporating an additional 86 GHz observing band. In contrast to 230 or 345 GHz, weather conditions at the ngEHT sites are reliably good enough for 86 GHz to enable year-round observations. Multi-frequency imaging that incorporates 86 GHz observations would sufficiently augment the (u,v) coverage at 230 and 345 GHz to permit detection of the M87 jet structure without requiring EHT stations to join the array. The general calibration and sensitivity of the ngEHT would also be enhanced by leveraging frequency phase transfer techniques, whereby simultaneous observations at 86 GHz and higher-frequency bands have the potential to increase the effective coherence times from a few seconds to tens of minutes. When observation at the higher frequencies is not possible, there are opportunities for standalone 86 GHz science, such as studies of black hole jets and spectral lines. Finally, the addition of 86 GHz capabilities to the ngEHT would enable it to integrate into a community of other VLBI facilities—such as the GMVA and ngVLA—that are expected to operate at 86 GHz but not at the higher ngEHT observing frequencies.