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 Physics A Mathematical and General

A high-temperature approximation for the path-integral quantum Monte Carlo method

Miroslav Kolář, S. F. O’Shea

Published Jul 7, 1996• 9 citations
Publisher

Abstract

A high-temperature approximation for the discretized path-integral quantum Monte Carlo (PIQMC) method is formulated. At higher temperatures, all P fictitious classical particles representing a single quantum particle stay close together, and an efficient approximation is obtained when, in the primitive short-time propagator, an essentially local harmonic approximation is used for the external potential at the common centre of mass of P fictitious particles - the integration over P - 1 dimensions can then be carried out analytically, and a classical formula of the effective-potential type is obtained for the partition function. Also discussed are the proper form and applicability of the virial total-energy estimator for finite systems, and the computation of the temperature dependence of the configurational partition function in a single PIQMC run.

DOI
10.1088/0305-4470/29/13/018

Related papers

Path-integral Monte Carlo

Richard M. Martin, Lucia Reining, David M. Ceperley

It is shown from first principles that, in spite of the large interatomic forces, liquid 4 He should exhibit a transition analogous to the transition in an ideal Bose–Einstein gas. The exact partition function is written as an integral over trajectories, using the space-time approach to quantum mechanics. R.P. Feynman, 1953 Summary In this chapter we discuss imaginary-time path integrals and the path-integral Monte Carlo method for the calculation of properties of quantum systems at non-zero temperature. We discuss how Fermi and Bose statistics enter, and how to generalize the fixed-node procedure to non-zero temperature. We then discuss an auxiliary-field method for the Hubbard model. The path-integral method can be used to perform ground-state calculations, allowing calculations of properties with less bias than the projector Monte Carlo method. We also discuss the problem of estimating real-time response functions using information from imaginary-time correlation functions. In previous chapters we described two QMC methods, namely variational QMC and projector (diffusion) QMC. Both of these methods are zero temperature, or, more properly, are formulated for single states. In this chapter we discuss the path-integral Monte Carlo (PIMC) method, which is explicitly formulated at non-zero temperature. Directly including temperature is important because many, if not most, measurements and practical applications involve significant thermal effects. One might think that to do calculations at a non-zero temperature we would have to explicitly sum over excited states. Such a summation would be difficult to accomplish once the temperature is above the energy gap, because there are so many possible many-body excitations. In addition, the properties for each excitation are more difficult to calculate than those for the ground state. As we will see, path-integral methods do not require an explicit sum over excitations. As an added bonus, they provide an interesting and enlightening window through which to view quantum systems. However, the sign problem, introduced in the previous chapter, is still present for fermion systems. The fixed-node approximation is used again. An advantage of PIMC is the absence of a trial wavefunction. As a result, quantum expectation values, including ones not involving the energy, can be computed directly. For the expert, the lack of an importance function may seem a disadvantage; without it one cannot push the simulation in a preferred direction.

Cambridge University Press eBooks•May 31, 2016•1 citation

The harmonic oscillator at finite temperature using path integrals

Åsmund Larsen, F. Ravndal

The statistical mechanics of a quantum harmonic oscillator can easily be formulated in terms of Feynman path integrals in imaginary time. By straightforward integrations, the partition function, the internal energy, and correlation functions are found. When time is discretized, the same results can be obtained directly using Monte Carlo methods on a computer. In the zero temperature limit, the usual quantum mechanical results are recovered.

American Journal of Physics•Dec 1, 1988•4 citations

Path-integral Monte Carlo method in quantum statistics for a system ofNidentical fermions

Alexander P. Lyubartsev, P. N. Vorontsov‐Velyaminov

A rigorous expression for the partition function of N identical fermions with spin 1/2 is obtained in a form suitable for a path-integral Monte Carlo (PIMC) simulation of a many-electron system at a finite temperature. A sum over permutations is reduced to a sum over classes with readily determinable factors. A distribution over projections and values of the total spin is also obtained, enabling PIMC calculations of spin-dependent quantities. An MC algorithm is proposed which accounts simultaneously for all types of trajectory linkages. A test calculation, four electrons in a spherical cavity with a radius of 1 nm at a temperature of 300--2000 K, is performed.

Physical Review A•Dec 1, 1993•38 citations

Path integral Monte Carlo on a lattice: Extended states

Mark O'Callaghan, Bruce N. Miller

The equilibrium properties of a single quantum particle (qp) interacting with a classical gas for a wide range of temperatures that explore the system's behavior in the classical as well as in the quantum regime is investigated. Both the qp and atoms are restricted to the sites of a one-dimensional lattice. A path integral formalism is developed within the context of the canonical ensemble in which the qp is represented by a closed, variable-step random walk on the lattice. Monte Carlo methods are employed to determine the system's properties. For the case of a free particle, analytical expressions for the energy, its fluctuations, and the qp-qp correlation function are derived and compared with the Monte Carlo simulations. To test the usefulness of the path integral formalism, the Metropolis algorithm is employed to determine the equilibrium properties of the qp for a periodic interaction potential, forcing the qp to occupy extended states. We consider a striped potential in one dimension, where every other lattice site is occupied by an atom with potential ε, and every other lattice site is empty. This potential serves as a stress test for the path integral formalism because of its rapid site-to-site variation. An analytical solution was determined in this case by utilizing Bloch's theorem due to the periodicity of the potential. Comparisons of the potential energy, the total energy, the energy fluctuations, and the correlation function are made between the results of the Monte Carlo simulations and the analytical calculations.

arXiv:1402.1820•Physical Review E•Apr 15, 2014•6 citations

Rapid mixing of path integral Monte Carlo for 1D stoquastic Hamiltonians

Elizabeth Crosson, Aram W. Harrow

Path integral quantum Monte Carlo (PIMC) is a method for estimating thermal equilibrium properties of stoquastic quantum spin systems by sampling from a classical Gibbs distribution using Markov chain Monte Carlo. The PIMC method has been widely used to study the physics of materials and for simulated quantum annealing, but these successful applications are rarely accompanied by formal proofs that the Markov chains underlying PIMC rapidly converge to the desired equilibrium distribution. In this work we analyze the mixing time of PIMC for 1D stoquastic Hamiltonians, including disordered transverse Ising models (TIM) with long-range algebraically decaying interactions as well as disordered XY spin chains with nearest-neighbor interactions. By bounding the convergence time to the equilibrium distribution we rigorously justify the use of PIMC to approximate partition functions and expectations of observables for these models at inverse temperatures that scale at most logarithmically with the number of qubits. The mixing time analysis is based on the canonical paths method applied to the single-site Metropolis Markov chain for the Gibbs distribution of 2D classical spin models with couplings related to the interactions in the quantum Hamiltonian. Since the system has strongly nonisotropic couplings that grow with system size, it does not fall into the known cases where 2D classical spin models are known to mix rapidly.

arXiv:1812.02144•Quantum•Feb 11, 2021•10 citations

Rapid mixing of path integral Monte Carlo for 1D stoquastic Hamiltonians

Elizabeth Crosson, Aram W. Harrow

Path integral quantum Monte Carlo (PIMC) is a method for estimating thermal equilibrium properties of stoquastic quantum spin systems by sampling from a classical Gibbs distribution using Markov chain Monte Carlo. The PIMC method has been widely used to study the physics of materials and for simulated quantum annealing, but these successful applications are rarely accompanied by formal proofs that the Markov chains underlying PIMC rapidly converge to the desired equilibrium distribution. In this work we analyze the mixing time of PIMC for 1D stoquastic Hamiltonians, including disordered transverse Ising models (TIM) with long-range algebraically decaying interactions as well as disordered XY spin chains with nearest-neighbor interactions. By bounding the convergence time to the equilibrium distribution we rigorously justify the use of PIMC to approximate partition functions and expectations of observables for these models at inverse temperatures that scale at most logarithmically with the number of qubits. The mixing time analysis is based on the canonical paths method applied to the single-site Metropolis Markov chain for the Gibbs distribution of 2D classical spin models with couplings related to the interactions in the quantum Hamiltonian. Since the system has strongly nonisotropic couplings that grow with system size, it does not fall into the known cases where 2D classical spin models are known to mix rapidly.

•Jan 1, 2021•7 citations
A high-temperature approximation for the discretized path-integral quantum Monte Carlo (PIQMC) method is formulated. At higher temperatures, all P fictitious classical particles representing a single quantum particle stay close together, and an efficient approximation is obtained when, in the primitive short-time propagator, an essentially local harmonic approximation is used for the external potential at the common centre of mass of P fictitious particles - the integration over P - 1 dimensions can then be carried out analytically, and a classical formula of the effective-potential type is obtained for the partition function. Also discussed are the proper form and applicability of the virial total-energy estimator for finite systems, and the computation of the temperature dependence of the configurational partition function in a single PIQMC run.
It is shown from first principles that, in spite of the large interatomic forces, liquid 4 He should exhibit a transition analogous to the transition in an ideal Bose–Einstein gas. The exact partition function is written as an integral over trajectories, using the space-time approach to quantum mechanics. R.P. Feynman, 1953 Summary In this chapter we discuss imaginary-time path integrals and the path-integral Monte Carlo method for the calculation of properties of quantum systems at non-zero temperature. We discuss how Fermi and Bose statistics enter, and how to generalize the fixed-node procedure to non-zero temperature. We then discuss an auxiliary-field method for the Hubbard model. The path-integral method can be used to perform ground-state calculations, allowing calculations of properties with less bias than the projector Monte Carlo method. We also discuss the problem of estimating real-time response functions using information from imaginary-time correlation functions. In previous chapters we described two QMC methods, namely variational QMC and projector (diffusion) QMC. Both of these methods are zero temperature, or, more properly, are formulated for single states. In this chapter we discuss the path-integral Monte Carlo (PIMC) method, which is explicitly formulated at non-zero temperature. Directly including temperature is important because many, if not most, measurements and practical applications involve significant thermal effects. One might think that to do calculations at a non-zero temperature we would have to explicitly sum over excited states. Such a summation would be difficult to accomplish once the temperature is above the energy gap, because there are so many possible many-body excitations. In addition, the properties for each excitation are more difficult to calculate than those for the ground state. As we will see, path-integral methods do not require an explicit sum over excitations. As an added bonus, they provide an interesting and enlightening window through which to view quantum systems. However, the sign problem, introduced in the previous chapter, is still present for fermion systems. The fixed-node approximation is used again. An advantage of PIMC is the absence of a trial wavefunction. As a result, quantum expectation values, including ones not involving the energy, can be computed directly. For the expert, the lack of an importance function may seem a disadvantage; without it one cannot push the simulation in a preferred direction.
A rigorous expression for the partition function of N identical fermions with spin 1/2 is obtained in a form suitable for a path-integral Monte Carlo (PIMC) simulation of a many-electron system at a finite temperature. A sum over permutations is reduced to a sum over classes with readily determinable factors. A distribution over projections and values of the total spin is also obtained, enabling PIMC calculations of spin-dependent quantities. An MC algorithm is proposed which accounts simultaneously for all types of trajectory linkages. A test calculation, four electrons in a spherical cavity with a radius of 1 nm at a temperature of 300--2000 K, is performed.
The equilibrium properties of a single quantum particle (qp) interacting with a classical gas for a wide range of temperatures that explore the system's behavior in the classical as well as in the quantum regime is investigated. Both the qp and atoms are restricted to the sites of a one-dimensional lattice. A path integral formalism is developed within the context of the canonical ensemble in which the qp is represented by a closed, variable-step random walk on the lattice. Monte Carlo methods are employed to determine the system's properties. For the case of a free particle, analytical expressions for the energy, its fluctuations, and the qp-qp correlation function are derived and compared with the Monte Carlo simulations. To test the usefulness of the path integral formalism, the Metropolis algorithm is employed to determine the equilibrium properties of the qp for a periodic interaction potential, forcing the qp to occupy extended states. We consider a striped potential in one dimension, where every other lattice site is occupied by an atom with potential ε, and every other lattice site is empty. This potential serves as a stress test for the path integral formalism because of its rapid site-to-site variation. An analytical solution was determined in this case by utilizing Bloch's theorem due to the periodicity of the potential. Comparisons of the potential energy, the total energy, the energy fluctuations, and the correlation function are made between the results of the Monte Carlo simulations and the analytical calculations.
Path integral quantum Monte Carlo (PIMC) is a method for estimating thermal equilibrium properties of stoquastic quantum spin systems by sampling from a classical Gibbs distribution using Markov chain Monte Carlo. The PIMC method has been widely used to study the physics of materials and for simulated quantum annealing, but these successful applications are rarely accompanied by formal proofs that the Markov chains underlying PIMC rapidly converge to the desired equilibrium distribution. In this work we analyze the mixing time of PIMC for 1D stoquastic Hamiltonians, including disordered transverse Ising models (TIM) with long-range algebraically decaying interactions as well as disordered XY spin chains with nearest-neighbor interactions. By bounding the convergence time to the equilibrium distribution we rigorously justify the use of PIMC to approximate partition functions and expectations of observables for these models at inverse temperatures that scale at most logarithmically with the number of qubits. The mixing time analysis is based on the canonical paths method applied to the single-site Metropolis Markov chain for the Gibbs distribution of 2D classical spin models with couplings related to the interactions in the quantum Hamiltonian. Since the system has strongly nonisotropic couplings that grow with system size, it does not fall into the known cases where 2D classical spin models are known to mix rapidly.
Path integral quantum Monte Carlo (PIMC) is a method for estimating thermal equilibrium properties of stoquastic quantum spin systems by sampling from a classical Gibbs distribution using Markov chain Monte Carlo. The PIMC method has been widely used to study the physics of materials and for simulated quantum annealing, but these successful applications are rarely accompanied by formal proofs that the Markov chains underlying PIMC rapidly converge to the desired equilibrium distribution. In this work we analyze the mixing time of PIMC for 1D stoquastic Hamiltonians, including disordered transverse Ising models (TIM) with long-range algebraically decaying interactions as well as disordered XY spin chains with nearest-neighbor interactions. By bounding the convergence time to the equilibrium distribution we rigorously justify the use of PIMC to approximate partition functions and expectations of observables for these models at inverse temperatures that scale at most logarithmically with the number of qubits. The mixing time analysis is based on the canonical paths method applied to the single-site Metropolis Markov chain for the Gibbs distribution of 2D classical spin models with couplings related to the interactions in the quantum Hamiltonian. Since the system has strongly nonisotropic couplings that grow with system size, it does not fall into the known cases where 2D classical spin models are known to mix rapidly.