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:1107.0976 · Journal of Physics Condensed Matter

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, Shiwei Zhang

Published Nov 30, 2011• 39 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1088/0953-8984/23/50/505601

Related papers

Spin and Charge Order in the Doped Hubbard Model: Long-Wavelength Collective Modes

Chia‐Chen Chang, Shiwei Zhang

Determining the properties of the two-dimensional Hubbard model is an outstanding problem in physics. Applying recent advances in constrained path auxiliary-field quantum Monte Carlo techniques and simulating large rectangular supercells, we characterize the magnetic and charge properties in the ground state as a function of doping. At intermediate interaction strengths, an incommensurate spin density wave (SDW) state is found, with antiferromagnetic order and essentially homogeneous charge correlation. The wavelength of the collective mode decreases with doping, as does its magnitude. The SDW order vanishes beyond a critical doping. As the interaction is increased, the holes go from a wavelike (delocalized) to a particlelike (localized) state, and charge ordering develops which eventually evolves into stripelike states.

arXiv:0910.0712•Physical Review Letters•Mar 18, 2010•93 citations

Magnetic order in the Hubbard model in three dimensions and the crossover to two dimensions

Jie Xu, Simone Chiesa, Eric J. Walter +1 more

Systems of fermions described by the three-dimensional (3D) repulsive Hubbard model on a cubic lattice have recently attracted considerable attention due to their possible experimental realization via cold atoms in an optical lattice. Analytical and numerical results are limited away from half-filling. We study the ground state of the doped system from weak to intermediate interaction strengths within the generalized Hartree-Fock approximation. The exact solution to the self-consistent-field equations in the thermodynamic limit is obtained and the ground state is shown to exhibit antiferromagnetic order and incommensurate spin-density waves (SDW). At low interaction strengths, the SDW has unidirectional character with a leading wavevector along the [100]-direction, and the system is metallic. As the interaction increases, the system undergoes a simultaneous structural and metal-to-insulator transition to a unidirectional SDW state along the [111]-direction, with a different wavelength. We systematically determine the real and momentum space properties of these states. The crossover from 3D to two dimensions (2D) is then studied by varying the inter-plane hopping amplitude, which can be experimentally realized by tuning the distance between a stack of square-lattice layers. Detailed comparisons are made between the exact numerical results and predictions from the pairing model, a variational ansatz based on the pairing of spins in the vicinity of the Fermi surface. Most of the numerical results can be understood quantitatively from the ansatz, which provides a simple picture for the nature of the SDW states.

arXiv:1301.6186•Journal of Physics Condensed Matter•Sep 18, 2013•4 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 zoo of states in the Hubbard model

Robin Scholle

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.

•Jan 1, 2026•0 citations

Group-theoretical and resonating Hartree-Fock studies for spin fluctuations in two-dimensional Hubbard model on triangular lattice

Fumikazu Satoh, Masa-aki Ozaki, Takaaki Maruyama +1 more

Magnetic structures of the two-dimensional Hubbard model on a uniform triangular lattice are investigated from the viewpoint of quantum fluctuations. First, possible spin density wave (SDW) solutions in the unrestricted Hartree-Fock approximation are group theoretically classified for the ordering vectors of ${\mathbf{Q}}_{1}=({\mathbf{G}}_{a}+{\mathbf{G}}_{b})/4$, ${\mathbf{Q}}_{2}={C}_{3}^{+}{\mathbf{Q}}_{1}$, and ${\mathbf{Q}}_{3}={C}_{3}^{\ensuremath{-}}{\mathbf{Q}}_{1}$. Then, to include the electron correlation effects efficiently, we apply the resonating Hartree-Fock method. It is shown that the spin fluctuations in the metallic ground state at $U/t\ensuremath{\sim}4.5$ are reasonably described as hybridization of two different SDW states obtained by the group theoretical classification and their low-energy excited states. These fluctuations make the long-range magnetic order significantly reduced in the correlated ground state.

Physical Review B•Dec 7, 2011•4 citations

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
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.
Determining the properties of the two-dimensional Hubbard model is an outstanding problem in physics. Applying recent advances in constrained path auxiliary-field quantum Monte Carlo techniques and simulating large rectangular supercells, we characterize the magnetic and charge properties in the ground state as a function of doping. At intermediate interaction strengths, an incommensurate spin density wave (SDW) state is found, with antiferromagnetic order and essentially homogeneous charge correlation. The wavelength of the collective mode decreases with doping, as does its magnitude. The SDW order vanishes beyond a critical doping. As the interaction is increased, the holes go from a wavelike (delocalized) to a particlelike (localized) state, and charge ordering develops which eventually evolves into stripelike states.
Systems of fermions described by the three-dimensional (3D) repulsive Hubbard model on a cubic lattice have recently attracted considerable attention due to their possible experimental realization via cold atoms in an optical lattice. Analytical and numerical results are limited away from half-filling. We study the ground state of the doped system from weak to intermediate interaction strengths within the generalized Hartree-Fock approximation. The exact solution to the self-consistent-field equations in the thermodynamic limit is obtained and the ground state is shown to exhibit antiferromagnetic order and incommensurate spin-density waves (SDW). At low interaction strengths, the SDW has unidirectional character with a leading wavevector along the [100]-direction, and the system is metallic. As the interaction increases, the system undergoes a simultaneous structural and metal-to-insulator transition to a unidirectional SDW state along the [111]-direction, with a different wavelength. We systematically determine the real and momentum space properties of these states. The crossover from 3D to two dimensions (2D) is then studied by varying the inter-plane hopping amplitude, which can be experimentally realized by tuning the distance between a stack of square-lattice layers. Detailed comparisons are made between the exact numerical results and predictions from the pairing model, a variational ansatz based on the pairing of spins in the vicinity of the Fermi surface. Most of the numerical results can be understood quantitatively from the ansatz, which provides a simple picture for the nature of the SDW states.
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.
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.
Magnetic structures of the two-dimensional Hubbard model on a uniform triangular lattice are investigated from the viewpoint of quantum fluctuations. First, possible spin density wave (SDW) solutions in the unrestricted Hartree-Fock approximation are group theoretically classified for the ordering vectors of Q1=(Ga+Gb)/4{\mathbf{Q}}_{1}=({\mathbf{G}}_{a}+{\mathbf{G}}_{b})/4Q1​=(Ga​+Gb​)/4, Q2=C3+Q1{\mathbf{Q}}_{2}={C}_{3}^{+}{\mathbf{Q}}_{1}Q2​=C3+​Q1​, and ${\mathbf{Q}}_{3}={C}_{3}^{\ensuremath{-}}{\mathbf{Q}}_{1}$. Then, to include the electron correlation effects efficiently, we apply the resonating Hartree-Fock method. It is shown that the spin fluctuations in the metallic ground state at $U/t\ensuremath{\sim}4.5$ are reasonably described as hybridization of two different SDW states obtained by the group theoretical classification and their low-energy excited states. These fluctuations make the long-range magnetic order significantly reduced in the correlated ground state.
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.