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

Journal of Mathematical Physics

Approximate Functional Integral Methods in Statistical Mechanics. I. Moment Expansions

Armand Siegel, Terence Burke

Published Nov 1, 1972• 13 citations
Publisher

Abstract

In this paper, four distinct ideas are combined, which under a wide range of circumstances can give very rapidly converging series expansions for functional integrals. (1) Expansion of the functional being integrated in functional Taylor series. In the familiar case arising in quantum statistical mechanics, that of the Wiener integral of exp[−∫V(x(t))dt],V(x) being the perturbing potential, this is equivalent to expanding the characteristic functional of the probability functional of V(x(t)) in central moments of V(x(t)). The lowest-order term of the series is the approximation obtained by Feynman and Hibbs through a variational method. (2) Transfer of the harmonic term of the potential, when the functional integral is the quantum-statistical density matrix (Green's function of the Bloch equation), to the weighting function. This transforms the functional integral from a Wiener to an Uhlenbeck-Ornstein integral. The formal expressions for the terms of the expansion are somewhat more complicated, but they can be worked out, and the result is a great improvement in the speed of convergence of the series with decreasing temperature and/or decreasing relative magnitude of the anharmonic part of the potential. (3) ``Reservation of variables'' in the integration. This amounts to breaking the averaging process down into an average over subsets of the distributed random function (conditional average), followed by an averaging of these averages. Any step of this kind (it may be repeated within the subsets, etc.) gives an improvement of accuracy. (4) When the quantity being evaluated through functional integration is the partition function, the device introduced by Feynman and Hibbs, of interchanging the functional integration with the integration of the Green's function over the equated initial and final configuration-space points, may be combined with the above techniques. This eliminates one integration in the terms of the expansion and seems to improve accuracy at the same time. The general series obtained is correlated with more conventional operator techniques of quantum-mechanical perturbation theory, in order to answer the perennial question, does the path-integral method bring with it anything that could not be derived by other methods? It is in some sense a Feynman-Dyson expansion of the Green's function, but one that is further modified mathematically in a way characteristic only of the path-integral point of view, and which, moreover, improves its accuracy. It thus appears unlikely that the result is merely one of standard type disguised as a functional integral result. Sample numerical calculations are given to assay the accuracy of the methods, which are shown to compare very favorably with the traditional approximation of finite subdivision of the time interval.

DOI
10.1063/1.1665892

Related papers

Thermal Green functions

C. Di Castro, Roberto Raimondi

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.

Cambridge University Press eBooks•Aug 27, 2015•0 citations

The Feynman path integral

Mark E. Tuckerman

Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations

•Aug 1, 2023•0 citations

Space-Temperature Correlations in Quantum-Statistical Mechanics

William Kunkin, Jerome K. Percus

The method of functional differentiation, used in classical statistical mechanics to obtain approximate integral equations for the pair distribution function, is extended to quantum systems obeying Maxwell-Boltzmann statistics. The grand partition function is written as a path integral, and space-temperature distributions are generated from it by successive functional differentiation. Distributions can be expanded in powers of an external potential via a functional Taylor series. When the series is truncated so that no distributions higher than second order enter and the external potential is specialized to one arising from a fixed particle in the system, coupled integral equations result for the two kinds of pair distributions. The first-order Ursell function calculated from these equations yields the Montroll-Ward ring summation for the grand potential. This approximate Ursell function determines also the Fourier transform of the momentum density (useful in calculating the p = 0 occupation number). Although divergences appear in the low-temperature limit, they can be removed by a simple modification of the independent functional. An extension to Fermi-Dirac and Böse-Einstein statistics is indicated at appropriate places in the text.

Journal of Mathematical Physics•Mar 1, 1970•0 citations

Series expansions and their summation

R. J. Rivers

Let us recapitulate. Although the form (4.32) for the path integral was motivated by analogy with finite-dimensional integrals, in no way is it a definition. Rather, as we have stressed repeatedly, it is an algorithm for handling the Feynman perturbation series. Strong coupling apart, all the manipulations of the path integral that we have performed (integration by parts, the introduction of auxiliary fields, the inclusion of background fields) stand in one-to-one correspondence to rearrangements of the perturbation series. In fact, path integrals are unnecessary in the derivation of expressions like (4.50). See Fried (1972), for example. The path integrals are as rigorous – or non-rigorous – as the series, irrespective of whether they can be given proper mathematical definition or not. At the same time, the analogy with finite-dimensional integrals suggests approaches that go beyond this purely operational stance. In particular, because of the dominant role played by the classical action in the path integral formalism we would like to use it to establish the relation of quantum dynamics to classical dynamics. This is the most compelling reason for introducing the path integral formalism, requiring a distinct relaxation in our treatment of it. By formal analogy with finite-dimensional saddle-point and stationary-phase methods, elevating the classical action to the phase of the integrand promotes series expansions in ħ as the natural way to proceed.

Cambridge University Press eBooks•Nov 12, 1987•0 citations

Summation of Perturbation Series in Quantum Mechanics

Anil C. Rae

This is the first in a series of papers. It deals with the linked-cluster summation of the perturbation series in quantum mechanics and quantum statistical physics. In particular, a formalism is developed for calculating the matrix element of the operator $\mathrm{exp}(+sH)$, where $H$ is the quantum-mechanical Hamiltonian operator. It is found that the matrix element may be written as an ordinary $c$-number integral of a $c$-number function over the phase space. This function is found by doing a linked-cluster summation of the usual perturbation series. It is in the form of an exponential of an effective Hamiltonian, ${H}_{\mathrm{eff}}$, which is a function of the parameter $s$ and the $c$-number momentum and position coordinates. We have assumed that the Hamiltonian has no spin-dependent terms. The matrix element is written in such a form that it may be used for calculating the transition matrix elements in time-dependent perturbation theory as well as for calculating the canonical partition function ${Z}_{N}$ in quantum statistical mechanics. The former is obtained by substituting $s=\ensuremath{-}\frac{\mathrm{it}}{\ensuremath{\hbar}}$ and the latter by substituting $s=\ensuremath{-}\ensuremath{\beta}$. The series for ${H}_{\mathrm{eff}}$ depends only on the potential and is independent of the statistics. The effects of the statistics (Bose-Einstein, Fermi-Dirac, or Maxwell-Boltzmann) are all included in a certain factor multiplying $\mathrm{exp}(s{H}_{\mathrm{eff}})$. Thus we have succeeded in separating the statistical and quantum corrections. These latter arise from the noncommutativity of the position and momentum operators. The effective Hamiltonian is an infinite series. A diagrammatic technique is given for evaluating this series. It is an expansion in powers of the coupling constant and can also be expanded in powers of Planck's constant. In general, it is not possible to give a closed-form expression for ${H}_{\mathrm{eff}}$. But for most practical cases it is sufficient to keep the first few powers of Planck's constant $\ensuremath{\hbar}$, because it is very small. We therefore have an entirely new method of approximation. The statistical effect may always be included exactly. No numerical estimates are made in this paper.

Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Feb 15, 1973•0 citations

On the approximation of Feynman-Kac path integrals for quantum statistical mechanics

Stephen Bond, Brian B. Laird, Benedict Leimkuhler

Discretizations of the Feynman-Kac path integral representation of the quantum mechanical density matrix are investigated. Each infinite-dimensional path integral is approximated by a Riemann integral over a finite-dimensional function space, by restricting the integration to a subspace of all admissible paths. Using this process, a wide class of methods can be derived, with each method corresponding to a different choice for the approximating subspace. The traditional ``short-time'' approximation and ``Fourier discretization'' can be recovered from this approach, using linear and spectral basis functions respectively. As an illustration, a novel method is formulated using cubic elements and is shown to have improved convergence properties when applied to a simple model problem.

arXiv:cond-mat/0007112•Jul 6, 2000•0 citations
In this paper, four distinct ideas are combined, which under a wide range of circumstances can give very rapidly converging series expansions for functional integrals. (1) Expansion of the functional being integrated in functional Taylor series. In the familiar case arising in quantum statistical mechanics, that of the Wiener integral of exp[−∫V(x(t))dt],V(x) being the perturbing potential, this is equivalent to expanding the characteristic functional of the probability functional of V(x(t)) in central moments of V(x(t)). The lowest-order term of the series is the approximation obtained by Feynman and Hibbs through a variational method. (2) Transfer of the harmonic term of the potential, when the functional integral is the quantum-statistical density matrix (Green's function of the Bloch equation), to the weighting function. This transforms the functional integral from a Wiener to an Uhlenbeck-Ornstein integral. The formal expressions for the terms of the expansion are somewhat more complicated, but they can be worked out, and the result is a great improvement in the speed of convergence of the series with decreasing temperature and/or decreasing relative magnitude of the anharmonic part of the potential. (3) ``Reservation of variables'' in the integration. This amounts to breaking the averaging process down into an average over subsets of the distributed random function (conditional average), followed by an averaging of these averages. Any step of this kind (it may be repeated within the subsets, etc.) gives an improvement of accuracy. (4) When the quantity being evaluated through functional integration is the partition function, the device introduced by Feynman and Hibbs, of interchanging the functional integration with the integration of the Green's function over the equated initial and final configuration-space points, may be combined with the above techniques. This eliminates one integration in the terms of the expansion and seems to improve accuracy at the same time. The general series obtained is correlated with more conventional operator techniques of quantum-mechanical perturbation theory, in order to answer the perennial question, does the path-integral method bring with it anything that could not be derived by other methods? It is in some sense a Feynman-Dyson expansion of the Green's function, but one that is further modified mathematically in a way characteristic only of the path-integral point of view, and which, moreover, improves its accuracy. It thus appears unlikely that the result is merely one of standard type disguised as a functional integral result. Sample numerical calculations are given to assay the accuracy of the methods, which are shown to compare very favorably with the traditional approximation of finite subdivision of the time interval.
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.
Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations
The method of functional differentiation, used in classical statistical mechanics to obtain approximate integral equations for the pair distribution function, is extended to quantum systems obeying Maxwell-Boltzmann statistics. The grand partition function is written as a path integral, and space-temperature distributions are generated from it by successive functional differentiation. Distributions can be expanded in powers of an external potential via a functional Taylor series. When the series is truncated so that no distributions higher than second order enter and the external potential is specialized to one arising from a fixed particle in the system, coupled integral equations result for the two kinds of pair distributions. The first-order Ursell function calculated from these equations yields the Montroll-Ward ring summation for the grand potential. This approximate Ursell function determines also the Fourier transform of the momentum density (useful in calculating the p = 0 occupation number). Although divergences appear in the low-temperature limit, they can be removed by a simple modification of the independent functional. An extension to Fermi-Dirac and Böse-Einstein statistics is indicated at appropriate places in the text.
Let us recapitulate. Although the form (4.32) for the path integral was motivated by analogy with finite-dimensional integrals, in no way is it a definition. Rather, as we have stressed repeatedly, it is an algorithm for handling the Feynman perturbation series. Strong coupling apart, all the manipulations of the path integral that we have performed (integration by parts, the introduction of auxiliary fields, the inclusion of background fields) stand in one-to-one correspondence to rearrangements of the perturbation series. In fact, path integrals are unnecessary in the derivation of expressions like (4.50). See Fried (1972), for example. The path integrals are as rigorous – or non-rigorous – as the series, irrespective of whether they can be given proper mathematical definition or not. At the same time, the analogy with finite-dimensional integrals suggests approaches that go beyond this purely operational stance. In particular, because of the dominant role played by the classical action in the path integral formalism we would like to use it to establish the relation of quantum dynamics to classical dynamics. This is the most compelling reason for introducing the path integral formalism, requiring a distinct relaxation in our treatment of it. By formal analogy with finite-dimensional saddle-point and stationary-phase methods, elevating the classical action to the phase of the integrand promotes series expansions in ħ as the natural way to proceed.
This is the first in a series of papers. It deals with the linked-cluster summation of the perturbation series in quantum mechanics and quantum statistical physics. In particular, a formalism is developed for calculating the matrix element of the operator exp(+sH)\mathrm{exp}(+sH)exp(+sH), where HHH is the quantum-mechanical Hamiltonian operator. It is found that the matrix element may be written as an ordinary ccc-number integral of a ccc-number function over the phase space. This function is found by doing a linked-cluster summation of the usual perturbation series. It is in the form of an exponential of an effective Hamiltonian, Heff{H}_{\mathrm{eff}}Heff​, which is a function of the parameter sss and the ccc-number momentum and position coordinates. We have assumed that the Hamiltonian has no spin-dependent terms. The matrix element is written in such a form that it may be used for calculating the transition matrix elements in time-dependent perturbation theory as well as for calculating the canonical partition function ZN{Z}_{N}ZN​ in quantum statistical mechanics. The former is obtained by substituting $s=\ensuremath{-}\frac{\mathrm{it}}{\ensuremath{\hbar}}$ and the latter by substituting $s=\ensuremath{-}\ensuremath{\beta}$. The series for Heff{H}_{\mathrm{eff}}Heff​ depends only on the potential and is independent of the statistics. The effects of the statistics (Bose-Einstein, Fermi-Dirac, or Maxwell-Boltzmann) are all included in a certain factor multiplying exp(sHeff)\mathrm{exp}(s{H}_{\mathrm{eff}})exp(sHeff​). Thus we have succeeded in separating the statistical and quantum corrections. These latter arise from the noncommutativity of the position and momentum operators. The effective Hamiltonian is an infinite series. A diagrammatic technique is given for evaluating this series. It is an expansion in powers of the coupling constant and can also be expanded in powers of Planck's constant. In general, it is not possible to give a closed-form expression for Heff{H}_{\mathrm{eff}}Heff​. But for most practical cases it is sufficient to keep the first few powers of Planck's constant $\ensuremath{\hbar}$, because it is very small. We therefore have an entirely new method of approximation. The statistical effect may always be included exactly. No numerical estimates are made in this paper.
Discretizations of the Feynman-Kac path integral representation of the quantum mechanical density matrix are investigated. Each infinite-dimensional path integral is approximated by a Riemann integral over a finite-dimensional function space, by restricting the integration to a subspace of all admissible paths. Using this process, a wide class of methods can be derived, with each method corresponding to a different choice for the approximating subspace. The traditional ``short-time'' approximation and ``Fourier discretization'' can be recovered from this approach, using linear and spectral basis functions respectively. As an illustration, a novel method is formulated using cubic elements and is shown to have improved convergence properties when applied to a simple model problem.