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

Avoiding ergodicity problems in lattice discretizations of the Hubbard model

Jan-Lukas Wynen, Evan Berkowitz, Christopher Körber, Timo A. Lähde, Thomas Luu

Published Aug 22, 2019• 22 citations
arXivHTMLPDF

Abstract

The Hubbard model arises naturally when electron-electron interactions are added to the tight-binding descriptions of many condensed matter systems. For instance, the two-dimensional Hubbard model on the honeycomb lattice is central to the ab initio description of the electronic structure of carbon nanomaterials, such as graphene. Such low-dimensional Hubbard models are advantageously studied with Markov chain Monte Carlo methods, such as hybrid Monte Carlo (HMC). HMC is the standard algorithm of the lattice gauge theory community, as it is well suited to theories of dynamical fermions. As HMC performs continuous, global updates of the lattice degrees of freedom, it provides superior scaling with system size relative to local updating methods. A potential drawback of HMC is its susceptibility to ergodicity problems due to so-called exceptional configurations, for which the fermion operator cannot be inverted. Recently, ergodicity problems were found in some formulations of HMC simulations of the Hubbard model. Here, we address this issue directly and clarify under what conditions ergodicity is maintained or violated in HMC simulations of the Hubbard model. We study different lattice formulations of the fermion operator and provide explicit, representative calculations for small systems, often comparing to exact results. We show that a fermion operator can be found which is both computationally convenient and free of ergodicity problems.

DOI
10.1103/physrevb.100.075141

Related papers

The Hubbard Model on the Honeycomb Lattice with Hybrid Monte Carlo

Johann Ostmeyer

We take advantage of recent improvements in the grand canonical Hybrid Monte Carlo (HMC) algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of $U_c/κ=3.835(14)$ and the critical exponent $ν=1.181(43)$ for the semimetal-antiferromagnetic Mott insulator quantum phase transition in the hexagonal Hubbard Model. Based on these results, we provide a unified, comprehensive treatment of all operators that contribute to the anti-ferromagnetic, ferromagnetic, and charge-density-wave structure factors and order parameters of the hexagonal Hubbard Model. We expect our findings to improve the consistency of Monte Carlo determinations of critical exponents. We perform a data collapse analysis and determine the critical exponent $β=0.898(37)$. We consider our findings in view of the $SU(2)$ Gross-Neveu, or chiral Heisenberg, universality class. We also discuss the computational scaling of the HMC algorithm. Our methods are applicable to a wide range of lattice theories of strongly correlated electrons. The Ising model, a simple statistical model for ferromagnetism, is one such theory. There are analytic solutions for low dimensions and very efficient Monte Carlo methods, such as cluster algorithms, for simulating this model in special cases. However most approaches do not generalise to arbitrary lattices and couplings. We present a formalism that allows one to apply HMC simulations to the Ising model, demonstrating how a system with discrete degrees of freedom can be simulated with continuous variables.

arXiv:2110.15432•Oct 28, 2021•0 citations

Fully ergodic simulations using radial updates

Finn L. Temmen, Evan Berkowitz, Anthony D. Kennedy +3 more

A sensible application of the hybrid Monte Carlo (HMC) method is often hindered by the presence of large—or even infinite—potential barriers. These potential barriers separate the configuration space into distinct sectors and can lead to ergodicity violations that bias measurements. In this work, we address this problem by augmenting HMC with a multiplicative Metropolis-Hastings update in a “radial direction” of the fields which enables crossing the potential barriers and ensures ergodicity of the sampling algorithm at comparably low computational cost. We demonstrate the algorithm on a simple toy model and show how it can be applied to the fermionic Hubbard model describing physics ranging from an exactly solvable two-site system to the C 20 H 12 perylene molecule. Our numerical results show that the radial updates successfully remove ergodicity violations, while simultaneously reducing autocorrelation times.

arXiv:2504.01575•Physical review. B./Physical review. B•Jun 24, 2025•4 citations

Hybrid Quantum Monte Carlo for Condensed Matter Models

Stefan Beyl

In this thesis we consider the hybrid quantum Monte Carlo method for simulations of the Hubbard and Su-Schrieffer-Heeger model. In the first instance, we discuss the hybrid quantum Monte Carlo method for the Hubbard model on a square lattice. We point out potential ergodicity issues and provide a way to circumvent them by a complexification of the method. Furthermore, we compare the efficiency of the hybrid quantum Monte Carlo method with a well established determinantal quantum Monte Carlo method for simulations of the half-filled Hubbard model on square lattices. One reason why the hybrid quantum Monte Carlo method loses the comparison is that we do not observe the desired sub-quadratic scaling of the numerical effort. Afterwards we present a formulation of the hybrid quantum Monte Carlo method for the Su-Schrieffer-Heeger model in two dimensions. Electron-phonon models like this are in general very hard to simulate using other Monte Carlo methods in more than one dimensions. It turns out that the hybrid quantum Monte Carlo method is much better suited for this model . We achieve favorable scaling properties and provide a proof of concept. Subsequently, we use the hybrid quantum Monte Carlo method to investigate the Su-Schrieffer-Heeger model in detail at half-filling in two dimensions. We present numerical data for staggered valence bond order at small phonon frequencies and an antiferromagnetic order at high frequencies. Due to an O(4) symmetry the antiferromagnetic order is connected to a superconducting charge density wave. Considering the Su-Schrieffer-Heeger model without tight-binding hopping reveals an additional unconstrained Z_2 gauge theory. In this case, we find indications for π-fluxes and a possible Z_2 Dirac deconfined phase as well as for a columnar valence bond ordered state at low phonon energies. In our investigations of the several phase transitions we discuss the different possibilities for the underlying mechanisms and reveal first insights into a rich phase diagram.

•Nov 5, 2019•11 citations

Revisiting the hybrid quantum Monte Carlo method for Hubbard and electron-phonon models

Stefan Beyl, Florian Goth, Fakher F. Assaad

A unique feature of the hybrid quantum Monte Carlo (HQMC) method is the potential to simulate negative sign free lattice fermion models with subcubic scaling in system size. Here we will revisit the algorithm for various models. We will show that for the Hubbard model the HQMC suffers from ergodicity issues and unbounded forces in the effective action. Solutions to these issues can be found in terms of a complexification of the auxiliary fields. This implementation of the HQMC that does not attempt to regularize the fermionic matrix so as to circumvent the aforementioned singularities does not outperform single spin flip determinantal methods with cubic scaling. On the other hand we will argue that there is a set of models for which the HQMC is very efficient. This class is characterized by effective actions free of singularities. Using the Majorana representation, we show that models such as the Su-Schrieffer-Heeger Hamiltonian at half filling and on a bipartite lattice belong to this class. For this specific model subcubic scaling is achieved.

arXiv:1708.03661•Physical review. B./Physical review. B•Feb 22, 2018•76 citations

Overcoming Ergodicity Problems of the Hybrid Monte Carlo Method using Radial Updates

Finn L. Temmen, Evan Berkowitz, Anthony D. Kennedy +3 more

Despite its many advantages, the sensible application of the Hybrid Monte Carlo (HMC) method is often hindered by the presence of large - or even infinite - potential barriers. These potential barriers partition the configuration space into distinct sectors, which leads to ergodicity violations and biased measurements of observables. In this work, we address this problem by augmenting the HMC method with a multiplicative Metropolis-Hastings update in a so-called "radial direction" of the fields, which enables jumps over the aforementioned potential barriers at comparably low computational cost. The effectiveness of this approach is demonstrated for the Hubbard model, formulated in a non-compact space by means of a continuous Hubbard-Stratonovich transformation. Our numerical results show that the radial updates successfully resolve the ergodicity violation, while simultaneously reducing autocorrelations.

arXiv:2410.19148•International Symposium on Lattice Field Theory (LATTICE)•Dec 5, 2024•3 citations

Hybrid Monte Carlo for the extended Kane-Mele-Hubbard model

Finn L. Temmen, Thomas Luu, Johann Ostmeyer +1 more

The Hubbard model offers a powerful minimal description of strongly correlated electrons in molecules and materials, yet it neglects effects observed in nature that give rise to a much richer variety of physics. For instance, including nonlocal interactions introduces a competition between strongly correlated phases, while incorporating spin-orbit coupling (SOC) can lead to the emergence of exotic topological phases. Seeking a more realistic platform for simulating graphene and, potentially, related materials such as transition metal dichalcogenides (TMDs), we therefore consider the extended Kane-Mele-Hubbard model, which features a Kane-Mele SOC term and nearest-neighbor interactions. Simulating this model is particularly challenging once the nonlocal interactions become dominant over the on-site Hubbard interaction, as a severe sign problem arises already at half filling, leaving this regime largely unexplored. In this work, we develop continuous auxiliary field formulations of the model, suitable for treatment with the Hybrid Monte Carlo (HMC) algorithm. Using these formulations, we characterize the sign problem across the model's parameter space and establish where simulations are feasible, as a first step toward the full phase diagram.

arXiv:2610.32458•Oct 23, 2026•0 citations
The Hubbard model arises naturally when electron-electron interactions are added to the tight-binding descriptions of many condensed matter systems. For instance, the two-dimensional Hubbard model on the honeycomb lattice is central to the ab initio description of the electronic structure of carbon nanomaterials, such as graphene. Such low-dimensional Hubbard models are advantageously studied with Markov chain Monte Carlo methods, such as hybrid Monte Carlo (HMC). HMC is the standard algorithm of the lattice gauge theory community, as it is well suited to theories of dynamical fermions. As HMC performs continuous, global updates of the lattice degrees of freedom, it provides superior scaling with system size relative to local updating methods. A potential drawback of HMC is its susceptibility to ergodicity problems due to so-called exceptional configurations, for which the fermion operator cannot be inverted. Recently, ergodicity problems were found in some formulations of HMC simulations of the Hubbard model. Here, we address this issue directly and clarify under what conditions ergodicity is maintained or violated in HMC simulations of the Hubbard model. We study different lattice formulations of the fermion operator and provide explicit, representative calculations for small systems, often comparing to exact results. We show that a fermion operator can be found which is both computationally convenient and free of ergodicity problems.
We take advantage of recent improvements in the grand canonical Hybrid Monte Carlo (HMC) algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of Uc/κ=3.835(14)U_c/κ=3.835(14)Uc​/κ=3.835(14) and the critical exponent ν=1.181(43)ν=1.181(43)ν=1.181(43) for the semimetal-antiferromagnetic Mott insulator quantum phase transition in the hexagonal Hubbard Model. Based on these results, we provide a unified, comprehensive treatment of all operators that contribute to the anti-ferromagnetic, ferromagnetic, and charge-density-wave structure factors and order parameters of the hexagonal Hubbard Model. We expect our findings to improve the consistency of Monte Carlo determinations of critical exponents. We perform a data collapse analysis and determine the critical exponent β=0.898(37)β=0.898(37)β=0.898(37). We consider our findings in view of the SU(2)SU(2)SU(2) Gross-Neveu, or chiral Heisenberg, universality class. We also discuss the computational scaling of the HMC algorithm. Our methods are applicable to a wide range of lattice theories of strongly correlated electrons. The Ising model, a simple statistical model for ferromagnetism, is one such theory. There are analytic solutions for low dimensions and very efficient Monte Carlo methods, such as cluster algorithms, for simulating this model in special cases. However most approaches do not generalise to arbitrary lattices and couplings. We present a formalism that allows one to apply HMC simulations to the Ising model, demonstrating how a system with discrete degrees of freedom can be simulated with continuous variables.
A sensible application of the hybrid Monte Carlo (HMC) method is often hindered by the presence of large—or even infinite—potential barriers. These potential barriers separate the configuration space into distinct sectors and can lead to ergodicity violations that bias measurements. In this work, we address this problem by augmenting HMC with a multiplicative Metropolis-Hastings update in a “radial direction” of the fields which enables crossing the potential barriers and ensures ergodicity of the sampling algorithm at comparably low computational cost. We demonstrate the algorithm on a simple toy model and show how it can be applied to the fermionic Hubbard model describing physics ranging from an exactly solvable two-site system to the C 20 H 12 perylene molecule. Our numerical results show that the radial updates successfully remove ergodicity violations, while simultaneously reducing autocorrelation times.
In this thesis we consider the hybrid quantum Monte Carlo method for simulations of the Hubbard and Su-Schrieffer-Heeger model. In the first instance, we discuss the hybrid quantum Monte Carlo method for the Hubbard model on a square lattice. We point out potential ergodicity issues and provide a way to circumvent them by a complexification of the method. Furthermore, we compare the efficiency of the hybrid quantum Monte Carlo method with a well established determinantal quantum Monte Carlo method for simulations of the half-filled Hubbard model on square lattices. One reason why the hybrid quantum Monte Carlo method loses the comparison is that we do not observe the desired sub-quadratic scaling of the numerical effort. Afterwards we present a formulation of the hybrid quantum Monte Carlo method for the Su-Schrieffer-Heeger model in two dimensions. Electron-phonon models like this are in general very hard to simulate using other Monte Carlo methods in more than one dimensions. It turns out that the hybrid quantum Monte Carlo method is much better suited for this model . We achieve favorable scaling properties and provide a proof of concept. Subsequently, we use the hybrid quantum Monte Carlo method to investigate the Su-Schrieffer-Heeger model in detail at half-filling in two dimensions. We present numerical data for staggered valence bond order at small phonon frequencies and an antiferromagnetic order at high frequencies. Due to an O(4) symmetry the antiferromagnetic order is connected to a superconducting charge density wave. Considering the Su-Schrieffer-Heeger model without tight-binding hopping reveals an additional unconstrained Z_2 gauge theory. In this case, we find indications for π-fluxes and a possible Z_2 Dirac deconfined phase as well as for a columnar valence bond ordered state at low phonon energies. In our investigations of the several phase transitions we discuss the different possibilities for the underlying mechanisms and reveal first insights into a rich phase diagram.
A unique feature of the hybrid quantum Monte Carlo (HQMC) method is the potential to simulate negative sign free lattice fermion models with subcubic scaling in system size. Here we will revisit the algorithm for various models. We will show that for the Hubbard model the HQMC suffers from ergodicity issues and unbounded forces in the effective action. Solutions to these issues can be found in terms of a complexification of the auxiliary fields. This implementation of the HQMC that does not attempt to regularize the fermionic matrix so as to circumvent the aforementioned singularities does not outperform single spin flip determinantal methods with cubic scaling. On the other hand we will argue that there is a set of models for which the HQMC is very efficient. This class is characterized by effective actions free of singularities. Using the Majorana representation, we show that models such as the Su-Schrieffer-Heeger Hamiltonian at half filling and on a bipartite lattice belong to this class. For this specific model subcubic scaling is achieved.
Despite its many advantages, the sensible application of the Hybrid Monte Carlo (HMC) method is often hindered by the presence of large - or even infinite - potential barriers. These potential barriers partition the configuration space into distinct sectors, which leads to ergodicity violations and biased measurements of observables. In this work, we address this problem by augmenting the HMC method with a multiplicative Metropolis-Hastings update in a so-called "radial direction" of the fields, which enables jumps over the aforementioned potential barriers at comparably low computational cost. The effectiveness of this approach is demonstrated for the Hubbard model, formulated in a non-compact space by means of a continuous Hubbard-Stratonovich transformation. Our numerical results show that the radial updates successfully resolve the ergodicity violation, while simultaneously reducing autocorrelations.
The Hubbard model offers a powerful minimal description of strongly correlated electrons in molecules and materials, yet it neglects effects observed in nature that give rise to a much richer variety of physics. For instance, including nonlocal interactions introduces a competition between strongly correlated phases, while incorporating spin-orbit coupling (SOC) can lead to the emergence of exotic topological phases. Seeking a more realistic platform for simulating graphene and, potentially, related materials such as transition metal dichalcogenides (TMDs), we therefore consider the extended Kane-Mele-Hubbard model, which features a Kane-Mele SOC term and nearest-neighbor interactions. Simulating this model is particularly challenging once the nonlocal interactions become dominant over the on-site Hubbard interaction, as a severe sign problem arises already at half filling, leaving this regime largely unexplored. In this work, we develop continuous auxiliary field formulations of the model, suitable for treatment with the Hybrid Monte Carlo (HMC) algorithm. Using these formulations, we characterize the sign problem across the model's parameter space and establish where simulations are feasible, as a first step toward the full phase diagram.