Search

We previously studied the inner crust and the pasta mantle of a neutron star within the 4th-order extended Thomas-Fermi (ETF) approach with consistent proton shell corrections added perturbatively via the Strutinsky integral (SI) theorem together with the contribution due to pairing. To speed up the computations and avoid numerical problems, we adopted parametrized nucleon density distributions. However, the errors incurred by the choice of the parametrization are expected to become more significant as the mean baryon number density is increased, especially in the pasta mantle where the differences in the energy per nucleon of the different phases are very small, typically a few keV. To improve the description of these exotic structures, we discuss the important features that a nuclear profile should fulfill and introduce two new parametrizations. Performing calculations using the BSk24 functional, we find that these parametrizations lead to lower ETF energy solutions for all pasta phases than the parametrization we adopted before and more accurately reproduce the exact equilibrium nucleon density distributions obtained from unconstrained variational calculations. Within the ETFSI method, all parametrizations predict the same composition in the region with quasi-spherical clusters. However, the two new parametrizations lead to a different mantle structure at mean baryon densities above about 0.07 fm^-3, at which point lasagna is energetically favored. Interestingly, spherical clusters reappear in the pasta region. The inverted pasta phases such as bucatini and Swiss cheese are still found in the densest region above the core in all cases., Comment: 12 pages, 9 figures, published in Phys. Rev. C

Our previous investigation of neutron-star crusts, based on the functional BSk24, led to a substantial reduction of the pasta mantle when Strutinsky integral and pairing corrections were added on top of the fourth-order extended Thomas-Fermi method (ETF). Here, our earlier calculations are widened to a larger set of functionals within the same family, and we find that the microscopic corrections weaken significantly the influence of the symmetry energy. In particular, the correlation observed at the pure ETF level between the density for the onset of pasta formation and the symmetry energy vanishes, not only for the $L$ coefficient but also for the symmetry-energy values at the relevant densities. Moreover, the inclusion of microscopic corrections results in a much lower abundance of pasta for all functionals., Comment: 11 pages, 4 figures

The accretion of matter onto the surface of a neutron star in a low-mass X-ray binary triggers X-ray bursts, whose ashes are buried and further processed thus altering the composition and the properties of the stellar crust. In this second paper of a series, the impact of accretion on the equation of state and on the global properties of neutron stars is studied in the framework of the nuclear energy-density functional theory. Considering ashes made of $^{56}$Fe, we calculated the equations of state using the same Brussels-Montreal nuclear energy-density functionals BSk19, BSk20, and BSk21, as those already employed for determining the crustal heating in our previous study for the same ashes. All regions of accreting neutron stars were treated in a unified and thermodynamically consistent way. With these equations of state, we determined the mass, radius, moment of inertia, and tidal deformability of accreted neutron stars and compared with catalyzed neutron stars for which unified equations of state based on the same functionals are available. The equation of state of accreted neutron stars is found to be significantly stiffer than that of catalyzed matter, with an adiabatic index $\Gamma \approx 4/3$ throughout the crust. For this reason, accreting neutron stars have larger radii. However, their crustal moment of inertia and their tidal deformability are hardly changed provided density discontinuities at the interface between adjacent crustal layers are properly taken into account. The enhancement of the stiffness of the equation of state of accreting neutron stars is mainly a consequence of nuclear shell effects, thus confirming the importance of a quantum treatment as stressed in our first study. With our previous calculations of crustal heating using the same functionals, we have thus obtained consistent microscopic inputs for simulations of accreting neutron stars., Comment: 10 pages, 7 figures. The tables of the equations of state are available at the CDS

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., Comment: 16 pages, 14 tables, 8 figures, accepted by Phys. Rev. C

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., Comment: 9 pages, 9 figures

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., Comment: 23 pages, 11 figures

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., Comment: 36 pages, 31 figures, 29 tables; corrected according to an erratum (2019, accepted)

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., Comment: 20 pages, 13 figures, submitted to A&A

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., Comment: 13 pages, 10 figures, 10 tables. In v.2, minor typos in the text and references are fixed. In v.3 and v.4, typos in Table 6 and Eq.(17), respectively, are fixed. In v.5, labels on Fig.7 are corrected

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., Comment: 10 pages, 9 figures, to appear in EPJA special volume on symmetry energy

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., Comment: Proceedings of IAUS 291 "Neutron Stars and Pulsars: Challenges and Opportunities after 80 years", J. van Leeuwen (ed.); 3 pages, 1 figure

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., Comment: Proceedings of IAUS 291 "Neutron Stars and Pulsars: Challenges and Opportunities after 80 years", J. van Leeuwen (ed.); 3 pages

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., Comment: 25 pages, 5 figures, revised version

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., Comment: 2 pages, to appear in the proceedings of the ERPM conference, Zielona Gora, Poland, April 2012

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., Comment: 4 pages, to appear in the proceedings of the ERPM conference, Zielona Gora, Poland, April 2012

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., Comment: 4 pages, to appear in the proceedings of the ERPM conference, Zielona Gora, Poland, April 2012

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., Comment: 35 pages, 13 figures, accepted for publication in Physical Review C

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., Comment: 9 pages, 8 figures

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., Comment: 16 pages, 3 figures. Contribution for the book "50 years of nuclear BCS", edited by R.A. Broglia and V. Zelevinsky

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., Comment: 11 pages, 4 figures; Physical Review C in press

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., Comment: 45 pages, 16 figures, accepted for publication in Physical Review C

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_{\tau}$ the value of -550 $\pm$ 100 MeV. Removing an approximation that they made, we find that the first-order estimate of $K_{\tau}$ 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., Comment: 6 pages, 1 figure

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., Comment: 15 pages, 31 figures, accepted for publication in Physical Review C

After a brief review of the Hartree-Fock-Bogoliubov method with Skyrme effective interactions, we show how it can be applied to the description of various nuclear systems, from finite nuclei to neutron-star crusts., Comment: 7 pages, 1 figure, talk presented at the 28th International Workshop on Nuclear Theory, Rila Mountains, Bulgaria, June 22-27, 2009

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., Comment: accepted for publication in Physical Review C

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., Comment: Proceedings of the Fifth International Conference on Exotic Nuclei and Atomic Masses, September 7-13 2008, Ryn (Poland). To appear in the European Physical Journal A

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.

In this latest of our series of Skyrme-HFB mass models, HFB-16, we introduce the new feature of requiring that the contact pairing force reproduce at each density the $^1S_0$ pairing gap of neutron matter as determined in microscopic calculations with realistic nucleon-nucleon forces. We retain the earlier constraints on the Skyrme force of reproducing the energy-density curve of neutron matter, and of having an isoscalar effective mass of $0.8M$ in symmetric infinite nuclear matter at the saturation density; we also keep the recently adopted device of dropping Coulomb exchange. Furthermore, the correction term for the spurious energy of collective motion has a form that is known to favour fission barriers that are in good agreement with experiment. Despite the extra constraints on the effective force, we have achieved a better fit to the mass data than any other mean field model, the rms error on the 2149 measured masses of nuclei with $N$ and $Z \ge$ 8 having been reduced to 0.632 MeV; the improvement is particularly striking for the most neutron-rich nuclei. Moreover, it turns out that even with no flexibility at all remaining for the pairing force, the spectral pairing gaps that we find suggest that level densities in good agreement with experiment should be obtained. This new force is thus particularly well-suited for astrophysical applications, such as stellar nucleosynthesis and neutron-star crusts., Comment: 38 pages, 9 figures accepted for publication in Nuclear Physics A

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., Comment: 41 pages, 6 figures

Starting from HFB-6, we have constructed a new mass table, referred to as HFB-8, including all the 9200 nuclei lying between the two drip lines over the range of Z and N > 6 and Z < 122. It differs from HFB-6 in that the wave function is projected on the exact particle number. Like HFB-6, the isoscalar effective mass is constrained to the value 0.80 M and the pairing is density independent. The rms errors of the mass-data fit is 0.635 MeV, i.e. better than almost all our previous HFB mass formulas. The extrapolations of this new mass formula out to the drip lines do not differ significantly from the previous HFB-6 mass formula., Comment: 9 pages, 7 figures, accepted for publication in Phys. Rev. C

We have constructed four new complete mass tables, referred to as HFB-4 to HFB-7, each one including all the 9200 nuclei lying between the two drip lines over the range of Z and N>8 and Z<120. HFB-4 and HFB-5 have the isoscalar effective mass M*_s$ constrained to the value 0.92 M, with the former having a density-independent pairing, and the latter a density-dependent pairing. HFB-6 and HFB-7 are similar, except that M*_s is constrained to 0.8 M. The rms errors of the mass-data fits are 0.680, 0.675, 0.686, and 0.676 MeV, respectively, almost as good as for the HFB-2 mass formula, for which M*_s was unconstrained. However, as usual, the single-particle spectra depend significantly on M*_s. This decoupling of the mass fits from the fits to the single-particle spectra has been achieved only by making the cutoff parameter of the delta-function pairing force a free parameter. An improved treatment of the center-of-mass correction was adopted, but although this makes a difference to individual nuclei it does not reduce the overall rms error of the fit. The extrapolations of all four new mass formulas out to the drip lines are essentially the same as for the original HFB-2 mass formula., Comment: 12 pages revtex, 9 eps figures, accepted for publication in Phys. Rev. C

Using the ETFSI (extended Thomas-Fermi plus Strutinsky integral) method, we have calculated the fission barriers of nearly 2000 exotic nuclei, including all the neutron-rich nuclei up to A=318 that are expected to be relevant to the r-process, and all the superheavy nuclei in the vicinity of N=184, with Z<=120. Our calculations were performed with the Skyrme force SkSC4, which was determined in the ETFSI-1 mass fit. For proton-deficient nuclei in the region of N=184 we find the barriers to be much higher than previously believed, which suggests that the r-process path might continue to mass numbers well beyond 300. For the superheavy nuclei we typically find barrier heights of 6-7 MeV., Comment: 33 pages, 6 figures, to appear in Nucl. Phys. A