Next Article in Journal
Cosmological Model Tests with JWST
Next Article in Special Issue
Expectations for Horizon-Scale Supermassive Black Hole Population Studies with the ngEHT
Previous Article in Journal
Hidden Spectra Treasures in the Foster Archive: A Pilot Study of the Be Stars α Eri, α Col, ω Car and η Cen
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relativistic Signatures of Flux Eruption Events near Black Holes

1
Center for Astrophysics, Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA
2
Black Hole Initiative, Harvard University, 20 Garden Street, Cambridge, MA 02138, USA
3
Department of Physics, Princeton University, Princeton, NJ 08544, USA
4
Department of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ 08544, USA
5
Center for Computational Astrophysics, Flatiron Institute, 162 Fifth Avenue, New York, NY 10010, USA
6
Institute for Theory and Computation (ITC), Harvard University, 60 Garden Street, Cambridge, MA 02138, USA
*
Authors to whom correspondence should be addressed.
John Harvard Distinguished Science and ITC Fellow.
Submission received: 12 October 2022 / Revised: 8 November 2022 / Accepted: 10 November 2022 / Published: 24 November 2022

Abstract

:
Images of supermassive black holes produced using very long baseline interferometry provide a pathway to directly observing effects of a highly curved spacetime, such as a bright “photon ring” that arises from strongly lensed emission. In addition, the emission near supermassive black holes is highly variable, with bright high-energy flares regularly observed. We demonstrate that intrinsic variability can introduce prominent associated changes in the relative brightness of the photon ring. We analyze both semianalytic toy models and GRMHD simulations with magnetic flux eruption events, showing that they each exhibit a characteristic “loop” in the space of relative photon ring brightness versus total flux density. For black holes viewed at high inclination, the relative photon ring brightness can change by an order of magnitude, even with variations in total flux density that are comparatively mild. We show that gravitational lensing, Doppler boosting, and magnetic field structure all significantly affect this feature, and we discuss the prospects for observing it in observations of M87 and Sgr A with the next-generation Event Horizon Telescope.

1. Introduction

The Event Horizon Telescope released the first images of M87 in 2019 and the first images of Sgr A in 2022 [1,2], enabling new measurements of black hole accretion flow properties directly from VLBI data and event-horizon-scale images. A theorized component of all black hole images is the photon ring: a thin annulus of light composed of photons traveling on nearly-bound geodesics. This photon ring splits into a series of self-similar subrings, each of which reflects a different degree of light-bending around the hole see, e.g., [3,4,5,6,7,8,9].
In particular, each subring is labeled by an index n, which counts the number of half-orbits that a photon completes on its trajectory from emitter to observer. In the optically thin limit, each successive subring has similar brightness but is exponentially demagnified, with the demagnification related to Lyapunov exponents that are governed by the properties of unstable spherical orbits of null geodesics see, e.g., [3,9,10,11].
One consequence of the demagnification is that the direct image ( n = 0 ) of a black hole tends to be the dominant source of observed flux, with indirect images ( n 1 ) appearing exponentially suppressed. Hence, an important quantity is the Photon ring Flux Ratio (PFR), defined as the fraction of total flux contained in a particular subimage. We focus on the n = 1 ratio, which we denote as f 1 :
f 1 F 1 F tot .
This quantity depends both on the spacetime (which entirely determines the relative demagnification of the subring), the emission geometry, and the magnetic field structure. Moreover, it can potentially be measured with the next-generation Event Horizon Telescope (ngEHT) by using modeling methods that isolate the contribution of the photon ring e.g., [12,13].
General relativistic magneto-hydrodynamic (GRMHD) simulations and general relativistic radiative transfer (GRRT) are important tools for connecting observations of the photon ring to the underlying plasma and emission physics of the accretion disk e.g., [14,15]. GRMHD simulations used to interpret the 2017 EHT observations favor strongly magnetized gas accreting onto Sgr A and M87  [16,17]. In the magnetically arrested disk MAD [18,19] limit, the magnetic field near the BH becomes strong enough to vertically squeeze the accreting gas. These fields ultimately undergo magnetic reconnection, allowing a bundle of vertical fields to escape from the vicinity of the BH, and in the process, eject out a large portion of the disk e.g., [20]. Observable signatures of such MAD system behavior in Sgr A* were recently reported by Wielgus et al. [21]. These “magnetic flux eruption events” exhibit large gas temperatures, strong vertical fields and occur quasi-periodically in the MAD state e.g., [22], and hence, are a prime candidate for the origin of high-energy flux eruptions, e.g., in Sgr A and M87 . An example of such a flux eruption is depicted in Figure 1, and a more detailed description of the accretion flow’s response to these events is outlined in Appendix A.
During these eruptions, the flux density in the image can drop by a factor of ∼10 (see Jia et al. in prep for more details about the lightcurve). Observations of Sagittarius A* indicate a flux density drop of a factor of 2 at millimeter wavelengths following a high energy flaring event [23,24]. Furthermore, simulations indicate that individual eruptions may last for ∼100–300 M,1 with the flux bundle eventually dissipating in the disk. This timescale corresponds to a few weeks to months for M87 or an ∼hour for Sgr A , and is reasonably consistent with the analysis of Sgr A* X-ray flares population [25]. Such a long evolution period allows us to possibly capture eruption events with instruments such as the ngEHT and to resolve the flaring state of a SMBH. The presence of strong vertical fields and the ejection of gas present an exceptional opportunity to probe the detailed structure of the flux eruption through its observable features in polarization and variability e.g., [26,27].
In this work, we explore how the dynamics of flux eruptions are manifest in image morphology, with specific attention to the underlying factors that directly control the relative brightness of the photon ring during these events. In Section 2, we introduce a toy model to simplify the depiction of a flux eruption event, allowing us to separate the individual influences of Doppler boosting and gravitational lensing on the photon ring’s appearance. In Section 3, we repeat a similar analysis on a set of GRMHD data, showing that the results of the toy model are recovered provided that the emission profile is not changed. In Section 4, we show that the introduction of a magnetic field (with a thermal emission profile) dramatically alters the photon ring brightness. Finally, in Section 5, we discuss our results in the context of observable targets for future EHT and ngEHT science.

2. Toy Model

In this section, we develop a toy model to illuminate the effects of Doppler boosting on the PFR during flux eruption events.

2.1. Description of the Model

During a flux eruption event, a substantial portion of the high-density material in the disk is ejected, leaving a low-density magnetosphere (with an equatorial current sheet) in the inner region. Our toy model therefore consists of a half-wedge (i.e., a wedge that spans 180 in ϕ ) surrounding a Kerr black hole, with emission extending down to the horizon. In practice, one could experiment with smaller wedges that span a narrower range of ϕ values, but we restrict our focus to the half-wedge for simplicity. The emission from the wedge represents the near-horizon sub-millimeter emission seen during a flux eruption event.
From here, we build off of the semianalytic models of Gold et al. [15], modifying the prescriptions so that emission is confined to exactly one half of the spacetime. Among these semianalytic test models, we adapt “Test 5”, which entails a simple, isotropic emissivity function (i.e., independent of magnetic field direction) and a thin scale height for the disk:
n ( r ) exp 1 2 r 10 2 + 100 3 cos θ 2 , j ν ( r ) n ( r ) , α ν ( r ) = 10 6 ν 230 GHz 2.5 × j ν ( r ) .
Here, n is the number density, j ν is the emissivity, and α ν is the absorptivity, with θ and r taking on their Boyer-Lindquist coordinate values. In addition to cutting out half of the emitting region, we also modify the original prescription of Gold et al. [15] so that distance of the black hole matches that of M87 , the spin of the black hole is a = + 0.94 , and the overall density is rescaled so F 230 0.5 Jy for one of the snapshots.

2.2. Photon Ring Flux Ratios

During a flux eruption event, there is a disruption of the approximate axisymmetry in the accretion structure. One expects the reconnection layer powering the eruption to emit strongly in the X-ray, leaving the cool gas in the disk to emit in the sub-millimeter and radio. Hence, we expect the PFR to depend strongly on the position of the observer relative to the flux eruption itself. This orientation is encoded through both the observer’s polar inclination angle ( θ cam ), which is measured from the black hole’s spin axis, as well as the observer’s azimuthal angle ( ϕ cam ). We orient the azimuthal coordinate so that ϕ cam = 0 when the region of highest density (and hence highest emissivity) is positioned directly in front of the hole. In a physical black hole, such a scenario would take place when the flux tube is positioned directly behind the black hole.
We ray-trace the toy model for a range of values of ϕ cam as a proxy for tracking the flux eruption over time, as increasing ϕ cam is equivalent to rotating the accretion disk with respect to the observer. We ray-trace the disk for both a “static” fluid rotation profile ( u ϕ = 0 ) as well as an asymptotically Keplerian rotation profile ( u ϕ ( r sin θ ) 1 / 2 as r ), corresponding to 0 = 0 and 0 = 1 in Gold et al. [15], respectively. This allows us to isolate the effects of Doppler boosting, which is relevant only when objects are in motion. We further repeat this procedure for θ cam = 17 and θ cam = 80 , representing both the low-inclination and high-inclination limits of the viewing geometry.
The ray-tracing is performed using the adaptive sampling scheme of ipole [28,29]. The field of view is taken to be 160 μ as with an effective pixel size of 0.15 μ as, sufficient to fully resolve the n = 1 subring. Each image is decomposed into its individual subrings using the procedure described by Gelles et al. [29]. In this scheme, the intensity of a pixel in the n th subimage is computed by performing radiative transfer along the corresponding geodesic’s n th pass around the black hole.
In Figure 2, we plot the ratio f 1 for each of these configurations and rotation profiles as a function of the total image flux. We normalize the total image flux F tot with a mean flux F ¯ tot taken over all values of ϕ cam , which effectively represents a time-averaged flux. Since we may be able to measure the ratio f 1 along with the total compact flux using the ngEHT, we choose to plot f 1 against F tot / F ¯ tot for each camera azimuthal position ϕ cam . The flux eruption induces a strong correlation between the total image flux and the fractional flux contained in the subring. As the half-disk orbits around the black hole, the PFR changes dramatically, more than doubling over the course of one revolution in the edge-on limit.
The specific shape of the curves traced out in Figure 2 is due to a confluence of numerous factors. In the low inclination limit ( θ cam = 17 ), the PFR remains relatively constant over azimuth, regardless of whether or not the disk is rotating. This is because the fluid velocity is confined to the midplane, so with a low-inclination viewing geometry, the half-disk cannot acquire a large velocity tangent to null geodesics that reach the observer.
The specific value of the PFR is consistent with theoretical predictions as well. For an a = + 0.94 Kerr black hole viewed at θ cam = 17 , the Lyapunov exponent ranges from γ = 2.36 to γ = 2.76 , depending on the azimuthal screen coordinate. Following Johnson et al. [9], in the asymptotic limit of large n (as well as optical transparency and axisymmetry), one expects
F n + 1 F n e γ F 1 F tot = F 1 F 0 1 1 + F 1 F 0 + . . . e γ 1 n = 0 e γ n = e γ e 2 γ ,
which ranges from 0.059 to 0.085 for this particular black hole. This range falls not far from the PFR of the low inclination snapshots in Figure 2.
On the other hand, in the high-inclination case ( θ cam = 80 ), the brightness of the photon ring depends heavily on the location of the observer with respect to the flux eruption. In the case of a static rotation profile (Figure 2; left panel), the PFR is primarily determined by the relative magnification from gravitational lensing. As the half-disk passes behind the black hole, more of the received flux is bent around the hole, leading to a magnification of the n = 1 image and a resultant brightening of the photon ring. Indeed, the PFR is largest when ϕ cam 180 and the emission wedge is located behind the black hole. The phases of increasing magnification and demagnification are demarcated with arrows in Figure 2.
However, when the disk is rotating (Figure 2; right panel), Doppler effects become important. The maximum and minimum PFR’s still occur when ϕ cam 180 and ϕ cam 0 respectively, as the gravitational lensing is identical to the non-rotating case. However, Doppler boosting stretches the PFR curve out horizontally, as the recessional speed of the eruption is largest when ϕ cam 90 and ϕ cam 270 . Indeed, in the bottom left snapshot of Figure 2, one can see the Doppler boosted direct image of the flux eruption on the left, leading to a Doppler deboosted indirect image on the right. In the bottom right snapshot, one sees the Doppler deboosted direct image of the flux eruption on the right, leading to a Doppler boosted indirect image on the left. The Doppler boosting and deboosting phases (which refer to the indirect image) are also demarcated with arrows in Figure 2. 2
We next return to the flaring GRMHD simulations to evaluate whether the gross behavior of the toy model is seen under more physically plausible circumstances.

3. GRMHD

In this section, we describe the GRMHD simulation and ray-tracing techniques that we use to measure subring fluxes in various accreting environments.

3.1. Procedure

To further investigate the connection between photon ring brightness and flux eruption events, we ray-trace the ideal GRMHD simulation from Ripperda et al. [20], performed using the h-amr code [30]. The simulation shows MAD accretion cycles, separated by prominent plasmoid-mediated magnetic reconnection events through which magnetic flux is expelled from the event horizon. The dimensionless black hole spin parameter is a = 15 / 16 and the effective grid resolution is N r × N θ × N ϕ = 5376 × 2304 × 2304 defined for logarithmic Kerr-Schild spherical polar coordinates. The simulation was evolved to t = 10,000 M.
We ray-trace the simulation using the adaptive sampling scheme of ipole. We used the mass and distance of M87 , as for the toy model discussed in Section 2, but we reduce the FOV to 80 μ as. The GRMHD scale factor is calibrated so that the average flux density is F 230 0.5  Jy. In generating the images, we rotated our azimuthal coordinates clockwise by 150 from the ipole default to align the region of highest synchrotron emissivity with the observer at ϕ cam = 0 , following the conventions of the toy model described in Section 2.
Over the course of the simulation, the flow exhibits a quiescent state (wherein material accretes steadily onto the black hole) and a flaring state (wherein magnetic flux is expelled outward as described in Section 1). We isolate two GRMHD time slices of the simulation that directly showcase these different states: t = 8858 M is identified with quiescent accretion, and t = 9553 M is identified with the flux eruption event. During the latter state, roughly half the fluid cells have radially inward velocities (hence accreting) while half have radially outward velocities (hence ejecting). The structure of the flux eruption in the GRMHD simulation is thus broadly consistent with our choice of toy model in Section 2.
We note that not only are there many smaller-scale eruption events that occur during the GRMHD simulation, but the size of the ejection region for any one specific eruption event also changes over time. We used our toy model to represent the peak of a particularly large flux eruption event from this simulation so as to capture the near-horizon image structure during a potentially bright high-energy flare.

3.2. Photon Ring Flux Ratios

As with the toy model, to a leading order approximation, we can ray-trace a single time-slice of data for a range of values of ϕ cam as a proxy for tracking the eruption over time. In particular, this eliminates the need to account for the time-dependence of the eruption shape, which would introduce additional non-linearities into our analysis.
Unlike the toy model, however, the GRMHD data contains information about magnetic fields. We expect the orientation of these fields in the accretion flow to directly control the brightness of the photon ring; for synchrotron processes, the plasma emissivity depends on k × B , where k is the spatial 3-vector of the null geodesic and B is the magnetic field, both of which are measured in the local Minkowski frame of the source [31]. Indeed, for M87 and Sgr A at millimeter wavelengths, the specific intensity will be roughly proportional to sin 2 ζ , where ζ is the “pitch angle” between the emitted wavevector and the magnetic field see, e.g., [31].
To compare the GRMHD case to Figure 2, we first ray-trace the data for a range of ϕ cam values with sin ζ 1 everywhere, thus eliminating magnetic field directional dependence and causing the GRMHD emissivity prescription to resemble that of the toy model. We then construct curves of the GRMHD PFR as a function of ϕ cam for both θ cam = 17 and θ cam = 80 , and the results are plotted in Figure 3.
Qualitatively, the high-inclination curves of Figure 3 match those of semianalytic case in Figure 2. All four phases (magnification, boosting, demagnification, and deboosting) are evident. Also evident are the Doppler boosted direct image of the eruption in the lower left panel of Figure 3, as well as the Doppler boosted indirect image of the eruption in the lower right panel Figure 3.
However, the low-inclination curve now rotates in the opposite direction (i.e., clockwise corresponds to increasing  ϕ cam ) and has been stretched out dramatically, indicating a break in axisymmetry; the region emitting in the sub-millimeter has been restricted to a narrower range of ϕ values. The effects of magnetic field magnitude are primarily responsible for this change.
Next, we will investigate the role of the magnetic field direction by re-introducing the sin ζ dependence in the synchrotron emissivity.

4. Magnetic Fields

In this section, we describe the specific role that magnetic field direction plays in altering the relative intensity of the direct vs. indirect images.

4.1. Background

The pitch angle ζ , which encodes the directional dependence of synchrotron emissivity, can dramatically alter the brightness of the image in a way that cannot be predicted from properties of the spacetime alone. In particular, ζ is different for the direct and indirect images and therefore directly influences the PFR. For many field configurations, the pitch angle will be larger for the strongly lensed geodesics and can hence artificially inflate the relative brightness of the photon ring.
For M87 , it is believed that the accretion flow has a strongly poloidal (i.e., a mix of vertical and radial) magnetic field and that the hole is viewed from Earth at a relatively low inclination [1,32]. For instance, for a stationary emitter at the ISCO ( r = 6 M) of a Schwarzschild black hole viewed at face-on inclination in a purely vertical magnetic field, the pitch angle is ζ = 19 . 4 for the n = 0 image and ζ = 41 . 8 for the n = 1 image. In this case, sin 2 ζ increases by more than a factor of 4 in the indirect image, leading to a corresponding increase in the relative brightness of the photon ring.
Such differences in field configurations are particularly relevant to our discussion of flux eruption events. During the quiescent state, the magnetic field is less ordered and accretes together with the infalling gas. However, during the flaring state, the magnetic field in the remaining magnetosphere (after material has been ejected) is connected to the event horizon and the jet, rather than the disk. This causes the magnetic field to be nearly equatorial and contain a magnetic null (i.e., a current sheet where reconnection takes place) separating the northern and southern jets. These contrasting magnetic field configurations may result in observable differences in direct emission that may be more prominent in polarized emission.

4.2. Effects of Magnetic Field Direction on GRMHD PFR’s

To demonstrate the effects of magnetic field direction on the observed PFR’s, we ray-trace the GRMHD simulation anew with a synchrotron emissivity profile that depends appropriately on sin ζ (as would normally be done). The resultant PFR’s are shown in Figure 4. Qualitatively, Figure 4 matches Figure 3. All four phases are once again identifiable in the high-inclination PFR curve, although the Doppler deboosting phase has shrunk to a smaller region of the phase space.
The low-inclination curve, while similar in orientation to Figure 3, has been dilated and translated upward, indicating the presence of a brighter photon ring. This transformation is consistent with a predominantly vertical magnetic field configuration in the accretion flow. As explained in Section 4.1, the contrast in pitch angle between direct/indirect images is particularly strong when the magnetic field threading the disk points directly at (or directly away from) the observer, as is the case for a vertical field viewed nearly face-on. In this case, we expect the PFR to increase, and we subsequently expect the behavior seen in Figure 4.
The correspondence between Figure 2, Figure 3 and Figure 4 suggest that the half-disk is an appropriate (albeit simplistic) toy model for broadly representing the shape of a flux eruption event, as well as its signatures of gravitational lensing and Doppler effects. However, a dynamically important magnetic field is necessary to fully capture the correct emissivity profile and resultant PFR.

5. Discussion

We have analyzed sub-millimeter images of black holes from semi-analytic toy models and GRMHD simulations with plasmoid-mediated magnetic reconnection events. We have shown that these flux eruption events illuminate clear signatures of light bending near a black hole, and that these specific signatures depend on a multitude of factors stemming from both the spacetime and astrophysics at hand. In particular, at sub-millimeter wavelengths, the relative brightness of the n = 1 photon ring, f 1 , is sensitive to the orbital phase of the flux eruption and can exceed 40% in systems viewed at high inclination. Furthermore, we have shown that during these flux eruption events, the PFR traces out a generic looping pattern over time, which is due to a combination of gravitational lensing, Doppler boosting, and changing synchrotron emissivity. Each of these factors is significant.
These striking features in the time-variable images of black holes are an exceptional opportunity for studies with the next generation EHT (ngEHT). In particular, the ngEHT will enable time-resolved images of both Sgr A and M87 over hundreds-to-thousands of gravitational timescales, sufficient to catch rare events and to monitor their evolution. The enhanced angular resolution, baseline coverage, and sensitivity of the ngEHT will enable detailed studies of the photon ring see, e.g., [13], allowing measurements of the PFR during flux eruptions and flares 3. Finally, simultaneous multi-frequency capabilities of the ngEHT will allow for spectral index and polarization measurements during flares see, e.g., [33], providing additional information that can resolve the degeneracies between effects from the curved spacetime and emitting plasma. Along with the ngEHT, simultaneous multiwavelength coverage at higher energies such as near-infrared and X-rays for Sgr A e.g., [26,34] or TeV for M87 e.g., [35] could constrain both the flow structure and particle acceleration mechanisms during eruption events.
We emphasize that our model does not capture several intricacies of the problem at hand. While our moving camera does serve as an effective time-proxy for the rotation of the accretion flow, we employ a “fast-light” ray tracing algorithm that does not take into account the time delay between direct and indirect images. While “slow-light” simulations are significantly more computationally expensive, time delays can be integral to our study of hotspot lightcurves and images e.g., [36,37], so they should be employed in future studies of the PFR. Qualitatively, we suspect that the effects of time delay on the PFR curves are two-fold. First, the curve should rotate, as the point of maximal/minimal lensing will be shifted. Second, the curve should shrink, as direct/indirect images will no longer be diametrically opposed on the observer’s screen, leading to a smaller range of possible values for f 1 . In any case, we do not expect the time delay to modify our conclusions, as the generic shape of the PFR curve will remain consistent.
This study has provided a first glimpse of the scientific opportunities that may be possible with time-resolved studies of the photon ring during flux eruptions of magnetically arrested disks. Additional crucial topics for future studies include the accessible signatures in polarization, the effects of optical depth, and reconstructed movies with the ngEHT.

Author Contributions

Conceptualization, Z.G., K.C. and M.J.; methodology, Z.G. and K.C.; software, Z.G., K.C. and M.L.; formal analysis, Z.G. and K.C.; writing—original draft preparation, Z.G., K.C. and M.J.; writing—review and editing, Z.G., K.C., M.J., B.R. and M.L.; visualization, Z.G., K.C. and M.J.; supervision, K.C. and M.J. All authors have read and agreed to the published version of the manuscript.

Funding

We thank the National Science Foundation (AST-1716536 and AST-1935980) and the Gordon and Betty Moore Foundation (GBMF-10423) for financial support of this work. This work was supported in part by the Black Hole Initiative, which is funded by grants from the John Templeton Foundation and the Gordon and Betty Moore Foundation to Harvard University. K.C. is also supported in part by the Black Hole PIRE program (NSF grant OISE-1743747). B.R. is supported by a Joint Princeton/Flatiron Postdoctoral Fellowship. Research at the Flatiron Institute is supported by the Simons Foundation. M.L. is supported by John Harvard Distinguished Science Fellowship and ITC Fellowship.

Data Availability Statement

Simulation and raytraced data available on request to authors. Software: eht-imaging library [38], h-amr [30], ipole [28], Numpy [39], Matplotlib [40].

Acknowledgments

We thank Maciek Wielgus and Dominic Chang for useful discussions. This research was enabled by support provided by grant no. NSF PHY-1125915 along with a INCITE program award PHY129, using resources from the Oak Ridge Leadership Computing Facility, Summit, which is a US Department of Energy office of Science User Facility supported under contract DE-AC05- 00OR22725, as well as Calcul Quebec (http://www.calculquebec.ca, accessed on 9 November 2022) and Compute Canada (http://www.computecanada.ca, accessed on 9 November 2022). The computational resources and services used in this work were partially provided by facilities supported by the Scientific Computing Core at the Flatiron Institute, a division of the Simons Foundation. This research is part of the Frontera (Stanzione et al. [41]) computing project at the Texas Advanced Computing Center (LRAC-AST20008). Frontera is made possible by National Science Foundation award OAC-1818253.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BHBlack Hole
EHTEvent Horizon Telescope
FOVField Of View
GRMHDGeneral Relativistic Magneto Hydro Dynamic
GRRTGeneral Relativistic Ray Tracing
ISCOInner Stable Circular Orbit
MADMagnetically Arrested Disk
ngEHTnext generation Event Horizon Telescope
PFRPhoton ring Flux Ratio
SMBHSuper Massive Black Hole
VLBIVery Long Baseline Interferometry

Appendix A. Magnetic Flux Eruptions in GRMHD

From Figure 1, we see that the gas inspirals towards the black hole during the quiescent state ( t = 8858 M), while roughly half of the disk is ejected during the flux eruption event ( t = 9553 M). Here we briefly discuss the changes in other properties of the accretion flow brought about by magnetic flux eruptions. Figure A1 shows the midplane cross-section of the gas density ρ , plasma- β (i.e., the ratio of the thermal and magnetic pressures, p gas / p mag ), and gas temperature T gas = p gas / ρ . We also show the proxy for the 230 GHz thermal synchrotron emissivity j syn given by the EHT code comparison project Porth et al. [14]:
j syn = ρ 3 p gas 2 exp C ρ 2 B p gas 2 1 / 3 .
This emissivity prescription is designed to resemble the true synchrotron emissivity fitting function given in Leung et al. [42]. Following Porth et al. [14], we assume C = 0.2 such that the emission drops exponentially beyond a few gravitational radii. We further normalize j syn such that the volume-integrated total synchrotron emissivity is 1.
The quiescent state is characterized by gas-rich spiral features that interact with strong magnetic fields (indicated by β < 1 ) near the black hole. There are occasional reconnecting sites where magnetic dissipation leads to relativistic gas temperatures (i.e., T gas > 1 ) . overall, we see a roughly axisymmetric disk structure that produces an azimuthally-uniform emissivity profile.
On the other hand, the flux eruption state exhibits highly relativistic temperatures in the evacuated region. This occurs due to the formation of a long thin current sheet that ultimately destabilizes and undergoes reconnection, allowing the escape of a magnetic flux-tube [20]. The large gas temperatures could potentially produce high-energy emission that may be able to explain X-ray/ γ ray flares seen in low-luminosity super-massive black holes such as Sgr A and M87 . For the 230 GHz image, the synchrotron emissivity map predicts that the bulk of the emission originates in the highly-dense relatively low-temperature accreting region. The evacuated region, despite the high gas temperatures, produces little to no flux, leading to the unique azimuthally-dependent image morphologies shown in Figure 4.
Figure A1. Magnetic flux eruptions produce regions of relativistically hot gas that could potentially produce high-energy flares in Sgr A and M87 . This figure shows midplane cross-sections of the quiescent (top) and flux eruption (bottom) states from the GRMHD simulation. From left-to-right, panels show the gas density ρ , plasma- β , gas temperature T gas (in relativistic units), and a proxy for the 230 GHz synchrotron emissivity j syn .
Figure A1. Magnetic flux eruptions produce regions of relativistically hot gas that could potentially produce high-energy flares in Sgr A and M87 . This figure shows midplane cross-sections of the quiescent (top) and flux eruption (bottom) states from the GRMHD simulation. From left-to-right, panels show the gas density ρ , plasma- β , gas temperature T gas (in relativistic units), and a proxy for the 230 GHz synchrotron emissivity j syn .
Galaxies 10 00107 g0a1

Notes

1
We use natural units, i.e., G = c = 1 , and the length and time units are both defined only in terms of the black hole mass M.
2
The effects of black hole spin on these conclusions are minimal, and these results are very similar for a Schwarzschild ( a = 0 ) black hole. In the Schwarzschild case, Doppler effects vanish completely, as the zero-angular momentum frame is motionless everywhere.
3
By flares, we refer to transient bright emission usually observed at wavelengths much smaller than sub-millimeter, e.g., X-ray flares in Sgr A [25].

References

  1. Akiyama, K. et al. [Event Horizon Telescope Collaboration] The Shadow of the Supermassive Black Hole. Astrophys. J. Lett. 2019, 875, L1. [Google Scholar] [CrossRef]
  2. Akiyama, K. et al. [Event Horizon Telescope Collaboration] First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way. Astrophys. J. Lett. 2022, 930, L12. [Google Scholar] [CrossRef]
  3. Luminet, J.P. Image of a spherical black hole with thin accretion disk. Astron. Astrophys. 1979, 75, 228–235. [Google Scholar]
  4. de Vries, A. The apparent shape of a rotating charged black hole, closed photon orbits and the bifurcation set A4. Class. Quantum Gravity 2000, 17, 123–144. [Google Scholar] [CrossRef]
  5. Takahashi, R. Shapes and Positions of Black Hole Shadows in Accretion Disks and Spin Parameters of Black Holes. Astrophys. J. 2004, 611, 996–1004. [Google Scholar] [CrossRef]
  6. Beckwith, K.; Done, C. Extreme gravitational lensing near rotating black holes. Mon. Not. R. Astron. Soc. 2005, 359, 1217–1228. [Google Scholar] [CrossRef] [Green Version]
  7. Johannsen, T.; Psaltis, D. Testing the No-hair Theorem with Observations in the Electromagnetic Spectrum. II. Black Hole Images. Astrophys. J. 2010, 718, 446–454. [Google Scholar] [CrossRef] [Green Version]
  8. Gralla, S.E.; Holz, D.E.; Wald, R.M. Black hole shadows, photon rings, and lensing rings. Phys. Rev. D 2019, 100, 024018. [Google Scholar] [CrossRef] [Green Version]
  9. Johnson, M.D.; Lupsasca, A.; Strominger, A.; Wong, G.N.; Hadar, S.; Kapec, D.; Narayan, R.; Chael, A.; Gammie, C.F.; Galison, P.; et al. Universal interferometric signatures of a black hole’s photon ring. Sci. Adv. 2020, 6, eaaz1310. [Google Scholar] [CrossRef] [Green Version]
  10. Darwin, C. The Gravity Field of a Particle. Proc. R. Soc. Lond. Ser. A 1959, 249, 180–194. [Google Scholar] [CrossRef]
  11. Ohanian, H.C. The black hole as a gravitational “lens”. Am. J. Phys. 1987, 55, 428–432. [Google Scholar] [CrossRef]
  12. Broderick, A.E.; Pesce, D.W.; Tiede, P.; Pu, H.Y.; Gold, R. Hybrid Very Long Baseline Interferometry Imaging and Modeling with THEMIS. Astrophys. J. 2020, 898, 9. [Google Scholar] [CrossRef]
  13. Tiede, P.; Broderick, A.E.; Palumbo, D.C.M.; Chael, A. Measuring the Ellipticity of M 87* Images. arXiv 2022, arXiv:2210.13499. [Google Scholar]
  14. Porth, O.; Chatterjee, K.; Narayan, R.; Gammie, C.F.; Mizuno, Y.; Anninos, P.; Baker, J.G.; Bugli, M.; Chan, C.k.; Davelaar, J.; et al. The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project. Astrophys. J. Suppl. Ser. 2019, 243, 26. [Google Scholar] [CrossRef] [Green Version]
  15. Gold, R.; Broderick, A.E.; Younsi, Z.; Fromm, C.M.; Gammie, C.F.; Mościbrodzka, M.; Pu, H.Y.; Bronzwaer, T.; Davelaar, J.; Dexter, J.; et al. Verification of Radiative Transfer Schemes for the EHT. Astrophys. J. 2020, 897, 148. [Google Scholar] [CrossRef]
  16. Akiyama, K. et al. [Event Horizon Telescope Collaboration] First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring. Astrophys. J. Lett. 2019, 875, L5. [Google Scholar] [CrossRef]
  17. Akiyama, K. et al. [Event Horizon Telescope Collaboration] First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole. Astrophys. J. Lett. 2022, 930, L16. [Google Scholar] [CrossRef]
  18. Narayan, R.; Igumenshchev, I.V.; Abramowicz, M.A. Magnetically Arrested Disk: An Energetically Efficient Accretion Flow. Publ. Astron. Soc. Jpn. 2003, 55, L69–L72. [Google Scholar] [CrossRef] [Green Version]
  19. Igumenshchev, I.V.; Narayan, R.; Abramowicz, M.A. Three-dimensional Magnetohydrodynamic Simulations of Radiatively Inefficient Accretion Flows. Astrophys. J. 2003, 592, 1042–1059. [Google Scholar] [CrossRef]
  20. Ripperda, B.; Liska, M.; Chatterjee, K.; Musoke, G.; Philippov, A.A.; Markoff, S.B.; Tchekhovskoy, A.; Younsi, Z. Black Hole Flares: Ejection of Accreted Magnetic Flux through 3D Plasmoid-mediated Reconnection. Astrophys. J. Lett. 2022, 924, L32. [Google Scholar] [CrossRef]
  21. Wielgus, M.; Moscibrodzka, M.; Vos, J.; Gelles, Z.; Martí-Vidal, I.; Farah, J.; Marchili, N.; Goddi, C.; Messias, H. Orbital motion near Sagittarius A*. Constraints from polarimetric ALMA observations. Astron. Astrophys. 2022, 665, L6. [Google Scholar] [CrossRef]
  22. Narayan, R.; Chael, A.; Chatterjee, K.; Ricarte, A.; Curd, B. Jets in magnetically arrested hot accretion flows: Geometry, power, and black hole spin-down. Mon. Not. R. Astron. Soc. 2022, 511, 3795–3813. [Google Scholar] [CrossRef]
  23. Akiyama, K. et al. [Event Horizon Telescope Collaboration] First Sagittarius A* Event Horizon Telescope Results. II. EHT and Multiwavelength Observations, Data Processing, and Calibration. Astrophys. J. Lett. 2022, 930, L13. [Google Scholar] [CrossRef]
  24. Wielgus, M.; Marchili, N.; Martí-Vidal, I.; Keating, G.K.; Ramakrishnan, V.; Tiede, P.; Fomalont, E.; Issaoun, S.; Neilsen, J.; Nowak, M.A.; et al. Millimeter Light Curves of Sagittarius A* Observed during the 2017 Event Horizon Telescope Campaign. Astrophys. J. Lett. 2022, 930, L19. [Google Scholar] [CrossRef]
  25. Haggard, D.; Nynka, M.; Mon, B.; de la Cruz Hernandez, N.; Nowak, M.; Heinke, C.; Neilsen, J.; Dexter, J.; Fragile, P.C.; Baganoff, F.; et al. Chandra Spectral and Timing Analysis of Sgr A*’s Brightest X-ray Flares. Astrophys. J. 2019, 886, 96. [Google Scholar] [CrossRef] [Green Version]
  26. Abuter, R. et al. [GRAVITY Collaboration]. Constraining particle acceleration in Sgr A with simultaneous GRAVITY, Spitzer, NuSTAR, and Chandra observations. Astron. Astrophys. 2021, 654, A22. [Google Scholar] [CrossRef]
  27. Gelles, Z.; Himwich, E.; Johnson, M.D.; Palumbo, D.C.M. Polarized image of equatorial emission in the Kerr geometry. Phys. Rev. D 2021, 104, 044060. [Google Scholar] [CrossRef]
  28. Mościbrodzka, M.; Gammie, C.F. IPOLE—Semi-analytic scheme for relativistic polarized radiative transport. Mon. Not. R. Astron. Soc. 2018, 475, 43–54. [Google Scholar] [CrossRef] [Green Version]
  29. Gelles, Z.; Prather, B.S.; Palumbo, D.C.M.; Johnson, M.D.; Wong, G.N.; Georgiev, B. The Role of Adaptive Ray Tracing in Analyzing Black Hole Structure. Astrophys. J. 2021, 912, 39. [Google Scholar] [CrossRef]
  30. Liska, M.; Chatterjee, K.; Tchekhovskoy, A.; Yoon, D.; van Eijnatten, D.; Hesp, C.; Markoff, S.; Ingram, A.; van der Klis, M. H-AMR: A New GPU-accelerated GRMHD Code for Exascale Computing With 3D Adaptive Mesh Refinement and Local Adaptive Time-stepping. arXiv 2019, arXiv:1912.10192. [Google Scholar] [CrossRef]
  31. Narayan, R.; Palumbo, D.C.M.; Johnson, M.D.; Gelles, Z.; Himwich, E.; Chang, D.O.; Ricarte, A.; Dexter, J.; Gammie, C.F.; Chael, A.A.; et al. The Polarized Image of a Synchrotron-emitting Ring of Gas Orbiting a Black Hole. Astrophys. J. 2021, 912, 35. [Google Scholar] [CrossRef]
  32. Akiyama, K. et al. [Event Horizon Telescope Collaboration] First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon. Astrophys. J. Lett. 2021, 910, L13. [Google Scholar] [CrossRef]
  33. Ricarte, A.; Gammie, C.; Narayan, R.; Prather, B.S. Probing Plasma Physics with Spectral Index Maps of Accreting Black Holes on Event Horizon Scales. arXiv 2022, arXiv:2202.02408. [Google Scholar]
  34. Boyce, H.; Haggard, D.; Witzel, G.; Fellenberg, S.v.; Willner, S.P.; Becklin, E.E.; Do, T.; Eckart, A.; Fazio, G.G.; Gurwell, M.A.; et al. Multiwavelength Variability of Sagittarius A* in 2019 July. Astrophys. J. 2022, 931, 7. [Google Scholar] [CrossRef]
  35. Aharonian, F.; Akhperjanian, A.G.; Bazer-Bachi, A.R.; Beilicke, M.; Benbow, W.; Berge, D.; Bernlöhr, K.; Boisson, C.; Bolz, O.; Borrel, V.; et al. Fast Variability of Tera-Electron Volt γ Rays from the Radio Galaxy M87. Science 2006, 314, 1424–1427. [Google Scholar] [CrossRef] [Green Version]
  36. Broderick, A.E.; Loeb, A. Imaging bright-spots in the accretion flow near the black hole horizon of Sgr A*. Mon. Not. R. Astron. Soc. 2005, 363, 353–362. [Google Scholar] [CrossRef] [Green Version]
  37. Broderick, A.E.; Loeb, A. Imaging optically-thin hotspots near the black hole horizon of Sgr A* at radio and near-infrared wavelengths. Mon. Not. R. Astron. Soc. 2006, 367, 905–916. [Google Scholar] [CrossRef]
  38. Chael, A.A.; Johnson, M.D.; Narayan, R.; Doeleman, S.S.; Wardle, J.F.C.; Bouman, K.L. High-resolution Linear Polarimetric Imaging for the Event Horizon Telescope. Astrophys. J. 2016, 829, 11. [Google Scholar] [CrossRef]
  39. Harris, C.R.; Millman, K.J.; van der Walt, S.J.; Gommers, R.; Virtanen, P.; Cournapeau, D.; Wieser, E.; Taylor, J.; Berg, S.; Smith, N.J.; et al. Array programming with NumPy. Nature 2020, 585, 357–362. [Google Scholar] [CrossRef]
  40. Hunter, J.D. Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
  41. Stanzione, D.; West, J.; Evans, R.T.; Minyard, T.; Ghattas, O.; Panda, D.K. Frontera: The Evolution of Leadership Computing at the National Science Foundation. In Proceedings of the Practice and Experience in Advanced Research Computing, PEARC’20, Portland, OR, USA,, 26–30 July 2020; Association for Computing Machinery: New York, NY, USA, 2020; pp. 106–111. [Google Scholar]
  42. Leung, P.K.; Gammie, C.F.; Noble, S.C. Numerical Calculation of Magnetobremsstrahlung Emission and Absorption Coefficients. Astrophys. J. 2011, 737, 21. [Google Scholar] [CrossRef]
Figure 1. Magnetic flux eruptions can remove over half of the disk from near the black hole, significantly changing the resultant horizon-scale image. Top: Snapshots of simulation in quiescent state (left) and flux eruption event (right), ray-traced with equal mass scale units at a viewing inclination of 17 . Bottom: 3D rendering of gas density (in GRMHD code units) for the quiescent state (left) and flux eruption event (right) within the inner 15 M, viewed at 17 . The white region in the density plot shows the evacuation of the disk and the formation of a low density magnetospheric region near the black hole.
Figure 1. Magnetic flux eruptions can remove over half of the disk from near the black hole, significantly changing the resultant horizon-scale image. Top: Snapshots of simulation in quiescent state (left) and flux eruption event (right), ray-traced with equal mass scale units at a viewing inclination of 17 . Bottom: 3D rendering of gas density (in GRMHD code units) for the quiescent state (left) and flux eruption event (right) within the inner 15 M, viewed at 17 . The white region in the density plot shows the evacuation of the disk and the formation of a low density magnetospheric region near the black hole.
Galaxies 10 00107 g001
Figure 2. PFR in half-disk toy model, along with a schematic depicting the disk orientation on the right. Snapshots are for “static” disk (top) and rotating disk (bottom), both with θ cam = 80 . Both curves show a steep change in the PFR with ϕ cam because of gravitational lensing, but the addition of rotation introduces Doppler effects that modulate the relative flux, spreading the curve horizontally. The Doppler effects are also seen to a small degree for the “static” case because of the angular velocity of the zero-angular momentum frame.
Figure 2. PFR in half-disk toy model, along with a schematic depicting the disk orientation on the right. Snapshots are for “static” disk (top) and rotating disk (bottom), both with θ cam = 80 . Both curves show a steep change in the PFR with ϕ cam because of gravitational lensing, but the addition of rotation introduces Doppler effects that modulate the relative flux, spreading the curve horizontally. The Doppler effects are also seen to a small degree for the “static” case because of the angular velocity of the zero-angular momentum frame.
Galaxies 10 00107 g002
Figure 3. GRMHD PFR’s after removing magnetic field directional dependence. Snapshots (bottom) are for θ cam = 80 . Here, “quiescent” and “flux eruption” refer to time slices t = 8858 M and t = 9553 M respectively. During the flux eruption, the shape of the high-inclination ( θ cam = 80 ) PFR curve closely resembles that of the toy model, driven by the effects of Doppler boosting and magnification. Furthermore, the direct and indirect images of the flux eruption can each be seen in the snapshots below.
Figure 3. GRMHD PFR’s after removing magnetic field directional dependence. Snapshots (bottom) are for θ cam = 80 . Here, “quiescent” and “flux eruption” refer to time slices t = 8858 M and t = 9553 M respectively. During the flux eruption, the shape of the high-inclination ( θ cam = 80 ) PFR curve closely resembles that of the toy model, driven by the effects of Doppler boosting and magnification. Furthermore, the direct and indirect images of the flux eruption can each be seen in the snapshots below.
Galaxies 10 00107 g003
Figure 4. PFR curves for GRMHD ray-traced with full synchrotron emissivities (i.e., including magnetic field dependence). The results are similar to Figure 3, showing that the effects of magnetic field direction are insignificant in shaping the PFR for this example both in the quiescent state and during the flux eruption.
Figure 4. PFR curves for GRMHD ray-traced with full synchrotron emissivities (i.e., including magnetic field dependence). The results are similar to Figure 3, showing that the effects of magnetic field direction are insignificant in shaping the PFR for this example both in the quiescent state and during the flux eruption.
Galaxies 10 00107 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gelles, Z.; Chatterjee, K.; Johnson, M.; Ripperda, B.; Liska, M. Relativistic Signatures of Flux Eruption Events near Black Holes. Galaxies 2022, 10, 107. https://doi.org/10.3390/galaxies10060107

AMA Style

Gelles Z, Chatterjee K, Johnson M, Ripperda B, Liska M. Relativistic Signatures of Flux Eruption Events near Black Holes. Galaxies. 2022; 10(6):107. https://doi.org/10.3390/galaxies10060107

Chicago/Turabian Style

Gelles, Zachary, Koushik Chatterjee, Michael Johnson, Bart Ripperda, and Matthew Liska. 2022. "Relativistic Signatures of Flux Eruption Events near Black Holes" Galaxies 10, no. 6: 107. https://doi.org/10.3390/galaxies10060107

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop