A general nonperturbation problem, described by the Hamiltonian
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
T4 correction to the spontaneous magnetization. For arbitrary temperatures, the zeroth-order relations give the Tyablikov equations for
〈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{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, ${\ensuremath{\theta}}_{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 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
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.
The theory of many-Boson systems was considered. A canonical transformation of the creation and annihilation operators was introduced in order to eliminate the singleparticle zero-momentam state and allow the application of field-theoretic penturbation techniques. These techniques were used to give a simpler derivation of a result due to Belyaev and to prove a theorem about the excitation spectrum of the system. General techniques for calculating the Green's functions are discussed. The excitation spectrum was calculated to second order for small wave number for the case of a repulsive square-well potential.
As a continuation of a recent paper, using density and current density as coordinates for interacting bosons at absolute zero, this work deals with the perturbation theory. The zerothorder Hamiltonian is the one shown to yield the Bogoliubov energy spectrum. First- or secondorder perturbation calculations have been made of the three-phonon and the four-phonon vertices for the energy corrections of the low-lying excited states. The results in the long-wavelength limit differ, although only slightly, from the quantum hydrodynamic calculations, and indicate that the perturbed excitation energy ${\ensuremath{\epsilon}}_{\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}}$ is an odd-power series in
k if the Fourier coefficients of the two-body potential concerned can be approximated by a constant or an evenpower series in
k. However, whether this remains so in all higher orders of various perturbations is an open question. Two methods were used in the calculation. The first involves functional representations of the coordinates in the ${\ensuremath{\rho}}_{\stackrel{\ensuremath{\rightarrow}}{\mathrm{k}}}$ representation recently reported. The second method is new, and employs "occupation-number" representations of the current algebra by introducing creation and annihilation operators for "phonons" from density fluctuations. These latter representations enable us to benefit from the advantages of the canonical field theory.