Next Article in Journal
Uncovering the First AGN Jets with AXIS
Next Article in Special Issue
Strange Quark Stars: The Role of Excluded Volume Effects
Previous Article in Journal
Detecting Wandering Intermediate-Mass Black Holes with AXIS in the Milky Way and Local Massive Galaxies
Previous Article in Special Issue
Measuring the Lense–Thirring Orbital Precession and the Neutron Star Moment of Inertia with Pulsars
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nuclear Matter Equation of State in the Brueckner–Hartree–Fock Approach and Standard Skyrme Energy Density Functionals

by
Isaac Vidaña
1,†,
Jérôme Margueron
2,3,*,† and
Hans-Josef Schulze
1,†
1
Sezione di Catania, Dipartimento di Fisica e Astronomia “Ettore Majorana“, Istituto Nazionale di Fisica Nucleare, Universitá di Catania, Via Santa Sofa 64, I-95123 Catania, Italy
2
Institut de Physique Nucléaire, Université Lyon 1 and CNRS/IN2P3, 4 rue Enrico Fermi, 69622 Villeurbanne, France
3
International Research Laboratory on Nuclear Physics and Astrophysics, Michigan State University and CNRS, East Lansing, MI 48824, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Universe 2024, 10(5), 226; https://doi.org/10.3390/universe10050226
Submission received: 29 March 2024 / Revised: 1 May 2024 / Accepted: 5 May 2024 / Published: 17 May 2024
(This article belongs to the Special Issue Studies in Neutron Stars)

Abstract

:
The equation of state of asymmetric nuclear matter as well as the neutron and proton effective masses and their partial-wave and spin–isospin decomposition are analyzed within the Brueckner–Hartree–Fock approach. Theoretical uncertainties for all these quantities are estimated by using several phase-shift-equivalent nucleon–nucleon forces together with two types of three-nucleon forces, phenomenological and microscopic. It is shown that the choice of the three-nucleon force plays an important role above saturation density, leading to different density dependencies of the energy per particle. These results are compared to the standard form of the Skyrme energy density functional, and we find that it is not possible to reproduce the BHF predictions in the ( S , T ) channels in symmetric and neutron matter above saturation density, already at the level of the two-body interaction, and even more including the three-body interaction.

1. Introduction

A new era of multi-messenger astronomy began with the first detection of gravitational wave signals, also known as GW170817 or AT2017gfo for the electromagnetic emission originating from the kilonova formed by the binary neutron star merger [1,2]. This new astronomy emphasizes the crucial role of neutron stars (NSs), which are prototypes for the most extreme phases of matter where the strong, the weak, the gravitational and the electromagnetic interactions can be studied in regimes that cannot be explored on Earth. In this paper, we are analyzing the properties of the dense matter equation of state (EOS) from microscopic calculations with a focus on the strong three-body interaction, and we discuss the link with empirical energy density functionals (EDFs).
In recent years, the determination of nuclear EDFs has incorporated more and more information deduced from microscopic calculations of nuclear matter. Usually, one specific or a very limited number of microscopic calculations are considered in the fitting protocols. However, the impact of the choices of the bare nuclear interaction, the regularization scheme, the experimental error in the phase shift data, and of the many-body framework is important, and has been cautiously addressed in the chiral effective field theory ( χ EFT) [3,4,5,6,7] in order to theoretically predict the nuclear EOS, including error estimates. There is, however, an upper limit on density (presently estimated to be (1–2) n sat with n sat 0.155 ± 0.005 fm 3 [8]), where the χ EFT breaks down and also the error estimates become meaningless [9,10,11].
At increasing density—and already around n sat —the role of nuclear three-body forces (3BF) in nuclear matter cannot be ignored. For traditional meson-exchange potential models, the compatible 3BF correction at saturation density is of limited magnitude, a few MeV, which represents (15–20)% of the total energy, in order to reproduce the saturation of nuclear matter. With the advent of regularized soft-core nuclear interactions, such as V low - k [12,13] as well as chiral effective field theory ( χ EFT) [3,4,5,6,7], a new paradigm has emerged. In short, any modification of the two-body interaction, by applying a unitary transformation, generates many-body interactions. The 3BF cannot be viewed as independent of the two-body force (2BF), but instead, it contains a part of the repulsive hard-core interaction. Regularized soft-core interactions need stronger many-body interactions to saturate, compared to the hard-core ones.
Building bridges between EDFs for the in-medium nuclear interaction and microscopic approaches of the nuclear many-body problem is, therefore, fundamental for a consistent understanding of nuclei, neutron stars, core-collapse supernovae, and kilonovae. Microscopic approaches, on one hand, have the advantage of being based on realistic 2BFs that reproduce with high precision the scattering phase shifts and the deuteron properties, and include the isospin asymmetry dependence in a natural way. However, the direct implementation of microscopic approaches in finite nuclei is not an easy task. A realistic 2BF cannot be used directly in nuclei. Its hard-core repulsion requires a resummation, as uniform matter and 3BFs are required to describe both the energy and charge radii of finite nuclei. EDF models, on the other hand, are usually based on effective density-dependent interactions with parameters frequently fitted to reproduce the global properties of nuclei and properties of symmetric and neutron matter [14,15,16,17,18]. These effective interactions capture the essence of the in-medium nuclear interaction as well as the effect of the complex correlations. Their predictions at high densities and isospin asymmetries, however, should be taken with care and discussed in view of the considered constraints. Combining phenomenological models with microscopic approaches would, therefore, help in setting up a nuclear model based on realistic 2BFs and 3BFs that can simultaneously describe the properties of nuclei and those of infinite nuclear matter in a large range of densities and isospin asymmetries.
Two major sources of uncertainty arise, however, when trying to build links between EDFs and microscopic approaches. The first is that there is a fairly large number of realistic 2BFs that reproduce the scattering data and the deuteron properties with equivalent high accuracy. Predictions employing these interactions are then compared to estimate the uncertainty due to the nuclear interaction (leading to an estimate of systematic uncertainties). The second source of uncertainties is due to the different methods employed to solve the nuclear many-body problem (which also contributes to systematical uncertainties). A critical comparison of various microscopic approaches using the same 2BF was performed in Ref. [19]. The aim of this work was to find the sources of discrepancies, and ultimately to determine a “systematic uncertainty” associated with the different microscopic approaches of the nuclear matter EOS. The approaches considered are the Brueckner–Hartree–Fock (BHF) [20,21], the Brueckner–Bethe–Goldstone (BBG) expansion up to the third order [22,23,24,25,26,27,28], the self-consistent Green’s function (SCGF) [29,30,31,32,33], the auxiliary field diffusion Monte Carlo (AFDMC) [34,35], the Green’s function Monte Carlo (GFMC) [36], and the Fermi hypernetted chain (FHNC) [37,38,39,40]. The properties of pure neutron matter (NM) and symmetric nuclear matter (SM) were computed with simplified versions [41] of the widely used Argonne Av18 potential [42], and a careful comparison of the results obtained with the different approaches was performed. The results of this work confirm that the tensor and spin–orbit components of the 2BF and their in-medium treatment are responsible for most of the observed discrepancies among these approaches [19]. A very similar study was performed recently in Refs. [43,44].
The EDF models adjusted to microscopic approaches shall incorporate these sources of uncertainties in the microscopic predictions for dense matter. Several authors have already determined EDF models from microscopic approaches (see, for example, Refs. [16,17,45,46,47,48,49,50,51,52,53,54,55]). The authors of Ref. [45], for instance, constructed a Skyrme-type EDF based on BHF results for infinite matter and applied it to a set of isospin-symmetric and -asymmetric nuclei through the mass table. A similar strategy to link the microscopic EOS in nuclear matter with finite nuclei properties was also considered in Refs. [46,47]. Their results provided a microscopic basis for a link between nuclear surface behavior and the NM EOS, as previously observed with phenomenological effective forces. The same authors also developed the so-called Barcelona–Catania–Paris–Madrid (BCPM) EDF [48,49,50,51] obtained from BHF calculations of nuclear matter within an approximation inspired by the Kohn–Sham formulation of density functional theory [56]. The BCPM EDF is built with a bulk part obtained directly from BHF results for NM and SM via a local-density approximation, and is supplemented by a phenomenological surface part together with the Coulomb, spin–orbit, and pairing contributions. The functional, with four adjustable parameters, is able to describe the ground-state properties of finite nuclei with an accuracy comparable to that of Skyrme and Gogny forces.
Later on, the so-called LNS Skyrme force was built in Ref. [52] based also on BHF calculations of infinite nuclear matter with consistent 2BFs and 3BFs. It was further refined in Ref. [54] to reproduce experimental binding energies and charge radii of some selected nuclei. A similar fitting procedure was also used in Refs. [16,17], where various Skyrme forces were adjusted to the ( S , T ) decomposition of the results of two different BHF calculations [57,58] 1. The comparison of the two microscopic calculations provides a kind of theoretical uncertainty based on the 2BFs and 3BFs employed, and suggests that the associated differences should be carefully analyzed. In Ref. [55], the tensor contribution to the ( S , T ) channels and different partial waves was analyzed. However, similarly to previous works, the uncertainty of the microscopic calculations was not included in the determination of the EDF. Finally, let us mention the analysis of Quantum Monte Carlo (QMC) calculations of spin-polarized neutron matter, which have been compared to predictions of several Skyrme EDFs [59]. Such analysis allowed to probe the time-odd part of the Skyrme EDFs independently of the time-even one and was found to be very constraining. Also, here an estimate of the uncertainty of the microscopic calculation would be very valuable, particularly as the density of the medium increases.
In the present paper, we aim to provide an estimation of the uncertainty of the nuclear EOS predicted by the BHF approach up to about 0.4 fm   3 2 n sat based on several nuclear two- and three-body interactions. Special attention is paid to the impact of the latter ones for which we consider two different types: phenomenological and microscopical 3BF. We then determine the corresponding parameter ranges of Skyrme-type nuclear energy density functionals, and more importantly, we show that the BHF calculations in the different ( S , T ) channels in SM and NM cannot be reproduced by the Skyrme EDF above n sat , already from the 2BF results.
The paper is organized as follows. In Section 2, a brief review of the BHF approach is presented. The ( S , T ) and partial-wave decomposition of the BHF correlation energy per particle for SM and NM is analyzed in Section 3. In Section 4, this decomposition is employed to constrain the parameters of Skyrme-type forces. A summary and the main conclusions are given in Section 5.

2. BHF Approach of Nuclear Matter

The BHF approach is the lowest-order realization of the BBG many-body theory of nuclear matter [60,61]. In this theory, the ground-state energy of nuclear matter is evaluated in terms of the so-called hole-line expansion, where the perturbative diagrams are grouped according to the number of independent hole-lines. The expansion is derived by means of the in-medium two-body scattering matrix G, which describes the effective interaction between two nucleons in the presence of a surrounding medium. The G matrix is obtained by solving the well-known Bethe–Goldstone equation
i j | G | i j = i j | V | i j + 1 Ω m l i j | V | m l Q m l ε i + ε j ε m ε l + i η m l | G | i j ,
where the multi-indices i , j , m , l , i , j indicate all the quantum numbers (momentum k, spin σ , and isospin τ projections) characterizing the two nucleons in the initial, intermediate, and final states; V denotes the bare N N interaction; Ω is the (large) volume enclosing the system; Q m l = ( 1 f m ) ( 1 f l ) with f i = θ ( k F ( τ i ) | k i | ) is the Pauli operator, taking into account the effect of the exclusion principle on the scattered holes; and k F ( τ i ) = ( 3 π 2 n τ i ) 1 / 3 is the Fermi momentum for particle i with density n τ i defined in the ground state. The constant η defines the width of the states (inverse of the lifetime). Since we consider the mean-field approximation, the lifetime of the states is infinite and η 0 . The single-particle (s.p.) energy of a nucleon with momentum k i and isospin τ i is
ε i ε τ i ( k i ) = k i 2 2 m τ i + R e U τ i ( k i ) ,
U τ i ( k i ) being the mean field felt by the nucleon i due to its interaction with the other nucleons of the medium. In the BHF approximation, it is calculated through the “on-shell” G-matrix,
U τ i ( k i ) = 1 Ω j f j i j | G | i j a ,
where the matrix elements are meant to be properly antisymmetrized. Note that Equations (1)–(3) should be solved self-consistently. The momentum dependence of the s.p. energy ε τ i can be characterized in terms of the effective mass [62],
m τ i ( k i ) m = k i m d ε τ i d k i 1 .
We consider in the following the isospin-asymmetric nuclear matter, and for clarity, we shall specify the isospin index τ . However, since we will discuss only spin-symmetric matter, the spin-up and -down potentials are equal and we do not specify the spin index in the following. For simplicity, we consider the same bare mass for neutrons and protons m τ = ( m n + m p ) / 2 m .
In the BHF approach, the total energy per particle of nuclear matter is given by the sum of only two-hole-line diagrams that include the effect of two-body correlations through the G matrix,
E tot = 1 n i f i k i 2 2 m + 1 2 n Ω i , j f i f j Re i j | G | i j a = 1 n i f i k i 2 2 m + E cor ,
where the first term on the right-hand side is the energy of a free Fermi gas and the second one is the so-called BHF correlation energy, E cor . We note here that the latter is usually referred to as BHF potential energy.
It has been shown in Refs. [22,23,24,25,26,27,28] for the V 18 , and recently confirmed in Ref. [63] for several modern 2BFs, that the contribution to the total energy from three-hole-line diagrams (which account for the effect of three-body correlations) is minimized when the so-called continuous prescription [62] is adopted for the BHF s.p. energy (2). This is a strong indication of the convergence of the hole-line expansion around saturation density and above. We adopt this prescription in our calculations and limit the exploration to densities less than 0.4 fm 3 , where the hole-line expansion parameter κ = ( c / d ) 3 (with c the interaction range and d the average distance between two nucleons) is still sufficiently small [63,64].
The present BHF calculations are carried out using a set of several phase-shift-equivalent N N potentials, namely, Av18 [39], NSC97a-f [65,66], a non-relativistic version of the Bonn B potential (Bonn) [67,68], and the charge-dependent Bonn potential (CD-Bonn) [69]. The predictions for the nuclear EOS based on these 2BFs are grouped together under the name 2BF. The group 2+3BFph combines the predictions of the same 2BF interactions supplemented with a 3BF of the Urbana type [70] consisting of the sum of the attractive two-pion-exchange Fujita–Miyazawa force with the excitation of an intermediate Δ -resonance [71] plus a phenomenological repulsive term. The group 2+3BFmic is a set of microscopic potentials [57] based on the Av18, Bonn B, or Nijmegen93 [72] potentials, which are also used to calculate consistently a three-body potential, including the virtual excitation of Δ (1232) and N (1440) resonances and anti-nucleons. In both cases, for the use in the BHF approach, these 3BF are reduced to an effective two-nucleon density-dependent force by averaging over the coordinates of the third nucleon [73,74,75]. The interested reader is referred to Refs. [57,58,76,77] for an extensive analysis of the use and effects of 3BF in nuclear matter.

3. Results of Brueckner–Hartree–Fock Microscopic Calculations

We discuss in this section the results of the BHF calculations of nuclear matter.

3.1. Correlation Energy Per Particle

In the BHF approach, the correlation energy per particle E cor , last term of Equation (5), can be easily linked to the spin–isospin and the partial-wave decomposition of the G-matrix, providing an interesting insight into the impact of the different ( S , T ) channels and partial waves contributing to the nuclear EOS. With the aim of providing more stringent constraints for phenomenological EDFs, we analyze in this section these decompositions for both SM and NM. A particular application of these constraints to a Skyrme force is discussed in the next section.

3.1.1. Spin–Isospin Decomposition of the Correlation Energy

The spin–isospin decomposition of the correlation energy per particle reads
E cor = S , T E cor ( S , T ) ,
where the contribution of each ( S , T ) channel,
E cor ( S , T ) = 1 2 n Ω i , j f i f j Re i j | G P S P T | i j a ,
can be obtained by using the projectors on spin and isospin singlet and triplet states, P S = 1 2 [ 1 + ( S 1 2 ) ( 1 + σ 1 · σ 2 ) ] and P T = 1 2 [ 1 + ( T 1 2 ) ( 1 + τ 1 · τ 2 ) ] . The partial waves contributing to the different channels are listed in Table 1.
The density dependence of the spin–isospin ( S , T ) decomposition of E cor in SM and NM is shown in Figure 1, where the bands correspond to the uncertainty associated with the different models used in our BHF calculation as previously discussed. In SM, the ( 0 , 0 ) channel is repulsive, while the ( 1 , 0 ) and ( 0 , 1 ) channels are attractive. The ( 1 , 1 ) channel is much weaker than the others and could be attractive or repulsive depending on the 3BF interactions (it is slightly repulsive with 2BF and 2+3BFph). The 2+3BFmic approach has a density dependence strong enough to change the attractive ( 0 , 1 ) channel for 2BF into a repulsive one at high density, or to turn the repulsive ( 1 , 1 ) channel into an attractive one around n sat . The ( 0 , 1 ) and ( 1 , 1 ) channels of NM are qualitatively similar to those of SM. The impact of the 3BF prescription can thus be opposite to the 2BF in the various channels. These features can be traced back to the general characteristics of the 3BF: the two-pion-exchange 3BF component is attractive in the ( 1 , 0 ) channel and repulsive in the ( 0 , 1 ) channel, and scalar repulsive 3BF components become dominant at high-density in all channels.

3.1.2. Partial-Wave Decomposition of the Correlation Energy

Let us now further continue this analysis by looking in more detail at the partial-wave decomposition of E cor for both SM and NM. The different partial waves contribute to the correlation energy as
E cor = L S J E cor ( L S J ) ,
where
E cor ( L S J ) = ( 2 J + 1 ) 4 n ( 2 π ) 6 τ τ d 3 k τ f τ ( k τ ) d 3 k τ f τ ( k τ ) f τ τ K q L S J | G τ τ τ τ | K q L S J
with K = | k τ + k τ | , q = | ( k τ k τ ) | / 2 and f τ τ = 1 + ( 1 ) L + S δ τ τ .
We first show in Figure 2 (for SM) and Figure 3 (for NM) the decomposition into partial waves with fixed angular momentum L, E cor ( L ) S J E cor ( L S J ) , i.e., all partial waves with the same value of L are summed together. Contributions up to a total angular momentum J = 8 are considered in the calculation, although only results for the S , P , D , and F partial waves are shown since the contributions of higher partial waves are small. It is indeed clear from the figures that the contributions of the different E cor ( L ) decrease as L increases for the range of densities considered here. Note that the P and D channels are found to be of the same order of magnitude in absolute value. Since the standard Skyrme EDF contributes only to the S and P channels, the relatively large importance of the D channel advocates for the extension of the standard Skyrme EDF as suggested in Ref. [55]. The F channel is found to be repulsive and less important than the others for 2BF and 2+3BFph but strongly attractive and comparable to the D channel for 2+3BFmic. This corresponds to the ( 0 , 0 ) and ( 1 , 1 ) channels in Figure 1.
The most significant difference between SM and NM is in the S channel, which is less attractive in NM compared to SM, due to the fact that the dominant attractive   3 S 1 channel contributes only to SM. The other channels are qualitatively similar in SM and NM. At low density n < n sat , one can clearly observe a dominance of the L = 0 contribution compared to the L > 0 terms, but as the density increases, those terms become increasingly important. As for the ( S , T ) decomposition, it is observed that the dispersion between the different predictions of the partial waves is generally larger in SM compared to NM.
Finally, we present individual partial wave contributions to E cor in Figure 4, for the   1 S 0 ,   1 P 1 ,   3 P 0 ,   3 P 1 ,   3 S D 1 and   3 P F 2 waves. It is only in the   1 S 0 channel that the 3BF are similar, while in the other channels, they act in an opposite way compared to the 2BF results. The coupled   3 S D 1 and   3 P F 2 waves are predicted to be attractive in all cases, caused by the tensor force. Note that the standard Skyrme EDFs contribute only to the S and P channels, while in Ref. [55], the partial waves   3 P 0 and   3 P 1 are used to calibrate the spin–orbit coupling W 0 and the tensor parameter t o ( n ) of extended EDFs.

3.2. BHF Single-Particle Energy

In the BHF approach, the s.p. energy ε τ ( k ) is defined by Equation (2). In some phenomenological models, however, like the Skyrme one, the momentum dependence of ε τ is simply quadratic and can be entirely incorporated in a modification of the mass appearing in the kinetic term,
ε τ ( k ) = k 2 2 m τ + u τ ,
where the effective mass m τ is a momentum-independent quantity and the s.p. potential u τ is the value of the s.p. potential at k = 0 : u τ = Re [ U τ ( k = 0 ) ] . Since we want to establish a link between microscopical and phenomenological approaches, using the BHF mean field potential U τ ( k ) , one can define a global effective mass:
m τ m = 1 + 2 m k F τ 2 Re [ U τ ( k F τ ) U τ ( 0 ) ] 1 .
Note that a local momentum-independent effective mass can also be determined from Equation (4) as m τ = m τ ( k = k F τ ) . This latter definition provides effective masses, which are typically 5–10% higher than using Equation (11) at saturation density, due to the presence of a wiggle in the BHF s.p. potential close to k F [78]. For the comparison with phenomenological models, we use Equation (11) for the global effective mass in the following.
The dependence of the effective mass on the isospin asymmetry is also an interesting result to extract from the BHF approach. One usually defines the effective-mass isospin splitting as
Δ β ( n ) m n ( n , β ) m p ( n , β ) m ,
where n = n n + n p is the total density and β = ( n n n p ) / n is the isospin asymmetry parameter. Δ β ( n ) is often defined for maximal asymmetry β = 1 (NM) and at saturation density n = n sat . Since this choice requires some technicalities for the definition of the proton effective mass in a pure neutron environment, we prefer in this work to employ Δ 1 / 2 ( n sat ) , using the approximation Δ 1 2 Δ 1 / 2 [79]. Note that this relation is exact for Skyrme-type interactions; see Section 4.
The density dependence of the effective mass in symmetric matter m ( n ) / m and of the isospin splitting Δ 1 / 2 ( n ) is shown in Figure 5(left,right), respectively. At normal density, the impact of 3BFs is small for both the effective mass and the isospin splitting. The effective mass is already well determined by 2BF only, which predicts a continuous decrease in the effective mass and a positive sign for Δ 1 / 2 ( n ) 0.03 ± 0.03 above n sat . While the differences between 2+3BFph and 2+3BFmic predictions are small below the saturation density, at larger density, the 2+3BFmic results remain close to the 2BF ones, while 2+3BFph predict a rise in the effective mass and a change in the sign of Δ 1 / 2 above 2 n sat .
We therefore conclude that the momentum dependence of the 3BFph results is more marked than the one for 3BFmic. In the next subsection, we show that this is due to compensation between two ( S , T ) contributions for 3BFmic, which is absent in the case of 3BFph. As a consequence, as the density increases, 3BFph tends to reduce the correction due to the effective mass and the difference between the momentum-dependent contributions from neutrons and protons in asymmetric matter.

Spin–Isospin Decomposition of the Effective Mass

In Figure 6, we show the spin–isospin decomposition of the inverse of the effective mass
m m τ = 1 + S , T m m τ ( S , T ) ,
where
m m τ ( S , T ) 2 m k F τ 2 Re U τ ( S , T ) ( k F τ ) U τ ( S , T ) ( 0 ) .
U τ ( S , T ) ( k ) being the contribution of the spin–isospin channel ( S , T ) to the s.p. mean field.
Figure 6. Spin−isospin decomposition of the effective mass ( m / m ) ( S , T ) , Equation (14), in SM and NM as a function of the density for various interactions.
Figure 6. Spin−isospin decomposition of the effective mass ( m / m ) ( S , T ) , Equation (14), in SM and NM as a function of the density for various interactions.
Universe 10 00226 g006
In the considered density range, the inverse of the effective mass is dominated by the S-wave contributions in the ( 0 , 1 ) and ( 1 , 0 ) channels, where there is a reasonable agreement between the two 3BF prescriptions. Note that in the SM ( 0 , 1 ) channel, the two 3BF predict an increase in the inverse of the effective mass compared to the 2BF above n sat : this originates from an attractive term in the 3BF, which moderates the repulsive contribution of the 2BF.
In the SM ( 1 , 1 ) (and NM ( S = 1 ) ) channels, the 2+3BFmic has a stronger density dependence than 2+3BFph. The larger contribution to m / m indicates that there is an additional repulsive momentum dependence in 3BFmic that is absent from 3BFph. For 3BFmic, it happens that this additional repulsive contribution in the ( 1 , 1 ) channel compensates for the additional attraction in ( 0 , 1 ) , such that the sum of all contributions to m / m is similar for 2BF and 2+3BFmic as remarked in the discussion of Figure 5. For 3BFph, the absence of the repulsive contribution in the ( 1 , 1 ) channel induces the marked difference between 2BF and 2+3BFph for the inverse of the effective mass shown in Figure 5. In the following, we show that the more repulsive 3BFmic leads also to more repulsion for the total energy per particle.

3.3. Total Energy

We show in Figure 7 the correlation and total energy per particle of SM and NM as a function of the density for the various 2BFs with and without the contribution of 3BFs. The SM results with only 2BF exhibit a saturation density well above the empirical value n sat . This is well known and due to the different off-shell behavior of the 2BF: the phase-shift-equivalent potentials lead to different saturation curves, with the corresponding saturation points lying in the so-called “Coester band“ [80]. The band defined by these curves can be used to estimate the contribution of 2BF to the uncertainty in the energy per particle. As expected, the error bar is small at low densities, where the interaction is mainly dominated by the S partial waves, but increases rapidly with density, where higher partial waves become important. Similar behavior is observed for the 3BFph and 3BFmic results, and also for NM, where the uncertainty bands are narrower due to the weaker interaction. Above n sat , we observe that 2+3BFmic predict more repulsion than 2+3BFph as discussed in the previous subsection.
It is also interesting to note that the 2+3BFmic approach has much wider error bars in the ( S , T ) decomposition of E cor shown in Figure 1 than in the total E cor shown in Figure 7 due to strong compensation between the different channels, which makes the total binding energy less dispersed than the ( S , T ) components.
At low density, our predictions are compared to χ -EFT predictions based on 2+3BF calculated within different orders and approaches:
(a)
N2LO Wlazłowski QMC 2014 [81]: Variational QMC calculation of NM with chiral nuclear forces at N2LO for 2BF and 3BF.
(b)
N2LO Tews QMC 2016 [82]: AFDMC calculation of NM with chiral nuclear forces at N2LO for 2BF. The 3BF is taken at leading order and local.
(c)
N3LO Drischler MBPT 2016 [83]: Many-body perturbation theory approach based on regularized chiral nuclear forces at N3LO for SM and NM.
Note that while there are no uncertainties estimated in (a), the estimation of the uncertainty is performed differently in (b) and (c) compared to our BHF results. In those cases, the uncertainty includes the unknown from the chiral 2BF fitted on phase shifts as well as the uncertainty in the truncation order, which is expected to decrease as the order increases. Thus, it is not surprising that the N3LO calculation [83] has slightly smaller error bars than the N2LO one [82]. There are also differences in the way the band widths are estimated among these calculations. It is, however, interesting to note that the band of uncertainty in our case is compatible with the most advanced N3LO MBPT calculation [83] in both SM and NM. This provides a solid ground for our extrapolation to higher density, even if one cannot exclude an additional uncertainty related to the theoretical approach itself.
The conclusion of this part is that there are differences between the 3BF considered in this paper, 3BFph and 3BFmic, due to their different construction. The 3BFmic predicts a more repulsive energy per particle than the 3BFph at high density.

4. Constraining the Skyrme Energy Density Functional from BHF Calculations

The commonly called standard Skyrme interaction is a contact interaction with momentum-dependent terms, which is expressed in the following form:
V ( r 1 , r 2 ) = t 0 ( 1 + x 0 P σ ) δ ( r ) + t 1 2 ( 1 + x 1 P σ ) k 2 δ ( r ) + δ ( r ) k 2 + t 2 ( 1 + x 2 P σ ) k · δ ( r ) k + t 3 6 ( 1 + x 3 P σ ) n ( R ) γ δ ( r ) + i W 0 ( σ 1 + σ 2 ) · [ k × δ ( r ) k ] ,
where r = r 1 r 2 , R = ( r 1 + r 2 ) / 2 , k = ( 1 2 ) / i is the relative momentum acting on the right, k is its conjugate acting on the left, and P σ = ( 1 + σ 1 · σ 2 ) / 2 is the spin exchange operator. The last term, proportional to W 0 , corresponds to the zero-range spin–orbit term. It does not contribute to the EOS in homogeneous systems and thus will be ignored for the rest of this article. In its standard form (15), the Skyrme interaction contributes only to the L = 0 and L = 1 channels, which implies some limitations, which will appear in the following discussion.
In the Hartree–Fock approximation, the total energy density of SM ( β = 0 ) and NM ( β = 1 ) is given respectively as
ε ( n , β = 0 ) = 1 2 m τ SM + C 0 τ n τ SM + C 0 n ( n ) n 2 ,
ε ( n , β = 1 ) = 1 2 m τ NM + [ C 0 τ + C 1 τ ] n τ NM + C 0 n ( n ) + C 1 n ( n ) n 2 ,
where the kinetic densities are τ SM = 3 5 ( 3 π 2 2 ) 2 / 3 n 5 / 3 and τ NM = 3 5 ( 3 π 2 ) 2 / 3 n 5 / 3 , and C t n ( n ) = C t 0 n + C t 3 n n γ for t = 0 , 1 . In this work, we employ the DFT coefficients, e.g., C t τ , C t 0 n , and C t 3 n , which can be defined in terms of the Skyrme parameters t i and x i ; see Refs. [84,85,86], where they were introduced.
Several authors (see, for example, Refs. [16,17,52,53,54]) have stressed the interest of using the spin−isospin decomposition of the BHF correlation energy as an additional constraint to better determine the parameters of energy density functionals based on effective interactions such as the Skyrme one. It indeed allows a more detailed determination of the parameters entering into the calculation of the energy per particle by increasing the number of equations to solve. From Equation (7), one can calculate the decomposition of the correlation energy density in terms of SE ε cor ( 0 , 1 ) , TE ε cor ( 1 , 0 ) , SO ε cor ( 0 , 0 ) , and TO ε cor ( 1 , 1 ) contributions as [53]
ε cor ( 0 , 0 ) = 1 16 C 0 τ 3 C 0 s T 3 C 1 τ + 9 C 1 s T n τ SM , ε cor ( 0 , 1 ) = 3 4 C 0 n ( n ) + C 1 n ( n ) n 2 + 3 16 C 0 τ 3 C 0 s T + C 1 τ 3 C 1 s T n τ SM , ε cor ( 1 , 0 ) = 1 4 C 0 n ( n ) 3 C 1 n ( n ) n 2 + 3 16 C 0 τ + C 0 s T 3 C 1 τ 3 C 1 s T n τ SM , ε cor ( 1 , 1 ) = 9 16 C 0 τ + C 0 s T + C 1 τ + C 1 s T n τ SM
in SM, and as
ε cor ( 0 ) = C 0 n ( n ) + C 1 n ( n ) n 2 + 1 4 C 0 τ 3 C 0 s T + C 1 τ 3 C 1 s T n τ NM ε cor ( 1 ) = 3 4 C 0 τ + C 0 s T + C 1 τ + C 1 s T n τ NM
in NM. The correlation energy in the channels ( 0 , 0 ) and ( 1 , 1 ) in SM, and ( 1 ) in NM scales only with the kinetic energy, τ SM and τ NM , because only the angular momentum L = 1 can contribute in these channels. In the channels ( 0 , 1 ) and ( 1 , 0 ) in SM, and ( 0 ) in NM, there is also a term that scales with the density n reflecting the L = 0 contribution of the Skyrme interaction. Depending on the sign of the combination of the coefficients, the energy per particle in the ( S , T ) channels can be positive or negative in the channels ( 0 , 0 ) and ( 1 , 1 ) in SM, and ( 1 ) in NM, or can change its sign as a function of the density in the other channels. This is a clear limitation of the standard Skyrme interaction (15) since BHF calculations predict a richer density dependence; see Figure 1, for instance.
In nuclear matter, the Skyrme interaction is determined by eight coefficients C t τ , C t s T , C t 0 ρ , and C t 3 ρ , as well as the exponent of the density dependence γ . We have six constraints in terms of the ( S , T ) decomposition and we decide to fix the value of the parameter γ , see hereafter. We therefore need two additional constraints to unambiguously determine all Skyrme parameters. These two other constraints will be given by the effective mass in SM and the isospin splitting, namely, from the BHF predictions, we fix m sat / m = 0.7 and Δ 1 , sat = 0.1 for 2BF, 2+3BFph, and 2+3BFmic.
For the Skyrme interaction, the nucleon effective mass is expressed as [53,84,85,86]
1 2 m τ ( n , β ) = 1 2 m + C 0 τ + τ 3 C 1 τ β n .
The parameters C 0 τ and C 1 τ can therefore be obtained as
C 0 τ = 1 n sat 1 2 m m m sat 1 ,
C 1 τ = 1 n sat 1 2 m 1 β Δ β 1 1 + Δ β m sat / m 2
1 n sat 1 2 m Δ β 2 β m m sat 2 ,
where the last approximation is valid for typical small values of Δ β m sat / m 2 0.02 .
We determine the eight Skyrme coefficients by imposing the global reproduction of the six ( S , T ) channels predicted by the BHF calculations in the density range going from 0.01 up to 0.4 fm   3 in SM and NM (including uncertainties), complemented by the effective mass m sat / m and its isospin splitting Δ 1 , sat . The fit is provided by the nonlinear least-squares Marquardt–Levenberg algorithm encoded in Gnuplot 2. The results are shown in Figure 8, where the parameter γ is explored from 1 to 1 / 10 . The error bars are estimated from the covariant matrix properties at the solution. There are four Skyrme coefficients C t τ and C t s T , with t = 0 , 1 , which are independent of γ , and four other coefficients C t 0 ρ and C t 3 ρ , which depend on the value taken for γ . We also observe a strong impact of the 3BF for the coefficient C 1 s T , and a noticeable but less strong impact for the coefficients C t 0 ρ and C t 3 ρ .
In order to evaluate the quality of the fits, we compare in Figure 9 the quantities targeted by our fits and the results provided by our fits. The Skyrme spin–isospin ( S , T ) channel decomposition of the correlation energy per particle in SM and NM are shown together with the constraints extracted from BHF calculations. We consider the centroid of the Skyrme coefficients obtained from Figure 8 and vary the parameter γ . The effect of γ is very weak since it is absorbed in the values of the coefficients C t 0 n and C t 3 n ; see Figure 8. The other four coefficients C t τ and C t s T are independent of γ . The quality of the fit is, however, not very good. This is particularly clear in the following channels: ( 0 , 0 ) , ( 1 , 1 ) in SM and ( 1 ) in NM. These channels are the ones where only the L = 1 contribution of the Skyrme interaction plays a role. The poor reproduction of the BHF results is therefore due to the lack of flexibility of the standard form (15) of the Skyrme interaction, where higher values of L are necessary.
Let us remark from Figure 9 that the 2BF is quite well reproduced by the Skyrme model in SM and for the ( 0 , 0 ) , ( 0 , 1 ) and ( 1 , 0 ) channels. The ( 1 , 1 ) channel has a wrong sign in the Skyrme model, although the strength of the interaction in this channel is quite weak. In NM, however, the fit of the 2BF is poor above n sat . We, therefore, conclude that already at the 2BF level, the Skyrme interaction cannot well reproduce the BHF predictions in both SM and NM above n sat .
One question that, however, remains—and which goes beyond the scope of the present study—is to what extent does the mismatch between the microscopic BHF predictions for SM and NM and the Skyrme functional impact the properties of finite nuclei? The present study cannot answer this question, which implies the study of a large number of finite nuclei properties. We can, however, refer to the existing studies focusing on a few properties, e.g., nuclear masses and charge radii.
In Refs. [16,17], a set of Skyrme interactions was obtained by employing the ( S , T ) decomposition in SM for two BHF calculations, and the NM channel was controlled globally by variational calculations. The difference with our approach is (i) that our uncertainty estimate is based on a larger number of microscopic predictions, and (ii) that we employ the ( S ) channel decomposition in NM to adjust the Skyrme EDF. In Ref. [54], only one BHF calculation is employed to fit the Skyrme EDF, and only the ( S , T ) decomposition in SM is employed. In Refs. [16,17,54], data from finite nuclei, e.g., binding energies and charge radii, are employed in addition to the BHF constraints. It is, therefore, difficult to perform detailed comparisons of our results with the ones presented in these papers. We, however, conclude from our analysis that the standard Skyrme EDF should be enriched, with higher L terms, for instance, in order to reproduce BHF calculations in the ( S , T ) channel in SM and NM up to at least 0.4 fm   3 . The results presented in this paper can be employed for such an aim.

5. Conclusions

In this paper, we analyzed the BHF results of SM and NM by decomposing the total energy and the effective mass into spin−isospin ( S , T ) channels and also into contributions from various partial waves. The two-body and three-body contributions were shown separately as bands representing the uncertainties after exploring a set of interactions. We found that 3BFmic induces more repulsion than 3BFph, visible in the energy per particle at high density. This repulsion originates from a stronger repulsive scalar meson exchange term of 3BFmic compared to 3BFph. Since the equation of state at high density is dominated by 3BF, we expect that the qualitative difference between 3BFph and 3BFmic that we discuss here has a noticeable impact: 3BFmic induces a stronger repulsion than 3BFph. The effective mass, however, seems to be more impacted by 3BFph than by 3BFmic. We show that this results from a compensation between two channels in 3BFmic, which is absent for 3BFph.
We then addressed the question of the effectiveness of the standard Skyrme model (15) to reproduce the ( S , T ) decompositions predicted by BHF calculations based on different 2BFs and 3BFs. Our results indicate that the standard Skyrme model cannot reproduce the ( S , T ) decomposition predicted by BHF calculations above n sat considering only 2BF. This conclusion is even stronger when 3BFs are added.
The detailed and numerous quantities obtained from BHF calculations and presented in this paper can, however, serve as guidelines for improvements in the Skyrme interaction with close links to BHF calculations. Having reliable models to describe the properties of baryonic matter at extreme densities, on the one hand, and fair estimates of the theoretical uncertainties, on the other hand, are important for the treatment of the astrophysical signals that will be collected in the future.

Author Contributions

All authors contributed to this work. All authors have read and agreed to the published version of the manuscript.

Funding

J.M. is partially supported by CNRS-IN2P3 MAC masterproject, as well as from the LABEX Lyon Institute of Origins (ANR-10-LABX-0066).

Data Availability Statement

Data shown in this paper are available on demand to the authors.

Acknowledgments

This work was also fostered by the scientific emulation created by the “NewCompStar” COST Action MP1304 and PHAROS COST Action MP16214.

Conflicts of Interest

The authors declare no conflicts of interest.

Notes

1
Note that in Figure 7 of Refs. [16,17], the BHF predictions for the ( S , T ) channels ( 0 , 1 ) and ( 1 , 1 ) in SM from catania2 [57] results are exchanged.
2

References

  1. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya V., B.; et al. GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef]
  2. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; Adya, V.B.; et al. Multi-messenger Observations of a Binary Neutron Star Merger. Astrophys. J. Lett. 2017, 848, L13. [Google Scholar]
  3. Weinberg, S. Phenomenological Lagrangians. Physica 1979, A96, 327. [Google Scholar] [CrossRef]
  4. Weinberg, S. Nuclear forces from chiral lagrangians. Phys. Lett. 1990, B251, 288. [Google Scholar] [CrossRef]
  5. Weinberg, S. Effective chiral lagrangians for nucleon-pion interactions and nuclear forces. Nucl. Phys. B 1991, 363, 3. [Google Scholar] [CrossRef]
  6. Epelbaum, E.; Hammer, H.-W.; Meißner, U.-G. Modern theory of nuclear forces. Rev. Mod. Phys. 2009, 81, 1773. [Google Scholar] [CrossRef]
  7. Machleidt, R.; Entem, D.R. Chiral effective field theory and nuclear forces. Phys. Rep. 2011, 503, 1. [Google Scholar] [CrossRef]
  8. Margueron, J.; Casali, R.H.; Gulminelli, F. Equation of state for dense nucleonic matter from metamodeling. I. Foundational aspects. Phys. Rev. 2018, C97, 025805. [Google Scholar] [CrossRef]
  9. Hebeler, K.; Lattimer, J.M.; Pethick, C.J.; Schwenk, A. Equation of State and Neutron Star Properties Constrained By Nuclear Physics and Observation. Astrophys. J. 2013, 773, 11. [Google Scholar] [CrossRef]
  10. Tews, I.; Carlson, J.; Gandolfi, S.; Reddy, S. Constraining the Speed of Sound inside Neutron Stars with Chiral Effective Field Theory Interactions and Observations. Astrophys. J. 2018, 860, 149. [Google Scholar] [CrossRef]
  11. Drischler, C.; Hebeler, K.; Schwenk, A. Chiral Interactions up to Next-to-Next-to-Next-to-Leading Order and Nuclear Saturation. Phys. Rev. Lett. 2019, 2019, 042501. [Google Scholar] [CrossRef] [PubMed]
  12. Bogner, S.K.; Kuo, T.T.S.; Coraggio, L. Low momentum nucleon-nucleon potentials with half-on-shell T-matrix equivalence. Nucl. Phys. 2001, A 684, 432c. [Google Scholar]
  13. Bogner, S.K.; Kuo, T.T.S.; Schwenk, A. Model-independent low momentum nucleon interaction from phase shift equivalence. Phys. Rep. 2003, 386, 1. [Google Scholar] [CrossRef]
  14. Chabanat, E.; Bonche, P.; Haensel, P.; Meyer, J.; Schaeffer, R. A Skyrme parametrization from subnuclear to neutron star densities. Nucl. Phys. A 1997, 627, 710. [Google Scholar] [CrossRef]
  15. Chabanat, E.; Bonche, P.; Haensel, P.; Meyer, J.; Schaeffer, R. A Skyrme parametrization from subnuclear to neutron star densities Part II. Nuclei far from stabilities. Nucl. Phys. A 1998, 635, 231. [Google Scholar] [CrossRef]
  16. Goriely, S.; Chamel, N.; Pearson, J.M. Further explorations of Skyrme-Hartree-Fock-Bogoliubov mass formulas. XII. Stiffness and stability of neutron-star matter. Phys. Rev. C 2010, 82, 035804. [Google Scholar] [CrossRef]
  17. Goriely, S.; Chamel, N.; Pearson, J.M. Hartree-Fock-Bogoliubov nuclear mass model with 0.50 MeV accuracy based on standard forms of Skyrme and pairing functionals. Phys. Rev. C 2013, 88, 024308. [Google Scholar] [CrossRef]
  18. Fantina, A.; Chamel, N.; Pearson, J.M.; Goriely, S. Neutron star properties with unified equations of state of dense matter. Astron. Astrophys. 2013, 559, A128. [Google Scholar] [CrossRef]
  19. Baldo, M.; Polls, A.; Rios, A.; Schulze, H.-J.; Vidaña, I. Comparative study of neutron and nuclear matter with simplified Argonne nucleon-nucleon potentials. Phys. Rev. 2012, C86, 064001. [Google Scholar] [CrossRef]
  20. Baldo, M. The Many-Body theory of the nuclear equation of state. In Nuclear Methods and the Nuclear Equation of State; International Review of Nuclear Physics; World Scientific: Singapore, 1999; Volume 8. [Google Scholar]
  21. Schulze, H.-J.; Cugnon, J.; Lejeune, A.; Baldo, M.; Lombardo, U. Status of the Brueckner-Hartree-Fock approximation to the nuclear matter binding energy with the Paris potential. Phys. Rev. C 1995, 52, 2785. [Google Scholar] [CrossRef] [PubMed]
  22. Song, H.Q.; Baldo, M.; Giansiracusa, G.; Lombardo, U. Bethe-Brueckner-Goldstone expansion in nuclear matter. Phys. Rev. Lett. 1998, 81, 1584. [Google Scholar] [CrossRef]
  23. Baldo, M.; Giansiracusa, G.; Lombardo, U.; Song, H.Q. Bethe-Brueckner-Goldstone expansion in neutron matter. Phys. Lett. 2000, B473, 1. [Google Scholar] [CrossRef]
  24. Baldo, M.; Fiasconaro, A.; Song, H.Q.; Giansiracusa, G.; Lombardo, U. High density symmetric nuclear matter in the BBG approach. Phys. Rev. 2001, C65, 017303. [Google Scholar]
  25. Sartor, R. Solution of the Bethe-Faddeev equations within the continuous version of the hole-line expansion. Phys. Rev. 2006, C73, 034307. [Google Scholar] [CrossRef]
  26. Baldo, M.; Maieron, C. Equation of state of nuclear matter at high baryon density. J. Phys. 2007, G34, R243. [Google Scholar] [CrossRef]
  27. Baldo, M.; Maieron, C. Nonlocality in the NN interaction and nuclear matter saturation. Phys. Rev. 2005, C72, 034005. [Google Scholar]
  28. Fukukawa, K.; Baldo, M.; Burgio, G.F.; Monaco, L.L.; Schulze, H.-J. Nuclear matter EOS from a Quark-Model nucleon-nucleon interaction. Phys. Rev. 2015, C92, 065802. [Google Scholar]
  29. Dickhoff, W.H.; Barbieri, C. Self-consistent Green’s function method for nuclei and nuclear matter. Prog. Part. Nucl. Phys. 2004, 52, 377. [Google Scholar] [CrossRef]
  30. Frick, T.; Müther, H.; Rios, A.; Polls, A.; Ramos, A. Correlated nucleons in configuration space. Phys. Rev. C 2005, 71, 014313. [Google Scholar] [CrossRef]
  31. Somà, V.; Bozèk, P. Diagrammatic calculation of thermodynamical quantities in nuclear matter. Phys. Rev. C 2006, 74, 045809. [Google Scholar] [CrossRef]
  32. Somà, V.; Bozèk, P. In-medium T matrix for nuclear matter with three-body forces: Binding energy and single-particle properties. Phys. Rev. C 2008, 78, 054003. [Google Scholar] [CrossRef]
  33. Rios, A.; Polls, A.; Vidan ña, I. Hot neutron matter from a self-consistent Green’s-functions approach. Phys. Rev. C 2009, 79, 025802. [Google Scholar] [CrossRef]
  34. Gandolfi, S.; Pederiva, F.; Fantoni, S.; Schmidt, K.E. Quantum Monte Carlo Calculations of Symmetric Nuclear Matter. Phys. Rev. Lett. 2007, 98, 102503. [Google Scholar] [CrossRef]
  35. Gandolfi, S.; Illarionov, A.Y.; Schmidt, K.E.; Pederiva, F.; Fantoni, S. Quantum Monte Carlo calculation of the equation of state of neutron matter. Phys. Rev. C 2009, 79, 054005. [Google Scholar] [CrossRef]
  36. Carlson, J.; Morales, J.; Pandharipande, V.R.; Ravenhall, D.G. Quantum Monte Carlo calculations of neutron matter. Phys. Rev. C 2003, 68, 025802. [Google Scholar] [CrossRef]
  37. Pandharipande, V.R.; Wiringa, R.B. Variations on a theme of nuclear matter. Rev. Mod. Phys. 1979, 51, 821. [Google Scholar] [CrossRef]
  38. Fantoni, S.; Fabrocini, A. Correlated basis function theory for fermion systems. In Microscopic Quantum Many- Body Theories and Their Applications; Lecture Notes in Physics; Navarro, J., Polls, A., Eds.; Springer-Verlag: Berlin, Germany, 1998; Volume 510. [Google Scholar]
  39. Lovato, A.; Benhar, O.; Fantoni, S.; Illarionov, A.Y.; Schmidt, K.E. Density-dependent nucleon-nucleon interaction from three-nucleon forces. Phys. Rev. C 2011, 83, 054003. [Google Scholar] [CrossRef]
  40. Wiringa, R.B.; Fiks, V.; Fabrocini, A. Equations of state for dense nucleon matter. Phys. Rev. C 1988, 38, 1010. [Google Scholar] [CrossRef] [PubMed]
  41. Wiringa, R.B.; Pieper, S.C. Evolution of Nuclear Spectra with Nuclear Forces. Phys. Rev. Lett. 2002, 89, 182501. [Google Scholar] [CrossRef]
  42. Wiringa, R.B.; Stoks, V.G.J.; Schiavilla, R. Accurate nucleon-nucleon potential with charge-independance breaking. Phys. Rev. C 1995, 51, 38. [Google Scholar] [CrossRef] [PubMed]
  43. Piarulli, M.; Bombaci, I.; Logoteta, D.; Lovato, A.; Wiringa, R.B. Benchmark calculations of pure neutron matter with realistic nucleon-nucleon interactions. Phys. Rev. C 2020, 101, 045801. [Google Scholar] [CrossRef]
  44. Lovato, A.; Bombaci, I.; Logoteta, D.; Piarulli, M.; Wiringa, R.B. Benchmark calculations of infinite neutron matter with realistic two- and three-nucleon potentials. Phys. Rev. C 2022, 105, 055808. [Google Scholar] [CrossRef]
  45. Baldo, M.; Maieron, C.; Schuck, P.; Viñas, X. Low densities in nuclear and neutron matters and in the nuclear surface. Nucl. Phys. A 2004, 736, 241. [Google Scholar] [CrossRef]
  46. Fayans, S.A. Towards a universal nuclear density functional. JETP Lett. 1998, 68, 169. [Google Scholar] [CrossRef]
  47. Steiner, A.W.; Prakash, M.; Lattimer, J.M.; Ellis, P.J. Isospin asymmetry in nuclei and neutron stars. Phys. Rep. 2005, 411, 325. [Google Scholar] [CrossRef]
  48. Baldo, M.; Schuck, P.; Viñas, X. Kohn–Sham density functional inspired approach to nuclear binding. Phys. Lett. B 2008, 663, 390. [Google Scholar] [CrossRef]
  49. Baldo, M.; Robledo, L.M.; Schuck, P.; Viñas, X. Energy density functional on a microscopic basis. J. Phys. G Nucl. Phys. 2010, 37, 064015. [Google Scholar] [CrossRef]
  50. Baldo, M.; Robledo, L.M.; Schuck, P.; Viñas, X. New Kohn-Sham density functional based on microscopic nuclear and neutron matter equations of state. Phys. Rev. C 2013, 87, 064305. [Google Scholar] [CrossRef]
  51. Baldo, M.; Robledo, L.M.; Schuck, P.; Viñas, X. Barcelona-Catania-Paris-Madrid functional with a realistic effective mass. Phys. Rev. C 2017, 95, 014318. [Google Scholar] [CrossRef]
  52. Cao, L.G.; Lombardo, U.; Shen, C.W.; Giai, N.V. From Brueckner approach to Skyrme-type energy density functional. Phys. Rev. C 2006, 73, 014313. [Google Scholar] [CrossRef]
  53. Lesinski, T.; Bennaceur, K.; Duguet, T.; Meyer, J. Isovector splitting of nucleon effective masses, ab initio benchmarks and extended stability criteria for Skyrme energy functionals. Phys. Rev. C 2006, 74, 044315. [Google Scholar] [CrossRef]
  54. Gambacurta, D.; Li, L.; Colò, G.; Lombardo, U.; Giai, N.V.; Zuo, W. Determination of local energy density functionals from BHF calculations. Phys. Rev. C 2011, 84, 024301. [Google Scholar] [CrossRef]
  55. Davesne, D.; Navarro, J.; Becker, P.; Jodon, R.; Meyer, J.; Pastore, A. Extended Skyrme pseudopotential deduced from infinite nuclear matter properties. Phys. Rev. C 2015, 91, 064303. [Google Scholar] [CrossRef]
  56. Khon, W.; Sham, L.J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 1965, 140, 1133. [Google Scholar] [CrossRef]
  57. Zhou, X.R.; Burgio, G.F.; Lombardo, U.; Schulze, H.-J.; Zuo, W. Three-body forces and neutron star structure. Phys. Rev. 2004, C69, 018801. [Google Scholar] [CrossRef]
  58. Li, Z.H.; Schulze, H.-J. Neutron star structure with modern nucleonic three-body forces. Phys. Rev. 2008, C78, 028801. [Google Scholar] [CrossRef]
  59. Roggero, A.; Mukherjee, A.; Pederiva, F. Constraining the Skyrme energy density functional with quantum Monte Carlo calculations. Phys. Rev. C 2015, 92, 054303. [Google Scholar] [CrossRef]
  60. Day, B.D. Elements of the Brueckner-Goldstone theory of nuclear matter. Rev. Mod. Phys. 1967, 39, 719. [Google Scholar] [CrossRef]
  61. Baldo, M. Nuclear Methods and the Nuclear Equation of State; World Scientific: Singapore, 1999. [Google Scholar]
  62. Jekenne, J.P.; Lejeunne, A.; Mahaux, C. Many-body theory of nuclear matter. Phys. Rep. 1976, 25, 83. [Google Scholar] [CrossRef]
  63. Lu, J.J.; Li, Z.H.; Chen, C.Y.; Baldo, M.; Schulze, H.-J. Convergence of the hole-line expansion with modern nucleon-nucleon potentials. Phys. Rev. 2017, C96, 044309. [Google Scholar] [CrossRef]
  64. Li, Z.H.; Schulze, H.-J. Correlation strength with modern NN potentials in the BHF approach. Phys. Rev. 2016, C94, 024322. [Google Scholar]
  65. Rijken, T.A.; Stoks, V.G.J.; Yamamoto, Y. Soft-core hyperon-nucleon potentials. Phys. Rev. C 1999, 59, 21. [Google Scholar] [CrossRef]
  66. Stoks, V.G.J.; Rijken, T.A. Soft-core baryon-baryon potentials for the complete baryon octet. Phys. Rev. C 1999, 59, 3009. [Google Scholar] [CrossRef]
  67. Machleidt, R.; Holinde, K.; Elster, C. The bonn meson-exchange model for the nucleon—nucleon interaction. Phys. Rep. 1987, 149, 1. [Google Scholar] [CrossRef]
  68. Machleidt, R. The Bonn meson-exchange model for the nucleon-nucleon interaction. Ad. Nucl. Phys. 1989, 19, 189. [Google Scholar] [CrossRef]
  69. Machleidt, R.; Sammarruca, F.; Song, Y. Nonlocal nature of the nuclear force and its impact on nuclear structure. Phys. Rev. C 1996, 53, 1483. [Google Scholar] [CrossRef] [PubMed]
  70. Pudliner, B.S.; Pandharipande, V.R.; Carlson, J.; Schiavilla, R.B. Quantum Monte Carlo Calculations of A≤6 Nuclei. Phys. Lett. B 1995, 74, 4396. [Google Scholar] [CrossRef] [PubMed]
  71. Fujita, J.; Miyazawa, H. Pion Theory of Three-Body Forces. Prog. Theor. Phys. 1957, 17, 360. [Google Scholar] [CrossRef]
  72. Stoks, V.G.J.; Klomp, R.A.M.; Terheggen, C.P.F.; de Swart, J.J. Construction of high-quality NN potential models. Phys. Rev. C 1994, 49, 2950. [Google Scholar] [CrossRef]
  73. Loiseau, B.A.; Nogami, Y.; Ross, C.K. Nucleon-nucleon correlation and two-pion-exchange three-body force in nuclear matter. Nucl. Phys. A 1971, 165, 601. [Google Scholar] [CrossRef]
  74. Grangé, P.; Martzolff, M.; Nogami, Y.; Sprung, D.W.L.; Ross, C.K. Three-body force in nuclear matter. Phys. Lett. B 1976, 60, 237. [Google Scholar] [CrossRef]
  75. Baldo, M.; Ferreira, L. Nuclear liquid-gas phase transition. Phys. Rev. C 1999, 59, 682. [Google Scholar] [CrossRef]
  76. Li, Z.H.; Lombardo, U.; Schulze, H.-J.; Zuo, W. Consistent nucleon-nucleon potentials and three-body forces. Phys. Rev. C 2008, 77, 034316. [Google Scholar] [CrossRef]
  77. Logoteta, D.; Vidaña, I.; Bombaci, I.; Kievsky, A. Comparative study of three-nucleon force models in nuclear matter. Phys. Rev. 2015, C91, 064001. [Google Scholar] [CrossRef]
  78. Shang, X.L.; Li, A.; Miao, Z.Q.; Burgio, G.F.; Schulze, H.-J. Nucleon effective mass in hot dense matter. Phys. Rev. 2020, C101, 065801. [Google Scholar] [CrossRef]
  79. Zuo, W.; Cao, L.G.; Li, B.A.; Lombardo, U.; Shen, C.W. Isospin splitting of the nucleon mean field. Phys. Rev. C 2005, 72, 014005. [Google Scholar] [CrossRef]
  80. Li, Z.H.; Lombardo, U.; Schulze, H.-J.; Zuo, W.; Chen, L.W.; Ma, H.R. Nuclear matter saturation point and symmetry energy with modern NN potentials. Phys. Rev. 2006, C74, 047304. [Google Scholar]
  81. Wlazłowski, G.; Holt, J.W.; Moroz, S.; Bulgac, A.; Roche, K.J. Auxiliary-Field Quantum Monte Carlo Simulations of Neutron Matter in Chiral Effective Field Theory. Phys. Rev. Lett. 2014, 113, 182503. [Google Scholar] [CrossRef]
  82. Tews, I.; Gandolfi, S.; Gezerlis, A.; Schwenk, A. Quantum Monte Carlo calculations of neutron matter with chiral three-body forces. Phys. Rev. 2016, C93, 024305. [Google Scholar] [CrossRef]
  83. Drischler, C.; Hebeler, K.; Schwenk, A. Neutron matter from chiral two- and three-nucleon calculations up to N 3 LO. Phys. Rev. 2016, C93, 054314. [Google Scholar]
  84. Bender, M.; Heenen, P.-H.; Reinhard, P.-G. Self-consistent mean-field models for nuclear structure. Rev. Mod. Phys. 2003, 75, 121. [Google Scholar] [CrossRef]
  85. Stone, J.R.; Reinhard, P.-G. The Skyrme interaction in finite nuclei and nuclear matter. Prog. Part. Nucl. Phys. 2007, 58, 587. [Google Scholar] [CrossRef]
  86. Erler, J.; Klüpfel, P.; Reinhard, P.-G. Self-consistent nuclear mean-field models: Example Skyrme-Hartree-Fock. J. Phys. 2011, G38, 033101. [Google Scholar] [CrossRef]
Figure 1. Spin−isospin ( S , T ) channel decomposition of the correlation energy per particle in MeV for SM and NM as a function of the density and for the various interactions considered. The bands quantify the internal accuracy of the BHF calculation associated with the different 2BFs and 3BFs used.
Figure 1. Spin−isospin ( S , T ) channel decomposition of the correlation energy per particle in MeV for SM and NM as a function of the density and for the various interactions considered. The bands quantify the internal accuracy of the BHF calculation associated with the different 2BFs and 3BFs used.
Universe 10 00226 g001
Figure 2. Decomposition in fixed L partial waves of the correlation energy per particle in SM as a function of the density for the various interactions.
Figure 2. Decomposition in fixed L partial waves of the correlation energy per particle in SM as a function of the density for the various interactions.
Universe 10 00226 g002
Figure 3. Same as Figure 2 in NM.
Figure 3. Same as Figure 2 in NM.
Universe 10 00226 g003
Figure 4. Contributions to the correlation energy per particle from various partial waves, as a function of the density for the different interactions.
Figure 4. Contributions to the correlation energy per particle from various partial waves, as a function of the density for the different interactions.
Universe 10 00226 g004
Figure 5. (left) Effective mass in SM and (right) isospin splitting of the effective mass for β = 0.5 as a function of density for the various forces.
Figure 5. (left) Effective mass in SM and (right) isospin splitting of the effective mass for β = 0.5 as a function of density for the various forces.
Universe 10 00226 g005
Figure 7. BHF correlation energy (left panels) and total energy (right panels) in SM (upper panels) and NM (lower panels) as a function of the density for the various N N potentials only (2BF) and including 3BF (2+3BFph, 2+3BFmic). Several χ EFT results [81,82,83] are also shown for comparison, see text for more details.
Figure 7. BHF correlation energy (left panels) and total energy (right panels) in SM (upper panels) and NM (lower panels) as a function of the density for the various N N potentials only (2BF) and including 3BF (2+3BFph, 2+3BFmic). Several χ EFT results [81,82,83] are also shown for comparison, see text for more details.
Universe 10 00226 g007
Figure 8. Skyrme coefficients obtained from the fit of ε cor ( S , T ) and the effective mass as function of 1 / γ . The models are calibrated to the following quantities; m sat / m = 0.7 and Δ 1 , sat = 0.1 , and to the BHF correlation energies in the ( S , T ) channels in SM and ( S ) channels in NM.
Figure 8. Skyrme coefficients obtained from the fit of ε cor ( S , T ) and the effective mass as function of 1 / γ . The models are calibrated to the following quantities; m sat / m = 0.7 and Δ 1 , sat = 0.1 , and to the BHF correlation energies in the ( S , T ) channels in SM and ( S ) channels in NM.
Universe 10 00226 g008
Figure 9. Spin−isospin ( S , T ) channel decomposition of the correlation energy per particle in SM and NM as a function of the density for the various interactions considered. The bands quantify the internal accuracy of the BHF calculation associated with the different 2BFs and 3BFs used. The curves are the result of the fit of the Skyrme EDF: the dashed lines represent fits of 2BF; dotted lines, 2+3BFph; and dashed-dotted lines, 2+3BFmic.
Figure 9. Spin−isospin ( S , T ) channel decomposition of the correlation energy per particle in SM and NM as a function of the density for the various interactions considered. The bands quantify the internal accuracy of the BHF calculation associated with the different 2BFs and 3BFs used. The curves are the result of the fit of the Skyrme EDF: the dashed lines represent fits of 2BF; dotted lines, 2+3BFph; and dashed-dotted lines, 2+3BFmic.
Universe 10 00226 g009
Table 1. Partial-wave contributions to the different ( S , T ) channels using the spectroscopic notation   2 S + 1 L J .
Table 1. Partial-wave contributions to the different ( S , T ) channels using the spectroscopic notation   2 S + 1 L J .
Name ( S , T ) LPartial WavesLegend in Figure 2 and Figure 3
Singlet Odd(0,0)1   1 P 1 P
SO 3   1 F 3 F
5   1 H 5 H
7   1 J 7 J
Singlet Even(0,1)0   1 S 0 S
SE 2   1 D 2 D
4   1 G 4 G
6   1 I 6 I
Triplet Even(1,0)0   3 S 1 S
TE 2   3 D 1 , 2 , 3 D
4   3 G 3 , 4 , 5 G
6   3 I 5 , 6 , 7 I
Triplet Odd(1,1)1   3 P 0 , 1 , 2 P
TO 3   3 F 2 , 3 , 4 F
5   3 H 4 , 5 , 6 H
7   3 J 6 , 7 , 8 J
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vidaña, I.; Margueron, J.; Schulze, H.-J. Nuclear Matter Equation of State in the Brueckner–Hartree–Fock Approach and Standard Skyrme Energy Density Functionals. Universe 2024, 10, 226. https://doi.org/10.3390/universe10050226

AMA Style

Vidaña I, Margueron J, Schulze H-J. Nuclear Matter Equation of State in the Brueckner–Hartree–Fock Approach and Standard Skyrme Energy Density Functionals. Universe. 2024; 10(5):226. https://doi.org/10.3390/universe10050226

Chicago/Turabian Style

Vidaña, Isaac, Jérôme Margueron, and Hans-Josef Schulze. 2024. "Nuclear Matter Equation of State in the Brueckner–Hartree–Fock Approach and Standard Skyrme Energy Density Functionals" Universe 10, no. 5: 226. https://doi.org/10.3390/universe10050226

APA Style

Vidaña, I., Margueron, J., & Schulze, H. -J. (2024). Nuclear Matter Equation of State in the Brueckner–Hartree–Fock Approach and Standard Skyrme Energy Density Functionals. Universe, 10(5), 226. https://doi.org/10.3390/universe10050226

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop