Search for: All records

Abstract A spinning black hole (BH) accreting from a disk of strongly magnetized plasma via a magnetically arrested disk is known to produce an efficient electromagnetic jet powered by the BH’s spin energy. We present general relativistic radiative magnetohydrodynamic simulations of magnetically arrested systems covering a range of sub- to super-Eddington accretion rates. Using the numerical results from these simulations, we develop formulae to describe the magnetization, jet efficiency, and spin evolution of an accreting BH as a function of its spin and accretion rate. A BH with near-Eddington accretion experiences a mild degree of spin-down because of angular momentum loss through the jet, leading to an equilibrium spin of 0.8 rather than 1.0 at the Eddington limit. As the accretion rate increases above Eddington, the spin-down effect becomes progressively stronger, ultimately converging on previous predictions based on nonradiative simulations. In particular, spin evolution drives highly super-Eddington systems toward a BH spin near zero. The formulae developed in this letter may be applied to galaxy- and cosmological-scale simulations that include BHs. If magnetically arrested disk accretion is common among supermassive BHs, the present results have broad implications for active galactic nucleus feedback and cosmological spin evolution. 

ABSTRACT We present two general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of magnetically arrested discs (MADs) around non-spinning (a* = 0) and spinning (a* = 0.9) supermassive black holes (BHs). In each simulation, the mass accretion rate is decreased with time such that we sample Eddington-scaled rates over the range $$3 \gtrsim \dot{M}/\dot{M}_{\rm {Edd}}\gtrsim 0.3$$. For the non-spinning BH model, the total and radiative efficiencies increase as the accretion rate decreases, varying over the range $$\eta _{\rm {tot}}\sim 9\!-\!16{{\ \rm per\ cent}}$$ and $$\eta _{\rm {rad}}\sim 6{-}12{{\ \rm per\ cent}}$$, respectively. This model shows very little jet activity. In contrast, the spinning BH model has a strong relativistic jet powered by spin energy extracted from the BH. The jet power declines with accretion rate such that $$\eta _{\rm {jet}}\sim 18{-}39{{\ \rm per\ cent}}$$ while the total and radiative efficiencies are $$\eta _{\rm {tot}}\sim 64{-}100{{\ \rm per\ cent}}$$ and $$\eta _{\rm {rad}}\sim 45{-}79{{\ \rm per\ cent}}$$, respectively. We confirm that mildly sub-Eddington discs can extract substantial power from a spinning BH, provided they are in the MAD state. The jet profile out to $$100\, GM/c^2$$ is roughly parabolic with a power-law index of k ≈ 0.43−0.53 during the sub-Eddington evolution. Both models show significant variability in the outgoing radiation which is likely associated with episodes of magnetic flux eruptions. The a* = 0.9 model shows semiregular variations with a period of $$\sim 2000\, GM/c^3$$ over the final $$\sim 10\, 000\, GM/c^3$$ of the simulation, which suggests that magnetic flux eruptions may be an important source of quasi-periodic variability. For the simulated accretion rates, the a* = 0 model is spinning up while the a* = 0.9 model is spinning down. Spinup–spindown equilibrium of the BH will likely be achieved at 0.5 < a*, eq < 0.6, assuming continuous accretion in the MAD state. 

Tidal disruption events (TDEs) around supermassive black holes (SMBHs) are a potential laboratory to study super-Eddington accretion disks and sometimes result in powerful jets or outflows which may shine in the radio and sub-millimeter bands. In this work, we modeled the thermal synchrotron emission of jets by general relativistic radiation magneto-hydrodynamics (GRRMHD) simulations of a BH accretion disk/jet system which assumed the TDE resulted in a magnetized accretion disk around a BH accreting at ∼12–25 times the Eddington accretion rate. Through synthetic observations with the Next Generation Event Horizon Telescope (ngEHT) and an image reconstruction analysis, we demonstrate that TDE jets may provide compelling targets within the context of the models explored in this work. In particular, we found that jets launched by a SANE super-Eddington disk around a spin a*=0.9 reach the ngEHT detection threshold at large distances (up to 100 Mpc in this work). A two-temperature plasma in the jet or weaker jets, such as a spin a*=0 model, requires a much closer distance, as we demonstrate detection at 10 Mpc for limiting cases of a*=0,R=1 or a*=0.9,R=20. We also demonstrate that TDE jets may appear as superluminal sources if the BH is rapidly rotating and the jet is viewed nearly face on. 

ABSTRACT We present general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of super-Eddington accretion flows around supermassive black holes (SMBHs), which may apply to tidal disruption events (TDEs). We perform long duration ($$t\ge 81,200\, GM/c^3$$) simulations that achieve mass accretion rates ≳11 times the Eddington rate and produce thermal synchrotron spectra and images of their jets. Gas flowing beyond the funnel wall expands conically and drives a strong shock at the jet head while variable mass ejection and recollimation, along the jet axis, results in internal shocks and dissipation. Assuming the ion temperature (Ti) and electron temperature (Te) in the plasma are identical, the radio/submillimetre spectra peak at >100 GHz and the luminosity increases with BH spin, exceeding $$\sim 10^{41} \, \rm {erg\, s^{-1}}$$ in the brightest models. The emission is extremely sensitive to Ti/Te as some models show an order-of-magnitude decrease in the peak frequency and up to four orders-of-magnitude decline in their radio/submillimetre luminosity as Ti/Te approaches 20. Assuming a maximum VLBI baseline distance of 10 Gλ, 230 GHz images of Ti/Te = 1 models shows that the jet head may be bright enough for its motion to be captured with the EHT (ngEHT) at D ≲ 110 (180) Mpc at the 5σ significance level. Resolving emission from internal shocks requires D ≲ 45 Mpc for both the EHT or ngEHT. 

ABSTRACT We use the general relativistic radiation magnetohydrodynamics code KORAL to simulate the accretion disc formation resulting from the tidal disruption of a solar mass star around a supermassive black hole (BH) of mass 106 M⊙. We simulate the disruption of artificially more bound stars with orbital eccentricity e ≤ 0.99 (compared to the more realistic case of parabolic orbits with e = 1) on close orbits with impact parameter β ≥ 3. We use a novel method of injecting the tidal stream into the domain, and we begin the stream injection at the peak fallback rate in this study. For two simulations, we choose e = 0.99 and inject mass at a rate that is similar to parabolic TDEs. We find that the disc only becomes mildly circularized with eccentricity e ≈ 0.6 within the 3.5 d that we simulate. The rate of circularization is faster for pericenter radii that come closer to the BH. The emitted radiation is mildly super-Eddington with $$L_{\rm {bol}}\approx 3{-}5\, L_{\rm {Edd}}$$ and the photosphere is highly asymmetric with the photosphere being significantly closer to the inner accretion disc for viewing angles near pericenter. We find that soft X-ray radiation with Trad ≈ 3–5 × 105 K may be visible for chance viewing angles. Our simulations suggest that TDEs should be radiatively inefficient with η ≈ 0.009–0.014. 

ABSTRACT We present the results of nine simulations of radiatively inefficient magnetically arrested discs (MADs) across different values of the black hole spin parameter a*: −0.9, −0.7, −0.5, −0.3, 0, 0.3, 0.5, 0.7, and 0.9. Each simulation was run up to $$t \gtrsim 100\, 000\, GM/c^3$$ to ensure disc inflow equilibrium out to large radii. We find that the saturated magnetic flux level, and consequently also jet power, of MAD discs depends strongly on the black hole spin, confirming previous results. Prograde discs saturate at a much higher relative magnetic flux and have more powerful jets than their retrograde counterparts. MADs with spinning black holes naturally launch jets with generalized parabolic profiles whose widths vary as a power of distance from the black hole. For distances up to 100GM/c2, the power-law index is k ≈ 0.27–0.42. There is a strong correlation between the disc–jet geometry and the dimensionless magnetic flux, resulting in prograde systems displaying thinner equatorial accretion flows near the black hole and wider jets, compared to retrograde systems. Prograde and retrograde MADs also exhibit different trends in disc variability: accretion rate variability increases with increasing spin for a* > 0 and remains almost constant for a* ≲ 0, while magnetic flux variability shows the opposite trend. Jets in the MAD state remove more angular momentum from black holes than is accreted, effectively spinning down the black hole. If powerful jets from MAD systems in Nature are persistent, this loss of angular momentum will notably reduce the black hole spin over cosmic time. 

We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models, including a constant electron to magnetic pressure (constant βe model) with a population of non-thermal electrons as well as an R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations are enhanced with every pair added. On the contrary, we show that, at higher frequencies, a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of a large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function. 

The Event Horizon Telescope (EHT) observation of M87^∗in 2018 has revealed a ring with a diameter that is consistent with the 2017 observation. The brightest part of the ring is shifted to the southwest from the southeast. In this paper, we provide theoretical interpretations for the multi-epoch EHT observations for M87^∗by comparing a new general relativistic magnetohydrodynamics model image library with the EHT observations for M87^∗in both 2017 and 2018. The model images include aligned and tilted accretion with parameterized thermal and nonthermal synchrotron emission properties. The 2018 observation again shows that the spin vector of the M87^∗supermassive black hole is pointed away from Earth. A shift of the brightest part of the ring during the multi-epoch observations can naturally be explained by the turbulent nature of black hole accretion, which is supported by the fact that the more turbulent retrograde models can explain the multi-epoch observations better than the prograde models. The EHT data are inconsistent with the tilted models in our model image library. Assuming that the black hole spin axis and its large-scale jet direction are roughly aligned, we expect the brightest part of the ring to be most commonly observed 90 deg clockwise from the forward jet. This prediction can be statistically tested through future observations.