pith. sign in

arxiv: 2605.10092 · v2 · pith:GKB2N2V3new · submitted 2026-05-11 · 🌌 astro-ph.GA

Tracing the kinematic perturbations of the Milky Way spiral arms with APOGEE DR17 and Gaia DR3

Xi-Can Tang , Zhi Li , Iulia T. Simion , Hao Tian , Zhijian Luo , Shuting Fan , Zi-Qi Li This is my paper

Pith reviewed 2026-06-30 22:48 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords Milky Wayspiral armskinematicsradial velocitypitch anglepattern speedAPOGEEGaia
0
0 comments X

The pith

Revised Milky Way spiral model incorporating both sine and cosine radial velocity terms reproduces observed stellar streaming motions to the 2 percent level.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a revised steady-state response model for how a two-armed logarithmic spiral potential perturbs the radial velocities of stars in the galactic disk. The model adds the previously neglected cosine component to the sine term used in earlier work and accounts for Lindblad and corotation resonances when fitting the smoothed velocity field of RGB stars from APOGEE DR17 and Gaia DR3. If the approach is valid, it delivers a pitch angle of roughly 10 degrees, a local density contrast of 5 to 18 percent, and a pattern speed between 10 and 20 km s^{-1} kpc^{-1}. A reader would care because these parameters describe the dynamical influence of the Milky Way's spiral arms on stellar orbits and set the stage for more accurate maps of the disk's non-axisymmetric structure.

Core claim

The revised model that includes both the V_R,sin and the dynamically important V_R,cos components for a two-armed logarithmic spiral potential reproduces the phase and amplitude of mock radial-velocity fields to the ~2 percent level. When applied to the observational data while subtracting resonance effects, it yields a robust pitch angle p ≃ 10° and a local surface density contrast ξ ≃ 5--18 percent at the solar radius, with the spiral pattern speed constrained to Ω_p ≈ 10--20 km s^{-1} kpc^{-1}.

What carries the argument

Revised steady-state radial-velocity response model for a two-armed logarithmic spiral potential that incorporates both V_R,sin and V_R,cos components and subtracts Lindblad and corotation resonance contributions.

If this is right

  • The spiral pitch angle is determined to be approximately 10 degrees.
  • The local surface density contrast at the solar radius lies between 5 and 18 percent.
  • The spiral pattern speed is limited to the range 10-20 km s^{-1} kpc^{-1}.
  • Resonance locations strongly shape the observed velocity field, producing large amplitudes near Lindblad resonances and near-zero amplitudes near corotation.
  • The radial scale length of the perturbation remains only weakly constrained owing to parameter covariance.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same kinematic fitting procedure could be applied to external galaxies with resolved stellar velocities to compare spiral arm properties across different systems.
  • Discrepancies between the model and data in specific radial zones could be used to test whether transient arms or bar-driven perturbations contribute measurably.
  • Higher-precision velocity surveys could tighten the resonance modeling and reduce the covariance that currently limits the radial scale length.
  • The derived density contrast range provides a concrete target for hydrodynamic simulations that aim to reproduce Milky Way-like spiral structure.

Load-bearing premise

Stellar streaming motions are dominated by a steady-state two-armed logarithmic spiral potential whose resonances can be accurately modeled and subtracted, without dominant contributions from transient arms, the galactic bar, or other unmodeled perturbations.

What would settle it

A high-resolution map of the radial velocity field that shows velocities near the expected Lindblad resonances remaining modest rather than becoming extremely large, or a pattern speed derived from independent tracers falling outside the 10-20 km s^{-1} kpc^{-1} interval.

Figures

Figures reproduced from arXiv: 2605.10092 by Hao Tian, Iulia T. Simion, Shuting Fan, Xi-Can Tang, Zhijian Luo, Zhi Li, Zi-Qi Li.

Figure 1
Figure 1. Figure 1: Radial velocity and corresponding uncertainty maps of RGB stars in APOGEE DR17 and Gaia DR3, smoothed over a spatial scale of 2.5 times the bin size. The magenta star marks the position of the Sun. Here, A ′ G is a fixed orthogonal matrix that transforms equatorial coordinates to the Galactic frame, as defined in equation 4.62 of the Gaia EDR3 online documentation. The matrix A represents the local orthono… view at source ↗
Figure 2
Figure 2. Figure 2: Comparison between the AGAMA orbit-integration results and our VR models. Panels (a-d) show, respectively, the mock VR field from AGAMA, the corrected VR model, their residuals (δVR = VR,AG−VR VR ), and the VR,sin only model. Black spirals mark the loci of minimum gravitational potential [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Results from the VR model with identical parameters but varying pattern speeds from Ωp = 2 to 24 km s−1 kpc−1 in steps of 2 km s−1 kpc−1 . Green dashed circles mark the theoretical corotation resonance positions (RCR) for a flat rotation curve, while light-green dot-dashed circles indicate the Lindblad resonance positions (RILR and ROLR). Black spirals denote the loci of the minimum gravitational potential… view at source ↗
Figure 4
Figure 4. Figure 4: Schematic illustration of the relationship between the radial ve￾locity field and resonance locations. The blue shaded region represents the distribution of VR, vertical dashed lines indicate the resonance radii, and the horizontal dashed line marks VR = 0 km s−1 . The amplitude of VR is also modulated with the resonances, as illustrated schematically in [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Radial velocity distribution as a function of Galactocentric radius after spatial binning. Panel (a) displays the observational distribution, while Panels (b-d) present the spatially resampled results derived from orbit integrations with Ωp = 2, 12 and 17.77 km s−1 kpc−1 , respectively. The gray points represent the binned data points, the blue shaded regions denote the 16th and 84th percentile intervals. … view at source ↗
Figure 6
Figure 6. Figure 6: Spiral-arm loci traced by the RGB sample kinematics. The first and second rows present the results for fixed pattern speeds of Ωp = 2 and 12 km s−1 kpc−1 , respectively, while the third row shows the result for a free pattern speed with Ωp = 17.77 km s−1 kpc−1 . From left to right, the columns display the observational velocity maps, the best-fit model maps, the residual maps, and the stellar surface densi… view at source ↗
Figure 7
Figure 7. Figure 7: Comparison between the kinematically predicted spiral-arm features and previous studies. Panel (a) shows the three spiral-arm patterns derived in this work. Panels (b-d) compare our results with the spiral-arm loci from Eilers et al. (2020) (thin green dashed line), Reid et al. (2019) (thin blue dashed line), and Sun et al. (2024) (thin magenta dashed line), overlaid on the background map of mean radial ve… view at source ↗
Figure 8
Figure 8. Figure 8: Corresponding frequency curves and resonances in our model. The shaded region indicates the range of spiral pattern speeds inferred from our VR model. onance effects are explicitly taken into account. To this end, we perform a simple exploratory test. Since an extremely slow pat￾tern speed, such as Ωp = 2 km s−1 kpc−1 , is difficult to reconcile with most theoretical expectations, we restrict the tested ra… view at source ↗
read the original abstract

Aims. We constrain the dynamical perturbations of the spiral arms in the Milky Way disk, based on the non-axisymmetric streaming motions of RGB stars revealed by APOGEE and \textit{Gaia}. Methods. We develop a revised steady-state radial-velocity response model that incorporates both the \(V_{R,\sin}\) and the dynamically important \(V_{R,\cos}\) components for a two-armed logarithmic spiral potential. The model is validated using orbit integrations with \texttt{AGAMA} and Bayesian parameter recovery with \texttt{dynesty}, and is applied to the smoothed two-dimensional radial-velocity field of RGB stars while accounting for Lindblad and corotation resonances. Results. The revised model reproduces the phase and amplitude of the mock radial-velocity field to the \(\sim2\%\) level, substantially improving upon earlier \(V_{R,\sin}\)-only formulations. Applied to the observational data, it yields a robust pitch angle of \(p \simeq 10^\circ\) and a local surface density contrast of \(\xi \simeq 5\)--\(18\%\) at the solar radius. The radial scale length is less well-constrained (\(h_{R,1} \simeq 40\)--\(50\,\mathrm{kpc}\)) due to intrinsic parameter covariance. Resonance effects strongly shape the velocity field, thus affecting the fitting: the radial velocity becomes extremely large near the Lindblad resonances, whereas it vanishes close to the corotation resonance. Conclusions. Our results demonstrate that including both the \(V_{R,\sin}\) and \(V_{R,\cos}\) terms is essential for a physically consistent interpretation of stellar streaming motions induced by a spiral potential. The observed kinematics constrain the spiral pattern speed to \(\Omega_{p} \approx 10\)--\(20\,\mathrm{km\,s}^{-1}\mathrm{kpc}^{-1}\).

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper develops a revised steady-state radial-velocity response model for a two-armed logarithmic spiral potential that includes both the V_{R,sin} and V_{R,cos} components. The model is validated via AGAMA orbit integrations and dynesty Bayesian recovery on mocks generated from the same model (reproducing the mock radial-velocity field to ~2%), then applied to the smoothed 2D radial-velocity field of APOGEE DR17 + Gaia DR3 RGB stars while accounting for Lindblad and corotation resonances. It reports a pitch angle p ≃ 10°, local surface density contrast ξ ≃ 5--18%, radial scale length h_{R,1} ≃ 40--50 kpc (with noted covariance), and pattern speed Ω_p ≈ 10--20 km s^{-1} kpc^{-1}.

Significance. If the central modeling assumption holds, the work strengthens constraints on Milky Way spiral structure by showing that inclusion of the dynamically important V_{R,cos} term is required for physical consistency and yields improved mock recovery over prior V_{R,sin}-only formulations. The use of independent AGAMA integrations and dynesty recovery on mocks is a methodological strength that supports the internal validity of the revised model.

major comments (2)
  1. [Methods] Methods (model development and resonance handling): The revised model and its application to the observational data presuppose that the observed radial-velocity field is generated by a steady-state two-armed logarithmic spiral potential whose resonances can be accurately modeled and subtracted. No quantitative test or discussion is provided of whether the galactic bar or transient spiral arms could contribute comparably to the streaming motions; if they do, the recovered parameters (p ≃ 10°, ξ ≃ 5--18%, Ω_p ≈ 10--20) would not reflect true spiral properties. This assumption is load-bearing for the data-application step and the final constraints.
  2. [Results] Results (parameter constraints): The reported ranges for Ω_p and ξ are stated to be robust, yet the text notes that resonance effects strongly shape the velocity field (radial velocity becomes extremely large near Lindblad resonances and vanishes near corotation). No sensitivity analysis is shown for how uncertainties in resonance locations or unmodeled perturbations propagate into the posterior ranges obtained with dynesty.
minor comments (1)
  1. [Abstract] The abstract and Results section could more explicitly state the selection criteria and smoothing procedure applied to the APOGEE+Gaia RGB sample before fitting.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive report and recommendation for major revision. We address the two major comments point-by-point below, agreeing that additional discussion and analysis are warranted while defending the core scope of the work.

read point-by-point responses

  1. Referee: [Methods] Methods (model development and resonance handling): The revised model and its application to the observational data presuppose that the observed radial-velocity field is generated by a steady-state two-armed logarithmic spiral potential whose resonances can be accurately modeled and subtracted. No quantitative test or discussion is provided of whether the galactic bar or transient spiral arms could contribute comparably to the streaming motions; if they do, the recovered parameters (p ≃ 10°, ξ ≃ 5--18%, Ω_p ≈ 10--20) would not reflect true spiral properties. This assumption is load-bearing for the data-application step and the final constraints.

    Authors: The manuscript explicitly frames the work as constraining perturbations from a steady-state two-armed logarithmic spiral (see Aims and Methods sections), with validation performed exclusively on mocks generated from the same model. We agree that no quantitative test of bar or transient-arm contributions is included, as such a decomposition would require an entirely different modeling framework. We will add a new subsection in the Discussion explicitly stating this assumption, its load-bearing nature, and the caveat that the reported parameters apply under the spiral-only hypothesis. revision: yes

  2. Referee: [Results] Results (parameter constraints): The reported ranges for Ω_p and ξ are stated to be robust, yet the text notes that resonance effects strongly shape the velocity field (radial velocity becomes extremely large near Lindblad resonances and vanishes near corotation). No sensitivity analysis is shown for how uncertainties in resonance locations or unmodeled perturbations propagate into the posterior ranges obtained with dynesty.

    Authors: Resonance locations are already incorporated by masking regions near the Lindblad and corotation resonances during the dynesty fit, as described in the Methods. We acknowledge that an explicit sensitivity analysis (varying resonance radii within observational uncertainties and re-running the posterior sampling) is absent. In the revision we will add this analysis, including a table or figure showing the resulting variation in the Ω_p and ξ posteriors. revision: yes

standing simulated objections not resolved
  • Quantitative test or decomposition of whether the galactic bar or transient spiral arms contribute comparably to the observed streaming motions.

Circularity Check

0 steps flagged

No significant circularity; derivation uses independent validation tools and external survey data

full rationale

The paper formulates a revised radial-velocity response model incorporating both V_R,sin and V_R,cos terms for a two-armed logarithmic spiral, validates it via AGAMA orbit integrations and dynesty Bayesian recovery on mock data generated from the model itself (standard forward-model testing), then fits parameters to independent APOGEE DR17 and Gaia DR3 observations. No load-bearing steps reduce by construction to self-citations, fitted inputs renamed as predictions, or self-definitional relations; the central claims (pitch angle, density contrast, pattern speed) are outputs of fitting external data under stated assumptions rather than tautological re-derivations of inputs. The model is presented as newly revised in the manuscript with no cited prior self-work serving as the sole justification for uniqueness or ansatz.

Axiom & Free-Parameter Ledger

4 free parameters · 3 axioms · 0 invented entities

Central claim rests on a new model formulation whose parameters are fitted to data and on standard domain assumptions about the form of the spiral potential; no new entities are postulated.

free parameters (4)
  • pitch angle p = ~10°
    Fitted to observational radial velocity field
  • surface density contrast ξ = 5-18%
    Fitted amplitude parameter at solar radius
  • radial scale length h_R,1 = 40-50 kpc
    Fitted but poorly constrained due to covariance
  • pattern speed Ω_p = 10-20 km s^{-1} kpc^{-1}
    Constrained from resonance-affected velocity field
axioms (3)
  • domain assumption Steady-state assumption for the spiral potential
    Invoked to derive the radial-velocity response model
  • domain assumption Two-armed logarithmic spiral potential
    Basis for the revised V_R response formulation
  • domain assumption Resonance locations (Lindblad and corotation) can be accurately identified and modeled
    Used when applying the model to the smoothed velocity field

pith-pipeline@v0.9.1-grok · 5902 in / 1692 out tokens · 42814 ms · 2026-06-30T22:48:43.981685+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

2 extracted references · 1 canonical work pages

  1. [1]

    & Palicio, P

    Alves, J., Zucker, C., Goodman, A. A., et al. 2020, Nature, 578, 237 Antoja, T., Figueras, F., Romero-Gómez, M., et al. 2011, MNRAS, 418, 1423 Antoja, T., Ramos, P., López-Guitart, F., et al. 2022, A&A, 668, A61 Baba, J., Kawata, D., Matsunaga, N., Grand, R. J. J., & Hunt, J. A. S. 2018, ApJ, 853, L23 Barbillon, M., Recio-Blanco, A., de Laverny, P., & Pal...

    work page arXiv 2020

  2. [2]

    Surface density In a gravitational field, the surface density of a galaxy can be derived from the Poisson equation (see Binney & Tremaine 2008, Sec- tion 2.3)

    The corresponding radial velocity is then given by: VR,f ull =  2Ω0A R0(Ω0 −Ω p) + dA dR ! R0  m(Ω0 −Ω p) ∆ sinχ+ mA R0 tanp m(Ω0 −Ω p) ∆ cosχ,(B.7) B.2. Surface density In a gravitational field, the surface density of a galaxy can be derived from the Poisson equation (see Binney & Tremaine 2008, Sec- tion 2.3). Assuming that the 3D perturbed ...

    2008