Globular Cluster CMD Explorer

A colour-magnitude diagram is the most-used observational tool in stellar astrophysics. Read the cluster's age from the main-sequence turnoff; see how metallicity shifts the whole sequence. Switch to multi-population mode to see why Omega Centauri's CMD doesn't fit a single isochrone β€” the fingerprint of its stripped-dwarf-galaxy origin.

πŸ”¬ Stellar evolution ⚠ OC multi-pop structure
This is a parametric isochrone, not a fit to data. The branches are smooth analytic approximations of MIST/PARSEC-style outputs β€” accurate enough that age, metallicity, distance, and reddening responses are correct in direction and within ~0.5 mag in magnitude, but not a substitute for a real stellar-evolution code if you need photometric-precision work. The shape of the OC multi-population fan, however, is qualitatively right and matches the published JWST + HST data.
Cluster preset:
Cluster parameters
11 Gyr
βˆ’0.75
13.5 mag
0.04
0.05
Display
Colour-magnitude diagram
main sequence subgiant + RGB horizontal branch AGB white dwarfs
MS turnoff M_V
β€”
at current age
Turnoff mass
β€”
Mβ˜‰
HB type
β€”
β€”
Implied age from turnoff
β€”
verification

Why colour-magnitude diagrams matter

A CMD plots each star's brightness (V-band magnitude) against its colour (Bβˆ’V), and stars at different evolutionary stages live in different regions. Reading a CMD is the bread and butter of stellar astrophysics: the position of the main-sequence turnoff fixes the cluster age (more massive stars exhaust their hydrogen first, so the turnoff moves down-and-redward over time); the shape of the red-giant branch encodes metallicity; the colour of the horizontal branch (the post-helium-flash core-burning phase) is jointly set by metallicity, age, and mass loss. A two-parameter (age, metallicity) isochrone fitted to a single-population cluster's CMD typically nails the age to ~10% and the metallicity to ~0.1 dex.

How this tool's isochrones are computed

Each branch (MS, MSTO, SGB, RGB, HB, AGB, WD) is a smooth analytic curve in (Bβˆ’V, M_V) space, parameterised by age (which sets the turnoff mass via the simplified mass-lifetime relation M_turn β‰ˆ (10/Ο„_Gyr)^0.5 Mβ˜‰ and the turnoff luminosity via L ∝ M^3.5) and [Fe/H] (which shifts colours redward as metallicity rises and rotates the HB morphology from blue at [Fe/H] < βˆ’1.5 to red-clump at [Fe/H] > βˆ’0.7). Reddening shifts colour by +E(Bβˆ’V) and dims V by +3.1 E(Bβˆ’V) per the standard R_V = 3.1 law. Distance modulus is a pure vertical shift.

This parametric approach is intentional: it works under file:// with no external data, captures the physics correctly enough for teaching, and lets you isolate which parameter affects which branch. For research-grade fitting, use PARSEC (Bressan et al. 2012), MIST (Choi et al. 2016), or BaSTI (Pietrinferni et al. 2004).

Omega Centauri's multi-population fingerprint

Most globular clusters fit a single isochrone reasonably well β€” they formed in a single starburst from chemically uniform gas. Omega Centauri does not. The classic Lee, Demarque & Zinn (1994) and later HST work (Bedin et al. 2004) found that the cluster's main sequence is split into multiple parallel tracks at different colours. The two- or three-population analyses of the past decade (oMEGACat VI / Sommer et al. 2025; Zhuang et al. 2025 JWST+HST; A&A 695:A12 2025 Gaia FPR + HST kinematics) decompose this into populations at [Fe/H] β‰ˆ βˆ’1.8, βˆ’1.4 and βˆ’0.9 with progressively enhanced helium fractions. The metal-rich population sits as a redder, brighter MS off to the right of the dominant metal-poor sequence. This is the colour-magnitude evidence that OC is not a simple cluster. Multiple distinct populations means multiple star-forming episodes or merger, both of which point to OC being the stripped core of a former dwarf galaxy rather than a coeval globular.

Why this matters for the IMBH question

The IMBH hypothesis (see Tool 2) relies on a well-defined kinematic centre and a believable enclosed-mass profile. If OC has three populations with distinct kinematics β€” and the A&A 695:A12 (2025) paper using Gaia FPR + HST shows it does β€” then mass profiles fitted to "all stars" mix together populations with different orbital structure. The HΓ€berle (2024) fast-star analysis already had to argue that the seven fast stars belong to a kinematically coherent sample. Future IMBH constraints will increasingly need to be done per population, with this CMD as the population assignment tool. The CMD is upstream of every IMBH measurement.

Classroom uses

HS / intro-astro: set up a young cluster (Pleiades preset) and an old cluster (M67), and ask students to estimate the age difference from the turnoff alone. Undergrad: compare a metal-poor globular (M92) to a metal-rich one (47 Tuc) at the same age; have students explain the HB morphology difference and the RGB shift. Graduate: turn on multi-population mode for Omega Centauri and decompose the synthetic data into individual [Fe/H] components; discuss what additional observations (high-resolution spectroscopy, asteroseismology) would tighten the decomposition.

Real-world reference numbers (as of May 2026)

πŸ“Έ The HST UV-Legacy Survey of GCs

The HST UV-Legacy Survey (Piotto et al. 2015 AJ 149:91, Nardiello et al. 2018 MNRAS 481:3382) produced ultra-deep multi-band photometry of 56 Galactic globular clusters using HST/WFC3-UVIS. For Omega Centauri specifically, the survey identified ~15 distinct stellar populations from the UV-blue colour distribution β€” the most spectacularly complex of any Milky Way GC, consistent with OC being the stripped nuclear cluster of an accreted dwarf galaxy. The oMEGACat collaboration (Sommer et al. 2025) extends this to ~1.4 million stars with photometry + 3D kinematics.

🌈 Multiple populations β€” the consensus

By 2025 essentially every well-studied Milky Way GC shows at least two distinct chemical populations (first-generation "P1" with primordial Big-Bang abundance pattern; second-generation "P2" enhanced in He, Na, Al and depleted in O, Mg). OC has at least 3–5 populations and possibly more. Gratton, Carretta & Bragaglia 2012 (A&ARv 20:50) is the canonical review. The mechanism that produces P2 stars (AGB pollution? Fast-rotating massive stars? Supermassive stars?) remains debated.

πŸ“ Isochrone fitting precision

Modern isochrone-fitting (PARSEC, MIST, BaSTI sets) achieves age precision of ~0.5 Gyr for the oldest GCs at ~12-13 Gyr β€” close to the age of the universe (13.8 Gyr). For OC the best-fit age is 11.5 Β± 0.5 Gyr (VandenBerg et al. 2013). The distance modulus is also fitted from the CMD: OC distance β‰ˆ 5.43 Β± 0.05 kpc (Soltis, Casertano & Riess 2021, ApJL 908:L5). These precision values are the foundation of all derived quantities β€” IMBH mass, structural parameters, dynamical models all flow from them.

v1.0 β€” 2026-05-16 Β· Code MIT Β· Prose CC BY 4.0 Β· Sommer et al. 2025 (oMEGACat VI, ApJ); Zhuang et al. 2025 (arXiv:2506.21187); A&A 695:A12 2025; Bressan et al. 2012 (PARSEC); Choi et al. 2016 (MIST)