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

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

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

Published Feb 17, 2022• 0 citations
Publisher

Abstract

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.

DOI
10.21203/rs.3.rs-1335580/v1

Related papers

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

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

Randall W. Hall

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.

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

Lagrange Approach in Stochastic Quantization of the 2D Hubbard Model

А. Н. Калинкин, Vitaliy M. Skorikov

A new method of the generation of Green's function in external local and bilocal fields is proposed for the stochastic quantization of the 2D Hubbard model. The differential equations for Green's functions transform to the system of integral equations and it is solved numerically for a local Hubbard-Stratanovich field. The Trotter-Suzuki expansion and the matrix elements of the evolution operator are not used in this approach. For a 4 × 4 lattice with half-filled band the total energy E, the density n, the density–density 〈 ni↑ni↓〉 and spin–spin 〈SziSzj〉 correlation functions were calculated by stochastic quantization as the functions of one-site repusion U and inverse temperature β. The obtained results are in agreement with the previous Monte Carlo calculations.

physica status solidi (b)•May 1, 2000•2 citations

A unified theory of variational and quantum Monte Carlo methods and beyond

Mohammad-Sadegh Vaezi, Abolhassan Vaezi

We present a unified theory of the variational Monte Carlo (VMC) and determinant quantum Monte Carlo (DQMC) methods using a novel density matrix formulation of VMC. We introduce an efficient algorithm for VMC to compute correlation functions and expectation values based on the auxiliary field Hirsch-Hubbard-Stratonovic transformation. We show that this new approach to VMC converges significantly faster than its traditional implementations. Furthermore, we generalize the Trotter-Suzuki decomposition to finite imaginary time steps $τ\sim O(1)$ and develop a variational quantum Monte Carlo (VQMC) method accordingly, which is more accurate than VMC and can incorporate quantum fluctuations more efficiently. The two extreme limits of the VQMC method, namely infinitesimal and infinite imaginary time steps, correspond to the DQMC and VMC techniques, respectively. We demonstrate that our VQMC allows us to access lower temperatures in comparison with the conventional DQMC before the sign problem comes into play. We finally show that our VQMC can also enhance the accuracy of the projector Monte Carlo methods by providing better and less biased candidates for the trial wave functions, requiring shorter projection times for a given accuracy and alleviating the sign problem further.

arXiv:1810.00864•Oct 1, 2018•2 citations
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.
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.
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 new method of the generation of Green's function in external local and bilocal fields is proposed for the stochastic quantization of the 2D Hubbard model. The differential equations for Green's functions transform to the system of integral equations and it is solved numerically for a local Hubbard-Stratanovich field. The Trotter-Suzuki expansion and the matrix elements of the evolution operator are not used in this approach. For a 4 × 4 lattice with half-filled band the total energy E, the density n, the density–density 〈 ni↑ni↓〉 and spin–spin 〈SziSzj〉 correlation functions were calculated by stochastic quantization as the functions of one-site repusion U and inverse temperature β. The obtained results are in agreement with the previous Monte Carlo calculations.
We present a unified theory of the variational Monte Carlo (VMC) and determinant quantum Monte Carlo (DQMC) methods using a novel density matrix formulation of VMC. We introduce an efficient algorithm for VMC to compute correlation functions and expectation values based on the auxiliary field Hirsch-Hubbard-Stratonovic transformation. We show that this new approach to VMC converges significantly faster than its traditional implementations. Furthermore, we generalize the Trotter-Suzuki decomposition to finite imaginary time steps τ∼O(1)τ\sim O(1)τ∼O(1) and develop a variational quantum Monte Carlo (VQMC) method accordingly, which is more accurate than VMC and can incorporate quantum fluctuations more efficiently. The two extreme limits of the VQMC method, namely infinitesimal and infinite imaginary time steps, correspond to the DQMC and VMC techniques, respectively. We demonstrate that our VQMC allows us to access lower temperatures in comparison with the conventional DQMC before the sign problem comes into play. We finally show that our VQMC can also enhance the accuracy of the projector Monte Carlo methods by providing better and less biased candidates for the trial wave functions, requiring shorter projection times for a given accuracy and alleviating the sign problem further.