Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2507.03825 · The European Physical Journal A

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

E. Krotscheck, Jia-Wei Wang

Published Apr 14, 2026• 0 citations
arXivHTMLPDF

Abstract

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$$ 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\textrm{P}_2$$ 3 P 2 and $$^3\textrm{P}_2$$ 3 P 2 - $$^3\textrm{F}_2$$ 3 F 2 gaps. We also find pairing in $$^3\textrm{P}_0$$ 3 P 0 states; the strength of the pairing gap depends sensitively on the potential model employed.

DOI
10.1140/epja/s10050-026-01825-4

Related papers

Variational and parquet-diagram calculations for neutron matter. II. Twisted chain diagrams

E. Krotscheck, Jiawei Wang

We develop a manifestly microscopic method to deal with strongly interacting nuclear systems that have different interactions in spin-singlet and spin-triplet states. In a first step we analyze variational wave functions that have been suggested to describe such systems and demonstrate that the so-called commutator contributions can have important effects whenever the interactions in the spin-singlet and the spin-triplet states are very different. We then identify these contributions as terms that correspond, in the language of perturbation theory, to non-parquet diagrams. We include these diagrams in a way that is suggested by the Jastrow-Feenberg approach and show that the corrections from non-parquet contributions are, at short distances, larger than all other many-body effects.

arXiv:2009.10849•Physical Review C•Dec 3, 2020•5 citations

Variational and parquet-diagram calculations for neutron matter: Structure and energetics

E. Krotscheck, Jiawei Wang

We develop the variational/parquet diagram approach to the structure of nuclear systems with strongly state-dependent interactions. For that purpose, we combine ideas of the general Jastrow-Feenberg variational method and the local parquet-diagram theory for bosons with state-dependent interactions [R. A. Smith and A. D. Jackson, Nucl. Phys. A 476, 448 (1988)]. The most tedious aspect of variational approaches, namely the symmetrization of an operator dependent variational wave function, is thereby avoided. We carry out calculations for neutron matter interacting via the Reid and Argonne ${v}_{6}$ models of the nucleon-nucleon interaction. While the equation of state is a rather robust quantity that comes out reasonably well even in very simplistic approaches, we show that effective interactions, which are the essential input for calculating dynamic properties, depend sensitively on the quality of the treatment of the many-body problem.

arXiv:1912.05685•Physical Review C•Jun 25, 2020•7 citations

A variational theory of nuclear matter

Robert B. Wiringa

The hypotheticai system of nucleons known as nuclear matter provides a very interesting field of study with relevance to both nuclear physics and astrophysics. The task of finding nucleon-nucleon potentials that explain the empirical properties of nuclear matter, as well as the two-body data, has been hampered by the lack of an adequate many-body theory for treating the complicated forces involved. Realistic nucleon-nucleon potentials have strong spin, isospin, tensor and spin-orbit components which induce correspondingly strong correlations in nuclear matter. In this thesis we present a variational theory designed to cope with these complicated correlations. To place the theory in perspective we give an overview of several different many-body theories and compare their predictions in the somewhat simpler systems of the 4He and 3He liquids, and for a simple central potential model of nuclear matter. The variational theory is then developed, starting with the choice of wave function. A generalized Jastrow product of two-body correlation functions containing central, spin, isospin, tensor, and spin-orbit operators is used. The correlation function is generated by a series of two-body Schrodinger equations with boundary conditions which require the wave function to heal at a distance d. The wave function is parameterized by d and by the magnitudes of the non-central correlations. Expectation values with variational wave functions of the Jastrow type may be studied conveniently with generalized Mayer diagrams. Diagram rules are given and a general diagrammatic cluster expansion is derived. It is an expansion in powers of commutators of the non-central operators, and in their absence reduces to the usual Mayer cluster expansion. Evaluation of diagrams in the cluster expansion is complicated by the presence of the non-commuting operators. Some simple rules and useful methods for calculating the contribution due to the central, spin, isospin, and tensor operators are given. Spin-orbit operators are neglected beyond this point, but the approach used is extendible to them. Large classes of diagrams contributing to an expectation value can be summed by means of integral equations. The development of single chain, hypernetted chain, and Fermi hypernetted chain (FHNC) equations for central correlations is reviewed. The chain summation methods are extended to sum single operator chains (SOC) for the non-central operators. Methods are also given for treating the leading commutator corrections in the expansi0n by simple vertex factors. The energy expectation value for potentials with six operators are evaluated using the FHNC/SOC functions, with an exact treatment of the commutators involved. The energy is found to have a minimum with respect to variations in all parameters. Results of calculations with model potentials based on the Reid, Bethe-Johnson, Hamada-Johnston, and Gammel-Thaler potentials are reported. A crude estimate of the effect of the neglected spin-orbit potentials on the nuclear matter binding energy and saturation density indicates that it could be significant. A correct treatment of the spin-orbit potentials is the next major step that should be taken.

•Jan 1, 1978•0 citations

α -cluster matter reexamined

J. W. Clark, E. Krotscheck

We examine in detail two alternative descriptions of a system of $\ensuremath{\alpha}$ particles interacting via local interactions of different character, highlighting the fact that a faithful microscopic description of such systems demands a consistent treatment of both short- and long-range correlations. In preparation, we examine four different versions of modern microscopic many-body theory and conclude by emphasizing that these approaches, although a priori very different, actually lead to the same equations for their efficient application. The only quantity that depends on the formulation of many-body theory chosen is an irreducible interaction correction. In the language of Green's functions and Feynman diagrams, it is the set of both particle-particle and particle-hole irreducible diagrams, and in variational Jastrow-Feenberg theory it is determined by multipartite correlations and elementary diagrams. We apply these theoretical methods to the calculation of the energetics, structure, thermodynamics, and dynamics of $\ensuremath{\alpha}$ matter, as well as its condensate fraction. In dimensionless units, $\ensuremath{\alpha}$ matter appears to be remarkably similar to the much-studied $^{4}\mathrm{He}$ quantum fluid, its low-temperature properties now basically solved in the Jastrow-Feenberg framework. Accordingly, one can have confidence in the results of application of the same procedure to $\ensuremath{\alpha}$ matter. Even so, closer examination reveals significant differences between the physics of the two systems. Within an infinite nuclear medium, $\ensuremath{\alpha}$ matter is subject to a spinoidal instability. Extended mixtures of nucleons and $\ensuremath{\alpha}$ particles are yet to be given rigorous consideration in a corresponding theoretical framework.

Physical Review C•Mar 21, 2024•4 citations

Variational Theory of Quantum Fluids.

Kevin E. Schmidt

"A theory is developed to allow the systematic calculation of the properties of quantum fluids given the microscopic two-body potential. Physical arguments are used to construct a correlation operator that approximately transforms the states of a non-interacting system into the corresponding states of the interacting-system. The theory can be regarded as a microscopic realization of Landau's assumption of a one-to-one correspondence between the non-interacting and interacting states. It is variational since the correlation operator is determined by minimizing E(p) the expectation value of the ground-state energy. We first give a survey of various other many-body theories including the simpler Jastrow variational method. Differential equations are derived to calculate the correlation operator and the related Jastrow correlation"" function. A comparison of results in several Bose systems with various methods is given. The Fermi hypernetted chain integral equations are generalized to allow the calculation of E(p) in Fermi systems, using the correlation operator wave function. The method is applied to liquid 3He and the v2-model of nuclear matter. In 3He the correlation operator gives results in much closer agreement with experiment than that with the Jastrow wave function. A new set of hypernetted chain equations are then derived to allow the calculation of the excited st;ate energies in Bose systems. The theory is then applied to the calculation of the single quasi-particle excitation spectrum of liquid 4He. The spectrum is found to be a stringent test on the form of the correlation operator, deficiencies in the correlation operator giving corresponding errors in the spectrum. Good results are obtained for the density dependence of the roton energy. A variational principle for the free energy is used to generalize the correlation operator technique to finite temperatures. Results are given for the finite temperature properties of the v8-model of nuclear matter."

•Aug 1, 1979•0 citations

The nature of short-ranged correlations in nuclear systems

E. Krotscheck, Jiawei Wang

We apply ideas of the parquet-diagram and optimized Fermi-hypernetted chain methods to determine the short-range structure of the pair wave function in neutron matter and compare these with Bethe-Goldstone results and those of low-order variational calculations. It is shown that the induced interaction, describing the exchange of density, spin, and tensor fluctuations, has a profound influence on the short-ranged structure of the pair wave function and, hence, on effective interactions in neutron matter.

arXiv:2006.11849•Physics Letters B•Sep 9, 2020•3 citations
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 v8v_8v8​ 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 3P2^3\textrm{P}_23P2​ 3 P 2 and 3P2^3\textrm{P}_23P2​ 3 P 2 - 3F2^3\textrm{F}_23F2​ 3 F 2 gaps. We also find pairing in 3P0^3\textrm{P}_03P0​ 3 P 0 states; the strength of the pairing gap depends sensitively on the potential model employed.
We develop a manifestly microscopic method to deal with strongly interacting nuclear systems that have different interactions in spin-singlet and spin-triplet states. In a first step we analyze variational wave functions that have been suggested to describe such systems and demonstrate that the so-called commutator contributions can have important effects whenever the interactions in the spin-singlet and the spin-triplet states are very different. We then identify these contributions as terms that correspond, in the language of perturbation theory, to non-parquet diagrams. We include these diagrams in a way that is suggested by the Jastrow-Feenberg approach and show that the corrections from non-parquet contributions are, at short distances, larger than all other many-body effects.
We develop the variational/parquet diagram approach to the structure of nuclear systems with strongly state-dependent interactions. For that purpose, we combine ideas of the general Jastrow-Feenberg variational method and the local parquet-diagram theory for bosons with state-dependent interactions [R. A. Smith and A. D. Jackson, Nucl. Phys. A 476, 448 (1988)]. The most tedious aspect of variational approaches, namely the symmetrization of an operator dependent variational wave function, is thereby avoided. We carry out calculations for neutron matter interacting via the Reid and Argonne v6{v}_{6}v6​ models of the nucleon-nucleon interaction. While the equation of state is a rather robust quantity that comes out reasonably well even in very simplistic approaches, we show that effective interactions, which are the essential input for calculating dynamic properties, depend sensitively on the quality of the treatment of the many-body problem.
The hypotheticai system of nucleons known as nuclear matter provides a very interesting field of study with relevance to both nuclear physics and astrophysics. The task of finding nucleon-nucleon potentials that explain the empirical properties of nuclear matter, as well as the two-body data, has been hampered by the lack of an adequate many-body theory for treating the complicated forces involved. Realistic nucleon-nucleon potentials have strong spin, isospin, tensor and spin-orbit components which induce correspondingly strong correlations in nuclear matter. In this thesis we present a variational theory designed to cope with these complicated correlations. To place the theory in perspective we give an overview of several different many-body theories and compare their predictions in the somewhat simpler systems of the 4He and 3He liquids, and for a simple central potential model of nuclear matter. The variational theory is then developed, starting with the choice of wave function. A generalized Jastrow product of two-body correlation functions containing central, spin, isospin, tensor, and spin-orbit operators is used. The correlation function is generated by a series of two-body Schrodinger equations with boundary conditions which require the wave function to heal at a distance d. The wave function is parameterized by d and by the magnitudes of the non-central correlations. Expectation values with variational wave functions of the Jastrow type may be studied conveniently with generalized Mayer diagrams. Diagram rules are given and a general diagrammatic cluster expansion is derived. It is an expansion in powers of commutators of the non-central operators, and in their absence reduces to the usual Mayer cluster expansion. Evaluation of diagrams in the cluster expansion is complicated by the presence of the non-commuting operators. Some simple rules and useful methods for calculating the contribution due to the central, spin, isospin, and tensor operators are given. Spin-orbit operators are neglected beyond this point, but the approach used is extendible to them. Large classes of diagrams contributing to an expectation value can be summed by means of integral equations. The development of single chain, hypernetted chain, and Fermi hypernetted chain (FHNC) equations for central correlations is reviewed. The chain summation methods are extended to sum single operator chains (SOC) for the non-central operators. Methods are also given for treating the leading commutator corrections in the expansi0n by simple vertex factors. The energy expectation value for potentials with six operators are evaluated using the FHNC/SOC functions, with an exact treatment of the commutators involved. The energy is found to have a minimum with respect to variations in all parameters. Results of calculations with model potentials based on the Reid, Bethe-Johnson, Hamada-Johnston, and Gammel-Thaler potentials are reported. A crude estimate of the effect of the neglected spin-orbit potentials on the nuclear matter binding energy and saturation density indicates that it could be significant. A correct treatment of the spin-orbit potentials is the next major step that should be taken.
We examine in detail two alternative descriptions of a system of $\ensuremath{\alpha}$ particles interacting via local interactions of different character, highlighting the fact that a faithful microscopic description of such systems demands a consistent treatment of both short- and long-range correlations. In preparation, we examine four different versions of modern microscopic many-body theory and conclude by emphasizing that these approaches, although a priori very different, actually lead to the same equations for their efficient application. The only quantity that depends on the formulation of many-body theory chosen is an irreducible interaction correction. In the language of Green's functions and Feynman diagrams, it is the set of both particle-particle and particle-hole irreducible diagrams, and in variational Jastrow-Feenberg theory it is determined by multipartite correlations and elementary diagrams. We apply these theoretical methods to the calculation of the energetics, structure, thermodynamics, and dynamics of $\ensuremath{\alpha}$ matter, as well as its condensate fraction. In dimensionless units, $\ensuremath{\alpha}$ matter appears to be remarkably similar to the much-studied 4He^{4}\mathrm{He}4He quantum fluid, its low-temperature properties now basically solved in the Jastrow-Feenberg framework. Accordingly, one can have confidence in the results of application of the same procedure to $\ensuremath{\alpha}$ matter. Even so, closer examination reveals significant differences between the physics of the two systems. Within an infinite nuclear medium, $\ensuremath{\alpha}$ matter is subject to a spinoidal instability. Extended mixtures of nucleons and $\ensuremath{\alpha}$ particles are yet to be given rigorous consideration in a corresponding theoretical framework.
"A theory is developed to allow the systematic calculation of the properties of quantum fluids given the microscopic two-body potential. Physical arguments are used to construct a correlation operator that approximately transforms the states of a non-interacting system into the corresponding states of the interacting-system. The theory can be regarded as a microscopic realization of Landau's assumption of a one-to-one correspondence between the non-interacting and interacting states. It is variational since the correlation operator is determined by minimizing E(p) the expectation value of the ground-state energy. We first give a survey of various other many-body theories including the simpler Jastrow variational method. Differential equations are derived to calculate the correlation operator and the related Jastrow correlation"" function. A comparison of results in several Bose systems with various methods is given. The Fermi hypernetted chain integral equations are generalized to allow the calculation of E(p) in Fermi systems, using the correlation operator wave function. The method is applied to liquid 3He and the v2-model of nuclear matter. In 3He the correlation operator gives results in much closer agreement with experiment than that with the Jastrow wave function. A new set of hypernetted chain equations are then derived to allow the calculation of the excited st;ate energies in Bose systems. The theory is then applied to the calculation of the single quasi-particle excitation spectrum of liquid 4He. The spectrum is found to be a stringent test on the form of the correlation operator, deficiencies in the correlation operator giving corresponding errors in the spectrum. Good results are obtained for the density dependence of the roton energy. A variational principle for the free energy is used to generalize the correlation operator technique to finite temperatures. Results are given for the finite temperature properties of the v8-model of nuclear matter."
We apply ideas of the parquet-diagram and optimized Fermi-hypernetted chain methods to determine the short-range structure of the pair wave function in neutron matter and compare these with Bethe-Goldstone results and those of low-order variational calculations. It is shown that the induced interaction, describing the exchange of density, spin, and tensor fluctuations, has a profound influence on the short-ranged structure of the pair wave function and, hence, on effective interactions in neutron matter.