Pith. sign in

REVIEW 3 major objections 4 minor 31 references

Imprint of the galactic acceleration scale on globular cluster systems

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read Density profiles of globular cluster systems break at the radius where the gravitational acceleration of the stars equals the galactic acceleration scale $a_0$.

desk verdict A new, honestly reported coincidence between GC system break radii and the a0 acceleration radius, but the breaks' reality is not yet established because no smooth-profile null is tested. read the letter →

arxiv 1908.04783 v1 pith:T3DZROOM submitted 2019-08-13 astro-ph.GA

classification astro-ph.GA
keywords globularclustersystemsgalacticaccelerationscalemodifiedNewtoniandynamics(MOND)radialrelationearly-typegalaxiesdarkmatterhalosdensityprofilebreaksbrokenpowerlaw
topics Dark Matter
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper reports that the radial density profiles of globular cluster systems in 17 early-type galaxies do not follow a single power law: they steepen at a break radius $r_{\rm br}$, and that radius coincides with the point where the gravitational acceleration generated by the stars equals $a_0$, the galactic acceleration scale known from rotation curves and the radial acceleration relation. The average ratio is $r_{\rm br}/r_{\rm acc,N}=0.8\pm0.3$ under Newtonian gravity and $r_{\rm br}/r_{\rm acc,M}=1.2\pm0.5$ under MOND, while other characteristic radii of the galaxies match the breaks less well. The authors argue that, if the match holds, the break radius becomes a photometric probe of the gravitational field and, in the standard dark-matter picture, of the dark halo concentration. A sympathetic reader would care because this adds a new observable to the short list of phenomena that single out the acceleration scale $a_0$.

What carries the argument

The carrying object is a broken power law for the globular cluster volume density, $\rho(r)=\rho_0 r^a$ inside the break and $\rho(r)=\rho_0 r_{\rm br}^{a-b} r^b$ outside it, with $r_{\rm br}$ the break radius. The identity doing the work is $r_{\rm br}\approx r_{\rm acc}$, where $r_{\rm acc}$ is the radius at which the gravitational acceleration from the stars equals $a_0$ under Newtonian gravity ($r_{\rm acc,N}$) or under the MOND interpolation ($r_{\rm acc,M}$); the average ratios are $0.8\pm0.3$ and $1.2\pm0.5$, respectively. The proposed mechanism behind the identity is framework-dependent: in the Lambda-CDM reading the gravitational potential changes its slope at the radius where stellar and dark-halo accelerations are equal, so infalling satellites are tidally stripped at the same radius but deposit their globular clusters with different radial spreads inside and outside that point; in MOND the break sits at $r_{\rm acc}$ because the dynamics change regime there, with external-field effects and enhanced dynamical friction acting only beyond the transition.

What would settle it

Fit the same archival globular cluster profiles with a single power law and with a Sérsic profile and compare the fits statistically; if the broken power law is not clearly preferred, or if the recovered $r_{\rm br}$ moves outside its quoted errors when red and blue subpopulations are fitted separately or when incompleteness is modeled, the reported match would be an artifact rather than an imprint of $a_0$.

Watch

Extended reading notes

Core claim

In the sample of 17 early-type galaxies, plus three projections of one simulated galaxy, the number density profiles of the globular cluster systems are well described by a broken power law that steepens beyond a break radius $r_{\rm br}$. The paper's central discovery is that $r_{\rm br}$ is nearly equal to $r_{\rm acc}$, the radius where the gravitational acceleration generated by the stars equals $a_0=1.2\times10^{-10}\,\mathrm{m\,s^{-2}}$: the average ratio is $0.8\pm0.3$ for Newtonian gravity and $1.2\pm0.5$ for the MOND interpolation. The match with the galaxy's effective radius, halo scale radius, and the red-blue crossover radius is substantially worse. The paper interprets this as an imprint of the acceleration scale on the globular cluster system, proposing that in the Lambda-CDM picture the break marks the radius where the dark halo starts to dominate the potential, and in MOND it marks the Newtonian-to-deep-MOND transition. The authors also note tentative evidence that globular cluster systems can reveal halo concentration as well as halo mass.

Load-bearing premise

The claim rests on the break radii $r_{\rm br}$ being real, sharp features of the globular cluster density profiles rather than artifacts of fitting a broken power law to an intrinsically smooth profile, of survey incompleteness, or of the superposition of red and blue cluster subpopulations.

Editorial extensions

If this is right

  • Photometry alone would provide the acceleration-scale radius: a confirmed $r_{\rm br}\approx r_{\rm acc}$ relation lets the globular cluster break radius stand in for kinematic tracers in galaxies where spectroscopy is impractical.
  • In the Lambda-CDM interpretation, combining the break radius with the standard relation between globular cluster number and halo mass would estimate both the mass and the concentration of the dark halo from imaging.
  • In the MOND interpretation, the break at $r_{\rm acc}$ is a prediction, and the outer slopes of globular cluster density profiles should approach the predicted $r^{-\alpha_\infty}$ behavior with $\alpha_\infty$ between about 3.5 and 4.5.
  • The match provides a new observational constraint on galaxy formation, and current Lambda-CDM galaxy formation simulations do not reproduce it: the simulated galaxy's break radii differ from both $r_{\rm acc}$ and $r_{\rm sh}$ by factors of a few.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the correlation is confirmed on an independent sample spanning a wider mass range, the globular cluster break radius could become a distance-free observable that tracks $a_0$ in galaxies where the stellar mass profile is known, effectively turning globular cluster systems into a gravitational-field probe.
  • A clean discrimination test between the two explanations would compare break radii in galaxies with similar stellar masses but very different dark halo concentrations: the Lambda-CDM version predicts $r_{\rm br}$ to track $r_{\rm sh}$, while the MOND version predicts it to track $r_{\rm acc}$.
  • The red-versus-blue globular cluster dichotomy remains the most obvious alternative driver of a break; separating the two color populations in each galaxy and measuring whether the break tracks the red-blue crossover radius would directly test that alternative.
  • The same argument should extend to globular-cluster-rich ultra-diffuse galaxies, where imaging-based estimates of halo concentration could be checked against independent dynamical measurements.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. This paper reports that the break radii of the radial number density profiles of globular cluster (GC) systems in 17 early-type galaxies coincide with the radius at which the gravitational acceleration from the stellar mass equals the galactic acceleration scale a0, both under Newtonian and MOND gravity. The authors compare the break radii rbr with other characteristic radii and find that rbr matches racc more closely than the effective radius, NFW scale radii, or the red/blue GC equality radius. They propose possible explanations in the ΛCDM and MOND frameworks and suggest that GC breaks could be used to estimate dark halo concentrations. The analysis uses previously published break radii from Bílek et al. (2019) and independently computed a0 radii, with no constant fitted to force agreement.

Significance. If the reported coincidence is real, it would add a new empirical manifestation of the acceleration scale in early-type galaxies and potentially offer a purely photometric probe of dark halo structure. The paper's main strengths are the parameter-free nature of the rbr–racc comparison and the checks against survey incompleteness and red/blue GC segregation. However, the significance is conditional: the central claim requires that the breaks are genuine features of the density profiles rather than fitting artifacts, and the statistical support for the rbr–racc match is not yet demonstrated. The proposed ΛCDM interpretation of a link between rbr and halo concentration also lacks quantitative backing.

major comments (3)
  1. [Section 2, Eq. (1), Appendix A] The central claim that GC systems exhibit density breaks at rbr rests entirely on the broken power-law fits from Bílek et al. (2019), but the paper gives no evidence that the broken power-law model is statistically preferred over a single power law or a Sérsic profile for any of the 17 galaxies. A broken power law can absorb the smooth curvature of a Sérsic profile, placing an artificial break at a radius that scales with the effective radius Re; since racc also correlates with galaxy mass, the reported rbr–racc agreement could then be a by-product of the mass–size relation rather than a physical break. I ask the authors to add model-comparison tests (e.g., likelihood-ratio or ΔAIC) against single power-law and Sérsic models, and to show the distribution of recovered rbr/Re for Sérsic-only synthetic profiles for comparison.
  2. [Section 2, Table A.1] The average ratios rbr/racc,N = 0.8 ± 0.3 and rbr/racc,M = 1.2 ± 0.5 are quoted without any measure of statistical significance. With only 17 galaxies and many rbr uncertainties spanning factors of 2–3, the reader cannot tell whether the agreement with unity is meaningful; please provide a formal test (e.g., a bootstrap confidence interval for the mean ratio, the number of galaxies within 1σ or 2σ of equality, or a correlation coefficient between rbr and racc). The paper also excludes the three simulated galaxies from the averages because they disagree; since that disagreement is relevant to the proposed explanations, the exclusion should be justified and its impact on the conclusions discussed.
  3. [Section 3, Table 1] The tentative ΛCDM claim that GC break radii can reveal halo concentrations is not quantitatively supported. The average ratio of fitted halo scale radii to break radii, rsh,f/rbr = 0.8 ± 0.7, and the large scatter indicate that the relation is not established, particularly because the illustrative example NGC 3115 was selected for its good match. I recommend a quantitative test of whether the rbr values are consistent with the fitted rsh values within the factor-of-two uncertainties on rs,f, or a statement that this part remains speculative.
minor comments (4)
  1. [Footnote 3] The phrase 'The breaks in the GC system profiles at racc can thus be taken as a MOND prediction' overstates the predictive status, because the comparison was made after the break radii were already known; please clarify that the MOND interpretation is a post-hoc identification rather than an a priori prediction.
  2. [Appendix A heading and abstract] The Appendix A heading 'Fits of radial number density profiles the GC systems' is missing the word 'of'; also, the abstract contains 'Lambda cold dark matter ( ΛCDM)' with an extra space after 'Lambda'.
  3. [Section 2] The statement that the acceleration profiles have no local maxima near the break radii is not accompanied by any figure or table; please specify the method and stellar mass profile used, or refer to a specific analysis in Bílek et al. (2019).
  4. [Table A.1] The meaning of the quoted error bars as 1σ limits is stated in the Appendix, but the table itself does not indicate this; consider adding a note to the table caption.

Circularity Check

1 steps flagged · score 3.0 of 10

The rbr–racc correlation is independently measured and not fitted, but the paper retroactively labels the already-fitted break radii as a MOND prediction (footnote 3), partially compromising that framing.

  1. fitted input called prediction [Section 3, footnote 3 (MOND paragraph)]
    "This consideration actually led MB to compare the break radii with racc, and after that with all the other characteristic radii. The breaks in the GC system profiles at racc can thus be taken as a MOND prediction."

    The rbr values are fitted parameters from the broken-power-law fits of Eq. (1) in Bílek et al. (2019). The footnote explicitly states that the MOND-motivated comparison was made only after the breaks were already known. Presenting the fitted breaks as 'a MOND prediction' turns the confirming quantity into the same fitted input that motivated the hypothesis, so the prediction is post hoc rather than independent. The underlying rbr–racc correlation is not statistically forced, because no parameter is fitted to make them agree, but the paper's 'prediction' wording overstates the independence of the evidence.

full rationale

The core empirical claim is not circular: rbr comes from fits to GC number-density profiles in the companion paper, while racc,N and racc,M are computed from stellar masses, distances, and the external constant a0 using Newtonian gravity or the McGaugh et al. (2016) MOND interpolation. No free parameter is adjusted to force rbr = racc, so the comparison is not equivalent to its inputs by construction. The reliance on Bílek et al. (2019) is a legitimate same-author prior fit rather than an unverified self-citation, and the ΛCDM halo radii are independent fitting products, not definitions of the target relation. The one genuinely circular element is the MOND-prediction framing exposed by footnote 3: the authors admit the MOND comparison was motivated after seeing the breaks, yet call the breaks a MOND prediction. That is a fitted parameter relabeled as a prediction. Whether a broken power law is statistically preferred over a smooth Sérsic profile is a real correctness risk about the reality of the breaks, but it is not a circularity of the present derivation and is therefore not scored here.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

No new entities are introduced. The paper's central comparison rests on previously fitted break radii, adopted stellar masses and distances, the external constant a0, and an extrapolated MOND interpolation function. The Lambda-CDM interpretation additionally depends on NFW halo parameters from scaling relations and kinematic fits. No free parameters are tuned in this paper to force the agreement.

free parameters (3)
  • rbr (break radius) per galaxy from broken power-law fits = Table A.1 (e.g., N821: 2.3 +2 -0.8 arcmin)
    The central comparison uses these fitted radii; they come from the same authors' prior paper and carry asymmetrical uncertainties that are large (factors of 2 or more) for several galaxies.
  • Stellar masses (via M/L) and distances of the sample galaxies = Table 1 (log M, d)
    racc is computed from the stellar mass distribution; these inputs are adopted from Bilek et al. (2019) and other literature, and their systematic errors directly shift racc. They are not fitted in this paper but are load-bearing.
  • NFW halo scale radii rs,s and rs,f = Table 1
    Used in the Lambda-CDM interpretation (rsh comparison). rs,s comes from scaling relations (Behroozi et al. 2013; Diemer and Kravtsov 2015) and rs,f from Jeans fits in Bilek et al. (2019); both carry factor-of-two uncertainties.
assumptions (5)
  • domain assumption The globular cluster system behaves as a relaxed tracer population whose number density profile reflects the gravitational potential (Jeans equation applies).
    The paper invokes the Jeans equation to interpret the breaks and refers to the companion kinematics paper; if GCs are not an equilibrium tracer, the interpretation fails.
  • domain assumption The stellar mass distribution is known and can be approximated as spherical to compute the acceleration field.
    racc is computed from the stellar distribution using spherical symmetry and a constant M/L; most ETGs are not spherical.
  • domain assumption The MOND interpolation function of McGaugh et al. (2016) (their Eq. 4) is valid for these early-type galaxies.
    The paper uses this empirical interpolation to compute racc,M; applying a relation calibrated on spiral rotation curves to ETGs is an extrapolation.
  • domain assumption The external field effect is negligible for the MOND acceleration calculations.
    racc,M is computed as for isolated galaxies; several sample members reside in clusters or groups (e.g., NGC 1399 in Fornax, NGC 4486 in Virgo), where the external field effect could alter the acceleration profile. This is not discussed in the computation.
  • standard math Standard mathematical tools (Abel transform, hypergeometric functions, Gaussian likelihood) are used for profile fitting and uncertainty estimation.
    Used in Appendix A for converting volume to surface density profiles and computing parameter uncertainties.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Imprint of the galactic acceleration scale on globular cluster systems." pith.science (2026). https://pith.science/paper/T3DZROOM

@misc{pith2026190804783,
  author       = {Pith},
  title        = {Pith review of: Imprint of the galactic acceleration scale on globular cluster systems},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/T3DZROOM}},
  note         = {Machine review of arXiv:1908.04783}
}
abstract

We report that the density profiles of globular cluster (GC) systems in a sample of 17 early-type galaxies (ETGs) show breaks at the radii where the gravitational acceleration exerted by the stars equals the galactic acceleration scale $a_0$ known from the radial acceleration relation or MOND. The match with the other characteristic radii in the galaxy is not that close. We propose possible explanations in the frameworks of the $\Lambda$CDM model and MOND. We find tentative evidence that in the $\Lambda$CDM context, GCs reveal not only the masses of the dark halos through the richness of the GC systems but also the concentrations through the break radii of the GC systems.

Figures

Figures reproduced from arXiv: 1908.04783 by the authors.

Figure 1
Figure 1. Demonstration that the break radii of GC systems are [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Illustration of a possible explanation of the breaks of the [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

31 extracted references · 21 canonical work pages

  1. [1]

    S., Wechsler, R

    Behroozi, P. S., Wechsler, R. H., & Conroy, C. 2013, ApJ, 770, 57

  2. [2]

    & Milgrom, M

    Bekenstein, J. & Milgrom, M. 1984, ApJ, 286, 7

  3. [3]

    A., Beasley, M

    Bekki, K., Forbes, D. A., Beasley, M. A., & Couch, W. J. 2003, MNRAS, 344, 1334 Bílek, M., Samurovi´c, S., & Renaud, F. 2019, A&A, 625, A32 Bílek, M., Thies, I., Kroupa, P., & Famaey, B. 2018, A&A, 614, A59

  4. [4]

    Brodie, J. P. & Strader, J. 2006, ARA&A, 44, 193

  5. [5]
  6. [6]

    & Tiret, O

    Combes, F. & Tiret, O. 2010, in American Institute of Physics Conference Series, V ol. 1241, American Institute of Physics Conference Series, ed. J.-M. Alimi & A. Fuözfa, 154–161 Côté, P., Marzke, R. O., & West, M. J. 1998, ApJ, 501, 554

  7. [7]

    & Kravtsov, A

    Diemer, B. & Kravtsov, A. V . 2015, ApJ, 799, 108

  8. [8]

    Durazo, R., Hernandez, X., Cervantes Sodi, B., & Sánchez, S. F. 2017, ApJ, 837, 179

Show all 31 references
  1. [9]

    Durazo, R., Hernandez, X., Cervantes Sodi, B., & Sanchez, S. F. 2018, ApJ, 863, 107

  2. [10]

    & McGaugh, S

    Famaey, B. & McGaugh, S. S. 2012, Living Reviews in Relativity, 15, 10

  3. [11]

    Forbes, D. A. 2017, MNRAS, 472, L104

  4. [12]

    E., Harris, G

    Harris, W. E., Harris, G. L., & Hudson, M. J. 2015, ApJ, 806, 36

  5. [13]

    Hernandez, X., Cortés, R. A. M., & Scarpa, R. 2017, MNRAS, 464, 2930

  6. [14]

    & Jiménez, M

    Hernandez, X. & Jiménez, M. A. 2012, ApJ, 750, 9

  7. [15]

    2015, Canadian Journal of Physics, 93, 169

    Kroupa, P. 2015, Canadian Journal of Physics, 93, 169

  8. [16]

    S., Lelli, F., & Schombert, J

    McGaugh, S. S., Lelli, F., & Schombert, J. M. 2016, Physical Review Letters, 117, 201101

  9. [17]

    1983, ApJ, 270, 365

    Milgrom, M. 1983, ApJ, 270, 365

  10. [18]

    1984, ApJ, 287, 571

    Milgrom, M. 1984, ApJ, 287, 571

  11. [19]

    2012, Physical Review Letters, 109, 131101

    Milgrom, M. 2012, Physical Review Letters, 109, 131101

  12. [20]

    2014, MNRAS, 437, 2531

    Milgrom, M. 2014, MNRAS, 437, 2531

  13. [21]

    F., Benítez-Llambay, A., Fattahi, A., et al

    Navarro, J. F., Benítez-Llambay, A., Fattahi, A., et al. 2017, MNRAS, 471, 1841

  14. [22]

    F., Frenk, C

    Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563

  15. [23]

    2008, MNRAS, 386, 2194

    Nipoti, C., Ciotti, L., Binney, J., & Londrillo, P. 2008, MNRAS, 386, 2194

  16. [24]

    2007, MNRAS, 381, L104

    Nipoti, C., Londrillo, P., & Ciotti, L. 2007, MNRAS, 381, L104

  17. [25]

    W., Forbes, D

    Pota, V ., Graham, A. W., Forbes, D. A., et al. 2013, MNRAS, 433, 235

  18. [26]

    2017, MNRAS, 465, 3622 Samurovi´c, S

    Renaud, F., Agertz, O., & Gieles, M. 2017, MNRAS, 465, 3622 Samurovi´c, S. 2014, A&A, 570, A132

  19. [27]

    & Falomo, R

    Scarpa, R. & Falomo, R. 2010, A&A, 523, A43

  20. [28]

    2003, A&A, 405, L15

    Scarpa, R., Marconi, G., & Gilmozzi, R. 2003, A&A, 405, L15

  21. [29]

    2007, A&A, 462, L9

    Scarpa, R., Marconi, G., Gilmozzi, R., & Carraro, G. 2007, A&A, 462, L9

  22. [30]

    & Combes, F

    Tiret, O. & Combes, F. 2008, in Astronomical Society of the Pacific Conference

  23. [31]

    2013, ApJ, 762, 39 van Dokkum, P., Abraham, R., Brodie, J., et al

    Tonini, C. 2013, ApJ, 762, 39 van Dokkum, P., Abraham, R., Brodie, J., et al. 2016, ApJ, 828, L6 Article number, page 4 of 6 M. Bílek et al.: Imprint of the galactic acceleration scale on globular cluster systems Appendix A: Fits of radial number density profiles the GC systems...

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.