Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStatistical Methods for Thermonuclear Reaction Rates and Nucleosynthesis Simulations
115
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Rigorous statistical methods for estimating thermonuclear reaction rates and nucleosynthesis are becoming increasingly established in nuclear astrophysics. The main challenge being faced is that experimental reaction rates are highly complex quantities derived from a multitude of different measured nuclear parameters (e.g., astrophysical S-factors, resonance energies and strengths, particle and gamma-ray partial widths). We discuss the application of the Monte Carlo method to two distinct, but related, questions. First, given a set of measured nuclear parameters, how can one best estimate the resulting thermonuclear reaction rates and associated uncertainties? Second, given a set of appropriate reaction rates, how can one best estimate the abundances from nucleosynthesis (i.e., reaction network) calculations? The techniques described here provide probability density functions that can be used to derive statistically meaningful reaction rates and final abundances for any desired coverage probability. Examples are given for applications to s-process neutron sources, core-collapse supernovae, classical novae, and big bang nucleosynthesis.
rate research
Read More
A method based on Monte Carlo techniques is presented for evaluating thermonuclear reaction rates. We begin by reviewing commonly applied procedures and point out that reaction rates that have been reported up to now in the literature have no rigorous statistical meaning. Subsequently, we associate each nuclear physics quantity entering in the calculation of reaction rates with a specific probability density function, including Gaussian, lognormal and chi-squared distributions. Based on these probability density functions the total reaction rate is randomly sampled many times until the required statistical precision is achieved. This procedure results in a median (Monte Carlo) rate which agrees under certain conditions with the commonly reported recommended classical rate. In addition, we present at each temperature a low rate and a high rate, corresponding to the 0.16 and 0.84 quantiles of the cumulative reaction rate distribution. These quantities are in general different from the statistically meaningless minimum (or lower limit) and maximum (or upper limit) reaction rates which are commonly reported. Furthermore, we approximate the output reaction rate probability density function by a lognormal distribution and present, at each temperature, the lognormal parameters miu and sigma. The values of these quantities will be crucial for future Monte Carlo nucleosynthesis studies. Our new reaction rates, appropriate for bare nuclei in the laboratory, are tabulated in the second paper of this series (Paper II). The nuclear physics input used to derive our reaction rates is presented in the third paper of this series (Paper III). In the fourth paper of this series (Paper IV) we compare our new reaction rates to previous results.
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.
Background: Exotic non-spherical nuclear pasta shapes are expected in nuclear matter at just below saturation density because of competition between short range nuclear attraction and long range Coulomb repulsion. Purpose: We explore the impact of nuclear pasta on nucleosynthesis, during neutron star mergers, as cold dense nuclear matter is ejected and decompressed. Methods: We perform classical molecular dynamics simulations with 51200 and 409600 nucleons, that are run on GPUs. We expand our simulation region to decompress systems from an initial density of 0.080 fm^{-3} down to 0.00125 fm^{-3}. We study proton fractions of Y_P=0.05, 0.10, 0.20, 0.30, and 0.40 at T =0.5, 0.75, and 1.0 MeV. We calculate the composition of the resulting systems using a cluster algorithm. Results: We find final compositions that are in good agreement with nuclear statistical equilibrium models for temperatures of 0.75 and 1 MeV. However, for proton fractions greater than Y_P=0.2 at a temperature of T = 0.5 MeV, the MD simulations produce non-equilibrium results with large rod-like nuclei. Conclusions: Our MD model is valid at higher densities than simple nuclear statistical equilibrium models and may help determine the initial temperatures and proton fractions of matter ejected in mergers.
Radiative capture reactions play a crucial role in stellar nucleosynthesis but have proved challenging to determine experimentally. In particular, the large uncertainty ($sim$100%) in the measured rate of the $^{12}$C$(alpha,gamma)^{16}$O reaction is the largest source of uncertainty in any stellar evolution model. With development of new high current energy-recovery linear accelerators (ERLs) and high density gas targets, measurement of the $^{16}$O$(e,e^prime alpha)^{12}$C reaction close to threshold using detailed balance opens up a new approach to determine the $^{12}$C$(alpha,gamma)^{16}$O reaction rate with significantly increased precision ($<$20%). We present the formalism to relate photo- and electro-disintegration reactions and consider the design of an optimal experiment to deliver increased precision. Once the new ERLs come online, an experiment to validate the new approach we propose should be carried out. This new approach has broad applicability to radiative capture reactions in astrophysics.
The thermonuclear rate of the 50Fe(p,gamma)51Co reaction in the Type I X-ray bursts (XRBs) temperature range has been reevaluated based on a recent precise mass measurement at CSRe lanzhou, where the proton separation energy Sp=142+/-77 keV has been determined firstly for the 51Co nucleus. Comparing to the previous theoretical predictions, the experimental Sp value has much smaller uncertainty. Based on the nuclear shell model and mirror nuclear structure information, we have calculated two sets of thermonuclear rates for the 50Fe(p,gamma)51Co reaction by utilizing the experimental Sp value. It shows that the statistical-model calculations are not ideally applicable for this reaction primarily because of the low density of low-lying excited states in 51Co. In this work, we recommend that a set of new reaction rate based on the mirror structure of 51Cr should be incorporated in the future astrophysical network calculations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
