*
Cococubed.com


Thirteen Isotope Alpha-Chain

Home

Teaching materials
Astronomy research
Astronomy codes
... Stellar equation of states
... EOS with ionization
... EOS for supernovae
... Chemical potentials
... Stellar atmospheres
... Voigt Function
... Polytropic stars
... Cold white dwarfs
... Hotter white dwarfs
... Cold neutron stars
... Stellar opacities
... Neutrino energy loss rates
... Ephemeris routines
... Fermi-Dirac functions
... Galactic chemical evolution

... Nuclear reaction networks
... Nuclear statistical equilibrium
... Laminar deflagrations
... CJ detonations
... ZND detonations
... Fitting to conic sections
... Unusual linear algebra
... Derivatives on uneven grids
... Pentadiagonal solver
... Quadratics, Cubics, Quartics
... Supernova light curves
... Exact Riemann solutions
... 1D PPM Hydrodynamics
... Verification problems

... EZ stellar evolution
... FLASH code
... Mesa code
Astronomy talks
Astronomy images
Outreach
Family album
Bicycle adventures
Artwork


Contact us:
J.D. Maldonado
F.X.Timmes, my vitae

A definition of what I mean by an alpha chain reaction network is prudent. A 'strict' alpha-chain reaction network is only composed of (alpha,gamma) and (gamma,alpha) links among the 13 isotopes 4He, 12C, 16O, 20Ne, 24Mg, 28Si, 32S, 36Ar, 40Ca, 44Ti, 48Cr, 52Fe, and 56Ni. It is essential, however, to include (alpha,p)(p,gamma) and (gamma,p)(p,alpha) links in order to obtain reasonably accurate energy generation rates and abundance levels when the temperature exceeds 2.5e9 K. At these elevated temperatures the flows through the (alpha,p)(p,gamma) sequences are faster than the flows through the (alpha,gamma) channels.

An (alpha,p)(p,gamma) sequence is, effectively, an (alpha,gamma) reaction through an intermediate isotope. In this alpha-chain reaction network, I include 8 (alpha,p)(p,gamma) sequences plus the corresponding inverse sequences by assuming steady-state proton flows through the intermediate isotopes 27Al, 31P, 35Cl, 39K, 43Sc, 47V, 51Mn, and 55Co . This strategy permits inclusion of (alpha,p)(p,gamma) sequences without explicitly evolving the proton or intermediate isotope abundances. Thus, this alpha-chain reaction network includes not just (alpha,gamma) and (gamma,alpha) links, but also links through the (alpha,p)(p,gamma) and (gamma,p)(p,alpha) sequences.

This 13-isotope alpha chain plus heavy-ion (12C+12C, 12C+16O, 16O+16O) reaction network is suitable for many multi-dimensional simulations of stellar phenomena where having a reasonably accurate energy generation rate is of primary concern.


image
mild c+o burn
image
13 isotope jacobian


public_aprox13.f
plus these include files:
burn_common.dek
cjdet.dek
const.dek
implno.dek
network.dek
sparse_matrix.dek
tfactors.dek
timers.dek
vector_eos.dek
plus a 13 MB eos data file :
helm_table.dat


 



Please cite the relevant references if you publish a piece of work that use these codes, pieces of these codes, or modified versions of them. If you're nice, offer co-authorship of the publication. At best, you'll love these programs so much that you'll send great wads of cash to me.