Dynamical History of the Local Group in $Λ$CDM

Technical Summary: Dynamical History of the Local Group in ΛCDM

Problem Statement
The positions and radial velocities of galaxies within the Local Group (LG) encode the gravitational history of the system, dominated by the Milky Way (MW) and Andromeda (M31). While the "timing argument" has historically been used to estimate the total mass of the LG, previous applications often relied on simplified assumptions, such as spherical symmetry or the exclusion of external perturbers. Furthermore, recent high-precision observations reveal a significant tension: the observed radial velocities (RVs) of LG galaxies systematically exceed those predicted by standard $\Lambda$CDM models. Specifically, many non-satellite galaxies appear to be receding faster than a pure Hubble flow, a phenomenon difficult to reconcile with standard gravitational dynamics in a homogeneous universe. This paper investigates whether a more sophisticated dynamical model, incorporating axisymmetry and external tides, can resolve this discrepancy or if it points to a fundamental limitation of the $\Lambda$CDM paradigm for the Local Group.

Methodology
The authors employ a numerical approach to trace the dynamical history of the LG backward in time. The methodology involves the following key components:

• Simulation Framework: The study evolves a large number of test particles forward in time from an initial epoch ($z \approx 9$, scale factor $a_i = 0.1$) where they follow a pure Hubble flow with no peculiar velocities. The gravitational field is modeled using the equations of motion derived from General Relativity in the weak-field limit, incorporating the expansion of the universe (Hubble drag) and the gravitational potential of massive bodies.
• Massive Bodies: The primary gravitational sources are the MW and M31, treated as point masses (with modified potentials at small radii to mimic flat rotation curves) on a radial orbit. A third massive body, Centaurus A (Cen A), is explicitly included due to its proximity to the MW-M31 line of sight and significant mass. The model assumes axisymmetry, justified by the small tangential motion of M31 relative to the MW and the radial alignment of Cen A.
• Mass Accretion: Unlike static models, the MW and M31 are allowed to accrete mass during the simulation. Test particles passing within a defined "accretion radius" are considered captured, increasing the mass of the host galaxy.
• Parameter Space Exploration: The authors systematically vary the total initial mass of the MW+M31 ($M_i$) and the mass fraction residing in the MW ($q_1$). They also treat the Local Standard of Rest (LSR) speed ($v_{c,\odot}$) and an extra velocity dispersion term ($\sigma_{extra}$) as parameters.
• Observational Comparison: For each observed target galaxy (32 non-satellite galaxies), the simulation identifies a test particle trajectory that ends at the galaxy's current position. The model-predicted Galactocentric Radial Velocity (GRV) is projected onto the line of sight, corrected for solar motion, and compared to the observed Heliocentric Radial Velocity (HRV).
• Statistical Analysis: A $\chi^2$ analysis is performed to determine the posterior probability distributions of the model parameters. An additional variance term, $\sigma_{extra}$, is introduced to account for astrophysical noise (e.g., interactions between dwarf galaxies, tides from large-scale structures) not captured by the deterministic model.

Key Contributions

• Axisymmetric Modeling: The paper moves beyond the spherical symmetry assumption common in previous timing argument studies. By modeling the MW and M31 as a binary system on a radial orbit and including Cen A, the authors capture the non-spherical nature of the gravitational potential, which is crucial for galaxies at distances comparable to the MW-M31 separation (~0.8 Mpc).
• Explicit Inclusion of External Perturbers: The study directly incorporates the gravity of Centaurus A and estimates the tidal effects of M81, IC 342, and the Great Attractor. This allows for a more rigorous test of whether external tides can explain the observed velocity anomalies.
• Quantification of Discrepancy: The authors rigorously quantify the mismatch between $\Lambda$CDM predictions and observations, isolating a significant "extra" velocity dispersion that cannot be attributed to observational errors alone.

Results

• Mass Estimates: The best-fitting model yields a total initial LG mass of $M_{LG} = 4.33^{+0.37}_{-0.32} \times 10^{12} M_\odot$, with the MW contributing approximately $14 \pm 7\%$ of this mass ($q_1 \approx 0.14$).
• Systematic Velocity Discrepancy: Despite optimizing the mass parameters, the model fails to reproduce the observed radial velocities. The observed RVs systematically exceed the model predictions. The typical disagreement is $45.1^{+7.0}_{-5.7}$ km/s, with a maximum discrepancy of $110 \pm 13$ km/s for the galaxy DDO 99.
• Ineffectiveness of Standard Explanations:
  • External Tides: Including Cen A reduces the required $\sigma_{extra}$ by ~8 km/s, but the discrepancy remains large. Tides from M81, IC 342, and the Great Attractor were found to either have negligible effects or, in some cases, worsen the discrepancy by reducing predicted velocities further.
  • LMC Kinematics: Including the kinematic recoil of the MW due to the Large Magellanic Cloud (LMC) reduces $\sigma_{extra}$ by at most ~5 km/s, insufficient to resolve the tension.
  • Local Under-density/High $H_0$: Assuming a local under-density (empty rest of LG) or a higher local Hubble constant reduces $\sigma_{extra}$ by ~5–7 km/s, but not enough to match observations.
  • Dwarf Interactions: Interactions between LG dwarf galaxies are estimated to contribute only ~15 km/s to the velocity dispersion, far below the observed ~45 km/s.
• The "Flyby" Hypothesis: The authors note that the high velocities of distant galaxies (receding faster than the Hubble flow) could be explained by a past close flyby between the MW and M31, which would impart a gravitational slingshot impulse to passing material. However, such a high-speed encounter (~600 km/s relative speed) is not predicted by standard $\Lambda$CDM simulations, which show the MW and M31 have never approached closely. This scenario is, however, a natural outcome of Modified Newtonian Dynamics (MOND).

Significance and Claims
The paper concludes that no plausible set of initial conditions within the standard $\Lambda$CDM framework yields a good match to the radial velocities of the Local Group sample. The systematic excess of observed velocities over model predictions suggests that either:

1. There are significant unmodeled astrophysical processes (though the authors argue standard mechanisms like tides and dwarf interactions are insufficient).
2. The gravitational dynamics of the Local Group are not well-described by standard gravity in a $\Lambda$CDM universe.

The authors suggest that a past close flyby between the MW and M31, occurring roughly 9 Gyr ago, could plausibly explain the high velocities of non-satellite galaxies via gravitational slingshot encounters. While this scenario is difficult to reconcile with $\Lambda$CDM (where the galaxies have never interacted closely), it arises naturally in modified gravity theories like MOND. The paper posits that the observed kinematic anomalies, combined with the age of the MW's thick disc, may provide evidence for such a past interaction, challenging the standard cosmological model's description of the Local Group's dynamical history. The authors emphasize that their results highlight a tension that requires further investigation, potentially pointing toward modified gravity or a need for a more complex understanding of LG dynamics.