M–Sigma Relation

The M–sigma relation (written symbolically as the M_BH–σ or M–σ relation) is a tight empirical power-law correlation in extragalactic astrophysics between the mass of a central supermassive black hole (M_BH) and the line-of-sight random stellar velocity dispersion (σ) of its host galaxy's spheroidal bulge. Discovered independently in 2000 by Laura Ferrarese and David Merritt and by Karl Gebhardt and collaborators using spatially resolved Hubble Space Telescope (HST) stellar and gas kinematics, the relation follows M_BH ∝ σ^β with an exponent β ≈ 4.0 to 5.1 and an intrinsic logarithmic scatter of only ~0.29–0.34 dex across three orders of magnitude in black hole mass (10⁶ to 10¹⁰ M_☉). Because the gravitational sphere of influence of a central black hole (r_inf = G M_BH / σ²) spans only a few parsecs—roughly 10⁻³ of the effective half-light radius R_e of the host Sérsic profile or Hernquist profile bulge—the M–σ relation provides compelling observational evidence for coevolution and self-regulated active galactic nucleus (AGN) feedback linking nuclear black hole accretion to galaxy-scale star formation and the Faber–Jackson relation.

In logarithmic form normalized to a fiducial velocity dispersion of σ₀ = 200 km s⁻¹, the M–sigma relation is parameterized as:

log₁₀(M_BH / M_☉) = α + β log₁₀(σ_e / 200 km s⁻¹)

where σ_e is the luminosity-weighted line-of-sight stellar velocity dispersion integrated within the host bulge's effective half-light aperture R_e (modeled dynamically via the Jeans equations or Schwarzschild orbit-superposition methods). When first published in 2000, Ferrarese & Merritt obtained β = 4.80 ± 0.54 from 12 galaxies with reliable dynamical black hole masses, while Gebhardt et al. found β = 3.75 ± 0.30 from 26 galaxies using aperture-integrated dispersions. Subsequent resolution of sample selection effects, inclusion of dark-matter halo profiles (such as the Navarro–Frenk–White profile) in triaxial orbit modeling, and expansion of dynamical censuses to ~90–145 galaxies (McConnell & Ma, 2013; Kormendy & Ho, 2013) converged on the modern calibration for elliptical galaxies and classical bulges:

log₁₀(M_BH / M_☉) = (8.49 ± 0.05) + (4.38 ± 0.29) log₁₀(σ_e / 200 km s⁻¹) (Kormendy & Ho, 2013)

When late-type spiral galaxies and brightest cluster galaxies are fitted together across the widest dynamic range, McConnell & Ma (2013) obtain α = 8.32 ± 0.05 and a slightly steeper slope β = 5.64 ± 0.32. Direct dynamical measurements of M_BH underlying this calibration rely on four high-angular-resolution kinematic tracers that resolve the Keplerian velocity rise inside the sphere of influence r_inf = G M_BH / σ² ≈ 10.7 pc (M_BH / 10⁸ M_☉) (σ / 200 km s⁻¹)⁻²:

1. Proper Motions of Individual Stars: Used for Sagittarius A* in the center of the Milky Way (M_BH = 4.15 × 10⁶ M_☉, σ ≈ 105 km s⁻¹).
2. Circumnuclear Water (H₂O) Megamaser Disks: Very Long Baseline Interferometry (VLBI) at 22 GHz traces thin, sub-parsec Keplerian molecular disks (e.g., NGC 4258 with M_BH = 4.0 × 10⁷ M_☉), providing geometric distance and mass precision better than 3%.
3. Spatially Resolved Stellar & Ionized Gas Kinematics: HST and adaptive-optics integral-field spectroscopy of nearby ellipticals and bulges such as the Andromeda Galaxy (M31, M_BH ≈ 1.4 × 10⁸ M_☉) and M87 (M_BH = 6.5 × 10⁹ M_☉, confirmed by Event Horizon Telescope horizon-scale interferometry).
4. Reverberation Mapping of Active Galactic Nuclei: Measures the light-travel time delay Δt = R_BLR / c between continuum flares and broad emission-line responses in distant quasars, calibrated against the local M–σ zero point via the virial factor f.

Because the gravitational potential of a central black hole is negligible at kiloparsec radii where σ_e is established, the tightness of the M–σ relation is explained by quasar-mode (radiative) and radio-mode (mechanical) feedback. Even though M_BH is only ~0.15% to 0.5% of the host bulge stellar mass M_bulge (the Magorrian relation), the relativistic rest-mass energy released during black hole accretion (E_BH = ε_r M_BH c² with radiative efficiency ε_r ≈ 0.1) exceeds the gravitational binding energy of the entire host bulge (E_bind ≈ M_bulge σ²) by a factor of ~10² for σ = 200 km s⁻¹. If even a few percent of the Eddington luminosity L_Edd = 4π G M_BH m_p c / σ_T couples to the cold interstellar gas fueling star formation via the Kennicutt–Schmidt law, the resulting outflow expels the ambient gas, terminating both black hole growth and bulge star formation simultaneously.

Two limiting thermodynamic regimes of an AGN wind sweeping up an isothermal gas bulge with gas fraction f_g produce the observed power-law slopes β = 4 and β = 5:

1. Momentum-Conserving (Radiatively Cooled) Outflow (M_BH ∝ σ⁴): Proposed by Andrew King (2003) and Fabian (1999), within the inner cooling radius of the bulge, inverse-Compton cooling by the quasar radiation field drains thermal energy from the shocked wind faster than the flow expansion time. Only the outward radiation momentum flux L_Edd / c acts on the swept-up gas shell. Setting the outward wind thrust equal to the inward gravitational weight of the gas shell in an isothermal potential yields the analytical threshold mass:

M_BH,σ = (f_g κ_es / π G²) σ⁴ ≈ 1.5 × 10⁸ M_☉ (σ / 200 km s⁻¹)⁴ which reproduces both the slope β = 4 and the empirical normalization α ≈ 8.2–8.5 without adjustable free parameters.

1. Energy-Conserving (Adiabatic) Outflow (M_BH ∝ σ⁵): First derived by Joseph Silk and Martin Rees (1998), at larger radii where the shocked wind does not cool radiatively, the thermalized mechanical wind luminosity f_w L_Edd performs adiabatic p dV work on the ambient interstellar medium. Balancing the wind energy deposition rate against the rate work must be done to unbind the bulge gas (~f_g σ⁵ / G) yields the steeper scaling M_BH ∝ σ⁵, which also explains why high-mass galaxies above the knee of the Schechter luminosity function are quenched into red-sequence ellipticals obeying the Kormendy relation.

At the highest masses (M_BH > 10^9.5 M_☉) in Brightest Cluster Galaxies (BCGs), stellar velocity dispersion saturates near σ ≈ 300–380 km s⁻¹ because dry dissipationless mergers increase M_BH and R_e while leaving the virial velocity dispersion σ² ≈ G M / R_e nearly constant; consequently, ultramassive black holes sit systematically above the extrapolation of the M–σ relation while remaining aligned with the M_BH–L_bulge relation. Conversely, at the low-mass end (M_BH ≈ 10⁵–10⁷ M_☉), disk-grown pseudobulges and bulgeless dwarf galaxies exhibit lower black hole masses and broader scatter below the classical M–σ line, reflecting inefficient secular gas feeding rather than merger-triggered Eddington-limited feedback.

What is the sphere of influence of a supermassive black hole?

The gravitational sphere of influence r_inf ≡ G M_BH / σ² is the characteristic radius within which the orbital motion of stars and gas is dominated by the point-mass Keplerian potential of the central black hole rather than the distributed stellar mass of the galaxy bulge. Reliable dynamical black hole mass determinations require angular resolution comparable to or smaller than r_inf / D, where D is the galaxy distance.

How does the M–sigma relation relate to the Faber–Jackson relation?

The Faber–Jackson relation demonstrates that the total stellar luminosity (and stellar mass M_bulge) of an elliptical galaxy scales with velocity dispersion as L ∝ σ⁴. Combining M_bulge ∝ σ⁴ with the momentum-driven M–σ scaling M_BH ∝ σ⁴ directly implies a linear proportionality between black hole mass and host bulge mass (M_BH ≈ 0.002–0.005 M_bulge), known as the Magorrian relation.

Why do Brightest Cluster Galaxies deviate upward from the M–sigma relation?

When two gas-poor ("dry") elliptical galaxies merge at the center of a galaxy cluster, their central supermassive black holes inspiral via dynamical friction and coalesce so that M_BH doubles, while orbital energy puffs up the stellar effective radius R_e in proportion to mass M. Because σ² ∝ M / R_e stays nearly unchanged during dry mergers, M_BH grows upward at nearly constant σ.

1. Ferrarese, L., & Merritt, D. (2000). "A fundamental relation between supermassive black holes and their host galaxies." The Astrophysical Journal, 539(1), L9–L12. https://doi.org/10.1086/312838
2. Gebhardt, K., Bender, R., Bower, G., Dressler, A., Faber, S. M., Filippenko, A. V., Green, R., Grillmair, C., Ho, L. C., Kormendy, J., Lauer, T. R., Magorrian, J., Pinkney, J., Richstone, D., & Tremaine, S. (2000). "A relationship between nuclear black hole mass and galaxy velocity dispersion." The Astrophysical Journal, 539(1), L13–L16. https://doi.org/10.1086/312840
3. King, A. (2003). "Black holes, outflows, and the M–σ relation." The Astrophysical Journal, 596(1), L27–L29. https://doi.org/10.1086/379143
4. McConnell, N. J., & Ma, C.-P. (2013). "Revisiting the scaling relations of black hole masses and host galaxy properties." The Astrophysical Journal, 764(2), 184. https://doi.org/10.1088/0004-637X/764/2/184
5. Kormendy, J., & Ho, L. C. (2013). "Coevolution (or not) of supermassive black holes and host galaxies." Annual Review of Astronomy and Astrophysics, 51, 511–653. https://doi.org/10.1146/annurev-astro-082708-101811