Modules for Experiments in Stellar Astrophysics (MESA): Convective Boundaries, Element Diffusion, and Massive Star Explosions
read the original abstract
We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact.
This paper has not been read by Pith yet.
Forward citations
Cited by 60 Pith papers
-
Unveiling the nature of barium stars. I. Asteroseismic masses and the evolutionary link between Ba dwarfs and giants
Asteroseismic masses average 1.29 Msun for Ba dwarfs versus 1.96 Msun for Ba giants, supporting main-sequence accretion evolution from dwarfs to giants, though models fail to match the observed [hs/ls] ratio.
-
Can Asteroseismic Structure Inversions Be Performed in Structure-Dependent Coordinates?
Asteroseismic structure inversions can be performed in acoustic and buoyancy coordinates, overcoming breakdowns for non-solar stars and removing the need for prior mass and radius knowledge in p-mode cases.
-
Evolution of Stars During the Main Sequence and the Transition to the Red Giant Phase
Analytical relations for convective core mass, variable-mu structure, and main-sequence end in intermediate-mass stars, with MESA validation.
-
A systematic survey for hypervelocity runaways from thermonuclear supernovae
A systematic survey with complete spectroscopic classification finds 13 hypervelocity D6 and LP 40-365 stars and shows that intermediate-heating evolutionary models best reproduce the observed population at birth rate...
-
Revealing the Origin of Desert Dwellers via Stellar Obliquities
Roche lobe overflow during gas giant destruction aligns stellar spins with orbits within tens of degrees regardless of starting conditions, offering an observable to differentiate this mechanism from high-eccentricity...
-
Extending asteroseismic magnetometry across the diverse landscape of magnetic structures
Extends magnetogravity polarization formalism to arbitrary magnetic field geometries, revealing avoided crossings and mode conversion below a local field threshold.
-
Revealing mixed modes in compressible hydrodynamical simulations of red giant stars
2D hydro simulations of a red giant star find that low-frequency mixed modes below 50 microhertz carry the largest kinetic energies, unlike empirical predictions peaking near nu_max, with surface velocities matching i...
-
Precision constraints on stellar physics from main sequence detached eclipsing binaries
Bayesian framework with active-learning surrogate for MESA models constrains ages and α_MLT from 38 main-sequence DEBs, recovering some α_MLT values below the solar calibration.
-
Accretion of Primordial Black Holes in Stellar Interiors
Self-consistent spherical accretion simulations show cooling-enhanced growth of PBHs with radiative efficiency ~10^{-2} in the bremsstrahlung regime, yielding a critical seed mass of ~10^{-16} M_sun to consume a solar...
-
The Life and Death of Stars That Capture Primordial Black Holes
Primordial black holes captured by stars lead to either quiet consumption or explosive disruption via disk formation, producing transients and high-spin remnants with potentially observable event rates.
-
Hydrodynamic Response of Mildly Evolved Common Envelope Donors in Luminous Red Novae
3D hydro simulations demonstrate that donor central concentration (ρ_c/ρ_bar) controls inspiral morphology and mass-ejection history in common-envelope events for luminous red novae.
-
TESS Observations of Stochastic Low-frequency Variability in Extreme Helium Stars
TESS data show stochastic low-frequency variability dominates in most extreme helium stars, with characteristic timescales of 0.5-10 days correlating to stellar parameters and matching subsurface convection predictions.
-
Blue Straggler Stars in Old Open Clusters and the Kraft Break
Blue straggler stars in old open clusters exhibit a Kraft break in rotation, with rapid rotators above the break and slow rotators below, indicating their envelopes behave like those of single stars.
-
Binarity at LOw Metallicity (BLOeM): massive star variability revealed using a novel software tool for point-spread function fitting of TESS images
A new PSF-fitting tool extracts TESS light curves for 91 SMC massive stars, revealing binarity, pulsations, and SLF variability whose morphology tracks HR diagram position similarly to Galactic stars.
-
Binarity at LOw Metallicity (BLOeM): massive star variability revealed using a novel software tool for point-spread function fitting of TESS images
A PSF-fitting pipeline extracts variability light curves for 91 SMC massive stars and finds that stochastic low-frequency variability morphology tracks HR-diagram position independently of metallicity.
-
The Tale of a Hungry Subgiant and Its Brown Dwarf: Interior Radiative Damping Dominates the Tidal Evolution of TOI-5882
A coupled MESA-GYRE framework shows interior radiative damping of gravity waves dominates tidal evolution in TOI-5882, predicting a 2-6 fold reduction in engulfment timescale versus equilibrium tide models.
-
The Tale of a Hungry Subgiant and Its Brown Dwarf: Interior Radiative Damping Dominates the Tidal Evolution of TOI-5882
Radiative damping of internal gravity waves dominates tidal evolution in TOI-5882, shortening the brown dwarf's engulfment timescale by a factor of 2-6 relative to classical models.
-
Type Ib Supernovae are bluer than Type Ic Supernovae
Type Ib supernovae are systematically bluer than Type Ic supernovae in optical colors, likely due to helium-rich versus helium-poor progenitors.
-
Hyperaccreting Neutron Stars inside Massive Envelopes: The Implausibility of Thorne-\.Zytkow Objects
First coupled GRHD simulations with M1 neutrino transport and alpha-chain network show hyperaccreting NSs in envelopes collapse to BHs without forming stable TZOs or unbound ejecta.
-
Hyperaccreting Neutron Stars inside Massive Envelopes: The Implausibility of Thorne-\.Zytkow Objects
Hypercritical accretion onto neutron stars embedded in massive envelopes leads to rapid collapse into black holes rather than stable Thorne-Zytkow objects.
-
White dwarf + M dwarf Detached Binaries in Long Period Radio Transients: Observed Binary Parameters, Evolution, and Population Constraints
Two long period radio transients are detached white dwarf-M dwarf binaries with matching periods, massive cool crystallized white dwarfs, low inclinations, and an estimated population of 100-2000 such systems within 2 kpc.
-
Irregularly Sampled Time Series Interpolation for Binary Evolution Simulations Using Dynamic Time Warping
Dynamic Time Warping with a shared warping path across parameters aligns binary stellar tracks for accurate interpolation while preserving physical relationships such as the Stefan-Boltzmann law.
-
The Black Hole Mass Gap as a New Probe of Millicharged Particles
Millicharged particles weaken pulsational pair-instability in massive stars, shifting the lower edge of the black hole mass gap upward and turning gravitational wave observations into a probe for particles with masses...
-
pop-cosmos: Star formation over 12 Gyr from generative modelling of a deep infrared-selected galaxy catalogue
A score-based diffusion generative model on deep infrared galaxy photometry yields a star formation rate density peaking at z=1.3 and shows distinct non-parametric star formation histories plus AGN activity peaking du...
-
A grid of fast-rotating, chemically-homogeneous, supernova and/or long-GRB progenitors
Presents a grid of 113 fast-rotating, chemically-homogeneous massive star models at Z=0.001 reaching core collapse with high angular momentum for use as supernova and GRB progenitors.
-
Dynamical Tides during High-Eccentricity Migration produces the Hot Jupiter Pile-up, Neptune Ridge, and Neptune Desert
Dynamical tides exciting f-modes during high-eccentricity migration produce the hot Jupiter pile-up, Neptune ridge, and Neptune desert via orbital circularization and selective atmospheric mass loss.
-
Stellar Population Spectra Incorporating Detailed Binary Evolution using POSYDON
New spectral models of solar-metallicity stellar populations using POSYDON show binary-stripped stars dominating ionizing radiation after about 16 Myr.
-
A systematic survey for hypervelocity runaways from thermonuclear supernovae
A complete classification of Gaia-selected hypervelocity candidates yields 10 D6 and 3 LP 40-365 stars, with forward modeling favoring intermediate-heating evolutionary scenarios and low birth rates relative to the SN...
-
JWST's Little Red Dots as collapsed Supermassive Dark Stars
Supermassive dark stars powered by dark matter annihilation can collapse into quasi-stars whose envelopes expand and cool to match the observed properties of many JWST Little Red Dots while bypassing the restrictive c...
-
Escape of Water- and Metal-enriched Atmospheres from compact Hot mini-Neptunes with CHAIN
Water- and metal-rich atmospheres on compact hot mini-Neptunes lose mass more slowly than H/He cases at high enrichment levels due to enhanced cooling and higher mean molecular weight.
-
Temperature-resolved sensitivities of $^{56}{\rm Ni}$ production to helium-burning reactions in pair-instability supernovae
Temperature-resolved Monte Carlo analysis of PISNe finds peak sensitivity of 56Ni production to triple-alpha and 12C(alpha,gamma)16O rates at T~2.5e8 K with opposite signs, tied to pre-carbon C/O ratio.
-
Reconstruction of annual solar irradiance over the last three millennia
Physics-based annual TSI reconstruction over three millennia yields a maximum difference of 1.04 W/m² in 50-year running means.
-
Neutron star-companion interaction in core collapse supernovae. Population synthesis based on detailed binary evolution models
Population synthesis from binary evolution models predicts periodic neutron star-companion interactions in more than half of surviving hydrogen-poor core-collapse supernovae, with periods peaking at 20-50 days and las...
-
The Distribution of Blue Straggler Stars in the Color-Magnitude Diagrams of Old Open Clusters
Blue straggler stars in old open clusters predominantly appear near the terminal-age main sequence because mass transfer from asymptotic giant branch donors enriches their cores with helium.
-
Double Neutron Star Delay Times Across Cosmic Metallicities: The Role of Helium Star Progenitors
Simulations show double neutron star mergers peak 80-250 million years after star formation across metallicities, with 15% quick mergers and over 20% delayed over a billion years.
-
Hot blue progenitors of stellar-mass black holes
Black hole progenitors are predominantly hot and blue at pre-collapse, often Wolf-Rayet stars luminous in ultraviolet, with only a minority as red supergiants, and a direct-collapse rate of about 0.4 per century for a...
-
Evolution of Massive Main-sequence Stars in Rapid Population Synthesis. I. Framework and Implementation
A new main-sequence evolution framework implemented in COMPAS yields more massive helium cores, more compact stripped-star radii, and systematically higher black-hole masses than standard prescriptions.
-
Dark matter energy exchange in stars orbiting supermassive black holes
Orbit-averaged DM energy exchange for S4714 reaches stellar luminosity at σ_χp ~ 10^{-36} cm² (MeV-GeV) and σ_χe ~ 5×10^{-38} cm² (sub-MeV) for spiked DM profiles.
-
OCCAM X. Neutron Capture Abundances with Keck/HIRES & Magellan/MIKE
New high-resolution spectra yield abundances for 7 neutron-capture elements in open cluster stars, revealing flat Milky Way gradients for second-peak s- and r-process species and shallower slopes for first-peak s-process.
-
The location of the upper edge of the pair-instability supernovae black hole mass gap
The upper edge of the PISN black hole mass gap shifts by up to 30 solar masses due to nuclear reaction rate uncertainties, primarily the 12C(α,γ)16O rate, while remaining robust to resolution variations unlike the lower edge.
-
Elemental Abundances from Off-center Carbon Burning in Accreting CO White Dwarfs: Implications for SN 2021yfj-like events
Numerical simulations find that off-center carbon burning in accreting CO white dwarfs produces substantial Si and S abundances that depend on initial carbon content, supporting the double WD merger scenario for SN 20...
-
Twin Peaks: Resolving Features in the Binary Black Hole Mass Function with COSMIC-METISSE
New MESA stellar tracks with varied winds and convective mixing produce a primary black hole mass function with twin peaks near 8 and 13 solar masses in most variations, the higher peak dominated by mass-ratio-reversa...
-
Formation of Black Hole-White Dwarf X-ray Binaries in Globular Clusters
Simulations identify black hole-giant collisions as the main formation channel for ultracompact black hole-white dwarf binaries in globular clusters, with mock catalogs suggesting order-one LISA-detectable sources in ...
-
Kepler Image-Subtracted Light Curves and Variable Star Catalog of NGC 6819
Kepler image-subtracted photometry yields 81,498 light curves and a catalog of 87 periodic variable candidates in NGC 6819, including 26 newly discovered ones.
-
Kepler Image-Subtracted Light Curves and Variable Star Catalog of NGC 6819
Generates and publicly releases 81,498 detrended Kepler light curves plus a catalog of 87 periodic variables (26 new) in the 2.5 Gyr cluster NGC 6819 using Gaia DR3 for membership.
-
Observing a 542-day transiting giant with large TTVs: The 2025 transit of HIP 41378 f and new constraints on the outer system
New 2025 transit timing of HIP 41378 f shows a 7-hour early arrival consistent with TTVs; N-body modeling with TRADES refines ephemerides for planets d, e, and f.
-
Observing a 542-day transiting giant with large TTVs: The 2025 transit of HIP 41378 f and new constraints on the outer system
New 2025 transit timing for HIP 41378 f confirms large TTVs and is combined with prior data on planets d and e in an N-body model to update ephemerides and predict future transits.
-
The fate of Earth during the Sun's giant phases: New constraints from ab initio tidal modelling and AGB mass loss
Updated ab initio tidal models predict Earth survives the Sun's giant phases, but outcomes depend strongly on AGB mass-loss rates, with L2 Pup observations as proxy indicating survival while low rates lead to engulfment.
-
The effect of near-core mixing on rejuvenation and the asteroseismic properties of massive accretors
1D models show convective boundary mixing dominates the asteroseismic imprint of accretion in massive stars, robust to semiconvection changes but drastically altered without it, with thermal relaxation as key.
-
Oscillations of red giant stars with magnetic damping in the core. II. Mixed mode visibilities on the red-giant branch
Synthetic spectra show that observational biases cause dipole mode visibilities to be overestimated by up to 20 percent on the red-giant branch, while partial energy preservation under magnetic damping can produce bot...
-
Building three-dimensional giant stellar models for common envelope simulations
Depositing stellar luminosity in an inner shell and cooling low-density outer cells produces a stable pulsating 3D red supergiant model for common envelope simulations without relaxation.
-
Mass-Orbital Period Distribution of Massive White Dwarfs Formed Through Stable Mass Transfer
Simulations show that intermediate-mass progenitors with non-degenerate cores before helium burning produce a mass-orbital period relation for massive white dwarfs that accounts for long-period systems.
-
The Impact of the New $^{59}$Fe Decay Rates on $^{60}$Fe and $^{26}$Al Nucleosynthesis in Massive Stars
New 59Fe beta-decay rates in non-rotating MESA models suppress 60Fe yields by 0.28 dex, yielding a Salpeter-IMF-integrated 60Fe/26Al flux ratio of 0.18 that matches the observed 0.184 ± 0.042.
-
Mergers via failed common envelope as a route towards intermediate-mass stripped stars
Failed common envelope mergers yield 6-14 solar mass stripped stars consistent with long-lived core He-burning objects that appear single or in wide binaries from hierarchical triples.
-
Formation of extremely low-mass white dwarf binaries undergoing enhanced angular momentum loss
Enhanced AML via L2-point mass loss in the RLOF channel alters ELM WD internal structure and mass-radius relation, reproducing observed shorter orbital periods.
-
Spectroscopic Disentangling Revealed the Tertiary Component in the Multiple System EM Boo
Tertiary in EM Boo is A-F type with Teff=7000K; system distance ~300 pc indicates Gaia DR3 underestimates true distance due to multiplicity.
-
Double Neutron Star Delay Times Across Cosmic Metallicities: The Role of Helium Star Progenitors
Population synthesis of helium star-NS systems yields DNS delay time distributions that peak between 80-250 Myr across metallicities, with 15% merging within 80 Myr and over 20% after 1 Gyr.
-
A Theoretical Study of the Structure and Elemental Abundances of HD 20794
MESA grid models find HD 20794 is a 0.80 solar-mass star aged ~9 Gyr whose observed abundances match core-collapse supernova enrichment and are preserved over Gyr timescales.
-
Hot blue progenitors of stellar-mass black holes
Stellar evolution and atmosphere models predict black hole progenitors are predominantly hot and blue with a direct-collapse rate of ~0.4 per century in a 1 Msun/yr star-forming galaxy.
-
Deep Adaptive Optics Imaging Rules Out a Helium Star Companion to PSR J1928+1815
Deep Keck/NIRC2 adaptive optics imaging rules out a helium star companion to PSR J1928+1815, supporting a massive white dwarf with possible wind-driven eclipses.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.