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:2107.11284

Determinantal Quantum Monte Carlo solver for Cluster Perturbation Theory

Edwin W. Huang, Yao Wang

Published Jul 23, 2021• 0 citations
arXivHTMLPDF

Abstract

Cluster Perturbation Theory (CPT) is a technique for computing the spectral function of fermionic models with local interactions. By combining the solution of the model on a finite cluster with perturbation theory on intra-cluster hoppings, CPT provides access to single-particle properties with arbitrary momentum resolution while incurring low computational cost. Here, we introduce Determinantal Quantum Monte Carlo (DQMC) as a solver for CPT. Compared to the standard solver, exact diagonalization (ED), the DQMC solver reduces finite size effects through utilizing larger clusters, allows study of temperature dependence, and enables large-scale simulations of a greater set of models. We discuss the implementation of the DQMC solver for CPT and benchmark the CPT+DQMC method for the attractive and repulsive Hubbard models, showcasing its advantages over standard DQMC and CPT+ED simulations.

DOI
10.48550/arxiv.2107.11284

Related papers

Determinantal quantum Monte Carlo solver for cluster perturbation theory

Edwin W. Huang, Shuhan Ding, Jiarui Liu +1 more

Cluster perturbation theory involves decomposing a 2D square lattice of interacting electrons into a set of connected clusters. By simulating a single cluster using determinantal quantum Monte Carlo, the spectral function A(k, w) for the infinite lattice can be recovered.

Physical Review Research•Oct 26, 2022•11 citations

Dynamical variational Monte Carlo as a quantum impurity solver: Application to cluster dynamical mean field theory

Peter Rosenberg, David Sénéchal, A.–M. S. Tremblay +1 more

Two of the primary sources of error in the Cluster dynamical mean field theory (CDMFT) technique arise from the use of finite size clusters and finite size baths, which makes the development of impurity solvers that can treat larger systems an essential goal. In this work we introduce an impurity solver based on the recently developed dynamical variational Monte Carlo (dVMC) method. Variational Monte Carlo possesses a favorable scaling as a function of system size, which enables the treatment of systems beyond the reach of current exact diagonalization solvers. To benchmark the technique, we perform a systematic set of CDMFT calculations on the single-band one-dimensional Hubbard model. We compare to results obtained with an exact diagonalization solver for small clusters, and against the exact solution in the thermodynamic limit obtained by Lieb and Wu for larger clusters. The development of improved impurity solvers will help extend the reach of quantum cluster methods, which can be applied to a wide range of strongly correlated electron systems, promising new insights into their emergent behavior.

arXiv:2307.15738•Physical review. B./Physical review. B•Dec 8, 2023•4 citations

Boosting determinant quantum Monte Carlo with submatrix updates: Unveiling the phase diagram of the 3D Hubbard model

Fanjie Sun, Yongbiao Xiao

Determinant Quantum Monte Carlo (DQMC) provides numerically exact solutions for strongly correlated fermionic systems but faces significant computational challenges with increasing system size. While submatrix updates were originally developed for Hirsch-Fye QMC with onsite interactions at finite temperatures, their comprehensive application in DQMC has remained unexplored despite noted algorithmic similarities. We present the first comprehensive application of submatrix updates in DQMC, significantly extending beyond the original scope by enabling simulations with extended interactions and at zero temperature. Building upon conventional fast updates and delay updates, our generalized implementation achieves an order-of-magnitude improvement in computational efficiency, enabling simulations of the half-filled Hubbard model on lattices up to 8,000 sites - a scale previously challenging with standard DQMC implementations. This enhanced computational capability allows us to accurately determine the finite-temperature phase diagram of the 3D Hubbard model at half-filling. Our findings not only shed light on the phase transitions within these complex systems but also pave the way for more effective simulations of strongly correlated electrons, potentially guiding experimental efforts in cold atom simulations of the 3D Hubbard model.

arXiv:2404.09989•SciPost Physics•Feb 17, 2025•8 citations

Applications of Cluster Perturbation Theory Using Quantum Monte Carlo Data

Fei Lin, Erik Schwartz Sørensen, Catherine Kallin +1 more

We study cluster perturbation theory [Phys. Rev. Lett. 84, 522 (2000)] when auxiliary field quantum Monte Carlo method is used for solving the cluster hamiltonian. As a case study, we calculate the spectral functions of the Hubbard model in one and two dimensions and compare our results for the spectral functions to results obtained using exact diagonalization to solve the cluster hamiltonian. The main advantage of using quantum Monte Carlo results as a starting point is that the initial cluster size can be taken to be considerably larger and hence potentially capture more of the relevant physics. The drawback is that quantum Monte Carlo methods yield results at imaginary times with stochastic errors.

arXiv:cond-mat/0512406•International Symposium on High Performance Computing Systems and Applications (HPCS)•Jan 1, 2006•0 citations

Perturbative solution of fermionic sign problem in quantum Monte Carlo computations

Sergei Iskakov, Mikhail I. Katsnelson, Alexander I. Lichtenstein

Abstract We have developed a strong-coupling perturbation scheme for a general doped Hubbard model around a particle-hole-symmetric reference system, which is free from the fermionic sign problem. Our approach is based on the lattice determinantal Quantum Monte Carlo (QMC) method in both continuous and discrete time versions for large periodic clusters in a fermionic bath. By considering the first-order perturbation in the shift of the chemical potential and the second-neighbor hopping, we are able to obtain an accurate electronic spectral function for a range of parameters that correspond to optimally doped cuprate systems at temperatures of up to T = 0.1 t , which are challenging to access using straightforward lattice QMC calculations. We also discuss the formation of the pseudogap and the nodal-antinodal dichotomy for a doped Hubbard system with the interaction parameter U equal to the bandwidth and the optimal value of the next-nearest-neighbor hopping parameter $${t}^{{\prime} }$$ t ′ for high-temperature superconducting cuprates.

arXiv:2303.01607•npj Computational Materials•Feb 15, 2024•6 citations

Characterizing the three-orbital Hubbard model with determinant quantum Monte Carlo

Y. F. Kung, Cheng-Chien Chen, Yao Wang +6 more

The authors systematically characterize the three-orbital Hubbard model using state-of-the-art determinant quantum Monte Carlo (DQMC) simulations with parameters relevant to the cuprate high-temperature superconductors. The DQMC results agree well with those from cuprate experiments, such as photoemission, and enable the identification of orbital content in the bands. A comparison of DQMC results to those from exact diagonalization and cluster perturbation theory elucidates how these different numerical techniques complement one another to produce a more complete understanding of the model and the cuprates.

arXiv:1601.05421•Physical review. B./Physical review. B•Apr 29, 2016•71 citations
Cluster Perturbation Theory (CPT) is a technique for computing the spectral function of fermionic models with local interactions. By combining the solution of the model on a finite cluster with perturbation theory on intra-cluster hoppings, CPT provides access to single-particle properties with arbitrary momentum resolution while incurring low computational cost. Here, we introduce Determinantal Quantum Monte Carlo (DQMC) as a solver for CPT. Compared to the standard solver, exact diagonalization (ED), the DQMC solver reduces finite size effects through utilizing larger clusters, allows study of temperature dependence, and enables large-scale simulations of a greater set of models. We discuss the implementation of the DQMC solver for CPT and benchmark the CPT+DQMC method for the attractive and repulsive Hubbard models, showcasing its advantages over standard DQMC and CPT+ED simulations.
Two of the primary sources of error in the Cluster dynamical mean field theory (CDMFT) technique arise from the use of finite size clusters and finite size baths, which makes the development of impurity solvers that can treat larger systems an essential goal. In this work we introduce an impurity solver based on the recently developed dynamical variational Monte Carlo (dVMC) method. Variational Monte Carlo possesses a favorable scaling as a function of system size, which enables the treatment of systems beyond the reach of current exact diagonalization solvers. To benchmark the technique, we perform a systematic set of CDMFT calculations on the single-band one-dimensional Hubbard model. We compare to results obtained with an exact diagonalization solver for small clusters, and against the exact solution in the thermodynamic limit obtained by Lieb and Wu for larger clusters. The development of improved impurity solvers will help extend the reach of quantum cluster methods, which can be applied to a wide range of strongly correlated electron systems, promising new insights into their emergent behavior.
Determinant Quantum Monte Carlo (DQMC) provides numerically exact solutions for strongly correlated fermionic systems but faces significant computational challenges with increasing system size. While submatrix updates were originally developed for Hirsch-Fye QMC with onsite interactions at finite temperatures, their comprehensive application in DQMC has remained unexplored despite noted algorithmic similarities. We present the first comprehensive application of submatrix updates in DQMC, significantly extending beyond the original scope by enabling simulations with extended interactions and at zero temperature. Building upon conventional fast updates and delay updates, our generalized implementation achieves an order-of-magnitude improvement in computational efficiency, enabling simulations of the half-filled Hubbard model on lattices up to 8,000 sites - a scale previously challenging with standard DQMC implementations. This enhanced computational capability allows us to accurately determine the finite-temperature phase diagram of the 3D Hubbard model at half-filling. Our findings not only shed light on the phase transitions within these complex systems but also pave the way for more effective simulations of strongly correlated electrons, potentially guiding experimental efforts in cold atom simulations of the 3D Hubbard model.
We study cluster perturbation theory [Phys. Rev. Lett. 84, 522 (2000)] when auxiliary field quantum Monte Carlo method is used for solving the cluster hamiltonian. As a case study, we calculate the spectral functions of the Hubbard model in one and two dimensions and compare our results for the spectral functions to results obtained using exact diagonalization to solve the cluster hamiltonian. The main advantage of using quantum Monte Carlo results as a starting point is that the initial cluster size can be taken to be considerably larger and hence potentially capture more of the relevant physics. The drawback is that quantum Monte Carlo methods yield results at imaginary times with stochastic errors.
Abstract We have developed a strong-coupling perturbation scheme for a general doped Hubbard model around a particle-hole-symmetric reference system, which is free from the fermionic sign problem. Our approach is based on the lattice determinantal Quantum Monte Carlo (QMC) method in both continuous and discrete time versions for large periodic clusters in a fermionic bath. By considering the first-order perturbation in the shift of the chemical potential and the second-neighbor hopping, we are able to obtain an accurate electronic spectral function for a range of parameters that correspond to optimally doped cuprate systems at temperatures of up to T = 0.1 t , which are challenging to access using straightforward lattice QMC calculations. We also discuss the formation of the pseudogap and the nodal-antinodal dichotomy for a doped Hubbard system with the interaction parameter U equal to the bandwidth and the optimal value of the next-nearest-neighbor hopping parameter t′{t}^{{\prime} }t′ t ′ for high-temperature superconducting cuprates.
The authors systematically characterize the three-orbital Hubbard model using state-of-the-art determinant quantum Monte Carlo (DQMC) simulations with parameters relevant to the cuprate high-temperature superconductors. The DQMC results agree well with those from cuprate experiments, such as photoemission, and enable the identification of orbital content in the bands. A comparison of DQMC results to those from exact diagonalization and cluster perturbation theory elucidates how these different numerical techniques complement one another to produce a more complete understanding of the model and the cuprates.