An extension is presented of the slightly refined Matsubara method for the computation of the grand partition function for quantum statistical averages. This computation is extended for the case of the electron gas, from which the results of Montroll and Ward are obtained. A slight
In Chapter 9 we applied the second quantization to the Hartree–Fock method. In this chapter we provide the basic aspects of the thermal Green function method, which is convenient for systematic perturbative expansion beyond the mean-field approximation. Several problems in quantum statistical mechanics can be elegantly and compactly described in terms of the so-called Matsubara Green functions (Matsubara (1955); Abrikosov et al. (1963)). Hence, if on one hand, the density matrix description, often used in the previous chapters, has a more direct physical interpretation, on the other hand the Green function method is far more powerful and takes advantage of a handy diagrammatic perturbation expansion. In the following chapters, we rely on the Green function for a more complete treatment of some many-body problems. To be self contained we introduce here this technical chapter with a short presentation of the by now standard Green function method. For a comprehensive treatment we refer the reader to the many books available and in particular to Abrikosov et al. (1963). The Matsubara Green function The main aim of quantum statistical mechanics is the evaluation of the partition function (see Eqs. (6.31) and (6.32)), e.g. in the grand canonical ensemble In the case of quantum gases one can evaluate Ƶ exactly as we did in Chapter 7. In general, the evaluation of the partition function is a difficult task and it is useful to use perturbation theory to obtain approximate solutions. To set up a perturbation theory, one begins by splitting the Hamiltonian as H = H 0 + H i, where H 0 represents the “easy” part, which can be solved exactly. In the second step, one could naturally think of expanding the exponential of the interacting part H i in order to write the partition function as a perturbative series. However, such a procedure requires extreme care due to the non-commutative nature of operators in quantum mechanics. The way of solving this problem is via the “ time ” ordering of operators in each term of the perturbative expansion. It can be shown, as we shall see in the present chapter, that the quantum statistical average of any time-ordered sequence of operators with respect to the unperturbed Hamiltonian H 0 can be expressed in terms of elementary quantum statistical averages of two field operators. Such elementary averages are called Green functions and in this section we introduce the basic definitions concerning them.
A systematic generalization of the Mayer cluster integral theory has been developed to deal with the quantum statistics of interacting particles. The grand partition function appears in a natural way and the cluster integrals are integrals over propagators which are derived from the Green's function solution of the Bloch equation (which follows from the Schroedinger equation by replacing it/ℏ by β = 1/kT). Every cluster integral can be represented by a hybrid of a Mayer graph and a Feynman diagram in (β, r) space. The generalization of classical ring cluster integrals has been analyzed. It is shown that in the case of the electron gas the classical limit of the contribution of these integrals to the grand partition function yields the Debye-Huckel theory while the low temperature limit leads to the Gell-Mann—Brueckner equation for the correlation energy of the ground state. A prescription is given for the construction of the cluster integral associated with any given diagram.
A new method of calculating the grand partition function of many-body system is developed, adopting extensively the techniques of calculus in quantum field theory. It is shown that the grand partition function, which is a trace of the density matrix expressed in terms of field operators, can be evaluated in a way almost parallel with the evaluation of the vacuum expectation value of the S-matrix in quantum field theory, provided that appropriate modifications in notation and definitions are made. As an example, the method is applied to electron-phonon system. Further, basing on this new formalism, various non-perturbational methods are discussed.
A general method of calculation is described for quantum statistical mechanics. It is based on a simplification of the Laplace transform of the density matrix which follows from a theorem due to Hugenholtz. The basic result is that an element of the density matrix can be written as a sum over graphs, with the contribution of each graph factored into contributions from connected or linked graphs. Applied to the grand partition function, the exponential formula of Bloch and DeDominicis is obtained in a simple way. A similar formula is then derived for the canonical ensemble for the case of a nondegenerate gas. In this way the familiar result of Uhlenbeck and Beth is obtained for the second virial coefficient. Techniques are also introduced for evaluating ensemble averages of operators. In this connection, some care must be exercised in the case of diagonal operators. Finally, these methods are used to calculate the pair-correlation function for a system of fermions interacting through short-range forces.