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:1606.04785 · Physical review. B./Physical review. B

Computation of dynamical correlation functions for many-fermion systems with auxiliary-field quantum Monte Carlo

Ettore Vitali, Hao Shi, Mingpu Qin, Shiwei Zhang

Published Aug 23, 2016• 56 citations
arXivHTMLPDF

Abstract

Quantum Monte Carlo simulations provide a powerful tool to compute static properties of quantum many-particles models and realistic systems. A vast array of ground-state energies and properties defined by equal-time correlation functions can be computed with high accuracy. On the other hand, the extensions needed to study time-displaced correlation functions, which provide access to excited states, have been challenging, especially for fermionic strongly correlated systems. The authors present a methodology that allows a high-accuracy calculation of the dynamical Green functions of many-fermions systems in the framework of the auxiliary-field quantum Monte Carlo method. They extract the the charge gap of the repulsive two-dimensional Hubbard model at half-filling as a function of the interaction strength and propose strategies for applying their methodology to more realistic systems.

DOI
10.1103/physrevb.94.085140

Related papers

Hamiltonian symmetries in auxiliary-field quantum Monte Carlo calculations for electronic structure

Mário Motta, Shiwei Zhang, Garnet Kin‐Lic Chan

Quantum Monte Carlo simulations provide a powerful tool to compute static properties of quantum many-particle models and realistic systems. Here, the authors present a technique to reduce the computational cost of computing ground-state energies, properties, and correlation functions in the framework of the auxiliary-field quantum Monte Carlo method. This is based on a simple application of group theory, and can be readily transferred to other computational methodologies for many-particle systems.

arXiv:1905.00511•Physical review. B./Physical review. B•Jul 19, 2019•27 citations

Real-time dynamics of strongly correlated fermions using auxiliary field quantum Monte Carlo

Matthew S. Church, Brenda M. Rubenstein

Spurred by recent technological advances, there is a growing demand for computational methods that can accurately predict the dynamics of correlated electrons. Such methods can provide much-needed theoretical insights into the electron dynamics probed via time-resolved spectroscopy experiments and observed in non-equilibrium ultracold atom experiments. In this article, we develop and benchmark a numerically exact Auxiliary Field Quantum Monte Carlo (AFQMC) method for modeling the dynamics of correlated electrons in real time. AFQMC has become a powerful method for predicting the ground state and finite temperature properties of strongly correlated systems mostly by employing constraints to control the sign problem. Our initial goal in this work is to determine how well AFQMC generalizes to real-time electron dynamics problems without constraints. By modeling the repulsive Hubbard model on different lattices and with differing initial electronic configurations, we show that real-time AFQMC is capable of accurately capturing long-lived electronic coherences beyond the reach of mean field techniques. While the times to which we can meaningfully model decrease with increasing correlation strength and system size as a result of the exponential growth of the dynamical phase problem, we show that our technique can model the short-time behavior of strongly correlated systems to very high accuracy. Crucially, we find that importance sampling, combined with a novel adaptive active space sampling technique, can substantially lengthen the times to which we can simulate. These results establish real-time AFQMC as a viable technique for modeling the dynamics of correlated electron systems and serve as a basis for future sampling advances that will further mitigate the dynamical phase problem.

The Journal of Chemical Physics•May 10, 2021•14 citations

Efficient calculation of imaginary-time-displaced correlation functions in the projector auxiliary-field quantum Monte Carlo algorithm

M. Feldbacher, Fakher F. Assaad

The calculation of imaginary-time-displaced correlation functions with the auxiliary-field projector quantum Monte Carlo algorithm provides valuable insight (such as spin and charge gaps) into the model under consideration. One of the authors and M. Imada proposed a numerically stable method to compute those quantities [J. Phys. Soc. Jpn. 65, 189 (1996)]. Although precise, this method is expensive in CPU time. Here we present an alternative approach which is an order of magnitude quicker, just as precise, and very simple to implement. The method is based on the observation that for a given auxiliary field the equal-time Green-function matrix G is a projector: ${G}^{2}=G.$

arXiv:cond-mat/0009447•Physical review. B, Condensed matter•Jan 29, 2001•38 citations

Imaginary Time Correlations for a High-Density two-dimensional Electron Gas

Mário Motta, Davide Emilio Galli, Saverio Moroni +1 more

We evaluate imaginary time density-density correlation functions for a two-dimensional homogeneous electron gas using the phaseless auxiliary field quantum Monte Carlo method. We show that such methodology, once equipped with suitable numerical stabilization techniques necessary to deal with exponentials, products and inversions of large matrices, gives access to the calculation of imaginary time correlation functions for medium-sized systems; we present simulations of a number up to 42 correlated fermions in the continuum, using up to 300 plane waves as basis set elements. We discuss the numerical stabilization techniques and the computational complexity of the methodology. We perform the inverse Laplace transform of the obtained density-density correlation functions, assessing the ability of the phaseless auxiliary field quantum Monte Carlo method to evaluate dynamical properties of many-fermion systems.

arXiv:1509.01219•Jul 28, 2015•0 citations

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Mingpu Qin, Hao Shi, Shiwei Zhang

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of $U/t$ from 2 to 8, for a range of densities including half-filling and $n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75$, and $0.875$. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.

arXiv:1605.09421•Physical review. B./Physical review. B•Aug 4, 2016•100 citations

Computational Studies Of Strongly Correlated Quantum Matter

Hao Shi

The study of strongly correlated quantum many-body systems is an outstanding challenge. Highly accurate results are needed for the understanding of practical and fundamental problems in condensed-matter physics, high energy physics, material science, quantum chemistry and so on. Our familiar mean-field or perturbative methods tend to be ineffective. Numerical simulations provide a promising approach for studying such systems. The fundamental difficulty of numerical simulation is that the dimension of the Hilbert space needed to describe interacting systems increases exponentially with the system size. Quantum Monte Carlo (QMC) methods are one of the best approaches to tackle the problem of enormous Hilbert space. They have been highly successful for boson systems and unfrustrated spin models. For systems with fermions, the exchange symmetry in general causes the infamous sign problem, making the statistical noise in the computed results grow exponentially with the system size. This hinders our understanding of interesting physics such as high-temperature superconductivity, metal-insulator phase transition. In this thesis, we present a variety of new developments in the auxiliary-field quantum Monte Carlo (AFQMC) methods, including the incorporation of symmetry in both the trial wave function and the projector, developing the constraint release method, using the force-bias to drastically improve the efficiency in Metropolis framework, identifying and solving the infinite variance problem, and sampling Hartree-Fock-Bogoliubov wave function. With these developments, some of the most challenging many-electron problems are now under control. We obtain an exact numerical solution of two-dimensional strongly interacting Fermi atomic gas, determine the ground state properties of the 2D Fermi gas with Rashba spin-orbit coupling, provide benchmark results for the ground state of the two-dimensional Hubbard model, and establish that the Hubbard model has a stripe order in the underdoped region.

•Jan 1, 2017•0 citations
Quantum Monte Carlo simulations provide a powerful tool to compute static properties of quantum many-particles models and realistic systems. A vast array of ground-state energies and properties defined by equal-time correlation functions can be computed with high accuracy. On the other hand, the extensions needed to study time-displaced correlation functions, which provide access to excited states, have been challenging, especially for fermionic strongly correlated systems. The authors present a methodology that allows a high-accuracy calculation of the dynamical Green functions of many-fermions systems in the framework of the auxiliary-field quantum Monte Carlo method. They extract the the charge gap of the repulsive two-dimensional Hubbard model at half-filling as a function of the interaction strength and propose strategies for applying their methodology to more realistic systems.
Quantum Monte Carlo simulations provide a powerful tool to compute static properties of quantum many-particle models and realistic systems. Here, the authors present a technique to reduce the computational cost of computing ground-state energies, properties, and correlation functions in the framework of the auxiliary-field quantum Monte Carlo method. This is based on a simple application of group theory, and can be readily transferred to other computational methodologies for many-particle systems.
Spurred by recent technological advances, there is a growing demand for computational methods that can accurately predict the dynamics of correlated electrons. Such methods can provide much-needed theoretical insights into the electron dynamics probed via time-resolved spectroscopy experiments and observed in non-equilibrium ultracold atom experiments. In this article, we develop and benchmark a numerically exact Auxiliary Field Quantum Monte Carlo (AFQMC) method for modeling the dynamics of correlated electrons in real time. AFQMC has become a powerful method for predicting the ground state and finite temperature properties of strongly correlated systems mostly by employing constraints to control the sign problem. Our initial goal in this work is to determine how well AFQMC generalizes to real-time electron dynamics problems without constraints. By modeling the repulsive Hubbard model on different lattices and with differing initial electronic configurations, we show that real-time AFQMC is capable of accurately capturing long-lived electronic coherences beyond the reach of mean field techniques. While the times to which we can meaningfully model decrease with increasing correlation strength and system size as a result of the exponential growth of the dynamical phase problem, we show that our technique can model the short-time behavior of strongly correlated systems to very high accuracy. Crucially, we find that importance sampling, combined with a novel adaptive active space sampling technique, can substantially lengthen the times to which we can simulate. These results establish real-time AFQMC as a viable technique for modeling the dynamics of correlated electron systems and serve as a basis for future sampling advances that will further mitigate the dynamical phase problem.
The calculation of imaginary-time-displaced correlation functions with the auxiliary-field projector quantum Monte Carlo algorithm provides valuable insight (such as spin and charge gaps) into the model under consideration. One of the authors and M. Imada proposed a numerically stable method to compute those quantities [J. Phys. Soc. Jpn. 65, 189 (1996)]. Although precise, this method is expensive in CPU time. Here we present an alternative approach which is an order of magnitude quicker, just as precise, and very simple to implement. The method is based on the observation that for a given auxiliary field the equal-time Green-function matrix G is a projector: G2=G.{G}^{2}=G.G2=G.
We evaluate imaginary time density-density correlation functions for a two-dimensional homogeneous electron gas using the phaseless auxiliary field quantum Monte Carlo method. We show that such methodology, once equipped with suitable numerical stabilization techniques necessary to deal with exponentials, products and inversions of large matrices, gives access to the calculation of imaginary time correlation functions for medium-sized systems; we present simulations of a number up to 42 correlated fermions in the continuum, using up to 300 plane waves as basis set elements. We discuss the numerical stabilization techniques and the computational complexity of the methodology. We perform the inverse Laplace transform of the obtained density-density correlation functions, assessing the ability of the phaseless auxiliary field quantum Monte Carlo method to evaluate dynamical properties of many-fermion systems.
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/tU/tU/t from 2 to 8, for a range of densities including half-filling and n=0.3,0.5,0.6,0.75n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75n=0.3,0.5,0.6,0.75, and 0.8750.8750.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.
The study of strongly correlated quantum many-body systems is an outstanding challenge. Highly accurate results are needed for the understanding of practical and fundamental problems in condensed-matter physics, high energy physics, material science, quantum chemistry and so on. Our familiar mean-field or perturbative methods tend to be ineffective. Numerical simulations provide a promising approach for studying such systems. The fundamental difficulty of numerical simulation is that the dimension of the Hilbert space needed to describe interacting systems increases exponentially with the system size. Quantum Monte Carlo (QMC) methods are one of the best approaches to tackle the problem of enormous Hilbert space. They have been highly successful for boson systems and unfrustrated spin models. For systems with fermions, the exchange symmetry in general causes the infamous sign problem, making the statistical noise in the computed results grow exponentially with the system size. This hinders our understanding of interesting physics such as high-temperature superconductivity, metal-insulator phase transition. In this thesis, we present a variety of new developments in the auxiliary-field quantum Monte Carlo (AFQMC) methods, including the incorporation of symmetry in both the trial wave function and the projector, developing the constraint release method, using the force-bias to drastically improve the efficiency in Metropolis framework, identifying and solving the infinite variance problem, and sampling Hartree-Fock-Bogoliubov wave function. With these developments, some of the most challenging many-electron problems are now under control. We obtain an exact numerical solution of two-dimensional strongly interacting Fermi atomic gas, determine the ground state properties of the 2D Fermi gas with Rashba spin-orbit coupling, provide benchmark results for the ground state of the two-dimensional Hubbard model, and establish that the Hubbard model has a stripe order in the underdoped region.