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

The Journal of Chemical Physics

Path integral studies of the 2D Hubbard model using a new projection operator

Randall W. Hall

Published Jan 15, 1991• 5 citations
Publisher

Abstract

Feynman’s path integral formulation of quantum mechanics, supplemented by an approximate projection operator (exact in the case of noninteracting particles), is used to study the 2D Hubbard model. The projection operator is designed to study Hamiltonians defined on a finite basis set, but extensions to continuous basis sets are suggested. The projection operator is shown to reduce the variance by a significant amount relative to straightforward Monte Carlo integration. Approximate calculations are usually within one standard deviation of exact results and virtually always within two to three standard deviations. In addition, the algorithm scales with the number of discretization points P as either P or P2 (depending on the method of implementation), rather than the P3 of the Hubbard–Stratonovich transformation. Accuracy to about 5%–10% in energies and spin–spin correlation functions are found using moderate amounts of computer time.

DOI
10.1063/1.459988

Related papers

Formally exact path integral Monte Carlo calculations using approximate projection operators

Randall W. Hall

A class of approximate projection operators is used to reduce the variance in path integral Monte Carlo calculations in a formally exact manner. Paths are classified according to the projection operators, allowing the identification of paths whose contribution to the variance is negligible. The approach is applied to two canonical systems: Two like-spin electrons in a three dimensional harmonic well and the two dimensional Hubbard model, representing localized and extended electronic states, respectively. Time savings of 15–900 over straightforward Monte Carlo calculations are observed.

The Journal of Chemical Physics•Nov 1, 1992•7 citations

Efficient World-Line-Based Quantum Monte Carlo Method Without Hubbard-Stratonovich Transformation

Jianpan Wang, W. Pan, D. Y. Sun

Abstract By precisely writing down the matrix element of the local Boltzmann operator (e-th, where h is the Hermitian conjugate pairs of off-diagonal operators), we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. With current formula, the Hubbard–Stratonovich transformation is not necessary, and is not based on the determinant approach, which can improve the computational efficiency. The results show that, the simulation time has the square-law scaling with system sizes, which is comparable with the usual first-principle calculation. The current formula also improves the accuracy of the Suzuki–Trotter decomposition. As an example, we have studied the one-dimensional half-filled Hubbard model at finite temperature. The obtained results are in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm can be used in various studies.

•Feb 17, 2022•0 citations

Development, justification, and use of a projection operator in path integral calculations in continuous space

Randall W. Hall, Melissa R. Prince

A projection operator, similar to one previously used by us for problems with a finite set of basis functions, is suggested for use with continuous basis sets. This projection operator requires knowledge of the nodes of the density matrix at all temperatures. We show that a class of nodes, determined from the noninteracting density matrix and present at high temperatures in the interacting system are preserved to first order in the interaction at low temperatures. While we cannot show that the nodes are present at intermediate temperatures, we suspect they do exist and, as a test of this conjecture, we perform a calculation of two electrons confined in a harmonic well, using the projection operator. We find that accurate results are obtained at a range of temperatures, suggesting that our conjecture is indeed correct. We find that the error limits determined using the projection operator are 1–2 times smaller than those obtained with straightforward Monte Carlo integration (corresponding to a reduction in time of 1–4 in obtaining a desired level of accuracy).

The Journal of Chemical Physics•Oct 15, 1991•5 citations

Efficient world-line-based quantum Monte Carlo method without Hubbard–Stratonovich transformation

J Wang, W. Pan, D. Y. Sun

By precisely writing down the matrix element of the local Boltzmann operator ([Formula: see text], where [Formula: see text] is the Hermitian conjugate pairs of off-diagonal operators), we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. With the current formula, the Hubbard-Stratonovich transformation is not necessary, accordingly the determinant calculation is not needed, which can improve the computational efficiency. The results show that, the simulation time has the square-law scaling with system sizes, which is comparable with the usual first-principles calculations. The current formula also improves the accuracy of the Suzuki-Trotter decomposition. As an example, we have studied the one-dimensional half-filled Hubbard model at finite temperature. The obtained results are in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm could be applied to various studies in future.

Scientific Reports•May 17, 2022•0 citations

Efficient World-line-based Quantum Monte Carlo Method Without Hubbard-Stratonovich Transformation

J. Wang, W. Pan, D. Y. Sun

By precisely writing down the matrix element of the local Boltzmann operator, we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. With current formula, the Hubbard-Stratonovich transformation is not necessary, and is not based on the determinant approach, which can improve the computational efficiency. The results show that, the simulation time has the square-law scaling with system sizes, which is comparable with the usual first-principle calculation. The current formula also improves the accuracy of the Suzuki-Trotter decomposition. As an example, we have studied the one-dimensional half-filled Hubbard model at finite temperature. The obtained results are in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm can be used in various studies.

arXiv:2110.08742•Oct 17, 2021•0 citations

Quantum Monte Carlo Method without Hubbard-Stratonovich Transformation for Many-body Systems

J. Wang, Weiwei Pan, D. Y. Sun

By precisely writing down the matrix element of the local Boltzmann operator ($e^(-{\tau}h$), where h is the paired operator formed by two interacting terms conjugate to each other), we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. Our formula presents remarkable outcomes such that: 1) the Hubbard-Stratonovich transformation is not necessary, and the auxiliary field is no longer needed, which can improve the computational efficiency. 2) The new formula improves the accuracy of the Suzuki-Trotter decomposition. Our method also provides a potential theoretical framework for analyzing and solving the fermion sign problem. As an example, we studied the one-dimensional half-filled Hubbard model at a finite temperature. The obtained results were in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm could be used in various studies.

•Oct 17, 2021•0 citations
Feynman’s path integral formulation of quantum mechanics, supplemented by an approximate projection operator (exact in the case of noninteracting particles), is used to study the 2D Hubbard model. The projection operator is designed to study Hamiltonians defined on a finite basis set, but extensions to continuous basis sets are suggested. The projection operator is shown to reduce the variance by a significant amount relative to straightforward Monte Carlo integration. Approximate calculations are usually within one standard deviation of exact results and virtually always within two to three standard deviations. In addition, the algorithm scales with the number of discretization points P as either P or P2 (depending on the method of implementation), rather than the P3 of the Hubbard–Stratonovich transformation. Accuracy to about 5%–10% in energies and spin–spin correlation functions are found using moderate amounts of computer time.
A class of approximate projection operators is used to reduce the variance in path integral Monte Carlo calculations in a formally exact manner. Paths are classified according to the projection operators, allowing the identification of paths whose contribution to the variance is negligible. The approach is applied to two canonical systems: Two like-spin electrons in a three dimensional harmonic well and the two dimensional Hubbard model, representing localized and extended electronic states, respectively. Time savings of 15–900 over straightforward Monte Carlo calculations are observed.
Abstract By precisely writing down the matrix element of the local Boltzmann operator (e-th, where h is the Hermitian conjugate pairs of off-diagonal operators), we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. With current formula, the Hubbard–Stratonovich transformation is not necessary, and is not based on the determinant approach, which can improve the computational efficiency. The results show that, the simulation time has the square-law scaling with system sizes, which is comparable with the usual first-principle calculation. The current formula also improves the accuracy of the Suzuki–Trotter decomposition. As an example, we have studied the one-dimensional half-filled Hubbard model at finite temperature. The obtained results are in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm can be used in various studies.
A projection operator, similar to one previously used by us for problems with a finite set of basis functions, is suggested for use with continuous basis sets. This projection operator requires knowledge of the nodes of the density matrix at all temperatures. We show that a class of nodes, determined from the noninteracting density matrix and present at high temperatures in the interacting system are preserved to first order in the interaction at low temperatures. While we cannot show that the nodes are present at intermediate temperatures, we suspect they do exist and, as a test of this conjecture, we perform a calculation of two electrons confined in a harmonic well, using the projection operator. We find that accurate results are obtained at a range of temperatures, suggesting that our conjecture is indeed correct. We find that the error limits determined using the projection operator are 1–2 times smaller than those obtained with straightforward Monte Carlo integration (corresponding to a reduction in time of 1–4 in obtaining a desired level of accuracy).
By precisely writing down the matrix element of the local Boltzmann operator ([Formula: see text], where [Formula: see text] is the Hermitian conjugate pairs of off-diagonal operators), we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. With the current formula, the Hubbard-Stratonovich transformation is not necessary, accordingly the determinant calculation is not needed, which can improve the computational efficiency. The results show that, the simulation time has the square-law scaling with system sizes, which is comparable with the usual first-principles calculations. The current formula also improves the accuracy of the Suzuki-Trotter decomposition. As an example, we have studied the one-dimensional half-filled Hubbard model at finite temperature. The obtained results are in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm could be applied to various studies in future.
By precisely writing down the matrix element of the local Boltzmann operator, we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. With current formula, the Hubbard-Stratonovich transformation is not necessary, and is not based on the determinant approach, which can improve the computational efficiency. The results show that, the simulation time has the square-law scaling with system sizes, which is comparable with the usual first-principle calculation. The current formula also improves the accuracy of the Suzuki-Trotter decomposition. As an example, we have studied the one-dimensional half-filled Hubbard model at finite temperature. The obtained results are in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm can be used in various studies.
By precisely writing down the matrix element of the local Boltzmann operator (e(−τhe^(-{\tau}he(−τh), where h is the paired operator formed by two interacting terms conjugate to each other), we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. Our formula presents remarkable outcomes such that: 1) the Hubbard-Stratonovich transformation is not necessary, and the auxiliary field is no longer needed, which can improve the computational efficiency. 2) The new formula improves the accuracy of the Suzuki-Trotter decomposition. Our method also provides a potential theoretical framework for analyzing and solving the fermion sign problem. As an example, we studied the one-dimensional half-filled Hubbard model at a finite temperature. The obtained results were in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm could be used in various studies.