An iterative method for simultaneously computing the synthesis of light and heavy elements

We propose an iterative algorithm for computing the synthesis of heavy elements through the rapid capture of neutrons (r-process) and, at sufficiently high temperatures, protons by simultaneously using two distinct computer codes. One of the codes describes the kinetics of nuclear reactions between light and intermediate chemical elements, which are the source of free neutrons and protons used by the second code to synthesize heavy elements from seed nuclides (isotopes near the iron peak of the cosmic abundance curve). The two codes interact through the neutron and proton reaction channels. We demonstrate the efficiency of our method with an example of the nucleosynthesis in a supernova’s helium shell triggered by the evaporation of neutrons and protons from α particles exposed to the neutrino flux from a collapsing stellar core. In this case, three or four iterations are enough to obtain an almost exact self-consistent solution.