Investigation of Important Weak Interaction Nuclei in Presupernova Evolution

Abstract

Abstract This project aims to investigate the most important weak interaction nuclei in the presupernova evolution of massive stars. To achieve this goal, an ensemble containing 728 nuclei in the mass range of A = 1–100 was considered. We computed the mass fractions of these nuclei using Saha’s equation for predetermined values of T, ρ, and Y e and assuming nuclear statistical equilibrium. The nuclear partition functions were obtained using a newly introduced recipe where excited states, up to 10 MeV, were treated as discrete. The weak interaction rates (electron capture (ec) and β-decay (bd)) were calculated in a totally microscopic fashion using the proton–neutron quasiparticle random phase approximation model and without assuming the Brink–Axel hypothesis. The calculated rates were coupled with the computed mass fractions to investigate the time rate of change of lepton to baryon fraction of the stellar matter. We compare our results with the previous calculations reported in the literature. Noticeable differences up to orders of magnitude are reported with previous calculations. These differences may influence the evolution of the star in the later stages of presupernova. We present a list of the top 50 ec and bd nuclei, which have the largest effect on Y e for conditions after silicon core burning. The competition between the ec and bd rates in the stellar core was investigated and it was found that Y e = 0.424–0.455 is the interval where the bd results are bigger than the ec rates.