Astronomy research
  Software Infrastructure:
     My codes
  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
     26Al & 60Fe
     44Ti, 60Co & 56Ni
     Yields of radionuclides
     Effects of 12C +12C
     SN 1987A light curve
     Constraints on Ni/Fe ratios
     An r-process
     Compact object IMF
     Pop III JWST
     Pre-SN 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:
     Hypatia catalog
     Zone models H to Zn
     Mixing ejecta
     γ-rays within 100 Mpc
  Thermodynamics & Networks
     Stellar EOS
     Reaction networks
     Proton-rich NSE
     Bayesian reaction rates
  Verification Problems:
     Validating an astro code
Software instruments
Bicycle adventures

AAS Journals
2017 MESA Marketplace
2017 MESA Summer School
2017 ASU+EdX AST111x
Teaching materials
Education and Public Outreach

Contact: F.X.Timmes
my one page vitae,
full vitae,
research statement, and
teaching statement.
The Accuracy, Consistency, and Speed of Five Equations of State for Stellar Hydrodynamics

In this paper, 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 equation-of-state 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 equation-of-state 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.

Pressure differences
Thermodynamic consistency

The Accuracy, Consistency, and Speed Of An Electron–Positron Equation Of State
Based On Table Interpolation Of The Helmholtz Free Energy

In this paper, an electron-positron 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.

Specific internal energy

Equation of State Codes
Open-source codes are avaliable from this link and this link