skip to main content

Title: Skye: A Differentiable Equation of State

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.

; ; ; ;
Publication Date:
Journal Name:
The Astrophysical Journal
Page Range or eLocation-ID:
Article No. 72
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

    We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings. Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker–Planck limit of the Boltzmann equation. Our model quantifies the dominant gravitational physics, combining both strong and weak single–binary interactions, which drives the time evolution of binary orbital parameter distributions in dense stellar environments. We focus in particular the distributions of binary orbital energies and eccentricities. We find a novel steady-state distribution of binary eccentricities, featuring strong depletions of both the highest and the lowest eccentricity binaries. In energy space, we compare the predictions of our analytic model to the results of numerical N-body simulations, and find that the agreement is good for the initial conditions considered here. This work is a first step towards the development of a fully self-consistent semi-analytic model for dynamically evolving binary star populations in dense stellar environments due to direct few-body interactions.


    Phase transitions play an important role for the style of mantle convection. While observations and theory agree that a substantial fraction of subducted slabs and rising plumes can move through the whole mantle at present day conditions, this behaviour may have been different throughout Earth’s history. Higher temperatures, such as in the early Earth, cause different phase transitions to be dominant, and also reduce mantle viscosity, favouring a more layered style of convection induced by phase transitions. A period of layered mantle convection in Earth’s past would have significant implications for the secular evolution of the mantle temperature and the mixing of mantle heterogeneities. The transition from layered to whole mantle convection could lead to a period of mantle avalanches associated with a dramatic increase in magmatic activity. Consequently, it is important to accurately model the influence of phase transitions on mantle convection. However, existing numerical methods generally preclude modelling phase transitions that are only present in a particular range of pressures, temperatures or compositions, and they impose an artificial lower limit on the thickness of phase transitions. To overcome these limitations, we have developed a new numerical method that solves the energy equation for entropy instead of temperature. Thismore »technique allows for robust coupling between thermodynamic and geodynamic models and makes it possible to model realistically sharp phase transitions with a wide range of properties and dynamic effects on mantle processes. We demonstrate the utility of our method by applying it in regional and global convection models, investigating the effect of individual phase transitions in the Earth’s mantle with regard to their potential for layering flow. We find that the thickness of the phase transition has a bigger influence on the style of convection than previously thought: with all other parameters being the same, a thin phase transition can induce fully layered convection where a broad phase transition would lead to whole-mantle convection. Our application of the method to convection in the early Earth illustrates that endothermic phase transitions may have induced layering for higher mantle temperatures in the Earth’s past.

    « less
  4. 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
  5. Abstract

    We describe the public release of the Cluster Monte Carlo (CMC) code, a parallel, star-by-starN-body code for modeling dense star clusters.CMCtreats collisional stellar dynamics using Hénon’s method, where the cumulative effect of many two-body encounters is statistically reproduced as a single effective encounter between nearest-neighbor particles on a relaxation timescale. The star-by-star approach allows for the inclusion of additional physics, including strong gravitational three- and four-body encounters, two-body tidal and gravitational-wave captures, mass loss in arbitrary galactic tidal fields, and stellar evolution for both single and binary stars. The public release ofCMCis pinned directly to theCOSMICpopulation synthesis code, allowing dynamical star cluster simulations and population synthesis studies to be performed using identical assumptions about the stellar physics and initial conditions. As a demonstration, we present two examples of star cluster modeling: first, we perform the largest (N= 108) star-by-starN-body simulation of a Plummer sphere evolving to core collapse, reproducing the expected self-similar density profile over more than 15 orders of magnitude; second, we generate realistic models for typical globular clusters, and we show that their dynamical evolution can produce significant numbers of black hole mergers with masses greater than those produced from isolated binary evolution (such as GW190521, amore »recently reported merger with component masses in the pulsational pair-instability mass gap).

    « less