Dissipation in relativistic superfluid neutron stars

We analyze damping of oscillations of general relativistic superfluid neutron stars. To this aim we extend the method of decoupling of superfluid and normal oscillation modes first suggested in [Gusakov & Kantor PRD 83, 081304(R) (2011)]. All calculations are made self-consistently within the finite temperature superfluid hydrodynamics. The general analytic formulas are derived for damping times due to the shear and bulk viscosities. These formulas describe both normal and superfluid neutron stars and are valid for oscillation modes of arbitrary multipolarity. We show that: (i) use of the ordinary one-fluid hydrodynamics is a good approximation, for most of the stellar temperatures, if one is interested in calculation of the damping times of normal f-modes; (ii) for radial and p-modes such an approximation is poor; (iii) the temperature dependence of damping times undergoes a set of rapid changes associated with resonance coupling of neighboring oscillation modes. The latter effect can substantially accelerate viscous damping of normal modes in certain stages of neutron-star thermal evolution.

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.