Source author record

N. Chamel

N. Chamel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2026arXiv

Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh: V. The N2LO extension of the Skyrme EDF

We present BSkG5, the latest entry in the Brussels-Skyrme-on-a-Grid (BSkG) series and the first large-scale nuclear structure model based on next-to-next-to-leading order (N2LO) Skyrme energy density functional (EDF). By extending the traditional Skyrme EDF ansatz with central terms containing up to four gradients, we are able to combine an excellent global description of nuclear ground state properties with a stiff equation of state for pure neutron matter that is consistent with all astronomical observations of neutron stars. More precisely, the new model matches the accuracy of earlier BSkG models but with two parameters less: we achieve root-mean-square deviations of 0.649 MeV for 2457 atomic masses, 0.0267 fm for 810 charge radii, and 0.43 MeV for 45 primary fission barriers of actinide nuclei. We demonstrate that the complexities of N2LO EDFs are not insurmountable, even for demanding many-body calculations.

preprint2020arXiv

Superfluidity in Disordered Neutron Stars Crusts

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, $ρ_s/ρ$, at densities above that of neutron drip, $ρ_d \approx 4\times 10^{11}\,g/cm^3$. For an amorphous solid crust the suppression of $ρ_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.

preprint2020arXiv

Unified equations of state for cold non-accreting neutron stars with Brussels-Montreal functionals. II. Pasta phases in semi-classical approximation

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.

preprint2019arXiv

Crystallization of the inner crust of a neutron star and the influence of shell effects

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.

preprint2019arXiv

Unified equations of state for cold non-accreting neutron stars with Brussels-Montreal functionals. I. Role of symmetry energy

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.

preprint2016arXiv

Constraint on the internal structure of a neutron star from Vela pulsar glitches

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.

preprint2016arXiv

Giant Pulsar Glitches and the Inertia of Neutron-Star Crusts

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.

preprint2016arXiv

Landau quantization and neutron emissions by nuclei in the crust of a magnetar

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.

preprint2015arXiv

Analytical representations of unified equations of state for neutron-star matter

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 ρare given as functions of the baryon number density n_b. In the second version, P, ρ, 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.

preprint2015arXiv

Heat capacity of low density neutron matter: from quantum to classical regimes

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.

preprint2015arXiv

Neutron drip transition in accreting and nonaccreting neutron star crusts

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 $ρ_\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.

preprint2014arXiv

Linear response theory and neutrino mean free path using Brussels-Montreal Skyrme functionals

The Brussels-Montreal Skyrme functionals have been successful to describe properties of both finite nuclei and infinite homogeneous nuclear matter. In their latest version, these functionals have been equipped with two extra density-dependent terms in order to reproduce simultaneously ground state properties of nuclei and infinite nuclear matter properties while avoiding at the same time the arising of ferromagnetic instabilities. In the present article, we extend our previous results of the linear response theory to include such extra terms at both zero and finite temperature in pure neutron matter. The resulting formalism is then applied to derive the neutrino mean free path. The predictions from the Brussels-Montreal Skyrme functionals are compared with ab-initio methods.

preprint2014arXiv

On the Lie subalgebra of Killing-Milne and Killing-Cartan vector fields in Newtonian space-time

The Galilean (and more generally Milne) invariance of Newtonian theory allows for Killing vector fields of a general kind, whereby the Lie derivative of a field is not required to vanish but only to be cancellable by some infinitesimal Galilean (respectively Milne) gauge transformation. In this paper, it is shown that both the Killing-Milne vector fields, which preserve the background Newtonian space-time structure, and the Killing-Cartan vector fields, which in addition preserve the gravitational field, form a Lie subalgebra.

preprint2013arXiv

Collective excitations in neutron-star crusts

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.

preprint2013arXiv

Crustal Entrainment and Pulsar Glitches

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.

preprint2013arXiv

Does a hadron-quark phase transition in dense matter preclude the existence of massive 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.

preprint2013arXiv

How "free" are free neutrons in neutron-star crusts and what does it imply for pulsar glitches ?

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.

preprint2013arXiv

Low-energy collective excitations in the neutron star inner crust

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.

preprint2013arXiv

Phase transitions in dense matter and the maximum mass of neutron stars

The recent precise measurement of the mass of pulsar PSR J1614$-$2230, as well as observational indications of even more massive neutron stars, has revived the question of the composition of matter at the high densities prevailing inside neutron-star cores. We study the impact on the maximum possible neutron-star mass of an "exotic" core consisting of non-nucleonic matter. For this purpose, we study the occurrence of a first-order phase transition in nucleonic matter. Given the current lack of knowledge of non-nucleonic matter, we consider the stiffest possible equation of state subject only to the constraints of causality and thermodynamic stability. The case of a hadron-quark phase transition is discussed separately. The purely nucleonic matter is described using a set of unified equations of state that have been recently developed to permit a consistent treatment of both homogeneous and inhomogeneous phases. We then compute the mass-radius relation of cold nonaccreting neutron stars with and without exotic cores from the Tolman-Oppenheimer-Volkoff equations. We find that even if there is a significant softening of the equation of state associated with the actual transition to an exotic phase, there can still be a stiffening at higher densities closer to the center of the star that is sufficient to increase the maximum possible mass. However with quarks the maximum neutron-star mass is always reduced by assuming that the sound speed is limited by $c/\sqrt{3}$ as suggested by QCD calculations. In particular, by invoking such a phase transition, it becomes possible to support PSR J1614$-$2230 with a nucleonic equation of state that is soft enough to be compatible with the kaon and pion production in heavy-ion collisions.

preprint2013arXiv

Properties of the outer crust of strongly magnetized neutron stars from Hartree-Fock-Bogoliubov atomic mass models

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.

preprint2013arXiv

Symmetry energy: nuclear masses and neutron stars

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.

preprint2013arXiv

Unified description of dense matter in neutron stars and magnetars

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.

preprint2012arXiv

Equation of state of magnetar crusts from Hartree-Fock-Bogoliubov atomic mass models

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.

preprint2012arXiv

Inner crust of neutron stars with mass-fitted Skyrme functionals

The equation of state and composition of the inner crust of neutron stars at zero temperature are calculated, using the T = 0 version of the TETFSI (temperature-dependent extended Thomas-Fermi plus Strutinsky integral) method, for each of a family of three functionals based on Skyrme-type forces BSk19, BSk20 and BSk21, which are characterized by different degrees of symmetry-energy stiffness, and also for the SLy4 functional. We also solve the Tolman-Oppenheimer-Volkoff equations to calculate the distribution of mass within the inner crust. Qualitatively similar results are found for all four functionals, and in particular the number of protons per Wigner-Seitz cell is in all cases equal to 40 throughout the inner crust.

preprint2012arXiv

Structure of neutron stars with unified equations of state

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.

preprint2012arXiv

Superfluidity and entrainment in neutron-star crusts

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.

preprint2011arXiv

Masses of neutron stars and nuclei

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.

preprint2011arXiv

The decompression of the outer neutron star crust and r-process nucleosynthesis

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.

preprint2010arXiv

Breathing-mode measurements in Sn isotopes and isospin dependence of nuclear incompressibility

T. Li {\it et al.}[Phys. Rev. C {\bf 81}, 034309 (2010)] have analyzed their measured breathing-mode energies of some tin isotopes in terms of a first-order leptodermous expansion, and find for the symmetry-incompressibility coefficient $K_τ$ the value of -550 $\pm$ 100 MeV. Removing an approximation that they made, we find that the first-order estimate of $K_τ$ shifts to -661 $\pm$ 144 MeV. However, taking into account higher-order terms in the leptodermous expansion shows that the data are compatible with the significantly lower magnitudes indicated by both another experiment and some theoretical estimates.

preprint2010arXiv

Further explorations of Skyrme-Hartree-Fock-Bogoliubov mass formulas. XII: Stiffness and stability of neutron-star matter

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.

preprint2010arXiv

Self-interaction errors in nuclear energy density functionals

When applied to a single nucleon, nuclear energy density functionals may yield a non-vanishing internal energy thus implying that the nucleon is interacting with itself. It is shown how to avoid this unphysical feature for semi-local phenomenological functionals containing all possible bilinear combinations of local densities and currents up to second order in the derivatives. The method outlined in this Rapid Communication could be easily extended to functionals containing higher order terms, and could serve as a guide for constraining the time-odd part of the functional.

preprint2010arXiv

Spin and spin-isospin instabilities in asymmetric nuclear matter at zero and finite temperatures using Skyrme functionals

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.

preprint2010arXiv

Unified description of neutron superfluidity in the neutron-star crust with analogy to anisotropic multi-band BCS superconductors

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.

preprint2009arXiv

Recent breakthroughs in Skyrme-Hartree-Fock-Bogoliubov mass formulas

We review our recent achievements in the construction of microscopic mass tables based on the Hartree-Fock-Bogoliubov method with Skyrme effective interactions. In the latest of our series of HFB-mass models, we have obtained our best fit ever to essentially all the available mass data, by treating the pairing more realistically than in any of our earlier models. The rms deviation on the 2149 measured masses of nuclei with N and Z>8 has been reduced for the first time in a mean field approach to 0.581 MeV. With the additional constraint on the neutron-matter equation of state, this new force is thus very well-suited for the study of neutron-rich nuclei and for the description of astrophysical environments like supernova cores and neutron-star crusts.

preprint2009arXiv

Skyrme-Hartree-Fock-Bogoliubov nuclear mass formulas: Crossing the 0.6 MeV threshold with microscopically deduced pairing

We present a new Skyrme-Hartree-Fock-Bogoliubov nuclear-mass model in which the contact-pairing force is constructed from microscopic pairing gaps of symmetric nuclear matter and neutron matter calculated from realistic two- and three-body forces, with medium-polarization effects included. With the pairing being treated more realistically than in any of our earlier models, the rms deviation with respect to essentially all the available mass data falls to 0.581 MeV, the best value ever found within the mean-field framework. Since our Skyrme force is also constrained by the properties of pure neutron matter, this new model is particularly well suited for application to astrophysical problems involving a neutron-rich environment, such as the elucidation of the r process of nucleosynthesis, and the description of supernova cores and neutron-star crusts.