pith. sign in

arxiv: 2605.20503 · v1 · pith:ORNZEZQUnew · submitted 2026-05-19 · 🌌 astro-ph.GA

The chemo-dynamical complexity of {ω} Centauri: different kinematics for different populations

Pith reviewed 2026-05-21 06:21 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords omega Centauriglobular clustermultiple populationschemo-dynamicskinematicsaluminum abundanceAPOGEEGaia proper motions
0
0 comments X

The pith

In ω Centauri, aluminum-rich stars sit closer to the center and move on more radial orbits than aluminum-poor stars.

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

The paper tests whether chemically distinct stars in the globular cluster ω Centauri also occupy different positions and follow different orbits. It separates red-giant stars from APOGEE into five groups based on their abundances of eight elements, then groups those components into two broader families according to aluminum content. The aluminum-rich family is more concentrated toward the center and shows stronger radial anisotropy in its velocities, while the aluminum-poor family stays closer to isotropy out to four half-light radii. Both families rotate at similar rates. These kinematic distinctions supply direct dynamical evidence that the cluster experienced multiple formation channels rather than a single simple collapse.

Core claim

A Gaussian Mixture Model in eight-dimensional abundance space applied to APOGEE DR17 red-giant stars identifies five chemical components that naturally separate into an aluminum-rich family and an aluminum-poor family. The aluminum-rich stars are more centrally concentrated and display stronger radial velocity anisotropy, whereas the aluminum-poor stars remain nearer to isotropy across the probed radial range. The two families nevertheless share a common rotation pattern. These chemo-dynamical differences are measured from the inner regions out to approximately four half-light radii and represent the first such link between detailed chemical tagging and internal kinematics in ω Centauri.

What carries the argument

Gaussian Mixture Model in eight-dimensional chemical-abundance space that partitions stars into components subsequently grouped by aluminum enrichment to expose spatial and kinematic contrasts.

If this is right

  • The spatial concentration and orbital anisotropy of the aluminum-rich family imply that these stars formed in a denser, more centrally located environment than the aluminum-poor family.
  • Shared rotation between the families indicates they experienced at least one common dynamical phase after their separate chemical enrichment.
  • The radial profiles extending to four half-light radii show that the kinematic differences persist into the cluster outskirts.
  • A formation path combining hierarchical assembly inside a dwarf-galaxy potential with later centrally concentrated star formation accounts for both the chemical and dynamical observations.

Where Pith is reading between the lines

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

  • Similar chemical-tagging plus kinematic analysis could be applied to other massive clusters suspected of ex-situ origins to test whether the same two-family pattern appears.
  • If the aluminum-rich stars formed later and more centrally, their orbits may have been affected by subsequent dynamical friction or relaxation that the aluminum-poor stars avoided.
  • The observed chemo-dynamical split supplies a concrete observational target for simulations that model both merger-driven accretion and in-situ star formation within the same cluster.

Load-bearing premise

The five chemical groups returned by the mixture model reflect genuinely separate stellar populations with distinct formation histories rather than arising from measurement noise or model choice.

What would settle it

Re-running the mixture model on the identical APOGEE sample while varying the number of components or inflating abundance errors so that the radial anisotropy and concentration differences between the aluminum-rich and aluminum-poor families disappear.

Figures

Figures reproduced from arXiv: 2605.20503 by A. Mastrobuono-Battisti, G. Pagnini, P. Bianchini, P. Di Matteo.

Figure 1
Figure 1. Figure 1: Chemical clustering and spatial distribution of APOGEE stars in [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Chemical, spatial, and global kinematic properties of the Al-poor and Al-rich populations. Top left: [Al/ [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Radial anisotropy profiles of the two chemically defined populations. Left and middle panels: [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Kinematic profiles of the three velocity components. Top row: mean tangential PM, radial PM, and line-of-sight velocity [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
read the original abstract

The origin of {\omega} Centauri remains one of the key open problems in stellar dynamics and chemical evolution. Its large abundance spreads and multiple populations suggest a formation history more complex than that of a typical globular cluster. We investigate whether the chemical sub-populations identified in APOGEE DR17 also exhibit distinct spatial and kinematic signatures. We analyse a sample of APOGEE DR17 red-giant stars using a Gaussian Mixture Model in an eight-dimensional chemical-abundance space. The resulting chemical components are combined with Gaia proper motions and APOGEE line-of-sight velocities to derive intrinsic mean velocities and velocity dispersions in all three observable directions. We measure both global kinematic quantities and radial profiles for each chemically defined group, extending from the inner regions to ~ 4 half-light radii. The Gaussian Mixture Model identifies five chemical components, which, when examined through their radial cumulative distributions, naturally group into two broader families characterised by lower and higher aluminium enrichment. The two families differ significantly in their spatial and kinematic properties: the Al-rich stars are more centrally concentrated and exhibit stronger radial anisotropy than the Al-poor stars, which remain closer to isotropy over the radial range probed. Despite these significant differences, the two populations share a common rotation pattern. This work represents the first chemo-dynamical study of {\omega} Cen linking detailed chemical tagging to internal kinematics from the inner regions to the cluster outskirts. A formation path involving both hierarchical assembly within a dwarf-galaxy potential and centrally concentrated, chemically enriched star formation offers a natural explanation for the observed chemo-dynamical complexity.

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 / 2 minor

Summary. The manuscript applies a Gaussian Mixture Model (GMM) to eight-dimensional APOGEE DR17 abundance data for red-giant stars in ω Centauri, identifying five chemical components. These are manually grouped post-hoc into two Al-rich and Al-poor families on the basis of mean aluminium abundance. Combining the assignments with Gaia proper motions and APOGEE line-of-sight velocities, the authors derive mean velocities, dispersions, and anisotropy profiles in radial bins out to ~4 half-light radii. They report that the Al-rich family is more centrally concentrated and shows stronger radial anisotropy, while the Al-poor family remains closer to isotropy, although both families share a common rotation pattern. The results are interpreted as evidence for a complex formation history involving both hierarchical assembly and centrally concentrated, chemically enriched star formation.

Significance. If the five-component decomposition is shown to be robust against abundance uncertainties and model-complexity choices, the work would provide one of the first direct links between detailed chemical tagging and spatially resolved kinematics across the full radial extent of ω Centauri. This would strengthen the case for a dwarf-galaxy progenitor and supply quantitative constraints on the relative spatial and orbital properties of the chemically distinct populations.

major comments (2)
  1. [§3] §3 (GMM analysis): The manuscript reports a single run of the GMM with K=5 but provides no BIC/AIC values, cross-validation scores, or stability tests under realistic abundance perturbations. Because the subsequent manual grouping into Al-rich/Al-poor families and all kinematic contrasts are derived directly from these component labels, the absence of such validation leaves open whether the reported differences in central concentration and radial anisotropy are stable or sensitive to the particular decomposition.
  2. [§4] §4 (kinematic profiles): The radial anisotropy and cumulative-distribution profiles for the two families are computed after hard assignment to the grouped families. No propagation of the GMM posterior probabilities into the velocity-dispersion or anisotropy uncertainties is described; this omission is load-bearing for the claim that the Al-rich family exhibits 'stronger radial anisotropy' while the Al-poor family remains 'closer to isotropy'.
minor comments (2)
  1. [Abstract] The abstract states that the components 'naturally group' into two families; a brief statement of the quantitative criterion (e.g., a threshold in mean [Al/Fe]) would remove ambiguity.
  2. [§2] Notation for the three velocity components (radial, tangential, line-of-sight) should be defined once at first use and used consistently in all figures and tables.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed report. The comments highlight important points regarding the robustness of our Gaussian Mixture Model (GMM) decomposition and the treatment of assignment uncertainties in the kinematic analysis. We have revised the manuscript to incorporate additional validation metrics and to propagate posterior probabilities into the derived kinematic profiles and uncertainties. Our point-by-point responses follow.

read point-by-point responses
  1. Referee: [§3] §3 (GMM analysis): The manuscript reports a single run of the GMM with K=5 but provides no BIC/AIC values, cross-validation scores, or stability tests under realistic abundance perturbations. Because the subsequent manual grouping into Al-rich/Al-poor families and all kinematic contrasts are derived directly from these component labels, the absence of such validation leaves open whether the reported differences in central concentration and radial anisotropy are stable or sensitive to the particular decomposition.

    Authors: We agree that explicit validation of the GMM choice strengthens the analysis. In the revised manuscript we now report BIC and AIC values for K = 2 to 8, which indicate that K = 5 is preferred. We have also added a stability test in which the eight-dimensional abundance vectors are perturbed 100 times by drawing from the reported APOGEE uncertainties; each perturbed catalogue is re-fit with the GMM. The resulting component centroids, particularly in [Al/Fe], remain consistent to within 0.05 dex, preserving the same post-hoc grouping into Al-rich and Al-poor families. Five-fold cross-validation scores are likewise reported and support K = 5. These additions demonstrate that the reported spatial and kinematic differences are insensitive to the precise decomposition within the range of abundance uncertainties. revision: yes

  2. Referee: [§4] §4 (kinematic profiles): The radial anisotropy and cumulative-distribution profiles for the two families are computed after hard assignment to the grouped families. No propagation of the GMM posterior probabilities into the velocity-dispersion or anisotropy uncertainties is described; this omission is load-bearing for the claim that the Al-rich family exhibits 'stronger radial anisotropy' while the Al-poor family remains 'closer to isotropy'.

    Authors: We acknowledge the value of propagating the full GMM posteriors. The revised analysis now weights each star by its posterior probability of belonging to the Al-rich or Al-poor family when computing velocity dispersions, anisotropy parameters, and their uncertainties in each radial bin. Uncertainties are estimated via bootstrap resampling that respects the posterior weights. The updated profiles confirm that the Al-rich family remains more radially anisotropic than the Al-poor family at >2σ significance across the radial range, while both families share the same mean rotation. The relevant text, equations, and figures have been updated to reflect this weighted approach. revision: yes

Circularity Check

0 steps flagged

No significant circularity: independent chemistry and kinematics measurements

full rationale

The paper applies a standard Gaussian Mixture Model to external APOGEE DR17 abundances in eight-dimensional space to define chemical components, then measures spatial distributions and kinematics separately using Gaia proper motions combined with APOGEE line-of-sight velocities. The reported differences in central concentration and radial anisotropy between the manually grouped Al-rich and Al-poor families are direct empirical measurements on these independently assigned groups, not quantities that reduce by construction to the GMM fit parameters or any self-citation. No uniqueness theorems, ansatzes, or renamings of known results are invoked to force the chemo-dynamical distinctions. The derivation chain remains self-contained against external data benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the GMM clustering and the assumption that chemical groups map to distinct dynamical populations; the number of components is chosen to fit the data.

free parameters (1)
  • Number of GMM components = 5
    Set to five to capture the chemical substructure in the APOGEE sample.
axioms (1)
  • domain assumption Velocity distributions within each chemically defined group are Gaussian
    Invoked when deriving intrinsic mean velocities and dispersions from the combined Gaia and APOGEE data.

pith-pipeline@v0.9.0 · 5834 in / 1408 out tokens · 43942 ms · 2026-05-21T06:21:35.396645+00:00 · methodology

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

64 extracted references · 64 canonical work pages

  1. [1]

    2022, ApJS, 259, 35

    Abdurro’uf, Accetta, K., Aerts, C., et al. 2022, ApJS, 259, 35

  2. [2]

    2019, ApJ, 886, 57 Alvarez Garay, D

    Alfaro-Cuello, M., Kacharov, N., Neumayer, N., et al. 2019, ApJ, 886, 57 Alvarez Garay, D. A., Mucciarelli, A., Bellazzini, M., Lardo, C., & Ventura, P. 2024, A&A, 681, A54

  3. [3]

    & van der Marel, R

    Anderson, J. & van der Marel, R. P. 2010, ApJ, 710, 1032

  4. [4]

    I., Vesperini, E., & Dalessandro, E

    Aros, F. I., Vesperini, E., & Dalessandro, E. 2025, A&A, 699, A44

  5. [5]

    2010, in Society of Photo-Optical In- strumentation Engineers (SPIE) Conference Series, V ol

    Bacon, R., Accardo, M., Adjali, L., et al. 2010, in Society of Photo-Optical In- strumentation Engineers (SPIE) Conference Series, V ol. 7735, Ground-based and Airborne Instrumentation for Astronomy III, ed. I. S. McLean, S. K. Ram- say, & H. Takami, 773508

  6. [6]

    & Hilker, M

    Baumgardt, H. & Hilker, M. 2018, MNRAS, 478, 1520

  7. [7]

    & Vasiliev, E

    Baumgardt, H. & Vasiliev, E. 2021, MNRAS, 505, 5957

  8. [8]

    R., Piotto, G., Anderson, J., et al

    Bedin, L. R., Piotto, G., Anderson, J., et al. 2004, ApJ, 605, L125

  9. [9]

    & Freeman, K

    Bekki, K. & Freeman, K. C. 2003, MNRAS, 346, L11

  10. [10]

    & Tsujimoto, T

    Bekki, K. & Tsujimoto, T. 2016, ApJ, 831, 70

  11. [11]

    P., et al

    Bellini, A., Anderson, J., van der Marel, R. P., et al. 2014, ApJ, 797, 115

  12. [12]

    R., Piotto, G., et al

    Bellini, A., Bedin, L. R., Piotto, G., et al. 2010, AJ, 140, 631

  13. [13]

    R., et al

    Bellini, A., Libralato, M., Bedin, L. R., et al. 2018, ApJ, 853, 86

  14. [14]

    P., Anderson, J., et al

    Bellini, A., Milone, A. P., Anderson, J., et al. 2017, ApJ, 844, 164

  15. [15]

    J., & Epps, E

    Butler, D., Dickens, R. J., & Epps, E. 1978, ApJ, 225, 148

  16. [16]

    2020, ApJ, 891, 167

    Calamida, A., Zocchi, A., Bono, G., et al. 2020, ApJ, 891, 167

  17. [17]

    & Lia, C

    Carraro, G. & Lia, C. 2000, A&A, 357, 977

  18. [18]

    C., Dotter, A., et al

    Clontz, C., Seth, A. C., Dotter, A., et al. 2024, ApJ, 977, 14

  19. [19]

    C., Wang, Z., et al

    Clontz, C., Seth, A. C., Wang, Z., et al. 2025, ApJ, 984, 162

  20. [20]

    P., Marino, A

    Cordoni, G., Milone, A. P., Marino, A. F., et al. 2020, ApJ, 898, 147

  21. [21]

    P., Marino, A

    Dondoglio, E., Milone, A. P., Marino, A. F., et al. 2026, A&A, 705, A2

  22. [22]

    R., Sollima, A., Pancino, E., et al

    Ferraro, F. R., Sollima, A., Pancino, E., et al. 2004, ApJ, 603, L81

  23. [23]

    W., Lang, D., & Goodman, J

    Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2021, A&A, 649, A1

  24. [24]

    G., Carretta, E., & Bragaglia, A

    Gratton, R. G., Carretta, E., & Bragaglia, A. 2012, A&A Rev., 20, 50

  25. [25]

    2016, MNRAS, 461, 3620 Häberle, M., Neumayer, N., Clontz, C., et al

    Guillard, N., Emsellem, E., & Renaud, F. 2016, MNRAS, 461, 3620 Häberle, M., Neumayer, N., Clontz, C., et al. 2025, ApJ, 983, 95

  26. [26]

    2004, A&A, 422, L9

    Hilker, M., Kayser, A., Richtler, T., & Willemsen, P. 2004, A&A, 422, L9

  27. [27]

    Johnson, C. I. & Pilachowski, C. A. 2010, ApJ, 722, 1373

  28. [28]

    2018, Monthly Notices of the Royal Astronomical Society, 473, 5591

    Kamann, S., Husser, T.-O., Dreizler, S., et al. 2018, Monthly Notices of the Royal Astronomical Society, 473, 5591

  29. [29]

    W., Joo, J

    Lee, Y . W., Joo, J. M., Sohn, Y . J., et al. 1999, Nature, 402, 55

  30. [30]

    I., Baumgardt, H., Cabrera-Ziri, I., et al

    Leitinger, E. I., Baumgardt, H., Cabrera-Ziri, I., et al. 2025, A&A, 694, A184 Lloyd Evans, T. 1977, MNRAS, 181, 591

  31. [31]

    R., Patterson, R

    Majewski, S. R., Patterson, R. J., Dinescu, D. I., et al. 2000, in Liege Inter- national Astrophysical Colloquia, V ol. 35, Liege International Astrophysical Colloquia, ed. A. Noels, P. Magain, D. Caro, E. Jehin, G. Parmentier, & A. A. Thoul, 619

  32. [32]

    K., & Ferraro, F

    Marcolini, A., Sollima, A., D’Ercole, A., Gibson, B. K., & Ferraro, F. R. 2007, MNRAS, 382, 443

  33. [33]

    F., Milone, A

    Marino, A. F., Milone, A. P., Piotto, G., et al. 2012, ApJ, 746, 14

  34. [34]

    F., Milone, A

    Marino, A. F., Milone, A. P., Piotto, G., et al. 2011, ApJ, 731, 64

  35. [35]

    F., Milone, A

    Marino, A. F., Milone, A. P., Renzini, A., et al. 2019, MNRAS, 487, 3815

  36. [36]

    C., Schiavon, R

    Mason, A. C., Schiavon, R. P., Kamann, S., et al. 2026, MNRAS, 547, stag433

  37. [37]

    & Perets, H

    Mastrobuono-Battisti, A. & Perets, H. B. 2013, ApJ, 779, 85

  38. [38]

    & Perets, H

    Mastrobuono-Battisti, A. & Perets, H. B. 2021, MNRAS, 505, 2548

  39. [39]

    1997, AJ, 114, 1087 Mészáros, S., Masseron, T., Fernández-Trincado, J

    Mayor, M., Meylan, G., Udry, S., et al. 1997, AJ, 114, 1087 Mészáros, S., Masseron, T., Fernández-Trincado, J. G., et al. 2021, MNRAS, 505, 1645

  40. [40]

    P., Marino, A

    Milone, A. P., Marino, A. F., Piotto, G., et al. 2015, ApJ, 808, 51

  41. [41]

    2020, A&A Rev., 28, 4

    Neumayer, N., Seth, A., & Böker, T. 2020, A&A Rev., 28, 4

  42. [42]

    S., Neumayer, N., Clontz, C., et al

    Nitschai, M. S., Neumayer, N., Clontz, C., et al. 2023, ApJ, 958, 8

  43. [43]

    & Bessell, M

    Norris, J. & Bessell, M. S. 1975, ApJ, 201, L75

  44. [44]

    Norris, J. E. & Da Costa, G. S. 1995, ApJ, 447, 680

  45. [45]

    2025, A&A, 693, A155

    Pagnini, G., Di Matteo, P., Haywood, M., et al. 2025, A&A, 693, A155

  46. [46]

    R., Bellazzini, M., Piotto, G., & Zoccali, M

    Pancino, E., Ferraro, F. R., Bellazzini, M., Piotto, G., & Zoccali, M. 2000, ApJ, 534, L83

  47. [47]

    R., & Bellazzini, M

    Pancino, E., Galfo, A., Ferraro, F. R., & Bellazzini, M. 2007, ApJ, 661, L155

  48. [48]

    2024, MNRAS, 528, 4941

    Pechetti, R., Kamann, S., Krajnovi´c, D., et al. 2024, MNRAS, 528, 4941

  49. [49]

    E., Frogel, J

    Persson, S. E., Frogel, J. A., Cohen, J. G., Aaronson, M., & Matthews, K. 1980, ApJ, 235, 452

  50. [50]

    A., Seitzer, P., Arnold, R., et al

    Reijns, R. A., Seitzer, P., Arnold, R., et al. 2006, A&A, 445, 503

  51. [51]

    2007, MNRAS, 376, 405

    Romano, D., Matteucci, F., Tosi, M., et al. 2007, MNRAS, 376, 405

  52. [52]

    P., Phillips, S

    Schiavon, R. P., Phillips, S. G., Myers, N., et al. 2024, MNRAS, 528, 1393

  53. [53]

    C., et al

    Souza, S., Neumayer, N., Seth, A. C., et al. 2026, arXiv e-prints, arXiv:2603.23589

  54. [54]

    Suntzeff, N. B. & Kraft, R. P. 1996, AJ, 111, 1913

  55. [55]

    2016, MN- RAS, 457, 4525

    Tailo, M., Di Criscienzo, M., D’Antona, F., Caloi, V ., & Ventura, P. 2016, MN- RAS, 457, 4525

  56. [56]

    A., Vesperini, E., & Varri, A

    Tiongco, M. A., Vesperini, E., & Varri, A. L. 2019, MNRAS, 487, 5535

  57. [57]

    I., & Korchagin, V

    Tsuchiya, T., Dinescu, D. I., & Korchagin, V . I. 2003, ApJ, 589, L29

  58. [58]

    I., & Dinescu, D

    Tsuchiya, T., Korchagin, V . I., & Dinescu, D. I. 2004, MNRAS, 350, 1141 van de Ven, G., van den Bosch, R. C. E., Verolme, E. K., & de Zeeuw, P. T. 2006, A&A, 445, 513 van Leeuwen, F., Le Poole, R. S., Reijns, R. A., Freeman, K. C., & de Zeeuw, P. T. 2000, A&A, 360, 472

  59. [59]

    & Baumgardt, H

    Vasiliev, E. & Baumgardt, H. 2021, MNRAS, 505, 5978

  60. [60]

    2025, A&A, 696, A168

    Vernekar, N., Lucatello, S., Kuzma, P., & Spina, L. 2025, A&A, 696, A168

  61. [61]

    2021, MNRAS, 502, 4290

    Vesperini, E., Hong, J., Giersz, M., & Hypki, A. 2021, MNRAS, 502, 4290

  62. [62]

    G., & Cassisi, S

    Villanova, S., Geisler, D., Gratton, R. G., & Cassisi, S. 2014, ApJ, 791, 107

  63. [63]

    R., et al

    Villanova, S., Piotto, G., King, I. R., et al. 2007, ApJ, 663, 296

  64. [64]

    C., Clontz, C., et al

    Wang, Z., Seth, A. C., Clontz, C., et al. 2026, arXiv e-prints, arXiv:2603.01041 Article number, page 10 of 14 G. Pagnini et al.:ωCen 0.4 0.2 0.0 0.2 0.4 0.6 0.8 [Mg/Fe] 0.50 0.25 0.00 0.25 0.50 0.75 1.00 1.25 [Al/Fe] Components Component 0 Component 1 Component 2 Component 3 Component 4 1.8 1.6 1.4 1.2 1.0 0.8 [Fe/H] Fig. A.1: Alternative representation ...