---
title: MESA Isochrones and Stellar Tracks (MIST)
url: https://www.emergentmind.com/topics/mesa-isochrones-and-stellar-tracks-mist
type: topic
---

# MESA Isochrones and Stellar Tracks (MIST)

Searching arXiv for recent and foundational MIST papers to ground the article.
MESA Isochrones and Stellar Tracks (MIST) is a stellar-model library and interpolation framework built on the one-dimensional stellar evolution code MESA, with the aim of providing self-consistent evolutionary tracks and isochrones across broad ranges of mass, age, composition, and evolutionary phase. In its solar-scaled release, MIST covers $0.1 \leq M/M_\odot \leq 300$, $5 \leq \log(\mathrm{Age})\,[\mathrm{yr}] \leq 10.3$, and $-2.0 \leq [Z/H] \leq 0.5$, with an extension to lower metallicity for models evolved from the pre-main sequence to the end of core helium burning; later releases add self-consistent $\alpha$-enhanced compositions and an explicit white-dwarf cooling-sequence extension [1604.08592][2602.22012][2509.21717]. Methodologically, MIST is closely identified with the Equivalent Evolutionary Point (EEP) formalism for constructing isochrones from heterogeneous stellar tracks, allowing a single interpolation framework to span phases from the pre-main sequence to the white-dwarf cooling sequence [1601.05144].

## 1. Historical development and project architecture

MIST emerged as a sequence of closely related releases rather than a single static grid. Paper 0 formalized the EEP-based interpolation machinery used to transform raw stellar tracks onto a uniform basis for track interpolation and isochrone construction [1601.05144]. Paper I presented the solar-scaled model library computed with MESA, including both theoretical and observational isochrones and extensive comparisons to other model sets and empirical constraints [1604.08592]. Subsequent work used the MIST framework to analyze cluster ages in the Gaia era, emphasizing how different stellar-model parameters project onto different observables and how multi-observable fitting can break classical degeneracies [1807.03789]. The project was then extended to $\alpha$-enhanced compositions over a two-dimensional abundance grid in $[\mathrm{Fe/H}]$ and $[\alpha/\mathrm{Fe}]$ [2602.22012], and later to carbon–oxygen-core white dwarfs with hydrogen atmospheres descended from full progenitor calculations [2509.21717].

The underlying computational engine is MESA (Modules for Experiments in Stellar Astrophysics), which solves the coupled stellar-structure and composition equations from the pre-main sequence through advanced phases. In the cluster-age analysis, MIST is described as using non-rotating tracks for the cluster fits, with a rotating grid published separately; the solar-scaled release itself includes both non-rotating and rotating families [1807.03789][1604.08592]. The project therefore combines three layers: interior stellar-evolution calculations, EEP-based interpolation, and atmosphere/bolometric-correction machinery for observational predictions.

A central organizational feature is the distinction between model families defined by composition and physics assumptions. The solar-scaled grids adopt solar-scaled heavy-element mixtures and a helium-enrichment law. The later $\alpha$-enhanced database adds explicit abundance variation in O, Ne, Mg, Si, S, Ar, Ca, and Ti, with microphysics tables and synthetic spectra recomputed for each abundance pattern [2602.22012]. The white-dwarf extension adds a further endpoint-specific branch to the library rather than a separate stand-alone product [2509.21717].

## 2. Stellar-physics content

The solar-scaled MIST release adopts a composite equation of state joining OPAL, SCVH, MacDonald, HELM, and PC tables across different thermodynamic regimes, together with OPAL Type I and Type II high-temperature opacities, Ferguson low-temperature opacities, and Cassisi electron-conduction opacities [1604.08592]. Its default nuclear network is `mesa_49.net`, containing pp-chains, cold and hot CNO cycles, triple-$\alpha$, $\alpha$-captures up to $^{32}\mathrm{S}$, Ne-Na, Mg-Al, and C/O burning, with rates drawn from JINA REACLIB [1604.08592].

Convection is treated with Henyey MLT under the Ledoux criterion. In MIST I, the mixing length is
$$
l_{\rm MLT}=\alpha_{\rm MLT}H_P,\qquad \alpha_{\rm MLT}=1.82,
$$
and diffusive exponential overshoot is written
$$
D_{\rm ov}(z)=D_0\exp\!\Bigl(-\frac{2z}{f_{\rm ov}H_P}\Bigr),
$$
with $f_{\rm ov,core}=0.016$ and $f_{\rm ov,env}=0.0174$ [1604.08592]. The same exponential-overshoot formalism appears in the cluster-age analysis as the parameterization used for H-core and He-core overshoot [1807.03789]. Semiconvection and thermohaline mixing are also included in the solar-scaled release, with explicit diffusion prescriptions and adopted efficiencies listed in the model description [1604.08592].

Mass loss is versioned by evolutionary regime. For low-mass stars, the solar-scaled release uses Reimers on the RGB and Bloecker on the AGB,
$$
\dot M_R=4\times10^{-13}\,\eta_R\,\frac{L/L_\odot}{M/M_\odot}\,\frac{R}{R_\odot},
$$
with $\eta_R=0.1$, and
$$
\dot M_B=4.83\times10^{-9}\,\eta_B\,(L/L_\odot)^{2.7}(M/M_\odot)^{-2.1}\,\dot M_R/\eta_R,
$$
with $\eta_B=0.2$; for high-mass stars it adopts the Dutch scheme, including Vink, Nugis–Lamers, and de Jager prescriptions, with a mass-loss cap of $10^{-3}\,M_\odot\,\mathrm{yr}^{-1}$ [1604.08592]. The cluster-age analysis likewise treats the Reimers efficiency as a parameter relevant to age inference, noting that it primarily affects the mass difference between RGB and red-clump stars while leaving CMD loci nearly unchanged [1807.03789].

Composition is parameterized through a helium-enrichment law. In MIST I,
$$
Y=Y_p+\frac{\Delta Y}{\Delta Z}Z,\qquad Y_p=0.249,\qquad \frac{\Delta Y}{\Delta Z}=1.5,
$$
with $Z_\odot=0.0142$ and $X=1-Y-Z$ [1604.08592]. The cluster-age analysis presents the same linear form and notes that it ties the initial helium mass fraction to metallicity in the standard solar-scaled framework [1807.03789].

The $\alpha$-enhanced release modifies this baseline in several ways. It recomputes EOS, OPAL, Ferguson, and atmosphere boundary tables for each $[\alpha/\mathrm{Fe}]$, replaces the earlier atomic-diffusion turbulence prescription, adopts a hybrid atmosphere boundary condition with a $T$–$\tau$ integration at $\tau=10$ for $M<8\,M_\odot$, and uses a metallicity-dependent envelope overshoot parameter,
$$
f_{\rm ov,env}=0.016-0.027\,[\mathrm{Fe/H}],
$$
constrained to $0.01 \leq f_{\rm ov,env} \leq 0.08$ [2602.22012]. Its solar calibration on the GS98 mixture yields $X_0=0.7080$, $Y_0=0.2735$, $Z_0=0.0185$, and $\alpha_{\rm MLT}=2.03$ [2602.22012]. These release-to-release changes are part of the model content and should not be conflated with a single immutable “MIST physics” specification.

## 3. EEP formalism and isochrone construction

The methodological core of MIST is the EEP formalism introduced by Dotter. The basic reparameterization replaces direct interpolation in age or timestep with a mapping
$$
\phi:(M_{\rm init},t)\rightarrow \xi,
$$
where $\xi$ is an index labeling equivalent evolutionary phases across tracks of different initial mass [1601.05144]. This avoids phase mismatch when adjacent tracks have entered different evolutionary states.

Primary EEPs are physically defined landmarks such as the pre-main-sequence start, ZAMS, IAMS, TAMS, RGB tip, onset and termination of core-helium burning, TPAGB or carbon-burning onset, post-AGB crossing, and WD cooling sequence [1601.05144]. Between adjacent primary EEPs, fixed numbers of secondary EEPs are inserted using a positive-definite metric,
$$
D_{i+1}=D_i+\sqrt{\sum_j w_j(x_{j,i+1}-x_{j,i})^2},
$$
with coordinates such as $\log L$, $\log T_{\rm eff}$, or $\log t$ and tunable weights $w_j$ [1601.05144]. Typical tracks use roughly $1\,500$–$2\,000$ EEP points [1601.05144].

Isochrone construction proceeds in two stages. First, each raw track is converted to the EEP basis. Second, for a desired age, one solves for the initial mass corresponding to each EEP by interpolating the age–mass relation at fixed EEP, then reads off $L$, $T_{\rm eff}$, $R$, $\log g$, surface abundances, and derived photometric quantities [1601.05144]. The solar-scaled MIST paper summarizes the same logic as interpolation in the three-dimensional grid $(M_i,\mathrm{Age},[Z/H])$ on a uniform EEP mesh, with final isochrones containing more than 500 EEP points and more than 50 stellar properties [1604.08592].

A technically important feature is the treatment of non-monotonic age–mass behavior. Dotter identifies two physical transitions that can violate the usual monotonic assumption: the onset of a convective core near $\sim1.1\,M_\odot$ and the degenerate–nondegenerate He-ignition boundary near $\sim1.8\,M_\odot$ [1601.05144]. MIST handles incidental non-monotonicities by isotonic regression using the Pool Adjacent Violators algorithm and can alternatively export secondary branches when users wish to retain genuinely multi-valued isochrone structure [1601.05144].

The same EEP framework remains operative in later releases. The $\alpha$-enhanced models state explicitly that MESA tracks are converted to uniform EEP tracks via `iso`, after which EEP tracks are interpolated in mass for a given age to form isochrones [2602.22012]. The continuity of the interpolation formalism across versions is one of the project’s central design features.

## 4. Grid coverage and released model families

The solar-scaled MIST library spans initial masses from $0.1$ to $300\,M_\odot$ in roughly 100 models with finer spacing near critical masses, an age grid of $\log\mathrm{Age/yr}=5.0$ to $10.3$ in steps of $0.05$, and metallicities from $[Z/H]=-2.0$ to $+0.5$ in steps of $0.25$ dex, plus an extension from $[Z/H]=-4.0$ to $-2.0$ evolved from the pre-main sequence to the end of core-helium burning [1604.08592]. Evolutionary endpoints depend on mass: low-mass tracks may terminate at TAMS if they do not reach helium burning within the modeled timespan; intermediate-mass low-mass-type tracks are followed through TPAGB, post-AGB, and into WD cooling; high-mass tracks are followed through carbon burning [1604.08592].

Rotation is represented by two model families in the v1.2-era description: $v/v_{\rm crit}=0.0$ and $v/v_{\rm crit}=0.4$ [2208.04969]. In MIST I, shellular rotation is imposed with solid-body rotation on the ZAMS and a ramp from $M_i=1.2$ to $1.8\,M_\odot$ up to $v_{\rm ZAMS}/v_{\rm crit}=0.4$, with diffusive transport coefficients following Heger et al. and rotationally enhanced mass loss following Langer [1604.08592]. The later $\alpha$-enhanced release again provides non-rotating and rotating families, with initial $\Omega/\Omega_c=0.4$ for $M>1.8\,M_\odot$, tapered to zero at $1.2\,M_\odot$, and adds a Tayler–Spruit dynamo prescription for $M<7\,M_\odot$ [2602.22012].

The $\alpha$-enhanced grid expands the composition space to 74 compositions: $[\mathrm{Fe/H}]$ from $-3.0$ to $+0.5$ in steps of $0.25$ dex and $[\alpha/\mathrm{Fe}] = -0.2, 0.0, +0.2, +0.4, +0.6$, with the $[\mathrm{Fe/H}]=+0.5$, $[\alpha/\mathrm{Fe}]=+0.6$ point omitted [2602.22012]. The abundance changes are applied self-consistently to both interior models and atmosphere/synthetic-spectrum grids. Bolometric-correction tables carry the same $[\alpha/\mathrm{Fe}]$ dimension, and the atmosphere grid spans $T_{\rm eff}=3500$–$15\,000\,\mathrm{K}$ and $\log g=-1$ to $5$ [2602.22012].

Observationally, MIST distributes both theoretical and photometric isochrones. The solar-scaled release provides `.iso` and `.bc.iso` files with magnitudes in approximately 30 photometric systems, including UBVRI, SDSS, 2MASS, Gaia, HST, JWST, Spitzer, and WISE [1604.08592]. The later $\alpha$-enhanced release preserves the general file structure while adding columns for gravity darkening, convective turnover times, and apsidal-motion $k_2$ [2602.22012].

## 5. Empirical performance, validation, and known discrepancies

Validation has been a defining component of MIST since the solar-scaled release. Paper I reports solar calibration, comparisons to open clusters including M 67, Praesepe, the Pleiades, NGC 6791, and Ruprecht 106, detached eclipsing binaries, asteroseismic diagnostics, massive-star population ratios, and other model sets such as PARSEC, $Y^2$, DSEP, BaSTI, and Lyon [1604.08592]. In the open-cluster comparison, the models reproduce the MS, MSTO, RGB, and red-clump loci in multiband CMDs, but also exhibit the well-known mismatch below $0.6$–$0.7\,M_\odot$, where the models are too blue in $BV$ and $VI$, attributed there to missing line opacities in bolometric corrections [1604.08592].

The cluster-age analysis sharpened this by examining how stellar age and model parameters such as $Y$, $\alpha_{\rm MLT}$, $f_{\rm ov}$, and $\eta_R$ imprint differently on CMD morphology, mass–radius relations, surface abundances, asteroseismic $\log g$, and population ratios [1807.03789]. For NGC 6819, MIST fits were reported as showing excellent agreement in the turn-off, subgiant, and clump; in M 67, the Henyey-hook morphology constrained core overshoot near $f_{\rm ov,H}\approx0.016$ but mild RGB-color and [C/N] tensions remained; in NGC 6791, no single fit reproduced all CMDs simultaneously, and Gaia DR2 parallax zero-point uncertainties were identified as a leading systematic [1807.03789].

The Hyades benchmark provides a more differential test of the single-star sequence. Using Gaia EDR3, non-rotating MESA models at $[\mathrm{Fe/H}] = +0.25$ fit stars above $0.85\,M_\odot$ and below $0.25\,M_\odot$ well, but systematically underpredict the luminosity of stars between $0.25$ and $0.85\,M_\odot$, by up to $\Delta G \simeq 0.2$ mag [2208.04969]. That study points to the fixed solar-calibrated mixing-length treatment as a potential limitation for partially convective stars and notes increased scatter near the fully convective boundary, possibly related to the convective-kissing instability driven by $^3$He burning [2208.04969].

A related empirical issue concerns Gaia colors. An analysis based on Hyades, Pleiades, and Praesepe found that current MIST isochrones are systematically too blue at low masses, with typical deviations of $\Delta_{BP-RP}\approx +0.20$–$0.25$ mag and $\Delta_{G-RP}\approx +0.05$–$0.08$ mag for $M_G>9$ [2411.12987]. Polynomial color-correction functions were then derived for non-rotating, solar-scaled MIST isochrones, reducing the RMS color difference from $\sim0.02$–$0.04$ mag to $\sim0.005$ mag and removing a systematic age bias of about $0.075$ dex in $\log(\mathrm{Age/yr})$ when the low-mass sequence is included [2411.12987]. The same work explicitly states that the underlying cause of the offset remains unclear, listing model-atmosphere opacities, low-mass convection, and Gaia photometric calibration as possibilities [2411.12987].

MIST is also used as a forward model in downstream Bayesian inference. BRUTUS embeds MIST v1.2 isochrones in a joint photometric and astrometric likelihood with Galactic priors on distance, metallicity, age, dust, and $R_V$, and applies empirical corrections to $\log T_{\rm eff}$ and $\log R$ on the lower main sequence together with photometric zero-point offsets derived from nearby field stars [2503.02227]. In that framework, the calibration is reported to reduce lower-main-sequence systematic errors to $<2\%$, and tests on mock and real data recover distances to $\lesssim10\%$, $T_{\rm eff}$ to $\lesssim100\,\mathrm{K}$, and $A_V$ to $\lesssim0.1$ mag when calibrated models and Gaia astrometry are used [2503.02227]. These results are not changes to MIST itself, but they document how MIST grids are operationally calibrated in precision inference workflows.

## 6. White-dwarf cooling sequence extension and current data ecosystem

The white-dwarf extension substantially broadens the evolutionary endpoint coverage of MIST. Previous versions ended at the end of the AGB, around $T_{\rm eff}=10^4\,\mathrm{K}$; the updated library extends tracks and isochrones down the white-dwarf cooling sequence to $T_{\rm eff}=2\,000\,\mathrm{K}$ and cooling ages beyond $14$ Gyr [2509.21717]. The released models are carbon–oxygen-core DA/DC white dwarfs descended from full progenitor calculations. For each metallicity point $([\mathrm{Fe/H}],Y)$ in the 84-point MIST grid, approximately 100 WD tracks are produced with $M_{\rm WD}\approx0.5$–$1.05\,M_\odot$, descended from progenitors of $0.6$–$6.5\,M_\odot$ [2509.21717].

The WD physics goes beyond a simple Mestel-like cooling treatment. Core compositions are inherited from full progenitor evolution, with typical central oxygen mass fractions $X_{\rm O,c}\simeq0.55$–$0.75$ and $X_{\rm Ne}\approx Z_{\rm init}$ from CNO $\rightarrow {}^{22}\mathrm{Ne}$ processing [2509.21717]. Envelope hydrogen and helium masses are computed self-consistently, with representative starting values around $M_{\rm H}\sim10^{-4}\,M_\odot$ for a $0.6\,M_\odot$ WD and $M_{\rm He}\sim10^{-2}\,M_\odot$ [2509.21717]. Time-dependent element diffusion is solved in Lagrangian form,
$$
\frac{\partial X_i}{\partial t}
=-\frac{1}{\rho r^2}\frac{\partial}{\partial r}\bigl(r^2\rho D_i \frac{\partial X_i}{\partial r}\bigr),
$$
and heavy-element sedimentation, residual nuclear burning, crystallization, latent heat, and carbon–oxygen phase separation are all included self-consistently [2509.21717].

The cooling problem is summarized by
$$
L=-\frac{dU}{dt}+\int(\epsilon_{\rm nuc}+\dot{\epsilon}_{\rm sed}+\Delta\epsilon_{\rm lat}+\Delta\epsilon_{\rm ps})\,dm,
$$
with residual burning capable of delaying cooling by up to about $1$ Gyr for low-$Z$, low-mass white dwarfs, and C/O phase separation delaying cooling by $0.5$–$1$ Gyr around $\log L/L_\odot\sim -3$ [2509.21717]. For practical use, cooling-age contour files are distributed in $(M_{\rm WD},\log L)$ space, with both raw and LOWESS-smoothed versions and a stated bilinear interpolation recipe for retrieving $t_{\rm cool}$ from measured mass and luminosity [2509.21717].

Validation of the WD extension is both internal and external. The resulting initial–final mass relation emerges self-consistently from the full progenitor calculations and is reported to reproduce recent empirical constraints; the cooling tracks show good agreement with LPCODE, BaSTI, and STELUM to approximately $2$ Gyr in cooling age, with later differences attributed to more modern treatments of crystallization and phase separation in MESA/Skye EOS [2509.21717]. A $0.6\,M_\odot$ WD track aligns precisely with the A-branch peak in the Gaia HRD, which is presented as confirmation that the adopted DA/DC bolometric corrections are well calibrated [2509.21717].

Current access reflects the layered history of the project. Solar-scaled tracks and isochrones were originally distributed through the MIST project website at `http://waps.cfa.harvard.edu/MIST/` [1604.08592]. The $\alpha$-enhanced release states that data products are available from `https://mist.science`, with a direct data portal at the earlier WAPS site and full reproducibility files on Zenodo [2602.22012]. The WD extension distributes model tracks, isochrones, and cooling-age contours via the MIST website and Zenodo DOI `10.5281/zenodo.15242047`, with full MESA inlists and run directories at DOI `10.5281/zenodo.15196934` [2509.21717]. Across these releases, MIST functions both as a published stellar-model database and as a reproducible framework for generating new tracks, bolometric-correction tables, and isochrones under controlled physics assumptions.

Source: https://www.emergentmind.com/topics/mesa-isochrones-and-stellar-tracks-mist