Dynamical Friction

Dynamical friction (also called Chandrasekhar dynamical friction or gravitational drag) is a fundamental collective dissipation process in stellar dynamics, celestial mechanics, and cosmology in which a massive body of mass M (such as a supermassive black hole, a globular star cluster, or an orbiting satellite galaxy) moving with velocity v_M through a collisionless background sea of lighter particles of mass m ≪ M (stars or dark matter particles) experiences a systematic retarding force opposite to its direction of motion—without any physical contact or hydrodynamic gas collisions. First derived in 1943 by Indian-American theoretical astrophysicist Subrahmanyan Chandrasekhar, dynamical friction has a simple physical origin in two dimensions of view:

1. In the Two-Body Scattering Picture: As the massive body M plows through the background sea of field stars, its gravity deflects each passing star (m) along a hyperbolic Rutherford trajectory around M. In the center-of-mass frame, elastic gravitational scattering rotates the relative velocity vector without changing its magnitude; transformed back into the galaxy rest frame, these cumulative two-body deflections transfer forward momentum and kinetic energy from the heavy object M into random velocity dispersion of the light background stars.
2. In the Gravitational Wake Picture (Mulder, 1983): Because deflected background stars are pulled toward the trajectory of M as it passes by, they converge into a tear-drop-shaped gravitational density wake (δρ > 0) trailing directly behind the moving mass M. The backward gravitational pull of this overdense wake on M itself is dynamical friction (Chandrasekhar, 1943; Mulder, 1983; Binney & Tremaine, 2008).

Consider a test body of mass M moving at velocity v_M through an infinite, homogeneous, isotropic background of field particles of individual mass m and mass density ρ = n m. When a single field particle with initial velocity v_m passes M at impact parameter b with relative speed V₀ = |v_M − v_m|, hyperbolic Keplerian scattering deflects the relative velocity vector by angle θ_def, where tan(θ_def / 2) = G(M + m) / (b V₀²) ≡ b₉₀ / b (with b₉₀ = G(M + m) / V₀² denoting the 90°-deflection impact parameter).

Integrating the parallel velocity change Δv_M,∥ over all impact parameters from b = 0 up to a maximum cutoff scale b_max (the size of the host galaxy or dark matter halo) and over a Maxwellian velocity distribution of field stars with one-dimensional velocity dispersion σ (so that f(v_m) = n (2πσ²)^−3/2 exp[−v_m² / (2σ²)])—and invoking Newton's Shell Theorem in velocity space (whereby only field particles moving slower than M, v_m < v_M, contribute net drag)—yields Chandrasekhar's (1943) dynamical friction equation:

dv_M / dt = − [ 4π G² M ρ ln Λ / v_M³ ] [ erf(X) − (2X / √π) exp(−X²) ] v_M

where:

• X = v_M / (√2 σ) is the dimensionless ratio of the body's speed to the background velocity dispersion.
• erf(X) = (2 / √π) ∫₀^X exp(−u²) du is the Gauss error function, so the bracketed factor [erf(X) − (2X/√π)exp(−X²)] represents the fraction of background field particles moving slower than v_M.
• ln Λ = ln( b_max / b₉₀ ) ≈ ½ ln( 1 + Λ² ) is the dimensionless Coulomb logarithm, where Λ = b_max V_typ² / [ G(M + m) ] ∼ (R_sys / R_M) (M_sys / M) typically takes values of ln Λ ≈ 3 to 15 across astrophysical systems.

Key Scaling Laws of Chandrasekhar Drag

Examining Chandrasekhar's formula reveals four counterintuitive physical properties (Binney & Tremaine, 2008):

1. The Drag Force Scales as F_df = M |dv_M/dt| ∝ M² ρ: Because the gravitational wake amplitude δρ/ρ is proportional to the perturber's mass M, and the gravitational force between M and its wake is proportional to M × δρ, the retarding force scales with the square of the moving body's mass (F_df ∝ M²), and the deceleration |dv_M/dt| is proportional to M. Consequently, more massive satellites, star clusters, and black holes sink to the center of a galaxy much faster than less massive ones (mass segregation).
2. Independence from Background Particle Mass (m): So long as m ≪ M, the deceleration depends on the total background mass density ρ = n m, not on the individual mass m of the background particles. A 10⁹ M_☉ satellite galaxy experiences the same dynamical friction whether the host dark matter halo is made of 10⁻⁵ eV axions, 100 GeV WIMPs, or 1 M_☉ primordial black holes.
3. Subsonic (v_M ≪ σ) vs. Supersonic (v_M ≫ σ) Velocity Regimes:
   • Low-Speed Stokes Limit (X ≪ 1, v_M ≪ σ): Taylor-expanding erf(X) gives [erf(X) − (2X/√π)e^−X2] ≈ (4 / [3√π]) X³. The v_M³ factors cancel, and dynamical friction behaves like viscous Stokes drag linearly proportional to velocity: dv_M/dt ≈ − (4√2 π / 3) (G² M ρ ln Λ / σ³) v_M.
   • High-Speed Limit (X ≫ 1, v_M ≫ σ): The bracket approaches 1, and the deceleration drops off inversely with the square of the speed: |dv_M/dt| ≈ (4π G² M ρ ln Λ) / v_M². Fast-moving objects barely deflect background stars as they pass, generating a weak wake and minimal drag—explaining why high-speed galaxy collisions in rich galaxy clusters (v_rel ∼ 1,500 km s⁻¹ ≫ σ_gal ∼ 250 km s⁻¹) pass through each other without merging, whereas slow galaxy group encounters (v_rel ∼ 200 km s⁻¹) merge via dynamical friction.

For a satellite galaxy, globular cluster, or black hole of mass M on a circular orbit at initial radius r_i inside a singular isothermal dark matter halo with flat circular velocity v_c = √2 σ and density profile ρ(r) = v_c² / (4π*G r*²), losing angular momentum L = M r v_c at rate dL/dt = r F_df causes the orbit to spiral to the center (r → 0) in the dynamical friction inspiral time (t_df) (Binney & Tremaine, 2008):

t_df = [ 1.17 / ln Λ ] ( r_i² v_c ) / ( G M ) ≈ [ 1.65 / ln Λ ] ( M_halo(r_i) / M ) t_dyn(r_i)

where t_dyn = r_i / v_c is the local dynamical crossing time.

In 1984, Scott Tremaine and Martin D. Weinberg (and Read et al., 2006) proved that Chandrasekhar's local formula breaks down when a massive satellite or globular cluster sinks into a constant-density harmonic core (ρ(r) = ρ₀ = constant, where gravitational potential Φ(r) = Φ₀ + ½ Ω₀² r² has a constant orbital frequency Ω₀ for all radii). Inside a harmonic core, background stars execute closed elliptical resonances synchronized with the orbiting perturber, causing dynamical friction to shut off ("core stalling"). This phenomenon provides a crucial cosmological test of dark matter profiles: in the Fornax dwarf spheroidal galaxy, five ancient globular clusters orbit at r ≈ 0.5–1.5 kpc despite having a classical Chandrasekhar inspiral time of only t_df ≈ 1–2 Gyr in a cuspy Navarro–Frenk–White (NFW, ρ ∝ r⁻¹) halo—providing dynamical evidence that Fornax possesses a constant-density dark matter core (r_core ≈ 1 kpc) that stalled their orbital decay (Goerdt et al., 2006; Read et al., 2006).

1. Chandrasekhar, S. (1943). Dynamical friction. I. General considerations: The coefficient of dynamical friction. The Astrophysical Journal, 97, 255–262. https://doi.org/10.1086/144517
2. Tremaine, S. D., Ostriker, J. P., & Spitzer, L., Jr. (1975). The formation of the nuclei of galaxies. I. M31. The Astrophysical Journal, 196, 407–411. https://doi.org/10.1086/153422
3. Mulder, W. A. (1983). Dynamical friction on high-speed bodies: The gravitational wake. Astronomy and Astrophysics, 117(1), 9–16. https://ui.adsabs.harvard.edu/abs/1983A%26A...117....9M
4. Tremaine, S., & Weinberg, M. D. (1984). Dynamical friction in spherical systems. Monthly Notices of the Royal Astronomical Society, 209(4), 729–757. https://doi.org/10.1093/mnras/209.4.729
5. Ostriker, E. C. (1999). Dynamical friction in a gaseous medium. The Astrophysical Journal, 513(1), 252–258. https://doi.org/10.1086/306858
6. Read, J. I., Goerdt, T., Moore, B., Pontzen, A. P., Stadel, J., & Lake, G. (2006). Dynamical friction in constant density cores: A failure of the Chandrasekhar formula. Monthly Notices of the Royal Astronomical Society, 373(4), 1451–1460. https://doi.org/10.1111/j.1365-2966.2006.11022.x
7. Conroy, C., Naidu, R. P., Garavito-Camargo, N., Besla, G., Kravtsov, A., Gnedin, O. Y., & Bonaca, A. (2021). All-sky dynamical response of the Galactic halo to the Large Magellanic Cloud. Nature, 592(7855), 534–536. https://doi.org/10.1038/s41586-021-03385-7