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:0805.4831 · Physical Review B

Spatially inhomogeneous phase in the two-dimensional repulsive Hubbard model

Chia‐Chen Chang, Shiwei Zhang

Published Oct 2, 2008• 88 citations
arXivHTMLPDF

Abstract

Using recent advances in auxiliary-field quantum Monte Carlo techniques and the phaseless approximation to control the sign/phase problem, we determine the equation of state in the ground state of the two-dimensional repulsive single-band Hubbard model at intermediate interactions. Shell effects are eliminated and finite-size effects are greatly reduced by boundary-condition integration. Spin-spin correlation functions and structure factors are also calculated. In lattice sizes up to $16\ifmmode\times\else\texttimes\fi{}16$, the results show signals for phase separation. Upon doping, the system separates into one phase of density $n=1$ (hole free) and the other at density ${n}_{c}$ $(\ensuremath{\sim}0.9)$. The long-range antiferromagnetic order is coupled to this process and is lost below ${n}_{c}$.

DOI
10.1103/physrevb.78.165101

Related papers

Spatially homogeneous ground state of the two-dimensional Hubbard model

Federico Becca, Massimo Capone, Sandro Sorella

We investigate the stability with respect to phase separation or charge-density-wave formation of the two-dimensional Hubbard model for various values of the local Coulomb repulsion and electron densities using Green-function Monte Carlo techniques. The well known sign problem is particularly serious in the relevant region of small hole doping. We show that the difference in accuracy for different doping makes it very difficult to probe the phase separation instability using only energy calculations, even in the weak-coupling limit $(U=4t)$ where reliable results are available. By contrast, the knowledge of the charge correlation functions allows us to provide clear evidence of a spatially homogeneous ground state up to $U=10t.$

arXiv:cond-mat/0010165•Physical review. B, Condensed matter•Nov 15, 2000•56 citations

Monte Carlo Study of the Two-Dimensional Hubbard Model

J. E. Hirsch

A two-dimensional model of interacting electrons on a lattice, the Hubbard model, has been studied by means of a numerical simulation technique. Results for the half-filled-band sector are presented for lattices of up to 6\ifmmode\times\else\texttimes\fi{}6 spatial sites, as a function of temperature and electron-electron repulsion $U$. In the ground state, the system is found to exhibit antiferromagnetic long-range order for all values of $U$, although greatly reduced from the mean-field-theory predictions.

Physical Review Letters•Nov 14, 1983•77 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

Phase separation in the two-dimensional Hubbard model: A fixed-node quantum Monte Carlo study

Andrea C. Cosentini, Massimo Capone, Leonardo Guidoni +1 more

Fixed-node Green's-function Monte Carlo calculations have been performed for very large 16\ifmmode\times\else\texttimes\fi{}16 two-dimensional Hubbard lattices, large interaction strengths $U=10,20,$ and $40$, and many $(15\ensuremath{\sim}20)$ densities between empty and half-filling. The nodes were fixed by a simple Slater-Gutzwiller trial wave function. For each value of U we obtained a sequence of ground-state energies which is consistent with the possibility of a phase separation close to half-filling, with a hole density in the hole-rich phase which is a decreasing function of $U$. The energies suffer, however, from a fixed-node bias: more accurate nodes are needed to confirm this picture. Our extensive numerical results and their test against size, shell, shape, and boundary-condition effects also suggest that phase separation is quite a delicate issue, on which simulations based on smaller lattices than considered here are unlikely to give reliable predictions.

arXiv:cond-mat/9801299•Physical review. B, Condensed matter•Dec 1, 1998•35 citations

Magnetism in the two-dimensional Hubbard model (abstract)

J. E. Hirsch, Hsueh‐Ju Lin

The Hubbard model has been widely used to describe magnetism due to itinerant electrons. However, except for the one-dimensional case, no exact solution of the model exists. We have studied the two-dimensional Hubbard model on a square lattice using a Monte Carlo simulation technique. Simulations were performed on lattices of up to 8×8 in spatial size. We will discuss results for thermodynamic properties and correlation functions as a function of temperature and interaction strength. For the half-filled band case, our results indicate that the system exhibits long-range antiferromagnetic order in the ground state for any nonzero interaction, although substantially reduced from the classical Néel state due to both charge and spin fluctuations. For the non-half-filled band case, our results suggest that no long-range magnetic order exists. We also don’t find any indication of short-ranged ferromagnetic correlations in the parameter range studied. We have also studied the predictions of various approximate treatments for this model, particularly the static approximation and the coherent potential approximation (CPA). By comparison with the Monte Carlo results, we discuss the validity of these schemes. Finally, we discuss the effect of changing the band structure on the magnetic properties of the system.

Journal of Applied Physics•Apr 15, 1985•0 citations
Using recent advances in auxiliary-field quantum Monte Carlo techniques and the phaseless approximation to control the sign/phase problem, we determine the equation of state in the ground state of the two-dimensional repulsive single-band Hubbard model at intermediate interactions. Shell effects are eliminated and finite-size effects are greatly reduced by boundary-condition integration. Spin-spin correlation functions and structure factors are also calculated. In lattice sizes up to 16\ifmmode×\else\texttimes\fi1616\ifmmode\times\else\texttimes\fi{}1616\ifmmode×\else\texttimes\fi16, the results show signals for phase separation. Upon doping, the system separates into one phase of density n=1n=1n=1 (hole free) and the other at density nc{n}_{c}nc​ (\ensuremath∼0.9)(\ensuremath{\sim}0.9)(\ensuremath∼0.9). The long-range antiferromagnetic order is coupled to this process and is lost below nc{n}_{c}nc​.
We investigate the stability with respect to phase separation or charge-density-wave formation of the two-dimensional Hubbard model for various values of the local Coulomb repulsion and electron densities using Green-function Monte Carlo techniques. The well known sign problem is particularly serious in the relevant region of small hole doping. We show that the difference in accuracy for different doping makes it very difficult to probe the phase separation instability using only energy calculations, even in the weak-coupling limit (U=4t)(U=4t)(U=4t) where reliable results are available. By contrast, the knowledge of the charge correlation functions allows us to provide clear evidence of a spatially homogeneous ground state up to U=10t.U=10t.U=10t.
A two-dimensional model of interacting electrons on a lattice, the Hubbard model, has been studied by means of a numerical simulation technique. Results for the half-filled-band sector are presented for lattices of up to 6\ifmmode\times\else\texttimes\fi{}6 spatial sites, as a function of temperature and electron-electron repulsion UUU. In the ground state, the system is found to exhibit antiferromagnetic long-range order for all values of UUU, although greatly reduced from the mean-field-theory predictions.
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.
Fixed-node Green's-function Monte Carlo calculations have been performed for very large 16\ifmmode\times\else\texttimes\fi{}16 two-dimensional Hubbard lattices, large interaction strengths U=10,20,U=10,20,U=10,20, and 404040, and many $(15\ensuremath{\sim}20)$ densities between empty and half-filling. The nodes were fixed by a simple Slater-Gutzwiller trial wave function. For each value of U we obtained a sequence of ground-state energies which is consistent with the possibility of a phase separation close to half-filling, with a hole density in the hole-rich phase which is a decreasing function of UUU. The energies suffer, however, from a fixed-node bias: more accurate nodes are needed to confirm this picture. Our extensive numerical results and their test against size, shell, shape, and boundary-condition effects also suggest that phase separation is quite a delicate issue, on which simulations based on smaller lattices than considered here are unlikely to give reliable predictions.
The Hubbard model has been widely used to describe magnetism due to itinerant electrons. However, except for the one-dimensional case, no exact solution of the model exists. We have studied the two-dimensional Hubbard model on a square lattice using a Monte Carlo simulation technique. Simulations were performed on lattices of up to 8×8 in spatial size. We will discuss results for thermodynamic properties and correlation functions as a function of temperature and interaction strength. For the half-filled band case, our results indicate that the system exhibits long-range antiferromagnetic order in the ground state for any nonzero interaction, although substantially reduced from the classical Néel state due to both charge and spin fluctuations. For the non-half-filled band case, our results suggest that no long-range magnetic order exists. We also don’t find any indication of short-ranged ferromagnetic correlations in the parameter range studied. We have also studied the predictions of various approximate treatments for this model, particularly the static approximation and the coherent potential approximation (CPA). By comparison with the Monte Carlo results, we discuss the validity of these schemes. Finally, we discuss the effect of changing the band structure on the magnetic properties of the system.