For the first time nonradial oscillations of superfluid nonrotating stars are self-consistently studied at finite stellar temperatures. We apply a realistic equation of state and realistic density dependent model of critical temperature of neutron an
d proton superfluidity. In particular, we discuss three-layer configurations of a star with no neutron superfluidity at the centre and in the outer region of the core but with superfluid intermediate region. We show, that oscillation spectra contain a set of modes whose frequencies can be very sensitive to temperature variations. Fast temporal evolution of the pulsation spectrum in the course of neutron star cooling is also analysed.
The breaking stress (the maximum of the stress-strain curve) of neutron star crust is important for neutron star physics including pulsar glitches, emission of gravitational waves from static mountains, and flares from star quakes. We perform many mo
lecular dynamic simulations of the breaking stress at different coupling parameters (inverse temperatures) and strain rates. We describe our results with the Zhurkov model of strength. We apply this model to estimate the breaking stress for timescales ~1 s - 1 year, which are most important for applications, but much longer than can be directly simulated. At these timescales the breaking stress depends strongly on the temperature. For coupling parameter <200, matter breaks at very small stress, if it is applied for a few years. This viscoelastic creep can limit the lifetime of mountains on neutron stars. We also suggest an alternative model of timescale-independent breaking stress, which can be used to estimate an upper limit on the breaking stress.
Using the results of extensive Monte Carlo simulations we discuss corrections to the linear mixing rule in strongly coupled binary ionic mixtures. We analyze the plasma screening function at zero separation, H_{jk}(0), for two ions (of types j=1,2 an
d k=1,2) in a strongly coupled binary mixture. The function H_{jk}(0) is estimated by two methods: (1) from the difference of Helmholtz Coulomb free energies at large and zero separations; (2) by fitting the Widom expansion of H_{jk}(x) in powers of interionic distance x to Monte Carlo data on the radial pair distribution function g_{jk}(x). These methods are shown to be in good agreement. For illustration, we analyze the plasma screening enhancement of nuclear burning rates in dense stellar matter.
We analyze the effect of plasma screening on nuclear reaction rates in dense matter composed of atomic nuclei of one or two types. We perform semiclassical calculations of the Coulomb barrier penetrability taking into account a radial mean field pote
ntial of plasma ions. The mean field potential is extracted from the results of extensive Monte Carlo calculations of radial pair distribution functions of ions in binary ionic mixtures. We calculate the reaction rates in a wide range of plasma parameters and approximate these rates by an analytical expression that is expected to be applicable for multicomponent ions mixtures. Also, we analyze Gamow-peak energies of reacting ions in various nuclear burning regimes. For illustration, we study nuclear burning in C-O mixtures.
We analyze the thermal conductivity of ions (equivalent to the conductivity of phonons in crystalline matter) in a neutron star envelope. We calculate the ion/phonon thermal conductivity in a crystal of atomic nuclei using variational formalism and
performing momentum-space integration by Monte Carlo method. We take into account phonon-phonon and phonon-electron scattering mechanisms and show that phonon-electron scattering dominates at not too low densities. We extract the ion thermal conductivity in ion liquid or gas from literature. Numerical values of the ion/phonon conductivity are approximated by analytical expressions, valid for T>10^5 K and 10^5 g cm^-3 < rho < 10^14 g cm^-3. Typical magnetic fields B~10^12 G in neutron star envelopes do not affect this conductivity although they strongly reduce the electron thermal conductivity across the magnetic field. The ion thermal conductivity remains much smaller than the electron conductivity along the magnetic field. However, in the outer neutron star envelope it can be larger than the electron conductivity across the field, that is important for heat transport across magnetic field lines in cooling neutron stars. The ion conductivity can greatly reduce the anisotropy of heat conduction in outer envelopes of magnetized neutron stars.
We compare Path Integral Monte Carlo calculations by Militzer and Pollock (Phys. Rev. B 71, 134303, 2005) of Coulomb tunneling in nuclear reactions in dense matter to semiclassical calculations assuming WKB Coulomb barrier penetration through the rad
ial mean-field potential. We find a very good agreement of two approaches at temperatures higher than ~1/5 of the ion plasma temperature. We obtain a simple parameterization of the mean field potential and of the respective reaction rates. We analyze Gamow-peak energies of reacting ions in various reaction regimes and discuss theoretical uncertainties of nuclear reaction rates taking carbon burning in dense stellar matter as an example.
