Robust ellipticity measurements of 29 Galactic globular clusters
Pith reviewed 2026-05-15 14:40 UTC · model grok-4.3
The pith
A new robust method shows rotation aligns with the short axis in ten globular clusters, indicating it drives their flattening.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that a newly developed robust ellipticity estimator, tested on simulated data and applied to 29 Galactic globular clusters, yields reliable shape parameters that align with rotation axes in NGC 104, NGC 1261, NGC 2808, NGC 3201, NGC 5286, NGC 5904, NGC 5986, NGC 6205, NGC 6341 and NGC 7078, thereby confirming rotation as the primary driver of their flattening while leaving room for anisotropy or tidal effects in others.
What carries the argument
The robust ellipticity estimator that corrects biases in kernel-density ellipse fits and principal-component analysis, used together with the V/σ diagram inside the isotropic oblate-rotator framework.
If this is right
- Dynamical models and simulations of globular clusters must incorporate measured flattening rather than assume sphericity.
- The method enables reliable shape measurements even for the small stellar samples typical of multiple-population studies.
- Rotation is established as the leading cause of flattening for the ten clusters where axes align.
- Velocity anisotropy or tidal forces must be considered when interpreting shapes of the remaining clusters.
- Improved ellipticity values provide a cleaner test of formation and evolution scenarios for the Galactic globular-cluster system.
Where Pith is reading between the lines
- The technique could be extended to extragalactic globular clusters to test whether rotation remains the dominant flattening mechanism outside the Milky Way.
- If applied to clusters at different evolutionary stages, the measurements might reveal how tides progressively alter shapes over time.
- Accounting for these shapes in N-body models would tighten constraints on the initial conditions of globular-cluster formation.
- The same bias-correction approach might improve shape estimates for other low-density stellar systems such as open clusters or dwarf galaxies.
Load-bearing premise
The systematic errors measured on mock data fully represent the errors present in real photometric observations of globular clusters.
What would settle it
High-resolution kinematic maps of the ten listed clusters that show no statistical alignment between their measured rotation axes and the semi-minor axes of the fitted ellipses would falsify the claim that rotation is the main driver.
read the original abstract
Globular clusters (GCs) exhibit varying degrees of flattening (ellipticity), which may provide insight into their internal dynamics and evolution histories. Commonly used methods to measure ellipticity, such as ellipse fitting of density contours and principal component analysis, often produce biased results, especially for clusters that are nearly round or have few observable stars. Using a combination of ground-based and space-based photometry, we investigate the shapes of 29 Galactic GCs. To that end, we test two commonly used methods: an ellipse fit to a kernel density profile and a principal component analysis. We find that both methods suffer from bias arising when the number of stars is small or the cluster is close to round. To solve this issue, we develop a robust method to measure the ellipticity of GCs, test it extensively on mock data, and apply it to the 29 Milky Way GCs in our sample. Using the $V/\sigma$ diagram used in the isotropic oblate rotator framework, we examine potential causes for the flattening, including rotation and velocity anisotropy. For ten clusters: NGC~104, NGC~1261, NGC~2808, NGC 3201, NGC 5286, NGC 5904, NGC 5986, NGC 6205, NGC 6341, and NGC 7078 we identify a very good agreement between the rotation angle and semi-minor axis of the ellipse, further corroborating the findings that rotation is the main driver of the ellipticity. The $V/\sigma$ diagram reveals that velocity anisotropy or tides could also be important in shaping the GCs. The robust method developed provides reliable measurements of the ellipticity of GCs, emphasising the importance of taking into account the flattening in theoretical models and simulations. It also offers a promising way to investigate the shapes of multiple stellar populations within GCs, where only small samples are usually available.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a robust ellipticity estimator for globular clusters after demonstrating biases in standard ellipse-fitting to kernel density profiles and principal component analysis, particularly for low star counts or near-round systems. The new method is validated extensively on mock data before application to combined ground- and space-based photometry for 29 Galactic GCs. For the ten clusters NGC 104, NGC 1261, NGC 2808, NGC 3201, NGC 5286, NGC 5904, NGC 5986, NGC 6205, NGC 6341, and NGC 7078, the authors report very good agreement between the photometric semi-minor axis and the rotation angle, supporting rotation as the dominant driver of flattening; the V/σ diagram is used to assess possible contributions from velocity anisotropy or tides.
Significance. A reliable ellipticity measurement technique for GCs with limited star samples would be a useful addition to the field, and the explicit mock validation plus the reported alignment for ten clusters constitute concrete strengths. If the method remains unbiased once real photometric systematics are fully accounted for, the results would strengthen the case that rotation shapes GC ellipticities and would encourage inclusion of flattening in dynamical models.
major comments (2)
- [Mock validation section] Mock validation section: the tests on synthetic data do not incorporate spatially varying completeness, field-star contamination, or PSF-induced centroid shifts that are present in the actual ground+space photometry. These effects can shift recovered ellipticities and position angles by amounts comparable to the reported agreement for the ten listed clusters, directly affecting the claim that rotation is the main driver.
- [Results for the ten clusters] Results for the ten clusters (NGC 104 etc.): the central claim of 'very good agreement' between rotation angle and semi-minor axis is presented without error bars on the ellipticity or angle measurements, quantitative fit statistics, or explicit sample-exclusion criteria. This leaves the statistical robustness of the alignment difficult to evaluate.
minor comments (2)
- [Abstract] The abstract states that the method was 'tested extensively on mock data' but does not specify the range of star counts, ellipticities, or contamination levels explored.
- [V/σ analysis] In the V/σ analysis, clarify the precise definition of the isotropic oblate rotator framework and whether any clusters were excluded from the diagram.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed comments, which have helped us identify areas where the manuscript can be clarified and strengthened. We address each major comment below and outline the revisions we will make.
read point-by-point responses
-
Referee: [Mock validation section] Mock validation section: the tests on synthetic data do not incorporate spatially varying completeness, field-star contamination, or PSF-induced centroid shifts that are present in the actual ground+space photometry. These effects can shift recovered ellipticities and position angles by amounts comparable to the reported agreement for the ten listed clusters, directly affecting the claim that rotation is the main driver.
Authors: We agree that our mock tests focused primarily on the dominant biases from low star counts and near-round shapes and did not fully replicate all real-data systematics. In the revised manuscript we will expand the mock validation section to include additional simulations that incorporate spatially varying completeness and field-star contamination. We will also add a quantitative discussion of the possible influence of PSF-induced centroid shifts on the recovered position angles, drawing on the characteristics of our combined ground- and space-based photometry. These changes will allow us to assess whether such effects could materially affect the reported alignments for the ten clusters. revision: yes
-
Referee: [Results for the ten clusters] Results for the ten clusters (NGC 104 etc.): the central claim of 'very good agreement' between rotation angle and semi-minor axis is presented without error bars on the ellipticity or angle measurements, quantitative fit statistics, or explicit sample-exclusion criteria. This leaves the statistical robustness of the alignment difficult to evaluate.
Authors: We acknowledge that the current presentation lacks the quantitative details needed for a full statistical assessment. In the revised manuscript we will report formal uncertainties on both the ellipticity and position-angle measurements obtained with the new robust method. We will also provide quantitative measures of agreement (e.g., the Pearson correlation coefficient between the photometric semi-minor axis and the literature rotation angle, together with associated p-values) and will explicitly state the selection criteria used to identify the ten clusters, including any sample-exclusion rules applied. These additions will make the statistical robustness of the alignment transparent to readers. revision: yes
Circularity Check
Independent mock validation keeps ellipticity and rotation-angle comparison non-circular
full rationale
The paper develops a new ellipticity estimator after showing bias in standard ellipse-fitting and PCA methods, validates the estimator on separate mock catalogs that are independent of the real 29-GC photometry, and only then applies it to the observed clusters. The reported agreement between photometric semi-minor-axis position angle and kinematic rotation angle for the ten listed clusters is therefore a cross-comparison between two independently measured quantities rather than a quantity that is forced by the fit itself. No equation in the derivation chain reduces the final ellipticity values or the claimed rotation–shape agreement to a self-defined input; any self-citations are peripheral and not load-bearing for the central result.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.