Flux Accretion Model in Astrophysics

• Flux Accretion Model is a framework that explains the role and dynamics of large-scale poloidal magnetic flux in accretion systems, emphasizing global and topological processes.
• It employs GRMHD and GRPIC simulations to analyze how coronal transport, reconnection, and flux saturation govern angular momentum transfer and state transitions.
• These models impact our understanding of jet formation, luminosity variability, and the cyclic behavior seen in magnetically arrested versus standard accretion states.

“Flux accretion model” denotes a family of astrophysical descriptions in which accretion is controlled, diagnosed, or observationally inferred through a transported flux. In black-hole plasma theory the phrase most commonly refers to the inward transport, saturation, reconnection, and eruption of large-scale poloidal magnetic flux, especially in GRMHD descriptions of weakly magnetized and magnetically arrested flows; in other literatures it can denote heat-flux-driven magnetic seed generation, or a mapping from time-dependent accretion rate to observed electromagnetic flux in compact objects and young stellar objects [(0906.2784); (Chatterjee et al., 2022); (Villarroel-Sepúlveda et al., 2024); (Das et al., 19 May 2026)]. This suggests that the expression is not used in a single standardized sense, but its dominant usage concerns the dynamics of magnetic flux in accretion systems.

1. Magnetic flux as an accretion variable

In the magnetic-transport literature, the relevant quantity is large-scale poloidal magnetic flux. Its generic definition is

$\Phi=\int_S \mathbf{B}\cdot d\mathbf{S},$

and in relativistic accretion studies it is commonly evaluated through the radial magnetic field on surfaces of constant radius or represented through contours of the azimuthal vector potential $A_\phi$, whose differences encode enclosed poloidal flux (0906.2784). The corresponding induction law is written in ideal-MHD form as

$$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$

but the global evolution of flux depends not only on local advection but also on changes of field-line connectivity through reconnection (0906.2784).

A central conclusion of the modern black-hole accretion literature is that mass inflow and magnetic-flux inflow are distinct processes. Matter can accrete through an MRI-turbulent disk while net poloidal flux is transported mainly through coronal motions, reconnection, or episodic eruptions rather than by smooth comoving advection in the dense equatorial flow [(0906.2784); (Chatterjee et al., 2022)]. This distinction is explicit in long-duration GRMHD simulations of advection-dominated accretion onto non-spinning black holes, where the average outward angular momentum transport is primarily radial in both SANE and MAD states, but the vertical transport channel becomes strong in MADs because magnetic-flux-eruption-driven winds carry angular momentum away from the disk (Chatterjee et al., 2022).

The same emphasis on flux as a dynamical variable appears in kinetic models of magnetized spherical accretion. There the horizon-threading flux $\Phi_{\rm H}$ is treated alongside the mass accretion rate $\dot M$, and the basic transport equation for the flux function $\Phi(r,\theta,t)$ is written in ideal form as

$\left(\partial_t+V^r\partial_r\right)\Phi=0,$

showing that flux growth is tied to the radial inflow only in phases where reconnection remains subdominant (Figueiredo et al., 4 Feb 2026).

2. Coronal transport and the rejection of purely local advection–diffusion pictures

A foundational result for magnetic-flux accretion in black-hole disks is the “coronal mechanism” identified in a global three-dimensional GRMHD simulation of a torus embedded in a large-scale vertical magnetic field (0906.2784). In that picture, orbital shear stretches the imposed field into large-scale half-loops in the low-density corona and disk surface layers. These loops move inward largely outside the turbulent disk body, then reconnect near the inner flow or funnel boundary, causing discontinuous jumps in the location of net flux (0906.2784).

The same work reports little or no substantial direct inward transport of net poloidal flux through the disk proper. The dense equatorial flow accretes mass, but the corona is the principal conveyor of organized large-scale flux; repeated coronal stretching and reconnection establish a dipolar magnetic field in the evacuated funnel around the axis (0906.2784). This directly undermines local formulations in which flux transport is described as a competition between an “effective viscosity” that advects field inward and an “effective resistivity” that diffuses it outward. The paper’s argument is that flux is intrinsically global, nonlinear, and topological, so local transport coefficients do not capture the dominant mechanism (0906.2784).

Analytic quasi-local modeling later reformulated the same problem for discs threaded by net vertical magnetic flux. In a slab model, the enclosed vertical flux

$\Phi(r,t)\equiv \int_0^r r\,dr\,B_z$

obeys a transport law in which inward advection of $B_z$ competes with outward turbulent diffusion through the disc (Begelman, 2024). The sensitivity of flux accumulation is controlled by a single parameter

$q\equiv \frac{\sqrt{3}\,v_{{\rm A}z}H}{\eta},$

interpreted as the ratio of vertical diffusion time to Alfvén crossing time (Begelman, 2024). Under wide-ranging conditions, inflow is governed by large-scale magnetic stresses rather than internal viscous stress, leading to “magnetically boosted” discs, and the time-dependent version of the model permits steady or episodic MAD-like central flux concentrations (Begelman, 2024).

3. Magnetically arrested states, eruptions, and cyclic flux release

In the MAD literature, flux accretion is not monotonic. Long-duration GRMHD simulations of sub-Eddington accretion onto non-spinning black holes show that magnetic-flux eruptions in MADs push gas out of the disk and launch strong winds with outflow efficiencies at times reaching $A_\phi$0 of the incoming accretion power, while average mass outflow rates remain small out to $A_\phi$1, reaching only $A_\phi$2 of the horizon accretion rate (Chatterjee et al., 2022). These simulations also find that the MAD state is highly transitory and non-axisymmetric: following an eruption, the accretion mode often changes to a SANE-like state before reattaining magnetic-flux saturation (Chatterjee et al., 2022).

A related GRMHD flare model for Sgr A* interprets violent episodes of flux escape from the black-hole magnetosphere as the source of coherent magnetized flux bundles orbiting in the inner MAD flow (Porth et al., 2020). In that model, accretion advects large-scale poloidal flux inward until a saturated MAD state is reached; excess flux then escapes into the disk as vertically magnetized, MRI-suppressed structures. The tracked bundles circularize at $A_\phi$3, have magnetic energies up to $A_\phi$4 erg, and are dynamically sub-Keplerian, a point explicitly identified as tension with the GRAVITY period–radius relation for three Sgr A* flares (Porth et al., 2020).

A kinetic extension of the flux-eruption picture replaces fluid closure by first-principles GRPIC dynamics in magnetized spherical accretion onto a Schwarzschild black hole (Figueiredo et al., 4 Feb 2026). There each cycle passes through three stages: an ideal advection phase in which $A_\phi$5 increases quasi-linearly with time; a reconnection-regulated phase in which intermittent near-horizon reconnection slows further growth and drives $A_\phi$6 to a saturation value of about $A_\phi$7; and a flaring phase in which large-scale reconnection expels horizon-threading flux, reducing $A_\phi$8 by a factor of $A_\phi$9 over $$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$0 within a cycle of $$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$1 (Figueiredo et al., 4 Feb 2026). The eruption phase obeys an exponential decay law,

$$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$2

with a fitted effective decay rate $$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$3, implying $$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$4 for $$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$5 (Figueiredo et al., 4 Feb 2026).

These eruptive models establish a common picture: flux is accumulated by accretion, partially diffused or redistributed by reconnection, and then expelled catastrophically once current sheets or magnetospheric stresses reach a threshold. The principal difference between fluid and kinetic treatments lies in microphysics: the GRPIC study resolves current-sheet thinning, plasmoid onset, and species-dependent particle acceleration directly, whereas the GRMHD models treat reconnection numerically (Porth et al., 2020, Figueiredo et al., 4 Feb 2026).

4. Flux generation and analytic transport theories

Not all flux-accretion models begin with externally supplied magnetic flux. One line of work treats large-scale flux as generated in situ. The Cosmic Battery model proposes that aberrated radiation pressure on orbiting electrons in the inner luminous flow induces a toroidal electric current and thereby generates poloidal magnetic loops (Contopoulos, 2019). In its flux-transport form, the induction equation is written as

$$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$6

so inward advection, outward diffusion, and a radiation-driven battery term jointly determine the secular flux budget (Contopoulos, 2019). In that picture the inner polarity accumulates near the black hole while the return polarity diffuses outward through the resistive disk during the rise to outburst; after the inner flow approaches equipartition and the Battery ceases, the transport bias reverses and flux reconnects with the accumulated central bundle during decline (Contopoulos, 2019).

A distinct non-ideal generation mechanism replaces radiation drag by conductive heat flux. In a Schwarzschild thin disk with purely radial thermodynamic profiles, the Biermann and relativistic battery terms vanish because $$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$7, $$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$8, and $$\frac{\partial \mathbf{B}}{\partial t}=\nabla\times(\mathbf{v}\times\mathbf{B}),$$9 are parallel, but the curl of the heat-flux source $\Phi_{\rm H}$0 remains nonzero and becomes the only seed-field source (Villarroel-Sepúlveda et al., 2024). The model predicts linearly time-growing magnetic field generation from an initially unmagnetized state and identifies radii $\Phi_{\rm H}$1 and $\Phi_{\rm H}$2 where the heat-flux battery dominates over the evolving fluid-vorticity term for the fiducial parameters (Villarroel-Sepúlveda et al., 2024).

Other analytic frameworks focus on how already-supplied flux modifies the bulk disc structure. A self-similar model of “flux-frozen” hyper-magnetized quasar disks adopts the closure

$\Phi_{\rm H}$3

arguing that the field is advected inward from the ISM and compressed approximately by flux freezing (Hopkins et al., 2023). The resulting disks are midplane-pressure-dominated by toroidal magnetic fields, with $\Phi_{\rm H}$4, $\Phi_{\rm H}$5, and a large Toomre parameter $\Phi_{\rm H}$6, so they remain gravitationally stable even at hyper-Eddington accretion rates (Hopkins et al., 2023). In another analytic RIAF treatment, allowing quasi-stationary magnetic evolution and restoring terms of order $\Phi_{\rm H}$7 produces a nonzero electric field and a localized Poynting flux near the inner edge, suggesting that inwardly transported ordered flux can power an electrodynamic jet base (Kaburaki, 2011).

5. Radiative-flux implementations and observational usage

In several other subfields, “flux accretion model” refers not to magnetic-flux transport but to models that map accretion dynamics onto observed electromagnetic flux. This usage is explicit in relativistic torus models for high-frequency QPOs, where a slender, polytropic, perfect-fluid, non-self-gravitating torus with constant specific angular momentum orbits in Kerr geometry and produces significant observed flux variability through radial and vertical oscillation modes (Bakala et al., 2014). The mode pair can be taken as $\Phi_{\rm H}$8 for a direct resonance picture or

$\Phi_{\rm H}$9

for a modified relativistic precession picture, with the torus center chosen so the relevant frequencies satisfy a $\dot M$0 ratio (Bakala et al., 2014).

Long-baseline monitoring of Sgr A* uses the phrase in an observationally inferential sense. Re-analysis of Keck/NIRC2 data from 2005–2022 shows that in 2019 the mean luminosity at $\dot M$1 increased by a factor of $\dot M$2, the median $\dot M$3 flux doubled from $\dot M$4 to $\dot M$5 mJy, and the 2019 light curves had higher variance than in all other intervals examined (Weldon et al., 2023). Because the entire flux distribution shifted upward rather than only its bright tail, the paper interprets 2019 as a temporary enhanced-accretion episode, plausibly linked to delayed accretion of tidally stripped gas from G2 (Weldon et al., 2023). The same study finds flux variations over a factor of $\dot M$6, with no evidence for a quasi-steady NIR quiescent state (Weldon et al., 2023).

In blazar X-ray spectroscopy, the phrase can designate a disc–corona accretion model fitted to observed flux. For Mrk 421, the Two Component Advective Flow model yields disc accretion rates of $\dot M$7, halo accretion rates of $\dot M$8, a dynamic-corona size from $\dot M$9 to $\Phi(r,\theta,t)$0, and a Keplerian-disc viscosity parameter of $\Phi(r,\theta,t)$1 during moderate 2017 states (Mondal et al., 2022). The same paper concludes that accretion-disc-based models can contribute to the observed X-ray emission when X-ray data are considered in isolation, but are disfavoured relative to relativistic jet models when multiwavelength information is included (Mondal et al., 2022).

Young-stellar-object applications use the term in a yet different sense: a parameterized time-dependent mapping $\Phi(r,\theta,t)$2. A FU Ori outburst model couples stellar photospheric emission, magnetospheric accretion shocks, an irradiated dust disk, and a viscously heated gas disk, and finds that red optical and near-infrared light curves generally follow the same or very similar form as the input accretion profile, whereas mid-infrared light curves are more responsive to the location and heating of the innermost dust disk (Das et al., 19 May 2026). For the Class I protostar EC 53, radiative-transfer modeling shows that a factor of $\Phi(r,\theta,t)$3 enhancement in $\Phi(r,\theta,t)$4 flux corresponds to an internal luminosity increase by a factor of $\Phi(r,\theta,t)$5, demonstrating that long-wavelength flux changes are muted reprocessed responses rather than direct luminosity scalings (Baek et al., 2020).

6. Open questions, caveats, and contested points

A recurring controversy is whether large-scale flux transport can be reduced to local mean-field coefficients. The coronal-mechanism study explicitly rejects a purely local “effective viscosity versus effective resistivity” picture for large-scale poloidal flux, arguing that the decisive events are global, reconnection-mediated, and topological (0906.2784). Later quasi-local models recover useful transport equations and control parameters, but do so by imposing strong geometric closures, such as slab structure, height-independent $\Phi(r,\theta,t)$6 and $\Phi(r,\theta,t)$7, or $\Phi(r,\theta,t)$8 (Begelman, 2024, Hopkins et al., 2023). This suggests that analytic tractability is often bought at the cost of uncertain closure physics.

Another unresolved issue concerns reconnection. In global GRMHD flux-eruption models, reconnection is numerical rather than explicitly resistive, so the exact rates and locations of topology change are not quantitatively predictive even when the qualitative mechanism appears robust (Porth et al., 2020, Chatterjee et al., 2022). Kinetic GRPIC calculations remedy that deficiency for spherical accretion, but at the price of severe idealization: the published models are two-dimensional, axisymmetric, Schwarzschild, and zero-angular-momentum, with no self-consistent radiation transfer or cooling (Figueiredo et al., 4 Feb 2026). A plausible implication is that current fluid and kinetic models illuminate complementary rather than identical regimes.

Observational identifications remain tentative. The Sgr A* MAD-flare model yields flux bundles with promising magnetic energies and predominantly vertical fields, yet their orbital motion is substantially sub-Keplerian and therefore in tension with the inferred period–radius relation of the GRAVITY flares (Porth et al., 2020). The NIR statistical study of Sgr A* supports an accretion-based interpretation of the 2019 episode, but explicitly does not derive a unique quantitative $\Phi(r,\theta,t)$9 law from NIR flux alone (Weldon et al., 2023). Likewise, the Mrk 421 accretion interpretation is presented as viable for isolated X-ray spectra but not as the preferred broadband source model (Mondal et al., 2022).

Taken together, these caveats indicate that “flux accretion model” is best treated as a technical category rather than a single doctrine. Its strongest and most coherent meaning is the magnetic one: accretion systems transport, concentrate, generate, and sometimes explosively release large-scale flux, and those flux budgets regulate angular momentum transport, vertical support, eruptive behavior, and, in some cases, the observable high-energy output of the source (Chatterjee et al., 2022, Begelman, 2024).