No Arabic abstract
Nuclear reaction rates of astrophysical applications are traditionally determined on the basis of Hauser-Feshbach reaction codes. These codes adopt a number of approximations that have never been tested, such as a simplified width fluctuation correction, the neglect of delayed or multiple-particle emission during the electromagnetic decay cascade, or the absence of the pre-equilibrium contribution at increasing incident energies. The reaction code TALYS has been recently updated to estimate the Maxwellian-averaged reaction rates that are of astrophysical relevance. These new developments enable the reaction rates to be calculated with increased accuracy and reliability and the approximations of previous codes to be investigated. The TALYS predictions for the thermonuclear rates of relevance to astrophysics are detailed and compared with those derived by widely-used codes for the same nuclear ingredients. It is shown that TALYS predictions may differ significantly from those of previous codes, in particular for nuclei for which no or little nuclear data is available. The pre-equilibrium process is shown to influence the astrophysics rates of exotic neutron-rich nuclei significantly. For the first time, the Maxwellian-averaged (n,2n) reaction rate is calculated for all nuclei and its competition with the radiative capture rate is discussed. The TALYS code provides a new tool to estimate all nuclear reaction rates of relevance to astrophysics with improved accuracy and reliability.
Nuclear Astrophysics requires the knowledge of reaction rates over a wide range of nuclei and temperatures. In recent calculations the nuclear level density - as an important ingredient to the statistical model (Hauser-Feshbach) - has shown the highest uncertainties. In a back-shifted Fermi-gas formalism utilizing an energy-dependent level density parameter and employing microscopic corrections from a recent FRDM mass formula, we obtain a highly improved fit to experimental level densities. The resulting level density is used for determining criteria for the applicability of the statistical model on neutron-induced reactions.
The majority of nuclear reactions in astrophysics involve unstable nuclei which are not fully accessible by experiments yet. Therefore, there is high demand for reliable predictions of cross sections and reaction rates by theoretical means. The majority of reactions can be treated in the framework of the statistical model (Hauser-Feshbach). The global parametrizations of the nuclear properties needed for predictions far off stability probe our understanding of the strong force and take it to its limit. The sensitivity of astrophysical scenarios to nuclear inputs is illustrated in the framework of a detailed nucleosynthesis study in type II supernovae. Abundances resulting from calculations in the same explosion model with two different sets of reaction rates are compared. Key reactions and required nuclear information are identified.
The study of d+d reactions is of major interest since their reaction rates affect the predicted abundances of D, 3He, and 7Li. In particular, recent measurements of primordial D/H ratios call for reduced uncertainties in the theoretical abundances predicted by big bang nucleosynthesis (BBN). Different authors have studied reactions involved in BBN by incorporating new experimental data and a careful treatment of systematic and probabilistic uncertainties. To analyze the experimental data, Coc et al. (2015) used results of ab initio models for the theoretical calculation of the energy dependence of S-factors in conjunction with traditional statistical methods based on Chi-2 minimization. Bayesian methods have now spread to many scientific fields and provide numerous advantages in data analysis. Astrophysical S-factors and reaction rates using Bayesian statistics were calculated by Iliadis et al. (2016). Here we present a similar analysis for two d+d reactions, d(d,n)3He and d(d,p)3H, that has been translated into a total decrease of the predicted D/H value by 0.16%.
The $ u p$ process appears in proton-rich, hot matter which is expanding in a neutrino wind and may be realised in explosive environments such as core-collapse supernovae or in outflows from accretion disks. The impact of uncertainties in nuclear reaction cross sections on the finally produced abundances has been studied by applying Monte Carlo variation of all astrophysical reaction rates in a large reaction network. As the detailed astrophysical conditions of the $ u p$ process still are unknown, a parameter study was performed, with 23 trajectories covering a large range of entropies and $Y_mathrm{e}$. The resulting abundance uncertainties are given for each trajectory. The $ u p$ process has been speculated to contribute to the light $p$ nuclides but it was not possible so far to reproduce the solar isotope ratios. It is found that it is possible to reproduce the solar $^{92}$Mo/$^{94}$Mo abundance ratio within nuclear uncertainties, even within a single trajectory. The solar values of the abundances in the Kr-Sr region relative to the Mo region, however, cannot be achieved within a single trajectory. They may still be obtained from a weighted superposition of different trajectories, though, depending on the actual conditions in the production site. For a stronger constraint of the required conditions, it would be necessary to reduce the uncertainties in the 3$alpha$ and $^{56}$Ni(n,p)$^{56}$Co rates at temperatures $T>3$ GK.
Nuclear reaction rate ($lambda$) is a significant factor in the process of nucleosynthesis. A multi-layer directed-weighted nuclear reaction network in which the reaction rate as the weight, and neutron, proton, $^4$He and the remainder nuclei as the criterion for different reaction-layers is for the first time built based on all thermonuclear reactions in the JINA REACLIB database. Our results show that with the increase of the stellar temperature ($T_{9}$), the distribution of nuclear reaction rates on the $R$-layer network demonstrates a transition from unimodal to bimodal distributions. Nuclei on the $R$-layer in the region of $lambda = [1,2.5times10^{1}]$ have a more complicated out-going degree distribution than the one in the region of $lambda = [10^{11},10^{13}]$, and the number of involved nuclei at $T_{9} = 1$ is very different from the one at $T_{9} = 3$. The redundant nuclei in the region of $lambda = [1, 2.5times10^{1}]$ at $T_{9} = 3$ prefer $(gamma,p)$ and $({gamma,alpha})$ reactions to the ones at $T_{9}=1$, which produce nuclei around the $beta$ stable line. This work offers a novel way to the big-data analysis on nuclear reaction network at stellar temperatures.