

Home Astronomy research Software instruments Stellar equation of states EOS with ionization EOS for supernovae Chemical potentials Stellar atmospheres Voigt Function Jeans escape Polytropic stars Cold white dwarfs Adiabatic white dwarfs Cold neutron stars Stellar opacities Neutrino energy loss rates Ephemeris routines FermiDirac functions Polyhedra volume Plane  cube intersection Coating an ellipsoid 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 Hydrodynamic test cases Galactic chemical evolution Universal twobody problem Circular and elliptical 3 body The pendulum Phyllotaxis MESA MESAWeb FLASH Zingale's software Brown's dStar GR1D code Iliadis' STARLIB database Herwig's NuGRID Meyer's NetNuc Presentations Illustrations cococubed YouTube Bicycle adventures Public Outreach Education materials 2023 ASU Solar Systems Astronomy 2023 ASU Energy in Everyday Life AAS Journals 2023 AAS YouTube 2023 AAS Peer Review Workshops 2023 MESA VI 2023 MESA Marketplace 2023 MESA Classroom 2023 Neutrino Emission from Stars 2023 White Dwarfs & ^{12}C(α,γ)^{16}O 2022 Earendel, A Highly Magnified Star 2022 Black Hole Mass Spectrum 2022 MESA in Don't Look Up 2021 Bill Paxton, Tinsley Prize Contact: F.X.Timmes my one page vitae, full vitae, research statement, and teaching statement. 
Before using these reaction networks you should probably glance at my method of madness, Raph Hix's & Brad Meyer's excellent article, Brad Meyer's annual review article, George Wallerstein's review of modern physics article, and this National Nuclear Physics Summer School lectures on reaction networks. There is a certain irremovable complexity associated with stiff systems of ordinary differential equations $$ \dot {{\bf y}} = {\bf f} \ ({\bf y}) \label{eq1} \tag{1} $$ when the right hand side is a complicated function, the Jacobian matrix $\tilde{{\bf J}}$ is sparse, and one wants a highquality time integration. The tools below use an analytical Jacobian, a variableorder BaderDeuflhard integration method, and MA28 sparse linear algebra. The reaction network and thermodynamics are integrated simultaneously. That is, they are fully coupled. Hydrostatic, onestep, adiabatic expansion, selfheating at constant density, selfheating through constant pressure, and arbitrary thermodynamic histories are currently supported. These reaction networks are a snapshot of my current research efforts. If you want to put these reaction networks in a stellar evolution or hydrodynamics software instrument, and/or you want the networks to execute as efficiently as possible, feel free to contact me.



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. Offer coauthorship as appropriate. 
