pith. sign in

arxiv: 2604.19992 · v2 · submitted 2026-04-21 · 🌌 astro-ph.GA

The apparent Large Magellanic Cloud star cluster age gap

Jonathan H. Klos , Andr\'es E. Piatti This is my paper

Pith reviewed 2026-05-12 01:14 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords Large Magellanic Cloudstar clustersage gapstar formation historycluster initial mass functionobservational completenesscluster mass loss
0
0 comments X

The pith

The apparent age gap in Large Magellanic Cloud star clusters stems from lower past star formation rates and detection limits on faded objects.

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

This paper constructs a model connecting the Large Magellanic Cloud's star formation history to the masses of clusters that formed at each epoch. It uses a power-law relation between the global star formation rate and the maximum initial mass of clusters formed, together with a fixed 5 percent formation efficiency. The model then applies calibrated mass loss and an observational completeness cutoff to predict which clusters remain visible at the present day. It finds that the 4-to-11 Gyr interval featured only low-mass clusters, which have since faded below detection thresholds. The observed gap therefore reflects visibility rather than an absence of cluster formation.

Core claim

In our model, the age gap is a consequence of the star-forming history and current observational limits. The age gap corresponds to a period characterised by a lower star formation rate, whereby no clusters with an initial mass above approximately 2 to 5·10^5 M_⊙ were formed. In the present day, these clusters have become so faint that only few of them have been detected. The pattern of both young-and-bright and old-and-massive clusters being more easily detectable than clusters of intermediate ages might reflect a more general phenomenon and not necessarily one specific to the LMC.

What carries the argument

Power-law relation between maximum initial cluster mass and global star formation rate, with constant 5% cluster-forming efficiency, N-body calibrated mass loss, and observational completeness limit.

If this is right

  • Clusters formed 4-11 Gyr ago had initial masses too low to remain above the completeness limit after mass loss.
  • Reproducing the old globular cluster population requires a higher maximum initial mass during epochs of elevated star formation.
  • A smooth linear change in the maximum-mass relation between 8 and 12 Gyr suffices to match the observed distribution.
  • The same visibility bias may affect age distributions measured in other galaxies.

Where Pith is reading between the lines

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

  • Deep surveys targeting faint, low-mass clusters in the LMC could reveal whether any intermediate-age objects were missed.
  • The same formation-efficiency and mass-loss framework could be applied to other dwarf galaxies to test for similar apparent gaps.
  • If confirmed, the result removes the need to invoke special dynamical disruption mechanisms during the gap interval.

Load-bearing premise

A power-law ties the maximum mass of newly formed clusters directly to the global star formation rate, the efficiency stays fixed at 5 percent, and the maximum mass changes linearly from 8 to 12 Gyr ago.

What would settle it

Detection of a substantial population of clusters aged 4-11 Gyr whose present-day masses imply initial masses above roughly 2 times 10^5 solar masses would falsify the model.

Figures

Figures reproduced from arXiv: 2604.19992 by Andr\'es E. Piatti, Jonathan H. Klos.

Figure 1
Figure 1. Figure 1: CFR and both complete and observable CAF computed for our model using M0 = 105M⊙. present-day mass m(t) is greater than mlim(t). Mlim(R, t) is then the observational completeness limit in initial mass at radius R and age t, such that fobs(t; R) = Z Mmax(t) Mlim(R,t) f(M) dM (15) gives the observable fraction at orbital radius R. Assuming that clusters are homogeneously distributed in the disc between radii… view at source ↗
Figure 3
Figure 3. Figure 3: Observable cluster fraction in our model as a function of look￾back time for Mmax(t) parametrised by β = 1 and M0 in the range 5 · 104M⊙ to 5 · 105M⊙. 2 4 6 8 10 12 14 100 101 102 103 CAF ηobs [Gyr −1 ] M0 = 5 · 104M 2 4 6 8 10 12 14 100 101 102 103 M0 = 105M 2 4 6 8 10 12 14 Age t [Gyr] 100 101 102 103 CAF ηobs [Gyr −1 ] M0 = 2 · 105M 2 4 6 8 10 12 14 Age t [Gyr] 100 101 102 103 M0 = 5 · 105M [PITH_FULL_… view at source ↗
Figure 4
Figure 4. Figure 4: Observable CAF ηobs in our model as a function of cluster age for Mmax(t) parametrised by β = 1 and M0 in the range 5 · 104M⊙ to 5 · 105M⊙. Shaded areas represent uncertainty bands based on SFR un￾certainties. fobs(t; R) > 0. As a result, the smaller fobs(t) becomes, the more sensitive it is to Mmax(t) and thus to differences in M0. If the observable fraction is reduced to 0 in this way, a gap in the range… view at source ↗
Figure 5
Figure 5. Figure 5: CAF ηobs in our model as a function of cluster age for M0 chang￾ing linearly from 105M⊙ to 2·105M⊙ between ages 8 and 12 Gyr (dotted lines). 5.0 5.2 5.4 5.6 5.8 6.0 10−5 10−3 10−1 integrated CIMF Young clusters (t < 2.5 Gyr) 5.0 5.2 5.4 5.6 5.8 6.0 10−5 10−3 10−1 integrated CIMF Age gap clusters (t = 4–11 Gyr) 5.0 5.2 5.4 5.6 5.8 6.0 log M/M 10−5 10−3 10−1 integrated CIMF Old clusters (t > 11 Gyr) [PITH_F… view at source ↗
Figure 6
Figure 6. Figure 6: Time-integrated CIMFs for recently formed clusters, age gap clusters and old clusters in our model for linearly changing M0. For M < 105M⊙, the CIMFs follow a power law with α = 2. In [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
read the original abstract

In the Large Magellanic Cloud (LMC), there have been very few clusters observed with ages between 4 and 11 Gyr. This phenomenon is sometimes referred to as the `LMC age gap'. We constructed a model of the cluster age distribution aimed at reproducing this scenario. We linked the star formation history to the cluster initial mass function via a power-law relation between maximum initial cluster mass and global star formation rate. Using a constant cluster-forming efficiency of 5%, we obtained the cluster formation history. Applying a model of cluster mass loss calibrated using N-body simulations and an observational completeness limit, we computed the observable fraction of initially formed clusters. We were then able to model the cluster age distribution. For a maximum initial cluster mass below $10^5$M$_\odot$ at a star formation rate of 1 M$_\odot$pc$^{-2}$Gyr$^{-1}$, our model reproduced the observed lack of clusters with ages between 4 and 11 Gyr. However, our model required a maximum initial mass at 1 M$_\odot$pc$^{-2}$Gyr$^{-1}$ of at least $2\cdot 10^5$M$_\odot$ in order to reproduce the population of ancient globular clusters. A linear change between maximum initial cluster mass relations from 8 to 12 Gyr reproduced the age gap to a satisfactory extent. In our model, the age gap is a consequence of the star-forming history and current observational limits. The age gap corresponds to a period characterised by a lower star formation rate, whereby no clusters with an initial mass above approximately 2 to 5$\cdot 10^5$M$_\odot$ were formed. In the present day, these clusters have become so faint that only few of them have been detected. The pattern of both young-and-bright and old-and-massive clusters being more easily detectable than clusters of intermediate ages might reflect a more general phenomenon and not necessarily one specific to the LMC.

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 paper constructs a model of the LMC cluster age distribution by connecting the star formation history to the cluster initial mass function via a power-law relating maximum initial cluster mass to global SFR, assuming constant 5% cluster-forming efficiency. It applies N-body-calibrated mass loss and observational completeness limits to compute the observable cluster population, reproducing the 4-11 Gyr age gap when the maximum initial mass is below ~10^5 M⊙ at SFR = 1 M⊙ pc^{-2} Gyr^{-1}. However, reproducing the ancient globular clusters requires a higher normalization (~2×10^5 M⊙), which is achieved by introducing a linear change in the maximum-mass relation between 8 and 12 Gyr. The authors conclude that the gap is apparent, arising from low-SFR epochs that produced no clusters massive enough to remain detectable today.

Significance. If the modeling assumptions hold, the work offers a physically motivated explanation for the LMC age gap as a selection effect driven by SFH, cluster mass loss, and detectability thresholds rather than a true formation hiatus. It demonstrates how young bright clusters and old massive clusters are preferentially observable compared to intermediate-age ones, with potential applicability to other galaxies. The incorporation of N-body-derived mass loss and explicit completeness modeling provides a concrete, testable framework that strengthens the link between global SFR and cluster populations.

major comments (2)
  1. [model description and results (linear change paragraph)] The linear interpolation of the maximum initial cluster mass normalization between 8 and 12 Gyr (described in the model construction and results) is introduced specifically because a constant relation calibrated to the young population underproduces ancient globular clusters while one calibrated to the old population overproduces intermediate-age clusters. This adjustment is load-bearing for reproducing both the gap and the old population simultaneously, yet lacks independent theoretical or extragalactic justification, making the match partly by construction.
  2. [model parameters and abstract] The power-law relation between maximum initial cluster mass and SFR, together with the fixed 5% cluster-forming efficiency, are tuned such that no clusters above ~2-5×10^5 M⊙ form during the low-SFR interval; the paper should demonstrate via sensitivity tests or alternative functional forms whether the gap reproduction is robust or sensitive to these choices (see abstract and model parameters).
minor comments (2)
  1. [model construction] Provide the explicit functional form and any equations for the power-law relation between maximum initial cluster mass and SFR, including how the normalization is defined at SFR = 1 M⊙ pc^{-2} Gyr^{-1}.
  2. [results] Clarify the derivation of the specific range 'approximately 2 to 5·10^5 M⊙' for clusters absent during the gap period, and add a comparison plot of the age distribution with and without the linear ramp to quantify its effect.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed review. We address each major comment below and have revised the manuscript accordingly to improve the robustness and transparency of our modeling choices.

read point-by-point responses

  1. Referee: [model description and results (linear change paragraph)] The linear interpolation of the maximum initial cluster mass normalization between 8 and 12 Gyr (described in the model construction and results) is introduced specifically because a constant relation calibrated to the young population underproduces ancient globular clusters while one calibrated to the old population overproduces intermediate-age clusters. This adjustment is load-bearing for reproducing both the gap and the old population simultaneously, yet lacks independent theoretical or extragalactic justification, making the match partly by construction.

    Authors: We agree that the linear change in normalization is introduced to reconcile the constraints from the young LMC population with the need for sufficiently massive ancient globular clusters. This is a minimal parameterization chosen for simplicity rather than a claim of unique physical validity. In the revised manuscript we have expanded the relevant section to discuss possible physical motivations, including the expectation of higher gas surface densities and pressures at earlier epochs that could permit more massive cluster formation. We have also added sensitivity tests using a step-function change at 10 Gyr and a constant high normalization; these show that the age gap is reproduced whenever the normalization is higher at old ages, indicating that the precise linear form is not critical to the main result. revision: yes

  2. Referee: [model parameters and abstract] The power-law relation between maximum initial cluster mass and SFR, together with the fixed 5% cluster-forming efficiency, are tuned such that no clusters above ~2-5×10^5 M⊙ form during the low-SFR interval; the paper should demonstrate via sensitivity tests or alternative functional forms whether the gap reproduction is robust or sensitive to these choices (see abstract and model parameters).

    Authors: We thank the referee for this suggestion. The adopted power-law slope and 5% efficiency are taken from observational calibrations of young clusters and standard literature values, respectively. In the revised manuscript we have added a dedicated sensitivity analysis varying the power-law index by ±0.2, the normalization by factors of 0.5–2, and the efficiency between 1% and 10%. These tests confirm that the age gap remains a robust feature provided the maximum cluster mass during the low-SFR epoch stays below ~5×10^5 M⊙, consistent with young-cluster observations. We have also tested an alternative broken-power-law form and obtain qualitatively similar results for the gap. revision: yes

Circularity Check

1 steps flagged

Ad hoc linear ramp in max initial cluster mass-SFR relation (8-12 Gyr) fitted to simultaneously reproduce age gap and ancient globular clusters

specific steps
  1. fitted input called prediction [Abstract]
    "For a maximum initial cluster mass below 10^5 M⊙ at a star formation rate of 1 M⊙ pc^{-2} Gyr^{-1}, our model reproduced the observed lack of clusters with ages between 4 and 11 Gyr. However, our model required a maximum initial mass at 1 M⊙ pc^{-2} Gyr^{-1} of at least 2·10^5 M⊙ in order to reproduce the population of ancient globular clusters. A linear change between maximum initial cluster mass relations from 8 to 12 Gyr reproduced the age gap to a satisfactory extent."

    The maximum initial mass normalization is explicitly adjusted (from <10^5 to ≥2×10^5 M⊙) and given a linear time dependence between 8-12 Gyr solely so that the model simultaneously matches the observed gap and the ancient cluster population. The gap is thereby reproduced by construction through parameter tuning rather than as an independent output of the SFH plus completeness model.

full rationale

The derivation links SFH to cluster IMF via an assumed power-law M_max(SFR) and fixed 5% efficiency, then applies mass-loss and completeness models. However, a constant power-law calibrated to young clusters under-produces old globulars, while one calibrated to old clusters over-produces intermediate-age objects. The paper therefore introduces an explicit linear interpolation of the power-law normalization between 8 and 12 Gyr; only this time-dependent adjustment suppresses clusters above ~2-5e5 M⊙ during the low-SFR interval while still allowing the high-mass tail at >12 Gyr. The age-gap 'prediction' therefore reduces to the choice of this fitted ramp rather than emerging independently from the SFH and observational limits.

Axiom & Free-Parameter Ledger

2 free parameters · 3 axioms · 0 invented entities

The claim rests on assumed scaling relations and tuned parameters chosen to match observations rather than first-principles derivations or independent data.

free parameters (2)
  • cluster-forming efficiency = 5%
    Held constant at 5% to convert star formation history into cluster formation history.
  • maximum initial cluster mass at SFR = 1 M_⊙ pc^{-2} Gyr^{-1} = <10^5 M_⊙ (gap) / ≥2×10^5 M_⊙ (old clusters)
    Set below 10^5 M_⊙ to produce the gap and at least 2×10^5 M_⊙ to produce ancient clusters, with linear change between 8-12 Gyr.
axioms (3)
  • domain assumption Power-law relation between maximum initial cluster mass and global star formation rate
    Used to link star formation history to the cluster initial mass function.
  • domain assumption Cluster mass loss model calibrated using N-body simulations
    Applied to compute the observable fraction of clusters at different ages.
  • domain assumption Observational completeness limit
    Determines which clusters remain detectable today.

pith-pipeline@v0.9.0 · 5672 in / 1715 out tokens · 69088 ms · 2026-05-12T01:14:33.889658+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We linked the star formation history to the cluster initial mass function via a power-law relation between maximum initial cluster mass and global star formation rate... A linear change between maximum initial cluster mass relations from 8 to 12 Gyr reproduced the age gap to a satisfactory extent.

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    For a maximum initial cluster mass below 10^5 M_⊙ at a star formation rate of 1 M_⊙ pc^{-2} Gyr^{-1}... our model required a maximum initial mass... of at least 2·10^5 M_⊙... linear change... from 8 to 12 Gyr

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages

  1. [1]

    2008, MNRAS, 390, 759

    Bastian, N. 2008, MNRAS, 390, 759

    work page 2008

  2. [2]

    Baumgardt, H., Parmentier, G., Anders, P., & Grebel, E. K. 2013, MNRAS, 430, 676

    work page 2013

  3. [3]

    J., Beasley, M

    Bekki, K., Couch, W. J., Beasley, M. A., et al. 2004, ApJ, 610, L93

    work page 2004

  4. [4]

    C., Reina-Campos, M., Eadie, G., & Sills, A

    Berek, S. C., Reina-Campos, M., Eadie, G., & Sills, A. 2023, MNRAS, 525, 1902 Da Costa, G. S. 1991, in IAU Symposium, V ol. 148, The Magellanic Clouds, ed. R. Haynes & D. Milne, 183

    work page 2023

  5. [5]

    Ferreira, B. P. L., Dias, B., Santos, J. F. C., et al. 2025, A&A, 695, L9

    work page 2025

  6. [6]

    2020, MNRAS, 499, 4114

    Gatto, M., Ripepi, V ., Bellazzini, M., et al. 2020, MNRAS, 499, 4114

    work page 2020

  7. [7]

    2024, A&A, 690, A164

    Gatto, M., Ripepi, V ., Bellazzini, M., et al. 2024, A&A, 690, A164

    work page 2024

  8. [8]

    2022, A&A, 664, L12

    Gatto, M., Ripepi, V ., Bellazzini, M., et al. 2022, A&A, 664, L12

    work page 2022

  9. [9]

    1997, AJ, 114, 1920

    Geisler, D., Bica, E., Dottori, H., et al. 1997, AJ, 114, 1920

    work page 1997

  10. [10]

    & Zaritsky, D

    Harris, J. & Zaritsky, D. 2009, AJ, 138, 1243

    work page 2009

  11. [11]

    Kruijssen, J. M. D. & Cooper, A. P. 2012, MNRAS, 420, 340

    work page 2012

  12. [12]

    R., McKee, C

    Krumholz, M. R., McKee, C. F., & Bland-Hawthorn, J. 2019, ARA&A, 57, 227

    work page 2019

  13. [13]

    Lamers, H. J. G. L. M., Baumgardt, H., & Gieles, M. 2010, MNRAS, 409, 305

    work page 2010

  14. [14]

    2018, ApJ, 869, 122

    Li, G.-W., Yanny, B., & Wu, Y . 2018, ApJ, 869, 122

    work page 2018

  15. [15]

    Maia, F. F. S., Piatti, A. E., & Santos, J. F. C. 2014, MNRAS, 437, 2005

    work page 2014

  16. [16]

    Massana, P., Ruiz-Lara, T., Noël, N. E. D., et al. 2022, MNRAS, 513, L40

    work page 2022

  17. [17]

    Mateo, M., Hodge, P., & Schommer, R. A. 1986, ApJ, 311, 113

    work page 1986

  18. [18]

    Pagel, B. E. J. & Tautvaisiene, G. 1998, MNRAS, 299, 535

    work page 1998

  19. [19]

    Piatti, A. E. 2014, MNRAS, 437, 1646

    work page 2014

  20. [20]

    Piatti, A. E. 2021, AJ, 161, 199

    work page 2021

  21. [21]

    Piatti, A. E. 2022, MNRAS, 511, L72

    work page 2022

  22. [22]

    E., Alfaro, E

    Piatti, A. E., Alfaro, E. J., & Cantat-Gaudin, T. 2019, MNRAS, 484, L19

    work page 2019

  23. [23]

    E., Illesca, D

    Piatti, A. E., Illesca, D. M. F., Chiarpotti, M., & Butrón, R. 2025, A&A, 702, A108

    work page 2025

  24. [24]

    Piatti, A. E. & Mackey, A. D. 2018, MNRAS, 478, 2164

    work page 2018

  25. [25]

    2016, MNRAS, 461, 519

    Pieres, A., Santiago, B., Balbinot, E., et al. 2016, MNRAS, 461, 519

    work page 2016

  26. [26]

    2020, A&A, 639, L3

    Ruiz-Lara, T., Gallart, C., Monelli, M., et al. 2020, A&A, 639, L3

    work page 2020

  27. [27]

    W., Kroupa, P., & Dib, S

    Zhou, J. W., Kroupa, P., & Dib, S. 2025, MNRAS, 541, 1276 1 https://numpy.org 2 https://matplotlib.org Article number, page 7 of 7

    work page 2025