Performance Improvements for Nuclear Reaction Network Integration

Abstract

Aims: The aim of this work is to compare the performance of three reaction network integration methods used in stellar nucleosynthesis calculations. These are the Gear's backward differentiation method, Wagoner's method (a 2nd-order Runge-Kutta method), and the Bader-Deuflehard semi-implicit multi-step method. Methods: To investigate the efficiency of each of the integration methods considered here, a test suite of temperature and density versus time profiles is used. This suite provides a range of situations ranging from constant temperature and density to the dramatically varying conditions present in white dwarf mergers, novae, and x-ray bursts. Some of these profiles are obtained separately from full hydrodynamic calculations. The integration efficiencies are investigated with respect to input parameters that constrain the desired accuracy and precision. Results: Gear's backward differentiation method is found to improve accuracy, performance, and stability in integrating nuclear reaction networks. For temperature-density profiles that vary strongly with time, it is found to outperform the Bader-Deuflehard method (although that method is very powerful for more smoothly varying profiles). Wagoner's method, while relatively fast for many scenarios, exhibits hard-to-predict inaccuracies for some choices of integration parameters owing to its lack of error estimations.

13 pages, 12 figures, accepted to Astronomy and Astrophysics (section 15) - corrected units in Figs. 6-11

References in corpus (12)

• The Effects of Thermonuclear Reaction Rate Variations on Nova Nucleosynthesis: A Sensitivity Study
• STARLIB: A Next-Generation Reaction-Rate Library for Nuclear Astrophysics
• Nuclear uncertainties in the NeNa-MgAl cycles and production of 22Na and 26Al during nova outbursts
• The Effects of Variations in Nuclear Processes on Type I X-Ray Burst Nucleosynthesis
• Stellar GADGET: A smooth particle hydrodynamics code for stellar astrophysics and its application to Type Ia supernovae from white dwarf mergers
• Hydrodynamic Models of Type I X-Ray Bursts: Metallicity Effects
• Double white dwarf mergers and elemental surface abundances in extreme helium and R Coronae Borealis stars
• On the Dynamics of Proto-Neutron Star Winds and r-Process Nucleosynthesis
• The Impact of Uncertainties in Reaction Q-values on Nucleosynthesis in Type I X-Ray Bursts
• Smoothed Particle Hydrodynamics simulations of white dwarf collisions and close encounters
• Nucleosynthesis during the Merger of White Dwarfs and the Origin of R Coronae Borealis Stars
• Algebraic Stabilization of Explicit Numerical Integration for Extremely Stiff Reaction Networks

Cited by in corpus (6)

• Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions
• SkyNet: A modular nuclear reaction network library
• The evolution of CNO elements in galaxies
• ${}^{60}$Fe in core-collapse supernovae and prospects for X-ray and $\gamma$-ray detection in supernova remnants
• First models of s process in AGB stars of solar metallicity for the stellar evolutionary code ATON with a novel stable explicit numerical solver
• On the parallelization of stellar evolution codes