Galaxy Rotation Curves Explained

• Galaxy rotation curves are the radial profiles of circular velocity that reveal the mass distribution in galaxies, highlighting the discrepancy between luminous and dark matter.
• Observations using HI 21-cm and optical spectroscopy capture a steep central rise and a flat outer region, necessitating complex mass decomposition techniques.
• The near-universality of flat rotation curves and the Tully–Fisher relation underscore the critical role of dark matter or modified gravity in explaining galaxy dynamics.

Galaxy rotation curves are the radial profiles of circular velocity, $v(r)$, traced by stars, gas, and other kinematic tracers as a function of galactocentric radius. They serve as the principal astrophysical observables for constraining the radial mass distribution in galaxies—including contributions from baryons (stars, gas, stellar remnants) and other matter components not directly detectable through electromagnetic emission. A defining empirical result in spiral galaxy dynamics is the near-flatness of rotation curves beyond the optical disk, in conflict with the rapid Keplerian decline expected from the distribution of luminous matter alone. This discrepancy underpins the inference of dark matter halos, as well as motivates precision tests of gravity on kiloparsec scales.

1. Theoretical Underpinnings and Observational Techniques

Galaxy rotation curves are measured by mapping the velocity field of neutral hydrogen (HI 21 cm) or other emission lines (e.g., CO, optical nebular lines) using radio or optical spectroscopy. The observed rotational speed, $V(R)$, at galactocentric radius $R$ is related to the enclosed mass $M(<R)$ via the equilibrium condition for circular orbits in an axisymmetric potential: $V^2(R) = \frac{GM(<R)}{R}$ This equation can be inverted to determine the mass profile $M(<R)$—a direct dynamical probe of both baryonic and non-luminous matter.

Observationally, the shape of the rotation curve is characterized by a steep central rise (bulge/star-dominated), a maximum or plateau (disk-dominated region, typically peaking at $R \approx 2.2\,R_\mathrm{d}$ for an exponential disk), and a nearly constant velocity up to large radii (the “flat rotation curve” regime), often persisting well beyond the inferred edge of the stellar disk (Sofue, 2013). High-resolution measurements in the Milky Way combine stellar, gas, and maser observations to reconstruct a continuous grand rotation curve from sub-parsec to hundreds of kiloparsec scales (Sofue, 2013).

2. Mass Decomposition and Modeling Approaches

The dynamical mass inferred from rotation curves is customarily decomposed into one or more components, each associated with physically motivated mass profiles:

Component	Typical Density/Surface Profile	Rotation Curve Contribution
Bulge	de Vaucouleurs: $$\Sigma_{\mathrm{b}}(R) = \Sigma_{be}\exp(-\kappa [(R/R_b)^{1/4}-1])$$	$V_\mathrm{b}(R)$ via numerical integration
Disk	Exponential: $\Sigma_{\mathrm{d}}(R) = \Sigma_{dc}\exp(-R/R_d)$	Closed-form: involves Bessel functions
Halo	Isothermal, NFW: e.g. $\rho_\mathrm{NFW}= \rho_0/[(r/h)(1+r/h)^2]$ (Sofue, 2013)	$V_\mathrm{h}(R)$ via Poisson/integral solution
Central BH	Point mass $\sim 10^6-10^9\,M_\odot$	$V_\mathrm{BH}(r)\propto r^{-1/2}$

In one approach (“decomposition method”), the observed $V(R)$ is fitted with the sum in quadrature of model components: $$V^2(R) = V_b^2(R) + V_d^2(R) + V_h^2(R) + V_{\mathrm{BH}}^2(R)$$ with the functional form and parameters (mass-to-light ratios, scale lengths, central densities) optimized for best fit (Sofue, 2013).

Alternatively, the “direct method” infers the mass density from $V(R)$ without assuming a functional form: for a spherical case,

$$M(R) = \frac{RV^2(R)}{G};\qquad \rho(R) = \frac{1}{4\pi R^2} \frac{dM}{dR}$$

An analogous (though more complicated) formulation exists for a thin axisymmetric disk (Sofue, 2013).

3. Key Empirical Results and Universality

Compilation of rotation curves across disk galaxies (Sa–Scd) reveals several robust features:

• The inner parts are generally baryon-dominated, showing a rise in $V(R)$ consistent with the luminous mass traced by photometry under standard mass-to-light ratios.
• The plateau or gently rising region beyond the optical radius ($R_{25}$) is incompatible with a Keplerian falloff, indicating the presence of a massive, spatially extended component.
• The outer flatness is well described by the addition of a nearly isothermal, spherical (or slightly oblate) dark matter halo, with $\rho \sim r^{-2}$ at the relevant radii.

The universality of these features is strongly evidenced by the similarity of scaled rotation curve shapes across galaxies—even as the normalization ($V_\mathrm{max}$) varies (basis for the Tully–Fisher relation) (Sofue, 2013, Tiley et al., 2018).

4. Modified Gravity and Alternative Theories

The observed flatness of rotation curves has catalyzed extensive study of modifications to Newtonian dynamics and theories of gravity:

• Fourth Order Gravity (FOG): FOG introduces higher-order curvature invariants into the gravitational action, leading to a Newtonian potential with Yukawa-type corrections:

$$\Phi_\mathrm{FOG}(r) = -\frac{GM}{r} \left[1 + \frac{1}{3}e^{-\mu_1 r} - \frac{4}{3}e^{-\mu_2 r}\right]$$

The two mass scales $\mu_1$, $\mu_2$ correspond to new gravitational degrees of freedom associated with the Ricci scalar and Ricci/Riemann tensor terms. At intermediate radii, the Yukawa terms can enhance $V(r)$ over the GR prediction, but they decay exponentially at large $r$, returning to Keplerian behavior. Fits to the Milky Way and NGC 3198 show that, without a dark matter halo, FOG cannot sustain flat $V(r)$ at large radii; a DM component remains required for consistency with data (Stabile et al., 2011).

• Grumiller Gravity: Adds a linear “Rindler term” ($ar$) to the potential, leading to $v(r) = \sqrt{M/r + a r}$ and an effective flatness at large $r$ if $a$ is nonzero. The data for several galaxies are well fit with a universal $a \simeq 0.3 \times 10^{-10}$ m s$^{-2}$, about a quarter of MonD’s characteristic acceleration $a_0$, but deviations appear at larger radii where the predicted $V(r)$ increases above flat; some galaxies require modifications to the disk scale length or a radially varying $a$ for a satisfactory fit (Lin et al., 2012).
• Newtonian Baryonic Models: By systematically integrating Newton’s law over a disk with empirically motivated baryonic density profiles, flat to rising rotation curves are recovered for a broad class of mass distributions—even in the absence of a DM halo. For instance, with a density profile $\rho \sim 1/(x+1)^k$, nearly flat $V(r)$ appears when the peripheral density remains a few percent of the central (Pavlovich et al., 2014). Empirical fits to 47 galaxies yield velocity correlations $>0.995$ even without invoking DM or modified gravity.
• Relativistic Corrections: Special relativistic mass enhancement, when factored into the density distribution, produces only a minor increase (a small “dark-matter-like effect”) in $V(r)$ at a given radius. The discrepancy between disk-based and Keplerian curves remains dominated by geometry rather than relativity (Jaracz, 2023).

5. Constraints from the Milky Way and Extragalactic Systems

Milky Way rotation curves, compiled from kinematic tracers spanning parsec to $\sim$100 kpc, provide stringent multi-scale constraints. The innermost region ($r \lesssim$ few pc) is dominated by the central black hole, transitioning to bulge ($r \sim 10^2$ pc) and main disk ($r \sim 3$–10 kpc) before reaching the dark-matter-dominated regime at tens of kpc. Detailed analysis deconvolves the bulge into at least two components (inner core and main bulge) with exponential (not de Vaucouleurs) profiles—a Plummer-like exponential sphere model yielding a much sharper central peak in $V(r)$ than the standard $r^{1/4}$ law (Sofue, 2013).

Kinematic tracers such as classical Cepheids indicate that outer Milky Way rotation speeds may decline below the canonical “flat” profile, with the implied dark matter density at least 60% lower than that inferred from flat-curve fits (Gnaciński, 2018). However, non-circular motions, potential triaxiality, and observational uncertainties complicate strong conclusions in the outermost disk and halo.

6. The Tully–Fisher Relation and Its Theoretical Basis

A robust empirical relation links the total baryonic mass $M_\mathrm{bar}$ to the asymptotic (flat) rotation velocity, $V_f$, as $M_\mathrm{bar} \propto V_f^4$ (baryonic Tully–Fisher relation, BTFR). Both standard isothermal halo models and some modified gravity and metric theories naturally reproduce this scaling: $M = \left(\frac{1}{2\beta G c^2}\right) V^4$ where $\beta$ is a universal inverse length scale, identified in some metric approaches with the inverse size of the observable universe. This metric-based explanation links the origin of the BTFR to modifications of global spacetime structure, bypassing explicit dark matter (Bhattacharyay, 2020).

7. Open Issues, Systematic Effects, and the Diversity Problem

Careful modeling of rotation curves must confront non-equilibrium processes, non-circular motions, and additional stress contributions (magnetic, thermal, viscosity, cosmic rays). Analysis with state-of-the-art hydrodynamical simulations (FIRE suite) indicates that the “constructed” rotation curve (from observed gas motions under assumptions of equilibrium and axisymmetry) can underestimate or overestimate the true circular velocity by up to 50% in dwarfs subject to feedback or mergers (Sands et al., 2024). This effect can give rise to an “artificial rotation curve diversity,” mimicking true diversity in dark matter halo structure. Even for well-ordered disks, remaining systematics limit the accuracy to 10% (Sands et al., 2024).

In massive star-forming galaxies over the past $\sim$10 Gyr, rigorous stacking and normalization of rotation curves reveal that, when properly scaled (using stellar disk scale-length, not dynamical turnover), the average curve remains flat or slightly rising to six scale lengths (about 13 kpc) at all redshifts. Inferred dark matter fractions within $6R_d$ are consistently $\gtrsim$60%, with little evolution, consistent with expectations from $\Lambda$CDM cosmology and hydrodynamical simulations (e.g., EAGLE) (Tiley et al., 2018).

Alternate and modified gravity models, including those using higher-order actions, conformal metrics, or superfluid-vacuum-induced potentials, can reproduce a variety of rotation curve morphologies observed in the HI Nearby Galaxy Survey (THINGS). However, in all frameworks, a common theme emerges: for models to explain all observed rotation curves (including those with declining $V(r)$ in the outer disk, or very flat curves in low surface brightness systems) sometimes additional parameters or new physics are required, and the data generally support significant non-baryonic mass or modifications at large radii (Zloshchastiev, 2023, Wojnar et al., 2018, Marr, 2015).

Summary Table: Select Theoretical Rotation Curve Formulas

Model/Class	Formula for $\Phi$ or $v(r)$	Flatness at Large $r$	Parameters
GR (Newtonian)	$\Phi = -GM/r$; $v^2 = GM/r$	Declines ($\propto 1/\sqrt{r}$)	$G$, $M$
Fourth Order Gravity	$$\Phi = -\frac{GM}{r}[1+\frac{1}{3}e^{-\mu_1 r} - \frac{4}{3}e^{-\mu_2 r}]$$	Returns to Keplerian	$\mu_1$, $\mu_2$ (curvature invariants)
Grumiller’s Gravity	$\Phi = -GM/r + a r$; $v^2 = GM/r + a r$	Rises ($\sqrt{a r})$	$a$ (Rindler acceleration), $M$
Empirical (THINGS fit)	$v^2 = (GM/r)[1 + b(1 + r/R_0)]$	Flat as $r \to \infty$, $v \sim$ const	$b$, $R_0$
NFW Halo	Standard Poisson solution with NFW density	Flat ($v \sim$ const)	$c$, $M_\mathrm{vir}$
Superfluid Vacuum Theory	$v^2(R) \sim$ Newtonian $+$ log/linear/quad terms	Shape from dominant term	$b_0$, $k_1$, $k_2$, etc.
Metric w/ Global Term	$$v = \frac{c}{\sqrt{2}}/\sqrt{\beta r + \lambda / r}$$	Flat ($v$ const) at large $r$	$\beta \sim 1 / r_0$ (cosmological size)

Concluding Perspectives

Galaxy rotation curves provide a critical tool for mapping the mass distribution in disk galaxies and testing fundamental theories of gravitation and dark matter. The near-universality of rotation curve shapes, deviations from the Newtonian predictions at large radii, and the consistent emergence of the Tully–Fisher relation are robust empirical results. Model-dependent details—such as the efficacy of Yukawa corrections in higher-order gravity, the role of Rindler or logarithmic potential terms, and the direct modeling of baryonic distributions—can all capture aspects of the data, but current observations generally imply the need for non-luminous mass or significant modifications to gravitational law on galactic scales.

Nuanced interpretation of rotation curve data now requires careful accounting for kinematic systematics (especially for dwarfs), empirical diversity in rotation curve shapes, and theoretical considerations ranging from structure formation to the quantum properties of the vacuum. Further progress depends critically on high-resolution observations, detailed modeling of baryonic and non-baryonic contributions, and continued development of both standard and alternative gravitational theories.