Eur. Phys. J. A (2026) 62: 77
https://doi.org/10.1140/epja/s10050-026-01825-4

Review

Variational theory and parquet diagrams for nuclear systems: a comprehensive study of neutron matter

¹ Department of Physics, University at Buffalo, SUNY, 14260, Buffalo, NY, USA
² Institut für Theoretische Physik, Johannes Kepler Universität, 4040, Linz, Austria

^a This email address is being protected from spambots. You need JavaScript enabled to view it.

Received: 9 July 2025
Accepted: 8 February 2026
Published online: 14 April 2026

Abstract

The main task of microscopic many-body theory is to provide an understanding and an explanation of properties of macroscopic systems from no other information than the properties of the underlying Hamiltonian, the particle statistics, and the macroscopic geometry of the system. Variational wave functions have been a very successful technique for examining the ground state and dynamic properties of dense quantum fluids. Specifically the “optimized (Fermi)-hypernetted chain” ((F)HNC-EL) hierarchy of approximations reproduces, for simple interactions like electrons and quantum fluids, the most basic features of a self-bound many-body system, namely binding, saturation, and spinodal decomposition. The relationships between Green’s functions based perturbation theory and variational wave functions have been clarified in much detail in the past. In the language of Feynman diagrams, a self-consistent summation of ring- and ladder-diagrams is the minimum requirement for a satisfactory microscopic description of these basic features. Realistic nucleon–nucleon interactions pose further problems: The most difficult aspect of the nuclear many-body problem is caused by the form of the microscopic nucleon–nucleon interaction which depends at least on the spin, isospin, the relative orientation and angular momentum of the interacting particles. Simple variational methods as used for electrons and quantum fluids are inadequate. Correlation functions that have the same structure as the interactions lead to so-called “commutator diagrams” that have, so far, been mostly neglected. We have shown in the past that these corrections can be very important if the two-body interactions are very different in different spin- or isospin channels. To deal with the problem of realistic nuclear interactions we have combined techniques of the Jastrow–Feenberg variational method and the local parquet-diagram theory. In the language of diagrammatic perturbation theory, “commutator diagrams” can be identified with non-parquet diagrams. We examine the physical processes described by these terms and include the relevant diagrams in a way that is suggested by the Jastrow–Feenberg approach. We show that the corrections from non-parquet contributions are, at short distances, larger than all other many-body effects. We examine here neutron matter as a prototype of systems with state-dependent interactions. Calculations are carried out for neutrons interacting via the so-called $v_8$ version of four popular interactions. We determine the structure and effective interactions and apply the method to the calculation of the energetics, structure and dynamic properties such as the single-particle self-energy and the dynamic response functions as well as BCS pairing in both singlet and triplet states. We find that many-body correlations lead to a strong reduction of the spin-orbit interaction, and, therefore, to a suppression of the $_{}^3{\text{P}}_2$ and $_{}^3{\text{P}}_2$-$_{}^3{\text{F}}_2$ gaps. We also find pairing in $_{}^3{\text{P}}_0$ states; the strength of the pairing gap depends sensitively on the potential model employed.