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

arXiv:0704.3792 · Journal of the Physical Society of Japan

Gaussian-Basis Monte Carlo Method for Numerical Study on Ground States of Itinerant and Strongly Correlated Electron Systems

Takeshi Aimi, Masatoshi Imada

Published Aug 10, 2007• 33 citations
arXivHTMLPDF

Abstract

We examine Gaussian-basis Monte Carlo method (GBMC) introduced by Corney and Drummond. This method is based on an expansion of the density-matrix operator rho by means of the coherent Gaussian-type operator basis Lambda and does not suffer from the minus sign problem. The original method, however, often fails in reproducing the true ground state and causes systematic errors of calculated physical quantities because the samples are often trapped in some metastable or symmetry broken states. To overcome this difficulty, we combine the quantum-number projection scheme proposed by Assaad, Werner, Corboz, Gull and Troyer in conjunction with the importance sampling of the original GBMC method. This improvement allows us to carry out the importance sampling in the quantum-number-projected phase-space. Some comparisons with the previous quantum-number projection scheme indicate that, in our method, the convergence with the ground state is accelerated, which makes it possible to extend the applicability and widen the range of tractable parameters in the GBMC method. The present scheme offers an efficient practical way of computation for strongly correlated electron systems beyond the range of system sizes, interaction strengths and lattice structures tractable by other computational methods such as the quantum Monte Carlo method.

DOI
10.1143/jpsj.76.084709

Related papers

Symmetry projection schemes for Gaussian Monte Carlo methods

Fakher F. Assaad, Philipp Werner, Philippe Corboz +2 more

A sign-free Monte Carlo method for the Hubbard model has recently been proposed by Corney and Drummond. High-precision measurements on small clusters show that ground-state correlation functions are not correctly reproduced. We argue that the origin of this mismatch lies in the fact that the low-temperature density matrix does not have the symmetries of the Hamiltonian. Here we show that supplementing the algorithm with symmetry projection schemes provides reliable and accurate estimates of ground-state properties.

arXiv:cond-mat/0509149•Physical Review B•Dec 30, 2005•30 citations

Gaussian Quantum Monte Carlo methods with symmetry projection

Fakher F. Assaad

Recently Corney and Drummond introduced the sign‐free Gaussian Quantum Monte Carlo method (GQMC) for the Hubbard model. We generalize their ideas to derive stochastic differential equations (SDEs) for the general electronic structure problem. We show that the method fails to reproduce the correct symmetries of the ground state leading to systematic errors at low temperatures. We propose to restore these symmetries a posteriori by projecting the resulting density matrix onto the symmetry sectors of the ground state. The so produced results are accurate and reliable, even for parameters where an auxiliary field QMC approach fails due to a severe sign problem. Possible sources for the systematic errors are discussed.

AIP conference proceedings•Jan 1, 2006•3 citations

Systematic errors in Gaussian quantum Monte Carlo and a systematic study of the symmetry projection method

Philippe Corboz, Matthias Troyer, Adrian Kleine +3 more

Gaussian quantum Monte Carlo (GQMC) is a stochastic phase space method for fermions with positive weights. In the example of the Hubbard model close to half filling it fails to reproduce all the symmetries of the ground state leading to systematic errors at low temperatures. In a previous work [Phys. Rev. B 72, 224518 (2005)] we proposed to restore the broken symmetries by projecting the density matrix obtained from the simulation onto the ground state symmetry sector. For ground state properties, the accuracy of this method depends on a large overlap between the GQMC and exact density matrices. Thus, the method is not rigorously exact. We present the limits of the approach by a systematic study of the method for two and three leg Hubbard ladders for different fillings and on-site repulsion strengths. We show several indications that the systematic errors stem from nonvanishing boundary terms in the partial integration step in the derivation of the Fokker-Planck equation. Checking for spiking trajectories and slow decaying probability distributions provides an important test of the reliability of the results. Possible solutions to avoid boundary terms are discussed. Furthermore we compare results obtained from two different sampling methods: Reconfiguration of walkers and the Metropolis algorithm.

arXiv:0707.4394•Physical Review B•Feb 11, 2008•15 citations

Quantum Monte Carlo method for the ground state of many-boson systems

Wirawan Purwanto, Shiwei Zhang

We formulate a quantum Monte Carlo (QMC) method for calculating the ground state of many-boson systems. The method is based on a field-theoretical approach, and is closely related to existing fermion auxiliary-field QMC methods which are applied in several fields of physics. The ground-state projection is implemented as a branching random walk in the space of permanents consisting of identical single-particle orbitals. Any single-particle basis can be used, and the method is in principle exact. We illustrate this method with a trapped atomic boson gas, where the atoms interact via an attractive or repulsive contact two-body potential. We choose as the single-particle basis a real-space grid. We compare with exact results in small systems and arbitrarily sized systems of untrapped bosons with attractive interactions in one dimension, where analytical solutions exist. We also compare with the corresponding Gross-Pitaevskii (GP) mean-field calculations for trapped atoms, and discuss the close formal relation between our method and the GP approach. Our method provides a way to systematically improve upon GP while using the same framework, capturing interaction and correlation effects with a stochastic, coherent ensemble of noninteracting solutions. We discuss various algorithmic issues, including importance sampling and the back-propagation technique for computing observables, and illustrate them with numerical studies. We show results for systems with up to N approximately 400 bosons.

arXiv:physics/0403146•Physical Review E•Nov 9, 2004•75 citations

Path-Integral Renormalization Group Method for Numerical Study on Ground States of Strongly Correlated Electronic Systems

Tsuyoshi Kashima, Masatoshi Imada

A new efficient numerical algorithm for interacting fermion systems is proposed and examined in detail. The ground state is expressed approximately by a linear combination of numerically chosen basis states in a truncated Hilbert space. Two procedures lead to a better approximation. The first is a numerical renormalization, which optimizes the chosen basis and projects onto the ground state within the fixed dimension, $L$, of the Hilbert space. The second is an increase of the dimension of the truncated Hilbert space, which enables the linear combination to converge to a better approximation. The extrapolation $L\to\infty$ after the convergence removes the approximation error systematically. This algorithm does not suffer from the negative sign problem and can be applied to systems in any spatial dimension and arbitrary lattice structure. The efficiency is tested and the implementation explained for two-dimensional Hubbard models where Slater determinants are employed as chosen basis. Our results with less than 400 chosen basis indicate good accuracy within the errorbar of the best available results as those of the quantum Monte Carlo for energy and other physical quantities.

arXiv:cond-mat/0104140•Journal of the Physical Society of Japan•Aug 15, 2001•89 citations

Combining the Transcorrelated Method with Full Configuration Interaction Quantum Monte Carlo: Application to the Homogeneous Electron Gas

Hongjun Luo, Ali Alavi

We suggest an efficient method to resolve electronic cusps in electronic structure calculations through the use of an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from [Formula: see text] to [Formula: see text], where M is the total number of orbital basis functions.

Journal of Chemical Theory and Computation•Feb 12, 2018•83 citations
We examine Gaussian-basis Monte Carlo method (GBMC) introduced by Corney and Drummond. This method is based on an expansion of the density-matrix operator rho by means of the coherent Gaussian-type operator basis Lambda and does not suffer from the minus sign problem. The original method, however, often fails in reproducing the true ground state and causes systematic errors of calculated physical quantities because the samples are often trapped in some metastable or symmetry broken states. To overcome this difficulty, we combine the quantum-number projection scheme proposed by Assaad, Werner, Corboz, Gull and Troyer in conjunction with the importance sampling of the original GBMC method. This improvement allows us to carry out the importance sampling in the quantum-number-projected phase-space. Some comparisons with the previous quantum-number projection scheme indicate that, in our method, the convergence with the ground state is accelerated, which makes it possible to extend the applicability and widen the range of tractable parameters in the GBMC method. The present scheme offers an efficient practical way of computation for strongly correlated electron systems beyond the range of system sizes, interaction strengths and lattice structures tractable by other computational methods such as the quantum Monte Carlo method.
A sign-free Monte Carlo method for the Hubbard model has recently been proposed by Corney and Drummond. High-precision measurements on small clusters show that ground-state correlation functions are not correctly reproduced. We argue that the origin of this mismatch lies in the fact that the low-temperature density matrix does not have the symmetries of the Hamiltonian. Here we show that supplementing the algorithm with symmetry projection schemes provides reliable and accurate estimates of ground-state properties.
Recently Corney and Drummond introduced the sign‐free Gaussian Quantum Monte Carlo method (GQMC) for the Hubbard model. We generalize their ideas to derive stochastic differential equations (SDEs) for the general electronic structure problem. We show that the method fails to reproduce the correct symmetries of the ground state leading to systematic errors at low temperatures. We propose to restore these symmetries a posteriori by projecting the resulting density matrix onto the symmetry sectors of the ground state. The so produced results are accurate and reliable, even for parameters where an auxiliary field QMC approach fails due to a severe sign problem. Possible sources for the systematic errors are discussed.
Gaussian quantum Monte Carlo (GQMC) is a stochastic phase space method for fermions with positive weights. In the example of the Hubbard model close to half filling it fails to reproduce all the symmetries of the ground state leading to systematic errors at low temperatures. In a previous work [Phys. Rev. B 72, 224518 (2005)] we proposed to restore the broken symmetries by projecting the density matrix obtained from the simulation onto the ground state symmetry sector. For ground state properties, the accuracy of this method depends on a large overlap between the GQMC and exact density matrices. Thus, the method is not rigorously exact. We present the limits of the approach by a systematic study of the method for two and three leg Hubbard ladders for different fillings and on-site repulsion strengths. We show several indications that the systematic errors stem from nonvanishing boundary terms in the partial integration step in the derivation of the Fokker-Planck equation. Checking for spiking trajectories and slow decaying probability distributions provides an important test of the reliability of the results. Possible solutions to avoid boundary terms are discussed. Furthermore we compare results obtained from two different sampling methods: Reconfiguration of walkers and the Metropolis algorithm.
We formulate a quantum Monte Carlo (QMC) method for calculating the ground state of many-boson systems. The method is based on a field-theoretical approach, and is closely related to existing fermion auxiliary-field QMC methods which are applied in several fields of physics. The ground-state projection is implemented as a branching random walk in the space of permanents consisting of identical single-particle orbitals. Any single-particle basis can be used, and the method is in principle exact. We illustrate this method with a trapped atomic boson gas, where the atoms interact via an attractive or repulsive contact two-body potential. We choose as the single-particle basis a real-space grid. We compare with exact results in small systems and arbitrarily sized systems of untrapped bosons with attractive interactions in one dimension, where analytical solutions exist. We also compare with the corresponding Gross-Pitaevskii (GP) mean-field calculations for trapped atoms, and discuss the close formal relation between our method and the GP approach. Our method provides a way to systematically improve upon GP while using the same framework, capturing interaction and correlation effects with a stochastic, coherent ensemble of noninteracting solutions. We discuss various algorithmic issues, including importance sampling and the back-propagation technique for computing observables, and illustrate them with numerical studies. We show results for systems with up to N approximately 400 bosons.
A new efficient numerical algorithm for interacting fermion systems is proposed and examined in detail. The ground state is expressed approximately by a linear combination of numerically chosen basis states in a truncated Hilbert space. Two procedures lead to a better approximation. The first is a numerical renormalization, which optimizes the chosen basis and projects onto the ground state within the fixed dimension, LLL, of the Hilbert space. The second is an increase of the dimension of the truncated Hilbert space, which enables the linear combination to converge to a better approximation. The extrapolation L→∞L\to\inftyL→∞ after the convergence removes the approximation error systematically. This algorithm does not suffer from the negative sign problem and can be applied to systems in any spatial dimension and arbitrary lattice structure. The efficiency is tested and the implementation explained for two-dimensional Hubbard models where Slater determinants are employed as chosen basis. Our results with less than 400 chosen basis indicate good accuracy within the errorbar of the best available results as those of the quantum Monte Carlo for energy and other physical quantities.
We suggest an efficient method to resolve electronic cusps in electronic structure calculations through the use of an effective transcorrelated Hamiltonian. This effective Hamiltonian takes a simple form for plane wave bases, containing up to two-body operators only, and its use incurs almost no additional computational overhead compared to that of the original Hamiltonian. We apply this method in combination with the full configuration interaction quantum Monte Carlo (FCIQMC) method to the homogeneous electron gas. As a projection technique, the non-Hermitian nature of the transcorrelated Hamiltonian does not cause complications or numerical difficulties for FCIQMC. The rate of convergence of the total energy to the complete basis set limit is improved from [Formula: see text] to [Formula: see text], where M is the total number of orbital basis functions.