Elementary diagrams in nuclear and neutron matter

Abstract

Variational calculations of nuclear and neutron matter are currently performed using a diagrammatic cluster expansion with the aid of nonlinear integral equations for evaluating expectation values. These are the Fermi hypernetted chain (FHNC) and single-operator chain (SOC) equations, which are a way of doing partial diagram summations to infinite order. A more complete summation can be made by adding elementary diagrams to the procedure. The simplest elementary diagrams appear at the four-body cluster level; there is one such E{sub 4} diagram in Bose systems, but 35 diagrams in Fermi systems, which gives a level of approximation called FHNC/4. We developed a novel technique for evaluating these diagrams, by computing and storing 6 three-point functions, S{sub xyz}(r{sub 12}, r{sub 13}, r{sub 23}), where xyz (= ccd, cce, ddd, dde, dee, or eee) denotes the exchange character at the vertices 1, 2, and 3. All 35 Fermi E{sub 4} diagrams can be constructed from these 6 functions and other two-point functions that are already calculated. The elementary diagrams are known to be important in some systems like liquid {sup 3}He. We expect them to be small in nuclear matter at normal density, but they might become significant at higher densities appropriate for neutron star calculations. This year we programmed the FHNC/4 contributions to the energy and tested them in a number of simple model cases, including liquid {sup 3}He and Bethe`s homework problem. We get reasonable, but not exact agreement with earlier published work. In nuclear and neutron matter with the Argonne v{sub 14} interaction these contributions are indeed small corrections at normal density and grow to only 5-10 MeV/nucleon at 5 times normal density.