REVIEW 2 major objections 4 minor 13 references
Kinematics of Multiple Stellar Populations in Globular Clusters with Gaia
T0 review · 2 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Two chemically distinct stellar populations in the globular cluster NGC 5904 (M5) rotate on different curves, with a phase shift in their sky-plane motions significant at the 3-sigma level.
desk verdict A thin but honest proceedings summary of work already in Cordoni+19; the M5 differential rotation claim is the headline, but the paper doesn't rule out a radial-segregation artifact. 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 central tool is the pseudo-color index $C_{U,B,I} = (U-B)-(B-I)$, which separates first- and second-generation stars along the red-giant branch after correcting for differential reddening. The kinematics come from Gaia DR2 proper motions for the same stars: the rotation curve is the median proper motion in angular bins around the cluster, and the internal kinematics use median radial and tangential components in circular annuli to build $\sigma_\mathrm{RAD}$, $\sigma_\mathrm{TAN}$, and the anisotropy profile $\sigma_\mathrm{TAN}/\sigma_\mathrm{RAD} - 1$. This combination lets the paper track the two populations' sky-plane motion over the full cluster, beyond the small field of view of earlier HST-based studies.
What would settle it
Recompute the NGC 5904 rotation curves after classifying first- and second-generation stars from medium-resolution spectroscopy rather than the pseudo-color index; if the 3-$\sigma$ phase shift in $\mu_\delta$ disappears under an independent classification, the claimed rotation difference is an artifact of the photometric separation.
Extended reading notes
Core claim
The central discovery is that the two chemically distinct populations of NGC 5904 (M5) do not share the same present-day rotation: their median proper-motion rotation curves are offset in phase, and the offset in the $\mu_\delta$ component is significant at the $3\sigma$ level. In the same cluster, first-generation stars show higher radial proper motions than second-generation stars between roughly 2 and 5 half-light radii. Across the rest of the sample, only NGC 0104 (47 Tuc) shows measurable plane-of-sky rotation, with the two populations rotating together; NGC 0104 also shows the second-generation population to be more radially anisotropic than the first. These differences are read as present-day traces of the different initial spatial configurations of the two generations, as expected if second-generation stars formed from gas segregated in the cluster center.
Load-bearing premise
The whole comparison rests on the assumption that the pseudo-color index $C_{U,B,I}$, after correcting for differential reddening, cleanly separates first- and second-generation stars along the red-giant branch; if photometric errors, residual reddening, or field-star contamination blur that separation, the rotation curves and anisotropy differences could be biased or spurious.
Editorial extensions
If this is right
- If the M5 phase shift is real, the two populations carry different angular momentum today, so formation models must produce or preserve a kinematic difference between first- and second-generation stars at the level of a few km/s.
- The lack of detectable rotation in the other studied clusters sets upper limits on present-day rotation that any model of globular cluster formation must not exceed.
- Confirmation that second-generation stars in NGC 0104 are more radially anisotropic than first-generation stars strengthens self-enrichment scenarios in which second-generation stars formed from a more centrally concentrated gas reservoir.
- The small velocity differences in the other clusters, under about 1 km/s, are consistent with two-body relaxation having erased primordial kinematic differences, meaning non-detections do not rule out initial differences.
- Wide-field Gaia proper motions can now expose population-level kinematics over the full cluster, making this a reusable diagnostic for the multiple-population phenomenon.
Reading between the lines
- A testable extension the paper does not perform: sort a larger sample of globular clusters by relaxation time and check whether first- versus second-generation kinematic differences shrink in dynamically older clusters; if they do, relaxation, not formation, sets the present-day signal.
- If the M5 phase shift is a real angular-momentum offset, the two populations should also show a detectable difference in line-of-sight radial velocities in the same radial range, which future spectroscopic surveys could confirm independently of Gaia astrometry.
- The same pseudo-color plus proper-motion pipeline could be applied to clusters of different mass and age to map how rotation differences correlate with dynamical age, turning the signal into a formation clock for the multiple-population phenomenon.
- The hints of different rotation in NGC 0104's outer regions, if confirmed with more stars, would suggest that any initial kinematic segregation is best preserved at large radius where relaxation is slowest.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings paper reports an extension of the authors' Gaia DR2-based kinematic analysis of multiple stellar populations in globular clusters, following Milone et al. (2018) and Cordoni et al. (2019). Using the pseudo-color C_U,B,I = (U-B)-(B-I) to separate first-generation and second-generation RGB stars and Gaia DR2 proper motions for membership and kinematics, it presents rotation curves on the plane of the sky and radial profiles of radial/tangential velocity dispersion and anisotropy for NGC 0104 and NGC 5904. The main claim is that, among the analyzed clusters, only NGC 0104 and NGC 5904 show significant rotation on the plane of the sky, and that the two stellar populations in NGC 5904 exhibit different rotation curves, most notably a phase shift in the mu_delta component that is reported as significant at the 3 sigma level. The other clusters are reported as non-rotating and mostly kinematically similar between the two populations.
Significance. If the result holds, this is a valuable contribution because it uses Gaia DR2 to extend the kinematic comparison of multiple stellar populations to large cluster radii, where HST fields cannot reach. The key finding---a present-day difference in the rotation curves of 1G and 2G stars in NGC 5904---would provide a new observational constraint on formation and self-enrichment scenarios. The empirical nature of the measurement, with no indication that the claimed effects are reduced by construction from fitted parameters, is a strength. The paper also adds to the sparse set of clusters with measured velocity-dispersion and anisotropy differences between populations. However, the proceedings format leaves out critical detail, and one potential confound, the radial segregation of 2G relative to 1G stars, is not addressed in the rotation-curve analysis.
major comments (2)
- [Section 2.1, Figure 1] The central claim that NGC 5904's 1G and 2G stars have different rotation curves rests on comparing median proper motions in angular regions without any stated control for the known radial segregation of the two populations. Since the self-enrichment scenarios cited in Section 1 predict that 2G stars are more centrally concentrated, and since rotation amplitude or axis can vary with radius, azimuthal medians over angular bins that contain different radial distributions can produce an apparent phase shift between the 1G and 2G curves even when the intrinsic rotation is identical. Please add a radial-matching control (for example, compute the rotation curves in matched radial shells, or subsample one population to match the radial distribution of the other) or quantitatively argue that M5's rotation is independent of radius within the sampled region. This test is needed before the 3 sigma phase shift can be interpreted as a kinematic difference between populations.
- [Section 2.1 and Figure 2] The paper reports a 3 sigma phase shift in mu_delta for NGC 5904 and Anderson-Darling p-values in Figure 2, but it does not state the sample sizes per population and per bin, the membership selection thresholds, the proper-motion uncertainty model, or how the significance and confidence intervals are computed. Without these details, the significance claims cannot be checked. Please report the number of stars in each angular and radial bin, the error bars or bootstrap/Monte Carlo estimates for the median curves, and the exact statistical test used for the phase-shift significance and for the dispersion-profile comparisons.
minor comments (4)
- [Abstract and Section 1] The cluster sample is inconsistent between the abstract and the body: the abstract lists NGC 0288, NGC 5904, NGC 6121, NGC 6752, and NGC 6838 as the new clusters, with NGC 0104 explored further, while Section 1 names six clusters including NGC 6254, and Section 2.1 later includes NGC 6254 as 'NCG 6254'. Please reconcile the list and correct the typo.
- [Figure 1] The rotation curves would be more informative with error bars or shaded confidence regions on the median proper-motion curves; currently the 3 sigma claim cannot be visually assessed from the figure.
- [Section 2.1] The statement that there is no evidence for rotation in the remaining clusters should be accompanied by quantitative upper limits (for example, the maximum rotation amplitude consistent with the data) rather than only a visual inspection of the median curves.
- [Section 2.2] Please clarify whether the velocity-dispersion profiles are corrected for Gaia DR2 measurement errors and astrometric systematics, since these corrections can be important at large radii where proper-motion uncertainties grow.
Circularity Check
No significant circularity: the central rotation and internal-kinematics results are independent empirical measurements using Gaia DR2 proper motions and photometric population labels.
full rationale
The paper's central claims—that only NGC0104 and NGC5904 show significant rotation on the plane of the sky and that the two populations in NGC5904 have different rotation curves—are direct observational measurements. The 1G/2G labels are assigned via the pseudo-color index CU,B,I = (U−B) − (B−I), which is photometric, while the kinematics are measured from Gaia DR2 proper motions, an independent observable. No equation in the paper defines the rotation curves in terms of the population labels or vice versa, and no fitted parameter is later renamed as a prediction. The references to Milone et al. (2018) and Cordoni et al. (2019) are self-citations used to situate the method and the ongoing project, but the load-bearing result for NGC5904 is a new empirical finding, not a consequence of those papers' assumptions. Potential concerns about radial segregation or sample matching are validity risks rather than circularity: the text does not show that the claimed 3-sigma phase shift is forced by the construction of the sample or by the definition of any derived quantity. The manuscript is therefore self-contained as a measurement report, and no circular step can be exhibited from the quoted text.
Assumptions & free parameters
assumptions (3)
- domain assumption The pseudo-color index CU,B,I = (U−B) − (B−I) effectively separates first-generation (1G) and second-generation (2G) RGB stars in each cluster.
- domain assumption Gaia DR2 proper motions, parallaxes, and astrometric quality parameters are sufficient to define cluster membership and to measure internal kinematics without significant spatially correlated systematics.
- domain assumption The differential-reddening correction is accurate and does not introduce population-dependent spatial patterns.
Cite this review
Pith. "Pith review of Kinematics of Multiple Stellar Populations in Globular Clusters with Gaia." pith.science (2026). https://pith.science/paper/FSUNTWUV
@misc{pith2026190811692,
author = {Pith},
title = {Pith review of: Kinematics of Multiple Stellar Populations in Globular Clusters with Gaia},
year = {2026},
howpublished = {\url{https://pith.science/paper/FSUNTWUV}},
note = {Machine review of arXiv:1908.11692}
}
read the original abstract
The internal dynamics of multiple stellar populations in Globular Clusters (GCs) provides unique constraints on the physical processes responsible for their formation. Specifically, the present-day kinematics of cluster stars, such as rotation and velocity dispersion, seems to be related to the initial configuration of the system. In recent work (Milone et al. 2018), we analyzed for the first time the kinematics of the different stellar populations in NGC0104 (47Tucanae) over a large field of view, exploiting the Gaia Data Release 2 proper motions combined with multi-band ground-based photometry. In this paper, based on the work by Cordoni et al. (2019), we extend this analysis to six GCs, namely NGC0288, NGC5904 (M5), NGC6121 (M4), NGC6752, NGC6838 (M71) and further explore NGC0104. Among the analyzed clusters only NGC0104 and NGC5904 show significant rotation on the plane of the sky. Interestingly, multiple stellar populations in NGC5904 exhibit different rotation curves.
Figures
Reference graph
Works this paper leans on
-
[1]
Bellini, A., Vesperini, E., Piotto, G., et al.\ 2015, ApJL, 810, L13
work page 2015
-
[2]
Bianchini, P., van der Marel, R. P., del Pino, A. et al.\ 2018, mnras, 481, 2125
work page 2018
-
[3]
Cordoni, G., Milone, A. P., Mastrobuono-Battisti, A. et al.\ 2019, arXiv, e-prints, 1905.09908
arXiv 2019
-
[4]
H \'e nault-Brunet, V., Gieles, M., Agertz, O., & Read, J. I.\ 2015, MNRAS, 450, 1164
work page 2015
-
[5]
Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al.\ 2018, A&A, 616, A1
work page 2018
-
[6]
Gaia Collaboration, Helmi, A., van Leeuwen, F., et al.\ 2018, A&A, 616, 12G
work page 2018
- [7]
-
[8]
P., Piotto, G., Renzini, A., et al.\ 2017, MNRAS, 464, 3636
Milone, A. P., Piotto, G., Renzini, A., et al.\ 2017, MNRAS, 464, 3636
2017
Show all 13 references
-
[9]
P., Marino, A
Milone, A. P., Marino, A. F., Mastrobuono-Battisti, A., & Lagioia, E. P.\ 2018, MNRAS, 479, 5005
2018
-
[10]
B., Heyl, J., Anderson, J., et al.\ 2013, ApJL, 771, L15
Richer, H. B., Heyl, J., Anderson, J., et al.\ 2013, ApJL, 771, L15
2013
-
[11]
Sollima, A., Baumgardt, H., Hilker, M.,\ 2019, MNRAS, 485, 1460
2019
-
[12]
B., Pancino, E., Zocchi, A
Stetson, P. B., Pancino, E., Zocchi, A. et al.\ 2019, MNRAS, 485, 3042
2019
-
[13]
Vesperini, E., McMillan, S. L. W., D'Antona, F., & D'Ercole, A.\ 2013, MNRAS, 429, 1913
2013
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.