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.

The possible presence of amorphous and heterogeneous phases in the inner crust of a neutron star is expected to reduce the electrical conductivity of the crust, with potentially important consequences on the magneto-thermal evolution of the star. In cooling simulations, the disorder is quantified by an impurity parameter which is often taken as a free parameter. We aim to give a quantitative prediction of the impurity parameter as a function of the density in the crust,performing microscopic calculations including up-to-date microphysics of the crust. A multi-component approach is developed at finite temperature using a compressible liquid drop description of the ions with an improved energy functional based on recent microscopic nuclear models and optimized on extended Thomas-Fermi calculations. Thermodynamic consistency is ensured by adding a rearrangement term and deviations from the linear mixing rule are included in the liquid phase. The impurity parameter is consistently calculated at the crystallization temperature as determined in the one-component plasma approximation for the different functionals. Our calculations show that at the crystallization temperature the composition of the inner crust is dominated by nuclei with charge number around $Z \approx 40$, while the range of the $Z$ distribution varies from about 20 near the neutron drip to about 40 closer to the crust-core transition. This reflects on the behavior of the impurity parameter that monotonically increases with density up to around 40 in the deeper regions of the inner crust. Our study shows that the contribution of impurities is non-negligible, thus potentially having an impact on the transport properties in the neutron-star crust. The obtained values of the impurity parameter represent a lower limit; larger values are expected in the presence of non-spherical geometries and/or fast cooling dynamics.

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 impact of electron-capture (EC) cross sections on neutron-rich nuclei on the dynamics of core-collapse during infall and early post-bounce is studied performing spherically symmetric simulations in general relativity using a multigroup scheme for neutrino transport and full nuclear distributions in extended nuclear statistical equilibrium models. We thereby vary the prescription for EC rates on individual nuclei, the nuclear interaction for the EoS, the mass model for the nuclear statistical equilibrium distribution and the progenitor model. In agreement with previous works, we show that the individual EC rates are the most important source of uncertainty in the simulations, while the other inputs only marginally influence the results. A recently proposed analytic formula to extrapolate microscopic results on stable nuclei for EC rates to the neutron rich region, with a functional form motivated by nuclear-structure data and parameters fitted from large scale shell model calculations, is shown to lead to a sizable (16%) reduction of the electron fraction at bounce compared to more primitive prescriptions for the rates, leading to smaller inner core masses and slower shock propagation. We show that the EC process involves $\approx$ 130 different nuclear species around 86 Kr mainly in the N = 50 shell closure region, and establish a list of the most important nuclei to be studied in order to constrain the global rates.

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.

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.

The equation of state (EoS) of hot and dense matter is a fundamental input to describe static and dynamical properties of neutron stars, core-collapse supernovae and binary compact-star mergers. We review the current status of the EoS for compact objects, that have been studied with both ab-initio many-body approaches and phenomenological models. We limit ourselves to the description of EoSs with purely nucleonic degrees of freedom, disregarding the appearance of strange baryonic matter and/or quark matter. We compare the theoretical predictions with different data coming from both nuclear physics experiments and astrophysical observations. Combining the complementary information thus obtained greatly enriches our insights into the dense nuclear matter properties. Current challenges in the description of the EoS are also discussed, mainly focusing on the model dependence of the constraints extracted from either experimental or observational data (specifically, concerning the symmetry energy), the lack of a consistent and rigorous many-body treatment at zero and finite temperature of the matter encountered in compact stars (e.g. problem of cluster formation and extension of the EoS to very high temperatures), the role of nucleonic three-body forces, and the dependence of the direct URCA processes on the EoS.

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.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

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.

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.

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.

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.

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

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.

In stellar core-collapse events matter is heated and compressed to densities above nuclear matter saturation density. For progenitors stars with masses above about 25 solar masses, which eventually form a black hole, the temperatures and densities reached during the collapse are so high that a traditional description in terms of electrons, nuclei, and nucleons is no longer adequate. We present here an improved equation of state which contains in addition pions and hyperons. They become abundant in the high temperature and density regime. We study the different constraints on such an equation of state, coming from both hyperonic data and observations of neutron star properties. In order to test the zero-temperature versions, we perform numerical simulations of the collapse of a neutron star with such additional particles to a black hole. We discuss the influence of the additional particles on the thermodynamic properties within the hot versions of the equation of state and we show that in regimes relevant to core-collapse and black hole formation, the effects of pions and hyperons on pressure, internal energy and sound speed are not negligible.

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.

We investigate to what extent the temperature dependence of the nuclear symmetry energy can affect the neutronization of the stellar core prior to neutrino trapping during gravitational collapse. To this end, we implement a one-zone simulation to follow the collapse until beta equilibrium is reached and the lepton fraction remains constant. Since the strength of electron capture on the neutron-rich nuclei associated to the supernova scenario is still an open issue, we keep it as a free parameter. We find that the temperature dependence of the symmetry energy consistently yields a small reduction of deleptonization, which corresponds to a systematic effect on the shock wave energetics: the gain in dissociation energy of the shock has a small yet non-negligible value of about 0.4 foe (1 foe = 10^51 erg) and this result is almost independent from the strength of nuclear electron capture. The presence of such a systematic effect and its robustness under changes of the parameters of the one-zone model are significative enough to justify further investigations with detailed numerical simulations of supernova explosions.