pith. sign in

arxiv: 2605.23802 · v1 · pith:75NBTWZXnew · submitted 2026-05-22 · 🌌 astro-ph.GA

Astrophysical Parameters of 5056 Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones

Olcay Plevne , Furkan Akbaba This is my paper

Pith reviewed 2026-05-25 03:12 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords open star clustersGaia DR3Bayesian nested samplingPARSEC isochronescolour-magnitude diagramsphotometric metallicitycluster parametershomogeneous catalogue
0
0 comments X

The pith

Bayesian nested sampling with PARSEC isochrones yields a uniform catalogue of age, distance, metallicity and extinction for 5056 open clusters from Gaia DR3 data.

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

The paper establishes a single, automated procedure that fits the same set of stellar evolution models to the colour-magnitude diagrams of every cluster in a large public list. By leaving initial metallicity free and supplying only broad priors drawn from Gaia parallaxes, XP spectra and a dust map, the method produces photometric [Fe/H] values together with the other three parameters without any cluster-by-cluster manual tuning. A reader would care because earlier compilations mixed results from many different techniques and pipelines; a homogeneous set removes those hidden inconsistencies and therefore supports cleaner statistical studies of how star clusters trace the structure and chemical evolution of the Milky Way disk.

Core claim

The authors apply Bayesian nested sampling to PARSEC isochrones on Gaia DR3 colour-magnitude diagrams for all 5056 clusters drawn from the Unified Cluster Catalogue. Metallicity is treated as a free parameter throughout, priors come only from Gaia astrometry, XP spectra where available, and the SFD dust map, and the entire process runs without manual intervention on individual objects. The resulting catalogue reports age, [Fe/H], heliocentric distance and E(G_BP-G_RP) together with posterior chains; 3766 clusters (74.5 percent) pass a fit-quality threshold of eta_fit greater than or equal to 0.80.

What carries the argument

Bayesian nested sampling applied uniformly to PARSEC isochrones fitted to Gaia DR3 colour-magnitude diagrams, with priors from Gaia astrometry, XP spectra and the SFD dust map.

If this is right

  • The catalogue supplies photometric metallicities for every cluster because initial metallicity remains a free parameter in every fit.
  • Seventy-four point five percent of the sample meets the eta_fit greater than or equal to 0.80 quality cut and spans ages from a few million years to 5.5 Gyr and distances out to nearly 19 kpc.
  • The full posterior chains are released so that downstream studies can propagate uncertainties directly rather than adopting single-point estimates.
  • No external parameter catalogue is used to anchor the fits, so the results are independent of prior compilations.

Where Pith is reading between the lines

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

  • The released posterior distributions could be combined with kinematic data to study cluster dissolution timescales as a function of metallicity.
  • Repeating the identical pipeline on future Gaia releases would test whether the current distance and extinction values change systematically for the most distant clusters.
  • Cross-matching the photometric [Fe/H] values against high-resolution spectroscopic surveys for the overlapping clusters would quantify any residual model-dependent bias.

Load-bearing premise

PARSEC isochrones accurately represent the observed colour-magnitude diagrams of open clusters across the full range of ages, metallicities and extinctions encountered in the sample.

What would settle it

Systematic offsets larger than the reported uncertainties between the catalogue ages or distances and independent measurements obtained for a substantial subset of the same clusters by another technique would falsify the central claim.

Figures

Figures reproduced from arXiv: 2605.23802 by Furkan Akbaba, Olcay Plevne.

Figure 1
Figure 1. Figure 1: Aitoff projection of the sky distribution in Galactic coordinates (ℓ, b) of the 5 056 open clusters processed in this work. Galactic centre is at the projection centre; longitude increases to the left. Class AA clusters (N = 2 121, dark teal) have membership solutions confirmed by two or more independent studies; class AB clusters (N = 2 767, light teal) have a reliable but singly confirmed solution (see G… view at source ↗
Figure 2
Figure 2. Figure 2: Dereddened colour–magnitude diagrams for four representative clusters spanning the full age range of the catalogue. Stars are colour-coded by membership probability (Pmem, coolwarm scale, 0.5–1.0) with black outlines; grey points are stars below the p ≥ 0.75 threshold. The orange curve shows the best-fitting PARSEC isochrone at the posterior median parameters. Key parameters are annotated in each panel: lo… view at source ↗
Figure 3
Figure 3. Figure 3: Distributions of the four derived parameters for the full sample (grey) and the high-quality subset (ηfit ≥ 0.80, teal). Top left: Logarithmic age. Top right: Heliocentric distance. Bottom left: Metallicity [Fe/H]. Bottom right: Extinction in the G band, AG [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Distribution of the fit quality indicator ηfit for all 5 056 clusters (grey) and the high-quality subset (ηfit ≥ 0.80, teal; N = 3 766, 74.5 per cent). The dashed vertical line marks the adopted quality threshold. The distribution peaks sharply near ηfit ≈ 0.95, indicating that the majority of the sample is well-constrained, while the tail below 0.80 captures clusters with degenerate or poorly populated CM… view at source ↗
Figure 5
Figure 5. Figure 5: Posterior distributions (corner plots) for four clusters spanning the full quality range of ηfit. Each panel shows the joint and marginal posteriors for the four free parameters: log(Age/yr), initial metallicity Zini, heliocentric distance D (pc), and colour excess E(GBP−GRP); vertical dashed lines mark the 16th, 50th, and 84th percentiles. The inset CMD in each panel shows stars colour-coded by membership… view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of [Fe/H] posterior widths (∆[Fe/H] ≡ [Fe/H]84 − [Fe/H]16) for clusters with (teal, N = 755) and without (brown, N = 3 011) an informative Gaia XP spectrophotometric metallicity prior, restricted to the high-quality subset (ηfit ≥ 0.80). Left: Normalized histogram; vertical dashed lines mark the respective medians. The dotted vertical line at 0.68 dex indicates the 68% width of the broad uniform… view at source ↗
Figure 7
Figure 7. Figure 7: Dereddened colour–magnitude diagrams for 30 NGC open clusters with fit quality ηfit ≥ 0.88, arranged in order of increasing age from top to bottom (six rows, five columns per row). Each panel shows probable members (p ≥ 0.75) colour-coded by membership probability (coolwarm scale) with the best-fitting PARSEC isochrone overlaid (orange curve); grey points are stars below the membership threshold. The y-axi… view at source ↗
Figure 8
Figure 8. Figure 8: Scatter comparison of age (log(Age/yr)), distance, and AG between this work and three photometric catalogues. Rows (top to bottom): Hunt et al. (2023), Dias et al. (2021), Cantat-Gaudin et al. (2020). Grey points show the full sample; teal points show clusters with ηfit ≥ 0.80. Statistics (N, ⟨∆⟩, σ) for each panel are given in the lower right corner. Dashed lines indicate the 1:1 relation [PITH_FULL_IMAG… view at source ↗
Figure 9
Figure 9. Figure 9: One-to-one comparison between photometric [Fe/H] (this work) and MWM spectroscopic [Fe/H] from the OCCAM DR19 catalogue (J. M. Otto et al. 2026) for the ηfit ≥ 0.80 subset (teal) and the full sample (grey). The dashed line shows the 1:1 relation. Statistics (N, mean offset, σ, MAD) are shown in the panel. cated Gaussians centred on the Gaia DR3 me￾dian parallax of cluster members. Metallicity pri￾ors are i… view at source ↗
read the original abstract

We present a homogeneous catalogue of fundamental astrophysical parameters -- age, metallicity ([Fe/H]), heliocentric distance, and colour excess $E(G_{\mathrm{BP}}-G_{\mathrm{RP}})$ -- for 5,056 open star clusters drawn from the Unified Cluster Catalogue (UCC). All parameters are derived uniformly from Gaia Data Release 3 (DR3) colour-magnitude diagrams via Bayesian Nested Sampling with PARSEC stellar isochrones, with no manual intervention on individual clusters. Initial metallicity $Z_{\mathrm{ini}}$ is treated as a free parameter throughout, yielding a photometric [Fe/H] estimate for every cluster. Physically motivated priors -- parallax-based distances from Gaia DR3 astrometry, spectrophotometric metallicity constraints from Gaia XP spectra where available, and interstellar reddening from the Schlegel-Finkbeiner-Davis dust map -- reduce CMD degeneracies without anchoring the fit to any external parameter catalogue. Of the 5,056 clusters, 3,766 (74.5\%) satisfy the fit-quality criterion $\eta_{\mathrm{fit}} \ge 0.80$. This high-quality subset spans ages 0.003-5.5~Gyr ($\log(\mathrm{Age/yr})$ median $8.33 \pm 0.34$~dex), heliocentric distances 88-19,011 pc (median 2,150~pc), metallicities $-1.17 \le \mathrm{[Fe/H]} \le +0.42$~dex (median $+0.002$ dex), and extinctions up to $A_G = 7.37$~mag (median 1.07~mag). The catalogue is made publicly available via CDS/VizieR; the complete nested-sampling posterior chains are archived on Zenodo.

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

1 major / 0 minor

Summary. The manuscript presents a homogeneous catalogue of age, [Fe/H], heliocentric distance and E(G_BP-G_RP) for 5056 open clusters from the UCC. Parameters are obtained uniformly via Bayesian nested sampling applied to Gaia DR3 colour-magnitude diagrams using PARSEC isochrones, with Z_ini treated as a free parameter and priors drawn from Gaia astrometry, XP spectra and the SFD dust map. No manual intervention is applied; 3766 clusters (74.5 %) satisfy η_fit ≥ 0.80, and the full posterior chains are released publicly.

Significance. If the PARSEC isochrones are shown to be an adequate representation of the observed CMDs across the reported parameter ranges and the priors introduce no net bias, the resulting catalogue would constitute a substantial, uniformly derived resource for open-cluster and Galactic-disk studies. The public archiving of complete nested-sampling chains is a clear strength that enhances reproducibility and downstream use.

major comments (1)
  1. [Abstract] Abstract (and fitting procedure): The claim that the derived parameters are free of coherent systematic bias rests on the assumption that PARSEC isochrones faithfully reproduce Gaia DR3 CMDs for all 5056 clusters (age 0.003–5.5 Gyr, [Fe/H] −1.17 to +0.42, distances to 19 kpc). No quantitative test against alternative isochrone libraries, no comparison with well-studied clusters having independent parameters, and no assessment of possible offsets arising from convective overshooting, helium enrichment or bolometric corrections are described. Because Z_ini is a free parameter and the isochrones are invoked at every step of the nested sampling, any model mismatch would propagate directly into the reported ages, distances and metallicities even for the η_fit ≥ 0.80 subset.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive review and for underscoring the need to address potential model-dependent systematics. We respond point-by-point below and outline revisions that will strengthen the manuscript without altering its core methodology or results.

read point-by-point responses

  1. Referee: [Abstract] Abstract (and fitting procedure): The claim that the derived parameters are free of coherent systematic bias rests on the assumption that PARSEC isochrones faithfully reproduce Gaia DR3 CMDs for all 5056 clusters (age 0.003–5.5 Gyr, [Fe/H] −1.17 to +0.42, distances to 19 kpc). No quantitative test against alternative isochrone libraries, no comparison with well-studied clusters having independent parameters, and no assessment of possible offsets arising from convective overshooting, helium enrichment or bolometric corrections are described. Because Z_ini is a free parameter and the isochrones are invoked at every step of the nested sampling, any model mismatch would propagate directly into the reported ages, distances and metallicities even for the η_fit ≥ 0.80 subset.

    Authors: We agree that the manuscript does not present quantitative comparisons to alternative isochrone sets (e.g., MIST or BaSTI) or direct validations against clusters with independent spectroscopic or asteroseismic parameters. Our work focuses on delivering a uniformly derived catalogue using a single, widely adopted library (PARSEC) together with Gaia-informed priors; the η_fit threshold is intended only as an internal consistency metric, not as proof of absolute accuracy. We will revise the abstract to remove any implication of zero systematic bias and add a dedicated paragraph in the Discussion section that (i) cites existing literature on PARSEC versus other isochrones for open clusters, (ii) notes that offsets in age, distance or metallicity of order 0.1–0.3 dex or 10–20 % in distance remain possible, and (iii) recommends that users treat the catalogue as a homogeneous but model-dependent resource. No new fits or external validations will be performed, but the added text will make the limitation explicit. revision: yes

Circularity Check

0 steps flagged

No circularity: parameters derived by direct Bayesian fit to independent Gaia CMD data

full rationale

The paper derives cluster parameters (age, [Fe/H], distance, extinction) via Bayesian nested sampling applied to observed Gaia DR3 colour-magnitude diagrams, using PARSEC isochrones as the forward model and external priors from Gaia astrometry, XP spectra, and the SFD dust map. No equation or step reduces any output quantity to a fitted input by construction, renames a prior result, or relies on a self-citation chain for its central claim. The fitting procedure is self-contained against the supplied data and model; the choice of PARSEC is an external modelling assumption whose validity is outside the scope of circularity analysis.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the accuracy of PARSEC isochrones for open-cluster CMDs and on the priors being sufficient to break degeneracies without bias; Z_ini is explicitly treated as free.

free parameters (1)
  • Z_ini
    Initial metallicity is treated as a free parameter for each cluster to produce a photometric [Fe/H] estimate.
axioms (1)
  • domain assumption PARSEC isochrones accurately represent the colour-magnitude diagrams of the open clusters in the sample across the fitted age, metallicity and extinction range
    The entire derivation consists of matching observations to these models.

pith-pipeline@v0.9.0 · 5871 in / 1449 out tokens · 56323 ms · 2026-05-25T03:12:19.451777+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

28 extracted references · 28 canonical work pages · 4 internal anchors

  1. [1]

    2023, A&AS, 267, 8, doi: 10.3847/1538-4365/acd53e Astropy Collaboration, Robitaille, T

    Andrae, R., Rix, H.-W., & Chandra, V. 2023, A&AS, 267, 8, doi: 10.3847/1538-4365/acd53e Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068

    work page doi:10.3847/1538-4365/acd53e 2023

  2. [2]

    Estimating distances from parallaxes. V: Geometric and photogeometric distances to 1.47 billion stars in Gaia Early Data Release 3

    Demleitner, M., & Andrae, R. 2021, AJ, 161, 147, doi: 10.3847/1538-3881/abd806

    work page internal anchor Pith review doi:10.3847/1538-3881/abd806 2021

  3. [3]

    Onions, J., Ascasibar, Y ., Behroozi, P., Casado, J., Elahi, P., Han, J., Knebe, A., Lux, H., Merch´an, M

    Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127, doi: 10.1111/j.1365-2966.2012.21948.x

    work page doi:10.1111/j.1365-2966.2012.21948.x 2012

  4. [4]

    Campello, R. J. G. B., Moulavi, D., & Sander, J. 2013, Advances in Knowledge Discovery and Data Mining, 160

    work page 2013

  5. [5]

    2022, Universe, 8, 111, doi: 10.3390/universe8020111

    Cantat-Gaudin, T. 2022, Universe, 8, 111, doi: 10.3390/universe8020111

    work page doi:10.3390/universe8020111 2022

  6. [6]

    2020, A&A, 640, A1, doi: 10.1051/0004-6361/202038192

    Cantat-Gaudin, T., Anders, F., Castro-Ginard, A., et al. 2020, A&A, 640, A1, doi: 10.1051/0004-6361/202038192

    work page doi:10.1051/0004-6361/202038192 2020

  7. [7]

    2021, A&A, 652, A25, doi: 10.1051/0004-6361/202039951

    Casamiquela, L., Tarricq, Y., Soubiran, C., et al. 2021, A&A, 652, A25, doi: 10.1051/0004-6361/202039951

    work page doi:10.1051/0004-6361/202039951 2021

  8. [8]

    2014, MNRAS, 444, 2525, doi: 10.1093/mnras/stu1605 De Angeli, F., Weiler, M., Montegriffo, P., et al

    Chen, Y., Girardi, L., Bressan, A., et al. 2014, MNRAS, 444, 2525, doi: 10.1093/mnras/stu1605 De Angeli, F., Weiler, M., Montegriffo, P., et al. 2023, A&A, 674, A2, doi: 10.1051/0004-6361/202243680

    work page doi:10.1093/mnras/stu1605 2014

  9. [9]

    S., Alessi, B

    Dias, W. S., Alessi, B. S., Moitinho, A., & L´ epine, J. R. D. 2002, A&A, 389, 871, doi: 10.1051/0004-6361:20020668

    work page doi:10.1051/0004-6361:20020668 2002

  10. [10]

    S., Monteiro, H., Moitinho, A., et al

    Dias, W. S., Monteiro, H., Moitinho, A., et al. 2021, MNRAS, 504, 356, doi: 10.1093/mnras/stab770

    work page doi:10.1093/mnras/stab770 2021

  11. [11]

    M., Cunha, K., et al

    Donor, J., Frinchaboy, P. M., Cunha, K., et al. 2020, AJ, 159, 199, doi: 10.3847/1538-3881/ab77bc

    work page doi:10.3847/1538-3881/ab77bc 2020

  12. [12]

    Drimmel, R., & Spergel, D. N. 2001, ApJ, 556, 181, doi: 10.1086/321611 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1, doi: 10.1051/0004-6361/201833051 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2021, A&A, 649, A1, doi: 10.1051/0004-6361/202039657 Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 202...

    work page doi:10.1086/321611 2001

  13. [13]

    L., & Reffert, S

    Hunt, E. L., & Reffert, S. 2021, A&A, 646, A104, doi: 10.1051/0004-6361/202039341

    work page doi:10.1051/0004-6361/202039341 2021

  14. [14]

    L., & Reffert, S

    Hunt, E. L., & Reffert, S. 2023, A&A, 673, A114, doi: 10.1051/0004-6361/202346267

    work page doi:10.1051/0004-6361/202346267 2023

  15. [15]

    L., & Reffert, S

    Hunt, E. L., & Reffert, S. 2024, A&A, 686, A42, doi: 10.1051/0004-6361/202348929

    work page doi:10.1051/0004-6361/202348929 2024

  16. [16]

    V., Piskunov, A

    Kharchenko, N. V., Piskunov, A. E., R¨ oser, S., Schilbach, E., & Scholz, R.-D. 2013, A&A, 558, A53, doi: 10.1051/0004-6361/201322302

    work page doi:10.1051/0004-6361/201322302 2013

  17. [17]

    2020, Gaia Early Data Release 3: Parallax bias versus magnitude, colour, and position, doi: 10.1051/0004-6361/202039653

    Lindegren, L., Klioner, S. A., Hern´ andez, J., et al. 2021, A&A, 649, A2, doi: 10.1051/0004-6361/202039653

    work page doi:10.1051/0004-6361/202039653 2021

  18. [18]

    2018, A&A, 616, A2, doi: 10.1051/0004-6361/201832727

    Lindegren, L., et al. 2018, A&A, 616, A2, doi: 10.1051/0004-6361/201832727

    work page doi:10.1051/0004-6361/201832727 2018

  19. [19]

    2017, ApJ, 835, 77, doi: 10.3847/1538-4357/835/1/77

    Marigo, P., Girardi, L., Bressan, A., et al. 2017, ApJ, 835, 77, doi: 10.3847/1538-4357/835/1/77

    work page doi:10.3847/1538-4357/835/1/77 2017

  20. [20]

    Accelerated Hierarchical Density Clustering

    McInnes, L., & Healy, J. 2017, arXiv e-prints, arXiv:1705.07321, doi: 10.48550/arXiv.1705.07321 M´ esz´ aros, S., Jofr´ e, P., Johnson, J. A., et al. 2025, AJ, 170, 96, doi: 10.3847/1538-3881/ade4b9

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1705.07321 2017

  21. [21]

    M., Frinchaboy, P

    Otto, J. M., Frinchaboy, P. M., Myers, N. R., et al. 2026, AJ, 171, 91, doi: 10.3847/1538-3881/ae28d8 22Plevne & Akbaba

    work page doi:10.3847/1538-3881/ae28d8 2026

  22. [22]

    I., V´ azquez, R

    Perren, G. I., V´ azquez, R. A., Piatti, A. E., et al. 2023, MNRAS, 526, 4107, doi: 10.1093/mnras/stad2826

    work page doi:10.1093/mnras/stad2826 2023

  23. [23]

    F., & Finkbeiner, D

    Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103, doi: 10.1088/0004-637X/737/2/103

    work page internal anchor Pith review doi:10.1088/0004-637x/737/2/103 2011

  24. [24]

    Maps of Dust IR Emission for Use in Estimation of Reddening and CMBR Foregrounds

    Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525, doi: 10.1086/305772

    work page internal anchor Pith review doi:10.1086/305772 1998

  25. [25]

    2004, in American Institute of Physics Conference Series, Vol

    Skilling, J. 2004, in American Institute of Physics Conference Series, Vol. 735, Bayesian Inference and Maximum Entropy Methods in Science and Engineering, ed. R. Fischer, R. Preuss, & U. V. Toussaint, 395–405, doi: 10.1063/1.1835238

    work page doi:10.1063/1.1835238 2004

  26. [26]

    Speagle, J. S. 2020, MNRAS, 493, 3132, doi: 10.1093/mnras/staa278

    work page doi:10.1093/mnras/staa278 2020

  27. [27]

    2021, A&A, 647, A19, doi: 10.1051/0004-6361/202039388

    Tarricq, Y., Soubiran, C., Casamiquela, L., et al. 2021, A&A, 647, A19, doi: 10.1051/0004-6361/202039388

    work page doi:10.1051/0004-6361/202039388 2021

  28. [28]

    2019, ApJ, 877, 116, doi: 10.3847/1538-4357/ab1c61

    Wang, S., & Chen, X. 2019, ApJ, 877, 116, doi: 10.3847/1538-4357/ab1c61

    work page doi:10.3847/1538-4357/ab1c61 2019