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.
hub
Modules for Experiments in Stellar Astrophysics (MESA)
38 Pith papers cite this work. Polarity classification is still indexing.
abstract
Stellar physics and evolution calculations enable a broad range of research in astrophysics. Modules for Experiments in Stellar Astrophysics (MESA) is a suite of open source libraries for a wide range of applications in computational stellar astrophysics. A newly designed 1-D stellar evolution module, MESA star, combines many of the numerical and physics modules for simulations of a wide range of stellar evolution scenarios ranging from very-low mass to massive stars, including advanced evolutionary phases. MESA star solves the fully coupled structure and composition equations simultaneously. It uses adaptive mesh refinement and sophisticated timestep controls, and supports shared memory parallelism based on OpenMP. Independently usable modules provide equation of state, opacity, nuclear reaction rates, and atmosphere boundary conditions. Each module is constructed as a separate Fortran 95 library with its own public interface. Examples include comparisons to other codes and show evolutionary tracks of very low mass stars, brown dwarfs, and gas giant planets; the complete evolution of a 1 Msun star from the pre-main sequence to a cooling white dwarf; the Solar sound speed profile; the evolution of intermediate mass stars through the thermal pulses on the He-shell burning AGB phase; the interior structure of slowly pulsating B Stars and Beta Cepheids; evolutionary tracks of massive stars from the pre-main sequence to the onset of core collapse; stars undergoing Roche lobe overflow; and accretion onto a neutron star. Instructions for downloading and installing MESA can be found on the project web site (http://mesa.sourceforge.net/).
hub tools
citation-role summary
citation-polarity summary
authors
representative citing papers
2D disc simulations with vertical volatile transport produce stable CO snow surface equilibria and eliminate limit-cycle behavior seen in 1D models.
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.
Extends magnetogravity polarization formalism to arbitrary magnetic field geometries, revealing avoided crossings and mode conversion below a local field threshold.
FIRE-2 simulations show per-galaxy tidal disruption rates peak near z=2.5 at 4e-4 per year, correlate with SFR and central density, and remain high in satellite galaxies at early times.
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.
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.
Gas giants form sub-snowline in binaries via dust traps at the tidal truncation radius, with observed planet semi-major axes following a_planet = 0.569 r_t (R²=0.94).
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.
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}.
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.
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.
GD 1212 showed the largest recorded outburst in a DAV white dwarf, attributed to parametric instability, with refined 17-hour rotation from mode splittings and a failed outburst event.
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 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.
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.
Hot Jupiter occurrence in the Galactic halo is low at ~0.13% with no significant difference between in-situ and accreted populations, well below disk rates.
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}.
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.
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 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.
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.
citing papers explorer
-
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.
-
CO snow lines are stabilised by the vertical transport of volatiles
2D disc simulations with vertical volatile transport produce stable CO snow surface equilibria and eliminate limit-cycle behavior seen in 1D models.
-
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.
-
TDEs on FIRE: Illuminating the Cosmic Evolution of Tidal Disruption Rates
FIRE-2 simulations show per-galaxy tidal disruption rates peak near z=2.5 at 4e-4 per year, correlate with SFR and central density, and remain high in satellite galaxies at early times.
-
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.
-
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.
-
Sub-Snowline Formation of Gas-Giant Planets in Binary Systems
Gas giants form sub-snowline in binaries via dust traps at the tidal truncation radius, with observed planet semi-major axes following a_planet = 0.569 r_t (R²=0.94).
-
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.
-
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}.
-
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.
-
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.
-
Mode Instability and a Massive, Isolated Outburst in the Pulsating White Dwarf GD 1212
GD 1212 showed the largest recorded outburst in a DAV white dwarf, attributed to parametric instability, with refined 17-hour rotation from mode splittings and a failed outburst event.
-
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.
-
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.
-
Exoplanets in ancient stellar populations: occurrence constraints and hot-Jupiter candidates in the Galactic halo
Hot Jupiter occurrence in the Galactic halo is low at ~0.13% with no significant difference between in-situ and accreted populations, well below disk rates.
-
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}.
-
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.
-
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.
-
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.
-
On the stellar parameter dependence of the combined Fe I and Fe II chromospheric emission-line in the wings of the Ca II K line
Observational analysis of 21 giant stars shows Fe I+Fe II blend emission proportional to Ca II K with matching Teff power-law exponent and a slope break in the flux ratio at log g ≈ 2.5.
-
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 blue hook stars and puffed-up pre-HB configurations.
-
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 three systems.
-
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.
-
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.