Dwarf Galaxy Dynamics

• Dwarf galaxy dynamics are defined by the interplay of rotation, dispersion, and tidal forces shaping low-mass systems.
• Advanced modeling techniques, including Jeans analysis and HI rotation curve decomposition, yield precise mass estimates.
• Interactions such as stochastic heating, tidal stripping, and feedback processes drive morphological changes and scaling relations.

Dwarf galaxies are low-mass systems (typical $M_\ast\sim10^5$–$10^9\,M_\odot$) whose dynamics provide unique laboratories for studying dark-matter physics, baryonic feedback, environmental processes, and the evolutionary pathways of galaxies. The interplay between internal and external mechanisms governing their dynamical structure and temporal evolution leads to a range of observed morphologies, from rotationally supported HI-rich irregulars to purely pressure-supported dwarf spheroidals (dSphs). Dwarf galaxy dynamics are tightly linked to cosmological context, feedback-driven structural transformations, and stochastic processes, and have driven the development and application of an array of advanced dynamical modeling techniques.

1. Fundamental Dynamical Principles and Classification

Dwarf galaxy dynamics are dictated by the balance between random motions (velocity dispersion $\sigma$) and coherent rotation ($V_\mathrm{rot}$), the mass distribution in baryons and dark matter (DM), the degree of environmental perturbation, and the impact of stochastic heating mechanisms. Dynamical structure is commonly classified as:

• Rotation-supported systems: High angular-momentum, gas-rich disks (dIrrs, many BCDs) with $V/\sigma\gtrsim2$ in their HI or young stellar components. Rotation curves can be well measured from HI 21cm data and corrected for pressure support and asymmetric drift (Lelli, 2022, Lelli et al., 2012, Lelli et al., 2011).
• Dispersion-supported systems: Classical dSphs and extreme star-forming dwarfs at high redshift generally show $V/\sigma\lesssim1$ in their stellar populations, with flat or weakly-rising velocity-dispersion profiles indicating dominant random motions (Battaglia et al., 2022, Lokas, 2019, Ceverino, 2018, Kirby et al., 2012).
• Transition/compound systems: Many dwarfs in the isolated field (e.g., WLM, VV124) exhibit both rotation and substantial pressure support, often with spatially distinct kinematic subcomponents (Leaman et al., 2012, Kirby et al., 2012, Kolhe et al., 2 Feb 2026).

Morphologically, systems span disk-like, thick oblate, and spheroidal configurations; structural parameters such as the intrinsic flattening, scale length, and half-light radius $r_{1/2}$ or $r_{\rm half}$ are tightly coupled to dynamical state.

2. Dynamical Modeling Techniques and Mass Inference

The modeling of dwarf galaxy dynamics relies on kinematic tracers (usually RGB stars, HI gas), structural measurements, and solution of the steady-state Jeans equations (spherical or axisymmetric). Canonical dynamical frameworks include:

• Jeans Equation Modeling: Utilizes observed velocity-dispersion profiles, tracer number density $\nu(r)$, and models for anisotropy parameter $\beta$ to constrain the enclosed mass $10^9\,M_\odot$0 (Battaglia et al., 2022, Lokas, 2019). In the spherical case,

$10^9\,M_\odot$1

is solved and projected to compare with the observed line-of-sight velocity dispersion profile $10^9\,M_\odot$2. Both NFW and cored dark-matter profiles are commonly tested, as is the impact of orbital anisotropy.

• Orbit-superposition and DF approaches: Schwarzschild’s method or direct DF modeling allow for flexible, non-parametric anisotropy and multiple kinematic subpopulations (Lokas, 2019). Action-based or made-to-measure schemes are also increasingly used.
• HI Rotation Curve Decomposition: For gas-rich dwarfs, high-resolution HI kinematics are used to fit tilted-ring models and decompose the inferred circular velocity into stellar, gaseous, and DM contributions (Lelli, 2022, Lelli et al., 2012, Lelli et al., 2011).
• Mass Estimators: For dispersion-supported systems, robust mass estimators at the 3D half-light radius are widely adopted:

$10^9\,M_\odot$3

(Wolf et al. 2010; (Battaglia et al., 2022, Kirby et al., 2014)).

Modeling is subject to degeneracies—notably the mass–anisotropy degeneracy—which limits unique inference of the DM profile unless higher-order kinematic moments or multiple tracer populations are available (Lokas, 2019, Battaglia et al., 2022). Binary orbital motion, if unrecognized, can induce radius-dependent biases on inferred masses, especially in low-dispersion systems and in the outskirts, requiring careful multi-epoch velocity measurements (Wang et al., 2023).

3. Stochastic Heating, Tidal Effects, and the Dynamical Attractor

Recent $10^9\,M_\odot$4-body experiments demonstrate that both subhalo-driven stochastic heating and tidal stripping establish a robust attractor for dSph dynamics, regardless of their formation scenario (Peñarrubia et al., 27 Feb 2026). The key phenomenology can be summarized:

• Heating by dark subhaloes: Encounters with subhaloes induce stochastic energy diffusion, causing the orbits of stars to gain energy irreversibly, expand radially, and evolve toward a state where the half-light radius approaches the radius of maximum circular velocity of the host halo, $10^9\,M_\odot$5, and the velocity dispersion to $10^9\,M_\odot$6.
• Role of tides: For satellites orbiting galaxies like the Milky Way, the tidal field strips DM, reducing $10^9\,M_\odot$7 and $10^9\,M_\odot$8, and preferentially removes super-virial stars, thus accelerating convergence to the attractor.
• Empirical consequences: Under this “heating argument,” nearly all MW dSphs fall on a narrow sequence in $10^9\,M_\odot$9 vs $\sigma$0, with small satellites ($\sigma$1 kpc) following the tidal track of cuspy haloes and only the most extended systems showing evidence for cored DM profiles. The observed mass–luminosity relation is consistent in slope with abundance matching, but with peak halo masses offset low by $\sigma$21.5 dex, tracing mass loss via tides.
• Prediction: Isolated, early-quenched dwarfs should be systematically more extended at fixed velocity dispersion or luminosity than satellites, providing a test of the universality of the heating attractor with upcoming wide-field surveys such as LSST (Peñarrubia et al., 27 Feb 2026).

This attractor framework unifies the interpretation of dSph scaling relations and highlights the dominance of internal heating and tides over formation conditions in setting structural diversity.

4. Scaling Relations, Baryon–Dark Matter Coupling, and Formation Pathways

Extensive kinematic and structural surveys have delineated a tight set of scaling relations for dwarfs:

• Velocity dispersion–radius locus: Both isolated and satellite dwarfs populate a similar relation between $\sigma$3 and $\sigma$4, with little distinction due to environment, suggesting universality (Kirby et al., 2014, Battaglia et al., 2022). This challenges scenarios in which tidal or ram-pressure processes uniquely create dSphs and constrains solutions to the “too-big-to-fail” problem.
• Inner circular velocity gradient (d$\sigma$5): The gradient measured at one disk scale length, $\sigma$6, is the primary dynamical parameter correlating tightly with central surface brightness, gas surface density, and star-formation intensity. BCDs are distinguished by steep inner gradients and high central $\sigma$7, and form a continuous sequence with compact irregulars and rotating spheroidals as star formation and gas content decline (Lelli et al., 2013, Lelli et al., 2012, Lelli et al., 2011).
• Baryonic Tully–Fisher (BTFR) and radial acceleration relations (RAR): Star-forming dwarfs extend the BTFR and RAR of spirals to lower masses and accelerations, with remarkably low intrinsic scatter (Lelli, 2022). This enforces a “fine-tuning” or tight coupling between baryons and DM on sub-kpc scales, the dynamical origin of which remains under investigation.
• Morphological evolution: Simulations and observations indicate a plausible pathway: BCDs arise from high gas-concentration irregulars undergoing centralized starbursts; after the burst, these fade to compact irregulars, and, upon gas removal (e.g., ram pressure, harassment), transition into rotating spheroidals (Lelli et al., 2013, Lelli et al., 2012). In more massive systems or through major mergers, tidal stirring or direct dwarf–dwarf merging can create spheroids with prolate rotation (Lokas, 2019).

5. Environmental, Cosmological, and Feedback Processes

The dynamics of dwarfs in different environments and their cosmological context are critical:

• Cluster dwarfs and group infall: Infalling dwarf groups into clusters retain boosted velocity dispersions and peculiar velocity signatures for several Gyr, preserving dynamical memory due to negligible dynamical friction, and can be identified via LOS velocity–radius phase-space structure despite spatial mixing (Vijayaraghavan et al., 2014). Early-type dwarfs (dEs/dSphs) in clusters are thus sensitive probes of accretion and assembly histories.
• Ram-pressure and gas removal: Isolated dwarfs (e.g., WLM) can undergo measurable dynamical perturbations due to ram pressure from the intergalactic medium, producing rotation curve asymmetries and significant biases (up to 20–30%) in dynamical mass estimates if not modeled explicitly (Kolhe et al., 2 Feb 2026).
• Cosmological context: Simulations at different epochs (FirstLight, VELA, AGORA) show that dwarfs begin as turbulent, dispersion-dominated systems at high z, become dynamically “hot” and elongated at $\sigma$8, and only later ($\sigma$9) settle into cold, rotation-supported discs as gas accretion rates drop (Ceverino, 2018).
• Feedback-driven structural transformation: Stellar feedback can create central DM “cores” in sufficiently massive systems with extended SFHs, while the lowest-mass, ancient systems retain DM “cusps.” The mass–anisotropy degeneracy and lack of observables beyond the half-light radius constrain firm conclusions for individual objects (Battaglia et al., 2022, Lokas, 2019).

6. Theoretical Alternatives and Outstanding Issues

Novel dynamical signatures and alternative gravity theories are being investigated:

• Fuzzy dark matter (FDM): Tidal perturbations of FDM solitons in dSphs induce long-lived, coherent breathing modes in the stellar and DM densities, imprinting radial oscillations in the velocity-dispersion of tracers detectable at the $V_\mathrm{rot}$0–$V_\mathrm{rot}$1 level, sharply contrasting with CDM which phase-mixes rapidly. Upcoming astrometric missions may test the presence of such oscillations (Widmark et al., 2023).
• Modified gravity theories: $V_\mathrm{rot}$2 gravity models can fit dSph velocity-dispersion profiles without DM by tuning Yukawa-type parameters, but produce a bimodal parameter distribution incompatible with a universal modification of gravity, challenging their viability on dwarf-galaxy scales (Martino et al., 2022).
• Host tidal field dominance: In the Milky Way, the gravitational field of the host can set the observed $V_\mathrm{rot}$3 via impulsive heating, particularly for first-infall or disrupting satellites, leading to overestimated DM fractions if not accounted for (Hammer et al., 2018).
• Binary star contamination: Unrecognized binaries can systematically bias dynamical mass and density inferences, with the greatest effect in ultra-faint dwarfs; multi-epoch high-precision spectroscopy is necessary to mitigate these biases (Wang et al., 2023).

7. Outlook and Future Prospects

Ongoing and future wide-field surveys and simulations are poised to clarify key aspects of dwarf galaxy dynamics:

• Large HI surveys (ASKAP, MeerKAT, SKA) and multi-object spectroscopic campaigns will yield resolved kinematics for thousands of dwarfs, testing dynamical laws to unprecedented precision and expanding the low-mass baseline for scaling relations (Lelli, 2022).
• Deep imaging (LSST, Euclid) and tailored N-body/hydrodynamic simulations are increasingly able to disentangle the signatures of minor mergers, tidal infall, and early gas/stellar feedback at faint surface brightness; these will be crucial for probing the dynamical attractor paradigm and the signatures of alternative DM models (Pascale et al., 2024, Peñarrubia et al., 27 Feb 2026).
• Advanced modeling tools, including orbit-superposition, action-based DFs, and full 3D non-spherical Jeans solvers, promise more robust estimation of DM profiles, anisotropy, and the role of triaxiality.

Dwarf galaxy dynamics thus remain at the forefront of efforts to understand small-scale structure, DM particle physics, baryon–DM coupling, and the cosmic history of galaxy evolution, with robust observational and theoretical tests continuing to emerge from the synergy between detailed modeling, cosmological context, and multi-wavelength surveys.