The observed thermal relaxation of transiently accreting neutron stars during quiescence periods in low-mass X-ray binaries suggests the existence of unknown heat sources in the shallow layers of neutron-star crust. Making use of existing experimental nuclear data, we estimate the maximum possible amount of heat that can be deposited in the outer crust of an accreting neutron star due to electron captures and pycnonuclear fusion reactions triggered by the burial of X-ray burst ashes.

Mutual entrainment effects in hot neutron-proton superfluid mixtures are studied in the framework of the self-consistent nuclear energy-density functional theory. The mass currents, derived from the time-dependent Hartree-Fock-Bogoliubov equations, are shown to have the same formal expression as in the absence of pairing at zero temperature. Simple analytical expressions for the entrainment matrix are obtained for application to superfluid neutron-star cores. Results are compared to Fermi-liquid predictions.

The iterative method recently proposed for determining the internal constitution of the outer crust of a nonaccreted neutron star is extended to magnetars by taking into account the Landau-Rabi quantization of electron motion induced by the presence of a very high magnetic field. It is shown that in the strongly quantizing regime, the method can be efficiently implemented using new analytical solutions for the transitions between adjacent crustal layers. Detailed numerical computations are performed to assess the performance and precision of the method.

A very fast iterative method is presented to calculate the internal constitution of the outer crust of a cold nonaccreted neutron star, making use of very accurate analytical formulas for the transition pressures between adjacent crustal layers and their density. In addition to the composition of the different crustal layers, their depth and their baryonic mass content can be simultaneously estimated using an approximate solution of Einstein's general relativistic equations. The overall computing time is drastically reduced compared to the traditional approach, thus opening the door to large-scale statistical studies and sensitivity analyses.

The role of the crust on the tidal deformability of a cold nonaccreted neutron star is studied using the recent unified equation of state BSk24. This equation of state, which is based on the nuclear-energy density functional theory, provides a thermodynamically consistent description of all stellar regions. Results obtained with this equation of state are compared to those calculated for a putative neutron star made entirely of homogeneous matter. The presence of the crustal layers is thus found to significantly reduce the Love number $k_2$, especially for low-mass stars. However, this reduction mainly arises from the increase in the stellar radius almost independently of the equation of state. This allows for a simple analytic estimate of $k_2$ for realistic neutron stars using the equation of state of homogeneous matter only.

Nonequilibrium conditions imposed by neutrino cooling through the liquid-solid transition lead to disorder in the solid crust of neutron stars. Disorder reduces the superfluid fraction, $\rho_s/\rho$, at densities above that of neutron drip, $\rho_d \approx 4\times 10^{11}\,g/cm^3$. For an amorphous solid crust the suppression of $\rho_s$ is small, except in the highest density regions of the crust. In contrast to the strong reduction in neutron conduction predicted for coherent Bragg scattering in a crystalline crust, the disordered solid crust supports sufficient neutron superfluid density to account for pulsar glitches.

Timing of the Crab and Vela pulsars have recently revealed very peculiar evolutions of their spin frequency during the early stage of a glitch. We show that these differences can be interpreted from the interactions between neutron superfluid vortices and proton fluxoids in the core of these neutron stars. In particular, pinning of individual vortices to fluxoids is found to have a dramatic impact on the mutual friction between the neutron superfluid and the rest of the star. The number of fluxoids attached to vortices turns out to be a key parameter governing the global dynamics of the star. These results may have implications for the interpretation of other astrophysical phenomena such as pulsar-free precession or the r-mode instability.

The superfluid and superconducting core of a cold rotating neutron star is expected to be threaded by a tremendous number of neutron quantised vortices and proton fluxoids. Their interactions are unavoidable and may have important astrophysical implications. In this paper, the various contributions to the force acting on a single vortex to which fluxoids are pinned are clarified. The general expression of the force is derived by applying the variational multifluid formalism developed by Carter and collaborators. Pinning to fluxoids leads to an additional Magnus type force due to proton circulation around the vortex. Pinning in the core of a neutron star may thus have a dramatic impact on the vortex dynamics, and therefore on the magneto-rotational evolution of the star.

We generalize our earlier work on neutron stars (arXiv:1903.04981), which assumed spherical Wigner-Seitz cells in the inner crust, to admit the possibility of pasta phases, i.e., non-spherical cell shapes. Full fourth-order extended Thomas-Fermi calculations using the density functional BSk24 are performed for cylindrical and plate-like cells. Unlike in our spherical-cell calculations, we do not include shell and pairing corrections, but there are grounds for expecting these corrections for pasta to be significantly smaller. It is therefore meaningful to compare the ETF pasta results with the full spherical-cell results, i.e., with shell and pairing corrections included. However, in view of the many previous studies in which shell and pairing corrections were omitted entirely, it is of interest to compare our pasta results with the ETF part of the corresponding spherical calculations. Making this latter comparison we find that as the density increases the cell shapes pass through the usual sequence sphere $\to$ cylinder $\to$ plate before the transition to the homogeneous core. The filling fractions found at the phase transitions are in close agreement with expectations based on the liquid-drop model. On the other hand, when we compare with the full spherical-cell results, we find the sequence to be sphere $\to$ cylinder $\to$ sphere $\to$ cylinder $\to$ plate. In neither case do any "inverted", i.e., bubble-like, configurations appear. We provide accurate fitting formulas to all our essential numerical results for each of the three phases, designed especially for the density range where the nonspherical shapes are expected, which enable one to capture not only the general behavior of the fitted functions, but also the differences between them in different phases.

The interior of a neutron star (NS) is usually assumed to be made of cold catalyzed matter. However, the outer layers are unlikely to remain in full equilibrium during the formation of the star and its cooling, especially after crystallization. We study the cooling and equilibrium composition of the outer layers of a NS down to crystallization. Here the impurity parameter, usually a free parameter in cooling simulations, is calculated self-consistently using a microscopic model for which a unified equation of state has recently been determined. We follow the evolution of the nuclear distributions of the multi-component Coulomb liquid plasma (MCP) fully self-consistently, adapting a general formalism originally developed for the description of supernova cores. We calculate the impurity parameter at the crystallization as determined in the one-component plasma (OCP) approximation. Our analysis shows that the sharp changes in composition obtained in the OCP approximation are smoothed out when a full nuclear distribution is allowed. The Coulomb coupling parameter at melting is found to be reasonably close to the canonical value of 175, except for specific pressures for which supercooling occurs in the OCP approximation. Our MCP treatment leads to non-monotonic variations of the impurity parameter with pressure. Its values can change by several orders of magnitude reaching about 50, suggesting that the crust may be composed of an alternation of pure (highly conductive) and impure (highly resistive) layers. The results presented here complement the recent unified equation of state obtained with the same model. Our self-consistent approach to hot MCP shows that the presence of impurities in the outer crust of a NS is non-negligible and may have a sizeable impact on transport properties. In turn, this may have important implications for the cooling of NS and their magneto-rotational evolution.

Context. In the cooling process of a non-accreting neutron star, the composition and properties of the crust are thought to be fixed at the finite temperature where nuclear reactions fall out of equilibrium. A lower estimation for this temperature is given by the crystallization temperature, which can be as high as $\approx 7\times 10^9$ K in the inner crust, potentially leading to sizeable differences with respect to the simplifying cold-catalyzed matter hypothesis. Aims. We extend the recent work by Fantina et al. (2019) on the outer crust, to the study of the crystallization of the inner crust and the associated composition in the one-component plasma approximation. Methods. The finite temperature variational equations for non-uniform matter in both the liquid and the solid phases are solved using a compressible liquid-drop approach with parameters optimized on four different microscopic models which cover the present uncertainties in nuclear modeling. Results. We consider separately the effect of the different nuclear ingredients with their associated uncertainties, namely the nuclear equation of state, the surface properties in the presence of a uniform gas of dripped neutrons, and the proton shell effects arising from the ion single-particle structure. Our results suggest that the highest source of model dependence comes from the smooth part of the nuclear functional. Conclusions. We show that shell effects play an important role at the lowest densities close to the outer crust, but the most important physical ingredient to be settled for a quantitative prediction of the inner crust properties is the surface tension at extreme isospin values.

The role of the symmetry energy and the neutron-matter stiffness on the tidal deformability of a cold nonaccreted neutron star is studied using a set of unified equations of state. Based on the nuclear energy-density functional theory, these equations of state provide a thermodynamically consistent treatment of all regions of the star and were calculated using functionals that were precision fitted to experimental and theoretical nuclear data. Predictions are compared to constraints inferred from the recent detection of the gravitational-wave signal GW170817 from a binary neutron-star merger and from observations of the electromagnetic counterparts.

The role of the symmetry energy on the internal constitution and the global structure of a cold nonaccreted neutron star is studied using a set of unified equations of state. Based on the nuclear energy-density functional theory, these equations of state provide a thermodynamically consistent treatment of all regions of the star and were calculated using the four different Brussels-Montreal functionals BSk22, BSk24, BSk25 and BSk26. Our predictions are compared to various constraints inferred from astrophysical observations including the recent detection of the gravitational wave signal GW170817 from a binary neutron-star merger.

Magnetic fields of order $10^{15}$ G have been measured at the surface of some neutron stars, and much stronger magnetic fields are expected to be present in the solid region beneath the surface. The effects of the magnetic field on the equation of state and on the composition of the crust due to Landau-Rabi quantization of electron motion are studied. Both the outer and inner crustal regions are described in a unified and consistent way within the nuclear-energy density functional theory.

The theory of the nuclear energy-density functional is used to provide a unified and thermodynamically consistent treatment of all regions of cold non-accreting neutron stars. In order to assess the impact of our lack of complete knowledge of the density dependence of the symmetry energy on the constitution and the global structure of neutron stars, we employ four different functionals. All of them were precision fitted to essentially all the nuclear-mass data with the Hartree-Fock-Bogoliubov method and two different neutron-matter equations of state based on realistic nuclear forces. For each functional, we calculate the composition, the pressure-density relation, and the chemical potentials throughout the star. We show that uncertainties in the symmetry energy can significantly affect the theoretical results for the composition and global structure of neutron stars. To facilitate astrophysical applications, we construct analytic fits to our numerical results.

Formed in the aftermath of gravitational core-collapse supernova explosions, neutron stars are unique cosmic laboratories for probing the properties of matter under extreme conditions that cannot be reproduced in terrestrial laboratories. The interior of a neutron star, endowed with the highest magnetic fields known and with densities spanning about ten orders of magnitude from the surface to the centre, is predicted to exhibit various phases of dense strongly interacting matter, whose physics is reviewed in this chapter. The outer layers of a neutron star consist of a solid nuclear crust, permeated by a neutron ocean in its densest region, possibly on top of a nuclear "pasta" mantle. The properties of these layers and of the homogeneous isospin asymmetric nuclear matter beneath constituting the outer core may still be constrained by terrestrial experiments. The inner core of highly degenerate, strongly interacting matter poses a few puzzles and questions which are reviewed here together with perspectives for their resolution. Consequences of the dense-matter phases for observables such the neutron-star mass-radius relationship and the prospects to uncover their structure with modern observational programmes are touched upon.

Observations of soft X-ray transients in quiescence suggest the existence of heat sources in the crust of accreting neutron stars. Heat is thought to be released by electroweak and nuclear processes triggered by the burying of ashes of X-ray bursts. The heating is studied using a fully quantum approach taking consistently into account nuclear shell effects. We have followed the evolution of ashes made of $^{56}$Fe employing the nuclear energy-density functional theory. Both the outer and inner crusts are described using the same functional, thus ensuring a unified and thermodynamically consistent treatment. To assess the role of the neutron-matter constraint, we have employed the set of accurately calibrated Brussels-Montreal functionals BSk19, BSk20, and BSk21 and for comparison the SLy4 functional. Due to nuclear shell effects, the fully accreted crust is found to be much less stratified than in previous studies. In particular, large regions of the inner crust contain clusters with the magic number $Z=14$. The heat deposited in the outer crust is tightly constrained by experimental atomic mass data. The shallow heating we obtain does not exceed $0.2$~MeV and is therefore not enough to explain the cooling of some soft X-ray transients. The total heat released in the crust is very sensitive to details of the nuclear structure and is predicted to lie in the range from $1.5$~MeV to $1.7$~MeV. The evolution of an accreted matter element and therefore the location of heat sources are governed to a large extent by the existence of nuclear shell closures. Ignoring these effects in the inner crust, the total heat falls to $\sim 0.6$~MeV. The neutron-matter constraint is also found to play a key role. The large amount of heat obtained by Steiner et al. (2012) could thus be traced back to unrealistic neutron-matter equations of state.

Neutron stars, the compact stellar remnants of core-collapse supernova explosions, are unique cosmic laboratories for exploring novel phases of matter under extreme conditions. In particular, the occurrence of superfluidity and superconductivity in neutron stars will be briefly reviewed.

In spite of the absence of viscous drag, the neutron superfluid permeating the inner crust of a neutron star cannot flow freely, and is entrained by the nuclear lattice similarly to laboratory superfluid atomic gases in optical lattices. The role of entrainment on the neutron superfluid dynamics is reviewed. For this purpose, a minimal hydrodynamical model of superfluidity in neutron-star crusts is presented. This model relies on a fully four-dimensionally covariant action principle. The equivalence of this formulation with the more traditional approach is demonstrated. In addition, the different treatments of entrainment in terms of dynamical effective masses or superfluid density are clarified. The nuclear energy density functional theory employed for the calculations of all the necessary microscopic inputs is also reviewed, focusing on superfluid properties. In particular, the microscopic origin of entrainment and the different methods to estimate its importance are discussed.

In this paper, we study in detail the role of general relativity on the global dynamics of giant pulsar glitches as exemplified by Vela. For this purpose, we carry out numerical simulations of the spin up triggered by the sudden unpinning of superfluid vortices. In particular, we compute the exchange of angular momentum between the core neutron superfluid and the rest of the star within a two-fluid model including both (non-dissipative) entrainment effects and (dissipative) mutual friction forces. Our simulations are based on a quasi-stationary approach using realistic equations of state (EoSs). We show that the evolution of the angular velocities of both fluids can be accurately described by an exponential law. The associated characteristic rise time $\tau_{\text{r}}$, which can be precisely computed from stationary configurations only, has a form similar to that obtained in the Newtonian limit. However, general relativity changes the structure of the star and leads to additional couplings between the fluids due to frame-dragging effects. As a consequence, general relativity can have a large impact on the actual value of $\tau_{\text{r}}$: the errors incurred by using Newtonian gravity are thus found to be as large as $\sim 40 \%$ for the models considered. Values of the rise time are calculated for Vela and compared with current observational limits. Finally, we study the amount of gravitational waves emitted during a glitch. Simple expressions are obtained for the corresponding characteristic amplitudes and frequencies. The detectability of glitches through gravitational wave observatories is briefly discussed.

We develop a detailed and self-consistent numerical model for extremely-magnetised white dwarfs, which have been proposed as progenitors of overluminous Type Ia supernovae. This model can describe fully-consistent equilibria of magnetic stars in axial symmetry, with rotation, general-relativistic effects and realistic equations of state (including electron-ion interactions and taking into account Landau quantisation of electrons due to the magnetic field). We study the influence of each of these ingredients onto the white dwarf structure and, in particular, on their maximum mass. We perform an extensive stability analysis of such objects, with their highest surface magnetic fields reaching $\sim 10^{13}~G$ (at which point the star adopts a torus-like shape). We confirm previous speculations that although very massive strongly magnetised white dwarfs could potentially exist, the onset of electron captures and pycnonuclear reactions may severely limit their stability. Finally, the emission of gravitational waves by these objects is addressed, showing no possibility of detection by the currently planned space-based detector eLISA.

Magnetars are neutron stars endowed with surface magnetic fields of the order of $10^{14}-10^{15}$~G, and with presumably much stronger fields in their interior. As a result of Landau quantization of electron motion, the neutron-drip transition in the crust of a magnetar is shifted to either higher or lower densities depending on the magnetic field strength. The impact of nuclear uncertainties is explored considering the recent series of Brussels-Montreal microscopic nuclear mass models. All these models are based on the Hartree-Fock-Bogoliubov method with generalized Skyrme functionals. They differ in their predictions for the symmetry energy coefficient at saturation, and for the stiffness of the neutron-matter equation of state. For comparison, we have also considered the very accurate but more phenomenological model of Duflo and Zuker. Although the equilibrium composition of the crust of a magnetar and the onset of neutron emission are found to be model dependent, the quantum oscillations of the threshold density are essentially universal.

Pulsars are spinning extremely rapidly with periods as short as about $1.4$ milliseconds and delays of a few milliseconds per year at most, thus providing the most accurate clocks in the Universe. Nevertheless, sudden spin ups have been detected in some pulsars like the emblematic Vela pulsar. These abrupt changes in the pulsar's rotation period have long been thought to be the manifestation of a neutron superfluid permeating the inner crust of neutron stars. However, the neutron superfluid has been recently found to be so strongly coupled to the crust that it does not carry enough angular momentum to explain the Vela data. We explore the extent to which pulsar-timing observations can be reconciled with the standard glitch theory considering the lack of knowledge of the dense-matter equation of state.

Giant pulsar frequency glitches as detected in the emblematic Vela pulsar have long been thought to be the manifestation of a neutron superfluid permeating the inner crust of a neutron star. However, this superfluid has been recently found to be entrained by the crust, and as a consequence it does not carry enough angular momentum to explain giant glitches. The extent to which pulsar-timing observations can be reconciled with the standard vortex-mediated glitch theory is studied considering the current uncertainties on dense-matter properties. To this end, the crustal moment of inertia of glitching pulsars is calculated employing a series of different unified dense-matter equations of state.

Context. A unified equation of state (EoS) should describe the crust and the core of a neutron star using the same physical model. The Brussels-Montreal group has recently derived a family of such EoSs based on the nuclear energy-density functional theory with generalized Skyrme effective forces, fitted to the available mass data. At the same time, these forces were constrained to reproduce microscopic calculations of homogeneous neutron matter based on realistic two- and three-nucleon forces. Aims. We represent basic physical characteristics of the latest Brussels-Montreal EoS models by analytical expressions to facilitate their inclusion in astrophysical simulations. Methods. We consider three EoS models, which significantly differ by stiffness: BSk19, BSk20, and BSk21. For each of them we constructed two versions of the EoS parametrization. In the first version, pressure P and gravitational mass density \rho are given as functions of the baryon number density n_b. In the second version, P, \rho, and n_b are given as functions of pseudo-enthalpy, which is useful for two-dimensional calculations of stationary rotating configurations of neutron stars. In addition to the EoS, we derived analytical expressions for several related quantities that are required in neutron-star simulations: number fractions of electrons and muons in the stellar core, nucleon numbers per nucleus in the inner crust, and equivalent radii and shape parameters of the nuclei in the inner crust. Results. We obtain analytical representations for the basic characteristics of the models of cold dense matter, which are most important for studies of neutron stars. We demonstrate the usability of our results by applying them to calculations of neutron-star mass-radius relations, maximum and minimum masses, thresholds of direct Urca processes, and the electron conductivity in the neutron-star crust.

The neutron-drip transition in the dense matter constituting the interior of neutron stars generally refers to the appearance of unbound neutrons as the matter density reaches some threshold density $\rho_\textrm{drip}$. This transition has been mainly studied under the cold catalyzed matter hypothesis. However, this assumption is unrealistic for accreting neutron stars. After examining the physical processes that are thought to be allowed in both accreting and nonaccreting neutron stars, suitable conditions for the onset of neutron drip are derived and general analytical expressions for the neutron drip density and pressure are obtained. Moreover, we show that the neutron-drip transition occurs at lower density and pressure than those predicted within the mean-nucleus approximation. This transition is studied numerically for various initial composition of the ashes from X-ray bursts and superbursts using microscopic nuclear mass models.

The heat capacity of neutron matter is studied over the range of densities and temperatures prevailing in neutron-star crusts, allowing for the transition to a superfluid phase at temperatures below some critical temperature $T_{sf}$ and including the transition to the classical limit. Finite temperature Hartree-Fock-Bogoliubov equations (FTHFB) are solved and compared to existing approximate expressions. In particular, the formula given by Levenfish and Yakovlev is found to reproduce the numerical results with a high degree of accuracy for temperatures $T\leq T_{sf}$. In the non-superfluid phase, $T\geq T_{sf}$, the linear approximation is valid only at temperature $T\ll T_{{\rm F} n}$ ($T_{{\rm F} n}$ being the Fermi temperature of the neutron gas) which is rarely the case in the shallow layers of the neutron star's crust. A non-perturbative interpolation between the quantal and the classical regimes is proposed here. The heat capacity, conveniently parametrized solely in terms of $T_{sf}$, $T_{{\rm F} n}$, and the neutron number density $n_n$, can be easily implemented in neutron-star cooling simulations.

One of the most intringuing questions about neutron stars concerns their maximum mass. The answer is intimately related to the properties of matter at densities far beyond that found in heavy atomic nuclei. The current view on the internal constitution of neutron stars and on their maximum mass, both from theoretical and observational studies, are briefly reviewed.

We explore the spectrum of low-energy collective excitations in the crust of a neutron star, especially in the inner region where neutron-proton clusters are immersed in a sea of superfluid neutrons. The speeds of the different modes are calculated systematically from the nuclear energy density functional theory using a Skyrme functional fitted to essentially all experimental atomic mass data.

It has been recently proposed that very massive white dwarfs endowed with strongly quantizing magnetic fields might be the progenitors of overluminous type Ia supernovae like SN 2006gz and SN 2009dc. In this work, we show that the onset of electron captures and pycnonuclear reactions in these putative super-Chandrasekhar white dwarfs may severely limit their stability.

We describe the main features of our most recent Hartree-Fock-Bogoliubov nuclear mass models, based on 16-parameter generalized Skyrme forces. They have been fitted to the data of the 2012 Atomic Mass Evaluation, and favour a value of 30 MeV for the symmetry coefficient J, the corresponding root-mean square deviation being 0.549 MeV. We find that this conclusion is compatible with measurements of neutron-skin thickness. By constraining the underlying interactions to fit various equations of state of neutron matter calculated {\it ab initio} our models are well adapted to a realistic and unified treatment of all regions of neutron stars. We use our models to calculate the composition, the equation of state, the mass-radius relation and the maximum mass. Comparison with observations of neutron stars again favours a value of J = 30 MeV.

The equilibrium properties of the outer crust of cold nonaccreting magnetars (i.e. neutron stars endowed with very strong magnetic fields) are studied using the latest experimental atomic mass data complemented with a microscopic atomic mass model based on the Hartree-Fock-Bogoliubov method. The Landau quantization of electron motion caused by the strong magnetic field is found to have a significant impact on the composition and the equation of state of crustal matter. It is also shown that the outer crust of magnetars could be much more massive than that of ordinary neutron stars.

We study the impact of a hadron-quark phase transition on the maximum neutron-star mass. The hadronic part of the equation of state relies on the most up-to-date Skyrme nuclear energy density functionals, fitted to essentially all experimental nuclear mass data and constrained to reproduce the properties of infinite nuclear matter as obtained from microscopic calculations using realistic forces. We show that the softening of the dense matter equation of state due to the phase transition is not necessarily incompatible with the existence of massive neutron stars like PSR J1614-2230.

We study the low-energy collective excitations of the neutron star inner crust, where a neutron superfluid coexists with a Coulomb lattice of nuclei. The dispersion relation of the modes is calculated systematically from a microscopic theory including neutron band structure effects. These effects are shown to lead to a strong mixing between the Bogoliubov-Anderson bosons of the neutron superfluid and the longitudinal crystal lattice phonons. In addition, the speed of the transverse shear mode is greatly reduced as a large fraction of superfluid neutrons are entrained by nuclei. Not only does the much smaller velocity of the transverse mode increase the specific heat of the inner crust, but it also decreases its electron thermal conductivity. These results may impact our interpretation of the thermal relaxation in accreting neutron stars. Due to strong mixing, the mean free path of the superfluid mode is found to be greatly reduced. Our results for the collective mode dispersion relations and their damping may also have implications for neutron star seismology.

Large pulsar frequency glitches are generally interpreted as sudden transfers of angular momentum between the neutron superfluid permeating the inner crust and the rest of the star. Despite the absence of viscous drag, the neutron superfluid is strongly coupled to the crust due to non-dissipative entrainment effects. These effects are shown to severely limit the maximum amount of angular momentum that can possibly be transferred during glitches. In particular, it is found that the glitches observed in the Vela pulsar require an additional reservoir of angular momentum.

The neutron superfluid permeating the inner crust of mature neutron stars is expected to play a key role in various astrophysical phenomena like pulsar glitches. Despite the absence of viscous drag, the neutron superfluid can still be coupled to the solid crust due to non-dissipative entrainment effects. Entrainment challenges the interpretation of pulsar glitches and suggests that a revision of the interpretation of other observed neutron-star phenomena might be necessary.

We have recently developed a set of equations of state based on the nuclear energy density functional theory providing a unified description of the different regions constituting the interior of neutron stars and magnetars. The nuclear functionals, which were constructed from generalized Skyrme effective nucleon-nucleon interactions, yield not only an excellent fit to essentially all experimental atomic mass data but were also constrained to reproduce the neutron-matter equation of state as obtained from realistic many-body calculations.

The equation of state (EoS) of the outer crust of a cold non-accreting magnetar has been determined using the model of Lai and Shapiro (1991). For this purpose, we have made use of the latest experimental atomic mass data complemented with a Hartree-Fock-Bogoliubov (HFB) mass model. Magnetar crusts are found to be significantly different from the crusts of ordinary neutron stars.

Despite the absence of viscous drag, the neutron superfluid permeating the inner crust of a neutron star can still be strongly coupled to nuclei due to non-dissipative entrainment effects. Neutron superfluidity and entrainment have been systematically studied in all regions of the inner crust of a cold non-accreting neutron star in the framework of the band theory of solids. It is shown that in the intermediate layers of the inner crust a large fraction of "free" neutrons are actually entrained by the crust. The results suggest that a revision of the interpretation of many observable astrophysical phenomena might be necessary.

We present a set of three unified equations of states (EoSs) based on the nuclear energy-density functional (EDF) theory.These EoSs are based on generalized Skyrme forces fitted to essentially all experimental atomic mass data and constrained to reproduce various properties of infinite nuclear matter as obtained from many-body calculations using realistic two- and three-body interactions. The structure of cold isolated neutron stars is discussed in connection with some astrophysical observations.

Nuclear pairing is studied both in atomic nuclei and in neutron-star crusts in the unified framework of the energy-density functional theory using generalized Skyrme functionals complemented with a local pairing functional obtained from many-body calculations in homogeneous nuclear matter using realistic forces.

Even though the "free" neutrons in the inner crust of a neutron star are superfluid, they are still strongly coupled to nuclei due to non-dissipative entrainment effects. These effects have been systematically studied in all regions of the inner crust in the framework of the band theory of solids. Using concepts from solid-state physics, it is shown that the density of conduction neutrons, i.e. neutrons that are effectively "free", can be much smaller than the density of unbound neutrons (by an order of magnitude in some layers) due to Bragg scattering. These results suggest that a revision of the interpretation of various observable neutron-star phenomena may be necessary.

We calculate the maximum mass of neutron stars for three different equations of state (EOS) based on generalized Skyrme functionals that are simultaneously fitted to essentially all the 2003 nuclear mass data (the rms deviation is 0.58 MeV in all three cases) and to one or other of three different equations of state of pure neutron matter, each determined by a different many-body calculation using realistic two- and three-body interactions but leading to significantly different degrees of stiffness at the high densities prevailing in neutron-star interiors. The observation of a neutron star with mass 1.97 $\pm$ 0.04 $\mathcal{M}_{\odot}$ eliminates the softest of our models (BSk19), but does not discriminate between BSk20 and BSk21. However, nuclear-mass measurements that have been made since our models were constructed strongly favor BSk21, our stiffest functional.

The rapid neutron-capture process, or r-process, is known to be fundamental for explaining the origin of approximately half of the A>60 stable nuclei observed in nature. In recent years nuclear astrophysicists have developed more and more sophisticated r-process models, by adding new astrophysical or nuclear physics ingredients to explain the solar system composition in a satisfactory way. Despite these efforts, the astrophysical site of the r-process remains unidentified. The composition of the neutron star outer crust material is investigated after the decompression that follows its possible ejection. The composition of the outer crust of a neutron star is estimated before and after decompression. Two different possible initial conditions are considered, namely an idealized crust composed of cold catalyzed matter and a crust initially in nuclear statistical equilibrium at temperatures around 10 GK. We show that in this second case before decompression and at temperatures typically corresponding to 8 GK, the Coulomb effect due to the high densities in the crust leads to an overall composition of the outer crust in neutron-rich nuclei with a mass distribution close to the solar system r-abundance distribution. Such distributions differ, however, from the solar one due to a systematic shift in the second peak to lower values. After decompression, the capture of the few neutrons per seed nucleus available in the hot outer crust leads to a final distribution of stable neutron-rich nuclei with a mass distribution of 80 < A < 140 nuclei in excellent agreement with the solar distribution, provided the outer crust is initially at temperatures around 8 GK and all layers of the outer crust are ejected. The decompression of the neutron star matter from the outer crust provides suitable conditions for a robust r-processing of the light species, i.e., r-nuclei with A < 140.

Self-consistent mean field methods based on phenomenological Skyrme effective interactions are known to exhibit spurious spin and spin-isospin instabilities both at zero and finite temperatures when applied to homogeneous nuclear matter at the densities encountered in neutron stars and in supernova cores. The origin of these instabilities is revisited in the framework of the nuclear energy density functional theory and a simple prescription is proposed to remove them. The stability of several Skyrme parametrizations is reexamined.

Mature neutron stars are expected to contain various kinds of superfluids in their interiors. Modeling such stars requires the knowledge of the mutual entrainment couplings between the different condensates. We present a unified equation of state describing the different regions of a neutron star with superfluid neutrons and superconducting protons in its core.

We construct three new Hartree-Fock-Bogoliubov (HFB) mass models, labeled HFB-19, HFB-20, and HFB-21, with unconventional Skyrme forces containing $t_4$ and $t_5$ terms, i.e., density-dependent generalizations of the usual $t_1$ and $t_2$ terms, respectively. The new forces underlying these models are fitted respectively to three different realistic equations of state of neutron matter for which the density dependence of the symmetry energy ranges from the very soft to the very stiff, reflecting thereby our present lack of complete knowledge of the high-density behavior of nuclear matter. All unphysical instabilities of nuclear matter, including the transition to a polarized state in neutron-star matter, are eliminated with the new forces. At the same time the new models fit essentially all the available mass data with rms deviations of 0.58 MeV and give the same high quality fits to measured charge radii that we obtained in earlier models with conventional Skyrme forces. Being constrained by neutron matter, these new mass models, which all give similar extrapolations out to the neutron drip line, are highly appropriate for studies of the $r$-process and the outer crust of neutron stars. Moreover, the underlying forces, labeled BSk19, BSk20 and BSk21, respectively, are well adapted to the study of the inner crust and core of neutron stars. The new family of Skyrme forces thus opens the way to a unified description of all regions of neutron stars.

A simple description of superfluid hydrodynamics in the inner crust of a neutron star is given. Particular attention is paid to the effect of the lattice of nuclei on the properties of the superfluid neutrons, and the effects of entrainment, the fact that some fraction of the neutrons are locked to the motion of the protons in nuclei.

Non-empirical effective contact pairing forces to be used in self-consistent mean-field calculations are presented. These pairing forces, constructed so as to reproduce exactly any given microscopic pairing gaps in infinite homogeneous nuclear matter for any isospin asymmetry, are given in analytical form. As a by-product, this work provides an analytical solution of the BCS gap equations which could be applied to describe various many-body systems.

The neutron superfluidity in the inner crust of a neutron star has been traditionally studied considering either homogeneous neutron matter or only a small number of nucleons confined inside the spherical Wigner-Seitz cell. Drawing analogies with the recently discovered multi-band superconductors, we have solved the anisotropic multi-band BCS gap equations with Bloch boundary conditions, thus providing a unified description taking consistently into account both the free neutrons and the nuclear clusters. Calculations have been carried out using the effective interaction underlying our recent Hartree-Fock-Bogoliubov nuclear mass model HFB-16. We have found that even though the presence of inhomogeneities lowers the neutron pairing gaps, the reduction is much less than that predicted by previous calculations using the Wigner-Seitz approximation. We have studied the disappearance of superfluidity with increasing temperature. As an application we have calculated the neutron specific heat, which is an important ingredient for modeling the thermal evolution of newly-born neutron stars. This work provides a new scheme for realistic calculations of superfluidity in neutron-star crusts.

We construct a new Hartree-Fock-Bogoliubov (HFB) mass model, labeled HFB-18, with a generalized Skyrme force. The additional terms that we have introduced into the force are density-dependent generalizations of the usual $t_1$ and $t_2$ terms, and are chosen in such a way as to avoid the high-density ferromagnetic instability of neutron stars that is a general feature of conventional Skyrme forces, and in particular of the Skyrme forces underlying all the HFB mass models that we have developed in the past. The remaining parameters of the model are then fitted to essentially all the available mass data, an rms deviation $\sigma$ of 0.585 MeV being obtained. The new model thus gives almost as good a mass fit as our best-fit model HFB-17 ($\sigma$ = 0.581 MeV), and has the advantage of avoiding the ferromagnetic collapse of neutron stars.

The inner crust of neutron stars, formed of a crystal lattice of uclear clusters immersed in a sea of unbound neutrons, may be the nique example of periodic nuclear systems. We have calculated the neutron specific heat in the shallow part of the crust using the band theory of solids with Skyrme nucleon-nucleon interactions. We have also tested the validity of various approximations. We have found that the neutron specific heat is well described by that of a Fermi gas, while the motion of the unbound neutrons is strongly affected by the nuclear lattice. These apparently contradictory results are explained by the particular properties of the neutron Fermi surface.

The physics of neutron star crusts is vast, involving many different research fields, from nuclear and condensed matter physics to general relativity. This review summarizes the progress, which has been achieved over the last few years, in modeling neutron star crusts, both at the microscopic and macroscopic levels. The confrontation of these theoretical models with observations is also briefly discussed.

An approach to the equation of state for the inner crust of neutron stars based on Skyrme-type forces is presented. Working within the Wigner-Seitz picture, the energy is calculated by the TETF (temperature-dependent extended Thomas-Fermi) method, with proton shell corrections added self-consistently by the Strutinsky-integral method. Using a Skyrme force that has been fitted to both neutron matter and to essentially all the nuclear mass data, we find strong proton shell effects: proton numbers $Z$ = 50, 40 and 20 are the only values possible in the inner crust, assuming that nuclear equilibrium is maintained in the cooling neutron star right down to the ambient temperature. Convergence problems with the TETF expansion for the entropy, and our way of handling them, are discussed. Full TETF expressions for the specific heat of inhomogeneous nuclear matter are presented. Our treatment of the electron gas, including its specific heat, is essentially exact, and is described in detail.

Both relativistic and non-relativistic two-fluid models of neutron star cores are constructed, using the constrained variational formalism developed by Brandon Carter and co-workers. We consider a mixture of superfluid neutrons and superconducting protons at zero temperature, taking into account mutual entrainment effects. Leptons, which affect the interior composition of the neutron star and contribute to the pressure, are also included. We provide the analytic expression of the Lagrangian density of the system, the so-called master function, from which the dynamical equations can be obtained. All the microscopic parameters of the models are calculated consistently using the non-relativistic nuclear energy density functional theory. For comparison, we have also considered relativistic mean field models. The correspondence between relativistic and non-relativistic hydrodynamical models is discussed in the framework of the recently developed 4D covariant formalism of Newtonian multi-fluid hydrodynamics. We have shown that entrainment effects can be interpreted in terms of dynamical effective masses that are larger in the relativistic case than in the Newtonian case. With the nuclear models considered in this work, we have found that the neutron relativistic effective mass is even greater than the bare neutron mass in the liquid core of neutron stars.

For more than three decades, the inner crust of neutron stars, formed of a solid lattice of nuclear clusters coexisting with a gas of electrons and neutrons, has been traditionally studied in the Wigner-Seitz approximation. The validity of this approximation is discussed in the general framework of the band theory of solids, which has been recently applied to the nuclear context. Using this novel approach, it is shown that the unbound neutrons move in the crust as if their mass was increased.

The inner crust of neutron stars formed of nuclear clusters immersed in a neutron sea has been widely studied in the framework of the Wigner-Seitz approximation since the seminal work of Negele and Vautherin. In this article, the validity of this approximation is discussed in the framework of the band theory of solids. For a typical cell of $^{200}$Zr, present in the external layers of the inner crust, it is shown that the ground state properties of the neutron gas are rather well reproduced by the Wigner-Seitz approximation, while its dynamical properties depend on the energy scale of the process of interest or on the temperature. It is concluded that the Wigner-Seitz approximation is well suited for describing the inner crust of young neutron stars and the collapsing core of massive stars during supernovae explosions. However the band theory is required for low temperature fluid dynamics.

The structure of the crust of a neutron star is completely determined by the experimentally measured nuclear masses up to a density of the order of 10^11 g.cm^{-3}. At higher densities, the composition of the crust still remains uncertain, mainly due to the presence of ``free'' superfluid neutrons which affect the properties of the nuclear ``clusters''. After briefly reviewing calculations of the equilibrium structure of the crust, we point out that the current approach based on the Wigner-Seitz approximation does not properly describe the unbound neutrons. We have recently abandoned this approximation by applying the band theory of solids. We have shown that the dynamical properties of the free neutrons are strongly affected by the clusters by performing 3D calculations with Bloch boundary conditions.

Hydrodynamical simulations of neutron star cores, based on a two fluid description in terms of a neutron-proton superfluid mixture, require the knowledge of the Andreev-Bashkin entrainment matrix which relates the momentum of one constituent to the currents of both constituents. This matrix is derived for arbitrary nuclear asymmetry at zero temperature and in the limits of small relative currents in the framework of the energy density functional theory. The Skyrme energy density functional is considered as a particular case. General analytic formulae for the entrainment parameters and various corresponding effective masses are obtained. These formulae are applied to the liquid core of a neutron star, composed of an homogeneous plasma of nucleons, electrons and possibly muons in beta equilibrium.

The build up of the stress whose relaxation is presumed to account for pulsar frequency glitches can be attributed to various mechanisms, of which the most efficient involve differential rotation of the neutron superfluid in the inner layers of the (magnetically braked) solid crust of a rotating neutron star. In such a case it is usually supposed that the stress is attributable to pinning of superfluid vortices to crust nuclei, but it has been suggested that, even if the pinning effect is too weak, a comparably large stress might still arise just from the deficit of centrifugal buoyancy in the slowed down crust. The present work is a re-examination that investigates the way such processes may be affected by considerations that were overlooked in the previous work -- notably uncertainties about the ``effective'' masses that have to be attributed to the ``free'' superfluid neutrons to allow for their entrainment by the ionic crust material. Though restricted to a Newtonian formulation, this analysis distinguishes more carefully than has been usual between true velocities, which are contravariantly vectorial, and so called ``superfluid velocities'' that are proportional to momenta, which are essentially covectorial, a technicality that is important when more than one independent current is involved. The results include a Proudman type theorem to the effect that the superfluid angular velocity must be constant on slightly deformed Taylor cylinders in the force free case, and it is shown how to construct a pair of integral constants of the motion that determine the solution for the pinned case assuming beta equilibrium.

Phenomena such as pulsar frequency glitches are believed to be attributable to differential rotation of a current of ``free'' superfluid neutrons at densities above the ``drip'' threshold in the ionic crust of a neutron star. Such relative flow is shown to be locally describable by adaption of a canonical two fluid treatment that emphasizes the role of the momentum covectors constructed by differentiation of action with respect to the currents, with allowance for stratification whereby the ionic number current may be conserved even when the ionic charge number Z is altered by beta processes. It is demonstrated that the gauge freedom to make different choices of the chemical basis determining which neutrons are counted as ``free'' does not affect their ``superfluid'' momentum covector, which must locally have the form of a gradient (though it does affect the ``normal'' momentum covector characterising the protons and those neutrons that are considered to be ``confined'' in the nuclei). It is shown how the effect of ``entrainment'' (whereby the momentum directions deviate from those of the currents) is controlled by the (gauge independent) mobility coefficient K, estimated in recent microscopical quantum mechanical investigations, which suggest that the corresponding (gauge dependent) ``effective mass'' m* of the free neutrons can become very large in some layers. The relation between this treatment of the crust layers and related work (using different definitions of ``effective mass'') intended for the deeper core layers is discussed.

The inner layers of a neutron star crust, composed of a Coulomb lattice of neutron rich nuclear clusters immersed in a sea of ``free'' superfluid neutrons, are closely analogous to periodic condensed matter systems such as electronic, photonic or phononic crystals. Applying methods from solid state physics to the neutron star context, we study the transport properties of those ``free'' neutrons for the outermost layers of the inner crust, near the drip point $\rho_{\rm drip} \sim 4\times 10^{11}$ g.cm$^{-3}$. In particular, we evaluate the effective neutron mass resulting from Bragg scattering by a band structure calculation. Comparison is made with the case of electrons in solids. The observational consequences are briefly discussed.

The relative current density $n^i$ of conduction neutrons in a neutron star crust beyond the neutron drip threshold can be expected to be related to the corresponding particle momentum covector $p_i$ by a linear relation of the form $n^i=K^{ij}p_j$ in terms of a physically well defined mobility tensor $K^{ij}$. This result is describable as an ``entrainment'' whose effect - wherever the crust lattice is isotropic - will simply be to change the ordinary neutron mass m to an effective mass $m_\star$ such that in terms of the relevant number density n of unconfined neutrons we shall have $K^{ij}=n/m_\star\gamma^{ij}$. In a preceding analysis based on an independent particle treatment using Bloch boundary conditions to obtain the distribution of energy $E_k$ and associated group velocity $v_k^i=\partial E_k/\partial\hbar k_i$ as a function of wavenumber $k_i$, it was shown that the mobility tensor would be given by $K^{ij}\propto\int d^3 k v_k^i v_k^j\delta\{E_k-\mu\}$, where $\mu$ is the Fermi energy. Using the approach due to Bogoliubov, it is shown here that the effect of BCS pairing with an energy gap \Delta_F and corresponding quasiparticle energy function \EE_k=\sqrt{(E_k-\mu)^2+\Delta_F^2}$ will just be to replace the Dirac distributional integrand by a smoother distribution : $K^{ij}\propto \int d^3 k v_k^i v_k^j \Delta_F^2 /\EE^3_k$. It is also shown how the pairing condensation gives rise to superfluidity in the technical sense of providing (meta) stability against resistive perturbations for a current that is not too strong (its momentum $p_i$ must be small enough to give $2|p_i v_k^i|<\EE_k^2/|E_k-\mu|$ for all modes). It is concluded that the prediction of a very large effective mass enhancement in the middle layers of the crust will not be significantly affected by the pairing mechanism.

This work develops the dynamics of perfectly elastic solid model for application to the outer crust of a magnetised neutron star. Particular attention is given to the Noether identities responsible for energy-momentum conservation, using a formulation that is fully covariant, not only (as is usual) in a fully relativistic treatment but also (sacrificing accuracy and elegance for economy of degrees of gravitational freedom) in the technically more complicated case of the Newtonian limit. The results are used to obtain explicit (relativistic and Newtonian) formulae for the propagation speeds of generalised (Alfven type) magneto-elastic perturbation modes.

In the inner crust of a neutron star, at densities above the ``drip'' threshold, unbound ``conduction'' neutrons can move freely past through the ionic lattice formed by the nuclei. The relative current density $n^i= n \bar v^i$ of such conduction neutrons will be related to the corresponding mean particle momentum $p_i$ by a proportionality relation of the form $n^i= {\cal K}p^i$ in terms of a physically well defined mobility coefficient $\cal K$ whose value in this context has not been calculated before. Using methods from ordinary solid state and nuclear physics, a simple quantum mechanical treatment based on the independent particle approximation, is used here to formulate $\cal K$ as the phase space integral of the relevant group velocity over the neutron Fermi surface. The result can be described as an ``entrainment'' that changes the ordinary neutron mass m to a macroscopic effective mass per neutron that will be given -- subject to adoption of a convention specifying the precise number density n of the neutrons that are considered to be ``free'' -- by $m_\star=n/{\cal K}$. The numerical evaluation of the mobility coefficient is carried out for nuclear configurations of the ``lasagna'' and ``spaghetti'' type that may be relevant at the base of the crust. Extrapolation to the middle layers of the inner crust leads to the unexpected prediction that $m_\star$ will become very large compared with m.

The outer layers of a neutron star are supposed to be formed of a solid Coulomb lattice of neutron rich nuclei. At densities above neutron drip density (about one thousandth of nuclear saturation density), this lattice is immersed in a neutron fluid. Bragg scattering of those dripped neutrons by the nuclei which has been usually neglected is investigated, within a simple mean field model with Bloch type boundary conditions. The main purpose of this work is to provide some estimates for the entrainment coefficients, as required for hydrodynamical two fluid simulations of neutron star crust, which relate the momentum of one fluid to the particle currents of the other two fluids. The implications for the equilibrium neutron star crust structure are also briefly discussed.

As a follow up to articles dealing firstly with a convective variational formulation in a Milne-Cartan framework for non-dissipative multi fluid models, and secondly with various ensuing stress energy conservation laws and generalised virial theorems, this work continues a series showing how analytical procedures developed in the context of General Relativity can be usefully adapted for implementation in a purely Newtonian framework where they provide physical insights that are not so easy to obtain by the traditional approach based on a 3+1 space time decomposition. The present article describes the 4-dimensionally covariant treatment of various dissipative mechanisms, including viscosity in non-superfluid constituents, superfluid vortex drag, ordinary resistivity (mutual friction) between relatively moving non-superfluid constituents, and the transvective dissipation that occurs when matter is transformed from one constituent to another due to chemical disequilibrium such as may be produced by meridional circulation in neutron stars. The corresponding non dissipative limit cases of vortex pinning, convection, and chemical equilibrium are also considered.

This is the first of a series of articles showing how 4 dimensionally covariant analytical procedures developed in the context of General Relativity can be usefully adapted for application in a purely Newtonian framework where they provide physical insights (e.g. concerning helicity currents) that are not so easy to obtain by the traditional approach based on a 3+1 space time decomposition. After an introductory presentation of the relevant Milne spacetime structure and the associated Cartan connection, the essential principles are illustrated by application to the variational formulation of simple barotropic perfect fluid models. This variational treatment is then extended to conservative multiconstituent self gravitating fluid models of the more general kind that is needed for treating the effects of superfluidity in neutron stars.

The 4-dimensionally covariant approach to multiconstituent Newtonian fluid dynamics presented in the preceding article of this series is developed by construction of the relevant 4-dimensional stress energy tensor whose conservation in the non-dissipative variational case is shown to be interpretable as a Noether identity of the Milne spacetime structure. The formalism is illustrated by the application to homogeneously expanding cosmological models, for which appropriately generalised local Bernouilli constants are constructed. Another application is to the Iordanski type generalisation of the Joukowski formula for the Magnus force on a vortex. Finally, at a global level, a new (formally simpler but more generally applicable) version of the ``virial theorem'' is obtained for multiconsituent -- neutron or other -- fluid star models as a special case within an extensive category of formulae whereby the time evolution of variously weighted mass moment integrals is determined by corresponding space integrals of stress tensor components, with the implication that all such stress integrals must vanish for any stationary equilibrium configuration.