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./Physical review. B

Hole clustering and mutual interplay in three-band Hubbard model

Yan Peng, Mi Jiang, Yi‐feng Yang, Guang-Ming Zhang

Published Sep 19, 2025• 1 citation
Publisher

Abstract

Recent scanning tunneling spectroscopy (STS) experiments revealed a remarkable role of a supercell consisting of $4\ifmmode\times\else\texttimes\fi{}4 {\mathrm{CuO}}_{2}$ unit cells in the emergence of local nematic state and preformed local Cooper pairs and phase coherent cuprate superconductivity. By employing the numerically exact determinant Quantum Monte Carlo simulations, we mimic the effects of an experimental Ca vacancy by an external local potential to investigate the charge and spectral properties of the system hosting two-doped holes. The model numerically supports the role of the $4\ifmmode\times\else\texttimes\fi{}4$ supercell as the building block of hole-doped cuprates via the hole density distribution and local spectra around the local potential. Our results might provide a theoretical support for the experimental observations and a platform for investigating local charge order and local Cooper pairs on the $4\ifmmode\times\else\texttimes\fi{}4$ supercell as the plausible route to understanding unconventional cuprate superconductivity.

DOI
10.1103/b9sq-gyf4

Related papers

Hole clustering and mutual interplay in three-band Hubbard model

Mi Jiang, Yi‐feng Yang, Guang-Ming Zhang

Recent scanning tunnelling spectroscopy (STS) experiments revealed remarkable role of a supercell consisting $4\times4$ CuO$_2$ unit cells in the emergence of local nematic state and preformed local Cooper pairs and phase coherent cuprate superconductivity. By employing the numerically exact determinant Quantum Monte Carlo simulations, we mimic the effects of experimental Ca vacancy by an external local potential to investigate the charge and spectral properties of the system hosting two doped holes. The model numerically support the role of the $4\times4$ supercell as the building block of hole doped cuprates via the hole density distribution and local spectra around the local potential. Our results might provide a theoretical support on the experimental observations and a platform for investigating local charge order and local Cooper pairs on the $4\times4$ supercell as the plausible route to understanding unconventional cuprate superconductivity.

arXiv:2506.00787•Jun 1, 2025•0 citations

Nature of the effective interaction in electron-doped cuprate superconductors: A sign-problem-free quantum Monte Carlo study

Zixiang Li, Fa Wang, Hong Yao +1 more

Understanding the mechanism of Cooper pairing amounts to determining the effective interaction that operates at low energies. Efforts to achieve such a goal for superconducting materials, especially strongly correlated ones, from both bottom-up and top-down approaches, have been plagued by having to use uncontrolled approximations. Here, we perform large-scale, numerically exact, sign-problem-free zero-temperature quantum Monte Carlo simulations on an effective theory based on ``hot spots'' plus fluctuating collective modes. Because hot spots are clearly identified by angle-resolved photoemission spectroscopy for electron-doped cuprates, we focus our attention on such materials. Our goal is to determine the minimum effective action that can describe the observed superconductivity and charge-density wave. The results suggest that antiferromagnetic fluctuation alone is not sufficient---the effective action needs to be amended with nematic fluctuations. We believe that our results address the pairing mechanism of high-${T}_{c}$ superconductivity in electron-doped cuprates, and they shed light on the pairing mechanism of hole-doped cuprates.

arXiv:1512.04541•Physical review. B./Physical review. B•Jun 5, 2017•41 citations

Magnetic and charge orders in the ground state of the Emery model: Accurate numerical results

Adam Chiciak, Ettore Vitali, Shiwei Zhang

We perform extensive auxiliary-field quantum Monte Carlo (AFQMC) calculations for the three-band Hubbard (Emery) model in order to study the ground-state properties of copper-oxygen planes in the cuprates. Employing cutting-edge AFQMC techniques with a self-consistent gauge constraint in auxiliary-field space to control the sign problem, we reach supercells containing $\ensuremath{\approx}500$ atoms to capture collective modes in the charge and spin orders and characterize the behavior in the thermodynamic limit. The self-consistency scheme interfacing with generalized Hartree-Fock calculations allows high accuracy in AFQMC to resolve small energy scales, which is crucial for determining the complex candidate orders in such a system. We present detailed information on the charge order, spin order, momentum distribution, and localization properties as a function of charge-transfer energy for the the underdoped regime. In contrast with the stripe and spiral orders under hole doping, we find that the corresponding $1/8$ electron-doped system exhibits purely antiferromagnetic order in the three-band model, consistent with the asymmetry between electron and hole doping in the phase diagram of cuprates.

arXiv:2008.06527•Physical review. B./Physical review. B•Dec 24, 2020•13 citations

Pairing correlations in the cuprates: A numerical study of the three-band Hubbard model

Peizhi Mai, Giovanni Balduzzi, Steven Johnston +1 more

We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, interaction strength, and charge-transfer energy. Using the DCA, we also solve the Bethe-Salpeter equation for the two-particle Green's function in the particle-particle channel to determine the transition temperature to the superconducting phase on smaller clusters. Our calculations reproduce many aspects of the cuprate phase diagram and indicate that there is an ``optimal'' value of the charge-transfer energy for the model where ${T}_{c}$ is largest. These results have implications for our understanding of superconductivity in both the cuprates and other doped charge-transfer insulators.

arXiv:2101.10965•Physical review. B./Physical review. B•Apr 19, 2021•24 citations

Mott localisation, antiferromagnetism and nematic superconductivity in the extended Hubbard model with non-local pairing

ALESSANDRO GORI

High temperature superconductivity in copper oxides (cuprates) represents one of the most striking phenomena in condensed matter physics. These materials host superconducting states with the highest known critical temperatures, emerging upon doping antiferromagnetic insulator compounds. The antiferromagnetic (AF) insulating and superconducting (SC) phases defy standard descriptions in terms of band picture and conventional phonon-mediated superconductivity. In fact, these phenomena are genuine signatures of strong electronic correlations which hinder effective single-particle descriptions of the many-electron system. After 40 years since their discovery, it is widely accepted that the Hubbard model on the square lattice provides a minimal description of the low-energy electronic properties of cuprates. The model describes the competition between the electronic kinetic energy, modelled by nearest-neighbours hopping t, and the Coulomb repulsion, modelled by a purely local interaction potential U. At half-filling, this competition results in the so-called "Mott localisation" phenomenon, which distinguishes the metallic phase for U/t « 1 from an insulating phase for U/t » 1. The latter is characterised by the almost complete localisation of electrons due to the high energetic cost of hopping towards occupied sites. Interestingly, this basic competition is intimately related to AF ordering and, upon doping, to the formation of the SC state with characteristic d-wave symmetry. The fact that a purely repulsive interaction U is responsible for the formation of a SC state is remarkable, being superconductivity typically associated with an effective attractive interaction. In fact, a d-wave SC state can also be obtained, in general, by considering an attractive potential between electrons of neighbouring sites. Motivated by this observation, in this thesis, we study a modified version of the Hubbard model, dubbed as the "extended Hubbard model", in which we include a nearest-neighbours attractive potential -V. The general goal of this work is to explore the interplay between the purely local repulsion and the non-local attraction. This is not only motivated by the fact that both interactions are, in principle, responsible for d-wave superconductivity, but also by recent experimental observations which strongly suggest that the term -V needs to be included in order to accurately describe the magnetic spectra of some specific compounds. In this thesis, we study the extended Hubbard model using the Hartree-Fock method. Namely, we approximate the ground state by implementing a variational search in the subspace of the many-body Hilbert space of gaussian states. Despite its simplicity, the Hartree-Fock method provides a powerful framework to study competition among different phases, originating from different symmetry breaking mechanisms. In particular, we study the competition between three different phases: the normal phase, characterised by no symmetry breaking, the AF phase, characterised by a lowering of SU(2) symmetry, and the SC phase, characterised by the breaking of U(1) charge conservation. This thesis work led to a number of original results, which can be summarised by the general observation that the non-local attraction -V acts cooperatively with the local repulsion U for the stabilisation of the various phases hosted by the Hubbard model. For the normal phase, we find that the non-local attraction acts as a source of Mott localisation, leading to an effective reduction of the electronic bandwidth. In particular, for values of V larger than a critical value and close to half-filling, we find a strong Mott localisation characterised by the complete flattening of the electronic bands. Interestingly, this result shows an example of Mott localisation that can be described at the level of the Hartree-Fock approximation, whereas Mott localisation in the conventional Hubbard model usually requires more sophisticated approximation schemes. For the AF phase, we restrict our analysis to the half-filled case. We find that the pairing-induced Mott localisation leads to a significant enhancement of the magnetisation, primarily originating from the local repulsion U. Finally, we focus on the superconducting phase and derive self-consistency equations in both singlet and triplet Cooper pairing channels, showing how Cooper pairing can be mediated simultaneously by both the local repulsion and the non-local attraction. Therefore, we numerically investigate the singlet channel and find evidence of "symmetry mixing", namely, the formation of a superconducting gap that is neither symmetric nor anti-symmetric under rotations of π/2 in the plane. We show that symmetry mixing is connected to the emergence of a finite nematic order parameter, that signals the breaking of planar rotational symmetry. This nematic symmetry breaking results in the fact that the renormalisation of bands becomes anisotropic, producing a difference in electronic kinetics in the x and y directions. To summarise, the original results obtained in this thesis show that the non-local pairing provides a different route towards the formation of key strongly correlated phenomena, including Mott localisation and nematic superconductivity, and call for future investigation using more sophisticated methods such as, for example, the Dynamical Mean Field Theory.

•Jun 17, 2026•0 citations

Fluctuating charge-density-wave correlations in the three-band Hubbard model

Peizhi Mai, Benjamin Cohen-Stead, Thomas Maier +1 more

The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in nontrivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using nonperturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model's fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.

arXiv:2405.13164•Proceedings of the National Academy of Sciences•Dec 5, 2024•6 citations
Recent scanning tunneling spectroscopy (STS) experiments revealed a remarkable role of a supercell consisting of 4\ifmmode×\else\texttimes\fi4CuO24\ifmmode\times\else\texttimes\fi{}4 {\mathrm{CuO}}_{2}4\ifmmode×\else\texttimes\fi4CuO2​ unit cells in the emergence of local nematic state and preformed local Cooper pairs and phase coherent cuprate superconductivity. By employing the numerically exact determinant Quantum Monte Carlo simulations, we mimic the effects of an experimental Ca vacancy by an external local potential to investigate the charge and spectral properties of the system hosting two-doped holes. The model numerically supports the role of the 4\ifmmode×\else\texttimes\fi44\ifmmode\times\else\texttimes\fi{}44\ifmmode×\else\texttimes\fi4 supercell as the building block of hole-doped cuprates via the hole density distribution and local spectra around the local potential. Our results might provide a theoretical support for the experimental observations and a platform for investigating local charge order and local Cooper pairs on the 4\ifmmode×\else\texttimes\fi44\ifmmode\times\else\texttimes\fi{}44\ifmmode×\else\texttimes\fi4 supercell as the plausible route to understanding unconventional cuprate superconductivity.
Recent scanning tunnelling spectroscopy (STS) experiments revealed remarkable role of a supercell consisting 4×44\times44×4 CuO2_22​ unit cells in the emergence of local nematic state and preformed local Cooper pairs and phase coherent cuprate superconductivity. By employing the numerically exact determinant Quantum Monte Carlo simulations, we mimic the effects of experimental Ca vacancy by an external local potential to investigate the charge and spectral properties of the system hosting two doped holes. The model numerically support the role of the 4×44\times44×4 supercell as the building block of hole doped cuprates via the hole density distribution and local spectra around the local potential. Our results might provide a theoretical support on the experimental observations and a platform for investigating local charge order and local Cooper pairs on the 4×44\times44×4 supercell as the plausible route to understanding unconventional cuprate superconductivity.
Understanding the mechanism of Cooper pairing amounts to determining the effective interaction that operates at low energies. Efforts to achieve such a goal for superconducting materials, especially strongly correlated ones, from both bottom-up and top-down approaches, have been plagued by having to use uncontrolled approximations. Here, we perform large-scale, numerically exact, sign-problem-free zero-temperature quantum Monte Carlo simulations on an effective theory based on ``hot spots'' plus fluctuating collective modes. Because hot spots are clearly identified by angle-resolved photoemission spectroscopy for electron-doped cuprates, we focus our attention on such materials. Our goal is to determine the minimum effective action that can describe the observed superconductivity and charge-density wave. The results suggest that antiferromagnetic fluctuation alone is not sufficient---the effective action needs to be amended with nematic fluctuations. We believe that our results address the pairing mechanism of high-Tc{T}_{c}Tc​ superconductivity in electron-doped cuprates, and they shed light on the pairing mechanism of hole-doped cuprates.
We perform extensive auxiliary-field quantum Monte Carlo (AFQMC) calculations for the three-band Hubbard (Emery) model in order to study the ground-state properties of copper-oxygen planes in the cuprates. Employing cutting-edge AFQMC techniques with a self-consistent gauge constraint in auxiliary-field space to control the sign problem, we reach supercells containing $\ensuremath{\approx}500$ atoms to capture collective modes in the charge and spin orders and characterize the behavior in the thermodynamic limit. The self-consistency scheme interfacing with generalized Hartree-Fock calculations allows high accuracy in AFQMC to resolve small energy scales, which is crucial for determining the complex candidate orders in such a system. We present detailed information on the charge order, spin order, momentum distribution, and localization properties as a function of charge-transfer energy for the the underdoped regime. In contrast with the stripe and spiral orders under hole doping, we find that the corresponding 1/81/81/8 electron-doped system exhibits purely antiferromagnetic order in the three-band model, consistent with the asymmetry between electron and hole doping in the phase diagram of cuprates.
We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, interaction strength, and charge-transfer energy. Using the DCA, we also solve the Bethe-Salpeter equation for the two-particle Green's function in the particle-particle channel to determine the transition temperature to the superconducting phase on smaller clusters. Our calculations reproduce many aspects of the cuprate phase diagram and indicate that there is an ``optimal'' value of the charge-transfer energy for the model where Tc{T}_{c}Tc​ is largest. These results have implications for our understanding of superconductivity in both the cuprates and other doped charge-transfer insulators.
High temperature superconductivity in copper oxides (cuprates) represents one of the most striking phenomena in condensed matter physics. These materials host superconducting states with the highest known critical temperatures, emerging upon doping antiferromagnetic insulator compounds. The antiferromagnetic (AF) insulating and superconducting (SC) phases defy standard descriptions in terms of band picture and conventional phonon-mediated superconductivity. In fact, these phenomena are genuine signatures of strong electronic correlations which hinder effective single-particle descriptions of the many-electron system. After 40 years since their discovery, it is widely accepted that the Hubbard model on the square lattice provides a minimal description of the low-energy electronic properties of cuprates. The model describes the competition between the electronic kinetic energy, modelled by nearest-neighbours hopping t, and the Coulomb repulsion, modelled by a purely local interaction potential U. At half-filling, this competition results in the so-called "Mott localisation" phenomenon, which distinguishes the metallic phase for U/t « 1 from an insulating phase for U/t » 1. The latter is characterised by the almost complete localisation of electrons due to the high energetic cost of hopping towards occupied sites. Interestingly, this basic competition is intimately related to AF ordering and, upon doping, to the formation of the SC state with characteristic d-wave symmetry. The fact that a purely repulsive interaction U is responsible for the formation of a SC state is remarkable, being superconductivity typically associated with an effective attractive interaction. In fact, a d-wave SC state can also be obtained, in general, by considering an attractive potential between electrons of neighbouring sites. Motivated by this observation, in this thesis, we study a modified version of the Hubbard model, dubbed as the "extended Hubbard model", in which we include a nearest-neighbours attractive potential -V. The general goal of this work is to explore the interplay between the purely local repulsion and the non-local attraction. This is not only motivated by the fact that both interactions are, in principle, responsible for d-wave superconductivity, but also by recent experimental observations which strongly suggest that the term -V needs to be included in order to accurately describe the magnetic spectra of some specific compounds. In this thesis, we study the extended Hubbard model using the Hartree-Fock method. Namely, we approximate the ground state by implementing a variational search in the subspace of the many-body Hilbert space of gaussian states. Despite its simplicity, the Hartree-Fock method provides a powerful framework to study competition among different phases, originating from different symmetry breaking mechanisms. In particular, we study the competition between three different phases: the normal phase, characterised by no symmetry breaking, the AF phase, characterised by a lowering of SU(2) symmetry, and the SC phase, characterised by the breaking of U(1) charge conservation. This thesis work led to a number of original results, which can be summarised by the general observation that the non-local attraction -V acts cooperatively with the local repulsion U for the stabilisation of the various phases hosted by the Hubbard model. For the normal phase, we find that the non-local attraction acts as a source of Mott localisation, leading to an effective reduction of the electronic bandwidth. In particular, for values of V larger than a critical value and close to half-filling, we find a strong Mott localisation characterised by the complete flattening of the electronic bands. Interestingly, this result shows an example of Mott localisation that can be described at the level of the Hartree-Fock approximation, whereas Mott localisation in the conventional Hubbard model usually requires more sophisticated approximation schemes. For the AF phase, we restrict our analysis to the half-filled case. We find that the pairing-induced Mott localisation leads to a significant enhancement of the magnetisation, primarily originating from the local repulsion U. Finally, we focus on the superconducting phase and derive self-consistency equations in both singlet and triplet Cooper pairing channels, showing how Cooper pairing can be mediated simultaneously by both the local repulsion and the non-local attraction. Therefore, we numerically investigate the singlet channel and find evidence of "symmetry mixing", namely, the formation of a superconducting gap that is neither symmetric nor anti-symmetric under rotations of π/2 in the plane. We show that symmetry mixing is connected to the emergence of a finite nematic order parameter, that signals the breaking of planar rotational symmetry. This nematic symmetry breaking results in the fact that the renormalisation of bands becomes anisotropic, producing a difference in electronic kinetics in the x and y directions. To summarise, the original results obtained in this thesis show that the non-local pairing provides a different route towards the formation of key strongly correlated phenomena, including Mott localisation and nematic superconductivity, and call for future investigation using more sophisticated methods such as, for example, the Dynamical Mean Field Theory.
The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in nontrivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using nonperturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model's fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.