Eddington Ratios in Black Hole Accretion

• Eddington ratios are dimensionless parameters that quantify black hole accretion efficiency by comparing observed bolometric luminosity to theoretical Eddington luminosity.
• Measurement techniques involve multiwavelength estimators and single-epoch virial methods, with uncertainties typically around 0.3–0.5 dex.
• The ratio influences AGN feedback and evolution, distinguishing radiatively inefficient, thin-disk, and super-Eddington accretion regimes.

The Eddington ratio is a dimensionless parameter that quantifies the normalized accretion rate of a compact object, typically a black hole. Defined as the ratio of the object's bolometric (total radiative) luminosity to its theoretical Eddington luminosity, it captures how efficiently mass is being accreted relative to the point where outward radiation pressure balances inward gravitational attraction. In black hole astrophysics, Eddington ratios are foundational for interpreting accretion physics, feedback, AGN demographics, and cosmic evolution.

1. Definition, Formalism, and Physical Interpretation

The classical Eddington luminosity, $L_{\rm Edd}$, is derived for fully ionized hydrogen by balancing gravitational force against radiation pressure from Thomson scattering: $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$ where $M_{\rm BH}$ is the black hole mass, $m_p$ the proton mass, $c$ the speed of light, $\sigma_T$ the Thomson cross-section, and $G$ Newton's constant (Izumi, 2018).

The Eddington ratio is then defined as: $$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$ where $L_{\rm bol}$ is the observed bolometric luminosity. For AGN, $L_{\rm bol}$ is typically estimated via multiwavelength SED integration or via bolometric corrections to monochromatic fluxes (e.g., $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$0, $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$1, X-ray bands, or emission-line luminosities) (Gupta et al., 16 Jul 2025, Kim et al., 2024, Kong et al., 2018). $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$2 traces the instantaneous accretion vigor: $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$3 signals radiatively inefficient flows (e.g., ADAF/RIAF), $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$4 characterizes thin/disk accretion, and $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$5 marks super-Eddington regimes with powerful radiation-driven outflows and geometrically thick disks (Panda et al., 1 Oct 2025, Zhang et al., 2 Jun 2025).

2. Measurement Techniques and Multiwavelength Estimators

Accurate determination of $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$6 requires robust measurements of both $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$7 and $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$8. Methodologies include:

• Single-epoch virial methods: $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$9 estimated from the BLR radius-luminosity relation and broad-line FWHM (typically H$M_{\rm BH}$0, Mg II, or CIV), combined with empirical virial factors (Izumi, 2018, Farina et al., 2022).
• Stellar velocity dispersion ($M_{\rm BH}$1) methods: For narrow-line (type 2) AGN, $M_{\rm BH}$2 derived from $M_{\rm BH}$3–$M_{\rm BH}$4 relations, often using absorption-line spectra or kinematic proxies from forbidden lines (Kong et al., 2018).
• Bolometric corrections: $M_{\rm BH}$5 inferred from monochromatic optical, MIR ($M_{\rm BH}$6), X-ray, or [O III] emission, applying empirically calibrated scaling relations (Kim et al., 2024, Kim et al., 2015, Gupta et al., 16 Jul 2025).
• Dust-immune proxies: MIR-based estimators (e.g., Kim et al. 2023) allow measurement of $M_{\rm BH}$7 in highly obscured quasars, mitigating extinction biases and extending applicability to merger-driven, dusty environments (Kim et al., 2024).
• Spectroscopic indicators: Line ratios such as [N II]/H$M_{\rm BH}$8 provide empirical proxies for $M_{\rm BH}$9 in cases where direct mass estimators are unavailable, with calibration scatter $m_p$00.6 dex (Oh et al., 2016, Oh et al., 2019).
• Far-UV iron features: The Fe II+Fe III $m_p$11123 emission shows strong correlation with optical Fe II/H$m_p$2, itself a classical $m_p$3 indicator, allowing photometric inference at low $m_p$4 (Zheng, 2021).

Uncertainties typically range from 0.3–0.5 dex, driven by line width systematics, bolometric correction variance, virial factor calibration, and extinction corrections. Sample selection further complicates the observed distribution of $m_p$5, frequently biasing against both low and super-Eddington accretors (Shirakata et al., 2019, Schulze et al., 2010).

3. Eddington Ratio Distribution Functions (ERDF): Shape, Evolution, and Dependencies

Intrinsic ERDFs are well modeled as broad log-normal or Schechter-like functions, with a rising low-$m_p$6 tail and an exponential cutoff near unity:

• ERDF empirical fits: For local broad-line AGN, the intrinsic ERDF follows $m_p$7 below a cutoff at $m_p$8 (Schulze et al., 2010). For obscured AGN (X-ray selected), double-power-law forms best fit the distribution, with breakpoints and slopes sensitive to obscuration class (Ananna et al., 2022).
• Redshift trends: The ERDF "knee," $m_p$9, increases with increasing $c$0, indicating "downsizing": high-$c$1 AGN preferentially accrete at near- or super-Eddington rates, while local AGN display broader, lower-$c$2 distributions (Shirakata et al., 2019, Zhang et al., 2023, Carraro et al., 2022). Quantitatively, super-Eddington duty cycles ($c$3) climb from $c$4 at $c$5 to $c$6 at $c$7 for $c$8 (Shirakata et al., 2019).
• Mass dependence: Lower-mass black holes exhibit higher typical $c$9 at fixed redshift, intensifying downsizing effects. At $\sigma_T$0, the active fraction above $\sigma_T$1 falls steeply with $\sigma_T$2: virtually no $\sigma_T$3 black hole accretes at significant levels locally (Schulze et al., 2010).
• Host gas fraction correlation: Empirically, local quasars with higher $\sigma_T$4 preferentially reside in gas-rich hosts ($\sigma_T$5), indicating a direct link between accretion rate and molecular fueling availability (Izumi, 2018).
• Selection biases: Flux/luminosity thresholds in surveys severely distort the observed ERDF, suppressing both low-$\sigma_T$6 and high-$\sigma_T$7 wings, emphasizing the necessity of selection-function-aware modeling (Shirakata et al., 2019, Schulze et al., 2010).

ERDF Feature	Physical Interpretation	Evidence
Log-normal or Schechter form	Many stochastic fueling channels	(DeGraf et al., 2012, Shankar et al., 2011)
Redshift-dependent $\sigma_T$8	Cosmological gas fraction evolution	(Shirakata et al., 2019, Zhang et al., 2023)
Mass-dependent shifts	Anti-hierarchical AGN downsizing	(Shankar et al., 2011, Schulze et al., 2010)

4. AGN Physics and Feedback Regimes as a Function of Eddington Ratio

Accretion flow structure and observational properties of AGN vary systematically with $\sigma_T$9:

• Sub-Eddington ($G$0): Accretion flows are radiatively inefficient, geometrically thick, and may manifest as low/hard states with strong hard X-ray emission, weak UV bump, and compact radio jets. The X-ray bolometric correction remains flat with luminosity, indicative of corona-dominated emission (Gupta et al., 16 Jul 2025, Zhang et al., 2 Jun 2025).
• Thermal/Thin Disk ($G$1): Standard Shakura-Sunyaev disks dominate, with efficient radiative output, prominent big blue bump, and strong Fe II emission. Bolometric corrections increase with $G$2 (Panda et al., 1 Oct 2025, Kim et al., 2024).
• Super-Eddington ($G$3): Disks become radiation pressure supported ("slim disks"), drive dense winds and powerful outflows, often characterized by low radiative efficiencies ($G$4 for $G$5) (Zhang et al., 2 Jun 2025, Zhang et al., 2023). AGN with $G$6 dominate cosmic SMBH growth at early epochs and are correlated with high host gas fractions (Izumi, 2018, Shirakata et al., 2019). Feedback in super-Eddington regimes effectively curtails further fueling, linking accretion rates to AGN self-regulation.

A critical transition occurs near $G$7, distinguishing radiatively inefficient from efficient accretion and marking boundaries for X-ray spectral behavior, variability, and jet production (Gupta et al., 16 Jul 2025, Panda et al., 1 Oct 2025).

5. Obscuration, Host Properties, and Evolutionary Pathways

The Eddington ratio is central to AGN unification, feedback, and demographic evolution:

• Obscured vs. unobscured AGN: Dust-obscured quasars exhibit systematically higher $G$8 than unobscured ones across $G$9, corresponding to early, rapid-growth phases in merger-driven evolutionary scenarios. Geometric unification would not predict such an offset, confirming the evolutionary nature (Kim et al., 2024, Kim et al., 2015).
• Obscured fraction vs. $$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$0: Covering factors (fraction of obscured AGN) are high ($$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$10.8) at low $$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$2, but drop precipitously above $$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$3 due to radiation pressure clearing of dusty circumnuclear material (Ananna et al., 2022). Temporal "clearing" models accurately reproduce the observed covering fraction with exponential or sigmoid decays.
• Feedback signatures: High $$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$4 AGN exhibit strong blue-shifted emission-line wings linked to radiation-driven outflows and broad-line region scale heights that exceed depleted torus heights, making BLR emission visible from all sightlines at high accretion rates (Ananna et al., 2022).
• Host coevolution: The positive correlation between $$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$5 and host gas fraction ($$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$6) supports models of coeval black hole–galaxy evolutionary tracks, merger-driven fueling, and AGN downsizing (Izumi, 2018, Shirakata et al., 2019).

6. Broader Impact: Cosmology, Variability, and Multiwavelength Diagnostics

The Eddington ratio is pivotal in both fundamental AGN research and practical observational applications:

• Radius-luminosity standardization: Correcting the BLR $$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$7–$$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$8 relation for $$\lambda_{\rm Edd} \equiv \frac{L_{\rm bol}}{L_{\rm Edd}}$$9 reduces intrinsic scatter, enabling use of quasars as cosmological probes and informing on the Hubble tension (Panda et al., 1 Oct 2025).
• Variability: The anti-correlation between fractional optical/UV variability ($L_{\rm bol}$0) and $L_{\rm bol}$1 is robust, providing a photometric estimator for accretion state applicable to surveys (e.g. ZTF, LSST) (Panda et al., 1 Oct 2025).
• SED and bolometric corrections: $L_{\rm bol}$2 governs multiwavelength AGN SEDs, dictating soft X-ray excess strength, UV/X-ray ratios, and bolometric corrections. The X-ray bolometric correction, $L_{\rm bol}$3, is primarily regulated by $L_{\rm bol}$4 not $L_{\rm bol}$5, with a breakpoint at $L_{\rm bol}$6 indicating a transition in accretion physics (Gupta et al., 16 Jul 2025).
• High-redshift SMBH formation: ERDF evolution with $L_{\rm bol}$7 constrains formation pathways for massive early black holes, indicating that early quasars accrete at high $L_{\rm bol}$8, often requiring super-Eddington growth or heavy seeds ($L_{\rm bol}$9) to reach observed masses by $L_{\rm bol}$0 (Farina et al., 2022, Zhang et al., 2023).
• Changing-Look AGN and state transitions: Rapid state changes (on/off transitions) in CLAGN cluster around low $L_{\rm bol}$1, confirming the role of accretion rate shifts in dramatic spectral variability and multi-epoch emission-line changes (Panda et al., 1 Oct 2025).

7. Summary Table: Key Empirical Relations

Context	Equation/Fit	Reference
Eddington ratio	$L_{\rm bol}$2	(Izumi, 2018, Panda et al., 1 Oct 2025)
X-ray bolometric correction	$L_{\rm bol}$3	(Gupta et al., 16 Jul 2025)
[N II]/H$L_{\rm bol}$4 proxy	$L_{\rm bol}$5	(Oh et al., 2016)
Dust-obscured H$L_{\rm bol}$6	$L_{\rm bol}$7	(Kong et al., 2018)
MIR bolometric estimator	$L_{\rm bol}$8	(Kim et al., 2024)
Host gas fraction correlation	$L_{\rm bol}$9	(Izumi, 2018)

• (Izumi, 2018) Supermassive black holes with higher Eddington ratios preferentially form in gas-rich galaxies
• (Gupta et al., 16 Jul 2025) The Eddington Ratio as the Primary Regulator of the Fraction of X-ray Emission in Active Galactic Nuclei
• (Ananna et al., 2022) Probing the Structure and Evolution of BASS AGN through Eddington Ratios
• (Shirakata et al., 2019) Slowing down of cosmic growth of supermassive black holes: Theoretical prediction of the Eddington ratio distribution
• (Kim et al., 2024) Eddington Ratios of Dust-obscured Quasars at $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$00: Evidence Supporting Dust-obscured Quasars as Young Quasars
• (Panda et al., 1 Oct 2025) Feeding frenzy in the mighty black holes: what we could learn from them?
• (Zhang et al., 2 Jun 2025) Radiation GRMHD Models of Accretion onto Stellar-Mass Black Holes: I. Survey of Eddington Ratios
• (Zhang et al., 2023) TRINITY III: Quasar Luminosity Functions Decomposed by Halo, Galaxy, and Black Hole Masses and Eddington Ratios from z=0-10
• (Schulze et al., 2010) Low redshift AGN in the Hamburg/ESO Survey: II. The active black hole mass function and the distribution function of Eddington ratios
• (Kim et al., 2015) Accretion Rates of Red Quasars from the Hydrogen P$$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$01 line
• (DeGraf et al., 2012) Growth of Early Supermassive Black Holes and the High-Redshift Eddington Ratio Distribution
• (Oh et al., 2016) BAT AGN Spectroscopic Survey-III. An observed link between AGN Eddington ratio and narrow emission line ratios
• (Farina et al., 2022) The X-shooter/ALMA Sample of Quasars in the Epoch of Reionization. II. Black Hole Masses, Eddington Ratios, and the Formation of the First Quasars
• (Zheng, 2021) Far-UV Fe Emission as Proxy of Eddington Ratios
• (Kong et al., 2018) The Black Hole Masses and Eddington Ratios of Type 2 Quasars
• (Carraro et al., 2022) An Eddington ratio-driven origin for the $$L_{\rm Edd} = \frac{4\pi G M_{\rm BH} m_p c}{\sigma_T} \simeq 1.26 \times 10^{38} \left(\frac{M_{\rm BH}}{M_\odot}\right) \mathrm{erg~s}^{-1}$$02 relation in quiescent and star forming active galaxies
• (Shankar et al., 2011) Accretion-Driven Evolution of Black Holes: Eddington Ratios, Duty Cycles, and Active Galaxy Fractions

The Eddington ratio, through its central role in accretion physics, feedback, AGN demographics, and multiwavelength observables, remains a cornerstone quantity for deciphering the growth and coevolution of black holes and their host galaxies.