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

The zoo of states in the Hubbard model

Robin Scholle

Published Jan 1, 2026• 0 citations
Publisher

Abstract

The two-dimensional Hubbard model is one of the most fundamental models in condensed matter physics. It captures many properties of cuprate high-temperature superconductors. Despite its apparent simplicity, solving the Hubbard model is a challenging task. In this work, we employ the Hartree-Fock approximation, as well as functional renormalization group methods, to construct phase diagrams of the Hubbard model and to shed light on the zoo of magnetic, charge ordered, and superconducting states. The Hartree-Fock approximation can self-consistently determine the expectation values of local charge and magnetic order parameters. Mean-field theories like this can lead to a deeper understanding of the competition, coexistence and transition of different orders. We start by employing the Hartree-Fock approximation in a completely unrestricted form on finite-size lattices and complement the results with calculations in the thermodynamic limit. We construct mean-field phase diagrams of the two-dimensional Hubbard model in the thermodynamic limit for different values of the next-nearest neighbor hopping as a function of filling and temperature. We find that the dominant magnetic orders at moderate coupling strength are Néel order around half-filling, followed by spiral order, and stripe order as the hole doping is increased. We also clarify the nature of the transition from magnetic spiral to stripe order. We then extend the Hartree-Fock scheme by a functional renormalization group flow, which generates new effective interactions and allows us to also access the superconducting phase of the Hubbard model. We present an improved phase diagram, featuring Néel, spiral, and stripe order coexisting with d-wave superconductivity at low temperatures. We analyze the coexistence and competition between these orders. Finally, we perform unrestricted Hartree-Fock calculations of a Hubbard model on the graphene lattice and the Bernal bilayer graphene lattice. We also construct magnetic phase diagrams for these geometries, which feature Néel order around half-filling and various stripe orders in a secondary magnetic dome.

DOI
10.18419/opus-18665

Related papers

Comprehensive mean-field analysis of magnetic and charge orders in the two-dimensional Hubbard model

Robin Scholle, Pietro M. Bonetti, Demetrio Vilardi +1 more

We present an unbiased mean-field analysis of magnetic and charge orders in the two-dimensional Hubbard model on a square lattice, at both zero and finite temperatures. Unrestricted Hartree-Fock calculations on large finite lattices are complemented by solutions restricted to N\'eel and circular spiral order in the thermodynamic limit. The magnetic states are classified by a systematic scheme based on the dominant Fourier components of the spin texture. On finite lattices a whole zoo of ordering patterns appears. We show that many of these states are finite-size artifacts related to the limited choice of ordering wave vectors on a finite lattice. In the thermodynamic limit only three classes of states with a relatively simple structure survive: N\'eel, circular spiral, and stripe states. Stripes involve also charge order and can be unidirectional or bidirectional, with horizontal and/or vertical orientation. We present complete phase diagrams in the plane spanned by electron density and temperature, for a moderate Hubbard interaction and various choices of the next-nearest-neighbor hopping amplitude.

arXiv:2303.15358•Physical review. B./Physical review. B•Jul 19, 2023•33 citations

Spiral to stripe transition in the two-dimensional Hubbard model

Robin Scholle, Walter Metzner, Demetrio Vilardi +1 more

We obtain an almost complete understanding of the mean-field phase diagram of the two-dimensional Hubbard model on a square lattice with a sizable next-nearest-neighbor hopping and a moderate interaction strength. In particular, we clarify the nature of the transition region between the spiral and the stripe phase. Complementing previous [] real-space Hartree-Fock calculations on large finite lattices, we solve the mean-field equations for coplanar unidirectional magnetic order directly in the thermodynamic limit, and we determine the nature of the magnetic states right below the mean-field critical temperature T * by a Landau free-energy analysis. While the magnetic order for filling factors n ≥ 1 is always of Néel type, for n ≤ 1 the following sequence of magnetic states is found as a function of increasing hole-doping: Néel, planar circular spiral, multispiral, and collinear spin-charge stripe states. Multispiral states are superpositions of several spirals with distinct wave vectors, and lead to concomitant charge order. We finally point out that nematic and charge orders inherited from the magnetic order can survive even in the presence of fluctuations, and we present a corresponding qualitative phase diagram. Published by the American Physical Society 2024

arXiv:2403.09862•Physical review. B./Physical review. B•Jun 27, 2024•12 citations

Coexisting magnetic, charge, and superconducting orders in the two-dimensional Hubbard model

Robin Scholle, Pietro M. Bonetti, Walter Metzner +1 more

We perform a renormalized mean-field study of the two-dimensional repulsive Hubbard model, focusing on the intricate interplay and possible coexistence of magnetic, charge, and superconducting orders. We improve on conventional mean-field theory by utilizing a renormalization group framework that captures high-energy fluctuations. This method generates effective magnetic and d -wave pairing interactions, and allows for an unbiased exploration of coexisting phases at weak and moderate interaction strengths. Unrestricted mean-field calculations of the effective Hamiltonian on large finite lattices are combined with analyses in the thermodynamic limit, revealing a rich phase diagram with extensive regions of coexisting orders. We find that d -wave superconductivity coexists with Néel order on the electron-doped side. On the hole-doped side, superconductivity is found to coexist with spiral or stripe magnetic orders. Within the stripe-ordered region, the superconducting order parameter is spatially modulated, with a period that follows the charge modulation of the stripes. Below Van Hove filling, pairing provides the primary energy gain, while the stripe order yields only a small, and hence fragile, additional energy lowering.

arXiv:2602.20073•Physical review. B./Physical review. B•Jul 28, 2026•0 citations

Coexistence of superconductivity with partially filled stripes in the Hubbard model

Hao Xu, Chia-Min Chung, Mingpu Qin +3 more

The Hubbard model is an iconic model in quantum many-body physics and has been intensely studied, especially since the discovery of high-temperature cuprate superconductors. Combining the complementary capabilities of two computational methods, we found superconductivity in both the electron- and hole-doped regimes of the two-dimensional Hubbard model with next-nearest-neighbor hopping. In the electron-doped regime, superconductivity was weaker and was accompanied by antiferromagnetic Néel correlations at low doping. The strong superconductivity on the hole-doped side coexisted with stripe order, which persisted into the overdoped region with weaker hole-density modulation. These stripe orders varied in fillings between 0.6 and 0.8. Our results suggest the applicability of the Hubbard model with next-nearest hopping for describing cuprate high–transition temperature ( T c ) superconductivity.

arXiv:2303.08376•Science•May 9, 2024•144 citations

Spin- and charge-density waves in the Hartree–Fock ground state of the two-dimensional Hubbard model

Jie Xu, Chia-Chen Chang, Eric J Walter +1 more

The ground states of the two-dimensional repulsive Hubbard model are studied within the unrestricted Hartree-Fock (UHF) theory. Magnetic and charge properties are determined by systematic, large-scale, exact numerical calculations, and quantified as a function of electron doping, h. In the solution of the self-consistent UHF equations, multiple initial configurations and simulated annealing are used to facilitate convergence to the global minimum. New approaches are employed to minimize finite-size effects in order to reach the thermodynamic limit. At low to moderate interacting strengths and low doping, the UHF ground state is a linear spin-density wave (l-SDW), with antiferromagnetic order and a modulating wave. The wavelength of the modulating wave is 2/h. Corresponding charge order exists but is substantially weaker than the spin order, hence holes are mobile. As the interaction is increased, the l-SDW states evolve into several different phases, with the holes eventually becoming localized. A simple pairing model is presented with analytic calculations for low interaction strength and small doping, to help understand the numerical results and provide a physical picture for the properties of the SDW ground state. By comparison with recent many-body calculations, it is shown that, for intermediate interactions, the UHF solution provides a good description of the magnetic correlations in the true ground state of the Hubbard model.

arXiv:1107.0976•Journal of Physics Condensed Matter•Nov 30, 2011•39 citations

Stripes and spin-density waves in the doped two-dimensional Hubbard model: Ground state phase diagram

Hao Xu, Hao Shi, Ettore Vitali +2 more

We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near $1/8$ doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about $1/5$ doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.

arXiv:2112.02187•Physical Review Research•Mar 28, 2022•74 citations
The two-dimensional Hubbard model is one of the most fundamental models in condensed matter physics. It captures many properties of cuprate high-temperature superconductors. Despite its apparent simplicity, solving the Hubbard model is a challenging task. In this work, we employ the Hartree-Fock approximation, as well as functional renormalization group methods, to construct phase diagrams of the Hubbard model and to shed light on the zoo of magnetic, charge ordered, and superconducting states. The Hartree-Fock approximation can self-consistently determine the expectation values of local charge and magnetic order parameters. Mean-field theories like this can lead to a deeper understanding of the competition, coexistence and transition of different orders. We start by employing the Hartree-Fock approximation in a completely unrestricted form on finite-size lattices and complement the results with calculations in the thermodynamic limit. We construct mean-field phase diagrams of the two-dimensional Hubbard model in the thermodynamic limit for different values of the next-nearest neighbor hopping as a function of filling and temperature. We find that the dominant magnetic orders at moderate coupling strength are Néel order around half-filling, followed by spiral order, and stripe order as the hole doping is increased. We also clarify the nature of the transition from magnetic spiral to stripe order. We then extend the Hartree-Fock scheme by a functional renormalization group flow, which generates new effective interactions and allows us to also access the superconducting phase of the Hubbard model. We present an improved phase diagram, featuring Néel, spiral, and stripe order coexisting with d-wave superconductivity at low temperatures. We analyze the coexistence and competition between these orders. Finally, we perform unrestricted Hartree-Fock calculations of a Hubbard model on the graphene lattice and the Bernal bilayer graphene lattice. We also construct magnetic phase diagrams for these geometries, which feature Néel order around half-filling and various stripe orders in a secondary magnetic dome.
We present an unbiased mean-field analysis of magnetic and charge orders in the two-dimensional Hubbard model on a square lattice, at both zero and finite temperatures. Unrestricted Hartree-Fock calculations on large finite lattices are complemented by solutions restricted to N\'eel and circular spiral order in the thermodynamic limit. The magnetic states are classified by a systematic scheme based on the dominant Fourier components of the spin texture. On finite lattices a whole zoo of ordering patterns appears. We show that many of these states are finite-size artifacts related to the limited choice of ordering wave vectors on a finite lattice. In the thermodynamic limit only three classes of states with a relatively simple structure survive: N\'eel, circular spiral, and stripe states. Stripes involve also charge order and can be unidirectional or bidirectional, with horizontal and/or vertical orientation. We present complete phase diagrams in the plane spanned by electron density and temperature, for a moderate Hubbard interaction and various choices of the next-nearest-neighbor hopping amplitude.
We obtain an almost complete understanding of the mean-field phase diagram of the two-dimensional Hubbard model on a square lattice with a sizable next-nearest-neighbor hopping and a moderate interaction strength. In particular, we clarify the nature of the transition region between the spiral and the stripe phase. Complementing previous [] real-space Hartree-Fock calculations on large finite lattices, we solve the mean-field equations for coplanar unidirectional magnetic order directly in the thermodynamic limit, and we determine the nature of the magnetic states right below the mean-field critical temperature T * by a Landau free-energy analysis. While the magnetic order for filling factors n ≥ 1 is always of Néel type, for n ≤ 1 the following sequence of magnetic states is found as a function of increasing hole-doping: Néel, planar circular spiral, multispiral, and collinear spin-charge stripe states. Multispiral states are superpositions of several spirals with distinct wave vectors, and lead to concomitant charge order. We finally point out that nematic and charge orders inherited from the magnetic order can survive even in the presence of fluctuations, and we present a corresponding qualitative phase diagram. Published by the American Physical Society 2024
We perform a renormalized mean-field study of the two-dimensional repulsive Hubbard model, focusing on the intricate interplay and possible coexistence of magnetic, charge, and superconducting orders. We improve on conventional mean-field theory by utilizing a renormalization group framework that captures high-energy fluctuations. This method generates effective magnetic and d -wave pairing interactions, and allows for an unbiased exploration of coexisting phases at weak and moderate interaction strengths. Unrestricted mean-field calculations of the effective Hamiltonian on large finite lattices are combined with analyses in the thermodynamic limit, revealing a rich phase diagram with extensive regions of coexisting orders. We find that d -wave superconductivity coexists with Néel order on the electron-doped side. On the hole-doped side, superconductivity is found to coexist with spiral or stripe magnetic orders. Within the stripe-ordered region, the superconducting order parameter is spatially modulated, with a period that follows the charge modulation of the stripes. Below Van Hove filling, pairing provides the primary energy gain, while the stripe order yields only a small, and hence fragile, additional energy lowering.
The Hubbard model is an iconic model in quantum many-body physics and has been intensely studied, especially since the discovery of high-temperature cuprate superconductors. Combining the complementary capabilities of two computational methods, we found superconductivity in both the electron- and hole-doped regimes of the two-dimensional Hubbard model with next-nearest-neighbor hopping. In the electron-doped regime, superconductivity was weaker and was accompanied by antiferromagnetic Néel correlations at low doping. The strong superconductivity on the hole-doped side coexisted with stripe order, which persisted into the overdoped region with weaker hole-density modulation. These stripe orders varied in fillings between 0.6 and 0.8. Our results suggest the applicability of the Hubbard model with next-nearest hopping for describing cuprate high–transition temperature ( T c ) superconductivity.
The ground states of the two-dimensional repulsive Hubbard model are studied within the unrestricted Hartree-Fock (UHF) theory. Magnetic and charge properties are determined by systematic, large-scale, exact numerical calculations, and quantified as a function of electron doping, h. In the solution of the self-consistent UHF equations, multiple initial configurations and simulated annealing are used to facilitate convergence to the global minimum. New approaches are employed to minimize finite-size effects in order to reach the thermodynamic limit. At low to moderate interacting strengths and low doping, the UHF ground state is a linear spin-density wave (l-SDW), with antiferromagnetic order and a modulating wave. The wavelength of the modulating wave is 2/h. Corresponding charge order exists but is substantially weaker than the spin order, hence holes are mobile. As the interaction is increased, the l-SDW states evolve into several different phases, with the holes eventually becoming localized. A simple pairing model is presented with analytic calculations for low interaction strength and small doping, to help understand the numerical results and provide a physical picture for the properties of the SDW ground state. By comparison with recent many-body calculations, it is shown that, for intermediate interactions, the UHF solution provides a good description of the magnetic correlations in the true ground state of the Hubbard model.
We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near 1/81/81/8 doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about 1/51/51/5 doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.