The cumulant representation of the Fourier path integral method is examined to determine the asymptotic convergence characteristics of the imaginary-time density matrix with respect to the number of path variables N included. It is proved that when the cumulant expansion is truncated at order p, the asymptotic convergence rate of the density matrix behaves like N(-(2p+1)). The complex algebra associated with the proof is simplified by introducing a diagrammatic representation of the contributing terms along with an associated linked-cluster theorem. The cumulant terms at each order are expanded in a series such that the asymptotic convergence rate is maintained without the need to calculate the full cumulant at order p. Using this truncated expansion of each cumulant at order p, the numerical cost in developing Fourier path integral expressions having convergence order N(-(2p+1)) is shown to be approximately linear in the number of required potential energy evaluations making the method promising for actual numerical implementation.
The cumulant representation of the Fourier path integral method is examined to determine the asymptotic convergence characteristics of the imaginary-time density matrix with respect to the number of path variables N included. It is proved that when the cumulant expansion is truncated at order p, the asymptotic convergence rate of the density matrix behaves like N(-(2p+1)). The complex algebra associated with the proof is simplified by introducing a diagrammatic representation of the contributing terms along with an associated linked-cluster theorem. The cumulant terms at each order are expanded in a series such that the asymptotic convergence rate is maintained without the need to calculate the full cumulant at order p. Using this truncated expansion of each cumulant at order p, the numerical cost in developing Fourier path integral expressions having convergence order N(-(2p+1)) is shown to be approximately linear in the number of required potential energy evaluations making the method promising for actual numerical implementation.
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.
We develop an approach of calculating the many-body path integral based on the linked cluster expansion method. First, we derive a linked cluster expansion and we give the diagrammatic rules for calculating the free energy and the pair distribution function g(r) as a systematic power-series expansion in the particle density. We also present a structured Pad\'e approximation scheme in the momentum space to determine g(r). The calculated g(r) for distinguishable particles interacting with Lennard-Jones and hard-sphere potential in various attempted schemes of approximation of the diagrammatic series compares very well with the results of path-integral Monte Carlo simulation. Our method is applicable to a wide range of problems of current general interest and may be extended to the case of identical particles and, in particular, to the case of the many-fermion problem.
Two path integral representations for the T matrix in nonrelativistic potential scattering are derived and proved to produce the complete Born series when expanded to all orders. They are obtained with the help of ``phantom'' degrees of freedom which take away explicit phases that diverge for asymptotic times. In addition, energy conservation is enforced by imposing a Faddeev-Popov-like constraint in the velocity path integral. These expressions may be useful for attempts to evaluate the path integral in real time and for alternative multiple scattering expansions. Standard eikonal-type high-energy approximations and systematic expansions immediately follow.
The direct transformation from the polygonal to the Fourier series representation of Feynman path integrals, via a change of integration variables, is effected explicitly for cases where the Lagrangian is of the form L(x,ẋ,t)= 1/2 ẋ2−V(x,t). This transformation involves a ‘‘functional Jacobian’’ stemming solely from the velocity term in the Lagrangian; this is because N-segment polygonal paths of the type contributing significantly to the integral, and their N-term Fourier series approximants, coalesce together as N→∞, but not their derivatives. We also consider integrals over paths with fixed means T−1∫T0 dt x(t) =const. The usefulness of the Fourier representation is illustrated with the harmonic oscillator case V(x)= 1/2 ω2x2, in both the free and fixed means situations; in particular, the Fourier evaluation of the path integrals trivially determines the large time phases (Maslov indices), and the ranges of ω2 values for which the integrals are finite or infinite in the imaginary time cases.
A useful approximation for the thermal operator exp(-betaH[over ]) is based on its representation in terms of either frozen or thawed Gaussian states. Such approximate representations are leading-order terms in respective series representations of the thermal operator. A numerical study of the convergence properties of the frozen Gaussian series representation has been recently published. In this paper, we extend the previous study to include also the convergence properties of the more expensive thawed Gaussian series representation of the thermal operator. We consider three different formulations for the series representation and apply them to a quartic double-well potential to find that the thawed Gaussian series representation converges faster than the frozen Gaussian one. Further analysis is presented as to the convergence properties and the numerical efficiency of three different thawed Gaussian series representation. The unsymmetrized form converges most rapidly, however, the lower order approximations of the symmetrized forms are more accurate. Comparison with a standard discretized path-integral evaluation demonstrates that the Gaussian based perturbation series representation converges much faster.
The reweighted random series techniques provide finite-dimensional approximations to the quantum density matrix of a physical system that have fast asymptotic convergence. We study two special reweighted techniques that are based upon the Lévy–Ciesielski and Wiener–Fourier series, respectively. In agreement with the theoretical predictions, we demonstrate by numerical examples that the asymptotic convergence of the two reweighted methods is cubic for smooth enough potentials. For each reweighted technique, we propose some minimalist quadrature techniques for the computation of the path averages. These quadrature techniques are designed to preserve the asymptotic convergence of the original methods.