Correcting NIRISS charge migration reveals KELT-7 b as a metal-enriched (~92x solar) ultra-hot Jupiter with H2O, CO2 and TiO but no H- or clouds, orbiting a young metal-rich star.
hub Canonical reference
Modules for Experiments in Stellar Astrophysics (MESA): Planets, Oscillations, Rotation, and Massive Stars
Canonical reference. 80% of citing Pith papers cite this work as background.
abstract
We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion that is more easily implemented when many nuclei are present at non-negligible abundances. This impacts the way in which MESA Star calculates semi-convective and thermohaline mixing. We exhibit the evolution of 3-8 Msun stars through the end of core He burning, the onset of He thermal pulses, and arrival on the white dwarf cooling sequence. We implement diffusion of angular momentum and chemical abundances that enable calculations of rotating-star models, which we compare thoroughly with earlier work. We introduce a new treatment of radiation-dominated envelopes that allows the uninterrupted evolution of massive stars to core collapse. This enables the generation of new sets of supernovae, long gamma-ray burst, and pair-instability progenitor models. We substantially modify the way in which MESA Star solves the fully coupled stellar structure and composition equations, and we show how this has improved MESA's performance scaling on multi-core processors. Updates to the modules for equation of state, opacity, nuclear reaction rates, and atmospheric boundary conditions are also provided. We describe the MESA Software Development Kit (SDK) that packages all the required components needed to form a unified and maintained build environment for MESA. [Abridged]
hub tools
citation-role summary
citation-polarity summary
claims ledger
- abstract We substantially update the capabilities of the open source software package Modules for Experiments in Stellar Astrophysics (MESA), and its one-dimensional stellar evolution module, MESA Star. Improvements in MESA Star's ability to model the evolution of giant planets now extends its applicability down to masses as low as one-tenth that of Jupiter. The dramatic improvement in asteroseismology enabled by the space-based Kepler and CoRoT missions motivates our full coupling of the ADIPLS adiabatic pulsation code with MESA Star. This also motivates a numerical recasting of the Ledoux criterion t
- dataset (2005), [11] Soubiran & Girard (2005), [12] Valenti & Fischer (2005), [13] Reddy et al. (2006), [14] Takeda et al. (2007), [15] Ram'ırez et al. (2007), [16] Sousa et al. (2008), [17] Holmberg et al. (2009), [18] Mena et al. (2010), [19] Ghezzi et al. (2010), [20] Gonzalez et al. (2010), [21] Prugniel et al. (2011), [22] Adibekyan et al. (2012), [23] Ram 'ırez et al. (2012), [24] Maldonado et al. (2012), [25] Da Silva et al. (2012), [26] Ram 'ırez et al. (2013), [27] Tsantaki et al. (2013), [28]
- background Jaeckel,Probing Minicharged Particles with Tests of Coulomb's Law,Phys. Rev. Lett.103(2009) 080402 [0904.1547]. [81] J. Jaeckel and A. Ringwald,The Low-Energy Frontier of Particle Physics,Ann. Rev. Nucl. Part. Sci.60 (2010) 405 [1002.0329]. [82] M. D. Diamond and P. Schuster,Searching for Light Dark Matter with the SLAC Millicharge Experiment, Phys. Rev. Lett.111(2013) 221803 [1307.6861]. [83] H. Vogel and J. Redondo,Dark Radiation constraints on minicharged particles in models with a hidden pho
- background companion star,E heat, which is given by Eheat =p·E kin,ej · ˜Ωeff,(1) withp=0.08 an efficiency factor (Hirai et al. 2018). Here, ˜Ωeff is the fractional solid angle subtended by the companion as seen from the exploding progenitor, which is a function of the com- panion's radiusR2 and its separationa: ˜Ωeff= Ω(R,a) 4π = 1 2 1− s 1− R2 a 2 .(2) The intersected energyE heat drives the stellar expansion un- til the luminosity matches the local Eddington luminosity in the lay
- background The brighter the star and/or the lower its effective temperature, the higher the mass-loss rate. We vary the linear scaling parameterη, choosing values of 0.2, 0.4, 0.8, and 1.0. Modern simulations sometimes adopt values ofη <1.0 to account for wind inhomo- geneities ("clumping") and other evidence suggesting the standard "Dutch" prescription withη= 1.0 is too ag- gressive [47, 48]. Note that Renzo et al. [19] showed that the choice of scaling parameterηis generally much more influential than th
- background (2013), arXiv:1301.0319 [astro-ph.SR]. [48] B. Paxton, P. Marchant, J. Schwab, E. B. Bauer, L. Bild- sten, M. Cantiello, L. Dessart, R. Farmer, H. Hu, N. Langer, R. H. D. Townsend, D. M. Townsley, and F. X. Timmes, Modules for Experiments in Stellar Astro- physics (MESA): Binaries, Pulsations, and Explosions, ApJS220, 15 (2015), arXiv:1506.03146 [astro-ph.SR]. [49] B. Paxton, J. Schwab, E. B. Bauer, L. Bildsten, S. Blinnikov, P. Duffell, R. Farmer, J. A. Goldberg, P. Marchant, E. Sorokina, A. Th
authors
co-cited works
representative citing papers
TOI-2195 A b is an inflated hot Neptune that likely originated as a Jovian planet losing ~90% mass through Roche lobe overflow during EKL-driven high-eccentricity migration triggered by a wide binary companion.
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.
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.
Analytical relations for convective core mass, variable-mu structure, and main-sequence end in intermediate-mass stars, with MESA validation.
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 migration.
Gravity-wave polarization theory is extended to non-axisymmetric stellar magnetic fields, revealing avoided crossings and a new mixing threshold distinct from the magnetic-suppression threshold.
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 in the less-truncated run.
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.
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-mass star in a Hubble time.
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.
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 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 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.
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.
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.
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.
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 35-200 keV and charges 10^{-10} to 10^{-9}.
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 during the quenching transition of massive galaxies.
SN 2020bij and four other Type IIP SNe with slow-rising light curves and high velocities are modeled with weak to no CSM interaction, suggesting a new subclass linked to confined CSM.
Orbit-averaged elastic DM scattering in S4714 reaches stellar luminosity at σ_χp ∼ 10^{-36} cm² (MeV–GeV) and σ_χe ∼ 5×10^{-38} cm² (sub-MeV) for a spiked profile.
Theta Eridani's historical brightness is explained as a millenary common-envelope transient powered by orbital energy extraction in a previously more eccentric binary.
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 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.
citing papers explorer
-
Mitigating Charge Migration in JWST NIRISS Reveals That KELT-7 b is a Metal-enriched Ultra-hot Jupiter Orbiting a Young Metal-rich Star
Correcting NIRISS charge migration reveals KELT-7 b as a metal-enriched (~92x solar) ultra-hot Jupiter with H2O, CO2 and TiO but no H- or clouds, orbiting a young metal-rich star.
-
Discovery of an Inflated Hot Neptune and Its Formation from Jovian Mass Loss
TOI-2195 A b is an inflated hot Neptune that likely originated as a Jovian planet losing ~90% mass through Roche lobe overflow during EKL-driven high-eccentricity migration triggered by a wide binary companion.
-
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.
-
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 migration.
-
Extending asteroseismic magnetometry across the diverse landscape of magnetic structures
Gravity-wave polarization theory is extended to non-axisymmetric stellar magnetic fields, revealing avoided crossings and a new mixing threshold distinct from the magnetic-suppression 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 in the less-truncated run.
-
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-mass star in a Hubble time.
-
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.
-
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.
-
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.
-
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 35-200 keV and charges 10^{-10} to 10^{-9}.
-
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 during the quenching transition of massive galaxies.
-
SN 2020bij and a Possible Slow-Rise High-Velocity Subclass of Type IIP Supernovae
SN 2020bij and four other Type IIP SNe with slow-rising light curves and high velocities are modeled with weak to no CSM interaction, suggesting a new subclass linked to confined CSM.
-
Dark matter energy exchange in stars orbiting supermassive black holes
Orbit-averaged elastic DM scattering in S4714 reaches stellar luminosity at σ_χp ∼ 10^{-36} cm² (MeV–GeV) and σ_χe ∼ 5×10^{-38} cm² (sub-MeV) for a spiked profile.
-
The forgotten bright star: Theta Eridani as a millenary stellar transient observed by Hipparchus, Ptolemy and al-Sufi
Theta Eridani's historical brightness is explained as a millenary common-envelope transient powered by orbital energy extraction in a previously more eccentric binary.
-
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 Ia rate.
-
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 conditions of nuclear-powered supermassive star formation.
-
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.
-
Subluminous Type IIP SN 2024abfl as a Result of a Significantly Low-energy Fe-core Collapse
SN 2024abfl is a subluminous Type IIP event from a low-energy Fe-core collapse of a compact 9–10 M⊙ progenitor with only ~0.003 M⊙ of nickel.
-
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 lasting 0.5-10 years.
-
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.
-
Revisiting Turner Window Axions: The Untapped Potential of NaI Dark Matter Detectors
Revised supernova opacity modeling reopens the Turner window for hadronically coupled axions, enabling resonant absorption searches in NaI detectors for couplings |g_app| between 10^{-6.5} and 10^{-2}.
-
Impacts of Multidimensional Progenitor Perturbations on Core-Collapse Supernova Explosions
Multi-D progenitor perturbations leave shock revival and explosion energetics essentially unchanged in 2D CHIMERA runs; any effects sit below the level of numerical stochasticity.
-
Deep Adaptive Optics Imaging Rules Out a Helium Star Companion to PSR J1928+1815
Deep NIR adaptive-optics imaging rules out a helium-star companion to the eclipsing millisecond pulsar PSR J1928+1815, favoring a massive white dwarf.
-
On the Importance of the Convective Urca Process in 3D Simulations of a Simmering White Dwarf
3D simulations show the convective Urca process substantially reduces the convection zone size in a simmering white dwarf, though convection extends past the Urca shell.
-
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.
-
Effects of Rotation on 3D Core-Collapse Supernova Models for Low-Mass Progenitors
For a low-mass CCSN progenitor, rotation alone weakly and non-monotonically affects explosion energy and observables; only the fastest spin yields T/|W| spiral modes and spin-kick alignment, with core spin amplified by ~4000.
-
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 2021yfj-like events.
-
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-reversal systems, with rates varying by a factor of six.
-
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 the Milky Way.
-
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 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 both present and absent mixed-mode signatures.
-
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.