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

Physical Review Research

Determinantal quantum Monte Carlo solver for cluster perturbation theory

Edwin W. Huang, Shuhan Ding, Jiarui Liu, Yao Wang

Published Oct 26, 2022• 11 citations
Publisher

Abstract

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.
DOI
10.1103/physrevresearch.4.l042015

Related papers

Determinantal Quantum Monte Carlo solver for Cluster Perturbation Theory

Edwin W. Huang, Yao Wang

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.

arXiv:2107.11284•Jul 23, 2021•0 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

Spectral function of the two-dimensional Hubbard model: A density matrix renormalization group plus cluster perturbation theory study

Chun Yang, Adrian E. Feiguin

We study the spectral function of the two-dimensional Hubbard model using cluster perturbation theory, and a density matrix renormalization group as a cluster solver. We reconstruct the two-dimensional dispersion at and away from half-filling using $2\ifmmode\times\else\texttimes\fi{}L$ ladders, with $L$ up to 80 sites, yielding results with unprecedented resolution in excellent agreement with quantum Monte Carlo. The main features of the spectrum can be described with a mean-field dispersion, with kinks and pseudogap traced back to scattering between spin and charge degrees of freedom.

arXiv:1512.06102•Physical review. B./Physical review. B•Feb 24, 2016•30 citations

Quantum Monte Carlo determinantal algorithm without Hubbard-Stratonovich transformation: a general consideration

Alexey N. Rubtsov

Continuous-time determinantal algorithm is proposed for the quantum Monte Carlo simulation of the interacting fermions. The scheme does not invoke Hubbard-Stratonovich transformation. The fermionic action is divided into two parts. One of them contains the interaction and certain additional terms; another one is purely Gaussian. The first part is considered as a perturbation. Terms of the series expansion for the partition function are generated in a random walk process. The sign problem and the complexity of the algorithm are analyzed. We argue that the scheme should be useful particularly for the systems with non-local interaction.

arXiv:cond-mat/0302228•Feb 12, 2003•4 citations

Quantum Monte Carlo impurity solvers for multi-orbital problems and frequency-dependent interactions

Hiroshi Shinaoka, Fakher F. Assaad, Nils Blümer +1 more

The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.

The European Physical Journal Special Topics•Jul 1, 2017•11 citations

Cluster Perturbation Theory for 2d Ising Model

Ilya Ivantsov, Sergey G. Ovchinnikov

This paper deals with 2d Ising model in the scope of cluster perturbation theory.Ising model is defined on a two-dimensional square lattice, the amount of nearest neighbors z=4.Lattice is divided into clusters of a given size and a complete set of energy eigenvalues and eigenvectors of the cluster is defined by the diagonalization method.On the basis of this, Hubbard operators are constructed and Green function is calculated, taking into account intercluster interactions according to perturbation theory, it allows us to obtain the dependence of the magnetization on the temperature in the Hubbard-I approximation.Obtained results are compared with the exact solution of the two-dimensional Ising model.

Journal of Siberian Federal University Mathematics & Physics•Feb 1, 2015•0 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.
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.
We study the spectral function of the two-dimensional Hubbard model using cluster perturbation theory, and a density matrix renormalization group as a cluster solver. We reconstruct the two-dimensional dispersion at and away from half-filling using $2\ifmmode\times\else\texttimes\fi{}L$ ladders, with LLL up to 80 sites, yielding results with unprecedented resolution in excellent agreement with quantum Monte Carlo. The main features of the spectrum can be described with a mean-field dispersion, with kinks and pseudogap traced back to scattering between spin and charge degrees of freedom.
Continuous-time determinantal algorithm is proposed for the quantum Monte Carlo simulation of the interacting fermions. The scheme does not invoke Hubbard-Stratonovich transformation. The fermionic action is divided into two parts. One of them contains the interaction and certain additional terms; another one is purely Gaussian. The first part is considered as a perturbation. Terms of the series expansion for the partition function are generated in a random walk process. The sign problem and the complexity of the algorithm are analyzed. We argue that the scheme should be useful particularly for the systems with non-local interaction.
The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.
This paper deals with 2d Ising model in the scope of cluster perturbation theory.Ising model is defined on a two-dimensional square lattice, the amount of nearest neighbors z=4.Lattice is divided into clusters of a given size and a complete set of energy eigenvalues and eigenvectors of the cluster is defined by the diagonalization method.On the basis of this, Hubbard operators are constructed and Green function is calculated, taking into account intercluster interactions according to perturbation theory, it allows us to obtain the dependence of the magnetization on the temperature in the Hubbard-I approximation.Obtained results are compared with the exact solution of the two-dimensional Ising model.