

Home Astronomy research Software Infrastructure: MESA FLASH STARLIB MESAWeb starkillerastro My instruments White dwarf supernova: Remnant metallicities Colliding white dwarfs Merging white dwarfs Ignition conditions Metallicity effects Central density effects Detonation density effects Tracer particle burning Subsonic burning fronts Supersonic burning fronts W7 profiles Massive star supernova: Rotating progenitors 3D evolution ^{26}Al & ^{60}Fe ^{44}Ti, ^{60}Co & ^{56}Ni Yields of radionuclides Effects of ^{12}C +^{12}C SN 1987A light curve Constraints on Ni/Fe ratios An rprocess Compact object IMF Stars: Neutrino HR diagram Pulsating white dwarfs Pop III with JWST Monte Carlo massive stars Neutrinos from preSN PreSN variations Monte Carlo white dwarfs SAGB stars Classical novae He shell convection Presolar grains He burn on neutron stars BBFH at 40 years Chemical Evolution: Iron Pseudocarbynes Radionuclides in the 2020s Hypatia catalog Zone models H to Zn Mixing ejecta γrays within 100 Mpc Thermodynamics & Networks Stellar EOS 12C(α,γ)16O Rate Protonrich NSE Reaction networks Bayesian reaction rates Verification Problems: Validating an astro code SuOlson Cog8 Mader RMTV Sedov Noh Software instruments Presentations Illustrations cococubed YouTube Bicycle adventures Public Outreach Education materials AAS Journals AAS YouTube 2020 Celebration of Margaret Burbidge 2020 Digital Infrastructure 2021 MESA Marketplace 2021 MESA Summer School 2021 ASU Solar Systems 2021 ASU Energy in Everyday Life Contact: F.X.Timmes my one page vitae, full vitae, research statement, and teaching statement. 
Skye: A Differentiable Equation of State (2021)
Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. In this article we present Skye, a new equation of state for fullyionized 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 selfconsistent 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. The Skye EOS = an improved Helmholtz EOS for the noninteracting parts + an improved Potekhin & Chabrier EOS for the Coulomb plama parts + autodifferentiation. Its the bees knees for ionized plasmas as of 2021. Skye is avaliable at https://github.com/adamjermyn/Skye.
The Accuracy, Consistency, and Speed Of An Electron–Positron Equation Of State Based On Table Interpolation Of The Helmholtz Free Energy (2000) In this article, an electronpositron equation of state based on table interpolation of the Helmholtz free energy is developed and analyzed. The interpolation scheme guarantees perfect thermodynamic consistency, independent of the interpolating function. The choice of a biquintic Hermite polynomial as the interpolating function results in accurately reproducing the underlying Helmholtz free energy data in the table, and yields derivatives of the pressure, specific entropy and specific internal energy which are smooth and continuous. The execution speed – evaluated across several different machine architectures, compiler options, and mode of operation – suggest that the Helmholtz equation of state routine is faster than any of the five equation of state routines surveyed by Timmes & Arnett (1999). When an optimal balance of accuracy, thermodynamic consistency, and speed is desirable, then the tabular Helmholtz equation of state is an excellent choice, particularly for multidimensional models of stellar phenomena.
The Accuracy, Consistency, and Speed of Five Equations of State for Stellar Hydrodynamics (1999) In this article, we compare the thermodynamic properties and execution speed of five independent equations of state. A wide range of temperatures, densities, and compositions are considerebdconditions appropriate for modeling the collapse of a cloud of hydrogen gas (or an exploding supernova) to the outer layers of a neutron star. The pressures and specific thermal energies calculated by each equationofstate routine are reasonably accurate (typically 0.1% error or less) and agree remarkably well with each other, despite the different approaches and approximations used in each routine. The derivatives of the pressure and specific thermal energies with respect to the temperature and density generally show less accuracy (typically 1% error or less) and more disagreement with one another. Thermodynamic consistency, as measured by deviations from the appropriate Maxwell relations, shows that the Timmes equation of state and the Nadyozhin equation of state achieve thermodynamic consistency to a high degree of precision. The execution speed of the five equationofstate routinebsevaluated across several different machine architectures, compiler options, and modes of operatiobndiffer dramatically. The Arnett equation of state is the fastest of the five routines, with the Nadyozhin equation of state close behind. Equation of State Software Instruments Opensource codes are avaliable from this link and this link. There are times when a simpler cold fermi gas EOS is a wonderful thing. Such an EOS is in cold_fermi_gas.tbz. One can see this equation of state in action on this cold white dwarf page. 


