Bridging the Gap between Intermediate and Massive Stars I: Validation of MESA against the State-of-the-Art Monash Stellar Evolution Program for a 2$M_{\odot }$ AGB Star

Abstract

One--dimensional stellar structure and evolution programs are built using different physical prescriptions and algorithms, which means there can be variations between models' predictions even when using identical input physics. This leads to questions about whether such deviations are physical or numerical; code validation studies are important and necessary tools for studying these questions. We provide the first direct comparison between the Monash stellar evolution program and MESA for a $2M_{\odot }$ model evolved from the zero-age main sequence to the tip of the thermally pulsing asymptotic giant branch. We compare the internal structure of the two models at six critical evolutionary points and find that they are in excellent agreement in characteristics like central temperature, central density and the temperature at the base of the convective envelope during the thermally pulsing asymptotic giant branch. The hydrogen-exhausted core mass between the models differs by less than 4.2% throughout the entire evolution, the final values vary only by 1.5%. Surface quantities such as luminosity and radius vary by less than 0.2% prior to the asymptotic giant branch. During thermal pulses, the difference extends to 3.4%, largely due to uncertainties in mixing and the treatment of atmospheric boundary conditions. Given that the veteran Monash code is closed source, the present work provides the first fully open-source computational analog. This increases accessibility to precision modeling on the asymptotic giant branch and lays the groundwork for higher-mass calculations that are performed with MESA but preserve the standards of the Monash code during the AGB.

Submitted, ApJ

References in corpus (27)

• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
• Nucleosynthesis and Remnants in Massive Stars of Solar Metallicity
• The Dawes Review 2: Nucleosynthesis and stellar yields of low and intermediate-mass single stars
• The nucleosynthesis of Al26 and Fe60 in solar metallicity stars extending in mass from 11 to 120 Msun: the hydrostatic and explosive contributions
• Stellar Models and Yields of Asymptotic Giant Branch Stars
• YREC: The Yale Rotating Stellar Evolution Code
• The updated BaSTI stellar evolution models and isochrones: II. alpha-enhanced calculations
• Reevaluating the Mass-Radius Relation for Low-Mass, Main Sequence Stars
• Mass loss on the red giant branch: the value and metallicity dependence of Reimers' η in globular clusters
• The Helium abundance and Delta Y / Delta Z in Lower Main Sequence stars
• Helium enrichment and Carbon-star Production in Metal-rich Populations
• Empirical Tests of Pre-Main-Sequence Stellar Evolution Models with Eclipsing Binaries
• Low temperature Rosseland opacities with varied abundances of carbon and nitrogen
• On the Numerical Treatment and Dependence of Thermohaline Mixing in Red Giants
• The updated BaSTI stellar evolution models and isochrones. III. White Dwarfs
• Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning
• The treatment of mixing in core helium burning models -- III. Suppressing core breathing pulses with a new constraint on overshoot
• Evolutionary modelling of subdwarf B stars using MESA with the predictive mixing and convective premixing schemes
• Diagnostics of Stellar Modelling from Spectroscopy and Photometry of Globular Clusters
• Explaining the differences in massive star models from various simulations
• A phenomenological modification of thermohaline mixing in globular cluster red giants
• Schwarzschild and Ledoux are equivalent on evolutionary timescales
• The most metal-rich stars in the universe: chemical contributions of low and intermediate mass asymptotic giant branch stars with metallicities between $0.04\le Z\le \mathrm{0.10.}$
• The most metal-rich asymptotic giant branch stars
• Atomic diffusion in solar-like stars with MESA. Comparison with the Montreal/Montpellier and CESTAM stellar evolution codes
• The NuGrid AGB Evolution and Nucleosynthesis Data Set

Cited by in corpus (9)

• Modules for Experiments in Stellar Astrophysics (MESA): Time-Dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure
• Characterizing Observed Extra Mixing Trends in Red Giants using the Reduced Density Ratio from Thermohaline Models
• Modelling the evolution and influence of dust in cosmological simulations that include the cold phase of the interstellar medium
• Red giants evolutionary status determination: the complete Kepler catalog
• New Constraints on Sodium Production in Globular Clusters From the ${}^{23}$Na${(}^3$He$,\text{d}{)}^{24}$Mg Reaction
• Evidence for the Keplerian orbit of a close companion around a giant star
• Beyond MESA Defaults: The Impact of Structural Resolution Uncertainty in p-mode Asteroseismology
• White dwarf eccentricity fluctuation and dissipation by AGB convection
• Aluminium-26 production in low- and intermediate-mass binary systems