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

Physical Review

Compact-Cluster Expansion for the Nuclear Many-Body Problem

B. H. Brandow

Published Dec 16, 1966• 184 citations
Publisher

Abstract

Bethe's treatment of three-body clusters is extended to the entire linked-cluster perturbation series. The result is a completely rearranged expansion in terms of "compact clusters." The spatial correlations within these terms are all of quite short range. As a direct consequence, each of these new compact-cluster terms is proportional to ${\ensuremath{\kappa}}^{{h}^{\ensuremath{'}}\ensuremath{-}1}$, where ${h}^{\ensuremath{'}}$ is the number of momenta inside the Fermi sea which can be summed independently after allowing for momentum conservation. The "small parameter" of the expansion is $\ensuremath{\kappa}=\ensuremath{\rho}\ensuremath{\int}{|\ensuremath{\zeta}|}^{2}d\ensuremath{\tau}$, where $\ensuremath{\rho}$ is the ordinary nuclear density and $\ensuremath{\zeta}=\ensuremath{\varphi}\ensuremath{-}{\ensuremath{\psi}}_{\mathrm{BG}}$ is the "wound" in the correlated two-body (bethe-Goldstone) wave function. In nuclear matter this $\ensuremath{\kappa}$ is of order 10%. This result requires that the single-particle potentials be defined in terms of a certain subset of all the self-energy insertions. The allowed insertions are those which can be evaluated entirely on the energy shell, by means of a "generalized time-ordering" factorization. The totality of these insertions is said to constitute an "on-energy-shell mass operator" ${M}^{\mathrm{on}}$. The resulting occupied-state potentials are quantitatively very similar to those of most previous versions of nuclear-matter theory. However, the present intermediate-state potentials are small enough (being only of order +1 MeV) to be safely ignored in practical calculations. This simplification is offset by the need for a separate calculation of the three-body clusters. These single-particle potentials are "physically meaningful" in the sense that they lead to an optimum treatment of the short-range correlations. The bulk properties of nuclear systems are determined almost entirely by these short-range correlations. It is therefore proposed that similar potentials be used in a many-body theory of finite nuclei. True occupation numbers occur in a simple and natural way throughout the entire expansion. The formal and practical consequences of this feature are carefully examined. It is argued that this is significant for a theory of finite nuclei. Comparisons are made with several different "renormalized" formulations of quantum statistical mechanics, and also with the Green's-function theory of nuclear matter. The present formulation tends to artificially suppress the long-range correlations (e.g., pairing correlations) between particles near the Fermi surface. This defect can be eliminated by using a degenerate analog of the Goldstone expansion.

DOI
10.1103/physrev.152.863

Related papers

Linked-Cluster Expansions for the Nuclear Many-Body Problem

B. H. Brandow

The Goldstone expansion is rederived by elementary time-independent methods, starting from Brillouin-Wigner (BW) perturbation theory. Interaction energy terms $\ensuremath{\Delta}E$ are expanded out of the BW energy denominators, and the series is then rearranged to obtain the linked-cluster result. Similar algebraic methods lead to the linked expansions for the total wave function (Hugenholtz) and the expectation value of a general operator (Thouless). Starting again with a degenerate version of BW perturbation theory, these methods are used to obtain the Bloch-Horowitz energy expansion, as well as the corresponding wave function, expectation-value, and transition-amplitude expansions. A "reduced" form of the Block-Horowitz expansion is described, and also a "completely linked" version. The latter is suggested as a tool for investigating superfiuid phenomena in nuclear matter, and for establishing contact with the Landau theory of Fermi liquids. The physical interpretation of these expansions is carefully studied, especially with regard to nuclear applications, to determine how they handle such "physical" features as antisymmetry, self-energy effects, wave-function renormalizations, and the distinction between "true" and "model" single-particle occupation probabilities. The problem of a correct theoretical definition for the shell-model potential is carefully examined, and a specific theory is presented. These expansions are seen to form a convenient and very powerful set of tools for studying the structure of actual nuclei.

Reviews of Modern Physics•Oct 1, 1967•1,275 citations

Many-Body Problem for Strongly Interacting Particles. II. Linked Cluster Expansion

K. A. Brueckner

An approximation method developed previously to deal with many particles in strong interaction is examined in further detail. It is shown that the series giving the interaction energy is a development in a sequence of linked or irreducible cluster terms each of which gives a contribution to the energy proportional to the total number of particles. Consequently the convergence of the expansion is independent of the total number of particles. The origin of this simple feature is illustrated by showing that a similar situation exists in the expansion of standard perturbation theory. The numerical convergence of the expansion is quantitatively discussed for the nuclear problem where it is shown that the correction arising from the first cluster term involving three particles is less than the leading term by a factor of about ${10}^{\ensuremath{-}4}$. The smallness of the correction is largely a result of the action of the exclusion principle.

Physical Review•Oct 1, 1955•770 citations

Many-Body Perturbation Theory

William Tobocman

The usual perturbation treatment of the many-body problem is modified by eliminating the two-body interaction and introducing in its place a $t$ matrix. The Watson, Brueckner, and Bethe (W.B.B.) treatment of the many-body problem is shown to be a special case of the modified perturbation treatment. The modified perturbation treatment involves no assumption except that of supposing the convergence of the perturbation expansion. Thus to each order in the expansion the results must be meaningful and no unphysical "unlinked cluster" terms can appear. Within the framework of the modified perturbation treatment it is easy to interpret the concepts of the Pauli principle for intermediate states and the self-consistent choice of the comparison potential which play such an important role in the W.B.B. treatment. This work is similar to that of J. Goldstone in that we relate the W.B.B. treatment to a perturbation treatment of the many-body problem. We differ from Goldstone in that we use the usual time-independent perturbation theory instead of the time-dependent one. The time-independent approach permits a direct comparison with the W.B.B. treatment and facilitates discussion of the Pauli principle for intermediate states and of the self-consistent potential.

Physical Review•Jul 1, 1957•28 citations

Method of Correlated Basis Functions

J. W. Clark, Paul A. Westhaus

The method of correlated basis functions provides a means of systematically exploiting our intuitive understanding of real physical systems. Here we seek especially to develop a formalism for describing systems of strongly interacting fermions, for example, nuclei and liquid ${\mathrm{He}}^{3}$, using a complete set of correlated functions ${\ensuremath{\Psi}}_{\mathrm{m}}=\frac{F{\ensuremath{\Phi}}_{\mathrm{m}}}{{({\ensuremath{\Phi}}_{\mathrm{m}}, {F}^{2}{\ensuremath{\Phi}}_{\mathrm{m}})}^{\frac{1}{2}}}$. $F$ is taken as a product of two-body factors, $\ensuremath{\Pi}{i<j}^{}[f({r}_{\mathrm{ij}})]$, the $f({r}_{\mathrm{ij}})$ embodying at least the correlations arising from strong short-range repulsions; and ${{\ensuremath{\Phi}}_{\mathrm{m}}}$, a set of independent-particle-model functions which might serve as a basis in the absence of highly singular interactions. The matrix elements ${\mathfrak{H}}_{\mathrm{mn}}$, ${\mathfrak{N}}_{\mathrm{mn}}$ of the Hamiltonian and identity operators with respect to the correlated basis ${{\ensuremath{\Psi}}_{\mathrm{m}}}$ are evaluated by generalization of cluster-expansion techniques developed originally for the calculation of the partition function of a classical imperfect gas and adapted to the quantum-mechanical many-body problem by Iwamoto and Yamada. Contributions to our expansions for the ${\mathfrak{H}}_{\mathrm{mn}}$ are classified according to order of magnitude in the quantity $\ensuremath{\rho}\ensuremath{\omega}$, where $\ensuremath{\omega}=\ensuremath{\int}({f}^{2}(r)\ensuremath{-}1)d\mathrm{r}$ is a "correlation parameter" and $\ensuremath{\rho}$ the average particle density. To the extent that $|\ensuremath{\rho}\ensuremath{\omega}|$ is small, it provides a true expansion parameter. A detailed technical discussion of the cluster expansions is given, with special emphasis on their symmetry properties upon truncation. Although the many-body problem is soluble by simultaneous diagonalization of the matrices (${\mathfrak{H}}_{\mathrm{mn}}$) and (${\mathfrak{N}}_{\mathrm{mn}}$), it is preferable first to transform to an orthonormal basis of "unperturbed" correlated functions ${{\ensuremath{\Theta}}_{\mathrm{m}}}$ with respect to which the Hamiltonian and identity operators have the matrix representations (${H}_{\mathrm{mn}}$) and $({N}_{\mathrm{mn}})=({\ensuremath{\delta}}_{\mathrm{mn}})$, and then to attempt the separate diagonalization of (${H}_{\mathrm{mn}}$). In general, such a transformation is expedited by the L\"owdin procedure carried out to the desired order in $\ensuremath{\rho}\ensuremath{\omega}$. However, a Schmidt orthogonalization scheme might prove more advantageous in specific cases. Field-theoretic techniques (of both perturbative and nonperturbative nature) may be brought to bear in the diagonalization of (${H}_{\mathrm{mn}}$). The inverse of the unitary transformation which takes the basis ${{\ensuremath{\Phi}}_{\mathrm{m}}}$ into the basis ${{\ensuremath{\Theta}}_{\mathrm{m}}}$, applied to the given singular Hamiltonian, leads---for any $F$ in the class considered---to a non-singular "effective Hamiltonian." Following this formal cue, we set out to construct explicitly a second-quantized effective Hamiltonian which generates the same matrix elements, to prescribed order in $\ensuremath{\rho}\ensuremath{\omega}$, when operating between independent-particle kets of the model occupation-number representation, as does the given Hamiltonian in the configuration-space representation, between the corresponding correlated functions of the L\"owdin basis. The quasiparticles described by this effective Hamiltonian interact via well-behaved two-body, three-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}, $Q$-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{} potentials, in contrast to the highly singular two-body interactions of the real particles. For small $|\ensuremath{\rho}\ensuremath{\omega}|$, the effective two-body forces predominate. An avenue to deeper understanding of the applicability of the nuclear-shell model and the dressed-quasiparticle theory of liquid ${\mathrm{He}}^{3}$ has been opened. This article treats only the formalism for a uniform extended medium. Numerical calculations, as well as the extension of the method to finite systems, will follow.

Physical Review•Jan 21, 1966•160 citations

Three-Body Clusters in Nuclear Matter

João da Providência, Carl M. Shakin

It is shown that starting with a correlated wave function for nuclear matter of a general form one may make direct correspondence to the conventional diagrammatic theory of Brueckner, Bethe, and Goldstone. In particular, in this work we show that the expression for the three-body-cluster energy is the same as that obtained by Bethe. Also, it is shown that the current approximation which neglects the potential for particle states in the diagrammatic approach corresponds to a neglect of a coupling between the two-body (Bethe-Goldstone) equation and the equation for the three-body-cluster wave function (Bethe-Faddeev). The theory presented here provides a clarification of the relation of the $n$-body-cluster wave functions and the structure of the wave function of the entire system.

Physical Review C•Jan 1, 1972•13 citations

Cluster Expansions for Correlated Wave Functions in Nuclear Physics

Carl M. Shakin

The cluster expansion for the energy for a system with Bijl-Dingle-Jastrow (BDJ) correlation structure is derived in a new form. It is shown how the theory may be related to the Brueckner-Bethe-Goldstone theory with a particular choice for the correlation operator (or wave matrix) appearing in the BDJ wave function. The three-body cluster terms are discussed. It is shown that these clusters contain a limited class of dispersive terms which may be included in the analysis of the two-body clusters. The three-body clusters are also shown to contain terms which correct for violations of the Pauli principle in the correlation operator.

Physical Review C•Sep 1, 1971•16 citations
Bethe's treatment of three-body clusters is extended to the entire linked-cluster perturbation series. The result is a completely rearranged expansion in terms of "compact clusters." The spatial correlations within these terms are all of quite short range. As a direct consequence, each of these new compact-cluster terms is proportional to \ensuremathκh\ensuremath′\ensuremath−1{\ensuremath{\kappa}}^{{h}^{\ensuremath{'}}\ensuremath{-}1}\ensuremathκh\ensuremath′\ensuremath−1, where h\ensuremath′{h}^{\ensuremath{'}}h\ensuremath′ is the number of momenta inside the Fermi sea which can be summed independently after allowing for momentum conservation. The "small parameter" of the expansion is \ensuremathκ=\ensuremathρ\ensuremath∫∣\ensuremathζ∣2d\ensuremathτ\ensuremath{\kappa}=\ensuremath{\rho}\ensuremath{\int}{|\ensuremath{\zeta}|}^{2}d\ensuremath{\tau}\ensuremathκ=\ensuremathρ\ensuremath∫∣\ensuremathζ∣2d\ensuremathτ, where \ensuremathρ\ensuremath{\rho}\ensuremathρ is the ordinary nuclear density and \ensuremathζ=\ensuremathφ\ensuremath−\ensuremathψBG\ensuremath{\zeta}=\ensuremath{\varphi}\ensuremath{-}{\ensuremath{\psi}}_{\mathrm{BG}}\ensuremathζ=\ensuremathφ\ensuremath−\ensuremathψBG​ is the "wound" in the correlated two-body (bethe-Goldstone) wave function. In nuclear matter this \ensuremathκ\ensuremath{\kappa}\ensuremathκ is of order 10%. This result requires that the single-particle potentials be defined in terms of a certain subset of all the self-energy insertions. The allowed insertions are those which can be evaluated entirely on the energy shell, by means of a "generalized time-ordering" factorization. The totality of these insertions is said to constitute an "on-energy-shell mass operator" Mon{M}^{\mathrm{on}}Mon. The resulting occupied-state potentials are quantitatively very similar to those of most previous versions of nuclear-matter theory. However, the present intermediate-state potentials are small enough (being only of order +1 MeV) to be safely ignored in practical calculations. This simplification is offset by the need for a separate calculation of the three-body clusters. These single-particle potentials are "physically meaningful" in the sense that they lead to an optimum treatment of the short-range correlations. The bulk properties of nuclear systems are determined almost entirely by these short-range correlations. It is therefore proposed that similar potentials be used in a many-body theory of finite nuclei. True occupation numbers occur in a simple and natural way throughout the entire expansion. The formal and practical consequences of this feature are carefully examined. It is argued that this is significant for a theory of finite nuclei. Comparisons are made with several different "renormalized" formulations of quantum statistical mechanics, and also with the Green's-function theory of nuclear matter. The present formulation tends to artificially suppress the long-range correlations (e.g., pairing correlations) between particles near the Fermi surface. This defect can be eliminated by using a degenerate analog of the Goldstone expansion.
The Goldstone expansion is rederived by elementary time-independent methods, starting from Brillouin-Wigner (BW) perturbation theory. Interaction energy terms $\ensuremath{\Delta}E$ are expanded out of the BW energy denominators, and the series is then rearranged to obtain the linked-cluster result. Similar algebraic methods lead to the linked expansions for the total wave function (Hugenholtz) and the expectation value of a general operator (Thouless). Starting again with a degenerate version of BW perturbation theory, these methods are used to obtain the Bloch-Horowitz energy expansion, as well as the corresponding wave function, expectation-value, and transition-amplitude expansions. A "reduced" form of the Block-Horowitz expansion is described, and also a "completely linked" version. The latter is suggested as a tool for investigating superfiuid phenomena in nuclear matter, and for establishing contact with the Landau theory of Fermi liquids. The physical interpretation of these expansions is carefully studied, especially with regard to nuclear applications, to determine how they handle such "physical" features as antisymmetry, self-energy effects, wave-function renormalizations, and the distinction between "true" and "model" single-particle occupation probabilities. The problem of a correct theoretical definition for the shell-model potential is carefully examined, and a specific theory is presented. These expansions are seen to form a convenient and very powerful set of tools for studying the structure of actual nuclei.
An approximation method developed previously to deal with many particles in strong interaction is examined in further detail. It is shown that the series giving the interaction energy is a development in a sequence of linked or irreducible cluster terms each of which gives a contribution to the energy proportional to the total number of particles. Consequently the convergence of the expansion is independent of the total number of particles. The origin of this simple feature is illustrated by showing that a similar situation exists in the expansion of standard perturbation theory. The numerical convergence of the expansion is quantitatively discussed for the nuclear problem where it is shown that the correction arising from the first cluster term involving three particles is less than the leading term by a factor of about ${10}^{\ensuremath{-}4}$. The smallness of the correction is largely a result of the action of the exclusion principle.
The usual perturbation treatment of the many-body problem is modified by eliminating the two-body interaction and introducing in its place a ttt matrix. The Watson, Brueckner, and Bethe (W.B.B.) treatment of the many-body problem is shown to be a special case of the modified perturbation treatment. The modified perturbation treatment involves no assumption except that of supposing the convergence of the perturbation expansion. Thus to each order in the expansion the results must be meaningful and no unphysical "unlinked cluster" terms can appear. Within the framework of the modified perturbation treatment it is easy to interpret the concepts of the Pauli principle for intermediate states and the self-consistent choice of the comparison potential which play such an important role in the W.B.B. treatment. This work is similar to that of J. Goldstone in that we relate the W.B.B. treatment to a perturbation treatment of the many-body problem. We differ from Goldstone in that we use the usual time-independent perturbation theory instead of the time-dependent one. The time-independent approach permits a direct comparison with the W.B.B. treatment and facilitates discussion of the Pauli principle for intermediate states and of the self-consistent potential.
The method of correlated basis functions provides a means of systematically exploiting our intuitive understanding of real physical systems. Here we seek especially to develop a formalism for describing systems of strongly interacting fermions, for example, nuclei and liquid He3{\mathrm{He}}^{3}He3, using a complete set of correlated functions ${\ensuremath{\Psi}}_{\mathrm{m}}=\frac{F{\ensuremath{\Phi}}_{\mathrm{m}}}{{({\ensuremath{\Phi}}_{\mathrm{m}}, {F}^{2}{\ensuremath{\Phi}}_{\mathrm{m}})}^{\frac{1}{2}}}$. FFF is taken as a product of two-body factors, $\ensuremath{\Pi}{i<j}^{}[f({r}_{\mathrm{ij}})]$, the f(rij)f({r}_{\mathrm{ij}})f(rij​) embodying at least the correlations arising from strong short-range repulsions; and ${{\ensuremath{\Phi}}_{\mathrm{m}}}$, a set of independent-particle-model functions which might serve as a basis in the absence of highly singular interactions. The matrix elements Hmn{\mathfrak{H}}_{\mathrm{mn}}Hmn​, Nmn{\mathfrak{N}}_{\mathrm{mn}}Nmn​ of the Hamiltonian and identity operators with respect to the correlated basis ${{\ensuremath{\Psi}}_{\mathrm{m}}}$ are evaluated by generalization of cluster-expansion techniques developed originally for the calculation of the partition function of a classical imperfect gas and adapted to the quantum-mechanical many-body problem by Iwamoto and Yamada. Contributions to our expansions for the Hmn{\mathfrak{H}}_{\mathrm{mn}}Hmn​ are classified according to order of magnitude in the quantity $\ensuremath{\rho}\ensuremath{\omega}$, where $\ensuremath{\omega}=\ensuremath{\int}({f}^{2}(r)\ensuremath{-}1)d\mathrm{r}$ is a "correlation parameter" and $\ensuremath{\rho}$ the average particle density. To the extent that $|\ensuremath{\rho}\ensuremath{\omega}|$ is small, it provides a true expansion parameter. A detailed technical discussion of the cluster expansions is given, with special emphasis on their symmetry properties upon truncation. Although the many-body problem is soluble by simultaneous diagonalization of the matrices (Hmn{\mathfrak{H}}_{\mathrm{mn}}Hmn​) and (Nmn{\mathfrak{N}}_{\mathrm{mn}}Nmn​), it is preferable first to transform to an orthonormal basis of "unperturbed" correlated functions ${{\ensuremath{\Theta}}_{\mathrm{m}}}$ with respect to which the Hamiltonian and identity operators have the matrix representations (Hmn{H}_{\mathrm{mn}}Hmn​) and $({N}_{\mathrm{mn}})=({\ensuremath{\delta}}_{\mathrm{mn}})$, and then to attempt the separate diagonalization of (Hmn{H}_{\mathrm{mn}}Hmn​). In general, such a transformation is expedited by the L\"owdin procedure carried out to the desired order in $\ensuremath{\rho}\ensuremath{\omega}$. However, a Schmidt orthogonalization scheme might prove more advantageous in specific cases. Field-theoretic techniques (of both perturbative and nonperturbative nature) may be brought to bear in the diagonalization of (Hmn{H}_{\mathrm{mn}}Hmn​). The inverse of the unitary transformation which takes the basis ${{\ensuremath{\Phi}}_{\mathrm{m}}}$ into the basis ${{\ensuremath{\Theta}}_{\mathrm{m}}}$, applied to the given singular Hamiltonian, leads---for any FFF in the class considered---to a non-singular "effective Hamiltonian." Following this formal cue, we set out to construct explicitly a second-quantized effective Hamiltonian which generates the same matrix elements, to prescribed order in $\ensuremath{\rho}\ensuremath{\omega}$, when operating between independent-particle kets of the model occupation-number representation, as does the given Hamiltonian in the configuration-space representation, between the corresponding correlated functions of the L\"owdin basis. The quasiparticles described by this effective Hamiltonian interact via well-behaved two-body, three-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}, QQQ-body, \ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{}\ifmmode\cdot\else\textperiodcentered\fi{} potentials, in contrast to the highly singular two-body interactions of the real particles. For small $|\ensuremath{\rho}\ensuremath{\omega}|$, the effective two-body forces predominate. An avenue to deeper understanding of the applicability of the nuclear-shell model and the dressed-quasiparticle theory of liquid He3{\mathrm{He}}^{3}He3 has been opened. This article treats only the formalism for a uniform extended medium. Numerical calculations, as well as the extension of the method to finite systems, will follow.
It is shown that starting with a correlated wave function for nuclear matter of a general form one may make direct correspondence to the conventional diagrammatic theory of Brueckner, Bethe, and Goldstone. In particular, in this work we show that the expression for the three-body-cluster energy is the same as that obtained by Bethe. Also, it is shown that the current approximation which neglects the potential for particle states in the diagrammatic approach corresponds to a neglect of a coupling between the two-body (Bethe-Goldstone) equation and the equation for the three-body-cluster wave function (Bethe-Faddeev). The theory presented here provides a clarification of the relation of the nnn-body-cluster wave functions and the structure of the wave function of the entire system.
The cluster expansion for the energy for a system with Bijl-Dingle-Jastrow (BDJ) correlation structure is derived in a new form. It is shown how the theory may be related to the Brueckner-Bethe-Goldstone theory with a particular choice for the correlation operator (or wave matrix) appearing in the BDJ wave function. The three-body cluster terms are discussed. It is shown that these clusters contain a limited class of dispersive terms which may be included in the analysis of the two-body clusters. The three-body clusters are also shown to contain terms which correct for violations of the Pauli principle in the correlation operator.