Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
read the original abstract
We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.
This paper has not been read by Pith yet.
Forward citations
Cited by 27 Pith papers
-
Cepheids with giant companions III. Evolutionary modeling of nine binary double Cepheids from the Milky Way and Magellanic Clouds
Extended q-PED modeling provides first mass estimates for SMC Cepheids and new period-mass-radius and mass-luminosity relations while indicating binary interactions in several systems.
-
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.
-
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.
-
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...
-
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.
-
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...
-
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.
-
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.
-
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 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.
-
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.
-
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.
-
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...
-
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.
-
Detection horizon for the neutrino burst from the stellar helium flash
The helium flash produces a neutrino burst with a 1.7 MeV line detectable up to almost 3 pc in future facilities, but asteroseismology remains the practical probe for now.
-
Comparative Study of Two Luminous Red Novae I. Progenitor Modeling and Dust Formation
Binary evolution modeling constrains donor masses of 14-23 solar masses for two luminous red novae and shows dust masses are 1-5 orders of magnitude below total ejected envelope masses.
-
Shaping the horizontal branch: The role of envelope mass in the evolution of stripped core-helium-burning stars
MESA models show residual hydrogen envelope mass sets effective temperature on the horizontal branch, with maximum values of 0.05-0.30 solar masses to avoid later thermally pulsing AGB evolution, plus explanations for...
-
Millicharged Particle Constraints from Asymptotic Giant Branch Stars
New upper bounds on millicharged particles (masses 10-100 keV, charges down to 5e-13) from the AGB-to-HB star ratio in globular clusters, improving prior limits by up to two orders of magnitude.
-
Modeling of the magnetic stellar wind braking of the ssrAp 33 Lib (HD137949)
Magnetic wind braking can slow ssrAp stars like 33 Lib to rotation periods of 80 years or more according to MESA evolution models that include magnetic field and wind changes.
-
Diagnosing the Properties and Evolutionary Fates of Black Hole and Wolf-Rayet X-ray Binaries as Potential Gravitational Wave Sources for the LIGO-Virgo-KAGRA Network
MESA binary evolution simulations with revised Bondi-Hoyle accretion efficiency and observational constraints yield lower BH mass upper limits for IC 10 X-1 and NGC 300 X-1 and predict Hubble-time BBH mergers for all ...
-
Study of photometric and spectral variability of the roAp star HD~210684
Observational analysis of HD210684 yields rotational period 5.02188 d, magnetic obliquity 77°, and main-sequence age ~1.45 Gyr from photometry, LSD magnetic measurements, and MESA/GYRE modeling.
-
Can current models predict the local black hole merger rate?
Theoretical predictions for local BBH merger rates exceed observations by a factor >10 under conservative SFRD and metallicity assumptions, indicating need for revisions in stellar evolution.
-
Red novae, their progenitors, and remnants
A synthesis of observational data on red novae as stellar merger events, including outburst properties, progenitor diversity, and long-term remnants.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.