Mutual entrainment effects in cold neutron-proton mixtures are studied in the framework of the self-consistent nuclear energy-density functional theory. Exact expressions for the mass currents, valid for both homogeneous and inhomogeneous systems, are directly derived from the time-dependent Hartree-Fock equations with no further approximation. The equivalence with the Fermi-liquid expression is also demonstrated. Focusing on neutron-star cores, a convenient and simple analytical formulation of the entrainment matrix in terms of the isovector effective mass is found, thus allowing one to relate entrainment phenomena in neutron stars to isovector giant dipole resonances in finite nuclei. Results obtained with different functionals are presented. These include the Brussels-Montreal functionals, for which unified equations of state of neutron stars have been recently calculated.

Mutual entrainment effects in hot neutron-proton superfluid mixtures are studied in the framework of the self-consistent nuclear energy-density functional theory. The local mass currents in homogeneous or inhomogeneous nuclear systems, which we derive from the time-dependent Hartree-Fock-Bogoliubov equations at finite temperatures, are shown to have the same formal expression as the ones we found earlier in the absence of pairing at zero temperature. Analytical expressions for the entrainment matrix are obtained for application to superfluid neutron-star cores. Results are compared to those obtained earlier using Landau's theory. Our formulas, valid for arbitrary temperatures and currents, are applicable to various types of functionals including the Brussels-Montreal series for which unified equations of state have been already calculated, thus, laying the ground for a fully consistent microscopic description of superfluid neutron stars.

Temperature and velocity-dependent 1S0 pairing gaps, chemical potentials and entrainment matrix in dense homogeneous neutron–proton superfluid mixtures constituting the outer core of neutron stars, are determined fully self-consistently by solving numerically the time-dependent Hartree–Fock–Bogoliubov equations over the whole range of temperatures and flow velocities for which superfluidity can exist. Calculations have been made for npeμ in beta-equilibrium using the Brussels–Montreal functional BSk24. The accuracy of various approximations is assessed and the physical meaning of the different velocities and momentum densities appearing in the theory is clarified. Together with the unified equation of state published earlier, the present results provide consistent microscopic inputs for modeling superfluid neutron-star cores.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

The interior of mature neutron stars is expected to contain superfluid neutrons and superconducting protons. The influence of temperature and currents on superfluid properties is studied within the self-consistent timedependent nuclear energy-density functional theory. We find that this theory predicts the existence of a regime in which nucleons are superfluid (the order parameter remains finite) even though the energy spectrum of quasiparticle excitations exhibits no gap. We show that the disappearance of the gap leads to a specific heat that is not exponentially suppressed at low temperatures as in the Bardeen-Cooper-Schrieffer regime but can be comparable to that in the normal phase. Introducing some dimensionless effective superfluid velocity, we show that the behavior of the specific heat is essentially universal and we derive general approximate analytical formulas for applications to neutron-star cooling simulations.