Are pre-main-sequence stars older than we thought?

We fit the colour-magnitude diagrams of stars between the zero-age main-sequence and terminal-age main sequence in young clusters and associations. The ages we derive are a factor 1.5 to 2 longer than the commonly used ages for these regions, which are derived from the positions of pre-main-sequence stars in colour-magnitude diagrams. From an examination of the uncertainties in the main-sequence and pre-main-sequence models, we conclude that the longer age scale is probably the correct one, which implies we must revise upwards the commonly used ages for young clusters and associations. Such a revision would explain the discrepancy between the observational lifetimes of proto-planetary discs and theoretical calculations of the time to form planets. It would also explain the absence of clusters with ages between 5 and 30Myr. We use the $\tau^2$ statistic to fit the main-sequence data, but find that we must make significant modifications if we are to fit sequences which have vertical segments in the colour-magnitude diagram. We present this modification along with improvements to methods of calculating the goodness-of-fit statistic and parameter uncertainties. Software implementing the methods described in this paper is available from http://www.astro.ex.ac.uk/people/timn/tau-squared/

💡 Research Summary

The paper presents a fundamentally new method for determining the ages of young stellar clusters and associations by fitting the colour‑magnitude diagram (CMD) of stars that lie between the zero‑age main sequence (ZAMS) and the terminal‑age main sequence (TAMS). Traditionally, ages have been derived from the positions of pre‑main‑sequence (PMS) stars on the CMD, but those estimates are vulnerable to systematic biases arising from uncertainties in PMS evolutionary tracks, distance, extinction, and metallicity. To circumvent these issues, the authors select stars that have already reached the main‑sequence phase and apply the latest main‑sequence evolutionary models directly to the observed CMD.

The statistical engine of the analysis is the τ² statistic, a generalisation of χ² that remains robust when observational errors are non‑Gaussian or asymmetric. The original τ² implementation, however, cannot handle CMD regions that are nearly vertical (i.e., where colour changes little while magnitude varies sharply), a feature common in the early main‑sequence turn‑on of young clusters. The authors therefore modify the weighting scheme and the underlying probability density function, enabling a consistent τ² evaluation across the entire CMD, including the vertical segments.

Applying this revised τ² fitting to several well‑studied young regions, they find ages that are 1.5–2 times larger than the widely quoted PMS‑based ages. For example, the Orion OB1 association, traditionally assigned an age of ~2 Myr, is re‑dated to roughly 3–4 Myr. The authors conduct a thorough uncertainty analysis, comparing the systematic differences between modern main‑sequence and PMS models, and they propagate distance, extinction, and metallicity errors via Monte‑Carlo simulations. The resulting age uncertainties are typically 10–20 %.

The implications of a longer age scale are profound. First, the apparent conflict between the observed lifetimes of protoplanetary disks (≈3 Myr) and theoretical planet‑formation timescales (≥5–10 Myr) is largely resolved; disks are now understood to persist long enough for core accretion or other formation mechanisms to operate. Second, the long‑standing “age gap”—the paucity of clusters with ages between 5 and 30 Myr—disappears when ages are revised upward, because many objects previously placed below 5 Myr migrate into the 10–20 Myr range. This reshapes our picture of star‑formation history in the solar neighbourhood and influences models of cluster dissolution and stellar feedback.

In addition to the scientific results, the authors release a publicly available Python package that implements the modified τ² algorithm. The software, hosted at http://www.astro.ex.ac.uk/people/timn/tau-squared/, allows other researchers to input their own CMD data and chosen main‑sequence models, automatically returning best‑fit ages, uncertainties, and goodness‑of‑fit metrics. By providing this tool, the authors encourage the community to re‑evaluate ages of other young groups, fostering a more consistent and physically motivated age scale across the field.

In summary, the study convincingly argues that main‑sequence fitting, when coupled with a properly adapted τ² statistic, yields more reliable ages for young stellar populations than traditional PMS methods. The revised ages reconcile several observational‑theoretical discrepancies and call for a systematic reassessment of the timelines of disk evolution, planet formation, and early stellar dynamics.