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:2006.10769 · Physical Review X

Tracking the Footprints of Spin Fluctuations: A MultiMethod, MultiMessenger Study of the Two-Dimensional Hubbard Model

Thomas Schäfer, Nils Wentzell, Fedor Šimkovic, Yuan-Yao He, Cornelia Hille, Marcel Klett, Christian Eckhardt, Behnam Arzhang, Viktor Harkov, François-Marie Le Régent, Alfred Kirsch, Yan Wang, Aaram J. Kim, Evgeny Kozik, E. A. Stepanov, Anna Kauch, Sabine Andergassen, P. Hansmann, Daniel Rohe, Yuri M. Vilk, J. P. F. LeBlanc, Shiwei Zhang, A.–M. S. Tremblay, Michel Ferrero, Olivier Parcollet, Antoine Georges

Published Mar 23, 2021• 222 citations
arXivHTMLPDF

Abstract

The Hubbard model represents the fundamental model for interacting quantum systems and electronic correlations.Using the two-dimensional half-filled Hubbard model at weak coupling as a testing ground, we perform a comparative study of a comprehensive set of state-of-the-art quantum many-body methods.Upon cooling into its insulating antiferromagnetic ground state, the model hosts a rich sequence of distinct physical regimes with crossovers between a high-temperature incoherent regime, an intermediate-temperature metallic regime, and a low-temperature insulating regime with a pseudogap created by antiferromagnetic fluctuations.We assess the ability of each method to properly address these physical regimes and crossovers through the computation of several observables probing both quasiparticle properties and magnetic correlations, with two numerically exact methods (diagrammatic and determinantal quantum Monte Carlo methods) serving as a benchmark.By combining computational results and analytical insights, we elucidate the nature and role of spin fluctuations in each of these regimes.Based on this analysis, we explain how quasiparticles can coexist with increasingly longrange antiferromagnetic correlations and why dynamical mean-field theory is found to provide a remarkably accurate approximation of local quantities in the metallic regime.We also critically

DOI
10.1103/physrevx.11.011058

Related papers

Antiferromagnetic pseudogap in the two-dimensional Hubbard model deep in the renormalized classical regime

Yu. M. Vil’k, C. Lahaie, A.–M. S. Tremblay

Long-wavelength spin fluctuations prohibit antiferromagnetic long-range order at finite temperature in two dimensions. Nevertheless, the correlation length starts to grow rapidly at a crossover temperature, leading to critical slowing down and to a renormalized-classical regime over a wide range of temperature, between a fraction of the mean-field transition temperature and the zero-temperature ordered state. This leads to a single-particle pseudogap of the kind observed in electron-doped cuprates. Very few theoretical methods can claim an accurate description of this regime. The challenge is that in this regime Fermi-liquid quasiparticles are already destroyed and new quasiparticles of the ordered state are not fully formed yet. Here, we study this problem for the two-dimensional Hubbard model by first generalizing the two-particle self-consistent approach. Using a special algorithm, spin fluctuations are treated in the thermodynamic limit even for large correlation lengths. The effects of Kanamori-Br\"uckner screening, of classical and of quantum fluctuations are taken into account. Results are presented at half filling for the one-band Hubbard model with nearest-neighbor hopping. They agree well with available benchmark diagrammatic quantum Monte Carlo at high temperature where the pseudogap opens up. In addition to temperature-dependent spectral properties, we find quantum corrections to the zero-temperature renormalized mean-field antiferromagnetic gap. Finally, analytic continuation of the Matsubara results for spectral functions show that the pseudogap opens up at significantly higher temperature than was previously identified based on the Matsubara data only.

arXiv:2407.11803•Physical review. B./Physical review. B•Sep 27, 2024•3 citations

Local and Nonlocal Electronic Correlations at the Metal-Insulator Transition in the Two-Dimensional Hubbard Model

Maria Chatzieleftheriou, Silke Biermann, E. A. Stepanov

Elucidating the physics of the single-orbital Hubbard model in its intermediate-coupling regime is a key missing ingredient to our understanding of metal-insulator transitions in real materials. Using recent nonperturbative many-body techniques that are able to interpolate between the spin-fluctuation-dominated Slater regime at weak coupling and the Mott insulator at strong coupling, we obtain the momentum-resolved spectral function in the intermediate regime and disentangle the effects of antiferromagnetic fluctuations and local electronic correlations in the formation of an insulating state. This allows us to identify the Slater and Heisenberg regimes in the phase diagram, which are separated by a crossover region of competing spatial and local electronic correlations. We identify the crossover regime by investigating the behavior of the local magnetic moment, shedding light on the formation of the insulating state at intermediate couplings.

arXiv:2312.03123•Physical Review Letters•Jun 5, 2024•20 citations

Two-dimensional Hubbard model at finite temperature: Weak, strong, and long correlation regimes

Fedor Šimkovic, Riccardo Rossi, Michel Ferrero

We investigate the momentum-resolved spin and charge susceptibilities, as well as the chemical potential and double occupancy in the two-dimensional Hubbard model as functions of doping, temperature, and interaction strength. Through these quantities, we identify a weak-coupling regime, a strong-coupling regime with short-range correlations and an intermediate-coupling regime with long magnetic correlation lengths. In the spin channel, we observe an additional crossover from commensurate to incommensurate correlations. In contrast, we find charge correlations to be only short ranged for all studied temperatures, which suggests that the spin and charge responses are decoupled. These findings were obtained by a connected determinant diagrammatic Monte Carlo algorithm for the computation of double expansions, which we introduce in this paper. This permits us to obtain numerically exact results at temperatures as low as $T=0.067$ and interactions up to $U=7$, while working on arbitrarily large lattices. Our method also allows us to gain physical insights from investigating the analytic structure of perturbative series. We connect to previous work by studying smaller lattice geometries and report substantial finite-size effects.

arXiv:2110.05863•Physical Review Research•Dec 20, 2022•45 citations

Extended Metal-Insulator Crossover with Strong Antiferromagnetic Spin Correlation in Half-Filled 3D Hubbard Model

Yufeng Song, Youjin Deng, Yuan-Yao He

The Hubbard model at temperatures above the Néel transition, despite being a paramagnet, can exhibit rich physics due to the interplay of Fermi surface, on-site interaction U and thermal fluctuations. Nevertheless, the understanding of the crossover physics remains only at a qualitative level, because of the intrinsically smooth behavior. Employing an improved variant of the numerically exact auxiliary-field quantum Monte Carlo algorithm equipped with numerical analytic continuation, we obtain a broad variety of thermodynamic and dynamical properties of the three-dimensional Hubbard model at half filling, quantitatively determine the crossover boundaries, and observe that the metal-insulator crossover regime, in which antiferromagnetic spin correlations appear strongest, exists over an extended regime in between the Fermi liquid for small U and the Mott insulator for large U. In particular, the location of the most rapid suppression of double occupancy as U increases, is found to fully reside in the metallic Fermi liquid regime, in contrast to the conventional intuition that it is a representative feature for entering the Mott insulator. Beside providing a reliable methodology for numerical study of crossover physics, our work can also serve as a timely and important guideline for the most recent optical lattice experiments.

arXiv:2404.08745•Physical Review Letters•Jan 10, 2025•17 citations

Spin and Charge Correlations across the Metal-to-Insulator Crossover in the Half-Filled 2D Hubbard Model

Aaram J. Kim, Fedor Šimkovic, Evgeny Kozik

The 2D Hubbard model with nearest-neighbor hopping on the square lattice and an average of one electron per site is known to undergo an extended crossover from metallic to insulating behavior driven by proliferating antiferromagnetic correlations. We study signatures of this crossover in spin and charge correlation functions and present results obtained with controlled accuracy using the diagrammatic Monte Carlo approach in the range of parameters amenable to experimental verification with ultracold atoms in optical lattices. The qualitative changes in charge and spin correlations associated with the crossover are observed at well-separated temperature scales, which encase the intermediary regime of non-Fermi-liquid character, where local magnetic moments are formed and nonlocal fluctuations in both channels are essential.

arXiv:1905.13337•Physical Review Letters•Mar 17, 2020•58 citations

Spin fluctuations and nearly isotropic quasiparticles in the doped three-dimensional Hubbard model

Pengcheng Hou, Kun Chen, Youjin Deng

Away from half filling and above its Neel temperature, the three-dimensional Hubbard model is a paramagnetic metal with antiferromagnetic spin fluctuations, where controlled momentum-resolved results have been lacking. Using connected-determinant diagrammatic Monte Carlo, we compute the momentum-resolved self-energy for n = 0.4-0.95, T/t = 0.05-0.33 and U/t = 4-8, and the full-zone spin susceptibility at twelve of these points, focusing on U/t = 4. At U/t = 4 and toward half filling the spin fluctuations become strong: for n = 0.7-0.9 at T/t = 0.14 they enhance the susceptibility peak up to sevenfold and raise the momentum-dependent part of the scattering rate to four to five times its second-order value. They are driven by an effective spin interaction screened to about 2U/3 and nearly uniform over the zone. Yet the peak height is set by the slightly smaller value of this interaction at the peak wave vector, and near magnetic order this small difference changes the height by a large factor. Two-particle self-consistent theory reproduces the peak because its sum rule fixes its single constant interaction near this value, but for n = 0.7-0.9 at T/t = 0.14 its self-energy overestimates the local scattering rate and underestimates the momentum-dependent part. Though strong, the fluctuations are short ranged, with a correlation length of at most about two lattice spacings, and for n = 0.7-0.9 and T/t >= 0.14 the scattering rate near the Fermi surface varies by at most about 1% at fixed band energy, without hot spots. At U/t = 4 dynamical mean-field theory (DMFT) reproduces the local scattering rate and quasiparticle weight within a few percent, and the nonlocal part it misses is approximately a function of the band energy. Strong spin fluctuations thus coexist with nearly isotropic quasiparticles, and the results provide a benchmark for extensions of DMFT and for spin-fluctuation theories.

arXiv:2610.27779•Oct 20, 2026•0 citations
The Hubbard model represents the fundamental model for interacting quantum systems and electronic correlations.Using the two-dimensional half-filled Hubbard model at weak coupling as a testing ground, we perform a comparative study of a comprehensive set of state-of-the-art quantum many-body methods.Upon cooling into its insulating antiferromagnetic ground state, the model hosts a rich sequence of distinct physical regimes with crossovers between a high-temperature incoherent regime, an intermediate-temperature metallic regime, and a low-temperature insulating regime with a pseudogap created by antiferromagnetic fluctuations.We assess the ability of each method to properly address these physical regimes and crossovers through the computation of several observables probing both quasiparticle properties and magnetic correlations, with two numerically exact methods (diagrammatic and determinantal quantum Monte Carlo methods) serving as a benchmark.By combining computational results and analytical insights, we elucidate the nature and role of spin fluctuations in each of these regimes.Based on this analysis, we explain how quasiparticles can coexist with increasingly longrange antiferromagnetic correlations and why dynamical mean-field theory is found to provide a remarkably accurate approximation of local quantities in the metallic regime.We also critically
Long-wavelength spin fluctuations prohibit antiferromagnetic long-range order at finite temperature in two dimensions. Nevertheless, the correlation length starts to grow rapidly at a crossover temperature, leading to critical slowing down and to a renormalized-classical regime over a wide range of temperature, between a fraction of the mean-field transition temperature and the zero-temperature ordered state. This leads to a single-particle pseudogap of the kind observed in electron-doped cuprates. Very few theoretical methods can claim an accurate description of this regime. The challenge is that in this regime Fermi-liquid quasiparticles are already destroyed and new quasiparticles of the ordered state are not fully formed yet. Here, we study this problem for the two-dimensional Hubbard model by first generalizing the two-particle self-consistent approach. Using a special algorithm, spin fluctuations are treated in the thermodynamic limit even for large correlation lengths. The effects of Kanamori-Br\"uckner screening, of classical and of quantum fluctuations are taken into account. Results are presented at half filling for the one-band Hubbard model with nearest-neighbor hopping. They agree well with available benchmark diagrammatic quantum Monte Carlo at high temperature where the pseudogap opens up. In addition to temperature-dependent spectral properties, we find quantum corrections to the zero-temperature renormalized mean-field antiferromagnetic gap. Finally, analytic continuation of the Matsubara results for spectral functions show that the pseudogap opens up at significantly higher temperature than was previously identified based on the Matsubara data only.
Elucidating the physics of the single-orbital Hubbard model in its intermediate-coupling regime is a key missing ingredient to our understanding of metal-insulator transitions in real materials. Using recent nonperturbative many-body techniques that are able to interpolate between the spin-fluctuation-dominated Slater regime at weak coupling and the Mott insulator at strong coupling, we obtain the momentum-resolved spectral function in the intermediate regime and disentangle the effects of antiferromagnetic fluctuations and local electronic correlations in the formation of an insulating state. This allows us to identify the Slater and Heisenberg regimes in the phase diagram, which are separated by a crossover region of competing spatial and local electronic correlations. We identify the crossover regime by investigating the behavior of the local magnetic moment, shedding light on the formation of the insulating state at intermediate couplings.
We investigate the momentum-resolved spin and charge susceptibilities, as well as the chemical potential and double occupancy in the two-dimensional Hubbard model as functions of doping, temperature, and interaction strength. Through these quantities, we identify a weak-coupling regime, a strong-coupling regime with short-range correlations and an intermediate-coupling regime with long magnetic correlation lengths. In the spin channel, we observe an additional crossover from commensurate to incommensurate correlations. In contrast, we find charge correlations to be only short ranged for all studied temperatures, which suggests that the spin and charge responses are decoupled. These findings were obtained by a connected determinant diagrammatic Monte Carlo algorithm for the computation of double expansions, which we introduce in this paper. This permits us to obtain numerically exact results at temperatures as low as T=0.067T=0.067T=0.067 and interactions up to U=7U=7U=7, while working on arbitrarily large lattices. Our method also allows us to gain physical insights from investigating the analytic structure of perturbative series. We connect to previous work by studying smaller lattice geometries and report substantial finite-size effects.
The Hubbard model at temperatures above the Néel transition, despite being a paramagnet, can exhibit rich physics due to the interplay of Fermi surface, on-site interaction U and thermal fluctuations. Nevertheless, the understanding of the crossover physics remains only at a qualitative level, because of the intrinsically smooth behavior. Employing an improved variant of the numerically exact auxiliary-field quantum Monte Carlo algorithm equipped with numerical analytic continuation, we obtain a broad variety of thermodynamic and dynamical properties of the three-dimensional Hubbard model at half filling, quantitatively determine the crossover boundaries, and observe that the metal-insulator crossover regime, in which antiferromagnetic spin correlations appear strongest, exists over an extended regime in between the Fermi liquid for small U and the Mott insulator for large U. In particular, the location of the most rapid suppression of double occupancy as U increases, is found to fully reside in the metallic Fermi liquid regime, in contrast to the conventional intuition that it is a representative feature for entering the Mott insulator. Beside providing a reliable methodology for numerical study of crossover physics, our work can also serve as a timely and important guideline for the most recent optical lattice experiments.
The 2D Hubbard model with nearest-neighbor hopping on the square lattice and an average of one electron per site is known to undergo an extended crossover from metallic to insulating behavior driven by proliferating antiferromagnetic correlations. We study signatures of this crossover in spin and charge correlation functions and present results obtained with controlled accuracy using the diagrammatic Monte Carlo approach in the range of parameters amenable to experimental verification with ultracold atoms in optical lattices. The qualitative changes in charge and spin correlations associated with the crossover are observed at well-separated temperature scales, which encase the intermediary regime of non-Fermi-liquid character, where local magnetic moments are formed and nonlocal fluctuations in both channels are essential.
Away from half filling and above its Neel temperature, the three-dimensional Hubbard model is a paramagnetic metal with antiferromagnetic spin fluctuations, where controlled momentum-resolved results have been lacking. Using connected-determinant diagrammatic Monte Carlo, we compute the momentum-resolved self-energy for n = 0.4-0.95, T/t = 0.05-0.33 and U/t = 4-8, and the full-zone spin susceptibility at twelve of these points, focusing on U/t = 4. At U/t = 4 and toward half filling the spin fluctuations become strong: for n = 0.7-0.9 at T/t = 0.14 they enhance the susceptibility peak up to sevenfold and raise the momentum-dependent part of the scattering rate to four to five times its second-order value. They are driven by an effective spin interaction screened to about 2U/3 and nearly uniform over the zone. Yet the peak height is set by the slightly smaller value of this interaction at the peak wave vector, and near magnetic order this small difference changes the height by a large factor. Two-particle self-consistent theory reproduces the peak because its sum rule fixes its single constant interaction near this value, but for n = 0.7-0.9 at T/t = 0.14 its self-energy overestimates the local scattering rate and underestimates the momentum-dependent part. Though strong, the fluctuations are short ranged, with a correlation length of at most about two lattice spacings, and for n = 0.7-0.9 and T/t >= 0.14 the scattering rate near the Fermi surface varies by at most about 1% at fixed band energy, without hot spots. At U/t = 4 dynamical mean-field theory (DMFT) reproduces the local scattering rate and quasiparticle weight within a few percent, and the nonlocal part it misses is approximately a function of the band energy. Strong spin fluctuations thus coexist with nearly isotropic quasiparticles, and the results provide a benchmark for extensions of DMFT and for spin-fluctuation theories.