## Magnetically Arrested Circumbinary Accretion Flows

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00757v2-figure-1952a3df-90c7-4d7a-b868-5ec8f13e864a.jpg)

1 TAPIR, Mailcode 350-17, California Institute of Technology, Pasadena, CA 91125, USA 2

Walter Burke Institute for Theoretical Physics, California Institute of Technology, Pasadena, CA 91125, USA

## ABSTRACT

Binary systems with comparable masses and a surrounding accretion disk can accrete gas through spiral accretion streams penetrating the central cavity formed by tidal interactions. Using threedimensional Newtonian magnetohydrodynamics simulations, we investigate the possibility of a magnetically arrested accretion flow through the cavity. Rather than solely continuously feeding the binary through spiral accretion streams, the accretion is regulated by the strong magnetic field inside the cavity. Transport of mass and angular momentum onto the binary then proceeds largely periodically in magnetic flux eruption episodes. The ejected flux tubes carry angular momentum outwards and away from the binary, inject hot plasma into the disk and can launch flares. This likely intermittent scenario could have potential implications for the emission signatures of supermassive black hole binaries, and shed light onto the role magnetic fields play in the binary's orbital evolution.

## 1. INTRODUCTION

Binary systems are a prime focus of study in astronomy and astrophysics, including protoplanetary disks (Williams &amp; Cieza 2011), binary stars (Kratter et al. 2008; Matsumoto 2024), and supermassive black hole binaries (Begelman et al. 1980; Milosavljevi´ &amp; Merritt c 2001; Milosavljevi´ &amp; Phinney 2005; Dotti et al. 2007; c Mayer et al. 2007). Accretion onto the binary has the potential to power electromagnetic transients and alter the dynamics of the orbit (see, e.g., Lai &amp; Mu˜oz (2023) n for a recent review). In the case of supermassive black hole binaries, accretion onto the system may help explain hardening the binary beyond the parsec separation scale (Begelman et al. 1980). As such, understanding gas accretion onto a supermassive black hole binary is of crucial importance for understanding the merging black hole population seen with the Pulsar Timing Array (PTA) (Agazie et al. 2023a,b) and the Laser Interferometer Space Antenna (LISA) (Amaro-Seoane et al. 2023).

resonance (Artymowicz &amp; Lubow 1994, 1996; Goldreich &amp; Tremaine 1979)). One of the main features of circumbinary accretion is that accretion proceeds mainly along two spiral arms penetrating the low-density cavity. A main driver behind circumbinary disk research is to understand the amount of angular momentum transferred onto the binary, and in particular whether the binary hardens or widens due to gas accretion (see, e.g., Mu˜ noz et al. 2019; Lai &amp; Mu˜oz 2023 and references n therein), as well as to unravel potential electromagnetic counterparts (e.g., Bogdanovi´ et al. 2022; Charisi et al. c 2022).

For binaries surrounded by an accretion disk, the timevarying gravitational potential of the binary exerts a torque on the surrounding circumbinary disk, balancing viscous torques and opens up a cavity at 2-3 times the binary separation (near the location of the inner Lindblad emost@caltech.edu

haiyangw@caltech.edu

Over the past decade, several numerical studies have clarified many of the hydrodynamical aspects of circumbinary accretion, including the dependence on binary properties (mass ratio (Duffell et al. 2020; Derdzinski et al. 2021; Siwek et al. 2023a,b), orbital eccentricity (D'Orazio &amp; Duffell 2021; Siwek et al. 2023a,b)) and the properties of the disk (aspect ratio (Tiede et al. 2020; Dittmann &amp; Ryan 2024), viscosity (Dittmann et al. 2023; Dittmann &amp; Ryan 2024), cooling (Sudarshan et al. 2022; Wang et al. 2023a,b; Pierens &amp; Nelson 2023), selfgravity of the gas (Franchini et al. 2021; Bourne et al. 2023) ). The stage when gravitational wave emission is non-negligible has also been studied by considering relativistic precession (DeLaurentiis et al. 2024) and decoupling on the final orbits (Tsokaros et al. 2022; Dittmann et al. 2023).

Due to the long integration times required to reach a quasi-steady state of angular momentum transfer be-

tween the binary and the disk (Miranda et al. 2017; Mu˜ noz et al. 2019), most simulations have been carried out in two-dimensions focusing on the orbital plane. Few studies have investigated the hydrodynamial torque evolution in full three dimensions (Moody et al. 2019; Bourne et al. 2023).

On black hole horizon scales several works have also investigated the dynamics of the black hole-gas interaction in numerical or post-Newtonian relativity (Noble et al. 2012; Avara et al. 2023; Farris et al. 2012; Gold et al. 2014), including jet launching and mini-disks (Bowen et al. 2017, 2018; Armengol et al. 2021; Combi et al. 2022; Guti´ errez et al. 2022), which form around the black holes within their respective Hill spheres (Paschalidis et al. 2021; Bright &amp; Paschalidis 2023). On arbitrary - and in the supermassive black context (sub)parsec - scales few simulations have been carried focusing mainly on the outer disk, and its magnetic torque evolution (Shi et al. 2012; Noble et al. 2021; Avara et al. 2023). Qualitatively, the inner cavity was either excised or remained in a state resembling broadly previous hydrodynamical works.

On the other hand, magnetic fields in single accreting black holes have been shown to cause drastic changes to the accretion flows (Yuan &amp; Narayan 2014; Davis &amp; Tchekhovskoy 2020). In particular, inferences from observations with the Event Horizon Telescope indicate, that black holes such as M87* and Sgr A* (Akiyama et al. 2019, 2021, 2022; Collaboration et al. 2024) might be in a magnetically arrested (MAD) state (Igumenshchev et al. 2003; Narayan et al. 2003; Tchekhovskoy et al. 2011). In contrast to advection dominated standard and normal accretion (SANE) flows, magnetically arrested flows are halted by a strong magnetosphere around the black hole, essentially establishing a cavity stabilized by magnetic pressure, which the flow can only cross via interchange-like instabilities pushing the field locally into the accretion disk (Tchekhovskoy et al. 2011; Porth et al. 2021; Ripperda et al. 2022). Accretion then proceeds in episodic bursts, with the initial development of Rayleigh-Taylor fingers (Zhdankin et al. 2023) ( see also, e.g., Kulkarni &amp; Romanova 2008; Takasao et al. 2022 for similar dynamics in T Tauri stars), before full magnetic flux tubes are ejected into the disk, and cause an asymmetric burst of accretion on the black hole (Porth et al. 2021; Ripperda et al. 2022; Chatterjee &amp; Narayan 2022). This inflow of magnetic flux re-establishes a magnetically dominated zone, pushing out the disk and resetting the magnetically arrested accretion cycle (Tchekhovskoy et al. 2011; Ripperda et al. 2022). Most importantly, angular momentum transport onto the black hole then mainly happens during the flux eruptions, and itself is episodic, magnetically driven, and in parts outward-directed (Chatterjee &amp; Narayan 2022), and can spin down the black hole significantly (Lowell et al. 2024). At the same time, it has been suggested that magnetic flux eruptions into the disk are mainly loaded with non-thermal electronpositron pairs and can potentially power X-ray transient emission (Dexter et al. 2020; Scepi et al. 2022; Zhdankin et al. 2023; Hakobyan et al. 2023; Vos et al. 2023) with implications for X-ray flares observed in the galactic center (GRAVITY Collaboration et al. 2018, 2021). It is important to point out that the picture painted here has been mainly investigated for accreting black holes, although it is generic and has also been found for accreting young proto-stars (Kulkarni &amp; Romanova 2008; Takasao et al. 2022) and for low-mass neutron star X-ray binaries (Parfrey &amp; Tchekhovskoy 2023; Murguia-Berthier et al. 2024) (see also Arons &amp; Lea 1976; Scharlemann 1978; Spruit et al. 1995).

Motivated by recent progress in supermassive black hole binary formation (e.g., Shi et al. 2024) and findings that large-scale accretion flows on parsec scales may well be very strongly magnetically dominated (Hopkins et al. 2024; Hopkins 2024; Guo et al. 2024), it seems possible that a strongly magnetized state inside the circumbinary disk cavity could form, with accretion properties similar to magnetically arrested accretion flows on black hole horizon scales.

In this work, using massively parallel GPU-enabled simulations, we demonstrate that the picture of magnetically arrested accretion flows around black holes naturally carries over to circumbinary disks, albeit with modifications introduced by the strong tidal potential of the binary. In particular, we show that the circumbinary disk cavity can be rapidly filled with vertical magnetic flux, establishing a magnetically arrested accretion state. We further identify Rayleigh-Taylor fingers, flares, and flux eruptions, which we find to reshape the cavity size and transport most of the angular momentum between the binary inside the cavity and the outer circumbinary disk. We also provide a detailed analysis of the angular momentum transport between the disk and the binary, finding that the magnetically arrested state has the potential to lead to an overall shrinking of the orbit.

## 2. METHODS

We study the dynamics of a circumbinary accretion disk endowed with a strong magnetic field. To this end, we model the binary as Newtonian point masses on a

fixed Keplerian orbit, and the dynamics of the gas flow using ideal magnetohydrodynamics (Stone et al. 2008). The fluid is described by a mass density, ρ , velocity, v , internal energy density e , pressure P , and magnetic field, B . The evolution of the system then follows (Stone et al. 2020),

$$\frac { \partial \rho } { \partial t } + \nabla \cdot ( \rho v ) = s _ { \rho } \,, \quad \quad ( 1 )$$

$$\frac { \partial ( \rho v ) } { \partial t } + \nabla \cdot ( \rho v v + \bar { P } \mathbb { I } - B B ) = s _ { p } - \rho \nabla \Phi \,, \quad ( 2 ) \quad \text{po} \\ \infty \, \pi$$

$$\frac { \partial E } { \partial t } + \nabla \cdot [ ( E + \bar { P } ) v - ( B \cdot v ) \, B ] = s _ { E } - \rho v \cdot \nabla \Phi \,,$$

(3)

$$\frac { \partial B } { \partial t } - \nabla \times [ v \times B ] = 0 \,, \quad \quad ( 4 )$$

where we have defined the conserved energy, E , and effective pressure ¯ , P

$$E = e + \frac { 1 } { 2 } \rho v ^ { 2 } + \frac { 1 } { 2 } B ^ { 2 } \,, \quad \quad ( 5 ) \quad \quad \, \quad \, |$$

$$\bar { P } = P + \frac { \bar { 1 } } { 2 } B ^ { 2 } \,. \, \text{ \quad \ \ } ( 6 )$$

The gravitational potential of the binary is modeled as

$$\Phi = - \frac { G M _ { 1 } } { ( | r - r _ { 1 } | ^ { 2 } + \epsilon ^ { 2 } ) ^ { 1 / 2 } } - \frac { G M _ { 2 } } { ( | r - r _ { 2 } | ^ { 2 } + \epsilon ^ { 2 } ) ^ { 1 / 2 } } \,, \quad ( 7 ) \quad \text{ wit} \\. \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ddots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

where M i and r i are the mass and location of the binary constituents respectively, G is the gravitational constant, r is the radial coordinate to the center of mass of the binary, and ϵ is a small number. We further adopt a unit convention of GM = a = 1 in our simulations, and focus in an equal mass system, M 1 = M 2 .

Since we model the binary dynamics on scales of the orbital separation, a , which are much larger than the physical extent of the binary constituents (e.g., the horizon scale in the case of supermassive black hole binaries), we require an effective boundary/sink prescription. While the type of sink can impact the formation of mini-disks and their dynamics (see, e.g., Dittmann &amp; Ryan (2021) for a detailed discussion), we here mainly focus on the evolution of the cavity, and adopt a simple prescription removing matter at the local Keplerian rate, s p = -Ω K,i ρ ( v -v i ), s ρ = -Ω K,i ρ , within a sink radius of r sink = 0 07 . a , where Ω K,i = ( GM i ) 1 / 2 / ( | r -r i | 2 + ϵ 2 ) 3 / 2 is the Keplerian frequency of each sink, ϵ = 0 05 . a is a small number acting as Plummer gravitational softening, and v i the velocity of the sink on a fixed Keplerian orbit. The form of s E follows from the density and momentum sink using the evolution equations. We further model the system as being locally isothermal, with the temperature being fixed to T = c 2 s = Φ | | / M 2 , where c 2 s is the speed of sound and M = 10 is the sonic Mach number (or having disk aspect ratio h = 1 / M = 0 1). . The surface density profile is initialized following Duffell et al. 2024; Mu˜oz &amp; Lithwick 2020, with the vertical n structure being given by assuming constant temperature along the vertical coordinate. We initialize the magnetic field as purely poloidal by adopting a vector potential A ϕ = R A 2 0 max( ρ -0 04 . ρ max , 0) 2 , where R is the cylindrical radius, ρ max the maximum mass density inside the initial disk, and A 0 is chosen to ensure an initial disk magnetization parameter β = 2 P/B 2 = 20. We further adopt a numerical floor value of β floor = 10 -3 .

We numerically integrate the above system using the public GPU-enabled version, AthenaK 1 , of the Athena++ code (Stone et al. 2020) based on the Kokkos library (Trott et al. 2021). In this work, we solve the Newtonian MHD system using piecewise parabolic reconstruction (Colella &amp; Woodward 1984), an HLLD Riemann solver (Miyoshi &amp; Kusano 2005), and a constraint transport algorithm (Gardiner &amp; Stone 2008) for the divergence-free magnetic field evolution. Timestepping is done using a second-order accurate Runge-Kutta stepper.

To perform the simulation, we adopt a Cartesian grid with 6 levels of static mesh refinement. The levels decrease in extent by a factor two each from the outer (outflow) domain boundary inward, which is located at ± 80 a . We adopt a finest resolution of ∆ x ≃ 0 0065 . a . We then numerically integrate the system on 1,260 V100 GPUs on the DOE OLCF Summit system for 230 binary orbits. Overall, the cost of the simulation presented here is about 135,000 V100 GPU-hours.

## 3. RESULTS

We present our results on magnetically arrested circumbinary disk accretion in the following way. First, we discuss the initial build-up of magnetization inside the cavity and its transition to a magnetically arrested state. We then describe the accretion cycle by means of magnetic flux eruptions, before concluding with an in-depth discussion on angular momentum transport.

## 3.1. Evolution towards a magnetically arrested state

We begin by detailing the overall evolution of the structures inside the cavity, which is shown in Fig. 1. Initially, the cavity is empty and begins to fill rapidly through tidal accretion streams, as can be seen after

t

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00757v2-figure-0222c835-6332-4355-a9b4-6a6a6e5db3e5.jpg)

= 8 0

.

t

t

= 51 0

.

t

t

= 110 0

.

t

Figure 1. Initial circumbinary evolution towards a magnetically arrested state over the first 110 binary orbits. (Top and bottom) Meridional slices along the orbital axis, (center) equatorial slices. Shown in color are the mass density, ρ , and magnetization parameter β . Times, t , are stated relative to the orbital period, t orbit , whereas distances are stated relative to the binary separation, a .

Figure 2. Circumbinary accretion cycle in the magnetically arrested state. (Left) Quiescent state, strong vertical magnetic fields (black streamlines) push out the accretion flow and support it. The cavity truncates well inside the inner Lindblad radius. We can also see a previous eruption epsiode to the right of the central cavity. (Center) Interchange instabilities trigger accretion streams through the cavity onto the mass sinks, with the cavity being highly distorted and expanded. We can also see that a reconnection-driven magnetic flare is launched vertically outwards. (Right) This leads to an ejection of vertical net flux from the cavity in the circumbinary disk. The ejected (magnetic pressure dominated) flux tube gets strongly sheared. Shown in color are the mass density, ρ , and magnetization parameter, β , with the top and bottom rows showing the meridional and equatorial planes, respectively. Times, t , are stated relative to the orbital time period, t orb .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00757v2-figure-c6c07780-c378-4d81-b9e6-d87227c8835c.jpg)

eight orbital periods. This phase is not unlike the early onset of two-dimensional circumbinary accretion (e.g., Mu˜ noz et al. 2019; Duffell et al. 2024), or recent threedimensional simulations with an outer accretion flow (Moody et al. 2019; Avara et al. 2023; Pierens &amp; Nelson 2023; Matsumoto 2024). In this stage, rapid minidisk formation sets in, with the outer disk subsequently becoming unstable to the magnetorotational instability (MRI) (Balbus &amp; Hawley 1991) and turbulent. The magnetization, β = 2 P/B 2 ≫ 1, in most parts of the cavity including inside the mini-disks. Overall, despite the more complex effective spatially varying viscosity (see Matsumoto (2024) for a detailed discussion), this stage of our evolution can be characterized as mainly hydrodynamical.

As the evolution proceeds, accretion from the circumbinary disk advects substantial amounts of magnetic flux into the cavity, leading to a build-up of net magnetic flux. We can see that already after 51 orbits, some parts of the cavity are entirely magnetically dominated, with β ≪ 1. While the tidal accretion streams are still predominantly present, the mini-disks start to exhibit turbulent features, with their magnetization level approaching equipartition values, β ≃ 1. The outer disk has become strongly turbulent at the inner cavity wall, with strong magnetically driven outflows present above the sink cells. These jet-like features are consistent with strong net vertical flux inside the cavity. In this transitional state, the flow is still largely hydrodynamical.

After about 100 orbits, however, the flow structure inside the cavity has drastically changed. The continued inflow of magnetic flux into the cavity has lead to a regime that is largely magnetically dominated, with β ≃ 10 -3 almost everywhere inside the cavity. Rather than having a turbulent accretion flow at the inner boundary wall, we can see a very clear and sharp delineation between the cavity and the accretion flow in the meridional plane. This state is reminiscent of a magnetically arrested accretion flow onto a single black hole (Tchekhovskoy et al. 2011; Begelman et al. 2022; Chatterjee &amp; Narayan 2022). There, the magnetic field pushes out the disk, with mass accretion in the quiescent state being confined to a thin current sheet in the equatorial plane of the black holes (Ripperda et al. 2022). Different from hydrodynamical simulations in which the radial location of the cavity is fully determined by the viscous and gravitational torque balance, here the joint torque from the magnetic field and gravitational interaction balances the torque from MRI turbulence in the circumbinary disk, causing a smaller cavity size. We ob-

Figure 3. Mass, ˙ M , electromagnetic energy, ˙ E EM , and angular momentum, ˙ , J fluxes into the cavity, normalized to the mass accretion rate, ˙ M out , at large distances. The time series is shown relative to the orbital period, t orbit ≡ 2 π/ Ω . b The red dashed line denotes the expected value of 〈 ˙ J/M ˙ 〉 ≃ 0 68 . a 2 Ω b for a converged non-eccentric, equal mass, vertically integrated hydrodynamic accretion flow (Mu˜oz et al. n 2019), the cyan dashed line the threshold value, 〈 ˙ J/M ˙ 〉 ≃ 3 / 8 a 2 Ω , for the binary orbit to shrink (Miranda et al. 2017; b Lai &amp; Mu˜ noz 2023).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00757v2-figure-0cc112c4-3fd6-4711-907c-f79bb5e83235.jpg)

serve the presence of a similar equatorial current sheet in our simulations (Fig. 2, top row).

## 3.2. Magnetic flux eruption cycle

One characteristic feature of magnetically arrested accretion flows around single black holes is an intermittent phenomenon regarding the exchange of mass and angular momentum (Tchekhovskoy et al. 2011; Ripperda et al. 2022; Chatterjee &amp; Narayan 2022). The system features two accretion states: a quiescent state and an eruption state. We show this cycle for our circumbinary disk simulation in Fig. 2.

In the quiescent state, strong magnetic fields prevent the cavity from accreting mass except through tidal accretion streams. However, our simulations indicate that the mass and angular momentum carried by these accretion streams might be insufficient to provide for permanent mini-disk formation. We caution that the size of the mini-disks (if any during this accretion phase) that we find in our simulation is subject to resolution effects and magnetic flux accumulation inside the sink, which likely will not influence the cavity wall dynamics we investigate here.

Figure 4. Radial profile of the mass accretion rate, ˙ M , and the normalized magnetic flux, ϕ B ( r ) = Φ B ( r / ) √ ˙ Mr 2 , where Φ ( ) is B r the magnetic flux at cylindrical radius R , given relative to the binary separation, a . Ω orb and M 0 are the binary orbital frequency and initial disk mass, respectively. Different times are shown with different colors, as denoted in the legend.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00757v2-figure-4f22bda6-9311-4ce6-b4dd-eaf5f1b06597.jpg)

At the end of the quiescent state, mass accumulates at the inner cavity wall on a viscous timescale of the outer accretion flow. Eventually the growth of an interchange instability (Spruit et al. 1995) leads to the inward accretion of Rayleigh-Taylor fingers into the cavity, which does not coincide with the tidal accretion streams. This appearance of Rayleigh-Taylor fingers is not unlike accretion states that have been found for T Tauri stars (Kulkarni &amp; Romanova 2008; Takasao et al. 2022), low mass X-ray binaries (Parfrey &amp; Tchekhovskoy 2023; Murguia-Berthier et al. 2024), and magnetically arrested black holes (Begelman et al. 2022; Zhdankin et al. 2023). These magnetic flux eruptions lead to outward-moving vertical magnetic flux tubes that get ejected into the accretion disk (Fig. 2, left and right column). In analogy with the accretion cycle for magnetically arrested black hole flows (Chatterjee &amp; Narayan 2022; Begelman et al. 2022), these low-density, high magnetic pressure flux eruption regions are asymmetric and will lead to net torques on the system, as we will detail in Sec. 3.3. In this state, mass accretion is enhanced, leading to a transient formation of a mini-disk as magnetic pressure

∣

∣

∣

∣

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00757v2-figure-c386afc7-dfa9-4b31-9581-3267f59284a2.jpg)

∣

∣

∣

∣

∣

∣

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00757v2-figure-2a44704b-06fa-47af-b98f-7798c4d6cff8.jpg)

∣

∣

Figure 5. Angular momentum transport in the circumbinary disk during a flux eruption. Shown is the ratio of the magnetic, ˙ J mag , to advected, ˙ J adv , angular momentum fluxes. (Top) Flux eruption in the equatorial plane. Angular momentum is largely transported outwards due to magnetic stresses (pink regions). (Bottom) Angular momentum fluxes inside a flux tube (right). As the flux tube gets ejected into the disk (1 &lt; x/a &lt; 1 5), it drives a wind carrying angular . momentum outwards.

support on the outer cavity wall is released. Concurrently, the rapid accretion of magnetic flux leads to reconnection events that can trigger the ejection of large plasmoids/flares (see Fig. 2, top center panel). These have also been found in magnetically arrested black hole disks, and have been suggested as an explanation for X-ray flares in the galactic center (Porth et al. 2021; Ripperda et al. 2020; Nathanail et al. 2022).

∣

∣

∣

∣ 〈

∣

∣

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00757v2-figure-40eaa4a1-59e4-4e7b-a63c-132a14c65ff7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00757v2-figure-43856a93-679b-4348-a8af-ad56d85ddd0c.jpg)

∣

∣

∣

〉

〈

〉

∣ 〈

∣

∣

Figure 6. Azimuthally and vertically integrated radial angular momentum fluxes as a spacetime diagram. (Left) Total radial angular momentum flux, 〈 ˙ J 〉 z,ϕ . s = +1(outwards) / -1 (inwards) indicates the sign of the flux. (Right) Ratio of the magnetic, 〈 ˙ J mag 〉 z,ϕ , to the advective, 〈 ˙ J adv 〉 z,ϕ , contribution of the radial angular momentum flux. The cylindrical radius, R , is shown in units of the binary separation, a . Times, t , are shown relative to the orbital period, t orbit .

After the magnetic flux tubes got ejected into the disk and sufficient magnetic flux has been re-accreted into the cavity, the cavity restores its magnetically dominated state, resetting the cycle and suppressing RayleighTaylor instabilities, at least intermittently. Meanwhile, the ejected flux tubes get sheared and mixed into the disk (Fig. 2, right column).

We can also quantitatively study the dynamics of the cavity. To this end, we compute the mass accretion rate, ˙ M , into the cavity, which we approximate by a sphere around the binary's center of mass with radius r ≃ 0 6 . a . The cavity we observe is generally smaller than in hydrodynamic simulations (Mu˜oz et al. 2019), although n we caution that our accretion state and cavity size may not yet be fully converged. As shown in Fig. 3, the magnetically arrested flow leads to periodic mass accretion patterns, which vary by a factor of 3-4, when normalized to the mass accretion rate, ˙ M out , at larger distances ( r ≃ 3 5 . a ). Indeed, we can spot clear accretion cycles matching the flux eruptions with a period of 3 -5 orbits as seen in Fig. 2. At the same time, we find that the mass accretion rate is already well converged out to radii of r ≃ 15 a (Fig. 4). This likely implies that the flow we observe is roughly in inflow equilibrium.

In the black hole accretion community, commonly the dimensionless flux ϕ B ( r ) = Φ B ( r / ) √ | ˙ M r | 2 , where Φ B ( r ) is the magnetic flux computed over a a sphere of radius r , commonly the horizon (Tchekhovskoy et al. 2011). In a magnetically arrested state, a large value of ϕ B ≫ 1 implies that energy in excess of the accreted mass energy can be extracted from the black hole (Blandford &amp; Znajek 1977; Tchekhovskoy et al. 2011). We therefore compare ϕ B as a function of radius, r , for different times (Fig. 4). We can see that ϕ B ≫ 10 locally inside the cavity. At the same time the outer disk region is at ϕ B ≲ 1. We also compare the net flux of electromagnetic energy ˙ E EM into the cavity (Fig. 3). Correlating this with the mass accretion rate, we find that spikes in the electromagnetic energy flux usually precede eruption events, consistent with the necessary requirement of replenishing the magnetic flux inside the cavity after a successful eruption.

## 3.3. Angular momentum transport via flux eruptions

One important aspect of circumbinary disk accretion is the question of angular momentum transport, or conversely, whether interactions with the gaseous disk make the binary orbit spin up or shrink for equal mass systems (Mu˜ noz et al. 2019; Duffell et al. 2020, 2024).

Magnetically arrested accretion flows can feature yet another angular momentum transport channel. As has been shown by Chatterjee &amp; Narayan (2022), flux eruptions themselves are able to transport angular momentum outwards. This process is intrinsically intermittent and quasi-periodic, mainly influenced by the magnetic accretion cycle. As we will show in the following, this

picture seems to largely carry over to magnetically arrested accretion flows. Conversely, the periodicity of angular momentum transport presented in hydrodynamical simulations is governed by tidal accretion streams and a slowly precessing eccentric cavity.

We quantify the amount of angular momentum transport by considering the transport equation for angular momentum in cylindrical coordinates, which can be obtained from the Euler equations (see the discussion of vertically integrated angular momentum transfer in Miranda et al. (2017), and the three-dimensional version in Sheikhnezami &amp; Fendt (2022)). The conserved angular momentum density, j = R ρv 2 ϕ , is then given by,

$$\text{mmanual} \, \text{conv} _ { y, \, J } - \imath \, \rho _ { v }, \, \text{as} \, \text{in} \, \text{given} \, \text{v} _ { y }, \quad \text{mu} \\ \partial _ { t } j + \partial _ { R } \left ( j v _ { R } - R ^ { 2 } B _ { R } B _ { \phi } \right ) + \partial _ { \phi } \left ( j - R B _ { \phi } B _ { \phi } + R \bar { P } \right ) \\ + \partial _ { z } \left ( j v _ { z } - R ^ { 2 } B _ { z } B _ { \phi } \right ) + R \rho \partial _ { \phi } \Phi = 0 \,, \quad ( 8 ) \quad \text{poi} \\ \text{where} \, R \, \text{as} \, the cylindrical radius, \, z \, \text{as} \, the \, \text{height} \, \text{coor} \quad \text{flux} \\ \frac { \cdot \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \dots \cdot \colon \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \, }$$

where R as the cylindrical radius, z as the height coordinate, and ϕ the azimuthal coordinate.

We compute the effective angular momentum fluxes as integrals on cylindrical surfaces located at the binary's center of mass, i.e.,

$$\text{sources}, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi,\, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \, \Phi, \,                                                                                                                                                                                                         \partial _ { t } \left \langle j \right \rangle _ { z, \phi } + \partial _ { R } \left ( \left \langle j \right \rangle _ { \text{adv} } + j \right \rangle _ { \text{mag} } + j \right \rangle _ { \text{grav} } \right \rangle _ { z, \phi } \right ) = 0 \,, \, \ ( 9 ) \, \quad \text{acc} \, \quad \text{lat} \, \quad \text{var} \, \quad \text{in} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \, \quad \text{in} \, \quad \text{var} \\ \text{where} \, \ \text{do} \, \ \text{do} \, \ \text{for} \, \ \text{do} \, \ \text{for} \, \ \text{do} \, \ \text{for} \, \ \text{do} \, \ \text{for} \, \ \text{for} \, \ \text{do} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{in} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{if} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{for} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{ try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{catch} \ \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \,\ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \,  \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \ \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try}\, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try} \, \ \text{try}$$

where we have defined the ⟨ . . . ⟩ z,ϕ as the z -and ϕ -integral over the cylinder, and

$$\dot { J } _ { \text{adv} } = R ^ { 2 } \rho v _ { \phi } v _ { R } \,,$$

$$j _ { \text{mag} } = - R ^ { 2 } B _ { R } B _ { \phi } \,, \quad \quad ( 1 1 ) \quad \$$

$$J _ { \text{mag} } & = - \mathcal { H } ^ { - } D _ { R } D _ { \phi } \,, & ( 1 1 ) & \quad \, \\ j _ { \text{grav} } & = \int _ { 0 } ^ { R } d r ^ { \prime } r ^ { \prime } \rho \, \partial _ { \phi } \Phi \,, & ( 1 2 ) & \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \$$

as the advective and magnetic angular momentum flux, respectively. We further define the total averaged angular momentum flux as ˙ J = ˙ J adv + ˙ J mag + ˙ J grav . Following Chatterjee &amp; Narayan (2022), we define the sign of the angular momentum flux such that ˙ J &gt; 0 corresponds to an outward flux of angular momentum, whereas ˙ J &lt; 0 corresponds to an inward one.

Using this decomposition, we can analyze the torque carried by a flux eruption event (Fig. 5). Focusing on the same accretion cycle shown in Fig. 2, we can see that the torque inside the flux eruption is entirely magnetically dominated. In magnetically dominated regions of the erupted flux tube, magnetic angular momentum fluxes are up to 100-times larger than the advected hydrodynamical fluxes. Outside of the eruption region, and also inside the Rayleigh-Taylor fingers, the angular momentum fluxes are largely hydrodynamical. The ejected flux tube carries angular momentum into the disk, which will eventually drive a wind/outflow due to the strong vertical net flux. This behavior is consistent with flux eruptions around single black holes (Chatterjee &amp; Narayan 2022; Lowell et al. 2024).

We can further quantify the flux transport behavior throughout the evolution of the system over several eruption cycles. We show the resulting temporal evolution of the angular momentum flux in Fig. 6, focusing on the evolution after the first 100 orbits when the system has become magnetically arrested. We can see that consistent with the appearance of magnetically arrested accretion cycles (over several binary orbits) the system features periodic outbursts of radial angular momentum flux. This flux starts roughly at or slightly outside the stationary cavity radius, r ≃ a , and extends outward to several binary distances. These fluxes are also outwardpointing. Indeed, when comparing the strength of the magnetically driven to the advected angular momentum flux (Fig. 6, right panel), we find that these outbursts are entirely magnetically driven, consistent with the flux eruption picture (Fig. 5), and with a similar analysis in magnetically arrested black hole disks (Chatterjee &amp; Narayan 2022).

It can be shown that for an equal mass non-eccentric accreting binary, such as the one consider here, the evolution of the orbital separation, a , should follow (Miranda et al. 2017; Lai &amp; Mu˜oz 2023), n

$$\cdots \, \infty, \, \infty \, \max \, \infty \, \max \, \infty _ { J }, \\ \frac { \dot { a } } { a } = 8 \left ( \ell _ { 0 } - \frac { 3 } { 8 } \right ) \frac { \langle \dot { M } \rangle } { M } \,, \quad \quad \quad ( 1 3 ) \\ \cdots \, \infty \, \cdots \, \infty \, \cdots \, \infty \, \cdots \,.$$

where ℓ 0 ≡ ⟨ ˙ J / ⟩ ( ⟨ ˙ M a ⟩ 2 Ω b ) , and Ω b is the Keplerian frequency of the binary. From this equation, we can see that for ℓ 0 &lt; 3 / 8, the binary separation would shrink. Hydrodynamical simulations in two-dimensions imply ℓ 0 ≃ 0 68 . (Mu˜ noz et al. 2019), whereas threedimensional simulations find ℓ 0 = 0 55 (Moody et al. . 2019). Both these values lead to an eventual growth of the binary orbit.

On the other hand, we find that during flux eruptions ℓ 0 can instantaneously drop below the threshold value, ℓ 0 &lt; 3 / 8, or even become negative (Fig. 3, bottom panel). Correlating these patterns with those of the mass accretion rate, ˙ M , we find that this behavior correlates directly with the magnetically arrested accretion cycle. In the quiescent phase, accretion of angular momentum onto the binary appears to be hydrodynamical and approaches approximately the equilibrium value, ℓ 0 ≃ 0 68, . found in two-dimensional studies. During flux eruptions when transport is entirely magnetically driven, angular momentum can move outwards.

We caution, that while our findings qualitatively match results for single black hole accretion flows (Chatterjee &amp; Narayan 2022; Lowell et al. 2024), our overall flow is

very likely not yet converged, and the cavity could still develop strong eccentricities (Mu˜oz et al. 2019). n A determination of which effect ultimately dominates total angular momentum transfer and budget in this system, will likely require substantially longer simulations than the ones presented here.

## 4. CONCLUSIONS

In this work, we have numerically investigated a magnetically arrested accretion state of a circumbinary disk.

Our simulations indicate, that just as for single accreting black holes, circumbinary accretion flows can become magnetically arrested. In this state, the inner cavity is filled with a strong magnetic flux regulating mass and angular momentum transport onto the binary (Chatterjee &amp; Narayan 2022). Magnetically arrested accretion states have been well studied in the black accretion literatures (Narayan et al. 2003; Igumenshchev et al. 2003; Igumenshchev 2008; Tchekhovskoy et al. 2011). Different from these works, we find that the presence of tidal forces from the binary induces constantly present tidal streams, which, however, are periodically disrupted and reformed during the magnetically arrested accretion cycle.

Periodic magnetic flux eruptions through magnetic interchange instabilities (Spruit et al. 1995; Kulkarni &amp; Romanova 2008), will then enhance mass transport inwards, and transiently drive angular momentum outwards. At the same time, this process triggers reconnection and the ejection of plasmoids and flares, akin to models for X-ray flares in the galactic center (Porth et al. 2021; Scepi et al. 2022; Ripperda et al. 2020, 2022). The ejection of flux tubes into the disk is likely also accompanied with potential non-thermal emission (Zhdankin et al. 2023; Ripperda et al. 2022). The mixing of the ejected flux tubes with strong vertical net flux with the disk may also affect wind launching from the system (Blandford &amp; Payne 1982; Zhu &amp; Stone 2018; Zhu et al. 2024). Furthermore, the strong magnetic field present in disk alters the cavity truncation radius, making it smaller than what is commonly observed in hydrodynamical flows (Mu˜oz et al. 2019). This may have n implications during the decoupling phase of the binary, potentially avoiding a drop of the mass accretion rate (Noble et al. 2021; Avara et al. 2023).

Most importantly for our context, magnetically arrested accretion flows around black holes have been shown to transport angular momentum outwards in magnetic flux eruption episodes (Chatterjee &amp; Narayan 2022), although for thick magnetically arrested disks the main spin-down mechanism may originate from the jet

(Lowell et al. 2024). We similarly find that magnetic flux eruption episodes carry angular momentum outwards, which can satisfy the criterion for the binary orbit to shrink. In this sense, magnetically arrested flows may offer a way to provide the necessary angular momentum loss for the binary's orbit to shrink. We, indeed, find prolonged binary cycles over which the orbit would be shrinking (see Fig. 3). However, we caution that our results are very likely not yet in a fully converged regime, with longer integration times being necessary to render a verdict on the actual net angular momentum transport rate.

While we find the period of the magnetically arrested circumbinary accretion cycle to be between 3 -5 orbits, the precise duration will likely depend on the amount of magnetic flux inside the accretion flow. Here we model the initial accretion disk as having largely poloidal flux, enabling the disk to easily replenish the amount of magnetic flux inside the cavity. Recent works on large-scale feeding onto supermassive black holes indicate that the magnetic field at large distances may be toroidal, with the disk being largely magnetically and radiation dominated (e.g., Hopkins et al. 2024; Hopkins 2024; see also Cho et al. 2023; Guo et al. 2024). When accreting such magnetic fluxes, further dynamo action either inside the circumbinary or mini-disks may likely be necessary to produce the required net vertical fluxes for the cavity to make the flux eruption happen (e.g., Liska et al. 2020; Jacquemin-Ide et al. 2024; Sadanari et al. 2024). This will potentially prolong the timescale between magnetically arrested cycles, likely making this feature intermittent.

The cavity structure and mini-disk formation could also be affected by changes in the magnetic field topology or large-scale fluxes. Additionally, we point out that our sink size is only 7% of the binary separation, and does not include an effective resistivity to regulate magnetic flux accumulation in the sink. Further investigations will be required as to how this will affect the magnetic flux filling on the cavity. Thermal effects in such disks may also naturally alter the circumbinary accretion state (Wang et al. 2023b). We plan to study the role of large scale magnetic fluxes in a future work.

Besides uncertainties in the (effectively initial) magnetic flux distribution, the development of an eccentric cavity (Mu˜ noz et al. 2019) may likely alter the onset of the interchange instability and the periodicity of the flux eruption. Quantifying this effect will require simulations such as the one presented here, over likely 10-times longer timescales, which is beyond the scope of this initial work, but well within in reach of modern GPU-based simulation approaches.

## ACKNOWLEDGMENTS

The authors are grateful for insightful discussions with Xue-Ning Bai, Manuela Campanelli, Luciano Combi, Alexander Dittmann, Philip F. Hopkins, Julian H. Krolik, Douglas N. C. Lin, Scott C. Noble, E. Sterl Phinney, Bart Ripperda, James M. Stone, Alexander Tchekhovskoy, and Zhaohuan Zhu. The simulations were performed on DOE OLCF Summit under allocation AST198. Additional simulations were done on DOE NERSC supercomputer Perlmutter under grant m4575. HYW thanks for the support from B. Thomas Soifer

Chair's Graduate Fellowship. Part of this project was completed during the second LISA Sprint meeting at Caltech. The meeting was supported by the Jet Propulsion Laboratory Astronomy and Physics Directorate. This research was supported in part by grant NSF PHY2309135 to the Kavli Institute for Theoretical Physics (KITP). This work was performed in part at Aspen Center for Physics, which is supported by National Science Foundation grant PHY-2210452.

Software: AthenaK (Stone et al. 2020), Kokkos (Trott et al. 2021), matplotlib (Hunter 2007), numpy (Harris et al. 2020), scipy (Virtanen et al. 2020)

## REFERENCES

| Agazie, G., Anumarlapudi, A., Archibald, A. M., et al. 2023a, ApJL, 951, L8, doi: 10.3847/2041-8213/acdac6 -. 2023b, ApJL, 952, L37, doi: 10.3847/2041-8213/ace18b                                                                |
|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Akiyama, K., et al. 2019, Astrophys. J. Lett., 875, L5, doi: 10.3847/2041-8213/ab0f43 -. 2021, Astrophys. J. Lett., 910, L13, doi: 10.3847/2041-8213/abe4de -. 2022, Astrophys. J. Lett., 930, L16, doi: 10.3847/2041-8213/ac6672 |
| Amaro-Seoane, P., Andrews, J., Arca Sedda, M., et al. 2023, Living Reviews in Relativity, 26, 2, doi: 10.1007/s41114-022-00041-y                                                                                                  |
| Astrophys. J., 913, 16, doi: 10.3847/1538-4357/abf0af                                                                                                                                                                             |
| Arons, J., & Lea, S. M. 1976, ApJ, 207, 914, doi: 10.1086/154562                                                                                                                                                                  |
| Artymowicz, P., & Lubow, S. H. 1994, ApJ, 421, 651, doi: 10.1086/173679 -. 1996, ApJL, 467, L77, doi: 10.1086/310200                                                                                                              |
| Avara, M. J., Krolik, J. H., Campanelli, M., et al. 2023. https://arxiv.org/abs/2305.18538                                                                                                                                        |
| Balbus, S. A., & Hawley, J. F. 1991, Astrophys. J., 376, 214, doi: 10.1086/170270                                                                                                                                                 |
| Begelman, M. C., Blandford, R. D., & Rees, M. J. 1980, Nature, 287, 307, doi: 10.1038/287307a0                                                                                                                                    |
| Begelman, M. C., Scepi, N., & Dexter, J. 2022, Mon. Not. Roy. Astron. Soc., 511, 2040, doi: 10.1093/mnras/stab3790                                                                                                                |
| Blandford, R. D., & Payne, D. G. 1982, Mon. Not. Roy. Astron. Soc., 199, 883, doi: 10.1093/mnras/199.4.883 Blandford, R. D., & Znajek, R. L. 1977, Mon. Not. Roy. Astron. Soc., 179, 433, doi: 10.1093/mnras/179.3.433            |
| Bogdanovi´, c T., Miller, M. C., & Blecha, L. 2022, Living Reviews in Relativity, 25, 3, doi: 10.1007/s41114-022-00037-8                                                                                                          |

Bourne, M. A., Fiacconi, D., Sijacki, D., Piotrowska, J. M., &amp; Koudmani, S. 2023, arXiv e-prints, arXiv:2311.17144, doi: 10.48550/arXiv.2311.17144

Bowen, D. B., Campanelli, M., Krolik, J. H., Mewes, V., &amp; Noble, S. C. 2017, Astrophys. J., 838, 42, doi: 10.3847/1538-4357/aa63f3

Bowen, D. B., Mewes, V., Campanelli, M., et al. 2018,

Astrophys. J. Lett., 853, L17, doi: 10.3847/2041-8213/aaa756

Bright, J. C., &amp; Paschalidis, V. 2023, Mon. Not. Roy.

Astron. Soc., 520, 392, doi: 10.1093/mnras/stad091

Charisi, M., Taylor, S. R., Runnoe, J., Bogdanovic, T., &amp;

Trump, J. R. 2022, MNRAS, 510, 5929, doi: 10.1093/mnras/stab3713

Chatterjee, K., &amp; Narayan, R. 2022, ApJ, 941, 30, doi: 10.3847/1538-4357/ac9d97

Cho, H., Prather, B. S., Narayan, R., et al. 2023, Astrophys.

J. Lett., 959, L22, doi: 10.3847/2041-8213/ad1048

Colella, P., &amp; Woodward, P. R. 1984, Journal of

Computational Physics, 54, 174, doi: 10.1016/0021-9991(84)90143-8

Collaboration, T. E. H. T., et al. 2024, Astrophys. J. Lett.,

964, L26, doi: 10.3847/2041-8213/ad2df1

Combi, L., Armengol, F. G. L., Campanelli, M., et al. 2022, Astrophys. J., 928, 187, doi: 10.3847/1538-4357/ac532a

Davis, S. W., &amp; Tchekhovskoy, A. 2020, Ann. Rev. Astron. Astrophys., 58, 407, doi: 10.1146/annurev-astro-081817-051905

DeLaurentiis, S., Haiman, Z., Westernacher-Schneider, J. R., et al. 2024, arXiv e-prints, arXiv:2405.07897, doi: 10.48550/arXiv.2405.07897

Derdzinski, A., D'Orazio, D., Duffell, P., Haiman, Z., &amp; MacFadyen, A. 2021, Mon. Not. Roy. Astron. Soc., 501, 3540, doi: 10.1093/mnras/staa3976

| Dexter, J., et al. 2020, Mon. Not. Roy. Astron. Soc., 497, 4999, doi: 10.1093/mnras/staa2288                                                                         |
|----------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Dittmann, A. J., & Ryan, G. 2021, Astrophys. J., 921, 71, doi: 10.3847/1538-4357/ac1bbd                                                                              |
| Dittmann, A. J., & Ryan, G. 2024, ApJ, 967, 12, doi: 10.3847/1538-4357/ad2f1e                                                                                        |
| Dittmann, A. J., Ryan, G., & Miller, M. C. 2023, ApJL, 949, L30, doi: 10.3847/2041-8213/acd183                                                                       |
| D'Orazio, D. J., & Duffell, P. C. 2021, ApJL, 914, L21, doi: 10.3847/2041-8213/ac0621                                                                                |
| Dotti, M., Colpi, M., Haardt, F., & Mayer, L. 2007, MNRAS, 379, 956, doi: 10.1111/j.1365-2966.2007.12010.x                                                           |
| Duffell, P. C., D'Orazio, D., Derdzinski, A., et al. 2020, ApJ, 901, 25, doi: 10.3847/1538-4357/abab95                                                               |
| Duffell, P. C., D'Orazio, D., Derdzinski, A., et al. 2020, Astrophys. J., 901, 25, doi: 10.3847/1538-4357/abab95                                                     |
| Duffell, P. C., Dittmann, A. J., D'Orazio, D. J., et al. 2024, arXiv e-prints, arXiv:2402.13039, doi: 10.48550/arXiv.2402.13039                                      |
| Duffell, P. C., et al. 2024. https://arxiv.org/abs/2402.13039                                                                                                        |
| Farris, B. D., Gold, R., Paschalidis, V., Etienne, Z. B., & Shapiro, S. L. 2012, Phys. Rev. Lett., 109, 221102, doi: 10.1103/PhysRevLett.109.221102                  |
| Franchini, A., Sesana, A., & Dotti, M. 2021, MNRAS, 507, 1458, doi: 10.1093/mnras/stab2234                                                                           |
| Gardiner, T. A., & Stone, J. M. 2008, Journal of Computational Physics, 227, 4123,                                                                                   |
| Gold, R., Paschalidis, V., Ruiz, M., et al. 2014, Phys. Rev. D, 90, 104030, doi: 10.1103/PhysRevD.90.104030                                                          |
| Goldreich, P. M., & Tremaine, S. 1979, Astrophys. J., 233, 857, doi: 10.1086/157448                                                                                  |
| GRAVITY Collaboration, Abuter, R., Amorim, A., et al. 2018, A&A, 618, L10, doi: 10.1051/0004-6361/201834294 -. 2021, A&A, 654, A22, doi: 10.1051/0004-6361/202140981 |
| Guo, M., Stone, J. M., Quataert, E., & Kim, C.-G. 2024, arXiv e-prints, arXiv:2405.11711, doi: 10.48550/arXiv.2405.11711                                             |
| Guti´rrez, e E. M., Combi, L., Noble, S. C., et al. 2022, Astrophys. J., 928, 137, doi: 10.3847/1538-4357/ac56de                                                     |
| Hakobyan, H., Ripperda, B., & Philippov, A. 2023, Astrophys. J. Lett., 943, L29, doi: 10.3847/2041-8213/acb264                                                       |
| Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2                                                    |
| Hopkins, P. F. 2024, arXiv e-prints, arXiv:2407.00160, doi: 10.48550/arXiv.2407.00160                                                                                |

| Hopkins, P. F., Squire, J., Su, K.-Y., et al. 2024, The Open Journal of Astrophysics, 7, 19, doi: 10.21105/astro.2310.04506                 |
|---------------------------------------------------------------------------------------------------------------------------------------------|
| Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, doi: 10.1109/MCSE.2007.55                                                    |
| Igumenshchev, I. V. 2008, Astrophys. J., 677, 317, doi: 10.1086/529025                                                                      |
| Igumenshchev, I. V., Narayan, R., & Abramowicz, M. A. 2003, Astrophys. J., 592, 1042, doi: 10.1086/375769                                   |
| Jacquemin-Ide, J., Rincon, F., Tchekhovskoy, A., & Liska, M. 2024, Mon. Not. Roy. Astron. Soc., 532, 1522, doi: 10.1093/mnras/stae1538      |
| Kratter, K. M., Matzner, C. D., & Krumholz, M. R. 2008, ApJ, 681, 375, doi: 10.1086/587543                                                  |
| Kulkarni, A. K., & Romanova, M. M. 2008, Mon. Not. Roy. Astron. Soc., 386, 673,                                                             |
| Lai, D., & Mu˜oz, n D. J. 2023, ARA&A, 61, 517, doi: 10.1146/annurev-astro-052622-022933                                                    |
| Liska, M. T. P., Tchekhovskoy, A., & Quataert, E. 2020, Mon. Not. Roy. Astron. Soc., 494, 3656,                                             |
| Lowell, B., Jacquemin-Ide, J., Tchekhovskoy, A., & Duncan, A. 2024, Astrophys. J., 960, 82, doi: 10.3847/1538-4357/ad09af                   |
| Matsumoto, T. 2024, ApJ, 964, 133, doi: 10.3847/1538-4357/ad25ee                                                                            |
| Mayer, L., Kazantzidis, S., Madau, P., et al. 2007, Science, 316, 1874, doi: 10.1126/science.1141858                                        |
| Milosavljevi´, c M., & Merritt, D. 2001, ApJ, 563, 34, doi: 10.1086/323830                                                                  |
| Milosavljevi´, c M., & Phinney, E. S. 2005, ApJL, 622, L93, doi: 10.1086/429618                                                             |
| T., & Kusano, K. 2005, Journal of Computational                                                                                             |
| Miyoshi, Physics, 208, 315, doi: 10.1016/j.jcp.2005.02.017                                                                                  |
| Moody, M. S. L., Shi, J.-M., & Stone, J. M. 2019, ApJ, 875, 66, doi: 10.3847/1538-4357/ab09ee                                               |
| Mu˜oz, n D. J., & Lithwick, Y. 2020, Astrophys. J., 905, 106, doi: 10.3847/1538-4357/abc74c                                                 |
| Mu˜oz, n D. J., Miranda, R., & Lai, D. 2019, ApJ, 871, 84, doi: 10.3847/1538-4357/aaf867                                                    |
| Murguia-Berthier, A., Parfrey, K., Tchekhovskoy, A., & Jacquemin-Ide, J. 2024, Astrophys. J. Lett., 961, L20, doi: 10.3847/2041-8213/ad16eb |
| Narayan, R., Igumenshchev, I. V., & Abramowicz, M. A. 2003, PASJ, 55, L69, doi: 10.1093/pasj/55.6.L69                                       |

| Nathanail, A., Mpisketzis, V., Porth, O., Fromm, C. M., & Rezzolla, L. 2022, Mon. Not. Roy. Astron. Soc., 513, 4267, doi: 10.1093/mnras/stac1118   |
|----------------------------------------------------------------------------------------------------------------------------------------------------|
| Noble, S. C., Krolik, J. H., Campanelli, M., et al. 2021, Astrophys. J., 922, 175, doi: 10.3847/1538-4357/ac2229                                   |
| Noble, S. C., Mundim, B. C., Nakano, H., et al. 2012, Astrophys. J., 755, 51, doi: 10.1088/0004-637X/755/1/51                                      |
| Parfrey, K., & Tchekhovskoy, A. 2023. https://arxiv.org/abs/2311.04291                                                                             |
| Paschalidis, V., Bright, J., Ruiz, M., & Gold, R. 2021, Astrophys. J. Lett., 910, L26, doi: 10.3847/2041-8213/abee21                               |
| Pierens, A., & Nelson, R. P. 2023, A&A, 670, A112, doi: 10.1051/0004-6361/202244828                                                                |
| Porth, O., Mizuno, Y., Younsi, Z., & Fromm, C. M. 2021, Mon. Not. Roy. Astron. Soc., 502, 2023, doi: 10.1093/mnras/stab163                         |
| Ripperda, B., Bacchini, F., & Philippov, A. 2020, Astrophys. J., 900, 100, doi: 10.3847/1538-4357/ababab                                           |
| Astrophys. J. Lett., 924, L32, doi: 10.3847/2041-8213/ac46a1 Sadanari, K. E., Omukai, K., Sugimura, K., Matsumoto,                                 |
| T., & Tomida, K. 2024, PASJ, doi: 10.1093/pasj/psae051                                                                                             |
| doi: 10.1086/155823                                                                                                                                |
| Scharlemann, E. T. 1978, ApJ, 219, 617,                                                                                                            |
| Sheikhnezami, S., & Fendt, C. 2022, ApJ, 925, 161, doi: 10.3847/1538-4357/ac3f31                                                                   |
| Shi, J.-M., Krolik, J. H., Lubow, S. H., & Hawley, J. F. 2012, Astrophys. J., 749, 118, doi: 10.1088/0004-637X/749/2/118                           |
| Shi, Y., Kremer, K., & Hopkins, P. F. 2024, ApJL, 969, L31, doi: 10.3847/2041-8213/ad5a95                                                          |
| Siwek, M., Weinberger, R., & Hernquist, L. 2023a, MNRAS, 522, 2707, doi: 10.1093/mnras/stad1131                                                    |
| Siwek, M., Weinberger, R., Mu˜oz, n D. J., & Hernquist, L. 2023b, MNRAS, 518, 5059, doi: 10.1093/mnras/stac3263                                    |
| Spruit, H. C., Stehle, R., & Papaloizou, J. C. B. 1995, Mon. Not. Roy. Astron. Soc., 275, 1223, doi: 10.1093/mnras/275.4.1223                      |

| Stone, J. M., Gardiner, T. A., Teuben, P., Hawley, J. F., & Simon, J. B. 2008, Astrophys. J. Suppl., 178, 137, doi: 10.1086/588755                 |
|----------------------------------------------------------------------------------------------------------------------------------------------------|
| Stone, J. M., Tomida, K., White, C. J., & Felker, K. G. 2020, ApJS, 249, 4, doi: 10.3847/1538-4365/ab929b                                          |
| Sudarshan, P., Penzlin, A. B. T., Ziampras, A., Kley, W., & Nelson, R. P. 2022, arXiv e-prints, arXiv:2206.07749. https://arxiv.org/abs/2206.07749 |
| Takasao, S., Tomida, K., Iwasaki, K., & Suzuki, T. K. 2022, Astrophys. J., 941, 73, doi: 10.3847/1538-4357/ac9eb1                                  |
| Mon. Not. Roy. Astron. Soc., 418, L79, doi: 10.1111/j.1745-3933.2011.01147.x                                                                       |
| Tiede, C., Zrake, J., MacFadyen, A., & Haiman, Z. 2020, ApJ, 900, 43, doi: 10.3847/1538-4357/aba432                                                |
| Trott, C., Berger-Vergiat, L., Poliakoff, D., et al. 2021, Computing in Science and Engineering, 23, 10, doi: 10.1109/MCSE.2021.3098509            |
| Tsokaros, A., Ruiz, M., Shapiro, S. L., & Paschalidis, V. 2022, Phys. Rev. D, 106, 104010, doi: 10.1103/PhysRevD.106.104010                        |
| Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261, doi: 10.1038/s41592-019-0686-2                                   |
| Vos, J., Davelaar, J., Olivares, H., Brinkerink, C., & Falcke, H. 2023. https://arxiv.org/abs/2310.16938                                           |
| Wang, H.-Y., Bai, X.-N., & Lai, D. 2023a, ApJ, 943, 175, doi: 10.3847/1538-4357/acac77                                                             |
| Wang, H.-Y., Bai, X.-N., Lai, D., & Lin, D. N. C. 2023b, MNRAS, 526, 3570, doi: 10.1093/mnras/stad2884                                             |
| Williams, J. P., & Cieza, L. A. 2011, ARA&A, 49, 67, doi: 10.1146/annurev-astro-081710-102548                                                      |
| Astrophys., 52, 529, doi: 10.1146/annurev-astro-082812-141003 Zhdankin, V., Ripperda, B., & Philippov, A. A. 2023,                                 |
| Phys. Rev. Res., 5, 043023, doi: 10.1103/PhysRevResearch.5.043023                                                                                  |
| Zhu, Z., & Stone, J. M. 2018, ApJ, 857, 34, doi: 10.3847/1538-4357/aaafc9                                                                          |
| Zhu, Z., Stone, J. M., & Calvet, N. 2024, MNRAS, 528, 2883, doi: 10.1093/mnras/stad3712                                                            |
Elias R. Most Hai-Yang Wang
Magnetically Arrested Circumbinary Accretion Flows
Binary systems with comparable masses and a surrounding accretion disk can
accrete gas through spiral accretion streams penetrating the central cavity
formed by tidal interactions. Using three-dimensional Newtonian
magnetohydrodynamics simulations, we investigate the possibility of a
magnetically arrested accretion flow through the cavity. Rather than solely
continuously feeding the binary through spiral accretion streams, the accretion
is regulated by the strong magnetic field inside the cavity. Transport of mass
and angular momentum onto the binary then proceeds largely periodically in
magnetic flux eruption episodes. The ejected flux tubes carry angular momentum
outwards and away from the binary, inject hot plasma into the disk and can
launch flares. This likely intermittent scenario could have potential
implications for the emission signatures of supermassive black hole binaries,
and shed light onto the role magnetic fields play in the binary's orbital
evolution.