skip to main content

Title: Modules for Experiments in Stellar Astrophysics (MESA): Time-dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The newauto_diffmodule implements automatic differentiation inMESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite-difference approximations. We significantly enhance the treatment of the growth and decay of convection inMESAwith a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron-degenerate ignition events. We strengthenMESA’s implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars inMESA, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator-split nuclear burning mode. We close by discussing major updates toMESA’s software infrastructure that enhance source code development and community more » engagement.

« less
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; « less
Award ID(s):
1663684 1663696
Publication Date:
Journal Name:
The Astrophysical Journal Supplement Series
Page Range or eLocation-ID:
Article No. 15
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Gravitational-wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernovae (PISNe). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections of merging BHs. We useMESAto evolve single, nonrotating, massive helium cores with a metallicity ofZ= 10−5, until they either collapse to form a BH or explode as a PISN, without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the ±3σuncertainty in our high-resolution tabulated12C(α,γ)16O reaction rate probability distribution function. We extensively test temporal and spatial resolutions for resolving the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower-resolution models reach convergence. We establish a new lower edge of the upper mass gap asMlower6014+32Mfrom the ±3σuncertainty inmore »the12C(α,γ)16O rate. We explore the effect of a larger 3αrate on the lower edge of the upper mass gap, findingMlower6918+34M. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog.

    « less
  2. Abstract

    Feedback from massive stars plays an important role in the formation of star clusters. Whether a very massive star is born early or late in the cluster formation timeline has profound implications for the star cluster formation and assembly processes. We carry out a controlled experiment to characterize the effects of early-forming massive stars on star cluster formation. We use the star formation software suiteTorch, combining self-gravitating magnetohydrodynamics, ray-tracing radiative transfer,N-body dynamics, and stellar feedback, to model four initially identical 104Mgiant molecular clouds with a Gaussian density profile peaking at 521.5 cm−3. Using theTorchsoftware suite through theAMUSEframework, we modify three of the models, to ensure that the first star that forms is very massive (50, 70, and 100M). Early-forming massive stars disrupt the natal gas structure, resulting in fast evacuation of the gas from the star-forming region. The star formation rate is suppressed, reducing the total mass of the stars formed. Our fiducial control model, without an early massive star, has a larger star formation rate and total efficiency by up to a factor of 3, and a higher average star formation efficiency per freefall time by up to a factor of 7. Early-forming massive stars promote the buildupmore »of spatially separate and gravitationally unbound subclusters, while the control model forms a single massive cluster.

    « less

    Planetary engulfment events can occur while host stars are on the main sequence. The addition of rocky planetary material during engulfment will lead to refractory abundance enhancements in the host star photosphere, but the level of enrichment and its duration will depend on mixing processes that occur within the stellar interior, such as convection, diffusion, and thermohaline mixing. We examine engulfment signatures by modelling the evolution of photospheric lithium abundances. Because lithium can be burned before or after the engulfment event, it produces unique signatures that vary with time and host star type. Using mesa stellar models, we quantify the strength and duration of these signatures following the engulfment of a 1, 10, or 100 M⊕ planetary companion with bulk Earth composition, for solar-metallicity host stars with masses ranging from 0.5 to 1.4 M⊙. We find that lithium is quickly depleted via burning in low-mass host stars ($\lesssim 0.7 \, {\rm M}_\odot$) on a time-scale of a few hundred Myrs, but significant lithium enrichment signatures can last for Gyrs in G-type stars ($\sim \! 0.9 \, {\rm M}_{\odot }$). For more massive stars (1.3−1.4 M⊙), engulfment can enhance internal mixing and diffusion processes, potentially decreasing the surface lithium abundance. Our predicted signatures frommore »exoplanet engulfment are consistent with observed lithium-rich solar-type stars and abundance enhancements in chemically inhomogeneous binary stars.

    « less
  4. Abstract

    Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. We present Skye, a new equation of state for fully ionized matter. Skye includes the effects of positrons, relativity, electron degeneracy, Coulomb interactions, nonlinear mixing effects, and quantum corrections. Skye determines the point of Coulomb crystallization in a self-consistent manner, accounting for mixing and composition effects automatically. A defining feature of this equation of state is that it uses analytic free energy terms and provides thermodynamic quantities using automatic differentiation machinery. Because of this, Skye is easily extended to include new effects by simply writing new terms in the free energy. We also introduce a novelthermodynamic extrapolationscheme for extending analytic fits to the free energy beyond the range of the fitting data while preserving desirable properties like positive entropy and sound speed. We demonstrate Skye in action in theMESAstellar evolution software instrument by computing white dwarf cooling curves.

  5. Abstract

    We present the discovery of a new double-detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with aPorb= 76.34179(2) minutes orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits, we find a low-mass sdB star,MsdB= 0.383 ± 0.028Mwith a massive white dwarf companion,MWD= 0.725 ± 0.026M. From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of ≈25 Myr whereas ourMESAmodel predicts an sdB age of ≈170 Myr. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion ≈25 Myr ago. Using theMESAstellar evolutionary code we find that the sdB star will start mass transfer in ≈6 Myr and in ≈60 Myr the white dwarf will reach a total mass of 0.92Mwith a thick helium layer of 0.17M. This will lead to a detonation that will likely destroymore »the white dwarf in a peculiar thermonuclear supernova. PTF1 J2238+7430 is only the second confirmed candidate for a double-detonation thermonuclear supernova. Using both systems we estimate that at least ≈1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions.

    « less