Evan B. Bauer

Experimental validation of equation-of-state (EOS) models for high-energy-density applications typically involves measuring the shock Hugoniot, which helps constrain the high-pressure, high-density EOS. However, EOS measurements at below ambient density in the expanded regime are less common. Recent experimental advances allow the isentrope to be inferred from measurements of the density profile of an adiabatically expanding proton-heated plasma, offering the potential to assess various theoretical EOS models over a wide range of below-ambient density conditions. On the theoretical side, determining the expanded-regime EOS poses a considerable challenge for ab initio methods like density functional theory--based molecular dynamics (MD) due to the big simulation cells and large number of orbitals required to simulate these low-density, high-temperature conditions. We thus take a more computationally friendly approach by leveraging the power of MD driven by machine-learned interatomic potentials (ML-IAPs) trained on ab initio simulations. We directly simulate the adiabatic expansion of heated diamond at conditions relevant to experiments using large cells with up to 1.3 million atoms. These simulations, using several different variants of carbon ML-IAPs, are used to extract release isentropes and to study possible microscopic chemistry and phase transitions at play. These insights gained from our MD simulations can serve to guide the design and interpret observations of experiments.

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense-matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. Here spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D_{2}O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense-matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

The hydrogen envelope is the outermost layer of a DA white dwarf; it makes up the entirety of the stellar photosphere, and yet its typical extent is difficult to model theoretically and remains poorly observationally constrained. As a result, hydrogen envelope mass is a substantial source of systematic uncertainty in the physical properties of white dwarfs, including overall masses and cooling ages. In this work, we fit a Gaussian mixture model to gravitational redshifts from high-resolution spectroscopy, paired with radius measurements from Gaia BP/RP spectra, to measure the mass–radius relation for a sample of 468 white dwarfs. Our results are in excellent agreement with the predicted mass–radius relations of state-of-the-art evolutionary models, including those from the MESA Isochrones and Stellar Tracks (MIST) library. We find that mass–radius relations such as those from MIST that assume a thick and mass-dependent hydrogen envelope are preferred by the observed probability density function over models that assume a hydrogen envelope of constant mass. Proper treatment of the evolution of white dwarf progenitors is thus important for accurately modeling the mass–radius relation. Our results indicate that gravitational redshift measurements of large samples of white dwarfs in wide binaries are promising probes of the hydrogen envelope masses of DA white dwarfs.

We update and expand the MESA Isochrones and Stellar Tracks (MIST) database to include variations in the α -capture elements, specifically [ α /Fe] = −0.2, 0, +0.2, +0.4, and +0.6 for the full range of −3 ≤ [Fe/H] ≤ +0.5 already covered by MIST. Variations in [ α /Fe] are included in a self-consistent manner from the stellar interior models to the synthetic spectra used to translate these models in the observational plane. We describe a number of updates to the physics utilized in these models as well as new information provided by the models. We validate the models with comparisons to other stellar evolution models including the previous generation of MIST and other models from the literature. MIST data products including stellar evolution tracks, isochrones, and bolometric correction tables can be obtained from the MIST project website, https://mist.science . All necessary files to reproduce MIST models are available from Zenodo.

We present a substantial update to the MESA Isochrones and Stellar Tracks (MIST) library, extending the MIST model grids and isochrones down the white dwarf (WD) cooling sequence with realistic physics for WD cooling timescales. This work provides a large grid of MESA models for carbon–oxygen core WDs with hydrogen atmospheres (spectral type DA/DC), descended from full prior stellar evolution calculations. The model tracks, isochrones, and WD cooling timescale contours are available on the MIST project website and at doi: 10.5281/zenodo.15242047 . Our WD models provide a very large, publicly available grid with detailed physics for WD cooling timescales: realistic interior and envelope compositions, with element diffusion and heavy-element sedimentation, nuclear burning at the base of the WD hydrogen envelope, core crystallization, and C/O phase separation. As a large grid of open-source stellar evolution models, these WD models provide both out-of-the-box model tracks for comparison with observations and a framework for building further WD models to investigate variations in WD physics.

A significant number of hypervelocity stars with velocities between 1500 − 2500 km s ⁻¹ have recently been observed. The only plausible explanation so far is that they were produced through thermonuclear supernovae in white dwarf binaries. Since these stars are thought to be surviving donors of Type Ia supernovae, a surprising finding was that these stars are inflated, with radii an order of magnitude higher than expected for Roche-lobe-filling donors. Recent attempts at explaining them have combined 3D hydrodynamical supernova explosion simulations with 1D stellar modelling to explain the impact of supernova shocks on runaway white dwarfs. However, only the hottest and most compact of those runaway stars can so far marginally be reproduced by detailed models of runaways from supernova explosions. In this and a companion paper, we introduce a new A REPO simulation of two massive CO white dwarfs that explode via a violent merger. During the merger, the primary white dwarf ignites when the secondary is on its last orbit and plunging towards the primary. In the corresponding aftermath, the core of the secondary white dwarf of 0.16 M _⊙ remains bound, moving at a velocity of ∼2800 km s ⁻¹ . We mapped this object into MESA and show that this runaway star can explain the observations of two hypervelocity stars that were dubbed D6-1 and D6-3 based on their original discovery motivated by the D6 scenario, though the violent merger scenario presented here is somewhat distinct from the D6 scenario.

Binary systems of two carbon-oxygen white dwarfs are one of the most promising candidates for the progenitor systems of Type Ia supernovae. Violent mergers, where the primary white dwarf ignites when the secondary white dwarf smashes into it while being disrupted on its last orbit, were the first double degenerate merger scenario proposed that ignites dynamically. However, violent mergers likely contribute only a few percent to the total Type Ia supernova rate and do not yield normal Type Ia supernova light curves. Here we revisit the scenario, simulating a violent merger with better methods and, in particular, a more accurate treatment of the detonation. We find good agreement with previous simulations but with one critical difference: The secondary white dwarf being disrupted and accelerated towards the primary white dwarf and impacted by its explosion does not fully burn, and its core survives as a bound object. The explosion leaves behind a 0.16 M _⊙ star travelling 2800 km/s, making it an excellent (and so far the only) candidate to explain the origin of the fastest observed hypervelocity stars. We also show that before the explosion, 5 × 10 ⁻³ M _⊙ of material predominantly consisting of helium, carbon, and oxygen had already been ejected at velocities above 1000 km/s. Finally, we argue that if a violent merger made the hypervelocity stars D6-1 and D6-3 and violent mergers require the most massive primary white dwarfs in binaries of two carbon-oxygen white dwarfs, there has to be a much larger population of white dwarf mergers with slightly lower mass primary white dwarfs. Because this population likely represents ≫10% of the Type Ia supernovae rate, it can essentially only give rise to normal Type Ia supernovae.

Type Ia supernovae (SNe Ia), critical for studying cosmic expansion, arise from thermonuclear explosions of white dwarfs, but their precise progenitor pathways remain unclear. Growing evidence supports the “double-degenerate scenario,” where two white dwarfs interact. The absence of nondegenerate companions capable of explaining the observed SN Ia rate, along with observations of hypervelocity white dwarfs, interpreted as surviving companions of such systems, provide compelling evidence for this scenario. Upcoming millihertz gravitational-wave observatories like the Laser Interferometer Space Antenna (LISA) are expected to detect thousands of double-degenerate systems, though the most compact known candidate SN Ia progenitors produce marginally detectable signals. Here, we report observations of ATLAS J1138-5139, a binary white dwarf system with an orbital period of just 28 minutes. Our analysis reveals a 1 M _☉ carbon–oxygen white dwarf accreting from a high-entropy helium-core white dwarf. Given its mass, the accreting carbon–oxygen white dwarf is poised to trigger a typical-luminosity SN Ia within a few million years, to evolve into a stably transferring AM Canum Venaticorum (or AM CVn) system, or undergo a merger into a massive white dwarf. ATLAS J1138-5139 provides a rare opportunity to calibrate binary evolution models by directly comparing observed orbital parameters and mass-transfer rates closer to merger than any known SN Ia progenitor. Its compact orbit ensures detectability by LISA, demonstrating the potential of millihertz gravitational-wave observatories to reveal a population of SN Ia progenitors on a Galactic scale, paving the way for multimessenger studies offering insights into the origins of these cosmologically significant explosions.