Nuclear Reaction Networks


Astronomy research
  Software Infrastructure:
     My instruments
  White dwarf pulsations:
     12C(α,γ) & overshooting
     Probe of 12C(α,γ)16O
     Impact of 22Ne
     Impact of ν cooling
     Variable white dwarfs
     MC reaction rates
  White dwarf supernova:
     Stable nickel production
     Remnant metallicities
     Colliding white dwarfs
     Merging white dwarfs
     Ignition conditions
     Metallicity effects
     Central density effects
     Detonation density
     Tracer particle burning
     Subsonic burning fronts
     Supersonic fronts
     W7 profiles
  Massive stars:
     Pop III with HST/JWST
     Rotating progenitors
     3D evolution to collapse
     MC reaction rates
     Pre-SN variations
  Massive star supernova:
     Yields of radionuclides
     26Al & 60Fe
     44Ti, 60Co & 56Ni
     SN 1987A light curve
     Constraints on Ni/Fe
     An r-process
     Effects of 12C +12C
  Neutron Stars and Black Holes:
     Black Hole spectrum
     Mass Gap with LVK
     Compact object IMF
     He burn neutron stars
  Neutrino Emission:
     Neutrino emission from stars
     Identifying the Pre-SN
     Neutrino HR diagram
     Pre-SN Beta Processes
     Pre-SN neutrinos
     Hypatia catalog
     SAGB stars
     Nugrid Yields I
     He shell convection
     BBFH at 40 years
     γ-rays within 100 Mpc
     Iron Pseudocarbynes
  Pre-Solar Grains:
     C-rich presolar grains
     SiC Type U/C grains
     Grains from massive stars
     Placing the Sun
     SiC Presolar grains
  Chemical Evolution:
     Radionuclides in 2020s
     Zone models H to Zn
     Mixing ejecta
  Thermodynamics, Opacities & Networks
     Radiative Opacity
     Skye EOS
     Helm EOS
     Five EOSs
     Equations of State
     12C(α,γ)16O Rate
     Proton-rich NSE
     Reaction networks
     Bayesian reaction rates
  Verification Problems:
     Validating an astro code
Software Instruments
AAS Journals
   2024 AAS YouTube
   2024 AAS Peer Review Workshops

2024 ASU Energy in Everyday Life
2024 MESA Classroom
Outreach and Education Materials

Other Stuff:
   Bicycle Adventures

Contact: F.X.Timmes
my one page vitae,
full vitae,
research statement, and
teaching statement.
Integration of Nuclear Reaction Networks for Stellar Hydrodynamics (2000)

In this article, methods for solving the stiff system of ordinary differential equations that constitute nuclear reaction networks are surveyed. Three semi-implicit time integration algorithms are examined; a traditional first-order-accurate Euler method, a fourth-order-accurate Kaps-Rentrop method, and a variable-order Bader-Deuflhard method. These three integration methods are coupled to eight different linear algebra packages. Four of the linear algebra packages operate on dense matrices (LAPACK, LUDCMP, LEQS, GIFT), three of them are designed for the direct solution of sparse matrices (MA28, UMFPACK, Y12M), and one uses an iterative method for sparse matrices (BiCG). The scaling properties and behavior of the 24 combinations (3 time integration methods times 8 linear algebra packages) are analyzed by running each combination on seven different nuclear reaction networks.

These reaction networks range from a hardwired 13 isotope $\alpha$-chain and heavy-ion reaction network, which is suitable for most multidimensional simulations of stellar phenomena, to a 489 isotope reaction network, which is suitable for determining the yields of isotopes lighter than technetium in spherically symmetric models of Type II supernovae. Each of the time integration methods and linear algebra packages are capable of generating accurate results, but the efficiency of the various methods -- evaluated across several different machine architectures and compiler options -- differ dramatically. If the execution speed of reaction networks that contain less than about 50 isotopes is an overriding concern, then the variable-order Bader-Deuflhard time integration method coupled with routines generated from the GIFT matrix package or LAPACK with vendor-optimized BLAS routines is a good choice. If the amount of storage needed for any reaction network is a concern, then any of the sparse matrix packages will reduce the storage costs by 70%-90%. When a balance between accuracy, overall efficiency, and ease of use is desirable, then the variable-order Bader-Deuflhard time integration method coupled with the MA28 sparse matrix package is a strong choice.

Kaps-Rentrop CPU time
Bader-Deuflard CPU time

An Inexpensive Nuclear Energy Generation Network for Stellar Hydrodynamics (2000)

In this article, we compare the nuclear energy generation rate and abundance levels given by an $\alpha$-chain nuclear reaction network that contains only seven isotopes with a standard 13 isotope α-chain reaction network. The energy generation rate of these two small networks are also compared to the energy generation rate given by a 489 isotope reaction network with weak reactions turned on and off. The comparison between the seven isotope and α-chain reaction networks indicate the extent to which one can be replaced by the other, and the comparison with the 489 isotope reaction network roughly indicates under what physical conditions it is safe to use the seven isotope and $\alpha$-chain reaction networks.

The seven isotope reaction network reproduces the nuclear energy generation rate of the standard $\alpha$-chain reaction network to within 30%, but often much better, during hydrostatic and explosive helium, carbon, and oxygen burning. It will also provide energy generation rates within 30% of an $\alpha$-chain reaction network for silicon burning at densities less than 10$^{7}$ g cm$^{-3}$. Provided there remains an equal number of protons and neutrons (Y$_e$ = 0.5) over the course of the evolution, and that flows through $\alpha$ particle channels dominate, then both of the small reaction networks return energy generation rates that are compatible with the energy generation rate returned by the 489 reaction network. If Y$_e$ is significantly different from 0.5, or if there are substantial flows through neutron and protons channels, then it is not generally safe to employ any $\alpha$-chain based reaction network. The relative accuracy of the 7 isotope reaction network, combined with its reduction in the computational cost, suggest that it is a suitable replacement for $\alpha$-chain reaction networks for parameter space surveys of a wide class of multidimensional stellar models.

Energy generation rate comparision
Ratio of energy generation rates
Isotope evolution

Nuclear Reaction Network Software Instruments

Open source software instruments are avaliable here.