Delayed phase mixing in the self-gravitating Galactic disc
Pith reviewed 2026-05-18 07:13 UTC · model grok-4.3
The pith
Self-gravity in the Milky Way disc delays phase spiral winding by about 300 million years after a perturbation.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the self-consistent N-body simulation of the MW-Sgr system, an initial coherent vertical oscillation lasts at least 300 Myr after pericentric passage before a clear phase spiral forms and begins winding, in contrast to test-particle models where winding starts right away. This self-gravitating response is qualitatively reproduced by Widrow's 2023 analytical shearing-box model and produces a lag of roughly 0.3 Gyr in the solar neighbourhood.
What carries the argument
The self-gravitating response of the disc that sustains a coherent vertical oscillation for hundreds of Myr before the phase spiral develops and winds.
If this is right
- The winding time measured from a phase spiral underestimates the true time since the perturbing event when self-gravity is active.
- Previous estimates of the Milky Way phase spiral age that ignored self-gravity are too young.
- The corrected age of 0.6-1.2 Gyr for the observed phase spiral aligns more closely with the Sagittarius pericentric passage.
- Analytical models that include self-gravity can be used to adjust future timing estimates.
Where Pith is reading between the lines
- The length of the delay may change with disc surface density or the mass and orbit of the perturber.
- Similar delays could affect phase-spiral dating in other disc galaxies that experience satellite impacts.
- Radial gradients in the disc might produce different lags at different distances from the Galactic centre.
Load-bearing premise
The self-gravitating delay measured in the N-body model of the MW-Sgr system applies to the real Milky Way disc near the Sun.
What would settle it
An N-body simulation with different initial conditions that produces immediate phase-spiral winding without a 300 Myr coherent-oscillation phase, or Gaia data showing that the observed phase spiral matches the Sagittarius impact time without any added delay.
read the original abstract
The Gaia phase spiral is considered to work as a dynamical clock for dating past perturbations, but some of the previous studies neglected the disc's self-gravity, potentially biasing estimates of the phase spiral's excitation time. We revisit the impact of self-gravitating effects on the evolution of vertical phase spirals and quantify the bias introduced in estimating their excitation time when such effects are ignored. We analysed a high-resolution, self-consistent $N$-body simulation of the MW-Sagittarius dwarf galaxy (Sgr) system, alongside four test particle simulations in potentials constructed from the $N$-body snapshots. In each case, we estimated the winding time of phase spirals by measuring the slope of the density contrast in the vertical angle-frequency space. In the test particle models, the phase spiral begins winding immediately after Sgr's pericentric passage, and the winding time closely tracks the true elapsed time since the Sgr impact. Adding the DM wake yields only a modest (< 100 Myr) reduction of the winding time relative to Sgr alone. By contrast, the self-consistent $N$-body simulation exhibits an initial, coherent vertical oscillation lasting $\gtrsim$ 300 Myr before a clear spiral forms, leading to systematic underestimation of excitation times. An analytical shearing-box model with self-gravity, developed by Widrow (2023), qualitatively reproduces this delay, supporting its origin in the disc's self-gravitating response. Assuming that self-gravity affects phase mixing in the MW to the same degree as the $N$-body model, the lag induced by self-gravity is estimated to be $\sim$ 0.3 Gyr in the solar neighbourhood. Accounting for this delay revises the inferred age of the MW's observed phase spiral to $\sim$0.6-1.2 Gyr, in better agreement with the Sgr's pericentric passage. (shortened for arXiv)
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that self-gravity in the Galactic disc induces a delay of ≳300 Myr in the onset of vertical phase-spiral winding after the Sagittarius pericentric passage. This is shown by comparing a self-consistent N-body simulation of the MW-Sgr system (where an initial coherent vertical oscillation precedes clear spiral formation) against test-particle integrations in snapshot potentials (where winding begins immediately) and against Widrow (2023) shearing-box analytics, which qualitatively recover the delay. The authors conclude that neglecting self-gravity leads to systematic underestimation of excitation times and, assuming the simulated lag applies to the real Milky Way, revise the age of the observed Gaia phase spiral to ~0.6-1.2 Gyr.
Significance. If the central extrapolation holds, the result would meaningfully revise dynamical clocks for Sagittarius-induced perturbations and highlight the importance of collective effects in phase mixing. The manuscript earns credit for its direct N-body versus test-particle comparison and for the qualitative match to an independent analytical model, which together isolate self-gravity without obvious circularity or free-parameter fitting.
major comments (1)
- [Abstract and Discussion] Abstract (final paragraph) and corresponding discussion: The headline revision of the phase-spiral age to 0.6-1.2 Gyr rests on the assumption that the ≳300 Myr coherent-oscillation delay seen in the single MW-Sgr N-body run occurs to the same degree in the real solar-neighbourhood disc. The manuscript provides no quantitative match of the simulated vertical frequency distribution, surface density, or velocity dispersion at R≈8 kpc to Gaia constraints; without such a match the lag duration cannot be assumed identical, weakening the direct application to observed MW data.
minor comments (1)
- [Methods] Methods section: The description of the N-body run lacks explicit values for particle number, gravitational softening, time-stepping, and the precise procedure used to fit the slope in vertical angle-frequency space; these details are needed for reproducibility even if they do not alter the qualitative conclusion.
Simulated Author's Rebuttal
We thank the referee for their constructive feedback and for recognizing the significance of our direct N-body versus test-particle comparison. We address the major comment below.
read point-by-point responses
-
Referee: The headline revision of the phase-spiral age to 0.6-1.2 Gyr rests on the assumption that the ≳300 Myr coherent-oscillation delay seen in the single MW-Sgr N-body run occurs to the same degree in the real solar-neighbourhood disc. The manuscript provides no quantitative match of the simulated vertical frequency distribution, surface density, or velocity dispersion at R≈8 kpc to Gaia constraints; without such a match the lag duration cannot be assumed identical, weakening the direct application to observed MW data.
Authors: We acknowledge that the manuscript does not include a detailed quantitative comparison of the simulated vertical frequency distribution, surface density, or velocity dispersion at R≈8 kpc against Gaia constraints. The N-body model was constructed as a representative Milky Way-like disc with a Sagittarius-like perturber, with parameters chosen to broadly capture the relevant dynamics rather than to fine-tune local observables. The central result—the ≳300 Myr delay—is isolated by the controlled comparison to test-particle integrations in the same snapshot potentials, and is further supported by the independent Widrow (2023) shearing-box analysis. We agree that the precise lag duration may depend on the exact local disc properties and that the extrapolation to the real Milky Way therefore rests on an assumption. We will revise the abstract and discussion to state this assumption and the lack of a detailed Gaia match more explicitly as a caveat, while retaining the main conclusion that neglecting self-gravity leads to underestimation of excitation times. This is a partial revision. revision: partial
Circularity Check
No significant circularity; derivation from independent simulations and external model
full rationale
The paper measures the ≳300 Myr delay directly from the evolution of vertical phase spirals in a self-consistent N-body simulation of the MW-Sgr system by tracking the slope of density contrast in vertical angle-frequency space, then contrasts this with test-particle runs in snapshot potentials and a qualitative match to the independent Widrow (2023) shearing-box model. The ~0.3 Gyr lag estimate for the real Milky Way is explicitly presented as an assumption rather than a derived necessity, with no fitted parameters renamed as predictions, no self-definitional loops, and no load-bearing self-citations or imported uniqueness theorems. The central claim therefore remains self-contained against the simulation outputs and external benchmark.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The N-body simulation correctly captures the self-gravitating dynamics of the galactic disc.
- domain assumption The analytical shearing-box model by Widrow (2023) applies qualitatively to the Milky Way context.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.