PARSEC v1.2S Isochrones: Stellar Evolution Models

• PARSEC v1.2S isochrones are a comprehensive grid of stellar evolution models covering a range of metallicities, masses, and evolutionary phases.
• They integrate advanced input physics including state-of-the-art opacities, nuclear reaction networks, and calibrated pre-main sequence treatments.
• Extensive empirical validation against clusters and eclipsing binaries ensures accurate synthetic CMDs for population synthesis and astrophysical analysis.

PARSEC v1.2S Isochrones are a foundational grid of stellar evolution models that deliver internally consistent evolutionary tracks and isochrone tables for a wide range of metallicities, initial masses, alpha-element enhancements, and evolutionary phases. Developed as an extension of the PAdova and TRieste Stellar Evolution Code (PARSEC), v1.2S incorporates significant advances in low-mass stellar physics, opacities, and boundary conditions and is calibrated to reproduce observed stellar populations in the Milky Way and nearby galaxies. These isochrones underpin synthetic color–magnitude diagrams, population synthesis, and precise inference of astrophysical parameters in studies ranging from resolved open clusters to galaxy-wide archaeology (Fu, 2017, Bressan et al., 2012, Brandner et al., 2023, Burgo et al., 2021, Wang et al., 2024, Marigo et al., 2017).

1. Input Physics, Mixtures, and Calibration

PARSEC v1.2S supersedes earlier Padova releases by introducing several advances in input microphysics:

• Equation of State: FreeEOS (Irwin), covering full ionization states and extending to very low temperatures (down to $T_\mathrm{eff} \sim 2500$ K), enabling robust modeling of low-mass, cool dwarfs (Fu, 2017).
• Opacities: High-temperature ($\log T > 4$) opacities from OPAL (Iglesias & Rogers), low-temperature opacities from AESOPUS (Marigo & Aringer) or Ferguson et al., with on-the-fly table generation for arbitrary metallicity, He- and $\alpha$-enhanced mixtures (Bressan et al., 2012).
• Nuclear Reaction Network: JINA Reaclib rates cover pp-chains, CNO cycles, Ne–Na and Mg–Al cycles, $\alpha$-captures, and light elements (Li, Be, B), enabling detailed pre-main sequence Li-depletion studies (Fu, 2017).
• Convection and Overshooting: Mixing-length theory, calibrated to solar values ($\alpha_\mathrm{MLT} = 1.74$), and diffusive overshooting treated with an exponentially decaying diffusion coefficient:

$$d_{\rm ov} = f_{\rm ov} H_P \, ,\quad D_{\rm ov}(r) = D_0 \exp\left[-2\frac{(r - r_0)}{f_{\rm ov} H_P}\right]$$

with $f_{\rm ov}$ mass-dependent (e.g., $\simeq 0.02$ for $1.1$–$1.7\,M_\odot$), smoothly vanishing for $M \lesssim 1.1\,M_\odot$ (Fu, 2017, Bressan et al., 2012).

• Pre-main Sequence Accretion and Photoevaporation: Incorporates residual disk accretion and EUV-driven outflows to address early Li evolution:

$$\dot{M}_\mathrm{acc}(t) = \dot{M}_\mathrm{ref}(t/t_\mathrm{ref})^{-\eta} \,,\;\; \dot{M}_\mathrm{pe} \simeq 4 \times 10^{-10} \left(\frac{\Phi_\mathrm{EUV}}{10^{41}\,\mathrm{s}^{-1}}\right)^{1/2} M_\odot\,\mathrm{yr}^{-1}$$

Isochrones are computed on a grid of metallicity ($Z$) and helium ($Y$):

• $Y = Y_p + (\Delta Y/\Delta Z) Z$\,, with $Y_p = 0.2485$ and $\Delta Y/\Delta Z = 1.40$ (or $1.78$ in some implementations), e.g., $Y \approx 0.273$ at solar $Z = 0.0152$ (Fu, 2017, Brandner et al., 2023).
• Three choices of $\alpha$-enhancement: $[\alpha/\mathrm{Fe}] = 0.0, 0.2, 0.4$ in v1.2S (Fu, 2017).
• Mass sampling from $0.1\,M_\odot$ to $15\,M_\odot$ (standard grid), up to $350\,M_\odot$ for massive-star extensions (Burgo et al., 2021).

Key calibration constants, such as $\alpha_\mathrm{MLT} = 1.74$, are set by matching the solar model, and overshooting parameters are tuned to open cluster and globular cluster color–magnitude diagrams.

2. Evolutionary Track Computation and Isochrone Extraction

Evolutionary tracks are computed from the pre-main sequence (Hayashi branch or ZAMS for $M \gtrsim 3\,M_\odot$) up to the onset of core carbon burning or helium flash:

• Low-mass ($0.1$–$2.0\,M_\odot$): Fine mass resolution ($\Delta M = 0.05\,M_\odot$) and detailed pre-main sequence evolution, with particular attention to PMS Li burning (Fu, 2017).
• Higher-mass stars: Coarser sampling, but consistently evolved up to post-MS phases.
• Advanced phases: Thermally-pulsing AGB (TP-AGB) treated with the COLIBRI code, appending detailed TPCs (Thermal Pulse Cycles) and surface abundance evolution (Marigo et al., 2017).

Isochrones are constructed by interpolating between equivalent evolutionary points (EEP) across tracks:

1. Select two tracks bracketing age $t$ at fixed EEP.
2. Interpolate observable (e.g., $L$, $T_\mathrm{eff}$, $g$) along each track to age $t$.
3. Interpolate between bracketing masses to yield desired isochrone points.
4. Repeat across desired age, metallicity, and $\alpha$-enhancement grid.

Isochrones cover ages $0.01$–$20$ Gyr with $\Delta\log t = 0.01$ dex, spanning pre-MS to TP-AGB, but not including white dwarfs (Fu, 2017, Marigo et al., 2017).

3. Output Structure, Bolometric Corrections, and Photometric Systems

Each v1.2S isochrone is provided as an ASCII table, with columns including:

• Age (Gyr)
• Initial mass $M_\mathrm{i}$; current mass $M$
• $\log L/L_\odot$, $\log T_\mathrm{eff}$, $\log g$
• Surface abundances: $X$, $Y$, $Z$, C, N, O (COLIBRI extension: individual CNO abundances through TP-AGB cycles)
• Bolometric corrections (BCs) for multiple systems: Johnson–Cousins UBVRI, 2MASS JHK, SDSS ugriz, Gaia $G_\mathrm{BP}$, $G$, $G_\mathrm{RP}$, and others

Model atmospheres for BCs:

• $T_\mathrm{eff} > 6500$ K: ATLAS9 (Castelli & Kurucz)
• $T_\mathrm{eff} < 5500$ K: PHOENIX BT-Settl (Allard et al.)
• $4000 < T_\mathrm{eff} < 5000$ K: Interpolation between ATLAS9 and PHOENIX (Wang et al., 2024)

TP-AGB phases incorporate dust radiative transfer, LPV pulsation periods, and synthetic photometry in the NIR/MIR.

4. Empirical Performance and Validation

PARSEC v1.2S isochrones have been extensively validated against key empirical datasets:

• Hyades open cluster: Fit to 600 single stars with $0.12 < M/M_\odot < 2.2$, achieving a median color-magnitude residual of $\leq 0.03$ mag for $0.4 < M/M_\odot < 2.0$ (Brandner et al., 2023).
• Detached eclipsing binaries: Main-sequence and core-helium burning masses inferred to within $4\%$ and $7\%$; subgiant and red giant masses underestimated on average by $12\%$ and $19\%$, respectively, indicating systematic biases in advanced phases (Burgo et al., 2021).
• Color-magnitude offsets: For low-mass stars ($BP-RP \gtrsim 1.2$), systematic color deviations occur in Gaia photometry: isochrones are $\sim$0.15 mag too blue in $BP-RP$ around $BP-RP \sim 2$–$3$ (Wang et al., 2024).
• TP-AGB benchmarking: Extended isochrones reproduce the observed number and properties of M- and C-rich TP-AGB stars, mean photometric trends, and the AGB-boosting effect near $\sim$1.6 Gyr (Marigo et al., 2017).

Limitations include lack of rotation (affecting turn-off morphology), possible He-enrichment deviations in specific clusters, and uncertain boundary conditions for the lowest-mass dwarfs.

5. Empirical Color Corrections and Practical Workflow

To address cluster-CMD discrepancies at low mass, empirical color correction functions have been introduced for Gaia bands (Wang et al., 2024). These take the form:

$$\Delta(\mathrm{BP} - \mathrm{RP})(x) = -0.089913\,x^9 + 0.75746\,x^8 - 3.4151\,x^7 + \ldots - 0.60593$$

$$\Delta(\mathrm{G} - \mathrm{RP})(x) = 1.1924 \times 10^3\,x^8 - 3.9313 \times 10^3\,x^7 + \ldots + 7.9764$$

with $x$ being the raw model color. Application workflow:

1. Download the uncorrected PARSEC v1.2S isochrone for desired age and $[\mathrm{M}/\mathrm{H}]$.
2. For each isochrone point, compute $\Delta$ from the above functions using its model color.
3. Add $\Delta$ to the respective color to obtain corrected $$(\mathrm{BP}-\mathrm{RP}), (\mathrm{G}-\mathrm{RP})$$.
4. Use the corrected isochrone for CMD fitting and parameter inference.

Corrected isochrones yield age estimates in close agreement with lithium depletion boundary measurements across a range of open clusters, eliminating systematic biases in low-mass age dating. The correction is restricted to $M \lesssim 0.8\,M_\odot$, $BP-RP \gtrsim 1.2$, and should not be extrapolated beyond these limits or to different metallicity regimes (Wang et al., 2024).

6. TP-AGB Extension: PARSEC v1.2S+COLIBRI

PARSEC v1.2S+COLIBRI isochrones comprise a detailed TP-AGB phase:

• Complete TPC reconstruction: Interpolate full COLIBRI tracks (quiescent and pulse phases) onto EEPs, reconstructing $L(\phi)$ and $T_\mathrm{eff}(\phi)$ variation within thermal pulse cycles (Marigo et al., 2017).
• Surface abundances: Tabulated for H, He, C, N, O including effects of diffusion, dredge-up events, and HBB.
• Dust and long-period variability: Self-consistent dust growth, NIR–MIR dust reprocessing, and variable star periods are included by combining COLIBRI and radiative transfer solutions.
• Output columns: Initial and current mass, $\log L$, $\log T_\mathrm{eff}$, $\log g$, core/envelope mass, mass fractions (X, Y, Z, C, N, O), C/O ratio, mass-loss rate, dust optical depth, pulsation periods, and evolutionary stage flags.
• User access: Available from http://stev.oapd.inaf.it/cmd and http://starkey.astro.unipd.it/cgi-bin/cmd with flexible photometric systems and time resolution.

This extension makes v1.2S+COLIBRI an end-to-end framework for population synthesis and modeling of resolved stellar populations, particularly in the near- and mid-IR, capable of supporting studies of chemically complex and dust-rich evolutionary phases (Marigo et al., 2017).

7. Implementation and Usage Recommendations

Retrieval of isochrones is performed via the CMD web interface (http://stev.oapd.inaf.it/cgi-bin/cmd):

• Specify model version (v1.2S), metallicity ($Z$ or $[\mathrm{M}/\mathrm{H}]$), age, $\alpha$-enhancement, photometric system, and advanced options (e.g., TP-AGB resolution for COLIBRI).
• Results are returned as ASCII files parsed by standard tools (e.g., Python, IDL, Fortran). A "parsec_interpolator" code (Python/C) enables fast multi-grid interpolation (Fu, 2017).
• Direct application to population synthesis, synthetic CMD generation, cluster sequence fitting, and parameter inference.
• Caveats: Do not extrapolate beyond published grid limits (mass, $Z$, $Y$, $\alpha$), nor apply empirical color corrections to uncalibrated domains. For stellar populations $t \lesssim 10$ Myr, use with caution due to evolving mass loss and accretion physics.
• Bayesian inference frameworks using PARSEC v1.2S benefit from precise main-sequence and core-He-burning parameter recovery, with uncertainties quantified for each evolutionary phase (Burgo et al., 2021). Polynomial fits to the zero-age main sequence are available for forward-modeling or analytic approximation.

PARSEC v1.2S thus establishes a rigorous and flexible isochrone suite for stellar and galactic astrophysics, with robust calibration, detailed microphysics, and broad empirical validation.