The variational approximation to the thermal properties of a quantum-mechanical system proposed by Giachetti and Tognetti [Phys. Rev. Lett. 55, 912 (1985); Phys. Rev. B 33, 7647 (1986)], Feynman and Kleinert [Phys. Rev. A 34, 5080 (1986)], and others is taken as the starting point, and a systematic perturbation series is developed for both the free energy and the particle density. Generalizations of the usual Feynman propagator appear naturally in the formalism. The anharmonic oscillator is treated as an example; the second-order correction to the free energy and the first-order correction to the particle density are evaluated explicitly and shown to lead to significant improvements over the variational solution.
New forms of the linked cluster expansion for the ground-state energy of a normal spherical system are given in terms of the exact single-particle propagator and of the associated self-energy operator. These forms are shown to be stationary with respect to variations of the latter quantities. A new concise proof of the Hugenholtz-Van Hove theorem on single-particle energies is then possible. The variational principle also explains the insensitivity of recent binding energy calculations to the choice of single-particle energies. It is shown how all results may be extended to systems with anisotropic features.
A variational perturbation theory based on the functional integral approach is formulated for many-particle systems. Using the variational action obtained through Jensen-Peierls' inequality, a perturbative expansion scheme for the thermodynamic potential is established. A modified Wick's theorem is obtained for the variational perturbation expansions. This theorem allows one to carry out systematic calculations of higher order terms without worrying about the double counting problem. A model numerical calculation was carried out on a nucleon gas system interacting through the Yukawa-type potential to test the efficiency of 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.
Variation principles for obtaining lower bounds by means of density matrices are applied to a number of many-body problems. This is done by means of restrictions derived from the study of exacly solvable models. Two classes of problems are considered: (a) obtaining a lower bound to the ground-state energy of a system of particles, and (b) obtaining a lower bound to the Helmholtz free energy. For the first case an exact lower bound is obtained for the ground-state energy of the particle-conserving Bogoliubov Hamiltonian. Of particular significance in the nonzero temperature case is a rigorous lower bound to the free energy of the Ising model with small external field at low temperatures.