Nuclear masses near $N=82$ that influence $r$-process abundances

Nuclear masses are one of the key ingredients of nuclear physics that go into astrophysical simulations of the $r$ process. Nuclear masses effect $r$-process abundances by entering into calculations of Q-values, neutron capture rates, photo-dissociation rates, beta-decay rates, branching ratios and the properties of fission. Most of the thousands of short-lived neutron-rich nuclei which are believed to participate in the $r$ process lack any experimental verification, thus the identification of the most influential nuclei is of paramount importance. We have conducted mass sensitivity studies near the $N=82$ closed shell in the context of a main $r$-process. Our studies take into account how an uncertainty in a single nuclear mass propagates to influence the relevant quantities of neighboring nuclei and finally to $r$-process abundances. We identify influential nuclei in various astrophysical conditions using the FRDM mass model. We show that our conclusions regarding these key nuclei are still retained when a superposition of astrophysical trajectories is considered.

Measurement of the reaction O-17(alpha,n)Ne-20 and its impact on the s process in massive stars

The ratio between the rates of the reactions O-17(alpha,n)Ne-20 and O-17(alpha,gamma)Ne-21 determines whether O-16 is an efficient neutron poison for the s process in massive stars, or if most of the neutrons captured by O-16(n,gamma) are recycled into the stellar environment. This ratio is of particular relevance to constrain the s process yields of fast rotating massive stars at low metallicity. Recent results on the (alpha,gamma) channel have made it necessary to measure the (alpha,n) reaction more precisely and investigate the effect of the new data on s process nucleosynthesis in massive stars. We present a new measurement of the O-17(alpha, n) reaction using a moderating neutron detector. In addition, the (alpha, n_1) channel has been measured independently by observation of the characteristic 1633 keV gamma-transition in Ne-20. The reaction cross section was determined with a simultaneous R-matrix fit to both channels. (alpha,n) and (alpha, gamma) resonance strengths of states lying below the covered energy range were estimated using their known properties from the literature. A new O-17(alpha,n) reaction rate was deduced for the temperature range 0.1 GK to 10 GK. It was found that in He burning conditions the (alpha,gamma) channel is strong enough to compete with the neutron channel. This leads to a less efficient neutron recycling compared to a previous suggestion of a very weak (alpha,gamma) channel. S process calculations using our rates confirm that massive rotating stars do play a significant role in the production of elements up to Sr, but they strongly reduce the s process contribution to heavier elements.

Large collection of astrophysical S-factors and its compact representation

Numerous nuclear reactions in the crust of accreting neutron stars are strongly affected by dense plasma environment. Simulations of superbursts, deep crustal heating and other nuclear burning phenomena in neutron stars require astrophysical S-factors for these reactions (as a function of center-of-mass energy E of colliding nuclei). A large database of S-factors is created for about 5000 non-resonant fusion reactions involving stable and unstable isotopes of Be, B, C, N, O, F, Ne, Na, Mg, and Si. It extends the previous database of about 1000 reactions involving isotopes of C, O, Ne, and Mg. The calculations are performed using the Sao Paulo potential and the barrier penetration formalism. All calculated S-data are parameterized by an analytic model for S(E) proposed before [Phys. Rev. C 82, 044609 (2010)] and further elaborated here. For a given reaction, the present S(E)-model contains three parameters. These parameters are easily interpolated along reactions involving isotopes of the same elements with only seven input parameters, giving an ultracompact, accurate, simple, and uniform database. The S(E) approximation can also be used to estimate theoretical uncertainties of S(E) and nuclear reaction rates in dense matter, as illustrated for the case of the 34Ne+34Ne reaction in the inner crust of an accreting neutron star.

A simple analytic model for astrophysical S-factors

We propose a physically transparent analytic model of astrophysical S-factors as a function of a center-of-mass energy E of colliding nuclei (below and above the Coulomb barrier) for non-resonant fusion reactions. For any given reaction, the S(E)-model contains four parameters [two of which approximate the barrier potential, U(r)]. They are easily interpolated along many reactions involving isotopes of the same elements; they give accurate practical expressions for S(E) with only several input parameters for many reactions. The model reproduces the suppression of S(E) at low energies (of astrophysical importance) due to the shape of the low-r wing of U(r). The model can be used to reconstruct U(r) from computed or measured S(E). For illustration, we parameterize our recent calculations of S(E) (using the Sao Paulo potential and the barrier penetration formalism) for 946 reactions involving stable and unstable isotopes of C, O, Ne, and Mg (with 9 parameters for all reactions involving many isotopes of the same elements, e.g., C+O). In addition, we analyze astrophysically important 12C+12C reaction, compare theoretical models with experimental data, and discuss the problem of interpolating reliably known S(E) values to low energies (E <= 2-3 MeV).

Astrophysical S-factors for fusion reactions involving C, O, Ne and Mg isotopes

Using the Sao Paulo potential and the barrier penetration formalism we have calculated the astrophysical factor S(E) for 946 fusion reactions involving stable and neutron-rich isotopes of C, O, Ne, and Mg for center-of-mass energies E varying from 2 MeV to 18-30 MeV (covering the range below and above the Coulomb barrier). We have parameterized the energy dependence S(E) by an accurate universal 9-parameter analytic expression and present tables of fit parameters for all the reactions. We also discuss the reduced 3-parameter version of our fit which is highly accurate at energies below the Coulomb barrier, and outline the procedure for calculating the reaction rates. The results can be easily converted to thermonuclear or pycnonuclear reaction rates to simulate various nuclear burning phenomena, in particular, stellar burning at high temperatures and nucleosynthesis in high density environments.