Josiah Schwab

Regular, automated testing is a foundational principle of modern software development. Numerous widely used continuous integration systems exist, but they are often not suitable for the unique needs of scientific simulation software. Here we describe the testing infrastructure developed for and used by the Modules for Experiments in Stellar Astrophysics (MESA) project. This system allows the computationally demanding MESA test suite to be regularly run on a heterogeneous set of computers and aggregates and displays the testing results in a form that allows for the rapid identification and diagnosis of regressions. Regularly collecting comprehensive testing data also enables longitudinal studies of the performance of the software and the properties of the models it generates.

Regular, automated testing is a foundational principle of modern software development. Numerous widely-used continuous integration systems exist, but they are often not suitable for the unique needs of scientific simulation software. Here we describe the testing infrastructure developed for and used by the Modules for Experiments in Stellar Astrophysics (MESA) project. This system allows the computationally-demanding MESA test suite to be regularly run on a heterogeneous set of computers and aggregates and displays the testing results in a form that allows for the rapid identification and diagnosis of regressions. Regularly collecting comprehensive testing data also enables longitudinal studies of the performance of the software and the properties of the models it generates.

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, 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 in MESA with 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 strengthen MESA’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 in MESA, 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 to MESA’s software infrastructure that enhance source code development and community engagement.

Helium star–carbon-oxygen white dwarf (CO WD) binaries are potential single-degenerate progenitor systems of thermonuclear supernovae. Revisiting a set of binary evolution calculations using the stellar evolution code MESA, we refine our previous predictions about which systems can lead to a thermonuclear supernova and then characterize the properties of the helium star donor at the time of explosion. We convert these model properties to near-UV/optical magnitudes assuming a blackbody spectrum and support this approach using a matched stellar atmosphere model. These models will be valuable to compare with pre-explosion imaging for future supernovae, though we emphasize the observational difficulty of detecting extremely blue companions. The pre-explosion source detected in association with SN 2012Z has been interpreted as a helium star binary containing an initially ultra-massive WD in a multiday orbit. However, extending our binary models to initial CO WD masses of up to 1.2 M _⊙, we find that these systems undergo off-center carbon ignitions and thus are not expected to produce thermonuclear supernovae. This tension suggests that, if SN 2012Z is associated with a helium star–WD binary, then the pre-explosion optical light from the system must be significantly modified by the binary environment and/or the WD does not have a carbon-rich interior composition.

Recently, a class of Roche-lobe-filling binary systems consisting of hot subdwarf stars and white dwarfs (WDs) with sub-hour periods has been discovered. At present, the hot subdwarf is in a shell He-burning phase and is transferring some of its remaining thin H envelope to its WD companion. As the evolution of the hot subdwarf continues, it is expected to detach, leaving behind a low-mass C/O-core WD secondary with a thick He layer. Then, on a timescale of ∼10 Myr, gravitational wave radiation will again bring the systems into contact. If the mass transfer is unstable and results in a merger and a catastrophic thermonuclear explosion is not triggered, it creates a remnant with a C/O-dominated envelope, but one still rich enough in He to support an R Corona Borealis-like shell-burning phase. We present evolutionary calculations of this phase and discuss its potential impact on the cooling of the remnant WD.

We present a set of ultramassive white dwarf models, focused on masses above 1.3 M _⊙. Given the uncertainties about the formation and compositions of such objects, we construct parameterized model sequences, guided by evolutionary calculations including both single star and double white dwarf merger formation channels. We demonstrate that the cooling of objects with central densities in excess of 10⁹ g cm⁻³ is dominated by neutrino cooling via the Urca process in the first ≈100 Myr after formation. Our models indicate that the recently discovered ultramassive white dwarf ZTF J190132.9+145808.7 is likely to have experienced this Urca-dominated cooling regime. We also show that the high densities imply that diffusion is unlikely to significantly alter the core compositions of these objects before they crystallize.

In this work we derive the minimum allowed orbital periods of H-rich bodies ranging in mass from Saturn’s mass to 1 M _⊙, emphasizing gas giants and brown dwarfs (BDs) over the range 0.0003–0.074 M _⊙. Analytic fitting formulae for as a function of the mass of the body and as a function of the mean density are presented. We assume that the density of the host star is sufficiently high so as not to limit the minimum period. In many instances this implies that the host star is a white dwarf. This work is aimed, in part, toward distinguishing BDs from planets that are found transiting the host white dwarf without recourse to near-infrared or radial velocity measurements. In particular, orbital periods of ≲100 minutes are very likely to be BDs. The overall minimum period over this entire mass range is ≃37 minutes.

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 novel thermodynamic extrapolation scheme 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 the MESA stellar evolution software instrument by computing white dwarf cooling curves.

We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M _⊙ carbon–oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of ²²Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6−7 significant figures of precision. Usually white dwarf models without ²²Ne are fit to the observed periods and other properties. The rms residuals to the ≃150−400 s low-order g-mode periods are typically in the range of σ _rms ≲ 0.3 s, for a fit precision of σ _rms/P ≲ 0.3%. We find average relative period shifts of ΔP/P ≃ ±0.5% for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of ΔP/P depending on the specific g-mode, abundance of ²²Ne, effective temperature, and the mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when ²²Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of ²²Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including ²²Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model-fitting processes.

We construct evolutionary models of the remnant of the merger of two carbon-oxygen (CO) core white dwarfs (WDs). With total masses in the range 1–2 M _⊙, these remnants may either leave behind a single massive WD or undergo a merger-induced collapse to a neutron star (NS). On the way to their final fate, these objects generally experience a ∼10 kyr luminous giant phase, which may be extended if sufficient helium remains to set up a stable shell-burning configuration. The uncertain, but likely significant, mass-loss rate during this phase influences the final remnant mass and fate (WD or NS). We find that the initial CO core composition of the WD is converted to oxygen-neon (ONe) in remnants with final masses ≳1.05 M _⊙. This implies that the CO core/ONe core transition in single WDs formed via mergers occurs at a similar mass as in WDs descended from single stars and thus that WD–WD mergers do not naturally provide a route to producing ultramassive CO-core WDs. As the remnant contracts toward a compact configuration, it experiences a “bottleneck” that sets the characteristic total angular momentum that can be retained. This limit predicts that single WDs formed from WD–WD mergers have rotational periods of ≈10–20 minutes on the WD cooling track. Similarly, it predicts remnants that collapse can form NSs with rotational periods ∼10 ms.