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.

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.

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.