Search for: All records

Abstract Accretion disks around both stellar-mass and supermassive black holes (BHs) are likely often warped. Whenever a disk is warped, its scale height varies with azimuth. Sufficiently strong warps cause extreme compressions of the scale height, which fluid parcels “bounce” off of twice per orbit to high latitudes. We study the dynamics of strong warps using (i) the nearly analytic “ring theory” of Fairbairn & Ogilvie, which we generalize to the Kerr metric, and (ii) three-dimensional general-relativistic hydrodynamic simulations of tori (“rings”) around BHs, using theH-AMRcode. We initialize a ring with a warp and study its evolution on tens of orbital periods. The simulations agree excellently with the ring theory until the warp amplitude,ψ, reaches a critical valueψ_c. Whenψ > ψ_c, the rings enter the bouncing regime. We analytically derive (and numerically validate) that ${\psi }_{\mathrm{c}}\approx {(r/r_{\mathrm{g}})}^{-1/2}$in the non-Keplerian regime, wherer_g = GM/c²is the gravitational radius, andMis the mass of the central object. Whenever the scale height bounces, the vertical velocity becomes supersonic, leading to “nozzle shocks” as gas collides at the scale height minima. Nozzle shocks damp the warp within ≈10–20 orbits, which is not captured by the ring theory. Nozzle shock dissipation leads to inflow timescales 1–2 orders of magnitude shorter than unwarpedαdisks, which may result in rapid variability, such as in changing-look active galactic nuclei or in the soft state of X-ray binaries. We propose that steady disks with strong warps may self-regulate to have amplitudes nearψ_c. 

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $$t_{\rm gen}$$, and one for the large-scale scale advection, $$t_{\rm adv}$$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $$\alpha$$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $$\alpha$$-dynamos. 

Abstract The conventional accretion disk lore is that magnetized turbulence is the principal angular momentum transport process that drives accretion. However, when dynamically important large-scale magnetic fields thread an accretion disk, they can produce mass and angular momentum outflows, known as winds,that also drive accretion. Yet, the relative importance of turbulent and wind-driven angular momentum transport is still poorly understood. To probe this question, we analyze a long-duration (1.2 × 10⁵r_g/c) simulation of a rapidly rotating (a= 0.9) black hole feeding from a thick (H/r∼ 0.3), adiabatic, magnetically arrested disk (MAD), whose dynamically important magnetic field regulates mass inflow and drives both uncollimated and collimated outflows (i.e., winds and jets, respectively). By carefully disentangling the various angular momentum transport processes within the system, we demonstrate the novel result that disk winds and disk turbulence both extract roughly equal amounts of angular momentum from the disk. We find cumulative angular momentum and mass accretion outflow rates of $\dot{L}\propto r^{0.9}$and $\dot{M}\propto r^{0.4}$, respectively. This result suggests that understanding both turbulent and laminar stresses is key to understanding the evolution of systems where geometrically thick MADs can occur, such as the hard state of X-ray binaries, low-luminosity active galactic nuclei, some tidal disruption events, and possibly gamma-ray bursts. 

Abstract Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius,R_B= 10³R_g. The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outsideR_Band suppress BH accretion to 1.5% of the Bondi rate, ${\dot{M}}_{\mathrm{B}}$. Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% of ${\dot{M}}_{\mathrm{B}}$reach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2% ${\dot{M}}_{\mathrm{B}}$. Because the BAD and RAD states tangle up the jets and destroy them well insideR_B, they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies. 

Abstract Accreting supermassive black holes (SMBHs) produce highly magnetized relativistic jets that tend to collimate gradually as they propagate outward. However, recent radio interferometric observations of the 3C 84 galaxy reveal a stunning, cylindrical jet already at several hundred SMBH gravitational radii,r≳ 350r_g. We explore how such extreme collimation emerges via a suite of 3D general relativistic magnetohydrodynamic simulations. We consider an SMBH surrounded by a magnetized torus immersed in a constant-density ambient medium that starts at the edge of the SMBH sphere of influence, chosen to be much larger than the SMBH gravitational radius,r_B= 10³r_g. We find that radiatively inefficient accretion flows (e.g., M87) produce winds that collimate the jets into parabolas near the black hole. After the disk winds stop collimating the jets atr≲r_B, they turn conical. Once outsider_B, the jets run into the ambient medium and form backflows that collimate the jets into cylinders some distance beyondr_B. Interestingly, for radiatively efficient accretion, as in 3C 84, the radiative cooling saps the energy out of the disk winds; at early times, they cannot efficiently collimate the jets, which skip the initial parabolic collimation stage, start out conical near the SMBH, and turn into cylinders already atr≃ 300r_g, as observed in 3C 84. Over time, the jet power remains approximately constant, whereas the mass accretion rate increases; the winds grow in strength and start to collimate the jets, which become quasi-parabolic near the base, and the transition point to a nearly cylindrical jet profile moves outward while remaining insider_B. 

ABSTRACT In general relativistic magnetohydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes in very high resolution GRMHD simulations. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating γ ≳ 100 (Te ≳ 1011 K) electrons. For models in which the electrons are heated to Te ∼ 1011 K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multiwavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization. 

Abstract The ongoing LIGO–Virgo–KAGRA observing run O4 provides an opportunity to discover new multimessenger events, including binary neutron star (BNS) mergers such as GW170817 and the highly anticipated first detection of a multimessenger black hole–neutron star (BH–NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH–NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH–NS mergers. Adopting optimal binary properties, a mass ratio ofq= 2, and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary’s evolution from premerger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of ∼−15 a few hours after the merger, longer than previous estimates, which did not consider the first principles–based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of ≳1 Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within ∼200 Mpc. The BH–NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals. 

Abstract A black hole (BH) traveling through a uniform, gaseous medium is described by Bondi–Hoyle–Lyttleton (BHL) accretion. If the medium is magnetized, then the black hole can produce relativistic outflows. We performed the first 3D, general-relativistic magnetohydrodynamic simulations of BHL accretion onto rapidly rotating black holes using theH-AMRcode, where we mainly varied the strength of a background magnetic field that threads the medium. We found that the ensuing accretion continuously drags the magnetic flux to the BH, which accumulates near the event horizon until it becomes dynamically important. Depending on the strength of the background magnetic field, the BHs can sometimes launch relativistic jets with high enough power to drill out of the inner accretion flow, become bent by the headwind, and escape to large distances. For stronger background magnetic fields, the jets are continuously powered, while at weaker field strengths they are intermittent, turning on and off depending on the fluctuating gas and magnetic flux distributions near the event horizon. We find that our jets reach extremely high efficiencies of ∼100%–300%, even in the absence of an accretion disk. We also calculated the drag forces exerted by the gas onto to the BH and found that the presence of magnetic fields causes the drag forces to be much less efficient than in unmagnetized BHL accretion. They can even sometimes become negative, accelerating the BH rather than slowing it down. Our results extend classical BHL accretion to rotating BHs moving through magnetized media, and demonstrate that accretion and drag are significantly altered in this environment. 

Abstract The angular momentum of gas feeding a black hole (BH) may be misaligned with respect to the BH spin, resulting in a tilted accretion disk. Rotation of the BH drags the surrounding spacetime, manifesting as Lense–Thirring torques that lead to disk precession and warping. We study these processes by simulating a thin (H/r= 0.02), highly tilted ( $\mathit{}=65°$) accretion disk around a rapidly rotating (a= 0.9375) BH at extremely high resolutions, which we performed using the general-relativistic magnetohydrodynamic codeH-AMR. The disk becomes significantly warped and continuously tears into two individually precessing subdisks. We find that mass accretion rates far exceed the standardα-viscosity expectations. We identify two novel dissipation mechanisms specific to warped disks that are the main drivers of accretion, distinct from the local turbulent stresses that are usually thought to drive accretion. In particular, we identify extreme scale height oscillations that occur twice an orbit throughout our disk. When the scale height compresses, “nozzle” shocks form, dissipating orbital energy and driving accretion. Separate from this phenomenon, there is also extreme dissipation at the location of the tear. This leads to the formation of low-angular momentum “streamers” that rain down onto the inner subdisk, shocking it. The addition of low-angular momentum gas to the inner subdisk causes it to rapidly accrete, even when it is transiently aligned with the BH spin and thus unwarped. These mechanisms, if general, significantly modify the standard accretion paradigm. Additionally, they may drive structural changes on much shorter timescales than expected inα-disks, potentially explaining some of the extreme variability observed in active galactic nuclei. 

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.