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.

Using Metropolis method to compute path integrals, the energy of quantum strongly-coupled Coulomb liquid ($1 \leq \Gamma \leq 175$) composed of distinguishable atomic nuclei and uniform incompressible electron background is calculated from first principles. The range of considered temperatures and densities represents fully-ionized layers of white dwarfs and neutron stars. In particular, the results allow one to determine reliably the heat capacity of ions in dense fluid stellar matter, which is a crucial ingredient for modelling thermal evolution of compact degenerate stars.

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

Although crystallized neutron star crust is responsible for many fascinating observational phenomena, its actual microscopic structure in tremendous gravitational and magnetic fields is not understood. Here we show that in a non-uniform magnetic field, three-dimensional ionic Coulomb crystals comprising the crust may stretch or shrink while their electrostatic pressure becomes anisotropic. The pressure depends non-linearly on the magnitude of the stretch, so that a continuous magnetic field evolution may result in an abrupt crystal elongation or contraction. This may provide a trigger for magnetar activity. A phonon mode instability is revealed, which sets the limits of magnetic field variation beyond which the crystal is destroyed. These limits sometimes correspond to surprisingly large deformations. It is not known what happens to crust matter subject to a pressure anisotropy exceeding these limits. We hypothesize that the ion system then possesses a long-range order only in one or two dimensions, that is becomes a liquid crystal.

Corrections to elastic moduli, including the effective shear modulus, of a solid neutron star crust due to electron screening are calculated. At any given mass density, the crust is modelled as a body-centred cubic Coulomb crystal of fully ionized atomic nuclei of a single type with a polarizable charge-compensating electron background. Motion of the nuclei is neglected. The electron polarization is described by a simple Thomas-Fermi model of exponential electron screening. The results of numerical calculations are fitted by convenient analytic formulae. They should be used for precise neutron star oscillation modelling, a rapidly developing branch of stellar seismology.

We study electrical and thermal conductivities of degenerate electrons emitting and absorbing phonons in a strongly magnetized crystalline neutron star crust. We take into account modification of the phonon spectrum of a Coulomb solid of ions caused by a strong magnetic field. Boltzmann transport equation is solved using a generalized variational method. The ensuing three-dimensional integrals over the transferred momenta are evaluated by two different numerical techniques, the Monte-Carlo method and a regular integration over the first Brillouin zone. The results of the two numerical approaches are shown to be in a good agreement. An appreciable growth of electrical and thermal resistivities is reported at quantum and intermediate temperatures $T \lesssim 0.1 T_{\rm p}$ ($T_{\rm p}$ is the ion plasma temperature) in a wide range of chemical compositions and mass densities of matter even for moderately magnetized crystals $\omega_{\rm B} \sim \omega_{\rm p}$ ($\omega_{\rm B}$ and $\omega_{\rm p}$ are the ion cyclotron and plasma frequencies). This effect is due to an appearance of a soft ($\omega \propto k^2$) phonon mode in the magnetized ion Coulomb crystal, which turns out to be easier to excite than acoustic phonons characteristic of the field-free case. These results are important for modelling magneto-thermal evolution of neutron stars.

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.

Shear modulus of solid neutron star crust is calculated by thermodynamic perturbation theory taking into account ion motion. At given density the crust is modelled as a body-centered cubic Coulomb crystal of fully ionized atomic nuclei of one type with the uniform charge-compensating electron background. Classic and quantum regimes of ion motion are considered. The calculations in the classic temperature range agree well with previous Monte Carlo simulations. At these temperatures the shear modulus is given by the sum of a positive contribution due to the static lattice and a negative $\propto T$ contribution due to the ion motion. The quantum calculations are performed for the first time. The main result is that at low temperatures the contribution to the shear modulus due to the ion motion saturates at a constant value, associated with zero-point ion vibrations. Such behavior is qualitatively similar to the zero-point ion motion contribution to the crystal energy. The quantum effects may be important for lighter elements at higher densities, where the ion plasma temperature is not entirely negligible compared to the typical Coulomb ion interaction energy. The results of numerical calculations are approximated by convenient fitting formulae. They should be used for precise neutron star oscillation modelling, a rapidly developing branch of stellar seismology.

The body-centered cubic Coulomb crystal of ions in the presence of a uniform magnetic field is studied using the rigid electron background approximation. The phonon mode spectra are calculated for a wide range of magnetic field strengths and for several orientations of the field in the crystal. The phonon spectra are used to calculate the phonon contribution to the crystal energy, entropy, specific heat, Debye-Waller factor of ions, and the rms ion displacements from the lattice nodes for a broad range of densities, temperatures, chemical compositions, and magnetic fields. Strong magnetic field dramatically alters the properties of quantum crystals. The phonon specific heat increases by many orders of magnitude. The ion displacements from their equilibrium positions become strongly anisotropic. The results can be relevant for dusty plasmas, ion plasmas in Penning traps, and especially for the crust of magnetars (neutron stars with superstrong magnetic fields $B \gtrsim 10^{14}$ G). The effect of the magnetic field on ion displacements in a strongly magnetized neutron star crust can suppress the nuclear reaction rates and make them extremely sensitive to the magnetic field direction.

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 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 analyze cooling of neutron stars, assuming the presence of localized protons in the densest region of their cores. Choosing a single threshold density for proton localization and adjusting neutron star mass, we reproduce the observational data on effective surface temperatures of Vela and PSR 0656+14, with or without an accreted hydrogen envelope. However, the presence of a tiny hydrogen envelope is mandatory, in this model, for reproducing the Geminga data.

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.

As pointed out by Kutschera and W{\'o}jcik, very low concentration of protons combined with a specific density dependence of effective neutron-proton interaction could lead to a localization of ``proton impurities'' in neutron medium at densities exceeding four times normal nuclear matter density. We study consequences of the localization of protons for transport processes in dense neutron star cores, assuming random distribution of proton impurities. Kinetic equations, relevant for the transport of charge, heat and momentum, are solved using variational method. Localization of protons removes a T^{-2} factor from the transport coefficients, which leads, at lower temperatures, to a strong decrease of thermal conductivity, electrical conductivity and shear viscosity of neutron star matter, as compared to the standard case, where protons form a Fermi liquid. Due to the localization of protons a number of conventional neutrino emission processes (including modified URCA process) become inoperative in neutron star cores. On the other hand, the energy loss rate from neutrino-antineutrino pair bremsstrahlung due to electron and neutron scattering off (localized) protons, will have a specific T^6 dependence, which could modify the cooling of the neutron star core, as compared to the standard case. Possible astrophysical implications of the localization of protons for neutron star evolution and dynamics are discussed.

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.

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.

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.

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.