REVIEW 3 major objections 4 minor 33 references
Tracing Hercules in Galactic azimuth with Gaia DR2
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Hercules' angular momentum decreases with azimuth, as expected for stars trapped at the co-rotation resonance of a slow, dynamically old Galactic bar.
desk verdict A short, honest letter proposing a genuinely new azimuthal diagnostic for Hercules, but the data-model comparison is visual and the young-bar OLR degeneracy means the conclusion is conditional, not unique. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the pendulum reduction of stellar motion near a bar resonance. In slow and fast action-angle variables, with slow angle $\theta_s = l\theta_R + m(\theta_\phi - \Omega_b t)$ and slow action $J_s = J_\phi/m$, the near-resonant Hamiltonian becomes a pendulum whose librating solutions are the trapped orbits. At the $m=2$ co-rotation resonance ($l=0$), $\theta_s = 2(\theta_\phi-\Omega_b t)$, so the angular momentum $J_\phi = 2J_s$ responds strongly to a change in azimuth toward the bar; at the outer Lindblad resonance ($l=1$), the radial angle varies faster than $J_s$, leaving $J_\phi$ nearly constant. The trapped distribution function is built by averaging the unperturbed quasi-isothermal distribution over the pendulum angle, the step that assumes the bar has been present long enough for orbits to phase-mix.
What would settle it
Simulate a self-consistent bar younger than about 2 Gyr with pattern speed near $55\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}^{-1}$ and check whether its outer-Lindblad ridge in the (azimuth, angular momentum) plane reproduces the observed $-8\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}^{-1}\,\mathrm{deg}^{-1}$ slope; if it does, the claimed uniqueness of the co-rotation interpretation fails. On the data side, extending the Hercules ridge over a much wider azimuth range and finding a clear deviation from the predicted monotonic decline would disfavour the co-rotation model.
Extended reading notes
Core claim
The central claim is that the key discriminant between different explanations of Hercules is geometric: resonant trapping zones leave different footprints in the plane of Galactocentric azimuth versus angular momentum measured at a fixed radius. For the $m=2$ co-rotation resonance of the bar, trapped orbits must change their angular momentum markedly with azimuth, producing a ridge with slope around $-8\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}^{-1}\,\mathrm{deg}^{-1}$ in the model; for the $m=2$ outer Lindblad resonance of a phase-mixed old bar, the same ridge is essentially flat. The authors show, using the second Gaia data release with Bayesian photometric distances, that the Hercules ridge in mean radial velocity versus (azimuth, angular momentum) has the negative slope predicted by co-rotation, over azimuths roughly from $-20^\circ$ to $+20^\circ$. They interpret this as evidence that Hercules is built from stars co-rotating with a large, dynamically old bar, and that the only remaining route to an outer-Lindblad origin requires a bar younger than about 2 Gyr, before the trapped orbits have phase-mixed, i.e. spread evenly around their orbits.
Load-bearing premise
The argument assumes the bar's gravitational perturbation has been present for more than about 2 Gyr, so trapped stars are spread evenly around their resonant orbits; if the bar is younger, the competing fast-bar explanation can mimic the observed azimuthal trend.
Editorial extensions
If this is right
- A successful co-rotation identification fixes the bar's pattern speed near $39\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{kpc}^{-1}$, placing the Sun just inside co-rotation and making the Milky Way's bar large and slow rather than small and fast.
- The (azimuth, angular momentum) slope becomes a standard diagnostic: ridges that tilt with azimuth are resonance-trapped structures, while flat ridges are consistent with phase-mixed outer-Lindblad or linear-deformation features.
- The argument demands a dynamically old bar, older than about 2 Gyr, so any independent evidence that the bar formed recently would undermine this particular interpretation.
- If the same modelling is applied to other ridges, the 'horn' at high positive azimuth should show a steep, high-order-resonance behaviour, giving a further prediction that future data can test.
Reading between the lines
- Because the diagnostic compares slopes rather than absolute zero-points, it should transfer directly to other Galactocentric radii and to action-space ridges, allowing the bar's resonant structure to be mapped across the disc once surveys cover more azimuths.
- A natural follow-up is to measure stellar ages inside the Hercules ridge: a phase-mixed co-rotation origin requires the trapped population to be old, whereas the young-bar loophole would predict a distinct age-dependent signature in the same ridge.
- The acknowledged young-bar loophole implies that the test's uniqueness is conditional on bar age; combining the azimuthal slope with independent bar-age and pattern-speed measurements could tighten the conclusion further, while a wider azimuth range would test whether the slope remains linear or curves under the influence of spiral arms and satellite encounters.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a dynamical test to distinguish between two competing explanations of the Hercules moving group in Gaia DR2: trapping at the bar's co-rotation resonance (CR) of a dynamically old, large bar (pattern speed ~39 km/s/kpc), versus a faster bar whose outer Lindblad resonance (OLR) lies near the Sun. Using a pendulum approximation for resonant orbits (Sect. 2), the authors argue heuristically and then demonstrate in a Portail et al. (2017) based model that trapped CR stars show a strong variation of angular momentum J_phi with azimuth, whereas phase-mixed OLR-trapped stars show almost constant J_phi. They predict a ridge slope of about -8 km/s/kpc/deg in the (azimuth, J_phi) plane for CR. In Gaia DR2 with StarHorse distances, selecting a 400 pc annulus around the Sun, they plot the mean radial velocity in (phi, J_phi) and claim that the Hercules ridge follows the CR slope, while the OLR prediction (zero slope) is disfavored. They explicitly acknowledge that an OLR origin could reproduce the trend if the bar is younger than ~2 Gyr, i.e. if the trapped distribution is not yet phase-mixed.
Significance. If the result is robust, it would be a valuable step toward resolving the long-standing degeneracy between CR and OLR explanations of Hercules and would support a slow, large bar with a corotation radius beyond the Sun. The analytic pendulum heuristic is elegant and transparent, and the prediction is derived from an independent dynamical model (Portail et al. 2017; M19) rather than fitted to the Hercules ridge, so the test is not circular. The paper is honest about its limitations, including the acknowledged young-bar OLR degeneracy. However, the central comparison between model and data is purely visual, with no error bars or statistical measure, and the diagnostic power is conditional on the bar being dynamically old, which is not independently established here. These issues, while not fatal, mean that the paper's central claim is not yet fully supported.
major comments (3)
- [Sect. 3, Fig. 2] The main observational evidence is a visual alignment of the Hercules ridge with an overplotted line of slope -8 km/s/kpc/deg, while the OLR zero-slope line is shown for contrast. There is no quantitative fit of the ridge slope, no error bars on the ridge location, and no significance statement that the zero-slope hypothesis is excluded. Because the entire Letter rests on this model-data comparison, the claim that 'the Hercules angular momentum changes significantly with azimuth' (Abstract) needs at least a simple quantitative assessment, such as fitting the ridge in (phi, J_phi) in both model and data and reporting the slope and its uncertainty from a bootstrap or similar procedure.
- [Sect. 3] Systematic uncertainties in the transformation from observed parallax/proper-motion/RVS quantities to Galactocentric (R, phi, v_R, J_phi) are not quantified. The assumed values of R0 = 8.2 kpc, the solar peculiar motion, and v_c = 233.1 km/s enter directly into J_phi and phi; StarHorse distance uncertainties and the selection function of the RVS sample can bias the mean v_R in each bin. The authors should either propagate these uncertainties into the ridge location or provide a robustness check (e.g., varying R0 and v_c within current uncertainties) to show that the observed slope is not an artifact of these choices.
- [Sect. 2.3 and Sect. 4] The paper's abstract describes the CR signature as 'a distinctive prediction of such a model,' but the paper itself notes that an OLR origin with a bar younger than ~2 Gyr can produce a significant azimuthal variation of angular momentum, citing Minchev et al. (2010) and Trick et al. (2019b). Since the age of the bar is not constrained in this work, the test is conditional rather than distinctive. The conclusions in Sect. 4 are appropriately hedged ('reinforce the case... in the case of a dynamically old bar'), but the framing in the Introduction and Abstract overstates the discriminating power. I recommend softening the 'distinctive prediction' language or explicitly stating in the Abstract that the test discriminates only under the assumption of a dynamically old bar.
minor comments (4)
- [Sect. 2.3 and Fig. 1 caption] There is an inconsistency in the velocity grid bin size: the text states 'bin-size Delta v = 5 km/s' while the Fig. 1 caption writes 'Delta v = 10 km/s'. Please correct one of them.
- [Sect. 2.2] Typo: 'whith' should be 'with' in the sentence '...will have evolved in magnitude by ~pi/8 = 22.5 deg, whith almost constant theta_s and J_s.'
- [Fig. 1 caption] The units of the slope are given inconsistently: the caption says 'slope of -8 km/s/deg' while the text says '-8 km/s/kpc/deg'. Since the horizontal axis is azimuth in degrees and the vertical axis is J_phi in km/s*kpc, the correct unit is km/s*kpc/deg (or km/s/deg if J_phi is already divided by R0). Please unify.
- [Sect. 3] The number of stars after quality cuts (6,350,087) and after the annulus cut (1,535,484) is given, but the initial RVS sample size (~7 million) is only mentioned in passing. A small clarification of the selection steps would help reproducibility.
Circularity Check
No circularity: the Hercules azimuth-slope prediction is computed from an independent bar model and compared to Gaia data, not fitted to the ridge.
full rationale
The model prediction of a negative ~ -8 km/s/kpc/deg slope for the Hercules ridge in (phi, J_phi) space is obtained by computing the trapped distribution function from pendulum perturbation theory in the Portail et al. (2017) bar potential, with pattern speed fixed at 39 km/s/kpc by independent bar/bulge, gas and stellar modelling; it is not fitted to the Gaia DR2 ridge. The Gaia ridge itself is constructed directly from Gaia DR2 astrometry and RVS velocities with StarHorse distances, using independent solar parameters, and the comparison is a qualitative overplot of a model-derived slope on the observed map. The paper explicitly flags the phase-mixed/old-bar assumption and concedes (Sect. 3, last paragraph; Sect. 4) that a bar younger than ~2 Gyr at the OLR could mimic the observed slope, citing Minchev et al. (2010) and Trick et al. (2019b); this limits uniqueness but is an acknowledged physical assumption, not a circular reduction. Self-citations to Monari et al. (2017a, 2019) transmit standard Hamiltonian perturbation methods and an externally constrained potential, not a result that already contains the target slope. Therefore no load-bearing step reduces to its own input by construction.
Assumptions & free parameters
free parameters (2)
- Bar pattern speed Omega_b =
39 km/s/kpc
- Solar azimuth relative to bar major axis =
28 degrees
assumptions (4)
- domain assumption Quasi-isothermal background distribution function f0(JR,Jphi) from Binney (2010) and M19 represents the unperturbed disk.
- domain assumption The epicyclic approximation gives reliable actions and angles for solar-neighborhood orbits.
- domain assumption Trapped orbits are phase-mixed along pendulum angles, so the bar perturbation is dynamically older than about 2 Gyr.
- standard math First-order resonant perturbation theory, reducing the resonance to a pendulum, is valid near the resonances considered.
Cite this review
Pith. "Pith review of Tracing Hercules in Galactic azimuth with Gaia DR2." pith.science (2026). https://pith.science/paper/IHGUDGNZ
@misc{pith2026190801318,
author = {Pith},
title = {Pith review of: Tracing Hercules in Galactic azimuth with Gaia DR2},
year = {2026},
howpublished = {\url{https://pith.science/paper/IHGUDGNZ}},
note = {Machine review of arXiv:1908.01318}
}
read the original abstract
The second data release of the Gaia mission has revealed, in stellar velocity and action space, multiple ridges, the exact origin of which is still debated. Recently, we demonstrated that a large Galactic bar with pattern speed 39 km/s/kpc does create most of the observed ridges. Among those ridges, the Hercules moving group would then be associated to orbits trapped at the co-rotation resonance of the bar. Here we show that a distinctive prediction of such a model is that the angular momentum of Hercules at the Sun's radius must significantly decrease with increasing Galactocentric azimuth, i.e. when getting closer to the major axis of the bar. We show that such a dependence of the angular momentum of trapped orbits on the azimuth would on the other hand not happen close to the outer Lindblad resonance of a faster bar, unless the orbital distribution is still far from phase-mixed, namely for a bar perturbation younger than ~ 2 Gyr. Using Gaia DR2 and Bayesian distances from the StarHorse code, and tracing the average Galactocentric radial velocity as a function of angular momentum and azimuth, we show that the Hercules angular momentum changes significantly with azimuth as expected for the co-rotation resonance of a dynamically old large bar.
Figures
Reference graph
Works this paper leans on
-
[1]
2019, arXiv e-prints, arXiv:1904.05721
Abuter, R., Amorim, A., Bauboeck, M., et al. 2019, arXiv e-prints, arXiv:1904.05721
arXiv 2019
-
[2]
Anders, F., Khalatyan, A., Chiappini, C., et al. 2019, A&A, 628, A94
work page 2019
- [3]
- [4]
- [5]
- [6]
-
[7]
P., Wegg, C., Gerhard, O., et al
Clarke, J. P., Wegg, C., Gerhard, O., et al. 2019, MNRAS
work page 2019
-
[8]
1998, AJ, 115, 2384
Dehnen, W. 1998, AJ, 115, 2384
1998
Show all 33 references
-
[9]
1999, ApJ, 524, L35
Dehnen, W. 1999, ApJ, 524, L35
1999
-
[10]
2000, AJ, 119, 800 D’Onghia, E
Dehnen, W. 2000, AJ, 119, 800 D’Onghia, E. & Aguerri, J. A. L. 2019, arXiv e-prints, arXiv:1907.08484
2000 arXiv
-
[11]
2005, A&A, 430, 165
Famaey, B., Jorissen, A., Luri, X., et al. 2005, A&A, 430, 165
2005
-
[12]
2019, MNRAS, 488, 3324
Fragkoudi, F., Katz, D., Trick, W., et al. 2019, MNRAS, 488, 3324
2019
- [13]
-
[14]
2001, A&A, 373, 511 Gaia Collaboration, Brown, A
Fux, R. 2001, A&A, 373, 511 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1
2001
-
[15]
Hunt, J. A. S. & Bovy, J. 2018, MNRAS, 477, 3945
2018
-
[16]
Hunt, J. A. S., Bub, M. W., Bovy, J., et al. 2019, arXiv e-prints [arXiv:1904.10968]
2019 arXiv
-
[17]
2018, A&A, 616, A11
Katz, D., Antoja, T., Romero-Gómez, M., et al. 2018, A&A, 616, A11
2018
-
[18]
Laporte, C. F. P., Minchev, I., Johnston, K. V ., & Gómez, F. A. 2019, MNRAS, 485, 3134
2019
-
[19]
2016, ApJ, 824, 13
Li, Z., Gerhard, O., Shen, J., Portail, M., & Wegg, C. 2016, ApJ, 824, 13
2016
- [20]
-
[21]
McMillan, P. J. 2017, MNRAS, 465, 76
2017
-
[22]
2010, MNRAS, 407, 2122
Minchev, I., Boily, C., Siebert, A., & Bienayme, O. 2010, MNRAS, 407, 2122
2010
-
[23]
Minchev, I., Nordhaus, J., & Quillen, A. C. 2007, ApJ, 664, L31
2007
-
[24]
2014, PhD thesis, Rijksuniversiteit Groningen
Monari, G. 2014, PhD thesis, Rijksuniversiteit Groningen
2014
-
[25]
2016, MNRAS, 457, 2569
Monari, G., Famaey, B., & Siebert, A. 2016, MNRAS, 457, 2569
2016
-
[26]
2019, A&A, 626, A41 (M19) Article number, page 4 of 5 G
Monari, G., Famaey, B., Siebert, A., Wegg, C., & Gerhard, O. 2019, A&A, 626, A41 (M19) Article number, page 4 of 5 G. Monari et al.: Tracing Hercules in Galactic azimuth with Gaia DR2
2019
-
[27]
Monari, G., Kawata, D., Hunt, J. A. S., & Famaey, B. 2017c, MNRAS, 466, L113 Pérez-Villegas, A., Portail, M., Wegg, C., & Gerhard, O. 2017, ApJ, 840, L2
2017
-
[28]
2017, MNRAS, 465, 1621 (P17)
Portail, M., Gerhard, O., Wegg, C., & Ness, M. 2017, MNRAS, 465, 1621 (P17)
2017
-
[29]
C., Dougherty, J., Bagley, M
Quillen, A. C., Dougherty, J., Bagley, M. B., Minchev, I., & Comparetta, J. 2011, MNRAS, 417, 762
2011
-
[30]
2018, A&A, 619, A72
Ramos, P., Antoja, T., & Figueras, F. 2018, A&A, 619, A72
2018
-
[31]
L., Smith, L., & Evans, N
Sanders, J. L., Smith, L., & Evans, N. W. 2019, arXiv e-prints, arXiv:1903.02009 Schönrich, R., Binney, J., & Dehnen, W. 2010, MNRAS, 403, 1829
2019 arXiv
-
[32]
C., Binney, J., & Magorrian, J
Sormani, M. C., Binney, J., & Magorrian, J. 2015, MNRAS, 454, 1818
2015
-
[33]
H., Fragkoudi, F., Hunt, J
Trick, W. H., Fragkoudi, F., Hunt, J. A. S., Mackereth, J. T., & White, S. D. M. 2019b, arXiv e-prints [arXiv:1906.04786] Article number, page 5 of 5
1906 arXiv
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.