Ages and thermal luminosities of neutron stars, inferred from observations, can be interpreted with the aid of the neutron star cooling theory to gain information on the properties of superdense matter in neutron-star interiors. We present a survey of estimated ages, surface temperatures and thermal luminosities of middle-aged neutron stars with relatively weak or moderately strong magnetic fields, which can be useful for these purposes. The catalogue includes results selected from the literature, supplemented with new results of spectral analysis of a few cooling neutron stars. The data are compared with the theory. We show that overall agreement of theoretical cooling curves with observations improves substantially for models where neutron superfluidity in stellar core is weak.

We present an accurate analytic approximation for the energy of a quantum one-component Coulomb liquid of ions in a uniform electron background which has been recently calculated from first principles (Baiko 2019). The approximation enables us to develop in an analytic form a complete thermodynamic description of quantum ions in a practically important range of mass densities at temperatures above crystallization. We show that ionic quantum effects in liquid cores of white dwarfs (WDs) affect heat capacity, cooling, thermal compressibility, pulsation frequencies and radii of sufficiently cold WDs, especially with relatively massive helium and carbon cores.

We study the bulk viscosity in two (anti-spaghetti and Swiss cheese) phases of non- spherical nuclei in the mantle of an oscillating neutron star near the boundary with the stellar core. The bulk viscosity is produced by non-equilibrium Urca neutrino emis- sion processes. In the mantle, the direct Urca process may be open (Gusakov et al., 2004) if neutrons and protons move in a periodic potential created by a lattice of non- spherical nuclei (which allows the nucleons to have large quasi-momenta and satisfy direct Urca momentum-conservation). This bulk viscosity can dominate over that due to the modified Urca process in the outer stellar core and over the shear viscosity. The bulk viscosity depends strongly on temperature, oscillation frequency and nucleon superfluidity. The enhanced bulk viscosity in the mantle can control propagation and damping of neutron star oscillations.

We study thermal relaxation in a neutron star after internal heating events (outbursts) in the crust. We consider thin and thick spherically symmetric heaters, superfluid and non-superfluid crusts, stars with open and forbidden direct Urca processes in their cores. In particular, we analyze long-term thermal relaxation after deep crustal heating produced by nuclear transformations in fully or partly accreted crusts of transiently accreting neutron stars. This long-term relaxation has a typical relaxation time and an overall finite duration time for the crust to thermally equilibrate with the core. Neutron superfluidity in the inner crust greatly affects the relaxation if the heater is located in the inner crust. It shortens and unifies the time of emergence of thermal wave from the heater to the surface. This is important for the interpretation of observed outbursts of magnetars and transiently accreting neutron stars in quasi-persistent low-mass X-ray binaries.

We derive analytic approximations for the neutrino luminosities and the heat capacities of neutron stars with isothernal nucleon cores as functions of the mass and radius of stars. The neutrino luminosities are approximated for the three basic neutrino emission mechanisms, and the heat capacities for the five basic combinations of the partial heat capacities. The approximations are valid for for a wide class of equations of state of dense nucleonmatter. The results significantly simplify the theoretical interpretation of observations of cooling neutron stars as well as of quasistationary thermal states of neutron stars in X-ray transients. For illustration, we present an analysis of the neutrino cooling functions of nine isolated neutron stars taking into account the effects of their magnetic fields and of the presence of light elements in their heat blanketing envelopes. These results allow one to investigate the superfluid properties of neutron star cores.

We derive analytic approximations of neutrino luminosities and heat capacities of neutron stars with nucleon cores valid for a wide class of equations of state of dense nucleon matter. The neutrino luminosities are approximated for the three cases in which they are produced by (i) direct Urca or (ii) modified Urca processes in non-superfluid matter, or (iii) neutrino-pair bremsstrahlung in neutron-neutron collisions (when other neutrino reactions are suppressed by strong proton superfluidity). The heat capacity is approximated for the two cases of (i) non-superfluid cores and (ii) the cores with strong proton superfluidity. The results can greatly simplify numerical simulations of cooling neutron stars with isothermal interiors at the neutrino and photon cooling stages as well as simulations of quasi-stationary internal thermal states of neutron stars in X-ray transients. For illustration, a model-independent analysis of thermal states of the latter sources is outlined.

We study the effects of heat blanketing envelopes of neutron stars on their cooling. To this aim, we perform cooling simulations using newly constructed models of the envelopes composed of binary ion mixtures (H--He, He--C, C--Fe) varying the mass of lighter ions (H, He or C) in the envelope. The results are compared with those calculated using the standard models of the envelopes which contain the layers of lighter (accreted) elements (H, He and C) on top of the Fe layer, varying the mass of accreted elements. The main effect is that the chemical composition of the envelopes influences their thermal conductivity and, hence, thermal insulation of the star. For illustration, we apply these results to estimate the internal temperature of the Vela pulsar and to study the cooling of neutron stars of ages of 0.1 - 1 Myr at the photon cooling stage. The uncertainties of the cooling models associated with our poor knowledge of chemical composition of the heat insulating envelopes strongly complicate theoretical reconstruction of the internal structure of cooling neutron stars from observations of their thermal surface emission.

Investigations of thermal evolution of neutron stars with hyperon cores require neutrino emissivities for many neutrino reactions involving strongly degenerate particles (nucleons, hyperons, electrons, muons). We calculate the angular integrals $I_n$ (over orientations of momenta of $n$ degenerate particles) for major neutrino reactions with $n$ =3, 4, 5 at all possible combinations of particle Fermi momenta. The integrals $I_n$ are necessary ingredients for constructing a uniform database of neutrino emissivities in dense nucleon-hyperon matter. The results can also be used in many problems of physical kinetics of strongly degenerate systems.

We construct new models of outer heat blanketing envelopes of neutron stars composed of binary ion mixtures (H - He, He - C, C - Fe) in and out of diffusive equilibrium. To this aim, we generalize our previous work on diffusion of ions in isothermal gaseous or Coulomb liquid plasmas to handle non-isothermal systems. We calculate the relations between the effective surface temperature Ts and the temperature Tb at the bottom of heat blanketing envelopes (at a density rhob= 1e8 -- 1e10 g/cc) for diffusively equilibrated and non-equilibrated distributions of ion species at different masses DeltaM of lighter ions in the envelope. Our principal result is that the Ts - Tb relations are fairly insensitive to detailed distribution of ion fractions over the envelope (diffusively equilibrated or not) and depend almost solely on DeltaM. The obtained relations are approximated by analytic expressions which are convenient for modeling the evolution of neutron stars.

The electron shear viscosity due to Coulomb scattering of degenerate electrons by atomic nuclei throughout a magnetized neutron star crust is calculated. The theory is based on the shear viscosity coefficient calculated neglecting magnetic fields but taking into account gaseous, liquid and solid states of atomic nuclei, multiphonon scattering processes, and finite sizes of the nuclei albeit neglecting the effects of electron band structure. The effects of strong magnetic fields are included in the relaxation time approximation with the effective electron relaxation time taken from the field-free theory. The viscosity in a magnetized matter is described by five shear viscosity coefficients. They are calculated and their dependence on the magnetic field and other parameters of dense matter is analyzed. Possible applications and open problems are outlined.

The central compact object in the supernova remnant HESS J1731-347 appears to be the hottest observed isolated cooling neutron star. The cooling theory of neutron stars enables one to explain observations of this star by assuming the presence of strong proton superfluidity in the stellar core and the existence of the surface heat blanketing envelope which almost fully consists of carbon. The cooling model of this star is elaborated to take proper account of the neutrino emission due to neutron-neutron collisions which is not suppressed by proton superfluidity. Using the results of spectral fits of observed thermal spectra for the distance of 3.2 kpc and the cooling theory for the neutron star of age 27 kyr, new constraints on the stellar mass and radius are obtained which are more stringent than those derived from the spectral fits alone.

We present a simple model for electron transport in a possible layer of exotic nuclear clusters (in the so called nuclear pasta layer) between the crust and liquid core of a strongly magnetized neutron star. The electron transport there can be strongly anisotropic and gyrotropic. The anisotropy is produced by different electron effective collision frequencies along and across local symmetry axis in domains of exotic ordered nuclear clusters and by complicated effects of the magnetic field. We also calculate averaged kinetic coefficients in case local domains are freely oriented. Possible applications of the obtained results and open problems are outlined.

We apply our recently suggested statistical approach to thermal evolution of isolated neutron stars and accreting quasistationary neutron stars in X-ray transients for constraining the position and relative broadening alpha of the direct Urca threshold of powerful neutrino emission in neutron star cores. We show that most likely explanation of observations corresponds to alpha = 0.08 - 0.10 and to the neutron star mass, at which the direct Urca process is open, M_D = (1.6 - 1.8) M_sun.

Thermal evolution of neutron stars is known to depend on the properties of superdense matter in neutron star cores. We suggest a statistical analysis of isolated cooling middle-aged neutron stars and old transiently accreting quasi-stationary neutron stars warmed up by deep crustal heating in low-mass X-ray binaries. The method is based on simulations of the evolution of stars of different masses and on averaging the results over respective mass distributions. This gives theoretical distributions of isolated neutron stars in the surface temperature--age plane and of accreting stars in the photon thermal luminosity--mean mass accretion rate plane to be compared with observations. This approach permits to explore not only superdense matter but also the mass distributions of isolated and accreting neutron stars. We show that the observations of these stars can be reasonably well explained by assuming the presence of the powerful direct Urca process of neutrino emission in the inner cores of massive stars, introducing a slight broadening of the direct Urca threshold (for instance, by proton superfluidity), and by tuning mass distributions of isolated and accreted neutron stars.

We consider models of cooling neutron stars with nucleon cores which possess moderately strong triplet-state superfluidity of neutrons. When the internal temperature drops below the maximum of the critical temperature over the core, $T_{\rm C}$, this superfluidity sets in. It produces a neutrino outburst due to Cooper pairing of neutrons which greatly accelerates the cooling. We show that the cooling of the star with internal temperature $T$ within $0.6\,T_{\rm C} \lesssim T \leq T_{\rm C}$ is described by analytic self-similar relations. A measurement of the effective surface temperature of the star and its decline, supplemented by assumptions on star's mass, radius and composition of heat-blanketing envelope, allows one to construct a family of cooling models parametrized by the value of $T_{\rm C}$. Each model reconstructs cooling history of the star including its neutrino emission level before neutron superfluidity onset and the intensity of Cooper pairing neutrinos. The results are applied to interpret the observations of the neutron star in the Cassiopeia A supernova remnant.

Based on the formalism by Kaminker et al. (Astron. Astrophys. 343 (1999) 1009) we derive an analytic approximation for neutrino-pair bremsstrahlung emissivity due to scattering of electrons by atomic nuclei in the neutron star crust of any realistic composition. The emissivity is expressed through generalized Coulomb logarithm which we fit by introducing an effective potential of electron-nucleus scattering. In addition, we study the conditions at which the neutrino bremsstrahlung in the crust is affected by strong magnetic fields. The results can be applied for modelling of many phenomena in neutron stars, such as thermal relaxation in young isolated neutron stars and in accreting neutron stars with overheated crust in soft X-ray transients.

Context: Central Compact Objects (CCOs) in supernova remnants are isolated thermally emitting neutron stars (NSs). They are most probably characterized by a magnetic field strength that is roughly two orders of magnitude lower than that of most of the radio and accreting pulsars. The thermal emission of CCOs can be modeled to obtain constraints on the physical parameters of the star such as its mass, radius, effective temperature, and chemical composition. Aims: The CCO in HESS, J1731-347 is one of the brightest objects of this class. We present our analysis of two new XMM-Newton observations of the source which increase the total exposure time of the by a factor of five compared to the analyses presented before. Methods: We use our numerical spectral models for carbon and hydrogen atmospheres to fit the spectrum of the CCO. From our fits, we derive constraints on the physical parameters of the emitting star. We also use the new data to derive new upper limits on pulsations and to confirm the absence of a long-term variability. Results: The analysis shows that atmosphere spectral models are clearly preferred by the fit over the blackbody function. Under the assumption that the X-ray emission is uniformly produced by the entire star surface, hydrogen atmosphere models lead to uncomfortably large distances of the CCO, above 7-8 kpc. On the other hand, the carbon atmosphere model formally excludes distances above 5-6 kpc and is compatible with the source located in the Scutum-Crux (3 kpc) or Norma-Cygnus (4.5 kpc) Galactic spiral arm. We provide and discuss the corresponding confidence contours in the NS mass-radius plane. The measured effective temperature indicates that the NS is exceptionally hot for the estimated age of 30 kyr. We discuss possible cooling scenarios to explain this property, as well as possible additional constraints on the star mass and radius from cooling theory.

We calculate interdiffusion coefficients in a two-component, weakly or strongly coupled ion plasma (gas or liquid, composed of two ion species immersed into a neutralizing electron background). We use an effective potential method proposed recently by Baalrud and Daligaut [PRL, 110, 235001, (2013)]. It allows us to extend the standard Chapman-Enskog procedure of calculating the interdiffusion coefficients to the case of strong Coulomb coupling. We compute binary diffusion coefficients for several ionic mixtures and fit them by convenient expressions in terms of the generalized Coulomb logarithm. These fits cover a wide range of plasma parameters spanning from weak to strong Coulomb couplings. They can be used to simulate diffusion of ions in ordinary stars as well as in white dwarfs and neutron stars.

Using 1D and 2D cooling codes we study thermal emission from neutron stars with steady state internal heaters of various intensities and geometries (blobs or spherical layers) located at different depths in the crust. The generated heat tends to propagate radially, from the heater down to the stellar core and up to the surface; it is also emitted by neutrinos. In local regions near the heater the results are well described with the 1D code. The heater's region projects onto the stellar surface forming a hot spot. There are two heat propagation regimes. In the first, conduction outflow regime (realized at heat rates $H_0 \lesssim 10^{20}$ erg cm$^{-3}$ s$^{-1}$ or temperatures $T_\mathrm{h} \lesssim 10^9$ K in the heater) the thermal surface emission of the star depends on the heater's power and neutrino emission in the stellar core. In the second, neutrino outflow regime ($H_0 \gtrsim 10^{20}$ erg cm$^{-3}$ s$^{-1}$ or $T_\mathrm{h} \gtrsim 10^9$ K) the surface thermal emission becomes independent of heater's power and the physics of the core. The largest (a few per cent) fraction of heat power is carried to the surface if the heater is in the outer crust and the heat regime is intermediate. The results can be used for modeling young cooling neutron stars (prior to the end of internal thermal relaxation), neutron stars in X-ray transients, magnetars and high-$B$ pulsars, as well as merging neutron stars.

We analyze enhancement of thermonuclear fusion reactions due to strong plasma screening in dense matter using a simple electron drop model. The model assumes fusion in a potential that is screened by an effective electron cloud around colliding nuclei (extended Salpeter ion-sphere model). We calculate the mean field screened Coulomb potentials for atomic nuclei with equal and nonequal charges, appropriate astrophysical S factors, and enhancement factors of reaction rates. As a byproduct, we study analytic behavior of the screening potential at small separations between the reactants. In this model, astrophysical S factors depend not only on nuclear physics but on plasma screening as well. The enhancement factors are in good agreement with calculations by other methods. This allows us to formulate the combined, pure analytic model of strong plasma screening in thermonuclear reactions. The results can be useful for simulating nuclear burning in white dwarfs and neutron stars.

We analyze diffusion equations in strongly coupled Coulomb mixtures of ions in dense stellar matter. Strong coupling of ions in the presence of gravitational forces and electric fields (induced by plasma polarization in the presence of gravity) produces a specific diffusion current which can separate ions with the same A/Z (mass to charge number) ratios but different Z. This Coulomb separation of ions can be important for the evolution of white dwarfs and neutron stars.

The thermal evolution of young neutron stars (NSs) reflects the neutrino emission properties of their cores. Heinke et al. (2010) measured a 3.6+/-0.6% decay in the surface temperature of the Cassiopeia A (Cas A) NS between 2000 and 2009, using archival data from the Chandra X-ray Observatory ACIS-S detector in Graded mode. Page et al. (2011) and Shternin et al. (2011) attributed this decay to enhanced neutrino emission from a superfluid neutron transition in the core. Here we test this decline, combining analysis of the Cas A NS using all Chandra X-ray detectors and modes (HRC-S, HRC-I, ACIS-I, ACIS-S in Faint mode, and ACIS-S in Graded mode) and adding a 2012 May ACIS-S Graded mode observation, using the most current calibrations (CALDB 4.5.5.1). We measure the temperature changes from each detector separately and test for systematic effects due to the nearby filaments of the supernova remnant. We find a 0.92%-2.0% decay over 10 years in the effective temperature, inferred from HRC-S data, depending on the choice of source and background extraction regions, with a best-fit decay of 1.0+/-0.7%. In comparison, the ACIS-S Graded data indicate a temperature decay of 3.1%-5.0% over 10 years, with a best-fit decay of 3.5+/-0.4%. Shallower observations using the other detectors yield temperature decays of 2.6+/-1.9% (ACIS-I), 2.1+/-1.0% (HRC-I), and 2.1+/-1.9% (ACIS-S Faint mode) over 10 years. Our best estimate indicates a decline of 2.9+/-0.9 (stat) +1.6/-0.3 (sys) % over 10 years. The complexity of the bright and varying supernova remnant background makes a definitive interpretation of archival Cas A Chandra observations difficult. A temperature decline of 1-3.5% over 10 years would indicate extraordinarily fast cooling of the NS that can be regulated by superfluidity of nucleons in the stellar core.

Free energy, internal energy, and specific heat for each of the three phonon spectrum branches of a magnetized Coulomb crystal with body-centered cubic lattice are calculated by numerical integration over the Brillouin zone in the range of magnetic fields $B$ and temperatures $T$, such that $0 \le \omega_{\rm B}/\omega_{\rm p}\le 10^3$ and $10^{-4} \le T/T_{\rm p} \le 10^4$. In this case, $\omega_{\rm B}$ is the ion cyclotron frequency, $\omega_{\rm p}$ and $T_{\rm p}$ are the ion plasma frequency and plasma temperature, respectively. The results of numerical calculations are approximated by simple analytical formulas. For illustration, these formulas are used to analyze the behavior of the heat capacity in the crust of a neutron star with strong magnetic field. Thermodynamic functions of magnetized neutron star crust are needed for modeling various observational phenomena in magnetars and high magnetic field pulsars.

We review the history of neutron star physics in the 1930s that is related to L. Landau. According to recollections of Rosenfeld (1974, Proc. 16th Solvay Conference on Physics, p. 174), Landau improvised the concept of neutron stars in a discussion with Bohr and Rosenfeld just after the news of the discovery of the neutron reached Copenhagen in February 1932. We present arguments that the discussion took place in March 1931, before the discovery of the neutron, and that they in fact discussed the paper written by Landau in Zurich in February 1931 but not published until February 1932 (Phys. Z. Sowjetunion, 1, 285). In his paper Landau mentioned the possible existence of dense stars which look like one giant nucleus; this can be regarded as an early theoretical prediction or anticipation of neutron stars, prior to the discovery of the neutron. The coincidence of the dates of the neutron's discovery and the paper's publication has led to an erroneous association of the paper with the discovery of the neutron. In passing, we outline the contribution of Landau to the theory of white dwarfs and to the hypothesis of stars with neutron cores.

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.

Contrary to what is claimed by Ferrer et al. [Phys. Rev. C 82, 065802 (2010)], the magnetic field of a neutron star cannot exceed 10^{19} G and the thermodynamic pressure of dense magnetized fermion gas is isotropic.

We present a new study of the X-ray spectral properties of the Crab Pulsar. The superb angular resolution of the Chandra X-ray Observatory enables distinguishing the pulsar from the surrounding nebulosity. Analysis of the spectrum as a function of pulse phase allows the least-biased measure of interstellar X-ray extinction due primarily to photoelectric absorption and secondarily to scattering by dust grains in the direction of the Crab Nebula. We modify previous findings that the line-of-sight to the Crab is under-abundant in oxygen and provide measurements with improved accuracy and less bias. Using the abundances and cross sections from Wilms, Allen & McCray (2000) we find [O/H] = $(5.28 \pm 0.28)\times10^{-4}$ ($4.9 \times10^{-4}$ is solar abundance). We also measure for the first time the impact of scattering of flux out of the image by interstellar grains. We find $\tau_{\rm scat} = 0.147 \pm 0.043$. Analysis of the spectrum as a function of pulse phase also measures the X-ray spectral index even at pulse minimum --- albeit with increasing statistical uncertainty. The spectral variations are, by and large, consistent with a sinusoidal variation. The only significant variation from the sinusoid occurs over the same phase range as some rather abrupt behavior in the optical polarization magnitude and position angle. We compare these spectral variations to those observed in Gamma-rays and conclude that our measurements are both a challenge and a guide to future modeling and will thus eventually help us understand pair cascade processes in pulsar magnetospheres. The data were also used to set new, and less biased, upper limits to the surface temperature of the neutron star for different models of the neutron star atmosphere.

According to recent results of Ho & Heinke (2009) and Heinke & Ho (2010), the Cassiopeia A supernova remnant contains a young neutron star which has carbon atmosphere and shows noticeable decline of the effective surface temperature. We report a new (November 2010) Chandra observation which confirms the previously reported decline rate. The decline is naturally explained if neutrons have recently become superfluid (in triplet-state) in the NS core, producing a splash of neutrino emission due to Cooper pair formation (CPF) process that currently accelerates the cooling. This scenario puts stringent constraints on poorly known properties of NS cores: on density dependence of the temperature $T_\mathrm{cn}(\rho)$ for the onset of neutron superfluidity [$T_\mathrm{cn}(\rho)$ should have a wide peak with maximum $\approx (7-9)\times 10^8$ K], on the reduction factor $q$ of CPF process by collective effects in superfluid matter ($q > 0.4$), and on the intensity of neutrino emission before the onset of neutron superfluidity (30--100 times weaker than the standard modified Urca process). This is serious evidence for nucleon superfluidity in NS cores that comes from observations of cooling NSs.

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

We explore the thermal state of the neutron star in the Cassiopeia A supernova remnant using the recent result of Ho & Heinke (Nature, 462, 71 (2009)) that the thermal radiation of this star is well-described by a carbon atmosphere model and the emission comes from the entire stellar surface. Starting from neutron star cooling theory, we formulate a robust method to extract neutrino cooling rates of thermally relaxed stars at the neutrino cooling stage from observations of thermal surface radiation. We show how to compare these rates with the rates of standard candles -- stars with non-superfluid nucleon cores cooling slowly via the modified Urca process. We find that the internal temperature of standard candles is a well-defined function of the stellar compactness parameter $x=r_g/R$, irrespective of the equation of state of neutron star matter ($R$ and $r_g$ are circumferential and gravitational radii, respectively). We demonstrate that the data on the Cassiopeia A neutron star can be explained in terms of three parameters: $f_\ell$, the neutrino cooling efficiency with respect to the standard candle; the compactness $x$; and the amount of light elements in the heat blanketing envelope. For an ordinary (iron) heat blanketing envelope or a low-mass ($\lesssim 10^{-13}\,M_\odot$) carbon envelope, we find the efficiency $f_\ell \sim 1$ (standard cooling) for $x \lesssim 0.5$ and $f_\ell \sim 0.02$ (slower cooling) for a maximum compactness $x\approx 0.7$. A heat blanket containing the maximum mass ($\sim 10^{-8}\,M_\odot$) of light elements increases $f_\ell$ by a factor of 50. We also examine the (unlikely) possibility that the star is still thermally non-relaxed.

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 study plasma effects on radiative transitions (e.g., decay of excited states of atoms or atomic nuclei) in a dense plasma at the transition frequencies $\omega \lesssim \omega_p$ (where $\omega_p$ is the electron plasma frequency). The decay goes through four channels -- the emission of real transverse and longitudinal plasmons as well as the emission of virtual transverse and longitudinal plasmons with subsequent absorption of such plasmons by the plasma. The emission of real plasmons dies out at $\omega \leq \omega_p$, but the processes with virtual plasmons strongly enhance the radiative decay. Applications of these results to radiative processes in white dwarf cores and neutron star envelopes are discussed.

We study the thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in an internal layer. We focus on the effect of magnetized (B > 10^{14} G) non-accreted and accreted outermost envelopes composed of different elements, from iron to hydrogen or helium. We discuss a combined effect of thermal conduction and neutrino emission in the outer neutron star crust and calculate the cooling of magnetars with a dipole magnetic field for various locations of the heat layer, heat rates and magnetic field strengths. Combined effects of strong magnetic fields and light-element composition simplify the interpretation of magnetars in our model: these effects allow one to interpret observations assuming less extreme (therefore, more realistic) heating. Massive magnetars, with fast neutrino cooling in their cores, can have higher thermal surface luminosity.

We calculate the shear viscosity $\eta = \eta_{e\mu}+\eta_{n}$ in a neutron star core composed of nucleons, electrons and muons ($\eta_{e\mu}$ being the electron-muon viscosity, mediated by collisions of electrons and muons with charged particles, and $\eta_{n}$ the neutron viscosity, mediated by neutron-neutron and neutron-proton collisions). Deriving $\eta_{e\mu}$, we take into account the Landau damping in collisions of electrons and muons with charged particles via the exchange of transverse plasmons. It lowers $\eta_{e\mu}$ and leads to the non-standard temperature behavior $\eta_{e\mu}\propto T^{-5/3}$. The viscosity $\eta_{n}$ is calculated taking into account that in-medium effects modify nucleon effective masses in dense matter. Both viscosities, $\eta_{e\mu}$ and $\eta_{n}$, can be important, and both are calculated including the effects of proton superfluidity. They are presented in the form valid for any equation of state of nucleon dense matter. We analyze the density and temperature dependence of $\eta$ for different equations of state in neutron star cores, and compare $\eta$ with the bulk viscosity in the core and with the shear viscosity in the crust.

We review current state of neutron star cooling theory and discuss the prospects to constrain the equation of state, neutrino emission and superfluid properties of neutron star cores by comparing the cooling theory with observations of thermal radiation from isolated neutron stars.

Strong (B >> 10^9 G) and superstrong (B > 10^{14} G) magnetic fields profoundly affect many thermodynamic and kinetic characteristics of dense plasmas in neutron star envelopes. In particular, they produce strongly anisotropic thermal conductivity in the neutron star crust and modify the equation of state and radiative opacities in the atmosphere, which are major ingredients of the cooling theory and spectral atmosphere models. As a result, both the radiation spectrum and the thermal luminosity of a neutron star can be affected by the magnetic field. We briefly review these effects and demonstrate the influence of magnetic field strength on the thermal structure of an isolated neutron star, putting emphasis on the differences brought about by the superstrong fields and high temperatures of magnetars. For the latter objects, it is important to take proper account of a combined effect of the magnetic field on thermal conduction and neutrino emission at densities \rho > 10^{10} g cm^{-3}. We show that the neutrino emission puts a B-dependent upper limit on the effective surface temperature of a cooling neutron star.

We simulate the cooling of the neutron star in the X-ray transient KS 1731-260 after the source returned to quiescence in 2001 from a long (>~ 12.5 yr) outburst state. We show that the cooling can be explained assuming that the crust underwent deep heating during the outburst stage. In our best theoretical scenario the neutron star has no enhanced neutrino emission in the core, and its crust is thin, superfluid, and has the normal thermal conductivity. The thermal afterburst crust-core relaxation in the star may be not over.

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 calculate the thermal conductivity of electrons and muons kappa_{e-mu} produced owing to electromagnetic interactions of charged particles in neutron star cores and show that these interactions are dominated by the exchange of transverse plasmons (via the Landau damping of these plasmons in nonsuperconducting matter and via a specific plasma screening in the presence of proton superconductivity). For normal protons, the Landau damping strongly reduces kappa_{e-mu} and makes it temperature independent. Proton superconductivity suppresses the reduction and restores the Fermi-liquid behavior kappa_{e-mu} ~ 1/T. Comparing with the thermal conductivity of neutrons kappa_n, we obtain kappa_{e-mu}> kappa_n for T>2 GK in normal matter and for any T in superconducting matter with proton critical temperatures T_c>3e9 K. The results are described by simple analytic formulae.

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.

We calculate the thermal conductivity of electrons produced by electron-electron Coulomb scattering in a strongly degenerate electron gas taking into account the Landau damping of transverse plasmons. The Landau damping strongly reduces this conductivity in the domain of ultrarelativistic electrons at temperatures below the electron plasma temperature. In the inner crust of a neutron star at temperatures T < 1e7 K this thermal conductivity completely dominates over the electron conductivity due to electron-ion (electron-phonon) scattering and becomes competitive with the the electron conductivity due to scattering of electrons by impurity ions.

We study thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in a spherical internal layer. We explore the location of this layer as well as the heating rate that could explain high observable thermal luminosities of magnetars and would be consistent with the energy budget of neutron stars. We conclude that the heat source should be located in an outer magnetar's crust, at densities rho < 5e11 g/cm^3, and should have the heat intensity of the order of 1e20 erg/s/cm^3. Otherwise the heat energy is mainly emitted by neutrinos and cannot warm up the surface.

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.

We study the "minimal" cooling scenario of superfluid neutron stars with nucleon cores, where the direct Urca process is forbidden and the enhanced cooling is produced by the neutrino emission due to Cooper pairing of neutrons. Extending our previous consideration (Gusakov et al. 2004a), we include the effects of accreted envelopes of light elements. We employ phenomenological density-dependent critical temperatures T_{cp}(\rho) and T_{cnt}(\rho) of singlet-state proton and triplet-state neutron pairing in a stellar core, as well as the critical temperature T_{cns}(\rho) of singlet-state neutron pairing in a stellar crust. We show that the presence of accreted envelopes simplifies the interpretation of observations of thermal radiation from isolated neutron stars in the scenario of Gusakov et al. (2004a) and widens the class of models for nucleon superfluidity in neutron star interiors consistent with the observations.

We have derived a set of equations to describe the thermal evolution of a neutron star which undergoes small-amplitude radial pulsations. We have taken into account, in the frame of the General Theory of Relativity, the pulsation damping due to the bulk and shear viscosity and the accompanying heating of the star. The neutrino emission of a pulsating non-superfluid star and its heating due to the bulk viscosity are calculated assuming that both processes are determined by the non-equilibrium modified Urca process. Analytical and numerical solutions to the set of equations of the stellar evolution are obtained for linear and strongly non-linear deviations from beta-equilibrium. It is shown that a pulsating star may be heated to very high temperatures, while the pulsations damp very slowly with time (a power law damping for 100-1000 years), as long as the damping is determined by the bulk viscosity. The contribution of the shear viscosity to the damping becomes important in a rather cool star with a low pulsation energy.

We study the cooling of superfluid neutron stars whose cores consist of nucleon matter with the Akmal-Pandharipande-Ravenhall equation of state. This equation of state opens the powerful direct Urca process of neutrino emission in the interior of most massive neutron stars. Extending our previous studies (Gusakov et al. 2004a, Kaminker et al. 2005), we employ phenomenological density-dependent critical temperatures T_{cp}(\rho) of strong singlet-state proton pairing (with the maximum T_{cp}^{max} \sim 7e9 K in the outer stellar core) and T_{cnt}(\rho) of moderate triplet-state neutron pairing (with the maximum T_{cnt}^{max} \sim 6e8 K in the inner core). Choosing properly the position of T_{cnt}^{max} we can obtain a representative class of massive neutron stars whose cooling is intermediate between the cooling enhanced by the neutrino emission due to Cooper pairing of neutrons in the absence of the direct Urca process and the very fast cooling provided by the direct Urca process non-suppressed by superfluidity.

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.

We discuss pycnonuclear burning of highly exotic atomic nuclei in deep crusts of neutron stars, at densities up to 1e13 g/cc. As an application, we consider pycnonuclear burning of matter accreted on a neutron star in a soft X-ray transient (SXT, a compact binary containing a neutron star and a low-mass companion). The energy released in this burning, while the matter sinks into the stellar crust under the weight of newly accreted material, is sufficient to warm up the star and initiate neutrino emission in its core. The surface thermal radiation of the star in quiescent states becomes dependent of poorly known equation of state (EOS) of supranuclear matter in the stellar core, which gives a method to explore this EOS. Four qualitatively different model EOSs are tested against observations of SXTs. They imply different levels of the enhancement of neutrino emission in massive neutron stars by (1) the direct Urca process in nucleon/hyperon matter; (2) pion condensates; (3) kaon condensates; (4) Cooper pairing of neutrons in nucleon matter with the forbidden direct Urca process. A low level of the thermal quiescent emission of two SXTs, SAX J1808.4-3658 and Cen X-4, contradicts model (4). Observations of SXTs test the same physics of dense matter as observations of thermal radiation from cooling isolated neutron stars, but the data on SXTs are currently more conclusive.

We investigate the combined effect of neutron and proton superfluidities on the cooling of neutron stars whose cores consist of nucleons and electrons. We consider singlet-state pairing of protons and triplet-state pairing of neutrons in the cores of neutron stars. The critical superfluid temperatures T_c are assumed to depend on the density of matter. We study two types of neutron pairing with different components of the total angular momentum of Cooper pairs along the quantization axis (|m_J| =0 or 2). Our calculations are compared with observations of thermal emission from isolated neutron stars. We show that the observations can be interpreted by using two classes of superfluidity models: (1) strong proton superfluidity with a maximum critical temperature in the stellar core T_c^{max} > 4 \times 10^9 K and weak neutron superfluidity of any type (T_c^{max} < 2 \times 10^8 K); (2) strong neutron superfluidity (pairing with |m_J|=0) and weak proton superfluidity. The two types of models reflect an approximate symmetry with respect to an interchange of the critical temperatures of neutron and proton pairing.

The impact of nuclear physics theories on cooling of isolated neutron stars is analyzed. Physical properties of neutron star matter important for cooling are reviewed such as composition, the equation of state, superfluidity of various baryon species, neutrino emission mechanisms. Theoretical results are compared with observations of thermal radiation from neutron stars. Current constraints on theoretical models of dense matter, derived from such a comparison, are formulated.

We briefly review recent theoretical studies of the effects of accreted envelopes, magnetic fields and crustal superfluidity on the cooling of neutron stars. These effects are especially important for slowly cooling low-mass neutron stars, where direct Urca process of neutrino emission is forbidden. The effects are useful for interpretation of observations of several isolated middle-aged neutron stars.

We show that the direct Urca process of neutrino emission is allowed in two possible phases of nonspherical nuclei (inverse cylinders and inverse spheres) in the mantle of a neutron star near the crust-core interface. The process is open because neutrons and protons move in a periodic potential created by inhomogeneous nuclear structures. In this way the nucleons acquire large quasimomenta needed to satisfy momentum-conservation in the neutrino reaction. The appropriate neutrino emissivity in a nonsuperfluid matter is about 2--3 orders of magnitude higher than the emissivity of the modified Urca process in the stellar core. The process may noticeably accelerate the cooling of low-mass neutron stars.

We simulate cooling of superfluid neutron stars with nucleon cores where direct Urca process is forbidden. We adopt density dependent critical temperatures $T_{cp}(\rho)$ and $T_{cn}(\rho)$ of singlet-state proton and triplet-state neutron pairing in a stellar core and consider a strong proton pairing (with maximum $T_{cp}^{max} \ga 5 \times 10^9$ K) and a moderate neutron pairing ($T_{cn}^{max} \sim 6 \times 10^8$ K). When the internal stellar temperature $T$ falls below $T_{cn}^{max}$, the neutrino luminosity $L_{CP}$ due to Cooper pairing of neutrons behaves $\propto T^8$, just as that produced by modified Urca process (in a non-superfluid star) but is higher by about two orders of magnitude. In this case the Cooper-pairing neutrino emission acts like an enhanced cooling agent. By tuning the density dependence $T_{cn}(\rho)$ we can explain observations of cooling isolated neutron stars in the scenario in which direct Urca process or similar process in kaon/pion condensed or quark matter are absent.

Observation of cooling neutron stars can potentially provide information about the states of matter at supernuclear densities. We review physical properties important for cooling such as neutrino emission processes and superfluidity in the stellar interior, surface envelopes of light elements due to accretion of matter and strong surface magnetic fields. The neutrino processes include the modified Urca process, and the direct Urca process for nucleons and exotic states of matter such as a pion condensate, kaon condensate, or quark matter. The dependence of theoretical cooling curves on physical input and observations of thermal radiation from isolated neutron stars are described. The comparison of observation and theory leads to a unified interpretation in terms of three characteristic types of neutron stars: high-mass stars which cool primarily by some version of the direct Urca process; low-mass stars, which cool via slower processes; and medium-mass stars, which have an intermediate behavior. The related problem of thermal states of transiently accreting neutron stars with deep crustal burning of accreted matter is discussed in connection with observations of soft X-ray transients.

We calculate the thermal structure and quiescent thermal luminosity of accreting neutron stars (warmed by deep crustal heating in accreted matter) in soft X-ray transients (SXTs). We consider neutron stars with nucleon and hyperon cores and with accreted envelopes. It is assumed that an envelope has an outer helium layer (of variable depth) and deeper layers of heavier elements, either with iron or with much heavier nuclei (of atomic weight A > 100) on the top (Haensel & Zdunik 1990, 2003, astro-ph/0305220). The relation between the internal and surface stellar temperatures is obtained and fitted. The quiescent luminosity of the hottest (low-mass) and coldest (high-mass) neutron stars is calculated, together with the ranges of its possible variations due to variable thickness of the helium layer. The results are compared with observations of SXTs, particularly, containing the coldest (SAX J1808.4-3658) and the hottest (Aql X-1) neutron stars. The observations of SAX J1808.4-3658 in a quiescent state on March 24, 2001 (Campana et al. 2002, astro-ph/0206376) can be explained only if this SXT contains a massive neutron star with a nucleon/hyperon core; a hyperon core with a not too low fraction of electrons is preferable. Future observations may discriminate between the various models of hyperon/nucleon dense matter. The thermal emission of SAX J1808.4-3658 is also sensitive to the models of plasma ionization in the outermost surface layers and can serve for testing such models.

We study the thermal structure of neutron stars with magnetized envelopes composed of accreted material, using updated thermal conductivities of plasmas in quantizing magnetic fields, as well as equation of state and radiative opacities for partially ionized hydrogen in strong magnetic fields. The relation between the internal and local surface temperatures is calculated and fitted by an analytic function of the internal temperature, magnetic field strength, angle between the field lines and the normal to the surface, surface gravity, and the mass of the accreted material. The luminosity of a neutron star with a dipole magnetic field is calculated for various values of the accreted mass, internal temperature, and magnetic field strength. Using these results, we simulate cooling of superfluid neutron stars with magnetized accreted envelopes. We consider slow and fast cooling regimes, paying special attention to very slow cooling of low-mass superfluid neutron stars. In the latter case, the cooling is strongly affected by the combined effect of magnetized accreted envelopes and neutron superfluidity in the stellar crust. Our results are important for interpretation of observations of isolated neutron stars hottest for their age, such as RX J0822-43 and PSR B1055-52.

We study thermal states of transiently accreting neutron stars (with mean accretion rates $\dot{M} \sim 10^{-14}-10^{-9}$ M$_\odot$ yr$^{-1}$) determined by the deep crustal heating of accreted matter sinking into stellar interiors. We formalize a direct correspondence of this problem to the problem of cooling neutron stars. Using a simple toy model we analyze the most important factors which affect the thermal states of accreting stars: a strong superfluidity in the cores of low-mass stars and a fast neutrino emission (in nucleon, pion-condensed, kaon-condensed, or quark phases of dense matter) in the cores of high-mass stars. We briefly compare the results with the observations of soft X-ray transients in quiescence. If the upper limit on the quiescent thermal luminosity of the neutron star in SAX J1808.4-3658 (Campana et al. 2002) is associated with the deep crustal heating, it favors the model of nucleon neutron-star cores with switched-on direct Urca process.

A generic toy model of a cooling neutron star (NS) is used to analyze cooling of NSs with nucleon and exotic composition of the cores. The model contains the parameters which specify the levels of slow and enhanced neutrino emission as well as the lower and upper densities of the layer where the slow emission transforms into the enhanced one. The prospects to constrain these parameters from the present and future observations of isolated middle-aged NSs are discussed.

The cooling theory of isolated neutron stars is reviewed. The main cooling regulators are discussed, first of all, operation of direct Urca process (or similar processes in exotic phases of dense matter) and superfluidity in stellar interiors. The prospects to constrain gross parameters of supranuclear matter in neutron-star interiors by confronting cooling theory with observations of isolated neutron stars are outlined. A related problem of thermal states of transiently accreting neutron stars with deep crustal heating of accreted matter is discussed in application to soft X-ray transients.

We outline the phenomenon of deep crustal heating in transiently accreting neutron stars. It is produced by nuclear transformations (mostly, by pycnonuclear reactions) in accreted matter while this matter sinks to densities rho > 10^{10} g/cc under the weight of freshly accreted material. We consider then thermal states of transiently accreting neutron stars (with mean mass accretion rates \dot{M}=(10^{-14}-10^{-9}) M_\odot/yr) determined by deep crustal heating. In a simplified fashion we study how the thermal flux emergent from such stars depends on the properties of superdense matter in stellar interiors. We analyze the most important regulators of the thermal flux: strong superfluidity in the cores of low-mass stars and fast neutrino emission (in nucleon, pion-condensed, kaon-condensed, or quark phases of dense matter) in the cores of high-mass stars. We compare the results with observations of soft X-ray transients in quiescent states.

We extend our modeling of cooling of superfluid neutron stars by including the production of muons in the core, in addition to neutrons, protons, and electrons. The results are confronted with observations of middle-aged isolated NSs. Muons have little effect on the hydrostatic structure of NSs, on the slow cooling of low-mass NSs (RX J0822--43 and PSR 1055--52 in our model) and on the rapid cooling of massive NSs. They affect, however, the moderately fast cooling of medium-mass NSs (1E 1207--52, RX J0002+62, PSR 0656+14, Vela, and Geminga) and shift appreciably the mass range of these NSs to lower masses, which is important for correct interpretation of the observations. Moreover, the effects of muons can accurately be reproduced by a simple renormalization of NS models with no muons in the NS cores.

The adiabatic index Gamma_1 for perturbations of dense matter is studied under various physical conditions which can prevail in neutron star cores. The dependence of Gamma_1 on the composition of matter (in particular, on the presence of hyperons), on the stellar pulsation amplitude, and on the baryon superfluidity is analyzed. Timescales of damping of stellar pulsations are estimated at different compositions, temperatures, and pulsation amplitudes. Damping of pulsations by bulk viscosity in the neutron-star cores can prevent the stars to pulsate with relative amplitudes > (1-15)% (depending on the composition of matter).

The upper limit of the effective surface temperature of the neutron star (NS) PSR J0205+6449 in the supernova remnant 3C 58 obtained recently by Slane et al. (2002) is analyzed using a modern theory of NS cooling (Kaminker et al. 2002). The observational limit can be explained by cooling of a superfluid NS with the core composed of neutrons, protons and electrons, where direct Urca process is forbidden. However, combined with the data on the surface temperatures of other isolated NSs, it gives evidence (emphasized by Slane et al.) that direct Urca process is open in the inner cores of massive NSs. This evidence turns out to be less stringent than the evidence provided by the well known observations of Vela and Geminga.

Cooling of neutron stars (NSs) with the cores composed of neutrons, protons, and electrons is analyzed. The main cooling regulators are discussed: opening of direct Urca process in a NS central kernel; superfluidity of nucleons in NS interiors; surface layers of light (accreted) elements; strong surface magnetic fields. An emphasis is paid on the cooling scenario with strong $^1$S$_0$ pairing of protons and weak $^3$P$_2$ pairing of neutrons in the NS core, as well as strong $^1$S$_0$ pairing of neutrons in the NS crust. The theory predicts three types of isolated cooling middle-aged NSs with distinctly different properties: low-mass, slowly cooling NSs; medium-mass, moderately cooling NSs; massive, rapidly cooling NSs. The theory is compared with observations of isolated NSs -- pulsars and radio quiet NSs in supernova remnants. The constraints on physical properties of NSs which can be inferred from such a comparison are outlined.

Cooling of neutron stars (NSs) with the cores composed of neutrons, protons, and electrons is simulated assuming $^1$S$_0$ pairing of neutrons in the NS crust, and also $^1$S$_0$ pairing of protons and weak $^3$P$_2$ pairing of neutrons in the NS core, and using realistic density profiles of the superfluid critical temperatures $T_{\rm c}(\rho)$. The theoretical cooling models of isolated middle-aged NSs can be divided into three main types. (I) {\it Low-mass}, {\it slowly cooling} NSs where the direct Urca process of neutrino emission is either forbidden or almost fully suppressed by the proton superfluidity. (II) {\it Medium-mass} NSs which show {\it moderate} cooling via the direct Urca process suppressed by the proton superfluidity. (III) {\it Massive} NSs which show {\it fast} cooling via the direct Urca process weakly suppressed by superfluidity. Confronting the theory with observations we treat RX J0822--43, PSR 1055--52 and RX J1856--3754 as slowly cooling NSs. To explain these sufficiently warm sources we need a density profile $T_{\rm c}(\rho)$ in the crust with a rather high and flat maximum and sharp wings. We treat 1E 1207--52, RX J0002+62, PSR 0656+14, Vela, and Geminga as moderately cooling NSs. We can determine their masses for a given model of proton superfluidity, $T_{\rm cp}(\rho)$, and the equation of state in the NS core. No rapidly cooling NS has been observed so far.

Bulk viscosity of neutron star cores containing hyperons is studied taking into account non-equilibrium weak process $n+n \rightleftharpoons p+\Sigma^-$. Rapid growth of the bulk viscosity within the neutron star core associated with switching on new reactions (modified Urca process, direct Urca process, hyperon reactions) is analyzed. The suppression of the bulk viscosity by superfluidity of baryons is considered and found out to be very important.

As shown recently by Kaminker et al. (2001), current observations of thermal emission of isolated middle-aged neutron stars (NSs) can be explained by cooling of NSs of different masses with the cores composed of neutrons, protons and electrons, assuming rather strong superfluidity (SF) of protons, weak triplet-state SF of neutrons and neglecting singlet-state SF of neutrons. We show that this explanation remains correct in the presence of singlet-state SF of neutrons in the NS crust and outermost core but under stringent constraints on the density profile of the SF critical temperature $T_{\rm cns}(\rho)$. In order to explain observations of (young and hot) RX J0822--43 and (old and warm) PSR 1055--52 and RX J1856--3754 as cooling not too massive NSs, the maximum $T_{\rm cns}^{\rm max}$ should be rather high ($\ga 5 \times 10^9$ K) and the decrease of $T_{\rm cns}(\rho)$ outside the maximum should be sharp. These results place important constraints on the models of nucleon SF in dense matter.

Phonon frequency moments and thermodynamic functions (electrostatic and vibrational parts of the free energy, internal energy, and heat capacity) are calculated for bcc and fcc Coulomb crystals in the harmonic approximation with a fractional accuracy within one part on 10^5. Temperature dependence of thermodynamic functions is fitted by analytic formulas with an accuracy of a few parts in 10^5. The static-lattice (Madelung) part of the free energy is calculated with an accuracy about 10^{-12}. The Madelung constant and frequency moments of hcp crystals are also computed.

The thermal structure of neutron stars with magnetized envelopes is studied using modern physics input. The relation between the internal (T_i) and local surface temperatures is calculated and fitted by analytic expressions for magnetic field strengths B from 0 to 10^{16} G and arbitrary inclination of the field lines to the surface. The luminosity of a neutron star with dipole magnetic field is calculated and fitted as a function of B, T_i, stellar mass and radius. In addition, we simulate cooling of neutron stars with magnetized envelopes. In particular, we analyse ultramagnetized envelopes of magnetars and also the effects of the magnetic field of the Vela pulsar on the determination of critical temperatures of neutron and proton superfluids in its core.

Cooling simulations of neutron stars (NSs) are performed assuming that stellar cores consist of neutrons, protons and electrons and using realistic density profiles of superfluid critical temperatures $T_{cn}(\rho)$ and $T_{cp}(\rho)$ of neutrons and protons. Taking a suitable profile of $T_{cp}(\rho)$ with maximum $\sim 5 \times 10^9$ K one can obtain smooth transition from slow to rapid cooling with increasing stellar mass. Adopting the same profile one can explain the majority of observations of thermal emission from isolated middle--aged NSs by cooling of NSs with different masses either with no neutron superfluidity in the cores or with a weak superfluidity, $T_{cn} < 10^8$ K. The required masses range from $\sim 1.2 M_\odot$ for (young and hot) RX J0822-43 and (old and warm) PSR 1055-52 and RX J1856-3754 to $\approx 1.45 M_\odot$ for the (colder) Geminga and Vela pulsars. Observations constrain the $T_{cn}(\rho)$ and $T_{cp}(\rho)$ profiles with respect to the threshold density of direct Urca process and maximum central density of NSs.

The diffusive thermal conductivity of neutrons in dense matter [$\rho \sim (1 - 8) \times 10^{14}$ g cm$^{-3}$] of neutron star cores is calculated. The contribution from neutron--neutron and neutron--proton collisions is taken into account. We use the transition probabilities calculated for symmetric dense nucleon matter on the basis of the Dirac--Brueckner approach to the in-medium effects and the Bonn model of bare nucleon--nucleon interaction. The diffusive thermal conductivity of neutrons in the presence of neutron and proton superfluidities is analyzed in a microscopic manner; the effects of superfluidity are shown to be significant. The low temperature behavior of the thermal conductivity appears to be extremely sensitive to the relation between critical temperatures of neutrons and protons. The results are fitted by simple analytic expressions. In combination with the formulae for the electron and muon thermal conductivities, obtained earlier, the present expressions provide a realistic description of the full diffusive thermal conductivity in the neutron star cores for normal and various superfluid phases.

We study the bulk viscosity in neutron star cores due to modified Urca processes involving nucleons, electrons and muons and analize its reduction by singlet-state or triplet-state superfluidity of nucleons. In combination with the results of our previous paper on the bulk viscosity due to direct Urca processes, a realistic description of the bulk viscosity in superfluid neutron star cores is obtained. Switching off direct Urca processes with decreasing density in a nonsuperfluid matter lowers the bulk viscosity by 3-5 orders of magnitude. The presense of muons opens additional source of bulk viscosity due to muon Urca processes and lowers the threshold density of the electron direct Urca process. The superfluidity may strongly reduce the bulk viscosity and affect thus damping of neutron star vibrations.

Matter of subnuclear density in the inner crusts of neutron stars consists of neutron-rich atomic nuclei immersed in strongly degenerate relativistic gas of electrons and strongly nonideal liquid of neutrons. We give a brief account of thermodynamic and kinetic properties of this matter which are greatly affected by Coulomb and nuclear interactions (a companion paper on the outer envelopes is astro-ph/0012316) and show how they can be studied from observations of thermal radiation of young (age < 100 yr) neutron stars (for a discussion of the latter possibility, see astro-ph/0012306).

Outer envelopes of neutron stars consist mostly of fully ionized, strongly coupled Coulomb plasmas characterized by typical densities about 10^4-10^{11} g/cc and temperatures about 10^4-10^9 K. Many neutron stars possess magnetic fields about 10^{11}-10^{14} G. Here we briefly review recent theoretical advances which allow one to calculate thermodynamic functions and electron transport coefficients for such plasmas with an accuracy required for theoretical interpretation of observations.

The internal properties of the neutron star crust can be probed by observing the epoch of thermal relaxation. After the supernova explosion, powerful neutrino emission quickly cools the stellar core, while the crust stays hot. The cooling wave then propagates through the crust, due to its finite thermal conductivity. When the cooling wave reaches the surface (age 10-100 yr), the effective temperature drops sharply from 250 eV to 30 or 100 eV, depending on the cooling model. The crust relaxation time is sensitive to the (poorly known) microscopic properties of matter of subnuclear density, such as the heat capacity, thermal conductivity, and superfluidity of free neutrons. We calculate the cooling models with the new values of the electron thermal conductivity in the inner crust, based on a realistic treatment of the shapes of atomic nuclei. Superfluid effects may shorten the relaxation time by a factor of 4. The comparison of theoretical cooling curves with observations provides a potentially powerful method of studying the properties of the neutron superfluid and highly unusual atomic nuclei in the inner crust.

We review the main neutrino emission mechanisms in neutron star crusts and cores. Among them are the well-known reactions such as the electron-positron annihilation, plasmon decay, neutrino bremsstrahlung of electrons colliding with atomic nuclei in the crust, as well as the Urca processes and neutrino bremsstrahlung in nucleon-nucleon collisions in the core. We emphasize recent theoretical achievements, for instance, band structure effects in neutrino emission due to scattering of electrons in Coulomb crystals of atomic nuclei. We consider the standard composition of matter (neutrons, protons, electrons, muons, hyperons) in the core, and also the case of exotic constituents such as the pion or kaon condensates and quark matter. We discuss the reduction of the neutrino emissivities by nucleon superfluidity, as well as the specific neutrino emission produced by Cooper pairing of the superfluid particles. We also analyze the effects of strong magnetic fields on some reactions, such as the direct Urca process and the neutrino synchrotron emission of electrons. The results are presented in the form convenient for practical use. We illustrate the effects of various neutrino reactions on the cooling of neutron stars. In particular, the neutrino emission in the crust is critical in setting the initial thermal relaxation between the core and the crust. Finally, we discuss the prospects of exploring the properties of supernuclear matter by confronting cooling simulations with observations of the thermal radiation from isolated neutron stars.

The bulk viscosity of the neutron star matter due to the direct Urca processes involving nucleons, electrons and muons is studied taking into account possible superfluidity of nucleons in the neutron star cores. The cases of singlet-state pairing or triplet-state pairing (without and with nodes of the superfluid gap at the Fermi surface) of nucleons are considered. It is shown that the superfluidity may strongly reduce the bulk viscosity. The practical expressions for the superfluid reduction factors are obtained. For illustration, the bulk viscosity is calculated for two models of dense matter composed of neutrons, protons,electrons and muons. The presence of muons affects the bulk viscosity due to the direct Urca reactions involving electrons and produces additional comparable contribution due to the direct Urca reactions involving muons. The results can be useful for studying damping of vibrations of neutron stars with superfluid cores.

The dynamic structure factor ${\tilde S}({\bf k},\omega)$ and the two-particle distribution function $g({\bf r},t)$ of ions in a Coulomb crystal are obtained in a closed analytic form using the harmonic lattice (HL) approximation which takes into account all processes of multi-phonon excitation and absorption. The static radial two-particle distribution function $g(r)$ is calculated for classical ($T \gtrsim \hbar \omega_p$, where $\omega_p$ is the ion plasma frequency) and quantum ($T \ll \hbar \omega_p$) body-centered cubic (bcc) crystals. The results for the classical crystal are in a very good agreement with extensive Monte Carlo (MC) calculations at $1.5 \lesssim r/a \lesssim 7$, where $a$ is the ion-sphere radius. The HL Coulomb energy is calculated for classical and quantum bcc and face-centered cubic crystals, and anharmonic corrections are discussed. The inelastic part of the HL static structure factor $S''(k)$, averaged over orientations of wave-vector {\bf k}, is shown to contain pronounced singularities at Bragg diffraction positions. The type of the singularities is different in classical and quantum cases. The HL method can serve as a useful tool complementary to MC and other numerical methods.

We study the heat capacity and neutrino emission reactions (direct and modified Urca processes, nucleon-nucleon bremsstrahlung, Cooper pairing of nucleons) in matter of supranuclear density of the neutron star cores with superfluid neutrons and protons. Various superfluidity types are analysed (singlet-state pairing and two types of triplet-state pairing, without and with nodes of the gap at a nucleon Fermi surface). The results are used for cooling simulations of isolated neutron stars. Both, the standard cooling and the cooling enhanced by the direct Urca process, are strongly affected by nucleon superfluidity. Comparison of cooling theory of isolated neutron stars with observations of their thermal radiation may give stringent constraints on the critical temperatures of the neutron and proton superfluidities in the neutron star cores.

New calculations of the thermal and electrical electron conductivities are performed for a broad range of physical parameters typical for envelopes of neutron stars and cores of white dwarfs. We consider stellar matter composed of astrophysically important chemical elements from H to Fe in the density range from 10^2-10^4 g/ccm up to 10^7-10^{10} g/ccm, where atoms are fully ionized and electrons are strongly degenerate. We have used modified ion structure factors suggested in physics/9811052. In the ion liquid, these modifications take into account, in an approximate way, instantaneous electron-band structures that reduce the electron-ion scattering rate. In crystallized matter, the new structure factors include multi-phonon processes important at temperatures not very much lower than the melting temperature T_m. The transport coefficients obtained differ significantly from those derived earlier in the important temperature range T_m/5 < T < 5 T_m. The results of our numerical calculations are fitted by analytical expressions convenient for astrophysical applications.

We calculate the neutrino energy emission rate due to singlet-state pairing of protons in the neutron star cores taking into account the relativistic correction to the non-relativistic rate. The non-relativistic rate is numerically small, and the relativistic correction appears to be about 10 -- 50 times larger. It plays thus the leading role, reducing great difference between the neutrino emissions due to pairing of protons and neutrons. The results are important for simulations of neutron star cooling.

Cooling of neutron stars (NSs) with superfluid cores is simulated taking into account neutrino emission produced by Cooper pairing of nucleons. The critical temperatures of neutron and proton superfluidities, $T_{cn}$ and $T_{cp}$, are assumed to be constant over the NS core, and treated as free parameters. They are constrained using the surface temperatures $T_s$ of isolated NSs (RX J0822-43, PSR 1055-52, 1E 1207-52, Vela, Geminga, PSR 0656+14, RX J0002+62), obtained by interpretation of observed thermal radiation either with black body spectrum or with hydrogen atmosphere models.

Neutrino-pair bremsstrahlung by relativistic degenerate electrons in a neutron-star crust at densities (10^9 - 1.5x10^{14}) g/cm^3 is analyzed. The processes taken into account are neutrino emission due to Coulomb scattering of electrons by atomic nuclei in a Coulomb liquid, and electron-phonon scattering and Bragg diffraction (the static-lattice contribution) in a Coulomb crystal. The static-lattice contribution is calculated including the electron band-structure effects for cubic Coulomb crystals of different types and also for the liquid crystal phases composed of rod- and plate-like nuclei in the neutron-star mantle (at 10^{14} - 1.5x10^{14} g/cm^3). The phonon contribution is evaluated with proper treatment of the multi-phonon processes which removes a jump in the neutrino bremsstrahlung emissivity at the melting point obtained in previous works. Below 10^{13} g/cm^3, the results are rather insensitive to the nuclear form factor, but results for the solid state near the melting point are affected significantly by the Debye-Waller factor and multi-phonon processes. At higher densities, the nuclear form factor becomes more significant. A comparison of the various neutrino generation mechanisms in neutron star crusts shows that electron bremsstrahlung is among the most important ones.

The neutrino energy emission rate due to formation of Cooper pairs of neutrons and protons in the superfluid cores of neutron stars is studied. The cases of singlet-state pairing with isotropic superfluid gap and triplet-state pairing with anisotropic gap are analysed. The neutrino emission due to singlet-state pairing of protons is found to be greatly suppressed with respect to the cases of singlet- and triplet-state pairings of neutrons. The neutrino emission due to pairing of neutrons is shown to be very important in the superfluid neutron-star cores with the standard neutrino luminosity and with the luminosity enhanced by the direct Urca process. It can greatly accelerate both, standard and enhanced, cooling of neutron stars with superfluid cores. This enables one to interpret the data on surface temperatures of six neutron stars, obtained by fitting the observed spectra with the hydrogen atmosphere models, by the standard cooling with moderate nucleon superfluidity.

The effect of the magnetic field on the energy loss rate in the direct Urca reactions is studied. The general expression for the neutrino emissivity at arbitrary magnetic field B is derived. The main emphasis is laid on a case, in which the field is not superstrong, and charged reacting particles (e and p) populate many Landau levels. The magnetic field keeps the process operative if Delta k / k_{Fn} < N_{Fp}^{-2/3} (N_{Fp} is the number of the Landau levels populated by protons and Delta k = k_{Fn}-k_{Fp}-k_{Fe}), that is beyond the well-known switch-on limit in the absence of the field, Delta k < 0. Cooling of magnetized neutron stars with strong neutron superfluid in the outer cores and nonsuperfluid inner cores is simulated. The magnetic field near the stellar center speeds up the cooling if the stellar mass M is slightly less than the minimum mass M_c at which the direct Urca reaction becomes allowed for B=0. If B=3x10^{16} G the affected mass range is M_c-M < 0.1M_c, while for B=3x10^{15} G the range is M_c-M < 0.015M_c. This may influence a theoretical interpretation of the observed thermal radiation as illustrated for the Geminga pulsar. The case of superstrong magnetic fields (B>10^{18} G), such that e and p populate only the lowest Landau levels is briefly outlined.

The dynamical structure factor of a Coulomb crystal of ions is calculated at arbitrary temperature below the melting point taking into account multi-phonon processes in the harmonic approximation. In a strongly coupled Coulomb ion liquid, the static structure factor is split into two parts, a Bragg-diffraction-like one, describing incipient long-range order structures, and an inelastic part corresponding to thermal ion density fluctuations. It is assumed that the diffractionlike scattering does not lead to the electron relaxation in the liquid phase. This assumption, together with the inclusion of multi-phonon processes in the crystalline phase, eliminates large discontinuities of the transport coefficients (jumps of the thermal and electric conductivities, as well as shear viscosity, reported previously) at a melting point.

The relationships between the effective surface (T_e) and internal temperatures of neutron stars (NSs) with and without accreted envelopes are calculated for T_e > 5\times10^4 K using new data on the equation of state and opacities in the outer NS layers. We examine various models of accreted layers (H, He, C, O shells produced by nuclear transformations in accreted matter). We employ new Opacity Library (OPAL) radiative opacities for H, He, and Fe. In the outermost NS layers, we implement the modern OPAL equation of state for Fe, and the Saumon-Chabrier equation of state for H and He. The updated thermal conductivities of degenerate electrons include the Debye-Waller factor for the electron-phonon scattering in solidified matter, while in liquid matter they include the contributions from electron-ion collisions (evaluated with non-Born corrections and with the ion structure factors in responsive electron background) and from the electron-electron collisions. For T_e < 10^{5.5} K, the electron conduction in non-degenerate layers of the envelope becomes important, reducing noticeably the temperature gradient. The accreted matter further decreases this gradient at T_e > 10^5 K. Even a small amount of accreted matter (with mass > 10^{-16} Msun) affects appreciably the NS cooling, leading to higher T_e at the neutrino cooling stage and to lower T_e at the subsequent photon stage.

We study the synchrotron emission of neutrino pairs by relativistic, degenerate electrons in strong magnetic fields. Particular attention is paid to the case in which the dominant contribution comes from one or several lowest cyclotron harmonics. Calculations are performed using the exact quantum formalism and the quasiclassical approach. Simple analytic fits to the neutrino synchrotron emissivity are obtained in the domain of magnetized, degenerate and relativistic electron gas provided the electrons populate either many Landau levels or the ground level alone. The significance of the neutrino synchrotron energy losses in the interiors of cooling neutron stars is discussed.

Practical expressions are derived for a rapid and accurate evaluation of electric and thermal conductivities and thermopower of degenerate relativistic electrons along quantizing magnetic fields in outer neutron star crusts. The electron Coulomb scattering on ions is considered in liquid matter, and on high-temperature phonons or charged impurities in solid matter. A semi-quantitative treatment of low-temperature phonons is also proposed. The transport properties are expressed through the energy dependent effective electron relaxation time, which is calculated using the theoretical formalism of the previous work (astro-ph/9603133). Accurate fitting expressions for the relaxation time are obtained, and an efficient energy averaging procedure proposed. A Fortran code, which calculates the longitudinal transport properties of degenerate electrons in strong magnetic fields for any parameters of dense stellar matter of practical interest, is available from the authors upon request.

We study the emission of neutrinos, resulting from the scattering of electrons off magnetic flux tubes (fluxoids) in the neutron star cores with superfluid (superconducting) protons. In the absence of proton superfluidity (T> T_{cp}), this process transforms into the well known electron synchrotron emission of neutrino pairs in a locally uniform magnetic field B, with the neutrino energy loss rate Q proportional to B^2 T^5. For temperatures T not much below T_{cp}, the synchrotron regime (Q \propto T^5) persists and the emissivity Q can be amplified by several orders of magnitude due to the appearance of the fluxoids and associated enhancement of the field within them. For lower T, the synchrotron regime transforms into the bremsstrahlung regime (Q \propto T^6) similar to the ordinary neutrino-pair bremsstrahlung of electrons which scatter off atomic nuclei. We calculate Q numerically and represent our results through a suitable analytic fit. In addition, we estimate the emissivities of two other neutrino-production mechanisms which are usually neglected -- neutrino-pair bremsstrahlung processes due to electron-proton and electron-electron collisions. We show that the electron-fluxoid and electron-electron scattering can provide the main neutrino production mechanisms in the neutron star cores with highly superfluid protons and neutrons at T < 5 10^8 K. The electron-fluxoid scattering is significant if the initial, locally uniform magnetic field B > 10^{13} G.

Recently, Bautista (1996, A&AS, 119, 105) reported new calculations of photoionization cross sections for Fe V, and compared them with the earlier results of Reilman & Manson (1979, ApJS, 40, 815) and Verner et al. (1993, Atomic Data Nucl. Data Tables, 55, 233). Bautista claimed that beyond 10 Ry the new cross sections ``converge well toward the results by Reilman and Manson'', whereas ``the calculations by Verner et al. still underestimate the cross section by almost factor of two''. We show that Bautista erroneously compared the total (summed over shells) cross sections from Reilman & Manson with the partial 3d-shell cross sections from Verner et al. Actually, the total cross sections from Verner et al. and Reilman & Manson agree within 3% at all energies.

Calculations are performed of the cooling of neutron stars with standard and enhanced neutrino energy losses in the presence of neutron and proton superfluidities in the stellar cores. The effects of superfluidity on the heat capacity and the neutrino luminosity produced by the direct and modified Urca processes and by the neutrino bremsstrahlung emission in nucleon-nucleon collisions are taken into account. The constraints are analyzed on the critical temperatures $\tn$ and $\tp$ of the transition of neutrons and protons to the superfluid state, which can be obtained from a comparison of observational data on thermal radiation from neutron stars with theoretical cooling curves for the standard and enhanced neutrino energy losses. Possible ranges of $\tn$ and $\tp$ for the pulsar Geminga are discussed.

Thermal and electric conductivities are calculated for degenerate electrons scattered by phonons in a crystal made of atomic nuclei. The exact phonon spectrum and the Debye--Waller factor are taken into account. Monte Carlo calculations are performed for body-centered cubic (bcc) crystals made of C, O, Ne, Mg, Si, S, Ca, and Fe nuclei in the density range from $10^3$ to $10^{11}$ g cm$^{-3}$ at temperatures lower than the melting temperature but higher than the temperature at which the Umklapp processes begin to be "frozen out". A simplified method of calculation is proposed, which makes it possible to describe the results in terms of simple analytic expressions, to extend these expressions to any species of nucleus, and to consider face-centered cubic (fcc) crystals. The kinetic coefficients are shown to depend tangibly on the lattice type. The results are applicable to studies of heat transfer and evolution of the magnetic field in the cores of white dwarfs and in the crusts of neutron stars. The thermal drift of the magnetic field in the crust of a neutron star is discussed.

Thermal and electric conductivities of relativistic degenerate electrons are calculated for the case when electrons scatter by phonons in Coulomb crystals made of spherical finite--size nuclei at densities $10^{11}$~g/cm$^3 < \rho < 10^{14}$~g/cm$^3$, corresponding to the inner crust of a neutron star. In combination with the results of the previous article (for lower $\rho$), simple unified fits are obtained which describe the kinetic coefficients in the range $10^3$~g/cm$^3 < \rho < 10^{14}$~g/cm$^3$, for matter with arbitrary nuclear composition. The results are valid for studying thermal evolution of neutron stars and evolution of their magnetic fields. The difference between the kinetic coefficients in the neutron star crust composed of ground state and accreted matters is analyzed. Thermal drift of the magnetic field in the neutron star crust is discussed.

Cooling of neutron stars with dipole magnetic fields is simulated using a realistic model of the anisotropic surface temperature distribution produced by magnetic fields. Suppression of the electron thermal conductivity of outer stellar layers across the field increases thermal isolation of these layers near the magnetic equator. Enhancement of the radiative and longitudinal electron thermal conductivities in quantizing magnetic fields reduces thermal isolation near the magnetic poles. The equatorial increase of the isolation is pronounced for $B \ga 10^{10}$ G, while the polar decrease -- for $B \ga 10^{12}$ G. The effects compensate partly each other, and the actual influence of the magnetic fields on the cooling is weaker than predicted by the traditional theories where the equatorial effects have been neglected.

Neutrino-pair bremsstrahlung radiation is considered due to electron--phonon scattering of degenerate, relativistic electrons in a lattice of spherical atomic nuclei in a neutron star crust. The neutrino energy generation rate is calculated taking into account exact spectrum of phonons, the Debye--Waller factor, and the nuclear form--factor in the density range from $10^7$~g~cm$^{-3}$ to $10^{14}$~g~cm$^{-3}$ at arbitrary nuclear composition for body-centered-cubic and face-centered-cubic Coulomb crystals. The results are fitted by a unified analytic expression. A comparison is given of the neutrino bremsstrahlung energy losses in a neutron star crust composed of ground state and accreted matter, in the solid and liquid phases.

Neutrino emissivity from the electron $\nu\bar\nu$ bremsstrahlung in the liquid layers of the neutron star crusts is studied. Nuclear composition of matter in neutron star crusts is considered for various scenarios of neutron star evolution. For the deep layers of the crust, the compositions of cold catalyzed matter, accreted matter and hot matter ($T \ga 5 \times 10^9$~K) are shown to be very different, and this implies differences in the neutrino emissivity at given density and temperature. Neutrino-pair bremsstrahlung, due to collisions of relativistic degenerate electrons in a Coulomb liquid of atomic nuclei, is considered. The neutrino energy loss rate is expressed in a simple form through a Coulomb logarithm $L$ -- a slowly varying function of density, temperature, and nuclear composition. Non-Born corrections, thermal width of electron Fermi sphere, and finite sizes of nuclei are taken into account. Implications for cooling of neutron stars are discussed.

We present a complete set of analytic fits to the non-relativistic photoionization cross sections for the ground states of atoms and ions of elements from H through Si, and S, Ar, Ca, and Fe. Near the ionization thresholds, the fits are based on the Opacity Project theoretical cross sections interpolated and smoothed over resonances. At higher energies, the fits reproduce calculated Hartree-Dirac-Slater photoionization cross sections.