On the radial velocity wave in the Galactic disk

Abstract: Stars in the Galactic disk have mean radial velocities $\overline{v}R$ that oscillate as a function of angular momentum $J\varphi$. This `$J_\varphi$-${\overline{v}}R$ wave' signal also exhibits a systematic phase shift when stars are binned by their dynamical temperatures. However, the origin of the wave is unknown. Here we use linear perturbation theory to derive a simple analytic formula for the $J\varphi$-$\overline{v}R$ signal that depends on the equilibrium properties of the Galaxy and the history of recent perturbations to it. The formula naturally explains the phase shift, but also predicts that different classes of perturbation should drive $J\varphi$-$\overline{v}R$ signals with very different morphologies. Ignoring the self-gravity of disk fluctuations, it suggests that neither a distant tidal kick (e.g., from the Sgr dwarf) nor a rigidly-rotating Galactic bar can produce a qualitatively correct $J\varphi$-$\overline{v}R$ wave signal. However, short-lived spiral arms can, and by performing an MCMC fit we identify a spiral perturbation that drives a $J\varphi$-${\overline{v}}_R$ signal in reasonable agreement with the data. We verify the analytic formula with test particle simulations, finding it to be highly accurate when applied to dynamically cold stellar populations. More work is needed to deal with hotter orbits, and to incorporate the fluctuations' self-gravity and the role of interstellar gas.

• The paper presents a linear perturbation theory that accurately reproduces the radial velocity wave seen in Gaia data through analytic derivations and test particle simulations.
• It compares various perturbation models, demonstrating that transient spiral structure best explains the observed multi-modal v wave morphology compared to bars or impulsive kicks.
• Numerical simulations and MCMC fitting validate the model for dynamically cold populations, while highlighting limitations for stars with higher velocity dispersions.

The Origin and Modeling of the Radial Velocity Wave in the Galactic Disk

Introduction

The paper "On the radial velocity wave in the Galactic disk" (2602.06182) presents a comprehensive theoretical and numerical analysis of the oscillatory mean radial velocity ($\overline{v}_R$) signals observed in the Galactic disk as a function of stellar angular momentum ($J_\varphi$), a phenomenon referred to as the "v wave". These features, uncovered by Gaia, challenge the equilibrium assumption for the Milky Way disk and point to recent perturbative events. The work builds a detailed linear perturbation theory formalism, tests its validity through high-statistics test particle simulations, and systematically assesses which classes of perturbations (external kicks, bars, spirals) can explain the observed data morphology.

Observational Context and Motivation

The "v wave" is an oscillatory feature in $\overline{v}_R(J_\varphi)$ observed for stars near the Sun and on dynamically cold orbits. Initial studies highlighted multiple comparable-amplitude wave-like modes with wavelengths of several hundred $\,\mathrm{kpc\,km\,s^{-1}}$ in $J_\varphi$. These features deviate from equilibrium models and indicate the disk's response to recent dynamical perturbations.

Figure 1: The observed v signal (black, panel a) for cold disk stars in a narrow solar-azimuth wedge, compared against analytic and simulation-based model fits and an agnostic multi-sinusoid fit.

Analytic Theory of the v Wave

The authors develop an explicit linear perturbation theory in action-angle variables for stars in a nearly-epicyclic approximation, focusing on the $J_\varphi$-$\overline{v}_R$ relation at fixed azimuth. The mean radial velocity is calculated as the response of a Schwarzschild DF to an externally prescribed potential perturbation.

The analytic results show:

• For a perturbation with $m$-fold symmetry, the response features two $J_\varphi$-space waves with frequencies $$\omega^\mathrm{eff}_{m,\pm} = m(\Omega + 2\kappa' \langle J_R\rangle) \pm \kappa$$.
• The $J_\varphi$-wavelengths increase with $J_\varphi$ and decrease with time since the perturbation.
• Dynamically hotter populations (larger $\langle J_R \rangle$) produce v waves with a systematic phase lag relative to colder populations.
• The analytic expressions are independent of the detailed $J_\varphi$-dependence of the background DF.

Numerical Validation with Test Particle Simulations

The analytic formalism is confronted with test particle simulations in various limits, for stars drawn from a Schwarzschild DF, evolving in a static galaxy potential under different perturbations.

Figure 2: Radial velocity signal for dynamically cold populations under an impulsive, quadrupolar external perturbation.

Figure 3: $J_\varphi$-$\overline{v}_R$ signal following an impulsive kick, for various perturbation strengths and velocity dispersions, showing strong agreement between theory and simulations in the dynamically cold, linear regime.

Results:

• For $$\epsilon = (\langle J_R\rangle/J_\varphi)^{1/2}\ll 0.1$$ and weak perturbation amplitude ($|\eta| \lesssim 0.05$), the linear analytic predictions are accurate.
• Deviations at larger velocity dispersion arise from breakdown of the long-wavelength approximation, not linearity.

Comparison of Perturbation Classes

Distant Impulsive Kicks

Single, impulsive, quadrupolar external perturbations (e.g., from the Sagittarius dwarf) induce a response dominated by a single low-frequency wave component, as expected from quadrupole selection rules. This is inconsistent with the observed multi-component structure in the Galactic v wave.

Rigidly Rotating Bar

A model bar perturbation leads to a v signal concentrated near Lindblad resonances, dominated by a single or (at most) two wave components depending on the bar switch-on timescale. The analytic and numerical results show pronounced features at the outer Lindblad resonance but fail to reproduce the observed multi-mode structure or the overall pattern for $J_\varphi$ outside resonance.

Figure 4: v wave signal in response to a rigidly-rotating bar perturbation under various dynamical conditions.

Transient Spiral Structure

A Gaussian-windowed, logarithmic spiral potential perturbation, either rigidly rotating or shearing, is shown to generate a v signal with multiple significant Fourier components, matching the observed multi-modal wave morphology.

Figure 5: Parameter space exploration of v wave responses to transient rigidly rotating and shearing spirals, showing how pitch angle, radial width, and time dependence modulate the v signal structure.

Figure 6: Comparison of analytic theory and simulations for the v wave induced by transient spirals, highlighting excellent agreement for cold, linear regimes.

Data Fitting and Model Inference

The analytic spiral model is fit to the Gaia v wave data using MCMC sampling, optimizing over spiral parameters (amplitude, pattern speed, width, lifetime, pitch angle, phase, corotation radius, offset). The best-fit parameters correspond to a transient, two-armed spiral with:

• $R_{\mathrm{CR}} \sim 11$ kpc
• $\alpha \sim 40^\circ$
• Short lifetime $\tau \sim 11$ Myr
• Dimensionless amplitude $\eta \sim 0.075$

  Figure 7: Dimensionless overdensity $\delta\Sigma/\Sigma_0$ at $t=0$ for the best-fit spiral, illustrating the physical plausibility of the underlying density fluctuation required to match the data.

  

  Figure 8: Model fit and simulation comparison for higher velocity dispersion; although the theory matches the data for this $\epsilon$, it fails to reproduce the simulation, indicating the breakdown of the underlying analytic assumptions.

  

  Figure 9: Posterior distributions from the MCMC parameter fit, indicating the sharp constraints achievable from the v wave data given the analytic model.

Strong numerical results: The best-fit analytic model, despite having two fewer free parameters than a multi-sinusoid fit, matches the amplitudes and phases of the Gaia v wave to within a few km/s and $\sim$50 $\,\mathrm{kpc\,km\,s^{-1}}$ in $J_\varphi$. However, the inferred spiral lifetime is implausibly short compared to theoretical expectations for swing-amplified or global unstable modes, highlighting possible deficiencies in the model or unaccounted systematic effects.

Theoretical and Practical Implications

Model Validity and Limitations

• The linear and long-wavelength analytic model is quantitatively accurate only for dynamically cold orbits $$(\epsilon\lesssim0.1,\,\sigma_{R}\lesssim30\,\mathrm{km\,s^{-1}})$$. For hotter orbits, omitted higher-order corrections become significant.
• Neglect of the self-gravity of the perturbation is justified only if the external potential applied accurately incorporates the self-consistent response of the disk (e.g., domination by stable Landau modes).

Origin of the v Wave

The study decisively rejects the hypothesis that the v wave is due to either a single external impulsive event or the bar alone (without collective self-gravity), as both predict qualitatively incorrect v wave morphologies. Only transient spiral structure—possibly excited internally or by mild perturbations—is dynamically consistent with the observed data.

Systematics and Future Prospects

• Systematic uncertainties in distances, potential models, and selection bias can alter v wave amplitudes and phase.
• To model hotter populations and include collective effects, the analytic framework should be extended beyond the current long-wavelength and linear truncations, or replaced by numerical linear response or N-body methods.
• The framework enables efficient parameter inference and is amenable to future integration with advanced dynamical modeling tools for the Milky Way, including chemical and vertical phase-space structure.

Conclusion

This work provides a detailed analytic and numerical synthesis of the origins of the radial velocity wave in the Galactic disk. The main conclusion is that only transient spiral structure, modeled as a short-lived, non-rigid spiral perturbation, successfully reproduces the observed $J_\varphi$-$\overline{v}_R$ v wave signal in the Milky Way. While the analytic theory is highly predictive in the cold limit, extensions are required to robustly infer disk perturbation histories for dynamically hotter populations or to incorporate self-gravity. The results guide future development of perturbative and simulation-based dynamical inference methods for Galactic structure and time-dependent phenomena.