X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models.

In a sufficiently hot and dense astrophysical environment the rate of the triple-alpha (3alpha) reaction can increase greatly over the value appropriate for helium burning stars owing to hadronically induced de-excitation of the Hoyle state. In this paper we use a statistical model to evaluate the enhancement as a function of temperature and density. For a density of $10^6 gm cm^{-3}$ enhancements can exceed a factor of one-hundred. In high temperature/density situations, the enhanced 3alpha rate is a better estimate of this rate and should be used in these circumstances. We then examine the effect of these enhancements on production of $^{12}$C in the neutrino wind following a supernova explosion and in an x-ray burster.

Neutrons produced by the carbon fusion reaction 12C(12C,n)23Mg play an important role in stellar nucleosynthesis. However, past studies have shown large discrepancies between experimental data and theory, leading to an uncertain cross section extrapolation at astrophysical energies. We present the first direct measurement that extends deep into the astrophysical energy range along with a new and improved extrapolation technique based on experimental data from the mirror reaction 12C(12C,p)23Na. The new reaction rate has been determined with a well-defined uncertainty that exceeds the precision required by astrophysics models. Using our constrained rate, we find that 12C(12C,n)23Mg is crucial to the production of Na and Al in Pop-III Pair Instability Supernovae. It also plays a non-negligible role in the production of weak s-process elements as well as in the production of the important galactic gamma-ray emitter 60Fe.

We have performed for the first time a comprehensive study of the sensitivity of $r$-process nucleosynthesis to individual nuclear masses across the chart of nuclides. Using the latest version (2012) of the Finite-Range Droplet Model, we consider mass variations of $\pm0.5$ MeV and propagate each mass change to all affected quantities, including $Q$-values, reaction rates, and branching ratios. We find such mass variations can result in up to an order of magnitude local change in the final abundance pattern produced in an $r$-process simulation. We identify key nuclei whose masses have a substantial impact on abundance predictions for hot, cold, and neutron star merger $r$-process scenarios and could be measured at future radioactive beam facilities.

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.

We have performed for the first time a complete $r$-process mass sensitivity study in the $N=82$ region. We take into account how an uncertainty in a single nuclear mass propagates to influence important quantities of neighboring nuclei, including Q-values and reaction rates. We demonstrate that nuclear mass uncertainties of $\pm0.5$ MeV in the $N=82$ region result in up to an order of magnitude local change in $r$-process abundance predictions. We identify key nuclei in the study whose mass has a substantial impact on final $r$-process abundances and could be measured at future radioactive beam facilities.

The temperature in the crust of an accreting neutron star, which comprises its outermost kilometer, is set by heating from nuclear reactions at large densities, neutrino cooling, and heat transport from the interior. The heated crust has been thought to affect observable phenomena at shallower depths, such as thermonuclear bursts in the accreted envelope. Here we report that cycles of electron capture and its inverse, $\beta^-$ decay, involving neutron-rich nuclei at a typical depth of about 150 m, cool the outer neutron star crust by emitting neutrinos while also thermally decoupling the surface layers from the deeper crust. This Urca mechanism has been studied in the context of white dwarfs and Type Ia supernovae, but hitherto was not considered in neutron stars, because previous models computed the crust reactions using a zero-temperature approximation and assumed that only a single nuclear species was present at any given depth. This thermal decoupling means that X-ray bursts and other surface phenomena are largely independent of the strength of deep crustal heating. The unexpectedly short recurrence times, of the order of years, observed for very energetic thermonuclear superbursts are therefore not an indicator of a hot crust, but may point instead to an unknown local heating mechanism near the neutron star surface.

Short-lived radionuclides (SLRs) with lifetimes \tau < 100 Ma are known to have been extant when the Solar System formed over 4.5 billion years ago. Identifying the sources of SLRs is important for understanding the timescales of Solar System formation and processes that occurred early in its history. Extinct 36Cl (t_1/2 = 0.301 Ma) is thought to have been produced by interaction of solar energetic particles (SEPs), emitted by the young Sun, with gas and dust in the nascent Solar System. However, models that calculate SLR production in the early Solar System (ESS) lack experimental data for the 36Cl production reactions. We present here the first measurement of the cross section of one of the main 36Cl production reactions, 33S(\alpha,p)36Cl, in the energy range 0.70 - 2.42 MeV/A. The cross section measurement was performed by bombarding a target and collecting the recoiled 36Cl atoms produced in the reaction, chemically processing the samples, and measuring the 36Cl/Cl ratio of the activated samples with accelerator mass spectrometry (AMS). The experimental results were found to be systematically higher than the cross sections used in previous local irradiation models and other Hauser-Feshbach calculated predictions. However, the effects of the experimentally measured cross sections on the modeled production of 36Cl in the early Solar System were found to be minimal. Reactions channels involving S targets dominate 36Cl production, but the astrophysical event parameters can dramatically change each reactions' relative contribution.

Despite much effort in the past decades, the C-burning reaction rate is uncertain by several orders of magnitude, and the relative strength between the different channels 12C(12C,alpha)20Ne, 12C(12C,p)23Na and 12C(12C,n)23Mg is poorly determined. Additionally, in C-burning conditions a high 12C+12C rate may lead to lower central C-burning temperatures and to 13C(alpha,n)16O emerging as a more dominant neutron source than 22Ne(alpha,n)25Mg, increasing significantly the s-process production. This is due to the rapid decrease of the 13N(gamma,p)12C with decreasing temperature, causing the 13C production via 13N(beta+)13C. Presented here is the impact of the 12C+12C reaction uncertainties on the s-process and on explosive p-process nucleosynthesis in massive stars, including also fast rotating massive stars at low metallicity. Using various 12C+12C rates, in particular an upper and lower rate limit of ~ 50000 higher and ~ 20 lower than the standard rate at 5*10^8 K, five 25 Msun stellar models are calculated. The enhanced s-process signature due to 13C(alpha,n)16O activation is considered, taking into account the impact of the uncertainty of all three C-burning reaction branches. Consequently, we show that the p-process abundances have an average production factor increased up to about a factor of 8 compared to the standard case, efficiently producing the elusive Mo and Ru proton-rich isotopes. We also show that an s-process being driven by 13C(alpha,n)16O is a secondary process, even though the abundance of 13C does not depend on the initial metal content. Finally, implications for the Sr-peak elements inventory in the Solar System and at low metallicity are discussed.

We consider the astrophysical reaction rates for radiative neutron capture reactions ($n,\gamma$) in the crust of a neutron star. The presence of degenerate neutrons at high densities (mainly in the inner crust) can drastically affect the reaction rates. Standard rates assuming a Maxwell-Boltzmann distribution for neutrons can underestimate the rates by several orders of magnitude. We derive simple analytical expressions for reaction rates at a variety of conditions with account for neutron degeneracy. We also discuss the plasma effects on the outgoing radiative transition channel in neutron radiative capture reactions and show that these effects can also increase the reaction rates by a few orders of magnitude. In addition, using detailed balance, we analyze the effects of neutron degeneracy and plasma physics on reverse ($\gamma,n$) photodisintegration. We discuss the dependence of the reaction rates on temperature and neutron chemical potential and outline the efficiency of these reactions in the neutron star crust.

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.

The location of electron capture heat sources in the crust of accreting neutron stars depends on the masses of extremely neutron-rich nuclei. We present first results from a new implementation of the time-of-flight technique to measure nuclear masses of rare isotopes at the National Superconducting Cyclotron Laboratory. The masses of 16 neutron-rich nuclei in the scandium -- nickel range were determined simultaneously, improving the accuracy compared to previous data in 12 cases. The masses of $^{61}${V}, $^{63}${Cr}, $^{66}${Mn}, and $^{74}${Ni} were measured for the first time with mass excesses of $-30.510(890)$ MeV, $-35.280(650)$ MeV, $-36.900(790)$ MeV, and $-49.210(990)$ MeV, respectively. With the measurement of the $^{66}$Mn mass, the locations of the two dominant electron capture heat sources in the outer crust of accreting neutron stars that exhibit superbursts are now experimentally constrained. We find that the location of the $^{66}$Fe$\rightarrow^{66}$Mn electron capture transition occurs significantly closer to the surface than previously assumed because our new experimental Q-value is 2.1 MeV (2.6$\sigma$) smaller than predicted by the FRDM mass model.

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).

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.

We analyze thermonuclear and pycnonuclear fusion reactions in dense matter containing atomic nuclei of different types. We extend a phenomenological expression for the reaction rate, proposed recently by Gasques et al. (2005) for the one-component plasma of nuclei, to the multi-component plasma. The expression contains several fit parameters which we adjust to reproduce the best microscopic calculations available in the literature. Furthermore, we show that pycnonuclear burning is drastically affected by an (unknown) structure of the multi-component matter (a regular lattice, a uniform mix, etc.). We apply the results to study nuclear burning in a carbon_12-oxygen_16 mixture. In this context we present new calculations of the astrophysical S-factors for carbon-oxygen and oxygen-oxygen fusion reactions. We show that the presence of a CO lattice can strongly suppress carbon ignition in white dwarf cores and neutron star crusts at densities > 3e9 g cm^{-3} and temperatures T<1e8 K.

In this paper we analyze the nuclear fusion rate between equal nuclei for all five different nuclear burning regimes in dense matter (two thermonuclear regimes, two pycnonuclear ones, and the intermediate regime). The rate is determined by Coulomb barrier penetration in dense environments and by the astrophysical S-factor at low energies. We evaluate previous studies of the Coulomb barrier problem and propose a simple phenomenological formula for the reaction rate which covers all cases. The parameters of this formula can be varied, taking into account current theoretical uncertainties in the reaction rate. The results are illustrated for the example of the ^{12}C+^{12}C fusion reaction. This reaction is very important for the understanding of nuclear burning in evolved stars, in exploding white dwarfs producing type Ia supernovae, and in accreting neutron stars. The S-factor at stellar energies depends on a reliable fit and extrapolation of the experimental data. We calculate the energy dependence of the S-factor using a recently developed parameter-free model for the nuclear interaction, taking into account the effects of the Pauli nonlocality. For illustration, we analyze the efficiency of carbon burning in a wide range of densities and temperatures of stellar matter with the emphasis on carbon ignition at densities rho > 10^9 g/cc.