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

Renormalization and Statistical Mechanics in Many-Particle Systems. II. Statistical Perturbation Method

Duane C. Wallace

Published Dec 2, 1966• 21 citations
Publisher

Abstract

A general nonperturbation problem, described by the Hamiltonian $\mathcal{H}$, is considered. Approximate creation operators ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ are defined as operators which satisfy the Hamiltonian commutator equations $[\mathcal{H},{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}]={\ensuremath{\omega}}_{i}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}+{{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, where ${\ensuremath{\omega}}_{i}$ are the creation energies, and the remainder operators ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ are small in the sense that statistical averages involving the ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ are small. A zeroth approximation is to neglect the ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$. In this case, a straight-forward derivation provides general relations between statistical averages such as $〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Omega}〉$ and $〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Omega}\ifmmode\pm\else\textpm\fi{}\ensuremath{\Omega}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}〉$, where $\ensuremath{\Omega}$ is any operator. The usual boson and fermion occupation numbers are a trivial result of these zeroth-order relations. Corrections to the general zeroth-order relations arise from the ${{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$; such corrections are considered in the spirit of perturbation theory. An explicit first-order correction is obtained for the averages $〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}{\ensuremath{\theta}}_{i}〉$, and this correction is related to a self-consistent energy-renormalization procedure. The present method is similar to the method of thermodynamic Green's functions in that the derivations given here do not require evaluation of any state vectors or of the partition function. The general results are applied to the Heisenberg ferromagnet problem. At low temperatures, the zeroth-order relations give the Bloch spin-wave results, and an approximate evaluation of the first-order contribution gives the leading term in Dyson's ${T}^{4}$ correction to the spontaneous magnetization. For arbitrary temperatures, the zeroth-order relations give the Tyablikov equations for $〈{S}^{z}〉$, and Callen's results are obtained from the first-order corrections. These examples illustrate the simplicity of calculating statistical averages of functions of the ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, ${\ensuremath{\theta}}_{i}$ operators by the present method.

DOI
10.1103/physrev.152.261

Related papers

Renormalization and Statistical Mechanics in Many-Particle Systems. I. Hamiltonian Perturbation Method

Duane C. Wallace

A general system composed of many weakly interacting bosons and/or fermions is considered. The object is to develop a procedure for renormalizing the single-particle creation operators, so as to remove the interactions to successively higher orders of perturbation. The basic idea is that creation operators ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ must satisfy the Hamiltonian commutator equations $[\mathcal{H},{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}]={\ensuremath{\omega}}_{i}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, where $\mathcal{H}$ is the Hamiltonian and ${\ensuremath{\omega}}_{i}$ are the particle energies. For the many-weakly-interacting-particle system, the zeroth-order boson and fermion creation operators satisfy these equations to zeroth order. It is shown that if the creation operators are renormalized so as to satisfy the Hamiltonian commutator equations to order $m$, and also to satisfy the appropriate boson commutators and fermion anticommutators to order $m$, then the problem is solved to order $m$. In particular, the vectors formed by operating on the ground state with the renormalized creation operators, according to the usual boson and fermion occupation-number representation, are eigenvectors of $\mathcal{H}$ and are orthonormal, all to order $m$. A procedure is otained for finding the ($m+1$)-order contributions to the particle creation energies, in terms of the $m$-order operators. Explicit first-order calculations of these general results are provided for a system of bosons and a system of fermions, and these first-order results are shown to include similar results of Rayleigh-Schr\"odinger perturbation theory, the random-phase approximation, and the method of thermodynamic Green's functions. The problem of anharmonic lattice dynamics is studied in detail, and a method of undetermined coefficients is used to renormalize the phonon creation operators to first order. The phonon energies are calculated to second order, and this calculation shows that the interactions between renormalized phonons cannot be removed in second order. Statistical averages of the phonon energies give the energy shifts and lifetimes which have been calculated previously by various propagator techniques. In addition, the renormalized energies are used to calculate the temperature-dependent part of the Helmholtz free energy correct to second order. As a final example, electron-phonon interactions in a normal metal are studied. The electron and phonon creation operators are renormalized to first order, the particle energies are calculated to second order, and statistical averages of the particle energies recover the usual thermodynamic Green's-function results for energy shifts and lifetimes. These examples show the simplicity by which the renormalization procedure obtains a great amount of detailed information about the single-particle nature of a many-particle system.

Physical Review•Dec 2, 1966•31 citations

Self-consistent many body theory for the standard basis operator Green's functions

R. Micnas

We propose a self-consistent many body theory for the standard basis operator Green 's functicms and obtain an exact Lyson type matriz equation for the interacting mcmy levei systems. A zeroth order crpproximation, which neglects ali the damping effects, is investigated in detail for the anisotropic Heisenberg model, the isotropic quadvopolar system and the Hubbard model. In case of the anisotropic Heisenberg ferromagnet with both exchange and single -ion anisotropy the low-temperature renomnalization of the apin-waves for the uniaxial ordering agrees with the Bloch-Dyson theory. For the spin -1 easy plane ferromagnet, the criticai parameters for the phase transition at zero temperature are determined and compared with other theories. The elementary excitation spectrum of the spin -1 isotropic quadrupolar system is calculated and compared with the random phase approximation and Callen's like deccupling schemes. Finally, the theory is applied to the studW of the single -particle excitation spectrum of the Hubbard model.

•Dec 1, 1980•0 citations

New approach to perturbation theory of many-particle systems

Colin John Lambert, W. E. Hagston

International Journal of Theoretical Physics•Feb 1, 1984•0 citations

Operator Method for Nonperturbative Calculation of the Thermodynamic Values in Quantum Statistics

И. Д. Феранчук, A. A. Ivanov

Operator method and cumulant expansion are used for nonperturbative calculation of the partition function and the free energy in quantum statistics. It is shown for some model systems that the zeroth approximation of the method proposed interpolates the thermodynamics values with rather good accuracy in the entire range of both the Hamiltonian parameters and temperature. The regular procedure for calculation of the corrections to the zeroth approximation is also considered.

•Sep 21, 2003•0 citations

Asymptotic nature of perturbation theory in the many-Fermion problem

R.J. Rennie

Physics Letters A•Feb 1, 1980•0 citations

Functional-derivative study of the Hubbard model. III. Fully renormalized Green's function

Tadashi Arai, Morrel H. Cohen

The functional-derivative method of calculating the Green's function developed earlier for the Hubbard model is generalized and used to obtain a fully renormalized solution. Higher-order functional derivatives operating on the basic Green's functions, $G$ and $\ensuremath{\Gamma}$, are all evaluated explicitly, thus making the solution applicable to the narrow-band region as well as the wide-band region. Correction terms $\ensuremath{\Phi}$ generated from functional derivatives of equal-time Green's functions of the type $\frac{{\ensuremath{\delta}}^{n}〈N〉}{\ensuremath{\delta}{\ensuremath{\epsilon}}^{n}}$, etc., with $n\ensuremath{\ge}2$. It is found that the $\ensuremath{\Phi}'\mathrm{s}$ are, in fact, renormalization factors involved in the self-energy $\ensuremath{\Sigma}$ and that the structure of the $\ensuremath{\Phi}'\mathrm{s}$ resembles that of $\ensuremath{\Sigma}$ and contains the same renormalization factors $\ensuremath{\Phi}$. The renormalization factors $\ensuremath{\Phi}$ are shown to satisfy a set of equations and can be evaluated self-consistently. In the presence of the $\ensuremath{\Phi}'\mathrm{s}$, all difficulties found in the previous results (papers I and II) are removed, and the energy spectrum $\ensuremath{\omega}$ can now be evaluated for all occupations $n$. The Schwinger relation is the only basic relation used in generating this fully self-consistent Green's function, and the Baym-Kadanoff continuity condition is automatically satisfied.

Physical review. B, Condensed matter•Apr 15, 1980•7 citations
A general nonperturbation problem, described by the Hamiltonian H\mathcal{H}H, is considered. Approximate creation operators \ensuremathθi\ifmmode†\else\textdagger\fi{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremathθi​\ifmmode†\else\textdagger\fi are defined as operators which satisfy the Hamiltonian commutator equations [H,\ensuremathθi\ifmmode†\else\textdagger\fi]=\ensuremathωi\ensuremathθi\ifmmode†\else\textdagger\fi+Ri\ifmmode†\else\textdagger\fi[\mathcal{H},{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}]={\ensuremath{\omega}}_{i}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}+{{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}[H,\ensuremathθi​\ifmmode†\else\textdagger\fi]=\ensuremathωi​\ensuremathθi​\ifmmode†\else\textdagger\fi+Ri​\ifmmode†\else\textdagger\fi, where \ensuremathωi{\ensuremath{\omega}}_{i}\ensuremathωi​ are the creation energies, and the remainder operators Ri\ifmmode†\else\textdagger\fi{{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}Ri​\ifmmode†\else\textdagger\fi are small in the sense that statistical averages involving the Ri\ifmmode†\else\textdagger\fi{{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}Ri​\ifmmode†\else\textdagger\fi are small. A zeroth approximation is to neglect the Ri\ifmmode†\else\textdagger\fi{{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}Ri​\ifmmode†\else\textdagger\fi. In this case, a straight-forward derivation provides general relations between statistical averages such as 〈\ensuremathθi\ifmmode†\else\textdagger\fi\ensuremathΩ〉〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Omega}〉〈\ensuremathθi​\ifmmode†\else\textdagger\fi\ensuremathΩ〉 and 〈\ensuremathθi\ifmmode†\else\textdagger\fi\ensuremathΩ\ifmmode±\else\textpm\fi\ensuremathΩ\ensuremathθi\ifmmode†\else\textdagger\fi〉〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremath{\Omega}\ifmmode\pm\else\textpm\fi{}\ensuremath{\Omega}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}〉〈\ensuremathθi​\ifmmode†\else\textdagger\fi\ensuremathΩ\ifmmode±\else\textpm\fi\ensuremathΩ\ensuremathθi​\ifmmode†\else\textdagger\fi〉, where \ensuremathΩ\ensuremath{\Omega}\ensuremathΩ is any operator. The usual boson and fermion occupation numbers are a trivial result of these zeroth-order relations. Corrections to the general zeroth-order relations arise from the Ri\ifmmode†\else\textdagger\fi{{R}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}Ri​\ifmmode†\else\textdagger\fi; such corrections are considered in the spirit of perturbation theory. An explicit first-order correction is obtained for the averages 〈\ensuremathθi\ifmmode†\else\textdagger\fi\ensuremathθi〉〈{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}{\ensuremath{\theta}}_{i}〉〈\ensuremathθi​\ifmmode†\else\textdagger\fi\ensuremathθi​〉, and this correction is related to a self-consistent energy-renormalization procedure. The present method is similar to the method of thermodynamic Green's functions in that the derivations given here do not require evaluation of any state vectors or of the partition function. The general results are applied to the Heisenberg ferromagnet problem. At low temperatures, the zeroth-order relations give the Bloch spin-wave results, and an approximate evaluation of the first-order contribution gives the leading term in Dyson's T4{T}^{4}T4 correction to the spontaneous magnetization. For arbitrary temperatures, the zeroth-order relations give the Tyablikov equations for 〈Sz〉〈{S}^{z}〉〈Sz〉, and Callen's results are obtained from the first-order corrections. These examples illustrate the simplicity of calculating statistical averages of functions of the \ensuremathθi\ifmmode†\else\textdagger\fi{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}\ensuremathθi​\ifmmode†\else\textdagger\fi, \ensuremathθi{\ensuremath{\theta}}_{i}\ensuremathθi​ operators by the present method.
A general system composed of many weakly interacting bosons and/or fermions is considered. The object is to develop a procedure for renormalizing the single-particle creation operators, so as to remove the interactions to successively higher orders of perturbation. The basic idea is that creation operators ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ must satisfy the Hamiltonian commutator equations $[\mathcal{H},{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}]={\ensuremath{\omega}}_{i}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, where H\mathcal{H}H is the Hamiltonian and ${\ensuremath{\omega}}_{i}$ are the particle energies. For the many-weakly-interacting-particle system, the zeroth-order boson and fermion creation operators satisfy these equations to zeroth order. It is shown that if the creation operators are renormalized so as to satisfy the Hamiltonian commutator equations to order mmm, and also to satisfy the appropriate boson commutators and fermion anticommutators to order mmm, then the problem is solved to order mmm. In particular, the vectors formed by operating on the ground state with the renormalized creation operators, according to the usual boson and fermion occupation-number representation, are eigenvectors of H\mathcal{H}H and are orthonormal, all to order mmm. A procedure is otained for finding the (m+1m+1m+1)-order contributions to the particle creation energies, in terms of the mmm-order operators. Explicit first-order calculations of these general results are provided for a system of bosons and a system of fermions, and these first-order results are shown to include similar results of Rayleigh-Schr\"odinger perturbation theory, the random-phase approximation, and the method of thermodynamic Green's functions. The problem of anharmonic lattice dynamics is studied in detail, and a method of undetermined coefficients is used to renormalize the phonon creation operators to first order. The phonon energies are calculated to second order, and this calculation shows that the interactions between renormalized phonons cannot be removed in second order. Statistical averages of the phonon energies give the energy shifts and lifetimes which have been calculated previously by various propagator techniques. In addition, the renormalized energies are used to calculate the temperature-dependent part of the Helmholtz free energy correct to second order. As a final example, electron-phonon interactions in a normal metal are studied. The electron and phonon creation operators are renormalized to first order, the particle energies are calculated to second order, and statistical averages of the particle energies recover the usual thermodynamic Green's-function results for energy shifts and lifetimes. These examples show the simplicity by which the renormalization procedure obtains a great amount of detailed information about the single-particle nature of a many-particle system.
We propose a self-consistent many body theory for the standard basis operator Green 's functicms and obtain an exact Lyson type matriz equation for the interacting mcmy levei systems. A zeroth order crpproximation, which neglects ali the damping effects, is investigated in detail for the anisotropic Heisenberg model, the isotropic quadvopolar system and the Hubbard model. In case of the anisotropic Heisenberg ferromagnet with both exchange and single -ion anisotropy the low-temperature renomnalization of the apin-waves for the uniaxial ordering agrees with the Bloch-Dyson theory. For the spin -1 easy plane ferromagnet, the criticai parameters for the phase transition at zero temperature are determined and compared with other theories. The elementary excitation spectrum of the spin -1 isotropic quadrupolar system is calculated and compared with the random phase approximation and Callen's like deccupling schemes. Finally, the theory is applied to the studW of the single -particle excitation spectrum of the Hubbard model.
The functional-derivative method of calculating the Green's function developed earlier for the Hubbard model is generalized and used to obtain a fully renormalized solution. Higher-order functional derivatives operating on the basic Green's functions, GGG and $\ensuremath{\Gamma}$, are all evaluated explicitly, thus making the solution applicable to the narrow-band region as well as the wide-band region. Correction terms $\ensuremath{\Phi}$ generated from functional derivatives of equal-time Green's functions of the type $\frac{{\ensuremath{\delta}}^{n}〈N〉}{\ensuremath{\delta}{\ensuremath{\epsilon}}^{n}}$, etc., with $n\ensuremath{\ge}2$. It is found that the $\ensuremath{\Phi}'\mathrm{s}$ are, in fact, renormalization factors involved in the self-energy $\ensuremath{\Sigma}$ and that the structure of the $\ensuremath{\Phi}'\mathrm{s}$ resembles that of $\ensuremath{\Sigma}$ and contains the same renormalization factors $\ensuremath{\Phi}$. The renormalization factors $\ensuremath{\Phi}$ are shown to satisfy a set of equations and can be evaluated self-consistently. In the presence of the $\ensuremath{\Phi}'\mathrm{s}$, all difficulties found in the previous results (papers I and II) are removed, and the energy spectrum $\ensuremath{\omega}$ can now be evaluated for all occupations nnn. The Schwinger relation is the only basic relation used in generating this fully self-consistent Green's function, and the Baym-Kadanoff continuity condition is automatically satisfied.