We study long-term thermal evolution of neutron stars in soft X-ray transients (SXTs), taking the deep crustal heating into account consistently with the changes of the composition of the crust. We collect observational estimates of average accretion rates and thermal luminosities of such neutron stars and compare the theory with observations. We perform simulations of thermal evolution of accreting neutron stars, considering the gradual replacement of the original nonaccreted crust by the reprocessed accreted matter, the neutrino and photon energy losses, and the deep crustal heating due to nuclear reactions in the accreted crust. We test and compare results for different modern theoretical models. We update a compilation of the observational estimates of the thermal luminosities in quiescence and average accretion rates in the SXTs and compare the observational estimates with the theoretical results. Long-term thermal evolution of transiently accreting neutron stars is nonmonotonic. The quasi-equilibrium temperature in quiescence reaches a minimum and then increases toward the final steady state. The quasi-equilibrium thermal luminosity of a neutron star in an SXT can be substantially lower at the minimum than in the final state. This enlarges the range of possibilities for theoretical interpretation of observations of such neutron stars. The updates of the theory and observations leave unchanged the previous conclusions that the direct Urca process operates in relatively cold neutron stars and that an accreted heat-blanketing envelope is likely present in relatively hot neutron stars in the SXTs in quiescence. The results of the comparison of theory with observations favor suppression of the triplet pairing type of nucleon superfluidity in the neutron-star matter.

We study equation of state (EOS) of an accreting neutron star crust. Usually, such EOS is obtained assuming (implicitly) that the free (unbound) neutrons and nuclei in the inner crust move together. We argue, that this assumption violates the condition $\mu_n^\infty={\rm const}$, required for hydrostatic (and diffusion) equilibrium of unbound neutrons ($\mu^\infty_n$ is the redshifted neutron chemical potential). We construct a new EOS respecting this condition, working in the compressible liquid-drop approximation. We demonstrate that it is close to the catalyzed EOS in most part of the inner crust, being very different from EOSs of accreted crust discussed in the literature. In particular, the pressure at the outer-inner crust interface does not coincide with the neutron drip pressure, usually calculated in the literature, and is determined by hydrostatic (and diffusion) equilibrium conditions within the star. We also find an instability at the bottom of fully accreted crust that transforms nuclei into homogeneous nuclear matter. It guarantees that the structure of fully accreted crust remains self-similar during accretion.

Observed temperatures of transiently accreting neutron stars in the quiescent state are generally believed to be supported by deep crustal heating, associated with non-equilibrium exothermic reactions in the crust. Traditionally, these reactions are studied by considering nuclear evolution governed by compression of the accreted matter. Here we show that this approach has a basic weakness, that is that in some regions of the inner crust the conservative forces, applied for matter components (nuclei and neutrons), are not in mechanical equilibrium. In principle the force balance can be restored by dissipative forces, however the required diffusion fluxes are of the same order as total baryon flux at Eddington accretion. We argue that redistribution of neutrons in the inner crust should be involved in realistic model of accreted crust.

Correct interpretation of X-ray observations of transiently accreting neutron stars requires modeling of nuclear-physical processes in these objects. We consider a chain of nuclear reactions that drives the crust composition in an accreting neutron star and heats up the star. We constructed multicomponent approach with the kinetics of nuclear reactions described in simplified stepwise manner. The redistribution of nucleons between nuclei by emission and capture of neutrons is shown to significantly affect the nuclear reaction chains and the composition of the inner crust. In particular, even if the outer crust has one-component composition, the appearance of free neutrons in the inner crust leads to branching of reaction chains and formation of the multicomponent composition. We apply the compressible liquid drop nuclear model, which includes effects of free neutrons on nuclear energies. It allows us to calculate the composition, the heating profile and the equation of state of matter up to densities \rho ~ 2 \times 10^{13} g cm^{-3}.

Transiently accreting neutron stars in low mass X-ray binaries are generally believed to be heated up by nuclear reactions in accreted matter during hydrostatic compression. Detailed modeling of these reactions is required for the correct interpretation of observations. In this paper, we construct a simplified reaction network, which can be easily implemented and depends mainly on atomic mass tables as nuclear physics input. We show that it reproduces results of the detailed network by Lau et al. (2018) very well, if one applies the same mass model. However, the composition and the heating power are shown to be sensitive to the mass table used and treatment of mass tables boundary, if one applies several of them in one simulation. In particular, the impurity parameter $Q_\mathrm{imp}$ at density $\rho=2\times 10^{12}$ g cm$^{-3}$ can differ for a factor of few, and even increase with density increase. The profile of integrated heat realize shown to be well confined between results by Fantina et al. (2018) and Lau et al. (2018).

With a spin frequency of 707 Hz, PSR J0952-0607 is the second fastest spinning pulsar known. It was discovered in radio by LOFAR in 2017 at an estimated distance of either 0.97 or 1.74 kpc and has a low-mass companion with a 6.42 hr orbital period. We report discovery of the X-ray counterpart of PSR J0952-0607 using XMM-Newton. The X-ray spectra can be well-fit by a single power law model (Gamma = 2.5) or by a thermal plus power law model (kTeff = 40 eV and Gamma = 1.4). We do not detect evidence of variability, such as that due to orbital modulation from pulsar wind and companion star interaction. Because of its fast spin rate, PSR J0952-0607 is a crucial source for understanding the r-mode instability, which can be an effective mechanism for producing gravitational waves. Using the high end of our measured surface temperature, we infer a neutron star core temperature of ~10^7 K, which places PSR J0952-0607 within the window for the r-mode to be unstable unless an effect such as superfluid mutual friction damps the fluid oscillation. The measured luminosity limits the dimensionless r-mode amplitude to be less than ~1x10^-9.

Dewitt et al. (1973) derived a useful relation between the plasma screening factor for a reaction of two fusing ions and their chemical potentials, based on the plasma pair distribution functions. We show that their result can be derived in a simpler, more straightforward way, by applying the principle of detailed balance, which also enables us to generalize the relation to reactions involving $N$ fusing ions. In order to demonstrate the usefulness of applying the principle of detailed balance, we calculate the screening factor for the $pep$ reaction, $p+e+p\rightarrow ^{2}$D $+\nu_{e}$. For the plasma conditions near the centre of the Sun, the reaction is suppressed by roughly the same amount by which the reaction $p+p\rightarrow ^{2}$D $+$ $e^{+}+\nu_{e}$ is enhanced. This effect may be measured in the near future.

I suggest a novel type of nuclear reactions in accreting neutron stars - neutron transfer, which is quantum tunneling of weakly bounded neutron from one nucleus to another. I estimate the rate of this process for fixed nuclei separation and then average the result over realistic distribution of nuclei to get the rate value for astrophysical conditions. The neutron transfer can modify reaction chains in accreting neutron stars, thus affecting their heating and cooling. In particular, it can suppress cooling by URCA pairs of nuclei, which is supposed to be crucial for the hottest neutron stars.

I consider differential rotation, associated with radiation-driven Chandrasekhar-Friedman-Schutz (CFS) instability and its possible observational evidences. I focus on the evolution of the apparent spin frequency, which is typically associated with the motion of a specific point on the stellar surface (e.g., polar cap). I start from long-term evolution (on the timescale when instability significantly changes the spin frequency). For this case, I reduce evolution equations to one differential equation and demonstrate that it can be directly derived from energy conservation law. This equation governs evolution rate through sequence of thermally equilibrium states and provides linear coupling for the cooling power and rotation energy losses via gravitational wave emission. In particular, it shows that differential rotation do not affect long-term spin-down. On the contrary, on the short timescales, differential rotation can significantly modify the apparent spin-down, if one examines a strongly unstable star with very small initial amplitude of unstable mode. This statement is confirmed by consideration of Newtonian non-magnetized perfect fluid and dissipative stellar models as well as magnetized stellar model. For example, despite that widely applied evolution equations predict effective spin to be constant in absence of dissipation, the CFS-unstable star should be observed as spinning-down. However, effects of differential rotation on apparent spin-down are negligible for realistic models of neutron star recycling, unless: (1) the neutron star is not magnetized, and (2) r-mode amplitude is modulated faster than the shear viscosity dissipation timescale, being large enough that spin-down can be measured on modulation timescale.

The strength of neutron star crust is crucial for modelling magnetar flares, pulsar glitches and gravitational wave emission. We aim to shed some light on this problem by analysing uniaxial stretch deformation (elongation and contraction) of perfect body-centered cubic Coulomb crystals, paying special attention to the inherent anisotropy of this process. Our analysis is based on the semi-analytical approach of Baiko and Kozhberov (2017), which, for any uniform deformation, allows one to calculate, in fully non-linear regime, critical deformation parameters beyond which the lattice loses its dynamic stability. We determine critical strain, pressure anisotropy and deformation energy for any stretch direction with respect to the crystallographic axes. These quantities are shown to be strongly anisotropic: they vary by a factor of almost 10 depending on the orientation of the deformation axis. For polycrystalline crust, we argue that the maximum strain for the stretch deformation sustainable elastically is 0.04. It is lower than the breaking strain of 0.1 obtained in molecular dynamic simulations of a shear deformation by Horowitz and Kadau (2009). The maximum pressure anisotropy of polycrystalline matter is estimated to be in the range from 0.005 to 0.014 $nZ^2e^2/a$, where $n$ is the ion number density, $Ze$ is the ion charge, and $a$ is the ion-sphere radius. We discuss possible mechanisms of plastic motion and formation of large crystallites in neutron star crust as well as analyse energy release associated with breaking of such crystallites in the context of magnetic field evolution and magnetar flaring activity.

At second order in perturbation theory, the unstable r-mode of a rotating star includes growing differential rotation whose form and growth rate are determined by gravitational-radiation reaction. With no magnetic field, the angular velocity of a fluid element grows exponentially until the mode reaches its nonlinear saturation amplitude and remains nonzero after saturation. With a background magnetic field, the differential rotation winds up and amplifies the field, and previous work where large mode amplitudes were considered suggests that the amplification may damp out the instability. A background magnetic field, however, turns the saturated time-independent perturbations corresponding to adding differential rotation into perturbations whose characteristic frequencies are of order the Alfv\'en frequency. As found in previous studies, we argue that magnetic- field growth is sharply limited by the saturation amplitude of an unstable mode. In contrast to previous work, however, we show that if the amplitude is small, i.e., of order 10^(-4), then the limit on the magnetic-field growth is stringent enough to prevent the loss of energy to the magnetic field from damping or significantly altering an unstable r-mode in nascent neutron stars with normal interiors and in cold stars whose interiors are type II superconductors. We show this result first for a toy model, and we then obtain an analogous upper limit on magnetic field growth using a more realistic model of a rotating neutron star. Our analysis depends on the assumption that there are no marginally unstable perturbations, and this may not hold when differential rotation leads to a magnetorotational instability.

We discuss magnetic field enhancement by unstable r-modes (driven by the gravitational radiation reaction force) in rotating stars. In the absence of a magnetic field, gravitational radiation exponentially increases the r-mode amplitude $\alpha$, and accelerates differential rotation (secular motion of fluid elements). For a magnetized star, differential rotation enhances the magnetic field energy. Rezzolla et al. (2000--2001) argued that if the magnetic energy grows faster than the gravitational radiation reaction force pumps energy into the r-modes, then the r-mode instability is suppressed. Chugunov (2015) demonstrated that without gravitational radiation, differential rotation can be treated as a degree of freedom decoupled from the r-modes and controlled by the back reaction of the magnetic field. In particular, the magnetic field windup does not damp r-modes. Here we discuss the effect of the back reaction of the magnetic field on differential rotation of unstable r-modes, and show that it limits the generated magnetic field and the magnetic energy growth rate preventing suppression of the r-mode instability by magnetic windup at low saturation amplitudes, $\alpha \ll 1$, predicted by current models.

I suggest a novel approach for deriving evolution equations for rapidly rotating relativistic stars affected by radiation-driven Chandrasekhar-Friedman-Schutz (CFS) instability. This approach is based on the multipolar expansion of gravitational wave emission and appeals to the global physical properties of the star (energy, angular momentum, and thermal state), but not to canonical energy and angular momentum, which is traditional. It leads to simple derivation of the CFS instability criterion for normal modes and the evolution equations for a star, affected by this instability. The approach also gives a precise form to simple explanation of the CFS instability: it occurs when two conditions met: (a) gravitational wave emission removes angular momentum from the rotating star (thus releasing the rotation energy) and (b) gravitational waves carry less energy, than the released amount of the rotation energy. Results are illustrated by r-mode instability in slowly rotating Newtonian stars. In resulting evolution equations emission of gravitational waves directly affects the spin frequency, being in line with the arguments by Levin & Ushomirsky (2001), but in contrast to widely accepted equations by Lindblom et al. (1998); Ho & Lai (2000), there effective spin frequency decrease is coupled with dissipation of unstable mode. This problem, initially stressed by Levin & Ushomirsky (2001), is shown to be superficial, and arises as a result of definition of the effective spin frequency by Lindblom et al. (1998); Ho & Lai (2000). Namely, it is shown, that if this definition is taken into account properly, the evolution equations coincide with obtained here in the leading order in mode amplitude. I also argue that the next-to-leading order terms in evolution equations (which differ for Owen et al. 1998 and Ho & Lai 2000) require clarification and thus it would be more self-consistent to omit them.

We combined Chandra ACIS observations of the globular cluster 47 Tucanae (hereafter, 47 Tuc) from 2000, 2002, and 2014-15 to create a deeper X-ray source list, and study some of the faint radio millisecond pulsars (MSPs) present in this cluster. We have detected 370 X-ray sources within the half-mass radius (2$'$.79) of the cluster, 81 of which are newly identified, by including new data and using improved source detection techniques. The majority of the newly identified sources are in the crowded core region, indicating cluster membership. We associate five of the new X-ray sources with chromospherically active BY Dra or W UMa variables identified by Albrow et al. (2001). We present alternative positions derived from two methods, centroiding and image reconstruction, for faint, crowded sources. We are able to extract X-ray spectra of the recently discovered MSPs 47 Tuc aa, 47 Tuc ab, the newly timed MSP 47 Tuc Z, and the newly resolved MSPs 47 Tuc S and 47 Tuc F. Generally, they are well fit by black body or neutron star atmosphere models, with temperatures, luminosities and emitting radii similar to those of other known MSPs in 47 Tuc, though 47 Tuc aa and 47 Tuc ab reach lower X-ray luminosities. We limit X-ray emission from the full surface of the rapidly spinning (542 Hz) MSP 47 Tuc aa, and use this limit to put an upper bound for amplitude of r-mode oscillations in this pulsar as $\alpha<2.5\times 10^{-9}$ and constrain the shape of the r-mode instability window.

To put new constraints on the r-mode instability window, we analyse the formation of millisecond pulsars (MSPs) within the recycling scenario, making use of three sets of observations: (a) X-ray observations of neutron stars (NSs) in low-mass X-ray binaries; (b) timing of millisecond pulsars; and (c) X-ray and UV observations of MSPs. As shown in previous works, r-mode dissipation by shear viscosity is not sufficient to explain observational set (a), and enhanced r-mode dissipation at the red-shifted internal temperatures $T^\infty\sim 10^8$ K is required to stabilize the observed NSs. Here, we argue that models with enhanced bulk viscosity can hardly lead to a self-consistent explanation of observational set (a) due to strong neutrino emission, which is typical for these models (unrealistically powerful energy source is required to keep NSs at the observed temperatures). We also demonstrate that the observational set (b), combined with the theory of internal heating and NS cooling, provides evidence of enhanced r-mode dissipation at low temperatures, $T^\infty\sim 2\times 10^7$ K. Observational set (c) allows us to set an upper limit on the internal temperatures of MSPs, $T^\infty<2\times 10^7$ K (assuming a canonical NS with the accreted crust). Recycling scenario can produce MSPs at these temperatures only if r-mode instability is suppressed in the whole MSP spin frequency range ($\nu\lesssim 750$ Hz) at temperatures $2\times 10^7\lesssim T^\infty\lesssim 3 \times 10^7$ K, providing thus a new constraint on the r-mode instability window. These observational constraints are analysed in more details in application to the resonance uplift scenario of Gusakov et al. [Phys. Rev. Lett., 112 (2014), 151101].

In the recent papers by Gusakov, Chugunov, and Kantor (2014) a new scenario describing evolution of rapidly rotating neutron stars in low-mass X-ray binaries was proposed. The scenario accounts for a resonant interaction of normal r modes with superfluid inertial modes at some specific internal stellar temperatures ("resonance temperatures"). This interaction results in an enhanced damping of r mode and appearance of the "stability peaks" in the temperature -- spin frequency plane, which split the r-mode instability window in the vicinity of the resonance temperatures. The scenario suggests that the hot and rapidly rotating NSs spend most of their life climbing up these peaks and, in particular, are observed there at the moment. We analyze in detail possible observational signatures of this suggestion. In particular, we show that these objects may exhibit `anti-glitches' -- sudden frequency jumps on a time scale of hours-months.

Rezzolla et al. [ApJ 531 (2000), L139; Phys. Rev. D 64 (2001), 104013; Phys. Rev. D 64 (2001), 104014] draw attention to the second order secular drift associated with r-modes and claimed that it should lead to magnetic field enhancement and suppression of r-mode instability in magnetized neutron stars. We critically revise these results. We present a particular second order r-mode solution with vanishing secular drift, thus refuting a widely believed statement that secular drift is an unavoidable feature of r-modes. This non-drifting solution is not affected by magnetic field $B$, if $B\ll B_{\mathrm{crit}}\approx 10^{17}\,(\nu/600\,\mathrm{Hz})$ G ($\nu$ is a spin frequency) and does not lead to secular evolution of magnetic field. For general second order r-mode solution the drift does not necessarily vanish, but the solution can be presented as a superposition of two solutions: one describes evolution of differential rotation in nonoscillating star (which describes secular drift; for nonmagnetized star it is arbitrary stationary rotation stratified on cylinders; for magnetized star differential rotation evolves on the Alfv\'{e}n timescale and may lead to magnetic energy enhancement), and another one is non-drifting r-mode solution mentioned above. This representation allows us to conclude that enhancement of magnetic field energy is limited by initial energy of differential rotation, which is much less (for a factor $\propto \alpha^2$, where $\alpha$ is mode amplitude) than the total energy of r-mode. Hence, magnetic field enhancement by drift cannot suppress r-mode instability. Results can be generalized for any oscillation mode in any medium, if this mode has non-drifting solution for $B=0$.

A new class of neutron stars (NSs) -- hot rapidly rotating non-accreting NSs, which we propose to call HOFNARs (HOt and Fast Non-Accreting Rotators) or "hot widows" (in analogy with "black widow" pulsars) -- is suggested. We argue that such stars should originate from the low-mass X-ray binaries (LMXBs) provided that they were unstable with respect to excitation of $r$-modes at the end of accretion epoch (when their low-mass companions ceased to fill the Roche lobe). High temperature of "hot widows"/HOFNARs is maintained by $r$-mode dissipation rather than by accretion. We analyse observational properties of "hot widows"/HOFNARs and demonstrate that these objects form a specific separate class of neutron stars. In particular, some of the most stable X-ray sources among the candidates to quiescent LMXB systems (qLMXBs), can, in fact, belong to that new class. We formulate observational criteria which allow to distinguish "hot widows"/HOFNARs from qLMXB systems, and argue that available observations of X-ray sources 47 Tuc X5 and X7 satisfy (or, at least, do not contradict) these criteria. In addition, we discuss indirect evidences in favor of "hot widows"/HOFNARs existence, following from the analysis of observations and predictions of population synthesis theories. If that new class of NSs does exist, it would prove the possibility to emit gravitational waves by mass-current multipole. Various applications of our results, such as prospects for constraining superdense matter properties with hot widows"/HOFNARs, are analyzed.

In a previous paper [M. E. Gusakov, A. I. Chugunov, and E. M. Kantor, Phys. Rev. Lett. 112, 151101 (2014)], we introduced a new scenario that explains the existence of rapidly rotating warm neutron stars (NSs) observed in low-mass X-ray binaries (LMXBs). Here it is described in more detail. The scenario takes into account the interaction between superfluid inertial modes and the normal (quadrupole) $m=2$ $r$-mode, which can be driven unstable by Chandrasekhar-Friedman-Schutz (CFS) mechanism. This interaction can only occur at some fixed "resonance" stellar temperatures; it leads to formation of the "stability peaks" which stabilize a star in the vicinity of these temperatures. We demonstrate that a NS in LMXB spends a substantial fraction of time on the stability peak, that is, in the region of stellar temperatures and spin frequencies, that has been previously thought to be CFS unstable with respect to excitation of $r$-modes. We also find that the spin frequencies of NSs are limited by the CFS instability of normal (octupole) $m=3$ $r$-mode rather than by $m=2$ $r$-mode. This result agrees with the predicted value of the cutoff spin frequency $\sim 730$ Hz in the spin distribution of accreting millisecond X-ray pulsars. In addition, we analyze evolution of a NS after the end of the accretion phase and demonstrate that millisecond pulsars can be born in LMXBs within our scenario. Besides millisecond pulsars, our scenario also predicts a new class of LMXB descendants - hot and rapidly rotating nonaccreting NSs ("hot widows"/HOFNARs). Further comparison of the proposed theory with observations of rotating NSs can impose new important constraints on the properties of superdense matter.

We study non-radial oscillations of neutron stars with superfluid baryons, in a general relativistic framework, including finite temperature effects. Using a perturbative approach, we derive the equations describing stellar oscillations, which we solve by numerical integration, employing different models of nucleon superfluidity, and determining frequencies and gravitational damping times of the quasi-normal modes. As expected by previous results, we find two classes of modes, associated to superfluid and non-superfluid degrees of freedom, respectively. We study the temperature dependence of the modes, finding that at specific values of the temperature, the frequencies of the two classes of quasi-normal modes show avoided crossings, and their damping times become comparable. We also show that, when the temperature is not close to the avoided crossings, the frequencies of the modes can be accurately computed by neglecting the coupling between normal and superfluid degrees of freedom. Our results have potential implications on the gravitational wave emission from neutron stars.

We consider an instability of rapidly rotating neutron stars in low-mass X-ray binaries (LMXBs) with respect to excitation of r-modes (which are analogous to Earth's Rossby waves controlled by the Coriolis force). We argue that finite temperature effects in the superfluid core of a neutron star lead to a resonance coupling and enhanced damping (and hence stability) of oscillation modes at certain stellar temperatures. Using a simple phenomenological model we demonstrate that neutron stars with high spin frequency may spend a substantial amount of time at these `resonance' temperatures. This finding allows us to explain puzzling observations of hot rapidly rotating neutron stars in LMXBs and to predict a new class of hot, non-accreting, rapidly rotating neutron stars, some of which may have already been observed and tentatively identified as quiescent LMXB (qLMXB) candidates. We also impose a new theoretical limit on the neutron star spin frequency, explaining the cut-off spin frequency ~730 Hz, following from the statistical analysis of accreting millisecond X-ray pulsars. Besides explaining the observations, our model provides a new tool to constrain superdense matter properties comparing measured and theoretically predicted resonance temperatures.

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.

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.

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 and 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 molecular 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 and 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.

Recently developed analytic approximation for the equation of state of fully ionized nonideal electron-ion plasma mixtures [Potekhin et al., Phys. Rev. E, 79, 016411 (2009); arXiv:0812.4344], which covers the transition between the weak and strong Coulomb coupling regimes and reproduces numerical results obtained in the hypernetted chain (HNC) approximation, is modified in order to fit the small deviations from the linear mixing in the strong coupling regime, revealed by recent Monte Carlo simulations. In addition, a mixing rule is proposed for the regime of weak coupling, which generalizes post-Debye density corrections to the case of mixtures and numerically agrees with the HNC approximation in that regime.

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 potential 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 radial 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.

We study acoustic oscillations (eigenfrequencies, velocity distributions, damping times) of normal crusts of strange stars. These oscillations are very specific because of huge density jump at the interface between the normal crust and the strange matter core. The oscillation problem is shown to be self-similar. For a low (but non-zero) multipolarity l the fundamental mode (without radial nodes) has a frequency ~300 Hz and mostly horizontal oscillation velocity; other pressure modes have frequencies >=20 kHz and almost radial oscillation velocities. The latter modes are similar to radial oscillations (have approximately the same frequencies and radial velocity profiles). The oscillation spectrum of strange stars with crust differs from the spectrum of neutron stars. If detected, acoustic oscillations would allow one to discriminate between strange stars with crust and neutron stars and constrain the mass and radius of the star.

We study eigenmodes of acoustic oscillations of high multipolarity l ~ 100 - 1000 and high frequency (~100 kHz), localized in neutron star envelopes. We show that the oscillation problem is self-similar. Once the oscillation spectrum is calculated for a given equation of state (EOS) in the envelope and given stellar mass M and radius R, it can be rescaled to a star with any M and R (but the same EOS in the envelope). For l>300 the modes can be subdivided into the outer and inner ones. The outer modes are mainly localized in the outer envelope. The inner modes are mostly localized near the neutron drip point, being associated with the softening of the EOS after the neutron drip. We calculate oscillation spectra for the EOSs of cold-catalyzed and accreted matter and show that the spectra of the inner modes are essentially different. A detection and identification of high-frequency pressure modes would allow one to infer M and R and determine also the EOS in the envelope (accreted or ground-state) providing thus a new and simple method to explore the main parameters and internal structure of neutron stars.

We calculate the electron shear viscosity (determined by Coulomb electron collisions) for a dense matter in a wide range of parameters typical for white dwarf cores and neutron star crusts. In the density range from ~10^3 g cm^-3 to 10^7-10^10 g cm^-3 we consider the matter composed of widely abundant astrophysical elements, from H to Fe. For higher densities, 10^10-10^14 g cm^-3, we employ the ground-state nuclear composition, taking into account finite sizes of atomic nuclei and the distribution of proton charge over the nucleus. Numerical values of the viscosity are approximated by an analytic expression convenient for applications. Using the approximation of plane-parallel layer we study eigenfrequencies, eigenmodes and viscous damping times of oscillations of high multipolarity, l~500-1000, localized in the outer crust of a neutron star. For instance, at l~500 oscillations have frequencies f >= 40 kHz and are localized not deeper than ~300 m from the surface. When the crust temperature decreases from 10^9 K to 10^7 K, the dissipation time of these oscillations (with a few radial nodes) decreases from ~1 year to 10-15 days.