The Hubbard Model on the Honeycomb Lattice with Hybrid Monte Carlo

Abstract

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/\kappa =3.835(14)$ and the critical exponent $\nu =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 $\beta =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.

Cumulative PhD thesis (HISKP, Bonn University) based on the publications arXiv:1912.03278, arXiv:2004.10472, arXiv:2005.11112, arXiv:2105.06936

References in corpus (17)

• Many-Body Physics with Ultracold Gases
• Theory of interacting electrons on the honeycomb lattice
• Light hadrons from lattice QCD with light (u,d), strange and charm dynamical quarks
• Quantum critical point in graphene approached in the limit of infinitely strong Coulomb interaction
• Fermionic quantum criticality in honeycomb and $\pi$-flux Hubbard models: Finite-size scaling of renormalization-group-invariant observables from quantum Monte Carlo
• Four-loop critical exponents for the Gross-Neveu-Yukawa models
• Machine Learning of Explicit Order Parameters: From the Ising Model to SU(2) Lattice Gauge Theory
• Quantum Critical Behaviour in a Graphene-like Model
• Structure of Lefschetz thimbles in simple fermionic systems
• Monte-Carlo study of the electron transport properties of monolayer graphene within the tight-binding model
• Momentum Distribution and Contact of the Unitary Fermi gas
• Deep Learning Beyond Lefschetz Thimbles
• Dirac electrons in the square lattice Hubbard model with a $d$-wave pairing field: chiral Heisenberg universality class revisited
• Quantum critical scaling in magnetic field near the Dirac point in graphene
• Hybrid Monte Carlo Simulation of Graphene on the Hexagonal Lattice
• Automated label flows for excited states of correlation functions in lattice gauge theory
• Data collapse in the critical region using finite-size scaling with subleading corrections