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

arXiv:2605.23802 · astro-ph.GA · Submitted 2026-05-22 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: I'm Vera.

Jocelyn: And I'm Jocelyn, here with Subrahmanyan.

Vera: Today's paper: "Astrophysical Parameters of 5056 Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones".

Jocelyn: The gist: We present a homogeneous catalogue of fundamental astrophysical parameters — age, metallicity (

Fe/H: ), heliocentric distance,

Vera: First, who's behind it and why it matters.

Title and authors: Vera: So we’ve looked at "Astrophysical Parameters of five thousand fifty-six Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones," and it gives us a super complete set of fundamental parameters for five thousand fifty-six open star clusters using Gaia data <ref:2605.23802#pg1>.

Jocelyn: Right, it’s not just about getting numbers like age or metallicity; it’s how they use that full sample of five thousand fifty-six clusters to build a real baseline for the Milky Way <ref:2605.23802#pg1>.

Jocelyn: And they treat metallicity as a free parameter initially, which means they get this photometric

Fe/H: posterior for every cluster without pre-supposing any abundance model <ref:2605.23802#pg1>.

Vera: So we’ve covered this paper "Astrophysical Parameters of five thousand fifty-six Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones," and it gives us a solid foundation <ref:2605.23802#pg1>.

Jocelyn: It really lays out how important it is to use these large samples to get those tight constraints on fundamental properties like age and metallicity <ref:2605.23802#pg1>.

Subrahmanyan: Definitely. Next time, we might look at how this type of stellar information connects to the larger cosmological models we’re working on <ref:2605.23802#pg1>.

The paper's summary: Vera: So we’ve looked at "Astrophysical Parameters of five thousand fifty-six Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones," and it gives us a super complete set of fundamental parameters for open clusters using Gaia data <ref:2605.23802#pg1>.

Subrahmanyan: That freedom is important because it lets the data drive the result rather than forcing it through an external assumption. It’s a way to see what the clusters actually tell us about their own chemistry <ref:2605.23802#pg1>.

Vera: So we’ve covered "Astrophysical Parameters of five thousand fifty-six Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones," and it gives us a solid foundation <ref:2605.23802#pg1>.

The paper's improvements: Vera: So we’ve just finished looking at "Astrophysical Parameters of five thousand fifty-six Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones," and basically, we’ve got a super complete set of fundamental parameters for five thousand fifty-six open star clusters using Gaia data and PARSEC models <ref:2605.23802#pg1,Astrophysical Parameters of 5056 Open Star Clusters from Bayesian Nested Sampling with>.

Subrahmanyan: And from a theoretical standpoint, this gives us something concrete to test our models on, you know? It helps constrain things like the disc structure and chemical enrichment history <ref:2605.23802#pg1>.

Vera: Exactly. The methodology is pretty solid because they’re using Bayesian Nested Sampling with PARSEC isochrones applied uniformly across all the data. That’s a big win for consistency <ref:2605.23802#pg1>.

Vera: And they do have some smart ways to check their work, too. They define this inlier fraction, eta fit, which has to be greater than zero point eighty for a cluster’s fit to be considered high quality <ref:2605.23802#pg1>.

Jocelyn: So if the data doesn't match the isochrone well enough—if it falls outside that plus or minus three sigma band—the parameters derived should probably be treated with caution <ref:2605.23802#pg1>.

Subrahmanyan: That filtering step is crucial for reliability. It prevents noisy fits from polluting our bigger picture calculations about structure formation in the disc <ref:2605.23802#pg1>.

Vera: They also compare their results against other catalogues and found a consistent distance offset across three different photometric sets, which gives us a real systematic error to account <ref:2605.23802#pg1>.

Jocelyn: That’s interesting because it means we can’t just assume one distance measurement is perfect; there's this shared bias across the literature <ref:2605.23802#pg1>.

Subrahmanyan: It shows that even when using different observational tools, there are common systematic effects in mapping our local environment. That needs to be factored into any cosmological reconstruction <ref:2605.23802#pg1>.

Vera: So what does this mean for us listening? It means we have a very robust, self-consistent reference set of parameters for open clusters <ref:2605.23802#pg1>.

Jocelyn: It’s like having a high-quality map of the Milky Way’s history built from the stars themselves. You can use it to track how chemical enrichment evolves over time <ref:2605.23802#pg1>.

Subrahmanyan: And that opens up avenues for testing our theories on galactic evolution, which is where you need those kinds of detailed constraints <ref:2605.23802#pg1>.

Conclusion: Vera: So we’ve just looked at "Astrophysical Parameters of five thousand fifty-six Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones," and basically, we’ve got a super complete set of fundamental parameters for five thousand fifty-six open star clusters using Gaia data and PARSEC models.

Jocelyn: Right, it’s not just about getting numbers like age or metallicity; it’s how they use that full sample of five thousand fifty-six clusters to build a real baseline for the Milky Way.

Subrahmanyan: And from a theoretical standpoint, this gives us something concrete to test our models on, you know? It helps constrain things like the disc structure and chemical enrichment history.

Vera: Exactly. The methodology is pretty solid because they’re using Bayesian Nested Sampling with PARSEC isochrones applied uniformly across all the data. That’s a big win for consistency.

Jocelyn: And they treat metallicity as a free parameter initially, which means they get this photometric

Fe/H: posterior for every cluster without pre-supposing any abundance model.

Subrahmanyan: That freedom is important because it lets the data drive the result rather than forcing it through an external assumption. It’s a way to see what the clusters actually tell us about their own chemistry.

Vera: And they do have some smart ways to check their work, too. They define this inlier fraction, eta fit, which has to be greater than zero point eighty for a cluster’s fit to be considered high quality.

Jocelyn: So if the data doesn't match the isochrone well enough—if it falls outside that plus or minus three sigma band—the parameters derived should probably be treated with caution.

Subrahmanyan: That filtering step is crucial for reliability. It prevents noisy fits from polluting our bigger picture calculations about structure formation in the disc.

Vera: They also compare their results against other catalogues and found a consistent distance offset across three different photometric sets, which gives us a real systematic error to account for.

Jocelyn: That’s interesting because it means we can’t just assume one distance measurement is perfect; there's this shared bias across the literature.

Subrahmanyan: It shows that even when using different observational tools, there are common systematic effects in mapping our local environment. That needs to be factored into any cosmological reconstruction.

Vera: So what does this mean for us listening? It means we have a very robust, self-consistent reference set of parameters for open clusters.

Jocelyn: It’s like having a high-quality map of the Milky Way’s history built from the stars themselves. You can use it to track how chemical enrichment evolves over time.

Subrahmanyan: And that opens up avenues for testing our theories on galactic evolution, which is where you need those kinds of detailed constraints.

Vera: So we’ve covered "Astrophysical Parameters of five thousand fifty-six Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones," and it gives us a solid foundation.

Jocelyn: It really lays out how important it is to use these large samples to get those tight constraints on fundamental properties like age and metallicity.

Subrahmanyan: Definitely. Next time, we might look at how this type of stellar information connects to the larger cosmological models we’re working on.

Olcay Plevne, Furkan Akbaba

Faculty of Science, Department of Astronomy and Space Sciences, Istanbul University · Institute of Graduate Studies in Science, Istanbul University

astro-ph.GA

Submitted: 2026-05-22

Updated: 2026-10-13

Comments: 23 pages, 11 figures. Accepted for publication in MNRAS. Posterior chains: https://doi.org/10.5281/zenodo.23349618 ; code: https://github.com/oplevne/ocnsfit

Code: https://github.com/oplevne/ocnsfit

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 89/100

The gist: The gist: We present a homogeneous catalogue of fundamental astrophysical parameters — age, metallicity ([Fe/H]), heliocentric distance, and colour excess E(GBP−GRP) — for 5 056 open star

Key concepts

Bayesian Nested Sampling
A statistical technique used to explore the complex parameter space of stellar cluster properties. It efficiently calculates the probability distributions (posteriors) of parameters like age and metallicity by systematically sampling the likelihood function, helping determine which physical values are most likely given the observational data.
PARSEC Stellar Isochrones
These are theoretical models that predict how stars evolve and appear on a color-magnitude diagram based on specific physical inputs. The PARSEC models provide the necessary theoretical framework to compare actual star cluster observations (color-magnitude diagrams) with predicted stellar evolution tracks.
Inlier Fraction ($\eta_{fit}$)
This metric assesses how well the best-fitting theoretical model (isochrone) matches the actual observed stars in a cluster. It calculates what percentage of probable member stars fall within a tight $\pm3\sigma$ band around the isochrone, providing a physically interpretable measure of fit quality.
Photometric [Fe/H] Posterior
This is the derived range of initial stellar metallicity ([Fe/H]) calculated for each cluster. Because the method treats metallicity as a free parameter, it yields a photometric estimate directly from the data, allowing comparison with spectroscopic surveys without needing an assumed abundance model.

Terminology

Summary

The gist: We present a homogeneous catalogue of fundamental astrophysical parameters — age, metallicity ([Fe/H]), heliocentric distance, and colour excess E(GBP−GRP) — for 5 056 open star clusters drawn from the Unified Cluster Catalogue (UCC) using Bayesian Nested Sampling with PARSEC stellar isochrones.

How it works

The analysis employs a Bayesian inference framework that compares extinction-corrected colour–magnitude diagrams (CMDs) of probable member stars with theoretical stellar isochrones, applied uniformly across the full sample with no manual intervention on individual clusters The model has four free parameters: logarithmic age θ1 = log(Age/yr), initial metallicity θ2 = Zini, heliocentric distance θ3 = d (pc), and reddening θ4 = E(GBP − GRP)

The input data includes the Unified Cluster Catalogue (UCC; G. I. Perren et al. 2023) which consolidates 16 588 Galactic stellar clusters into a single homogeneous database The Gaia DR3 photometry (Gaia Collaboration et al. 2021) provides the G, GBP, and GRP magnitudes used to construct the colour–magnitude diagrams fitted in this work

How it works

The pipeline treats initial metallicity Zini as a fully free parameter across the entire grid, yielding a photometric [Fe/H] posterior for every cluster that is not anchored to an assumed abundance model Where Gaia XP spectrophotometric metallicities are available (N = 1 346 clusters), we use them as an informative starting point for the prior rather than a hard constraint For the remaining clusters a broad, environment-dependent uniform prior is adopted

How it works

The distance prior is a truncated Gaussian centred on the distance inferred from the median Gaia DR3 parallax of probable members The uncertainty on the median parallax is propagated as a truncated Gaussian with a lower bound at zero distance

The reddening prior is a truncated Gaussian, p(E(GBP − GRP)) ∝ N (µE, σ squared E) (E(GBP − GRP) ≥ 0) The prior mean is calculated using the D. J. Schlegel et al. (1998) full-sky dust map recalibrated by E. F. Schlafly & D. P. Finkbeiner (2011)

How it works

The likelihood function models the CMD residuals as a weighted chi-squared statistic, where each probable member star i is assigned the perpendicular distance ∆mi to the nearest point on the dereddened, distance-shifted isochrone The log-likelihood is lnL(θ) = −1/2 X i wi / (∆mi σ i 2)

How it works

Posterior distributions are sampled with the dynesty package, which implements the nested sampling algorithm of J. Skilling (2004) We use Nlive = 400 live points and the MultiEllipsoid bounding method Sampling terminates when the estimated remaining log-evidence falls below ∆ lnZ = 0.01

How it works

A central motivation is to address a gap in existing large catalogues: T. Cantat-Gaudin et al. (2020) fix metallicity to the solar value; W. S. Dias et al. (2021) include metallicity but adopt a position-dependent abundance–radius relation as the metallicity prior Our pipeline treats initial metallicity Zini as a fully free parameter across the entire grid, yielding a photometric [Fe/H] posterior for every cluster that is not anchored to an assumed abundance model

How it works

To assess the goodness of fit in a physically interpretable way, we define the inlier fraction ηfit: the fraction of probable members that fall within a ±3σ band around the best-fitting isochrone A cluster is classified as a high-quality fit if ηfit ≥ 0.80

How it works

The resulting catalogue spans ages from ∼3 Myr to ∼5.5 Gyr, distances from 88 to 19 011 pc, and metallicities from −1.17 to +0.42 dex

How it works

The median [Fe/H] of +0.002 dex is indistinguishable from solar, consistent with a sample dominated by disc clusters at Galactocentric radii near the solar circle where the disc metallicity is close to the solar value

How it works

The most stringent test of our metallicity scale is provided by the OCCAM DR19 catalogue (J. M. Otto et al. 2026) which lists high-resolution Milky Way Mapper (MVM: S. M´esz´aros et al. 2025) spectroscopic [Fe/H] measurements for open clusters in the MW disc We find 127 clusters in common after applying the ηfit ≥ 0.80 quality cut

How it works

The catalogue provides the largest homogeneous open-cluster parameter set to date derived with fully free photometric metallicity — 5 056 clusters processed with a single, uniform Bayesian pipeline on Gaia DR3 data This constitutes a self-consistent reference baseline for a broad range of Galactic science: reconstruction of the age–metallicity and age–velocity-dispersion relations of the disc; chemical tagging and identification of co-natal stellar populations; mapping of the three-dimensional dust distribution; and calibration of stellar evolution models

How it works

The pipeline source code will be released on GitHub with a citable DOI The complete set of dynesty posterior samples in HDF5 format (one file per cluster, totalling ∼5 000 files) is archived separately on Zenodo

How it works

The catalogue contains 5 056 entries, one per cluster, and is made available in full through the CDS/VizieR service and as a supplementary machine-readable file accompanying this article

How it works

The mean offset of only +0.029 dex is well within the systematic uncertainty of isochrone-based methods (≤ 0.1 dex) and is consistent with zero at the ∼0.2σ level

The catalogue provides a homogeneous determination of fundamental astrophysical parameters — age, initial stellar metallicity (Zini, converted to [Fe/H]), heliocentric distance, and selective extinction E(GBP−GRP) — for 5 056 open star clusters drawn from the Unified Cluster Catalogue (UCC; G. I. Perren et al. 2023) These parameters are derived by Bayesian Nested Sampling (NS) with PARSEC isochrones fitted to Gaia DR3 colour–magnitude diagrams (CMDs) The principal findings are summarised below >

How it works

The catalogue provides the largest homogeneous open-cluster parameter set to date derived with fully free photometric metallicity — 5 056 clusters processed with a single, uniform Bayesian pipeline on Gaia DR3 data It constitutes a self-consistent reference baseline for a broad range of Galactic science: reconstruction of the age–metallicity and age–velocity-dispersion relations of the disc; chemical tagging and identification of co-natal stellar populations; mapping of the three-dimensional dust distribution; and calibration of stellar evolution models

The analysis successfully processed all 5 056 AA and AB quality clusters in the UCC (G. I. Perren et al.

Improvements for AI systems

  1. The improved AI system can perform model-independent metallicity estimation by utilizing the photometric [Fe/H] values are derived without assuming any external abundance model and being directly comparable with spectroscopic surveys.

  2. The system can provide robust uncertainty quantification for derived physical parameters by reporting the posterior median together with the 16th and 84th percentile credible intervals (±1σ equivalent).

  3. The system can rigorously assess the reliability of isochrone fits using a physically interpretable metric, specifically by calculating the inlier fraction ηfit: the fraction of probable members that fall within a ±3σ band around the best-fitting isochrone.

  4. The system can identify and flag unreliable fits based on quality criteria, as it allows for classifying clusters where the reported parameters should be treated with appropriate caution if they do not meet the ηfit ≥ 0.80 threshold.

  5. The system can compare its derived results against literature catalogues to quantify systematic offsets, such as finding that Across all three photometric catalogues a consistent positive distance offset is present: +178, +225, and +43 pc.

  6. The system can provide a measure of the impact of prior information on parameter precision by comparing posterior widths; for example, it can show that "Clusters with an informative Gaia XP spectrophotometric metallicity prior (nXP ≥ 3 member stars) have significantly narrower [Fe/H] posteriors (median ∆[Fe/H] = 0.42 dex) than clusters fitted with the broad uniform prior (median 0.65 dex)."

  7. The system can perform automated cross-matching with spectroscopic surveys to confirm photometric results, as demonstrated by finding that The mean offset of only +0.029 dex is well within the systematic uncertainty of isochrone-based methods (≤ 0.1 dex) when compared to OCCAM DR19.

Abstract

We present a homogeneous catalogue of ages, metallicities, distances, and reddenings for 5056 open clusters in the Unified Cluster Catalogue (UCC), derived uniformly from Gaia Data Release 3 (DR3) colour-magnitude diagrams via Bayesian nested sampling with PARSEC isochrones. Unlike previous large-scale catalogues, metallicity is fitted as a free parameter for every cluster rather than fixed or adopted from the literature: where available, Gaia XP spectrophotometric metallicities inform the fit (1201 clusters), while the remaining clusters are constrained by broad priors from Gaia parallaxes and the Schlegel-Finkbeiner-Davis dust map. No cluster is anchored to an external open-cluster parameter catalogue, and no per-cluster manual intervention is applied. A fit-quality metric, η fit, identifies 3807 clusters (75 per cent) with reliable colour-magnitude solutions, spanning ages of 0.004-5.9 Gyr and distances out to 17 kpc. We validate the resulting parameters against six literature catalogues and, more critically, against independent high-resolution spectroscopic metallicities from OCCAM DR19: for 132 clusters in common we find a median absolute deviation of only 0.054 dex (-0.021 plus or minus 0.132 dex mean offset), comparable to the intrinsic systematic uncertainty of isochrone-based metallicities. As a direct consequence of this fidelity, the radial metallicity gradient recovered from our XP-informed clusters (-0.064 dex kpc-1) steepens towards the spectroscopically calibrated value (-0.079 dex kpc-1) for the highest-quality subset - demonstrating that a purely photometric pipeline, built without reference to spectroscopy, recovers a physically meaningful abundance scale. The full catalogue and posterior chains are publicly released.

Sources

Related papers