Thermal radiation of neutron stars in soft X-ray transients (SXTs) in a quiescent state is believed to be powered by the heat deposited in the stellar crust due to nuclear reactions during accretion. Confronting observations of this radiation with simulations helps to verify theoretical models of the dense matter in neutron stars. We simulate the thermal evolution of the SXTs with theoretical models of the equation of state and composition of the accreted crust. The new family of such models were recently developed within a thermodynamically consistent approach by modeling the nuclear evolution of an accreted matter as it sinks toward the stellar center, starting from representative thermonuclear ash compositions. The crust cooling curves computed with the traditional and modern theory are compared with observations of SXTs MXB 1659-29 and IGR J17480-2446. We show that the new and traditional models of the accreted neutron star crusts are similar in their capability to explain the thermal evolution of neutron stars in SXTs. Both kinds of models require inclusion of additional ingredients not supplied by the current theory, such as the shallow heating and variation of thermal conductivity, to fit observations.

Historically, a sequence of nuclear pasta shapes was predicted to appear in the deepest region of the inner crust of a neutron star within the compressible liquid-drop picture, when the filling fraction u exceeds some threshold values. However, later calculations showed that these values depend on the details of the liquid-drop model. Here we investigate the existence of pasta in neutron stars within the semiclassical extended Thomas-Fermi approach using various generalized Skyrme functionals. The filling fractions for the different transitions are found to be quasi-universal, unlike the pasta density ranges governed by the symmetry energy at relevant densities. In particular, pasta emerge at u_sp≈0.13-0.15. By applying a simplified stability criterion within the liquid-drop framework, we show that these values of u_sp can be explained by the nuclear curvature correction. In this way, the abundance of pasta can be easily estimated. This criterion can also be used to optimize the search of pasta within the more realistic extended Thomas-Fermi approach.

Self-consistent mean-field methods with Skyrme-type effective interactions and semiclassical approximations, such as the Thomas-Fermi approach and its extensions are particularly well-suited for describing in a thermodynamically consistent way the various phases of the dense matter present in the interior of neutron stars. These methods have been applied to predict the composition of the different regions, including the inner crust constituted by nuclear clusters coexisting with free neutrons and electrons. Because of the computational cost, the energy is typically calculated for a few selected average baryon number densities, and the results are interpolated to obtain the pressure numerically. However, this may introduce systematic errors in the calculations of the global structure of a neutron star and its dynamical evolution. In this paper, we show how the full equation of state can be consistently calculated within the same framework by deriving exact formulas for the chemical potentials and for the pressure that can be easily implemented in existing computer codes. These formulas are applicable to both catalyzed and accreted crusts. We discuss in each case the suitable conditions to impose to determine the composition. Numerical examples are also presented and discussed. Results from refined calculations of the BSk24 equation of state for the inner crust of nonaccreted neutron stars and the corresponding adiabatic index are provided.

The elastic properties of the neutron star crust are thought to play a crucial role in various phenomena of neutron stars (glitches, oscillations, gravitational wave emission) and should be described quantitatively to model these phenomena. The fundamental problem of this description is associated with the polycrystalline nature of the crust: similar to terrestrial materials, the elastic moduli, strictly speaking, depend on the shape and orientation of crystallites, but for the crust, they are unknown. As a result, some assumptions are generally required to predict the elastic properties or constrain their possible range. In this paper, we follow the commonly believed assumption that the crust is (locally) isotropic, which allows us to describe elastic properties by two (effective) parameters: bulk and shear moduli. The bulk modulus is well determined by the Voigt-Reuss bounds, and we constrain the shear modulus by applying, for the first time in astrophysics of compact stars, the variational Hashin-Shtrikman approach, based on the additional assumption that there are no correlations in the orientation of crystallites. We analyse the Hashin-Shtrikman bounds for the one-component crust taking into account the electron screening and the motion of the nuclei, and for two-component static crystals. In particular, we demonstrate that within applied assumptions the effective shear modulus should be lower than the Voigt estimate, typically applied in the astrophysical literature.

Observations of accreting neutron stars are widely used to constrain the microphysical properties of superdense matter. A key ingredient in this analysis is the heating associated with nuclear reactions in the outer layers of the neutron star (crust), as well as the equation of state and composition of these layers. As recently shown, the neutron hydrostatic/diffusion (nHD) condition is valid in the inner part of the crust, where some of the neutrons are not bound to the nuclei, and this condition should be properly incorporated into crustal models. Here we construct models of the accreted crust of a neutron star, taking into account the nHD condition and proton shell effects in nuclei. For numerical illustration, we employ the recently proposed compressible liquid drop model, which incorporates shell effects. However, our approach is general and can also be used in future studies relying on more sophisticated nuclear physics models.

Parametrized nucleon density distributions are widely employed for the calculation of the properties of atomic nuclei and dense inhomogeneous matter in compact stars within the Thomas-Fermi method and its extensions. We show that the use of insufficiently smooth parametrizations may deteriorate the accuracy of this method. We discuss and clarify the smoothness condition using the example of the so-called "nuclear pasta" in the neutron star mantle.

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.

We model the nuclear evolution of an accreted matter as it sinks toward the stellar centre, in order to find its composition and equation of state. To this aim, we developed a simplified reaction network that allows for redistribution of free neutrons in the inner crust to satisfy the recently suggested neutron hydrostatic and diffusion equilibrium condition. We analyse the main reaction pathways for the three representative thermonuclear ash compositions: Superburst, Kepler, and Extreme rp. In contrast to the previous results, which neglect redistribution of free (unbound) neutrons in the inner crust, the most significant reactions in our calculations are neutron captures and electron emissions. The pycnonuclear fusion plays some role only for Kepler ashes. For the direct application of our results in astrophysical codes we present profiles of the average charge, & LANGBRAC;Z & rang;, impurity parameter, Q(imp), and equation of state for a set of models, parametrized by the pressure at the outer-inner crust interface. Typically, for Superburst ashes Q(imp) & AP; 1 - 4, while for Kepler ashes Q(imp) decreases from & AP;23 at the outer-inner crust interface to & AP;5 at the end of our simulation (the corresponding density equals & rho;(dc) & AP; 2 x 10(12) g cm(-3)). At the same time, for Extreme rp ashes Q(imp) remains large & AP;30 - 35 in the considered inner crust region. Our results are important for modelling the thermal relaxation of transiently accreting neutron stars after the end of the outburst.

ABSTRACT The elasticity of neutron star crust is important for adequate interpretation of observations. To describe elastic properties one should rely on theoretical models. The most widely used is Coulomb crystal model (system of point-like charges on neutralizing uniform background), in some works it is corrected for electron screening. These models neglect finite size of nuclei. This approximation is well justified except for the innermost crustal layers, where nuclei size becomes comparable with the inter-nuclear spacing. Still, even in those dense layers it seems reasonable to apply the Coulomb crystal result, if one assumes that nuclei are spherically symmetric: Coulomb interaction between them should be the same as interaction between point-like charges. This argument is indeed correct; however, as we point here, shear of crustal lattice generates (microscopic) quadrupole electrostatic potential in a vicinity of lattice cites, which induces deformation on the nuclei. We analyse this problem analytically within compressible liquid drop model. In particular, for ground state crust composition the effective shear modulus is reduced for a factor of $1-u^{5/3}/(2+3\, u-4\, u^{1/3})$, where u is the ratio of the nuclei volume to the volume of the cell. This result is universal, i.e. it does not depend on the applied nucleon interaction model within applied approach. For the innermost layers of inner crust u ∼ 0.2 leading to reduction of the shear modulus by $\sim 25{{\ \rm per\ cent}}$, which can be important for correct interpretation of quasi-periodic oscillations in the tails of magnetar flares.

Neutron stars are the densest objects in the Universe. They have a microscopically homogeneous core and heterogeneous crust. In particular, there may be a specific layer inside neutron stars, the mantle, which consists of substantially non-spherical nuclei immersed in a background of relativistic degenerate electrons and quasi-free neutrons. In this paper, we reconsider the transverse shear modulus for cylindrical phases of the mantle within the framework of the compressible liquid drop model. We demonstrate that transverse shearing affects the shape of nuclear clusters: their cross-section becomes elliptical. This effect reduces the respective elastic constant. Using a simple model, we perform all derivations analytically and obtain the expression for the transverse shear modulus, which can be useful for astrophysical applications.

Thermal emission of neutron stars in soft X-ray transients (SXTs) in a quiescent state is believed to be powered by the heat deposited in the stellar crust due to nuclear reactions during accretion (deep crustal heating paradigm). Confronting observations of SXTs with simulations helps to verify theoretical models of the dense matter in the neutron stars. Usually, such simulations were carried out assuming that the free neutrons and nuclei in the inner crust move together. A recently proposed thermodynamically consistent approach allows for independent motion of the free neutrons. We simulate the thermal evolution of the SXTs within the thermodynamically consistent approach and compare the results with the traditional approach and with observations. For the latter, we consider a collection of quasi-equilibrium thermal luminosities of the SXTs in quiescence and the observed neutron star crust cooling in SXT MXB 1659-29. We test different models of the equation of state and baryon superfluidity and take into account additional heat sources in the shallow layers of neutron-star crust (the shallow heating). We find that the observed quasi-stationary thermal luminosities of the SXTs can be equally well fitted using the traditional and thermodynamically consistent models, provided that the shallow heat diffusion into the core is taken into account. The observed crust cooling in MXB 1659-29 can also be fitted in the frames of both models, but the choice of the model affects the derived parameters responsible for the thermal conductivity in the crust and for the shallow heating.

ABSTRACT Neutron stars in low-mass X-ray binaries are thought to be heated up by accretion-induced exothermic nuclear reactions in the crust. The energy release and the location of the heating sources are important ingredients of the thermal evolution models. Here, we present thermodynamically consistent calculations of the energy release in three zones of the stellar crust: at the outer–inner crust interface, in the upper layers of the inner crust (up to the density ρ ≤ 2 × 1012 g cm−3), and in the underlying crustal layers. We consider three representative models of thermonuclear ashes (superburst, extreme rp, and Kepler ashes). The energy release in each zone is parametrized by the pressure at the outer–inner crust interface, Poi, which encodes all uncertainties related to the physics of the deepest inner-crust layers. Our calculations allow us to set new theoretical lower limits on the net energy release (per accreted baryon): Q ≳ 0.28 MeV for extreme rp ashes and Q ≳ 0.43–0.51 MeV for superburst and Kepler ashes. Our results can be directly incorporated into numerical codes and provide an opportunity to constrain Poi by comparing thermal evolution models of accreting neutron stars with observations.

Neutron stars are the densest objects in the Universe. In this paper, we consider the so-called inner crust—the layer where neutron-excess nuclei are immersed in the degenerate gas of electrons and a sea of quasi-free neutrons. It was generally believed that spherical nuclei become unstable with respect to quadrupole deformations at high densities, and here, we consider this instability. Within the perturbative approach, we show that spherical nuclei with equilibrium number density are, in fact, stable with respect to infinitesimal quadrupole deformation. This is due to the background of degenerate electrons and associated electrostatic potential, which maintain stability of spherical nuclei. However, if the number of atomic nuclei per unit volume is much less than the equilibrium value, instability can arise. To avoid confusion, we stress that our results are limited to infinitesimal deformations and do not guarantee strict thermodynamic stability of spherical nuclei. In particular, they do not exclude that substantially non-spherical nuclei (so-called pasta phase) represent a thermodynamic equilibrium state of the densest layers of the neutron star crust. Rather, our results point out that spherical nuclei can be metastable even if they are not energetically favourable, and the timescale of transformation of spherical nuclei to the pasta phases should be estimated subsequently.

Nuclear pasta phases in the neutron stars mantle can affect the mechanical and transport properties of superdense matter, thus playing an important role in the dynamics and evolution of neutron stars. In this paper, we compare results obtained by the Extended Thomas–Fermi (ETF) method with the compressible liquid drop model (CLDM), based on the thermodynamically consistent description of the surface properties calculated for the two-phase plane interface and the same energy-density functional (for numerical illustration, we applied the Skyrme-type functional SLy4). Our ETF calculations found that pasta phases in cylindrical form cover a significant crustal region (both normal and inverse phases, aka spaghetti and bucatini are presented). Meanwhile, within the applied CLDM framework, which includes the thermodynamically required effect of neutron adsorption on the cluster’s surface but neglects curvature corrections, only the spaghetti phase was found to be energetically favorable in the small density range prior to crust–core transition. On the other hand, the recent CLDM of Dinh Thi et al., 2021, which, on the contrary, accounts for curvature term but neglects neutron adsorption, predicts pasta phase onset in better agreement with the ETF. This fact highlights the importance of the curvature effects and allows counting on the potential validity of the CLDMs as a convenient, transparent and accurate tool for investigation of the pasta-phase properties.

ABSTRACT The main contribution to the effective shear modulus of neutron star crust can be calculated within Coulomb solid model and can be approximated by simple analytical expression for arbitrary (even multicomponent) composition. Here I consider correction associated with electron screening within Thomas–Fermi approximation. In particular, I demonstrate that for relativistic electrons (density ρ > 106 g cm−3) this correction can be estimated as $\delta \mu _\mathrm{eff}^\mathrm{V}= -9.4\times 10^{-4}\sum _Z n_Z Z^{7/3} e^2/a_\mathrm{e},$ where summation is taken over ion species, nZ is number density of ions with charge Ze, kTF is Thomas–Fermi screening wavenumber. Finally, ae = (4πne/3)−1/3 is electron sphere radius. Quasi-neutrality condition ne = ∑ZZnZ is assumed. This result holds true for arbitrary (even multicomponent and amorphous) matter and can be applied for neutron star crust and (dense) cores of white dwarfs. For example, the screening correction reduces shear modulus by ∼9 per cent for Z ∼ 40, which is typical for inner layers of neutron star crust.

Using path-integral Monte Carlo (PIMC) simulations, we have calculated energy of a crystal composed of atomic nuclei and uniform incompressible electron background in the temperature and density range, covering fully ionized layers of compact stellar objects, white dwarfs, and neutron stars, including the high-density regime, where ion quantization is important. We have approximated the results by convenient analytic formulae, which allowed us to integrate and differentiate the energy with respect to temperature and density to obtain various thermodynamic functions such as Helmholtz free energy, specific heat, pressure, entropy etc. In particular, we have demonstrated, that the total crystal specific heat can exceed the well-known harmonic lattice contribution by a factor of 1.5 due to anharmonic effects. By combining our results with the PIMC thermodynamics of a quantum Coulomb liquid, updated in the present work, we were able to determine density dependences of such melting parameters as the Coulomb coupling strength at melting, latent heat, and a specific heat jump. Our results are necessary for realistic modelling of thermal evolution of compact degenerate stars.

Neutron stars can appear as sources of different nature. In this paper we address the observability of a hypothetical class of neutron stars—HOt and Fast Non-Accreting Rotators, HOFNARs. These objects are heated due to the r-mode instability. With surface temperatures ∼106 K they are expected to be thermal soft X-ray emitters. We perform a population synthesis modeling of HOFNARs to predict the number of potentially detectable sources in the eROSITA all-sky survey. For surface temperatures ∼106 K we obtain ∼500 sources above the detection limit 0.01 cts s−1 and ∼100 easier identifiable sources with >0.1 cts s−1. Temperatures ≳1.2 × 106 K start to be in contradiction with non-detection of HOFNARs by ROSAT. Only for T ≲ 5 × 105 K numbers predicted for eROSITA turn out to be so low that identification does not look possible. We conclude that eROSITA has good chances to discover HOFNARs, if they exist. Non-detection will put very stringent limits on the properties of this type of neutron stars.

ABSTRACT The deep crustal heating, associated with exothermal nuclear reactions, is believed to be a key parameter for describing the thermal evolution of accreting neutron stars. In this paper, we present the first thermodynamically consistent calculations of the crustal heating for realistic compositions of thermonuclear ashes. In contrast to previous studies based on the traditional approach, we account for neutron hydrostatic/diffusion (nHD) equilibrium condition imposed by superfluidity of neutrons in a major part of the inner crust and rapid diffusion in the remaining part of the inner crust. We apply a simplified reaction network to model nuclear evolution of various multi-component thermonuclear burning ashes (superburst, KEPLER, and extreme rp-process ashes) in the outer crust and calculate the deep crustal heating energy release Q, parametrized by the pressure at the outer–inner crust interface, Poi. Using the general thermodynamic arguments, we set a lower limit on Q, Q ≳ 0.13−0.2 MeV per baryon (an actual value depends on the ash composition and the employed mass model).

Neutron stars are superdense compact astrophysical objects. The central region of the neuron star (the core) consists of locally homogeneous nuclear matter, while in the outer region (the crust) nucleons are clustered. In the outer crust these nuclear clusters represent neutron-rich atomic nuclei and all nucleons are bound within them. Whereas in the inner crust some neutrons are unbound, but nuclear clusters still keeps generally spherical shape. Here we consider the region between the crust and the core of the star, so-called mantle, where non-spherical nuclear clusters may exist. We apply compressible liquid drop model to calculate the energy density for several shape types of nuclear clusters. It allows us to identify the most energetically favorable configuration as function of baryon number density. Employing four Skyrme-type forces (SLy4 and BSk24, BSk25, BSk26), which are widely used in the neutron star physics, we faced with strong model dependence of the ground state composition. In particular, in agreement with previous works within liquid drop model, mantle is absent for SLy4 (nuclear spheres directly transit into homogeneous nuclear matter; exotic nuclear shapes do not appear).

We show that, in order to determine the equation of state of the inner crust of an accreting neutron star, one should minimize not the Gibbs free energy, as it is generally assumed in the literature, but a different thermodynamic potential Ψ, which tends to the minimum at fixed pressure and neutron chemical potential. Once this potential is specified, one can calculate the heat-release distribution in the stellar crust due to nonequilibrium nuclear reactions induced by accretion of matter onto the neutron-star surface. The results are important for adequate modeling of the accreted crust and interpretation of the observations of accreting neutron stars in low-mass X-ray binaries.