Supermassive Black Hole Mergers Research Articles

Viscous torque in turbulent magnetized active galactic nucleus accretion disks and its effects on the gravitational waves of extreme mass ratio inspirals

The merger of supermassive black holes produces millihertz gravitational waves (GWs), which are potentially detectable by the future Laser Interferometer Space Antenna (LISA). Such binary systems are usually embedded in an accretion disk environment at the center of an active galactic nucleus (AGN). Recent studies suggest the plasma environment imposes measurable imprints on the GW signal if the mass ratio of the binary is around q $ $. The effect of the gaseous environment on the GW signal is strongly dependent on the disk’s parameters; therefore, it is believed that future low-frequency GW detections will provide us with precious information about the physics of AGN accretion disks. We investigated this effect by measuring the viscous torque via modeling of the evolution of magnetized tori around the primary massive black hole. Using the general relativistic magnetohydrodynamic HARM-COOL code, we performed 2D and 3D simulations of weakly magnetized, thin accretion disks, with a possible truncation and transition to advection-dominated accretion flow. We studied the angular momentum transport and turbulence generated by the magnetorotational instability. We quantified the disk’s effective alpha viscosity and its evolution over time. We applied our numerical results to quantify the relativistic viscous torque on a hypothetical low-mass secondary black hole via a 1D analytical approach, and we estimated the GW phase shift due to the gas environment.

Supermassive primordial black holes from inflation

There is controversy surrounding the origin and evolution of our universe's largest supermassive black holes (SMBHs). In this study, we consider the possibility that some of these black holes formed from the direct collapse of primordial density perturbations. Since the mass of a primordial black hole is limited by the size of the cosmological horizon at the time of collapse, these SMBHs must form rather late, and are naively in conflict with constraints from CMB spectral distortions. These limits can be avoided, however, if the distribution of primordial curvature perturbations is highly non-Gaussian. After quantifying the departure from Gaussianity needed to evade these bounds, we explore a model of multi-field inflation — a non-minimal, self-interacting curvaton model — which has all the necessary ingredients to yield such dramatic non-Gaussianities. We leave the detailed model building and numerics to a future study, however, as our goal is to highlight the challenges associated with forming SMBHs from direct collapse and to identify features that a successful model would need to have. This study is particularly timely in light of recent observations of high-redshift massive galaxy candidates by the James Webb Space Telescope as well as evidence from the NANOGrav experiment for a stochastic gravitational wave background consistent with SMBH mergers.

Primordial Black Holes from Spatially Varying Cosmological Constant Induced by Field Fluctuations in Extra Dimensions

The origin and evolution of supermassive black holes (SMBHs) in our universe have sparked controversy. In this study, we explore the hypothesis that some of these black holes may have seeded from the direct collapse of dark energy domains with density significantly higher than the surrounding regions. The mechanism of the origin of such domains relies on the inflationary evolution of a scalar field acting in D dimensions, which is associated with the cosmological constant in our four-dimensional spacetime manifold. Inner space quantum fluctuations of the field during inflation are responsible for the spatial variations of the dark energy density in our space. This finding holds particular significance, especially considering recent evidence from pulsar timing array observations, which supports the existence of a stochastic gravitational wave background consisting of SMBH mergers.

RABBITS – I. The crucial role of nuclear star formation in driving the coalescence of supermassive black hole binaries

ABSTRACT In this study of the ‘Resolving supermAssive Black hole Binaries In galacTic hydrodynamical Simulations’ (RABBITS) series, we focus on the hardening and coalescing process of supermassive black hole (SMBH) binaries in galaxy mergers. For simulations including different galaxy formation processes (i.e. gas cooling, star formation, SMBH accretion, stellar, and AGN feedback), we systematically control the effect of stochastic eccentricity by fixing it to similar values during the SMBH hardening phase. We find a strong correlation between the SMBH merger time-scales and the presence of nuclear star formation. Throughout the galaxy merging process, gas condenses at the centre due to cooling and tidal torques, leading to nuclear star formation. These recently formed stars, which inherit low angular momenta from the gas, contribute to the loss cone and assist in the SMBH hardening via three-body interactions. Compared to non-radiative hydrodynamical runs, the SMBH merger time-scales measured from the runs including cooling, stellar, and SMBH physical processes tend to be shortened by a factor of ∼1.7. After fixing the eccentricity to the range of e ∼ 0.6–0.8 during the hardening phase, the simulations with AGN feedback reveal merger time-scales of ∼100–500 Myr for disc mergers and ∼1–2 Gyr for elliptical mergers. With a semi-analytical approach, we find that the torque interaction between the binary and its circumbinary disc has minimal impact on the shrinking of the binary orbit in our retrograde galaxy merger. Our results are useful in improving the modelling of SMBH merger time-scales and gravitational-wave event rates.

Analysis of Kozai cycles in equal-mass hierarchical triple supermassive black hole mergers in the presence of a stellar cluster

ABSTRACT Supermassive black holes (SMBHs) play an important role in galaxy evolution. Binary and triple SMBHs can form after galaxy mergers. A third SMBH may accelerate the SMBH merging process, possibly through the Kozai mechanism. We use N-body simulations to analyse oscillations in the orbital elements of hierarchical triple SMBHs with surrounding star clusters in galaxy centres. We find that SMBH triples spend only a small fraction of time in the hierarchical merger phase (i.e. a binary SMBH with a distant third SMBH perturber). Most of the time, the enclosed stellar mass within the orbits of the innermost or the outermost SMBH is comparable to the SMBH masses, indicating that the influence of the surrounding stellar population cannot be ignored. We search for Eccentric Kozai–Lidov (EKL) oscillations for which (i) the eccentricity of the inner binary and inclination are both oscillate and are antiphase or in-phase and (ii) the oscillation period is consistent with EKL time-scale. We find that EKL oscillations are short-lived and rare: the triple SMBH spends around 3 per cent of its time in this phase over the ensemble of simulations, reaching around 8 per cent in the best-case scenario. This suggests that the role of the EKL mechanism in accelerating the SMBH merger process may have been overestimated in previous studies. We follow-up with three-body simulations, using initial conditions extracted from the simulation, and the result can to some extent repeat the observed EKL-like oscillations. This comparison provides clues about why those EKL oscillations with perturbing stars are short-lived.

Anisotropic Star Clusters around Recoiling Supermassive Black Holes

Gravitational-wave recoil kicks from merging supermassive black hole binaries can have a profound effect on the surrounding stellar population. In this work, we study the dynamic and kinematic properties of nuclear star clusters following a recoil kick. We show that these postkick structures present unique signatures that can provide key insight to observational searches for recoiling supermassive black holes. In our previous paper, we showed that an in-plane recoil kick turns a circular disk into a lopsided, eccentric disk such as the one we observe in the Andromeda nucleus. Building on this work, here we explore many recoil kick angles as well as initial stellar configurations. For a circular disk of stars, an in-plane kick causes strong apsidal alignment with a significant fraction of the disk becoming retrograde at large radii. If the initial orbits are highly eccentric, an in-plane kick forms a bar-like structure made up of two antialigned lopsided disks. An out-of-plane kick causes clustering in the argument of periapsis, ω, regardless of the initial eccentricity distribution. Initially, isotropic configurations form anisotropies in the form of a torus of eccentric orbits oriented perpendicular to the recoil kick. Postkick surface density and velocity maps are presented in each case to highlight the distinct, observable structures of these systems.

High-redshift supermassive black hole mergers in simulations with dynamical friction modelling

ABSTRACT In the near future, projects like Laser Interferometer Space Antenna (LISA) and pulsar timing arrays are expected to detect gravitational waves from mergers between supermassive black holes, and it is crucial to precisely model the underlying merger populations now to maximize what we can learn from this new data. Here, we characterize expected high-redshift (z > 2) black hole mergers using the very large volume Astrid cosmological simulation, which uses a range of seed masses to probe down to low-mass black holes (BHs), and directly incorporates dynamical friction so as to accurately model the dynamical processes that bring black holes to the galaxy centre where binary formation and coalescence will occur. The black hole populations in Astrid include black holes down to $\sim 10^{4.5} \, \mathrm{M}_\odot$, and remain broadly consistent with the TNG simulations at scales $\gt 10^6 \, \mathrm{M}_\odot$ (the seed mass used in TNG). By resolving lower mass black holes, the overall merger rate is ∼5× higher than in TNG. However, incorporating dynamical friction delays mergers compared to a recentring scheme, reducing the high-z merger rate mass-matched mergers by a factor of ∼2×. We also calculate the expected LISA signal-to-noise values, and show that the distribution peaks at high SNR (>100), emphasizing the importance of implementing a seed mass well below LISA’s peak sensitivity ($\sim 10^6 \, \mathrm{M}_\odot$) to resolve the majority of LISA’s gravitational wave detections.

A Candidate Dual QSO at Cosmic Noon

We report the discovery of a candidate dual QSO at z = 1.889, a redshift that is in the era known as “cosmic noon” where most of the universe’s black hole and stellar mass growth occurred. The source was identified in Hubble Space Telescope WFC3/IR images of a dust-reddened QSO that showed two closely separated point sources at a projected distance of 0.″26, or 2.2 kpc. This red QSO was targeted for imaging to explore whether red QSOs are hosted by merging galaxies. We subsequently obtained a spatially resolved Space Telescope Imaging Spectrograph spectrum of the system, covering the visible spectral range, and verifying the presence of two distinct QSO components. We also obtained high-resolution radio continuum observations with the Very Long Baseline Array at 1.4 GHz (21 cm L band) and found two sources coincident with the optical positions. The sources have similar black hole masses, bolometric luminosities, and radio-loudness parameters. However, their colors and reddenings differ significantly. The redder QSO has a higher Eddington ratio, consistent with previous findings. We consider the possibility of gravitational lensing and find that it would require extreme and unlikely conditions. If confirmed as a bona fide dual QSO, this system would link dust reddening to galaxy and supermassive black hole mergers, opening up a new population in which to search for samples of dual active galactic nuclei.

Aligning Retrograde Nuclear Cluster Orbits with an Active Galactic Nucleus Accretion Disc

ABSTRACT Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination i < 90°) decrease inclination, as well as semimajor axis (a) and eccentricity (e) until orbital alignment with the gas disc (‘disc capture’). Captured stellar-origin black holes (sBH) add to the embedded AGN population that drives sBH–sBH mergers detectable in gravitational waves using LIGO–Virgo–KAGRA or sBH–SMBH mergers detectable with Laser Interferometer Space Antenna. Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive ‘immortal’ stars. Here, we investigate the behaviour of polar and retrograde orbiters (i ≥ 90°) interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation (i < 90°) during capture, and decrease a dramatically towards the SMBH. For sBH, we find a critical angle iret ∼ 113°, below which retrograde sBH decay towards embedded prograde orbits (i → 0°), while for io > iret sBH decay towards embedded retrograde orbits (i → 180°). sBH near polar orbits (i ∼ 90°) and stars on nearly embedded retrograde orbits (i ∼ 180°) show the greatest decreases in a. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial a. For sBH, disc capture time is longest for polar orbits, low-mass sBH, and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the spin distribution of embedded sBH.

Searching for clues of past binary supermassive black hole mergers in nuclear star clusters

ABSTRACTGalaxy mergers are common processes in the Universe. As a large fraction of galaxies hosts at their centres a central supermassive black hole (SMBH), mergers can lead to the formation of a supermassive black hole binary (SMBHB). The formation of such a binary is more efficient when the SMBHs are embedded in a nuclear star cluster (NSC). NSCs are dense and massive stellar clusters present in the majority of the observed galaxies. Their central densities can reach up to $10^7\, {\rm M_{\odot }}\,{\rm pc^{-3}}$ and their masses can be as large as a few $10^7\, {\rm M_{\odot }}$. The direct detection of an SMBHB is observationally challenging. In this work, we illustrate how the large-scale structural and dynamical properties of an NSC can help to identify nucleated galaxies that recently went through a merger that possibly led to the formation of a central SMBHB. Our models show that the merger can imprint signatures on the shape, density profile, rotation, and velocity structure of the NSC. The strength of the signatures depends on the mass ratio between the SMBHs and on the orbital initial conditions of the merger. In addition, the number of hypervelocity stars produced in the mergers is linked to the SMBHB properties. The merger can also contribute to the formation of the nuclear stellar disc of the galaxy.

Chasing supermassive black hole merging events withAthenaandLISA

ABSTRACTThe European Space Agency is studying two large-class missions bound to operate in the decade of the 30s, and aiming at investigating the most energetic and violent phenomena in the Universe. Athena is poised to study the physical conditions of baryons locked in large-scale structures from the epoch of their formation, as well as to yield an accurate census of accreting supermassive black holes down to the epoch of reionization; LISA will extend the hunt for Gravitational Wave (GW) events to the hitherto unexplored mHz regime. We discuss in this paper the science that their concurrent operation could yield, and present possible Athena observational strategies. We focus on Supermassive (M$\lesssim 10^7\, \rm {M_\odot }$) Black Hole Mergers (SMBHMs), potentially accessible to Athena up to z ∼ 2. The simultaneous measurement of their electromagnetic (EM) and GW signals may enable unique experiments in the domains of astrophysics, fundamental physics, and cosmography, such as the magnetohydrodynamics of fluid flows in a rapidly variable space–time, the formation of coronae and jets in Active Galactic Nuclei, and the measurement of the speed of GW, among others. Key to achieve these breakthrough results will be the LISA capability of locating a SMBHM event with an error box comparable to, or better than the field-of-view of the Athena Wide Field Imager ($\simeq 0.4\,$ deg2) and Athena capability to slew fast to detect the source during the inspiral phase and the post-merger phase. Together, the two observatories will open in principle the exciting possibility of truly concurrent EM and GW studies of the SMBHMs

Modelling the accretion and feedback of supermassive black hole binaries in gas-rich galaxy mergers

ABSTRACTWe introduce a new model for the accretion and feedback of supermassive black hole (SMBH) binaries to the ketju code, which enables us to resolve the evolution of SMBH binaries down to separations of tens of Schwarzschild radii in gas-rich galaxy mergers. Our subgrid binary accretion model extends the widely used Bondi–Hoyle–Lyttleton accretion into the binary phase and incorporates preferential mass accretion on to the secondary SMBH, which is motivated by results from small-scale hydrodynamical circumbinary disc simulations. We perform idealized gas-rich disc galaxy merger simulations using pure thermal or pure kinetic active galactic nuclei (AGNs) feedback. Our binary accretion model provides more physically motivated SMBH mass ratios, which are one of the key parameters for computing gravitational wave (GW) induced recoil velocities. The merger time-scales of our simulated SMBH binaries are in the range tmerge ∼ 10–400 Myr. Prograde in-plane equal-mass galaxy mergers lead to the shortest merger time-scales, as they experience the strongest starbursts, with the ensuing high stellar density resulting in a rapid SMBH coalescence. Compared to the thermal AGN feedback, the kinetic AGN feedback predicts longer merger time-scales and results in more core-like stellar profiles, as it is more effective in removing gas from the galaxy centre and quenching star formation. This suggests that the AGN feedback implementation plays a critical role in modelling SMBH coalescences. Our model will be useful for improving the modelling of SMBH mergers in gas-rich galaxies, the prime targets for the upcoming LISA GW observatory.

Gravitational waves from domain walls in Pulsar Timing Array datasets

We present a model-independent search for a gravitational wave background from cosmic domain walls (DWs) in the NANOGrav 12.5 years dataset and International PTA Data Release 2. DWs that annihilate at temperatures ∼ 20–50 MeV with tensions ∼ (40–100 TeV)3 provide as good a fit to both datasets as the astrophysical background from supermassive black hole mergers. DWs may decay into the Standard Model (SM) or a dark sector. In the latter case we predict an abundance Δ N eff of dark radiation well within the reach of upcoming CMB surveys. Complementary signatures at colliders and laboratories can arise if couplings to the SM are present. As an example, we discuss heavy axion scenarios, where DW annihilation may interestingly be induced by QCD confinement.

A unique, ring-like radio source with quadrilateral structure detected with machine learning

ABSTRACT We report the discovery of a unique object in the MeerKAT Galaxy Cluster Legacy Survey (MGCLS) using the machine learning anomaly detection framework astronomaly. This strange, ring-like source is 30′ from the MGCLS field centred on Abell 209, and is not readily explained by simple physical models. With an assumed host galaxy at redshift 0.55, the luminosity (1025 W Hz−1) is comparable to powerful radio galaxies. The source consists of a ring of emission 175 kpc across, quadrilateral enhanced brightness regions bearing resemblance to radio jets, two ‘ears’ separated by 368 kpc, and a diffuse envelope. All of the structures appear spectrally steep, ranging from −1.0 to −1.5. The ring has high polarization (25 per cent) except on the bright patches (<10 per cent). We compare this source to the Odd Radio Circles recently discovered in ASKAP data and discuss several possible physical models, including a termination shock from starburst activity, an end-on radio galaxy, and a supermassive black hole merger event. No simple model can easily explain the observed structure of the source. This work, as well as other recent discoveries, demonstrates the power of unsupervised machine learning in mining large data sets for scientifically interesting sources.

Ultramassive Black Holes Formed by Triple Quasar Mergers at z ∼ 2

The origin of rare and elusive ultramassive black holes (UMBH; with M BH > 1010 M ⊙) is an open question. Using the large volume cosmological hydrodynamic simulation ASTRID, we report on the formation of an extremely massive UMBH with M BH ∼ 1011 M ⊙ at z ∼ 2. The UMBH is assembled as a result of two successive mergers of massive galaxies each with stellar mass M * > 3 × 1011 M ⊙ that also produces a bright, rare triple quasar system powered by three ∼109 M ⊙ black holes. The second merger of supermassive black holes (SMBHs) follows the first after 150 Myr. The merger events lead to sustained Eddington accretion onto the central SMBH, forming a UMBH in the center of a massive compact stellar core with M * > 2 × 1012 M ⊙. The strong feedback of the UMBH quenches the surrounding star formation to <10 M ⊙ yr−1 in the inner 50 h −1 kpc region. There are two more UMBHs with M BH > 5 × 1010 M ⊙ at z > 2 in ASTRID that are also produced by major mergers of galaxies, and their progenitors can be observed as quasar triplets of lower luminosity. The rarely observed quasar multiples can be the cradle of UMBHs at high redshift, and likely end up in the center of the most massive clusters.

Extragalactic neutrino-emission induced by supermassive and stellar mass black hole mergers

ABSTRACT The recent detections of binary stellar mass black hole mergers by the LIGO and Virgo Collaborations suggest that such mergers are common occurrences. Galaxy mergers further indicate that supermassive black holes in centres of galaxies also merge and are typically expected to have had at least one merger in their lifetime, possibly many. In the presence of a jet, these mergers are almost always accompanied by a change of the jet direction and a connected jet precession motion, leading to interactions of the jet with ambient matter and producing high-energy particles, and consequently high-energy gamma-rays and neutrinos. In this work, we investigate the possibility under which conditions such mergers could be the sources of the diffuse astrophysical neutrino flux measured by the IceCube Neutrino Observatory. The main free parameters in the calculation concern the frequency of the mergers and the fraction of energy that is transferred from the gravitationally released energy to neutrinos. We show that the merger rate for SMBBHs must lie between ∼10−7 and 10−5 Gpc−3 yr−1. The ratio of energy going to neutrinos during such mergers lies then between ∼10−6 and 3 × 10−4. For stellar mass BBH mergers, the rate needs to be ∼10–100 Gpc−3 yr−1 and the expected ratio of neutrino to GW energy lies in a comparable range as for SMBBHs, ∼2 × 10−5–10−3. These values lie in a reasonable parameter range, so that the production of neutrinos at the level of the detected neutrino flux is a realistic possibility.

Mergers of Supermassive and Intermediate-mass Black Holes in Galactic Nuclei from Disruptions of Star Clusters

Gravitational waves (GWs) offer an unprecedented opportunity to survey the sky and detect mergers of compact objects. While intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt with either dynamical or accretion signatures, the GW landscape appears very promising. Mergers of an IMBH with a supermassive black hole (SMBH) will be primary sources for the planned space-based mission LISA and could be observed up to the distant universe. SMBH–IMBH binaries can be formed as a result of the migration and merger of stellar clusters at the center of galaxies, where an SMBH lurks. We build for the first time a semianalytical framework to model this scenario and find that the comoving merger rate of SMBH–IMBH binaries is ∼10−4 Gpc−3 yr−1 in the local universe for a unity IMBH occupation fraction, scales linearly with it, and has a peak at z ≈ 0.5–2. Our model predicts ∼0.1 events yr−1 within redshift z ≈ 3.5 if 10% of the inspiraled star clusters hosted an IMBH, while ∼1 event yr−1 for a unity occupation fraction. More than 90% of these systems will be detectable with LISA with a signal-to-noise ratio larger than 10, promising to potentially find a family of IMBHs.

Trinity I: self-consistently modelling the dark matter halo–galaxy–supermassive black hole connection from z = 0–10

ABSTRACT We present Trinity, a flexible empirical model that self-consistently infers the statistical connection between dark matter haloes, galaxies, and supermassive black holes (SMBHs). Trinity is constrained by galaxy observables from 0 &amp;lt; z &amp;lt; 10 [galaxies’ stellar mass functions, specific and cosmic star formation rates (SFRs), quenched fractions, and UV luminosity functions] and SMBH observables from 0 &amp;lt; z &amp;lt; 6.5 (quasar luminosity functions, quasar probability distribution functions, active black hole mass functions, local SMBH mass–bulge mass relations, and the observed SMBH mass distributions of high-redshift bright quasars). The model includes full treatment of observational systematics [e.g. active galactic nucleus (AGN) obscuration and errors in stellar masses]. From these data, Trinity infers the average SMBH mass, SMBH accretion rate, merger rate, and Eddington ratio distribution as functions of halo mass, galaxy stellar mass, and redshift. Key findings include: (1) the normalization and the slope of the SMBH mass–bulge mass relation increases mildly from z = 0 to z = 10; (2) The best-fitting AGN radiative+kinetic efficiency is ∼0.05–0.06, but can be in the range ∼0.035–0.07 with alternative input assumptions; (3) AGNs show downsizing, i.e. the Eddington ratios of more massive SMBHs start to decrease earlier than those of lower mass objects; (4) The average ratio between average SMBH accretion rate and SFR is ∼10−3 for low-mass galaxies, which are primarily star-forming. This ratio increases to ∼10−1 for the most massive haloes below z ∼ 1, where star formation is quenched but SMBHs continue to accrete.

X-ray view of a merging supermassive black hole binary candidate SDSS J1430+2303: Results from the first ∼200 days of observations

Context. Recently we discovered an unprecedented supermassive black hole binary (SMBHB) candidate in the nearby Seyfert galaxy SDSS J1430+2303, which is predicted to merge within three years. X-ray spectroscopy may bring unique kinematic evidence for the last inspiraling stage, when the binary is too close to allow each of them to hold an individual broad line region. Aims. We try to confirm the unique SMBHB merger event and understand the associated high-energy processes from a comprehensive X-ray view. Methods. We observed SDSS J1430+2303 with XMM-Newton, NuSTAR, Chandra, and Swift spanning the first ∼200 days since its discovery. Results. X-ray variability, up to a factor of 7, has been detected on a timescale of a few days. The broadband spectrum from 0.2–70 keV can be well fitted with a model consisting of a power law and a relativistic reflection covered by a warm absorber. The properties of the warm absorber changed dramatically, for example, with a decrease in the line-of-sight velocity from ∼0.2c to ∼0.02c, between the two XMM-Newton observations separated by only 19 days, which can be naturally understood in the context of the SMBHB; although, the clumpy wind scenario cannot be completely excluded. Broad Fe Kα emission has been robustly detected, though its velocity shift or profile change is not yet measurable. Further longer X-ray observations are highly encouraged to detect the expected orbital motion of the binary.

Bridging the Bondi and Event Horizon Scales: 3D GRMHD Simulations Reveal X-shaped Radio Galaxy Morphology

X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up ≲10% of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole (SMBH), SMBH merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic magnetohydrodynamic (GRMHD) simulation for the first time follows magnetized rotating gas from outside the SMBH sphere of influence of radius R B to the SMBH of gravitational radius R g at the largest scale separation, R B/R g = 103, to date. Initially, our axisymmetric system of constant-density hot gas contains a weak vertical magnetic field and rotates in the equatorial plane of a rapidly spinning SMBH. We seed the gas with small-scale 2% level pressure perturbations. Infalling gas forms an accretion disk, and the SMBH launches relativistically magnetized collimated jets reaching well outside R B. Under the pressure of the infalling gas, the jets intermittently turn on and off, erratically wobble, and inflate pairs of cavities in different directions, resembling an X-shaped jet morphology. Synthetic X-ray images reveal multiple pairs of jet-powered shocks and cavities. Large-scale magnetic flux accumulates on the SMBH, becomes dynamically important, and leads to a magnetically arrested disk state. The SMBH accretes at 2% of the Bondi rate ( for M87*) and launches twin jets at η = 150% efficiency. These jets are powerful enough (P jets ≃ 2 × 1044 erg s−1) to escape along the SMBH spin axis and end the short-lived intermittent jet state, whose transient nature can account for the rarity of XRGs.