In order to go over from a perturbation expansion of the grand partition function (the unlinked-cluster expansion) to an expansion of the thermodynamic potential (the linked-cluster expansion) in powers of the interaction, it is necessary to treat carefully those terms in which creation (or annihilation) operators for the same state occur twice or more. The unlinked- and linked-cluster expansions for a system of fermions are here shown to be equivalent by a direct comparison of the terms which occur in each. The relation between the two expansions is illustrated by the example of a system of fermions interacting only with a single-particle potential.
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.
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.
We developed a reasonably detailed and comprehensive understanding of boson and fermion statistics from the previous chapter. While this is interesting physics in itself, more important for us is an understand of how the symmetry properties of these particles impact on the entropy associated with them. In this chapter, we start by considering situations where the number of particles in a system is not fixed. Next, the theory developed is used to derive quantum statistics again, but this time without recourse to various approximations. Access to the requisite probability functions enables the Gibbs entropy to be calculated for fermions and bosons .