Semi-Analytic Inflow Framework

• Semi-analytic inflow frameworks are computationally tractable models that describe gas dynamics using a compact set of differential equations with tunable physical parameters.
• They incorporate key processes such as gravity, pressure gradients, cooling, and feedback, enabling rapid exploration of galaxy accretion, chemical enrichment, and black hole environments.
• These frameworks effectively reproduce observable signatures from star formation histories to horizon-scale imaging, linking theoretical models with astrophysical observations.

A semi-analytic inflow framework refers to a class of physically motivated, computationally tractable models that describe the dynamics, thermodynamics, or chemistry of inflowing gas (or particles) in an astrophysical or galactic context. Unlike fully numerical simulations, these frameworks encode the relevant physical processes into a compact set of equations—ordinary differential equations, or population-based rate equations for ensembles—often with explicit analytic or semi-analytic solutions. They are widely used to model radiative accretion flows, gas supply to galaxies and black holes, and circumgalactic media, with applications extending to interpreting observations from instruments such as the Event Horizon Telescope or large-scale spectroscopic galaxy surveys.

1. Basic Concepts and Governing Equations

A semi-analytic inflow framework aims to describe how gas (or clouds, or plasma) accretes into a physical system under the influence of key physical processes such as gravity, pressure gradients, drag, cooling, conduction, mixing, and feedback. The framework is defined by a set of governing equations, whose structure depends on the target system:

• Galaxy Accretion & Chemical Evolution: One-zone models track the gas mass $M_g(t)$, the stellar mass $M_*(t)$, and gas metallicity $Z(t)$ in a system. Typical ODE systems include:

$$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$

$\dot M_* = (1-R)\psi(t)$

$$\dot M_Z = Z_{\rm in}\dot M_{\rm in}(t) - Z(1-R)\psi(t) + (1-R)y_Z\psi(t) - Z\dot M_{\rm out}(t)$$

with inflow rates $\dot M_{\rm in}$ often parameterized as an exponential decay or as a fixed multiple of the SFR (Zhou et al., 2022, Yabe et al., 2014).

• Circumgalactic Clouds: The evolution of pressure-confined clouds is described by local pressure equilibrium, $n_cT_c=n_hT_h(r)$, with dynamics governed by gravity and drag:

$$M_c\frac{d^2r}{dt^2} = -\frac{GM_{\rm DM}(<r)M_c}{r^2} - F_{\rm ram}(r)$$

and mass-loss rates via conduction-evaporation models (Lan et al., 2018).

• Accretion Flows Near Compact Objects: In the context of black hole accretion, the dynamical variables include the four-velocity $u^\mu$, density profile $M_*(t)$0, and magnetic field geometry, often solved using mass continuity, stationary Euler equations (momentum balance), and a specified closure (e.g., polytropic or isothermal):

$M_*(t)$1

with additional prescriptions for $M_*(t)$2, $M_*(t)$3, and boundary conditions at the event horizon (Saurabh et al., 15 Aug 2025, Ogihara et al., 2024).

These equations are reduced to a computationally manageable system by adopting symmetry assumptions (e.g., axisymmetry, stationarity), closure relations (e.g., Schmidt law for star formation), or by introducing key dimensionless parameters (e.g., wind mass-loading factor, conduction suppression factor, inflow timescale).

2. Primary Physical Ingredients and Parameterization

Semi-analytic inflow frameworks are characterized by the incorporation of physically motivated but tunable parameters to account for complex processes unresolved by analytic theory:

• Inflow Rate Parameterization:
  • Exponential-onset: $M_*(t)$4 for $M_*(t)$5 (Zhou et al., 2022).
  • Proportional-to-SFR: $M_*(t)$6, with $M_*(t)$7 a free parameter (Yabe et al., 2014).
  • Fixed mass-flux for ensembles: $M_*(t)$8 for CGM clouds (Lan et al., 2018).
• Mass-loss and Feedback: Explicit models for outflow $M_*(t)$9, or evaporation/conduction terms with saturation and suppression corrections ($Z(t)$0, $Z(t)$1), are included to regulate the fate of inflowing matter (Lan et al., 2018, Zhou et al., 2022).
• Initial Condition Specification: Frameworks often define an initial cloud mass function (Schechter-like or power-law), radial velocity distributions, or delay times for the onset of inflow.
• Magnetic and Geometrical Configuration: For accretion flows, magnetic field geometry is incorporated via a vector potential $Z(t)$2, with configurations interpolating between monopole, parabolic, and vertical field topologies, and a geometry parameter $Z(t)$3 distinguishing toroidal and poloidal dominance (Saurabh et al., 15 Aug 2025).
• Gas Dynamics and Thermodynamics: The frameworks account for the physical conditions governing sound speed, disk thickness $Z(t)$4, electron temperature prescriptions, and advection-dominated energy transport as appropriate.

The table below illustrates a selected mapping of parameterizations across various implementations:

Framework Target	Inflow Parameter	Outflow/Feedback
One-zone galactic SFH (Zhou)	Exp. decline: $Z(t)$5	$Z(t)$6
"Leaky box" galaxy (Yabe)	$Z(t)$7	$Z(t)$8
CGM cloud ensemble (Lan & Mo)	Fixed $Z(t)$9, mass Fn	Evaporation $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$0, $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$1
RIAF/BH inflow (EHT)	Flow-mixing $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$2	N/A

These parameterizations enable the modeler to flexibly explore how variations in gas supply, retention, and loss shape global system properties and observables.

3. Methodologies for Implementation and Solution

The methodological workflow in semi-analytic inflow frameworks generally follows a staged architecture:

1. Specification of Host System Properties: Quantities such as halo mass, concentration, disk radius, input gravitational potential (e.g., NFW or Kerr metric), and gas/magnetic field distributions are set.
2. Parameter Selection and Grid Construction: Choices for inflow rates, mass functions, wind efficiencies, conduction suppression factors, and magnetic geometry are made according to the model context (Lan et al., 2018, Zhou et al., 2022).
3. Pre-computation of Local Physical Properties: For multi-phase or ensemble models, pre-computed grids (e.g., of ionization fractions, $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$3 from CLOUDY models) are established as a function of radius, cloud mass, or other relevant variables.
4. Numerical Integration or Ensemble Evolution: The governing ODEs or system of rate equations are advanced forward, typically employing explicit integrators. Individual inflow events/objects are created and tracked, updating positions, masses, and observables at each timestep.
5. Observables Extraction: Model snapshots are analyzed to generate mock observables:
   • Covering fractions, kinematics, and column densities for CGM (impact parameter analysis) (Lan et al., 2018).
   • SFH, chemical evolution, SED fitting for integrated galaxy properties (via stellar population synthesis and likelihood maximization) (Zhou et al., 2022).
   • Synthetic images and polarization metrics for accretion flows using GR ray-tracing and polarized radiative transfer (Saurabh et al., 15 Aug 2025, Ogihara et al., 2024).
6. Statistical Model Fitting: Likelihood or $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$4 minimization, often with nested sampling (e.g., MultiNest), is used to fit model predictions to observed data (spectra, gas fractions, metallicities, imaging metrics).

4. Principal Results and Observational Comparisons

Semi-analytic inflow models have produced key results across physical scales:

• Gas Accretion in Galaxies: The best-fit gas inflow rates in star-forming galaxies at $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$5 are constrained as $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$6 times the SFR, with outflow rates $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$7 (Yabe et al., 2014). Both $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$8 and $$\dot M_g = \dot M_{\rm in}(t) - (1-R)\psi(t) - \dot M_{\rm out}(t)$$9 decline with decreasing redshift, paralleling the decline in cosmological accretion and wind activity.
• Star Formation and Chemical Enrichment Histories: Application to MaNGA spiral galaxies yields:
  • Infall onset $\dot M_* = (1-R)\psi(t)$0 Gyr,
  • Infall timescale $\dot M_* = (1-R)\psi(t)$1 Gyr decreasing with stellar mass,
  • Wind parameter $\dot M_* = (1-R)\psi(t)$2 falling from $\dot M_* = (1-R)\psi(t)$310 at low mass to $\dot M_* = (1-R)\psi(t)$45 at high mass.
  • These results reproduce the mass--metallicity relation, present-day $\dot M_* = (1-R)\psi(t)$5, and "downsizing" in galaxy formation (Zhou et al., 2022).
• Circumgalactic Pressure-confined Clouds: The inflow model can reproduce observed covering fractions of high-$\dot M_* = (1-R)\psi(t)$6 gas in passive galaxies but fails to match the full range of low-$\dot M_* = (1-R)\psi(t)$7 systems (Lan et al., 2018). Outflow-driven models, in contrast, better match the combination of covering fraction, spatial distribution, and velocity dispersion in star-forming galaxies.
• Inflow in Black Hole Accretion Environments: Semi-analytic RIAF models in Kerr spacetime with varying magnetic field geometry, disk thickness, and inflow dynamics quantitatively reproduce synthetic EHT-observables such as ring diameter, brightness asymmetry, polarization fraction, and EVPA structure in M87*. Flows with a poloidal field dominated configuration and partially radial inflow best match the observed metrics (Saurabh et al., 15 Aug 2025). The impact of disk thickness on observables is minor; the magnetic geometry and inflow prescription are most significant.
• Jets and Magnetized Funnel Regions: The structure of radio images in highly magnetized jet funnels is dominated by the location of the inflow-outflow separation surface, leading to multiple ring-like features in synthetic images and providing a direct testable prediction for next-generation EHT observations (Ogihara et al., 2024).

5. Assumptions, Sensitivities, and Limitations

Every semi-analytic inflow framework entails a suite of simplifying assumptions, sensitivities to parameter choices, and modeling boundaries:

• Mixing and Homogenization: One-zone galaxy models assume instantaneous, homogeneous mixing, with no radial flows or environmental interactions (Zhou et al., 2022, Yabe et al., 2014).
• Metallicity of Inflow: Inflowing material is typically assumed to be pristine ($\dot M_* = (1-R)\psi(t)$8) (Yabe et al., 2014).
• Scaling Laws: The prescribed star formation law (e.g., linear Schmidt law or Kennicutt law) fixes the conversion of gas mass to SFR, with uncertainties in parameters such as return fraction $\dot M_* = (1-R)\psi(t)$9 and yield $$\dot M_Z = Z_{\rm in}\dot M_{\rm in}(t) - Z(1-R)\psi(t) + (1-R)y_Z\psi(t) - Z\dot M_{\rm out}(t)$$0 influencing the inferred inflow/outflow efficiencies.
• Feedback and Instabilities: For circumgalactic clouds, conduction suppression (fiducially $$\dot M_Z = Z_{\rm in}\dot M_{\rm in}(t) - Z(1-R)\psi(t) + (1-R)y_Z\psi(t) - Z\dot M_{\rm out}(t)$$1), the assumption of pressure equilibrium, and the neglect of explicit hydrodynamic instabilities (which are assumed suppressed by cooling or magnetic fields) systematically affect cloud lifetimes and observable covering fractions (Lan et al., 2018).
• Angular Momentum and Multiphase Effects: Most inflow models for galaxies disregard angular momentum transport except in analytic models of torque-driven inflow to black holes, where the amplitude of non-axisymmetric perturbations and the shock criterion set the inflow rates (Hopkins et al., 2010).
• Magnetic and Spin Sensitivities: In RIAF and GRMHD-based models, the geometry of the magnetic field and black hole spin directly modulate key image and polarization properties, with the separation surface between inflow and outflow fixed primarily by spin (Saurabh et al., 15 Aug 2025, Ogihara et al., 2024).

A plausible implication is that the predictive power of semi-analytic inflow frameworks, while high, is contingent on the calibration of these free parameters to observed galaxies, clouds, or accretion systems and on the fidelity of physical assumptions regarding feedback, phase structure, and angular momentum transport.

6. Applications and Extensions

Semi-analytic inflow frameworks serve as computational laboratories for rapid exploration of parameter spaces, efficient model fitting, and physical intuition building across scales. Their applications include:

• Reconstruction of Galactic Star Formation and Metallicity Histories from spectra and emission line data, providing insight into downsizing and feedback processes (Zhou et al., 2022).
• Interpretation of EHT Horizon-Scale Images for black hole environments, through parameter studies of magnetic field geometry, disk thickness, and inflow dynamics (Saurabh et al., 15 Aug 2025, Ogihara et al., 2024).
• Modeling of Circumgalactic and Halo Gas Accretion, Feedback, and Cloud Lifetimes, guiding the design and interpretation of CGM absorber surveys and cosmological hydrodynamic simulations (Lan et al., 2018).
• Analytic Sub-grid Prescriptions for Gas and BH Accretion in cosmological simulations, leveraging validated inflow/shock criteria from semi-analytic torque models (Hopkins et al., 2010).

This flexibility, combined with favorable computational performance, underpins the broad utility of semi-analytic inflow frameworks in modern astrophysical research. The approach provides vital links between simulation, analytic theory, and observation, particularly in regimes where full-scale simulations are either infeasible or lack the necessary physical detail.