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

Physical Review B

Insights into the structure of many-electron wave functions of Mott-insulating antiferromagnets: The three-band Hubbard model in full configuration interaction quantum Monte Carlo

Lauretta R. Schwarz, George H. Booth, Ali Alavi

Published Jan 30, 2015• 16 citations
Publisher

Abstract

We investigate the ground-state wave function of a prototypical strongly correlated system, a three-band $(p\ensuremath{-}d)$ Hubbard model of cuprates, using full configuration interaction quantum Monte Carlo. We show that the configuration interaction description of the exact ground state wave function is profoundly affected by the choice of single-particle representation, in a counterintuitive manner. Thus a broken-symmetry unrestricted Hartree-Fock basis, which at a single configuration level produces a qualitatively correct description of the antiferromagnet, results in a highly entangled exact solution consisting of high particle-hole excitations of the reference. This wave function is found to be very difficult to approximate using subspace diagonalizations. Conversely, a restricted Hartree-Fock basis, which yields at a single configuration level a qualitatively incorrect paramagnetic metal, results in a relatively rapidly converging configuration interaction expansion. Convergence can be further accelerated by adopting a natural orbital representation. Our results suggest that with the correct single-particle basis, such strongly correlated systems may be described by relatively compact wave functions.

DOI
10.1103/physrevb.91.045139

Related papers

Single-Particle Excitations under Coexisting Electron Correlation and Disorder: A Numerical Study of the Anderson–Hubbard Model

Hiroshi Shinaoka, Masatoshi Imada

Interplay of electron correlation and randomness is studied by using the Anderson-Hubbard model within the Hartree-Fock approximation. Under the coexistence of short-range interaction and diagonal disorder, we obtain the ground-state phase diagram in three dimensions, which includes an antiferromagnetic insulator, an antiferromagnetic metal, a paramagnetic insulator (Anderson-localized insulator) and a paramagnetic metal. Although only the short-range interaction is present in this model, we find unconventional soft gaps in the insulating phases irrespective of electron filling, spatial dimensions and long-range order, where the single-particle density of states (DOS) vanishes with a power-law scaling in one dimension (1D) or even faster in two dimensions (2D) and three dimensions (3D) toward the Fermi energy. We call it soft Hubbard gap. Moreover, exact-diagonalization results in 1D support the formation of the soft Hubbard gap beyond the mean-field level. The formation of the soft Hubbard gap cannot be attributed to a conventional theory by Efros and Shklovskii (ES) owing the emergence of soft gaps to the long-range Coulomb interaction. Indeed, based on a picture of multivalley energy landscape, we propose a phenomenological scaling theory, which predicts a scaling of the DOS in perfect agreement with the numerical results. We further discuss a correction of the scaling of the DOS by the long-range part of the Coulomb interaction, which modifies the scaling of Efros and Shklovskii. Furthermore, explicit formulae for the temperature dependence of the DC resistivity via variable-range hopping under the influence of the soft gaps are derived. Finally, we compare the present theory with experimental results of SrRu_{1-x}Ti_xO_3.

arXiv:0905.3248•Journal of the Physical Society of Japan•Aug 28, 2009•53 citations

Quantum Monte Carlo Simulations and Weak-Coupling Approximations for the Three-Band Hubbard Model

R. Putz, H. Endres, A. Muramatsu +1 more

We study the three-band Hubbard model for the cuprate superconductors using Quantum Monte Carlo simulations and the fluctuation exchange approximation including a self-consistent many-body renormalization of the single-particle propagator. The energy dispersion of the low-lying one-particle excitations is in very good agreement with spectroscopic measurements for a hole doping concentration δ = 0.25 of the CuO2 planes. At the same doping, we examine the spin susceptibility yielding a non-ordered phase with short-range incommensurate fluctuations near the antiferromagnetic wave vector. In the superconducting regime, the interaction vertex for the pairing correlation function in the extended s-wave channel shows a maximum for δ = 0.20, resembling the dependence of transition temperature Tc on doping.

Notes on numerical fluid mechanics•Jan 1, 1994•1 citation

Quantum Monte Carlo study of a single hole in a quantum antiferromagnet

Sandro Sorella

Using the standard quantum Monte Carlo technique for the Hubbard model, I present here a numerical investigation of the hole propagation in a quantum antiferromagnet. The calculation is very well stabilized, using selected-size systems and a special choice of the trial wave function that satisfies the ``closed-shell condition'' in the presence of an arbitrarily weak Zeeman magnetic field. First it is shown that the presence of this magnetic field does not affect the thermodynamic properties of the half-filled system. Then I use the same selected sizes for the one-hole ground state. I investigate the question of vanishing or nonvanishing quasiparticle weight, in order to clarify whether the Mott insulator behaves as a conventional insulator with an upper and lower Hubbard band. By comparing the present finite-size scaling with several techniques that predict a finite quasiparticle weight, the data seem more consistent with a vanishing quasiparticle weight, i.e., the Mott-Hubbard insulator should be characterized by nontrivial excitations that cannot be interpreted in a simple picture. However, one cannot exclude, based only on numerical grounds, a very small but nonvanishing quasiparticle weight. In the last part I give some theoretical arguments to explain the results of the Monte Carlo simulation.

Physical review. B, Condensed matter•Nov 1, 1992•29 citations

Anomalous low-doping phase of the Hubbard model

C. Gröber, R. Eder, Werner R. Hanke

We present results of a systematic quantum Monte Carlo study for the single-band Hubbard model. Thereby we evaluated single-particle spectra (PES and IPES), two-particle spectra (spin and density correlation functions), and the dynamical correlation function of suitably defined diagnostic operators, all as a function of temperature and hole doping. The results allow us to identify different physical regimes. Near half-filling we find an anomalous ``Hubbard-I phase,'' where the band structure is, up to some minor modifications, consistent with the Hubbard-I predictions. At lower temperatures, where the spin response becomes sharp, additional dispersionless ``bands'' emerge due to the dressing of electrons/holes with spin excitations. We present a simple phenomenological fit that reproduces the band structure of the insulator quantitatively. The Fermi surface volume in the low-doping phase, as derived from the single-particle spectral function, is not consistent with the Luttinger theorem, but qualitatively in agreement with the predictions of the Hubbard-I approximation. The anomalous phase extends up to a hole concentration of $\ensuremath{\approx}15%,$ i.e., the underdoped region in the phase diagram of high-${T}_{c}$ superconductors. We also investigate the nature of the magnetic ordering transition in the single-particle spectra. We show that the transition to a spin-density wave-like band structure is not accomplished by the formation of any resolvable ``precursor bands,'' but rather by a (spectroscopically invisible) band of spin-3/2 quasiparticles. We discuss implications for the ``remnant Fermi surface'' in insulating cuprate compounds and the shadow bands in the doped materials.

arXiv:cond-mat/0001366•Physical review. B, Condensed matter•Aug 15, 2000•137 citations

The Many Guises of Quasiparticles in Correlated Electron Systems

Ivan Pasqua

Many strongly correlated electron systems fall outside the conventional quasiparticle paradigm. Their low-energy fermionic excitations need not be adiabatically connected to the bare electrons or appear as poles of the physical electron single-particle Green's function. This thesis investigates to what extent an effective quasiparticle description can nevertheless survive, and which emergent degrees of freedom provide the appropriate low-energy structure. First, the interaction-driven evolution from a quantum spin Hall insulator to a Mott insulator is studied using the dynamical cluster approximation. Dispersive zeros of the single-particle Green's function emerge near the Mott transition and retain a nontrivial topological character. When symmetry breaking is allowed, a nontopological excitonic insulator is found to intrude between quantum spin Hall and Mott insulators. Encoding both poles and zeros within a low-energy quasiparticle Hamiltonian provides a continuous one-body description of the three phases. On the Mott side, its gap evolution supports interpreting the soft excitonic mode as a bound state of quasiparticle and quasihole excitations tied to the Green's function zeros, rather than to the widely separated Hubbard bands. Second, the intrinsic anomalous Hall conductivity of a topological metal is analyzed within a multiband extension of Landau Fermi-liquid theory. Residual interactions among quasiparticles at the Fermi surface dress the dynamic current vertex and generate corrections to the Berry-curvature contribution obtained from the quasiparticle bands alone. Our result supports recent claims that the correct expressions for topological observables include vertex corrections besides the topological invariants built just upon the Green's function. It also demonstrates that such corrections are naturally accounted for by Landau Fermi-liquid theory. Extracting the relevant ingredients for an effective description of material-relevant correlated systems requires computational tools that remain reliable at low temperatures. To this end, a deterministic strong-coupling impurity solver for dynamical mean-field theory is discussed. By providing static, dynamical, and thermodynamic observables without stochastic sampling, it offers a promising route toward multiorbital and spin-orbit coupled problems where conventional Monte Carlo methods currently struggle. Finally, the auxiliary quasiparticle framework underlying the ghost-Gutzwiller Ansatz is investigated. On a model for correlated quantum spin Hall insulators, we show that the topology of the interacting state is encoded in the auxiliary one-body Hamiltonian. The resulting band-structure reveals topological Hubbard bands whose topological character can be tuned through a finite magnetization. When applied to the single-band t-J model, the same framework yields a fractionalized Fermi liquid (FL*). In this phase, dispersive neutral spinons, with vanishing physical-electron spectral weight, coexist with conventional quasiparticles forming a small hole-like Fermi surface, thus violating Luttinger's theorem. The corresponding temperature--doping phase diagram contains a low-doping FL*, a d-wave superconducting dome, and an overdoped conventional Fermi liquid, thereby reproducing key qualitative features of cuprate phenomenology.

•Sep 14, 2026•0 citations

Optimized wave function by kinetic renormalization effect in strongly correlated region of the three-band Hubbard model

Takashi Yanagisawa, Izumi Hase, Mitake Miyazaki +1 more

We investigate the ground state of the two-dimensional single-band Hubbard model and three-band Hubbard model (d-p model) by using a variational Monte Carlo method. We employ an optimization variational Monte Carlo method. In this method the many-body wave function is improved and optimized by introducing variational parameters in the form: ψ = exp(- S ) ψ 0 where ψ 0 is a wave function, for example, the Gutzwiller function, and S is an operator used to include the many-body effect. The strong magnetic correlation is induced by the on-site Coulomb repulsive interaction U . The magnetic correlation plays a significant role in correlated electron systems as an origin of electron pairing as well as a main competitor of superconductivity. High-temperature superconductivity occurs in the strongly correlated region where U is greater than the bandwidth. We investigate the stability of antiferromagnetic state in the three-band d-p model.

Journal of Physics Conference Series•Jul 1, 2018•2 citations
We investigate the ground-state wave function of a prototypical strongly correlated system, a three-band (p\ensuremath−d)(p\ensuremath{-}d)(p\ensuremath−d) Hubbard model of cuprates, using full configuration interaction quantum Monte Carlo. We show that the configuration interaction description of the exact ground state wave function is profoundly affected by the choice of single-particle representation, in a counterintuitive manner. Thus a broken-symmetry unrestricted Hartree-Fock basis, which at a single configuration level produces a qualitatively correct description of the antiferromagnet, results in a highly entangled exact solution consisting of high particle-hole excitations of the reference. This wave function is found to be very difficult to approximate using subspace diagonalizations. Conversely, a restricted Hartree-Fock basis, which yields at a single configuration level a qualitatively incorrect paramagnetic metal, results in a relatively rapidly converging configuration interaction expansion. Convergence can be further accelerated by adopting a natural orbital representation. Our results suggest that with the correct single-particle basis, such strongly correlated systems may be described by relatively compact wave functions.
Interplay of electron correlation and randomness is studied by using the Anderson-Hubbard model within the Hartree-Fock approximation. Under the coexistence of short-range interaction and diagonal disorder, we obtain the ground-state phase diagram in three dimensions, which includes an antiferromagnetic insulator, an antiferromagnetic metal, a paramagnetic insulator (Anderson-localized insulator) and a paramagnetic metal. Although only the short-range interaction is present in this model, we find unconventional soft gaps in the insulating phases irrespective of electron filling, spatial dimensions and long-range order, where the single-particle density of states (DOS) vanishes with a power-law scaling in one dimension (1D) or even faster in two dimensions (2D) and three dimensions (3D) toward the Fermi energy. We call it soft Hubbard gap. Moreover, exact-diagonalization results in 1D support the formation of the soft Hubbard gap beyond the mean-field level. The formation of the soft Hubbard gap cannot be attributed to a conventional theory by Efros and Shklovskii (ES) owing the emergence of soft gaps to the long-range Coulomb interaction. Indeed, based on a picture of multivalley energy landscape, we propose a phenomenological scaling theory, which predicts a scaling of the DOS in perfect agreement with the numerical results. We further discuss a correction of the scaling of the DOS by the long-range part of the Coulomb interaction, which modifies the scaling of Efros and Shklovskii. Furthermore, explicit formulae for the temperature dependence of the DC resistivity via variable-range hopping under the influence of the soft gaps are derived. Finally, we compare the present theory with experimental results of SrRu_{1-x}Ti_xO_3.
We study the three-band Hubbard model for the cuprate superconductors using Quantum Monte Carlo simulations and the fluctuation exchange approximation including a self-consistent many-body renormalization of the single-particle propagator. The energy dispersion of the low-lying one-particle excitations is in very good agreement with spectroscopic measurements for a hole doping concentration δ = 0.25 of the CuO2 planes. At the same doping, we examine the spin susceptibility yielding a non-ordered phase with short-range incommensurate fluctuations near the antiferromagnetic wave vector. In the superconducting regime, the interaction vertex for the pairing correlation function in the extended s-wave channel shows a maximum for δ = 0.20, resembling the dependence of transition temperature Tc on doping.
Using the standard quantum Monte Carlo technique for the Hubbard model, I present here a numerical investigation of the hole propagation in a quantum antiferromagnet. The calculation is very well stabilized, using selected-size systems and a special choice of the trial wave function that satisfies the ``closed-shell condition'' in the presence of an arbitrarily weak Zeeman magnetic field. First it is shown that the presence of this magnetic field does not affect the thermodynamic properties of the half-filled system. Then I use the same selected sizes for the one-hole ground state. I investigate the question of vanishing or nonvanishing quasiparticle weight, in order to clarify whether the Mott insulator behaves as a conventional insulator with an upper and lower Hubbard band. By comparing the present finite-size scaling with several techniques that predict a finite quasiparticle weight, the data seem more consistent with a vanishing quasiparticle weight, i.e., the Mott-Hubbard insulator should be characterized by nontrivial excitations that cannot be interpreted in a simple picture. However, one cannot exclude, based only on numerical grounds, a very small but nonvanishing quasiparticle weight. In the last part I give some theoretical arguments to explain the results of the Monte Carlo simulation.
We present results of a systematic quantum Monte Carlo study for the single-band Hubbard model. Thereby we evaluated single-particle spectra (PES and IPES), two-particle spectra (spin and density correlation functions), and the dynamical correlation function of suitably defined diagnostic operators, all as a function of temperature and hole doping. The results allow us to identify different physical regimes. Near half-filling we find an anomalous ``Hubbard-I phase,'' where the band structure is, up to some minor modifications, consistent with the Hubbard-I predictions. At lower temperatures, where the spin response becomes sharp, additional dispersionless ``bands'' emerge due to the dressing of electrons/holes with spin excitations. We present a simple phenomenological fit that reproduces the band structure of the insulator quantitatively. The Fermi surface volume in the low-doping phase, as derived from the single-particle spectral function, is not consistent with the Luttinger theorem, but qualitatively in agreement with the predictions of the Hubbard-I approximation. The anomalous phase extends up to a hole concentration of $\ensuremath{\approx}15%,$ i.e., the underdoped region in the phase diagram of high-Tc{T}_{c}Tc​ superconductors. We also investigate the nature of the magnetic ordering transition in the single-particle spectra. We show that the transition to a spin-density wave-like band structure is not accomplished by the formation of any resolvable ``precursor bands,'' but rather by a (spectroscopically invisible) band of spin-3/2 quasiparticles. We discuss implications for the ``remnant Fermi surface'' in insulating cuprate compounds and the shadow bands in the doped materials.
Many strongly correlated electron systems fall outside the conventional quasiparticle paradigm. Their low-energy fermionic excitations need not be adiabatically connected to the bare electrons or appear as poles of the physical electron single-particle Green's function. This thesis investigates to what extent an effective quasiparticle description can nevertheless survive, and which emergent degrees of freedom provide the appropriate low-energy structure. First, the interaction-driven evolution from a quantum spin Hall insulator to a Mott insulator is studied using the dynamical cluster approximation. Dispersive zeros of the single-particle Green's function emerge near the Mott transition and retain a nontrivial topological character. When symmetry breaking is allowed, a nontopological excitonic insulator is found to intrude between quantum spin Hall and Mott insulators. Encoding both poles and zeros within a low-energy quasiparticle Hamiltonian provides a continuous one-body description of the three phases. On the Mott side, its gap evolution supports interpreting the soft excitonic mode as a bound state of quasiparticle and quasihole excitations tied to the Green's function zeros, rather than to the widely separated Hubbard bands. Second, the intrinsic anomalous Hall conductivity of a topological metal is analyzed within a multiband extension of Landau Fermi-liquid theory. Residual interactions among quasiparticles at the Fermi surface dress the dynamic current vertex and generate corrections to the Berry-curvature contribution obtained from the quasiparticle bands alone. Our result supports recent claims that the correct expressions for topological observables include vertex corrections besides the topological invariants built just upon the Green's function. It also demonstrates that such corrections are naturally accounted for by Landau Fermi-liquid theory. Extracting the relevant ingredients for an effective description of material-relevant correlated systems requires computational tools that remain reliable at low temperatures. To this end, a deterministic strong-coupling impurity solver for dynamical mean-field theory is discussed. By providing static, dynamical, and thermodynamic observables without stochastic sampling, it offers a promising route toward multiorbital and spin-orbit coupled problems where conventional Monte Carlo methods currently struggle. Finally, the auxiliary quasiparticle framework underlying the ghost-Gutzwiller Ansatz is investigated. On a model for correlated quantum spin Hall insulators, we show that the topology of the interacting state is encoded in the auxiliary one-body Hamiltonian. The resulting band-structure reveals topological Hubbard bands whose topological character can be tuned through a finite magnetization. When applied to the single-band t-J model, the same framework yields a fractionalized Fermi liquid (FL*). In this phase, dispersive neutral spinons, with vanishing physical-electron spectral weight, coexist with conventional quasiparticles forming a small hole-like Fermi surface, thus violating Luttinger's theorem. The corresponding temperature--doping phase diagram contains a low-doping FL*, a d-wave superconducting dome, and an overdoped conventional Fermi liquid, thereby reproducing key qualitative features of cuprate phenomenology.
We investigate the ground state of the two-dimensional single-band Hubbard model and three-band Hubbard model (d-p model) by using a variational Monte Carlo method. We employ an optimization variational Monte Carlo method. In this method the many-body wave function is improved and optimized by introducing variational parameters in the form: ψ = exp(- S ) ψ 0 where ψ 0 is a wave function, for example, the Gutzwiller function, and S is an operator used to include the many-body effect. The strong magnetic correlation is induced by the on-site Coulomb repulsive interaction U . The magnetic correlation plays a significant role in correlated electron systems as an origin of electron pairing as well as a main competitor of superconductivity. High-temperature superconductivity occurs in the strongly correlated region where U is greater than the bandwidth. We investigate the stability of antiferromagnetic state in the three-band d-p model.