# Dirac-Brueckner-Hartree-Fock versus chiral effective field theory

###### Abstract

We compare nuclear and neutron matter predictions based on two different ab initio approaches to nuclear forces and the nuclear many-body problem. The first consists of a realistic meson-theoretic nucleon-nucleon potential together with the relativistic counterpart of the Brueckner-Hartree-Fock theory of nuclear matter. The second is based on chiral effective field theory, with density-dependent interactions derived from leading order chiral three-nucleon forces. We find the results to be very close and conclude that both approaches contain important features governing the physics of nuclear and neutron matter.

###### pacs:

21.65.+f, 21.30.Fe## I Introduction

Nuclear matter is a convenient laboratory to test nuclear forces and many body theories. In particular, the equation of state (EoS) of extremely neutron-rich matter has attracted considerable attention lately because of its broad applications, ranging from the structure of rare isotopes to the properties of neutron stars.

Constraints on the properties of nuclear matter, symmetric or isospin-asymmetric, can be obtained from a variety of experiments, such as measurements of nuclear binding energies (including isobaric analog state energies), parity-violating electron scattering, neutron skin thickness measurements, nucleus-nucleus collisions, and astrophysical observations. For a recent review of available constraints, particularly on the nuclear symmetry energy, the reader is referred to Ref. Tsang+ ().

Diverse theoretical frameworks have been employed to describe the properties of nuclear and neutron matter. They include: phenomenological approaches, both relativistic and non-relativistic; the Brueckner-Hartree-Fock (BHF) method, typically implemented with three-nucleon forces (3NF); variational approaches; the relativistic Dirac-Brueckner-Hartree-Fock (DBHF) approach; and chiral effective field theories. Predictions are model-dependent, particularly at the higher densities, where constraints are scarse and less stringent.

In our previous work with both symmetric nuclear matter (SNM) and isospin-asymmetric nuclear matter (IANM), we have relied on the DBHF approach, together with a relativistic meson-theoretic nucleon-nucleon (NN) potential which uses the pseudovector coupling for the pion Sam10 (). The latter choice was motivated by the considerations we outline next. Already when QCD (and its approximate symmetries) were unknown, it was observed that the contribution from the nucleon-antinucleon pair diagram, Fig. 1, becomes unreasonably large if the pseudoscalar (ps) coupling is used for the pion, leading to very large pion-nucleon scattering lengths Wei68 (); GB79 (). We recall that the Lagrangian density for pseudoscalar coupling of the nucleon field () with a pseudoscalar meson field () is

(1) |

On the other hand, the same contribution shown in Fig. 1, is heavily suppressed when the pseudovector (pv) coupling is used instead (a mechanism which became known as “pair suppression”). The reason for the suppression is that the covariant derivative in the pseudovector Lagrangian,

(2) |

generates a vertex that is proportional to momentum (leading to a weak coupling for low momenta) and, thus, explains the small value of the pion-nucleon scattering length at threshold. Non-linear realizations of chiral symmetry Wei68 () can further motivate a preference for the pseudovector coupling.

Back to the many-body problem, the main strength of the DBHF approach is its inherent ability to account for important three-body forces through its density dependence. In Fig. 2 we show a three-body force originating from virtual excitation of a nucleon-antinucleon pair, known as “Z-graph”. Notice that the observations from the previous paragraph ensure that this diagram, as well as the corresponding diagram at the two-body level, Fig. 1, is moderate in size when the pv coupling, Eq. (2), is used. (Hence, the importance of using pv potentials in relativistic nuclear structure calculations.) The main feature of the DBHF method turns out to be closely related to the 3NF depicted in Fig. 2, as we will argue next. In the DBHF approach, one describes the positive energy solutions of the Dirac equation in the medium as

(3) |

where the nucleon effective mass, , is defined as , with an attractive scalar potential. (This will be derived below.) It can be shown that both the description of a single-nucleon via Eq. (3) and the evaluation of the Z-diagram, Fig. 2, generate a repulsive effect on the energy/particle in symmetric nuclear matter which depends on the density approximately as

(4) |

and provides the saturating mechanism missing from conventional Brueckner calculations.

The approximate equivalence of the effective-mass description of Dirac states and the contribution from the Z-diagram has a simple intuitive explanation in the observation that Eq. (3), like any other solution of the Dirac equation, can be written as a superposition of positive and negative energy solutions using free nucleon masses. On the other hand, the “nucleon” in the middle of the Z-diagram, Fig. 2, is precisely a superposition of positive and negative energy states. In summary, the DBHF method effectively takes into account a particular class of 3NF, which are crucial for nuclear matter saturation.

The conventional BHF theory, together with meson-exchange 3NF, can also describe the saturation of nuclear matter in a satisfactory manner, although consistency between the parameters of the two- and the three-body forces can be problematic Catania ().

Ideally, one wishes to derive nuclear forces from the fundamental theory of strong interactions, QCD. Such task is however not feasible, due to the non-perturbative nature of the theory in the low-energy regime typical of nuclear physics. An alternative way is to respect the symmetries of the QCD Lagrangian while retaining the traditional degrees of freedom relevant to nuclear physics, namely nucleons and pions Wein79 (). This is the philosophy of chiral effective field theorie (EFT), which has become popular in recent years.

Effective field theories allow for a systematic expansion in powers of the momentum (or the pion mass) known as chiral perurbation theory (ChPT) such that, at any order, the irreducible two- and many-body diagrams to be included are precisely defined. An extensive review on the subject, including a comprehensive list or references, can be found in Ref. ME11 ().

Chiral EFT has validity only up to the chiral symmetry breaking scale of GeV. Thus, the low-momentum expansion has a limited range of applicability. This is another reason why relativistic meson theory has been considered more suitable for applications to dense systems, where high momenta are involved due to the high Fermi momenta. To remedy this problem, ways to extend EFT predictions to higher densities are employed, such as parametrizing the available predictions with accurate fitting functions which are then used to predict the properties of dense astrophysical systems.

It is the purpose of this paper to conduct a comparison between nuclear and neutron matter (NM) predictions obtained with two ab-initio approaches:

We would like to explore how the DBHF phenomenology compares with the predictions of chiral EFT at NLO plus leading 3NF, in both symmetric and neutron matter.

We will start with a comparison at the two-body force level. The meson-exchange side of this part of the comparison will be represented by the conventional BHF model. In the meson-exchange sector, realistic saturation properties will be generated through the DBHF model, whereas, on the chiral side, the contribution from leading 3NF will be included. We wish to compare the size and density dependence of the saturating effects in each scheme. We will compare predictions of the equation of state and the symmetry energy in both cases. Neutron star masses and radii will also be addressed.

This paper is organized as follows: In Section II, we review the main points of the DBHF calculation leading to the energy per particle in nuclear matter with arbitrary degree of isospin asymmetry. A brief outline of chiral two- and three-nucleon forces is contained in Section III. Our findings and conclusions are presented in Sections IV and V, respectively.

## Ii Relativistic meson-exchange potential and DBHF approach

As stated in the Introduction, the starting point of our many-body calculation is a realistic NN interaction which is then applied in the nuclear medium without any additional free parameters.

Relativistic meson theory is an appropriate framework to deal with the high momenta encountered in dense matter. In particular, the one-boson-exchange (OBE) model has proven very successful in describing NN data in free space up to high energy and has a good theoretical foundation. The OBE potential is defined as a sum of one-particle-exchange amplitudes of certain bosons with given mass and coupling. In general, six non-strange bosons with masses below 1 GeV are used. Thus,

(5) |

with and pseudoscalar, and scalar, and and vector particles. For more details, see Ref. Mac89 ().

Among the many available OBE potentials, some being part of the “high-precision generation” pot1 (); pot2 (), we seek a momentum-space potential developed within a relativistic scattering equation, such as the one obtained through the Thompson Thom () three-dimensional reduction of the Bethe-Salpeter equation BS ().

Having summarized in the Introduction the main DBHF philosophy, for completeness we now proceed to review the main aspects of our approach and the various approximations we perform through the application of the DBHF procedure. The equations we present are those suitable for isospin-asymmetric nuclear matter (IANM), since they naturally contain both the cases of SNM and NM.

We start from the Thompson Thom () relativistic three-dimensional reduction of the Bethe-Salpeter equation BS (). The Thompson equation is applied to nuclear matter in strict analogy to free-space scattering and reads, in the nuclear matter rest frame,

(6) |

where is the in-medium reaction matrix (=, , or ), and the asterix signifies that medium effects are applied to those quantities. Thus the NN potential, , is constructed in terms of effective Dirac states (in-medium spinors) as explained above. In Eq. (6), , , and are the initial, final, and intermediate relative momenta, and . The momenta of the two interacting particles in the nuclear matter rest frame have been expressed in terms of their relative momentum and the center-of-mass momentum, , through

(7) |

and

(8) |

The energy of the two-particle system is

(9) |

and is the starting energy. The single-particle energy includes kinetic energy and potential energy contributions (see Eq. (23) below). The Pauli operator, , prevents scattering to occupied , , or states. To eliminate the angular dependence from the kernel of Eq. (6), it is customary to replace the exact Pauli operator with its angle-average. Detailed expressions for the Pauli operator and the average center-of-mass momentum in the case of two different Fermi seas can be found in Ref. AS03 ().

With the definitions

(10) |

and

(11) |

one can rewrite Eq. (6) as

(12) |

which is formally identical to its non-relativistic counterpart.

The goal is to determine self-consistently the nuclear matter single-particle potential which, in IANM, will be different for neutrons and protons. To facilitate the description of the procedure, we will use a schematic notation for the neutron/proton potential. We write, for neutrons,

(13) |

and for protons

(14) |

where each of the four pieces on the right-hand-side of Eqs. (13-14) signifies an integral of the appropriate -matrix elements (, , or ) obtained from Eq. (12). Clearly, the two equations above are coupled through the component and so they must be solved simultaneously. Furthermore, the -matrix equation and Eqs. (13-14) are coupled through the single-particle energy (which includes the single-particle potential, itself defined in terms of the -matrix). So we have a coupled system to be solved self-consistently.

Before proceeding with the self-consistency, one needs an ansatz for the single-particle potential. The latter is suggested by the most general structure of the nucleon self-energy operator consistent with all symmetry requirements. That is:

(15) |

where and are an attractive scalar field and a repulsive vector field, respectively, with the timelike component of the vector field. These fields are in general density and momentum dependent. We take

(16) |

which amounts to assuming that the spacelike component of the vector field is much smaller than both and . Furthermore, neglecting the momentum dependence of the scalar and vector fields and inserting Eq. (16) in the Dirac equation for neutrons/protons propagating in nuclear matter,

(17) |

naturally leads to rewriting the Dirac equation in the form

(18) |

with positive energy solutions as in Eq. (3), , and

(19) |

The subscript “” signifies that these parameters are different for protons and neutrons.

As in the symmetric matter case BM84 (), evaluating the expectation value of Eq. (16) leads to a parametrization of the single particle potential for protons and neutrons (Eqs.(13-14)) in terms of the constants and which is given by

(20) |

Also,

(21) |

which, along with Eq. (20), allows the self-consistent determination of the single-particle potential as explained below.

The kinetic contribution to the single-particle energy is

(22) |

and the single-particle energy is

(23) |

The constants and

(24) |

are convenient to work with as they facilitate the connection with the usual non-relativistic framework HT70 ().

Starting from some initial values of and , the -matrix equation is solved and a first approximation for is obtained by integrating the -matrix over the appropriate Fermi sea, see Eq. (21). This solution is again parametrized in terms of a new set of constants, determined by fitting the parametrized , Eq. (20), to its values calculated at two momenta, a procedure known as the “reference spectrum approximation”. The iterative procedure is repeated until satisfactory convergence is reached.

Finally, the energy per neutron or proton in nuclear matter is calculated from the average values of the kinetic and potential energies as

(25) |

The EoS, or energy per nucleon as a function of density, is then written as

(26) |

or

(27) |

Clearly, symmetric nuclear matter is obtained as a by-product of the calculation described above by setting =0, whereas =1 corresponds to pure neutron matter.

## Iii Chiral two- and three-nucleon forces

Ideally, one wishes to base a derivation of the nuclear force on QCD. However, the well-known problem with QCD is that it is non-perturbative in the low-energy regime characteristic for nuclear physics. For many years this fact was perceived as a great obstacle to a derivation of nuclear forces from QCD–impossible to overcome except with lattice QCD. The effective field theory concept has shown the way out of this dilemma. One has to realize that the scenario of low-energy QCD is characterized by pions and nucleons interacting via a force governed by spontaneously broken approximate chiral symmetry. This chiral EFT allows for a systematic low-momentum expansion known as chiral perturbation theory (ChPT). Contributions are analyzed in terms of powers of small external momenta over the large scale, , where is generic for an external momentum (nucleon three-momentum or pion four-momentum) or pion mass, and 1 GeV is the chiral symmetry breaking scale (‘hard scale’). (See Ref. ME11 () and references therein.)

The past fifteen years have seen great progress in applying ChPT to nuclear forces. As a result, NN potentials of high precision have been constructed, which are based on ChPT carried to next-to-next-to-next-to-leading order (NLO). We will apply here the chiral NN potential of Ref. EM03 () which uses a cutoff, , equal to 500 MeV.

A great advantage of the EFT approach to nuclear forces is that it creates two- and many-body forces on an equal footing. Three-nucleon forces make their appearance at the third order in the chiral power counting. These leading-order contributions are: the long-range two-pion exchange graph; the medium-range one-pion exchange diagram; and the short-range contact term.

In Ref. Holt (), density-dependent corrections to the in-medium NN interaction have been derived from the leading-order chiral 3NF. These are effective two-nucleon interactions that reflect the underlying three-nucleon forces and are therefore computationally very convenient, whereas realistic models of three-nucleon forces would be prohibitive.

A total of six one-loop diagrams contribute at this order. Three are generated by the two-pion exchange graph of the chiral three-nucleon interaction and depend on the low-energy constants , which are fixed in the NN system EM03 (). We use GeV, GeV, and GeV. Two are generated by the one-pion exchange diagram and depend on the low-energy constant . Finally, the short-range component depends on the constant . The constants and can be fixed by fitting properties of few-nucleon systems, such as the triton and He binding energies Nav07 (). We use and .

In pure neutron matter, the contributions proportional to the low-energy constants , and vanish Holt (). Analytical expressions for these corrections are provided in Ref. Holt () in terms of the well-known non-relativistic two-body nuclear force operators. These can be conveniently incorporated in the usual NN partial wave formalism and the conventional BHF theory.

## Iv Results

In Fig. 3, we display the EoS of SNM obtained with the two approaches outlined above. The (blue) solid lines denote the results from the meson-exchange Bonn potential and the conventional BHF approximation (“BnB BHF”) or the relativistic DBHF (“BnB DBHF”), while the chiral calculations are shown by the (red) dashed lines, with and without 3NF, as denoted.

As to be expected, both the conventional BHF calculation and the one using only the two-body chiral potential NLO display excessive attraction and are unable to produce saturation up to very high density. The “Dirac effect” brings in a powerful saturation mechanism as seen by comparing the two (blue) solid curves. A very similar effect, both qualitatively and quantitatively, is generated by the effective chiral 3NF.

In spite of the differences between the two approaches, the predictions are very close, not only in the value of the energy at saturation, but also the density dependence of the saturation mechanism. We notice, though, that the DBHF energies tend to grow faster at the higher densities.

The same comparison is done in Fig. 4 for neutron matter, yielding similar conclusions. Here, too, the DBHF energies show a faster growth as density increases.

In Fig. 5, we display the symmetry energy as obtained with the well-known parabolic approximation, . (Notice that at this point we only consider the versions of the BnB- and the N3LO-based calculations which have realistic saturation properties, see Fig. 3.)

The slope of the symmetry energy, usually defined in terms of the parameter, (), is a measure of the pressure gradient between neutron and symmetric matter. A variety of experiments are aimed at constraining this important quantity Tsang+ (), which correlates nearly linearly with the neutron skin of neutron-rich nuclei. At the present time there seems to be a consensus from various experiments that the acceptable range of values of the symmetry energy and its slope at saturation are centered around and equal to 32.5 MeV and 70 MeV, respectively Tsang+ (). The values of the parameter obtained with “BnB+DBHF” and “N3LO+3NF” are 69 and 61 MeV, respectively. Notice that the typical uncertainty on the parameter is as large as 20-25 MeV Tsang+ ().

How the differences/similarities noted above impact neutron star bulk properties is examined in Fig. 6. (Since our focal point is a comparison between two theoretical approaches, rather than a detailed calculation of -stable matter composition, we consider only neutron stars made of pure neutron matter.) We also note that, in order to extend the chiral predictions to the high densities probed by compact stars, we fit a three-parameter function, =, to the “N3LO+3NF” predictions of Fig. 4 and use this ansatz to obtain the high-density EoS needed for neutron star calculations.

The comparison in Fig. 6 shows a larger star maximum mass for “BnB DBHF”, most likely the result of larger repulsion of that model at high density. Overall, the predicted star bulk properties are only moderately different in the two approaches.

## V Conclusions

We have considered two very different methods to approach the study of nucleonic matter: one based on a meson-theoretic potential and the Dirac-Brueckner-Hartree-Fock approximation; the other based on a high-precision chiral NN potential and chiral effective three-nucleon forces at NNLO. The predictions we have considered include: the EoS of nuclear and neutron matter, the symmetry energy and its slope, and the mass-radius relation in a neutron star.

From our results, we conclude that the DBHF method is an excellent phenomenology capable of incorporating important many-body effects that are crucial to nuclear saturation.

In both approaches, the effective 3NF is generated by one nucleon interacting with the Fermi sea. That is, in both cases we have effective two-nucleon interactions that reflect the underlying three-nucleon forces. No matter if this interaction proceeds via relativistic meson exchange or via chiral EFT forces, the results are very similar. This is reassuring and confirms that the two ways of describing nuclear forces are complementary.

## Acknowledgments

Support from the U.S. Department of Energy under Grant No. DE-FG02-03ER41270 is acknowledged, and from the Italian Ministero dell’Istruzione dell’Universitá della Ricerca (MIUR) under PRIN 2009. We are grateful to F. Weber for providing his TOV code for the calculation of neutron star properties.

## References

- (1) M.B. Tsang et al., Phys. Rev. C 86, 015803 (2012).
- (2) F. Sammarruca, Int. J. Mod. Phys. E, 19, 1259 (2010).
- (3) S. Weinberg, Phys. Rev. 166, 1568 (1968).
- (4) G.E. Brown, in Mesons in Nuclei, eds. M. Rho and D.H. Wilkinson, Vol. I (North-Holland, Amsterdam, 1979), 330.
- (5) Z.H. Li, U. Lombardo, H.-J. Schulze, and W. Zuo, Phys. Rev. C 77, 034316 (2008).
- (6) S. Weinberg, Physica 96A, 327 (1979).
- (7) R. Machleidt and D.R. Entem, Phys. Rep. 503, 1 (2011).
- (8) R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989).
- (9) D.R. Entem and R. Machleidt, Phys. Rev. C 68, 041001 (2003).
- (10) R. Machleidt, Phys. Rev. C 63, 024001 (2001).
- (11) V.G.J. Stoks et al., Phys. Rev. C 49, 2950 (1994).
- (12) R.H. Thompson, Phys. Rev. D 1, 110 (1970).
- (13) E.E. Salpeter and H.A. Bethe, Phys. Rev. 84, 1232 (1951).
- (14) D. Alonso and F. Sammarruca, Phys. Rev. C 67, 054301 (2003).
- (15) R. Brockmann and R. Machleidt, Phys. Lett. B 149, 283 (1984); Phys. Rev. C 42, 1965 (1990).
- (16) M.I. Haftel and F. Tabakin, Nucl. Phys. A158, 1 (1970).
- (17) J.W. Holt, N. Keiser, and W. Weise, Phys. Rev. C 81, 024002 (2010).
- (18) P. Navratil, V.G. Gueorguiev, J.P. Vary, W.E. Ormand, and A. Nogga, Phys. Rev. Lett. 99, 042501 (2007).