Skye: A Differentiable Equation of State

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, non-linear 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.

Accepted in ApJ. 27 pages, 19 figures

References in corpus (15)

• The NumPy array: a structure for efficient numerical computation
• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
• Updated Electron-Conduction Opacities: The Impact on Low-Mass Stellar Models
• Relativistic effective interaction for nuclei, giant resonances, and neutron stars
• Nuclear equation of state for core-collapse supernova simulations with realistic nuclear forces
• Coulomb tunneling for fusion reactions in dense matter: Path integral Monte Carlo versus mean field
• Ab initio equations of state for hydrogen (H-REOS.3) and helium (He-REOS.3) and their implications for the interior of Brown Dwarfs
• A New Open-Source Nuclear Equation of State Framework based on the Liquid-Drop Model with Skyrme Interaction
• Phase diagram of hydrogen and a hydrogen-helium mixture at planetary conditions by Quantum Monte Carlo simulations
• Equation of state of classical Coulomb plasma mixtures
• Toward Precision Cosmochronology: A New C/O Phase Diagram for White Dwarfs
• Direct Evaluation of the Phase Diagrams of Dense Multicomponent Plasmas by Integration of the Clapeyron Equations
• Quantum ion thermodynamics in liquid interiors of white dwarfs
• Interpolation of equation-of-state data

Cited by in corpus (77)

• Modules for Experiments in Stellar Astrophysics (MESA): Time-Dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure
• Current challenges in the physics of white dwarf stars
• Axion Astrophysics
• Carbon-Oxygen Phase Separation in MESA White Dwarf Models
• White dwarfs as Physics laboratories: lights and shadows
• Giant planet engulfment by evolved giant stars: light curves, asteroseismology, and survivability
• An Atlas of Convection in Main-Sequence Stars
• The effects of surface fossil magnetic fields on massive star evolution: IV. Grids of models at Solar, LMC, and SMC metallicities
• Mass Transfer and Stellar Evolution of the White Dwarfs in AM CVn Binaries
• Reversing the verdict: Cataclysmic variables could be the dominant progenitors of AM CVn binaries after all
• Period bouncers as detached magnetic cataclysmic variables
• Magnetic field breakout from white dwarf crystallization dynamos
• Modeling contact binaries, II. The effect of energy transfer
• The Asteroseismic Imprints of Mass Transfer: A Case Study of a Binary Mass Gainer in the SPB Instability Strip
• A theory of mass transfer in binary stars
• A Short Intense Dynamo at the Onset of Crystallization in White Dwarfs
• Isochrone fitting to the open cluster M67 in the era of Gaia and improved model physics
• The implications of large binding energies of massive stripped core collapse supernova progenitors on the explosion mechanism
• Efficient Gravitational Wave Template Bank Generation with Differentiable Waveforms
• Observational predictions for Thorne-Żytkow objects
• $Abinitio$ thermodynamics of one-component plasma for astrophysics of white dwarfs and neutron stars
• Orbital decay in an accreting and eclipsing 13.7 minute orbital period binary with a luminous donor
• Magnetic field evolution for crystallization-driven dynamos in C/O white dwarfs
• Characterizing Observed Extra Mixing Trends in Red Giants using the Reduced Density Ratio from Thermohaline Models
• Long-duration Gamma-ray Burst Progenitors and Magnetar Formation
• Fast and Slow Crystallization-driven Convection in White Dwarfs
• Cooling Models for the Most Massive White Dwarfs
• The jittering jets explosion mechanism (JJEM) in electron capture supernovae
• Further evidence of saturated, boosted, and disrupted magnetic braking from evolutionary tracks of cataclysmic variables
• Theoretical investigation of the occurrence of tidally excited oscillations in massive eccentric binary systems
• The formation of the magnetic symbiotic star FN Sgr
• An upper limit on the spins of merging binary black holes formed through binary evolution
• Modeling contact binaries, III. Properties of a population of close, massive binaries
• The Influence of Disk Composition on the Evolution of Stars in the Disks of Active Galactic Nuclei
• The Impacts of Neutron-Star Structure and Base Heating on Type I X-Ray Bursts and Code Comparison
• Stellar Neutrino Emission Across The Mass-Metallicity Plane
• Supernova Shocks Cannot Explain the Inflated State of Hypervelocity Runaways from White Dwarf Binaries
• Phase diagrams of binary ionic mixtures and white dwarf cooling
• Collision-induced mass loss and mass gain on an extremely massive star. An analytical approach and a static proto-globular cluster test-case
• On the Anti-Correlation between Duration and Redshift in Gamma-ray Bursts
• The Secular Periodic Evolution of X-ray Quasi-periodic Eruptions Driven by Star-disc Collisions
• Asteroseismic Inversions for Internal Sound Speed Profiles of Main-sequence Stars with Radiative Cores
• Post-main sequence thermal evolution of planetesimals
• HIP 15429: A newborn Be star on an eccentric binary orbit
• The core degenerate scenario for the type Ia supernova SN 2020eyj
• Resolution of a paradox: SDSS J1257+5428 can be explained as a descendant of a cataclysmic variable with an evolved donor
• Terminating a common envelope jets supernova impostor event with a super-Eddington blue supergiant
• On excess entropy and latent heat in crystallizing white dwarfs
• Reinterpreting the Polluted White Dwarf SDSS J122859.93+104032.9 in Light of Thermohaline Mixing Models: More Polluting Material from a Larger Orbiting Solid Body
• Suggested magnetic braking prescription derived from field complexity fails to reproduce the cataclysmic variable orbital period gap
• Evolving massive stars to core collapse with GENEC: Extension of equation of state, opacities and effective nuclear network
• Asteroseismic Structure Inversions of Main-Sequence Solar-like Oscillators with Convective Cores
• The structure and evolution of a high-mass stellar merger in the Hertzsprung gap
• The efficiency of mixed modes for angular momentum transport
• Revisiting the extremely long-period cataclysmic variables V479 Andromedae and V1082 Sagitarii
• Evolutionary models for the Very Massive Stars in the R136 cluster of 30 Doradus in the Large Magellanic Cloud
• Three-component Phase Separation for Ultramassive White Dwarf Models
• Seismology and diffusion of ultramassive white dwarf magnetic fields
• Toward a Comprehensive Grid of Cepheid Models with MESA. III. Evolutionary and Pulsation Relations for Models with Core and Envelope Overshooting
• Metallicity Ceiling in Quasars from Recycled Stellar Winds
• Evolution of a Long-Lived Deep-Seated Main-Sequence Magnetic Field During White Dwarf Cooling
• Electrostatic energy of solid binary ionic mixtures
• From Main Sequence Binary to Blast: MESA Modeling of the Double-Detonation Progenitor PTF1~J2238+7430
• Toward a Comprehensive Grid of Cepheid Models with MESA II. Impact of Physical and Numerical Assumptions on Elemental Abundances
• Magnetic dynamos powered by white dwarf superficial convection
• A study of the electrostatic properties of the interiors of low-mass stars: Possible implications for the observed rotational properties
• Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy
• Three-dimensional equation of state extension of quark matter in Fermi-liquid theory
• Slow convection and fast rotation in crystallization-driven white dwarf dynamos
• Toward a Comprehensive Grid of Cepheid Models with MESA. IV. Modest Effects of Rotation on Blue Loops
• A Theoretical Study of the Structure and Elemental Abundances of HD 20794
• Gravitational-Wave Signatures of Massive Black Hole Formation
• Continuous Gravitational Waves from Supersoft X-ray Sources: Promising Targets for deci-Hz Detectors
• Resonance locking and tidal evolution in rotating γ-Doradus binaries
• Constraints on the Axion-Photon Coupling Using Stellar Modelling
• Evolution of the ZTF SLRN-2020 star-planet merger
• Ensemble seismic study of the properties of the core of Red Clump stars