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

Sep 30, 2013
13 pages
Published in:
  • Astron.Astrophys. 560 (2013) A48
e-Print:

Citations per year

201320162019202220250510152025
Abstract: (arXiv)
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.
Note:
  • 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
  • [1998]
    (c) Li & Schulze
  • V18

  • These NeuM constraints make the EDFs BSk19-21 suitable for application to the neutron-rich environments encountered in many different astrophysical situations, and our making three different such EDFs available reflects the current lack of knowledge of the high-density behavior of dense matter. In addition, these three EDFs were also constrained to reproduce several other properties of homogeneous nuclear matter as obtained from many-body calculations using realistic two- and threenucleon interactions
    • among those, the ratio of the isoscalar effective mass m∗ s to bare nucleon mass m in symmetric nuclear matter at saturation was set to the realistic value of 0.8, and all three EDFs predict a neutron effective mass that is higher than the proton effective mass in neutron-rich matter, as found both experimentally and from microscopic calculations. Various properties of the BSk19, BSk20, and BSk21 EDFs are summarized in Table 1 namely: the rms deviations to the 2149 measured atomic masses σ(matom) and to the 782 measured charge radii σ(Rch), the energy per nucleon of symmetric nuclear matter at saturation density av, the baryon number density at saturation nb,0, the incompressibility of symmetric nuclear matter at saturation Kv (which was required to fall in the experimental range 240 ± 10 MeV, according to Colò et al
      • Goriely
    • [2004]
      the symmetry energy coefficient J and its slope L, the isospin compressibility Kτ, the isoscalar and isovector effective masses m∗ s and m∗ v relative to the bare nucleon mass m, and the limiting baryonic density ncaus after which the EoSs of neutronstar matter violate causality. The last line indicates the NeuM
      • EoS to which each EDF was fitted
        • The Brussels-Montreal EDFs BSk19, BSk20, and BSk21 were used to compute the EoS of all regions of a neutron star
          • Following the BPS model the outer crust was assumed to consist of fully ionized atoms arranged in a bodycentered cubic lattice at zero temperature. The EoSs of the outer crust were calculated using either experimental atomic masses when available or theoretical masses obtained from the HFB mass models constructed with the BSk19, BSk20, and BSk21
            • Baym
          • EDFs, as appropriate (see for details). For the inner crust, where neutron-proton clusters coexist with free neutrons, the kinetic-energy part of the appropriate EDF was calculated using the semi-classical extended Thomas-Fermi method with proton quantum shell corrections added via the Strutinsky integral theorem
            • Pearson
          • neutron shell effects, which are known to be much smaller, were neglected (see for details). This method is a computationally very fast approximation to the full Hartree-Fock equations. In order to further reduce the computations, nuclear clusters were assumed to be spherical, and parametrized nucleon distributions were introduced. Finally, the electrostatic energy was calculated using the Wigner-Seitz approximation, and pairing effects were neglected. The overall resulting errors in the EoS of the inner crust were found to be about
            • Pearson
          • 5% at the neutron-drip point. The EoSs of the core, assumed to consist of homogeneous beta-equilibrated matter made of nucleons and leptons (electrons and muons), were calculated essentially analytically from the adopted EDF
            • Goriely
          • In a recent paper analyzed 240 Skyrme parametrizations, including BSk EDFs, by comparing their predicted properties of symmetric nuclear matter and pure NeuM to some empirical constraints. On this basis, these authors rejected
            • Dutra
          • SM3

            • Dutra
            • [2012]
              a constraint on the EoS of symmetric nuclear matter obtained from the analysis of particle-flow measurements in Au-Au collisions that can be represented by a band in the plot of pressure vs. density. Now the pressures calculated with the EDFs BSk19-21 fall within this band over 80% of the density range and never deviate by more than about 20% from those inferred by Danielewicz et al (see, e.g., Fig
              • Danielewicz
            • [3]
              of 2011). Nevertheless, Dutra et al (2012) reject these functionals on the grounds that the calculated EoSs do not fall within the band of Danielewicz et al over 95% of the density range. This criterion is quite arbitrary, and without any sound statistical foundation. In addition, it has to be noted that the interpretation by Danielewicz et al of the raw experimental data is subject to uncertainties of two different kinds of model dependence: i) the transport model that determines the particle flow in a given collision experiment
              • Chamel
            • [3]
              ii) extrapolation from the charge-asymmetricAu + Au system ((N-Z)/A ≈ 20%) to symmetric nuclear matter. As demonstrated by
              • Gale
              • [1990]
                the flow is generated during the early nonequilibrium stages of the collision, whereas the interpretation in terms of the
                • EoS by is an equilibrium version of a simplified momentum-dependent interaction. Viewed in this light, the deviations of the pressures obtained with BSk19-21 from the values of Danielewicz et al are altogether insignificant
                  • Danielewicz
                • PNM2

                  • Dutra
                • SM3

                  • MIX5

                    • Dutra
                  • However, the situation regarding the symmetry energy still remains a matter of debate, different experiments and/or theoretical calculations leading to different and sometimes contradicting predictions, as shown e.g. in Fig. 12 of Lattimer (see also
                    • Tsang
                  • Sect. IIIB). Different
                    • Goriely
                  • Skyrme EDFs can thus be selected depending on the experiments or many-body calculations. For example, the values of L obtained with BSk19-21 EDFs are all compatible with the range of values L = 55 ± 25 MeV deduced by and Warda et al. from measurements of the neutron-skin thickness in nuclei. The L coefficients obtained with BSk20 and
                    • Centelles
                  • BSk21 are also in agreement with the values 36 MeV< L <