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

Dual parquet scheme for the two-dimensional Hubbard model: Modeling low-energy physics of high- T c cuprates with high momentum resolution

Grigory V. Astretsov, G. Rohringer, Alexey N. Rubtsov

Published Feb 10, 2020• 48 citations
arXivHTMLPDF

Abstract

We present a method to treat the two-dimensional (2D) Hubbard model for parameter regimes which are relevant for the physics of the high-${T}_{c}$ superconducting cuprates. Unlike previous attempts to attack this problem, our approach takes into account all fluctuations in different channels on equal footing and is able to treat reasonable large lattice sizes up to $32\ifmmode\times\else\texttimes\fi{}32$. This is achieved by the following three-step procedure: (i) We transform the original problem to a representation (dual fermions) in which all purely local correlation effects from the dynamical mean field theory are already considered in the bare propagator and bare interaction of the new problem. (ii) The strong $1/{(i\ensuremath{\nu})}^{2}$ decay of the bare propagator allows us to integrate out all higher Matsubara frequencies aside from the lowest using low-order diagrams. The new effective action depends only on the two lowest Matsubara frequencies which allows us to (iii) apply the two-particle self-consistent parquet formalism, which takes into account the competition between different low-energy bosonic modes in an unbiased way, on much finer momentum grids than usual. In this way, we were able to map out the phase diagram of the 2D Hubbard model as a function of temperature and doping. Consistently with the experimental evidence for hole-doped cuprates and previous dynamical cluster approximation calculations, we find an antiferromagnetic region at low doping and a superconducting dome at higher doping. Our results also support the role of the van Hove singularity as an important ingredient for the high value of ${T}_{c}$ at optimal doping. At small doping, the destruction of antiferromagnetism is accompanied by an increase of charge fluctuations supporting the scenario of a phase-separated state driven by quantum critical fluctuations.

DOI
10.1103/physrevb.101.075109

Related papers

Superconductivity, antiferromagnetism, and phase separation in the two-dimensional Hubbard model: A dual-fermion approach

Junya Otsuki, Hartmut Hafermann, Alexander I. Lichtenstein

The dual-fermion approach offers a way to perform diagrammatic expansion around the dynamical mean field theory. Using this formalism, the influence of antiferromagnetic fluctuations on the self-energy is taken into account through ladder-type diagrams in the particle-hole channel. The resulting phase diagram for the (quasi-)two-dimensional Hubbard model exhibits antiferromagnetism and $d$-wave superconductivity. Furthermore, a uniform charge instability, i.e., phase separation, is obtained in the low-doping regime around the Mott insulator. We also examine spin/charge density wave fluctuations including $d$-wave symmetry. The model exhibits a tendency towards an unconventional charge density wave, but no divergence of the susceptibility is found.

arXiv:1410.1246•Physical Review B•Dec 18, 2014•98 citations

Two-particle level diagrammatic approaches for strongly correlated systems

Herbert F. Fotso

This thesis presents the development of new numerical methods for the treatment of strongly correlated electron systems based on self-consistent approaches at both the one and the two-particle level such as the parquet formalism. The parquet formalism was solved for the first time on a two-dimensional cluster. When the fully irreducible vertex is approximated by the bare vertex, we obtain the parquet approximation. Its validity was investigated by comparing results that it produces to those of other conserving approximations such as the FLuctuation EXchange (FLEX) approximation or the Second Order Perturbation Theory (SOPT). We found that it provides a significant improvement of FLEX or SOPT and a satisfactory agreement with Quantum Monte Carlo results despite instabilities in the self-consistency at low temperatures and for strong Coulomb interaction. We use the parquet formalism to study the Quantum Critical Point at finite doping in the Hubbard model by decomposing the vertex into its contributions from different channels. We apply this decomposition to the pairing channel and we find that the dominant contribution to the vertex originates in the spin channel even at the quantum critical doping. Furthermore, we explore the divergence of the two parts of the pairing matrix at optimal doping and observe that the irreducible vertex decreases monotonically as the doping is increased while the bare susceptibility exhibits an algebraic divergence at the quantum critical doping supporting the Quantum Critical BCS scenario proposed by She and Zaanen. To circumvent the instabilities in the iteration of the parquet formalism, we explored the dual fermion approach introduced by Rubtsov et al. Here, we extended the formalism to the Dynamical Cluster Approximation, in the process introducing a small parameter in the dual fermions perturbation theory. We demonstrate the quality of the resulting Dual Fermion DCA through a systematic study of the cluster size dependence and of the different perturbative approximations. These efforts represent the initial steps in the development of the Multi-Scale Many Body approach that appropriately treats correlations at different length scales.

Civil War Book Review•Jul 1, 2011•1 citation

Superfluid ground state phase diagram of the two-dimensional Hubbard model in the emergent Bardeen-Cooper-Schrieffer regime

Fedor Šimkovic, Youjin Deng, Evgeny Kozik

In nonperturbative regimes, the superfluid instability in the two-dimensional (2D) Hubbard model can be described by an emergent BCS theory with small effective pairing constants. We compute the effective couplings using a controlled bold-line diagrammatic Monte Carlo approach, which stochastically sums all skeleton Feynman diagrams dressed in the one- and two-particle channels to high expansion orders, and map out the resulting superfluid ground-state phase diagram in a range of next-nearest-neighbor hopping $0\ensuremath{\le}{t}^{\ensuremath{'}}\ensuremath{\le}0.3t$, interaction strength $0\ensuremath{\le}U\ensuremath{\le}3t$, and lattice filling $0\ensuremath{\le}n\ensuremath{\le}2$. The phase diagram is dramatically transformed in the hole-doped region and becomes particularly rich at larger doping and ${t}^{\ensuremath{'}}$. At ${t}^{\ensuremath{'}}=0.3$, the weak-coupling picture with the dominant triplet pairing sharply peaked at $n\ensuremath{\approx}0.73$ due to the Van Hove singularity is replaced by a plateau of the singlet ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ paring, while for $U\ensuremath{\gtrsim}3t$ the effective couplings are consistent with the ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ high-temperature superconductivity in the hole-doped region near cuprates' optimal doping.

arXiv:1912.13054•Physical review. B./Physical review. B•Jul 20, 2021•10 citations

Next generation multi-scale quantum simulations for strongly correlated materials

Shuxiang Yang

This thesis represents our effort to develop the next generation multi-scale quantum simulation methods suitable for strongly-correlated materials, where complicated phase-diagrams prevail, suggesting complicated underlying physics. We first give a detailed description of the parquet formalism. With its help, different approximate methods can be unified and a hierarchy of approximate methods with different accuracies and computational complexity can thus be designed. Next, we present a numerical solution of the parquet approximation. Results on the Hubbard model are compared to those obtained from Determinant Quantum Monte Carlo (DQMC), FLuctuation EXchange (FLEX), and self-consistent second-order approximation methods. The comparison shows a satisfactory agreement with DQMC and a significant improvement over the FLEX or the self-consistent second-order approximation. The parquet formalism can also be used to analyze the superconducting mechanism of the high-temperature superconductors. The dynamical cluster approximation (DCA) method is used to understand the proximity of the superconducting dome to the quantum critical point in the 2-D Hubbard model. At optimal doping, where Vd is revealed to be featureless, we find a power-law behavior of chi_0d(w=0), replacing the BCS logarithm behavior, and strongly enhanced T_c. After that we propose another multi-scale approach by combining the DCA and the recently introduced dual-fermion formalism. Within this approach, short and long length scale physics is addressed by the DCA cluster calculation, while intermediate length scale physics is addressed diagrammatically using dual fermions. The bare and dressed dual fermionic Green functions scale as O(1/Lc), so perturbation theory on the dual lattice converges very quickly. Lastly, we study the responses to dynamical modulation of the optical lattice potential by analyzing properties of the repulsive fermionic Hubbard model in an optical lattice. We provide numerical evidence showing the modulations by on-site local interaction cannot be ignored, and instead can even strongly contribute to the dynamical behaviors of the system in highly-doped cases.

Civil War Book Review•Aug 22, 2011•2 citations

Evolution of the superconductivity dome in the two-dimensional Hubbard model

K.-S. Chen, Zi Yang Meng, Shuxiang Yang +3 more

In a recent publication [Chen et al., Phys. Rev. B 86, 165136 (2012)], we identified a line of Lifshitz transition points separating the Fermi liquid and pseudogap regions in the hole-doped two-dimensional Hubbard model. Here, we extend the study to further determine the superconducting transition temperature in the phase diagram. By means of large-scale dynamical cluster quantum Monte Carlo simulations, we are able to identify the evolution of the $d$-wave superconducting dome in the hole-dope side of the phase diagram, with next-nearest-neighbor hopping (${t}^{\ensuremath{'}}$), chemical potential, and temperature as control parameters. To obtain the superconducting transition temperature ${T}_{c}$, we employ two-particle measurements of the pairing susceptibilities. As ${t}^{\ensuremath{'}}$ goes from positive to negative values, we find the $d$-wave projected irreducible pairing vertex function is enhanced, and the curvature of its doping dependence changes from convex to concave, which fixes the position of the maximum superconducting temperature at the same filling ($n\ensuremath{\approx}0.85$) and constraints the dome from precisely following the Lifshitz line. We furthermore decompose the irreducible vertex function into fully irreducible, charge and spin components via the parquet equations, and consistently find that the spin component dominates the pairing vertex function in the doping range where the dome is located. Our investigations deepen the understanding of the phase diagram of the two-dimensional Hubbard model and, more importantly, pose new questions to the field. For example, we found as ${t}^{\ensuremath{'}}$ goes from positive to negative values, the curvature of the pairing strength as a function of doping changes from convex to concave, and the nature of the dominant fluctuations changes from charge degree of freedom to spin degree of freedom. The study of these issues will lead to further understanding of the phase diagram of the two-dimensional Hubbard model and also the physics of the hole-doped cuprate high-temperature superconductors.

arXiv:1308.5946•Physical Review B•Dec 10, 2013•33 citations

Antiferromagnetic to superconducting phase transition in the hole- and electron-doped Hubbard model at zero temperature

Markus Aichhorn, Enrico Arrigoni, Michael Potthoff +1 more

The competition between $d$-wave superconductivity (SC) and antiferromagnetism (AF) in the high-${T}_{c}$ cuprates is investigated by studying the hole- and electron-doped two-dimensional Hubbard model with a recently proposed variational quantum-cluster theory. The approach is shown to provide a thermodynamically consistent determination of the particle number, provided that an overall shift of the on-site energies is treated as a variational parameter. The consequences for the single-particle excitation spectra and for the phase diagram are explored. By comparing the single-particle spectra with quantum Monte Carlo and experimental data, we verify that the low-energy excitations in a strongly correlated electronic system are described appropriately. The cluster calculations also reproduce the overall ground-state phase diagram of the high-temperature superconductors. In particular, they include salient features such as the enhanced robustness of the antiferromagnetic state as a function of electron doping and the tendency towards phase separation into a mixed antiferromagnetic-superconducting phase at low doping and a pure superconducting phase at high (both hole and electron) doping.

arXiv:cond-mat/0511460•Physical Review B•Jul 13, 2006•154 citations
We present a method to treat the two-dimensional (2D) Hubbard model for parameter regimes which are relevant for the physics of the high-Tc{T}_{c}Tc​ superconducting cuprates. Unlike previous attempts to attack this problem, our approach takes into account all fluctuations in different channels on equal footing and is able to treat reasonable large lattice sizes up to 32\ifmmode×\else\texttimes\fi3232\ifmmode\times\else\texttimes\fi{}3232\ifmmode×\else\texttimes\fi32. This is achieved by the following three-step procedure: (i) We transform the original problem to a representation (dual fermions) in which all purely local correlation effects from the dynamical mean field theory are already considered in the bare propagator and bare interaction of the new problem. (ii) The strong 1/(i\ensuremathν)21/{(i\ensuremath{\nu})}^{2}1/(i\ensuremathν)2 decay of the bare propagator allows us to integrate out all higher Matsubara frequencies aside from the lowest using low-order diagrams. The new effective action depends only on the two lowest Matsubara frequencies which allows us to (iii) apply the two-particle self-consistent parquet formalism, which takes into account the competition between different low-energy bosonic modes in an unbiased way, on much finer momentum grids than usual. In this way, we were able to map out the phase diagram of the 2D Hubbard model as a function of temperature and doping. Consistently with the experimental evidence for hole-doped cuprates and previous dynamical cluster approximation calculations, we find an antiferromagnetic region at low doping and a superconducting dome at higher doping. Our results also support the role of the van Hove singularity as an important ingredient for the high value of Tc{T}_{c}Tc​ at optimal doping. At small doping, the destruction of antiferromagnetism is accompanied by an increase of charge fluctuations supporting the scenario of a phase-separated state driven by quantum critical fluctuations.
The dual-fermion approach offers a way to perform diagrammatic expansion around the dynamical mean field theory. Using this formalism, the influence of antiferromagnetic fluctuations on the self-energy is taken into account through ladder-type diagrams in the particle-hole channel. The resulting phase diagram for the (quasi-)two-dimensional Hubbard model exhibits antiferromagnetism and ddd-wave superconductivity. Furthermore, a uniform charge instability, i.e., phase separation, is obtained in the low-doping regime around the Mott insulator. We also examine spin/charge density wave fluctuations including ddd-wave symmetry. The model exhibits a tendency towards an unconventional charge density wave, but no divergence of the susceptibility is found.
This thesis presents the development of new numerical methods for the treatment of strongly correlated electron systems based on self-consistent approaches at both the one and the two-particle level such as the parquet formalism. The parquet formalism was solved for the first time on a two-dimensional cluster. When the fully irreducible vertex is approximated by the bare vertex, we obtain the parquet approximation. Its validity was investigated by comparing results that it produces to those of other conserving approximations such as the FLuctuation EXchange (FLEX) approximation or the Second Order Perturbation Theory (SOPT). We found that it provides a significant improvement of FLEX or SOPT and a satisfactory agreement with Quantum Monte Carlo results despite instabilities in the self-consistency at low temperatures and for strong Coulomb interaction. We use the parquet formalism to study the Quantum Critical Point at finite doping in the Hubbard model by decomposing the vertex into its contributions from different channels. We apply this decomposition to the pairing channel and we find that the dominant contribution to the vertex originates in the spin channel even at the quantum critical doping. Furthermore, we explore the divergence of the two parts of the pairing matrix at optimal doping and observe that the irreducible vertex decreases monotonically as the doping is increased while the bare susceptibility exhibits an algebraic divergence at the quantum critical doping supporting the Quantum Critical BCS scenario proposed by She and Zaanen. To circumvent the instabilities in the iteration of the parquet formalism, we explored the dual fermion approach introduced by Rubtsov et al. Here, we extended the formalism to the Dynamical Cluster Approximation, in the process introducing a small parameter in the dual fermions perturbation theory. We demonstrate the quality of the resulting Dual Fermion DCA through a systematic study of the cluster size dependence and of the different perturbative approximations. These efforts represent the initial steps in the development of the Multi-Scale Many Body approach that appropriately treats correlations at different length scales.
In nonperturbative regimes, the superfluid instability in the two-dimensional (2D) Hubbard model can be described by an emergent BCS theory with small effective pairing constants. We compute the effective couplings using a controlled bold-line diagrammatic Monte Carlo approach, which stochastically sums all skeleton Feynman diagrams dressed in the one- and two-particle channels to high expansion orders, and map out the resulting superfluid ground-state phase diagram in a range of next-nearest-neighbor hopping $0\ensuremath{\le}{t}^{\ensuremath{'}}\ensuremath{\le}0.3t$, interaction strength $0\ensuremath{\le}U\ensuremath{\le}3t$, and lattice filling $0\ensuremath{\le}n\ensuremath{\le}2$. The phase diagram is dramatically transformed in the hole-doped region and becomes particularly rich at larger doping and ${t}^{\ensuremath{'}}$. At ${t}^{\ensuremath{'}}=0.3$, the weak-coupling picture with the dominant triplet pairing sharply peaked at $n\ensuremath{\approx}0.73$ due to the Van Hove singularity is replaced by a plateau of the singlet ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ paring, while for $U\ensuremath{\gtrsim}3t$ the effective couplings are consistent with the ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ high-temperature superconductivity in the hole-doped region near cuprates' optimal doping.
This thesis represents our effort to develop the next generation multi-scale quantum simulation methods suitable for strongly-correlated materials, where complicated phase-diagrams prevail, suggesting complicated underlying physics. We first give a detailed description of the parquet formalism. With its help, different approximate methods can be unified and a hierarchy of approximate methods with different accuracies and computational complexity can thus be designed. Next, we present a numerical solution of the parquet approximation. Results on the Hubbard model are compared to those obtained from Determinant Quantum Monte Carlo (DQMC), FLuctuation EXchange (FLEX), and self-consistent second-order approximation methods. The comparison shows a satisfactory agreement with DQMC and a significant improvement over the FLEX or the self-consistent second-order approximation. The parquet formalism can also be used to analyze the superconducting mechanism of the high-temperature superconductors. The dynamical cluster approximation (DCA) method is used to understand the proximity of the superconducting dome to the quantum critical point in the 2-D Hubbard model. At optimal doping, where Vd is revealed to be featureless, we find a power-law behavior of chi_0d(w=0), replacing the BCS logarithm behavior, and strongly enhanced T_c. After that we propose another multi-scale approach by combining the DCA and the recently introduced dual-fermion formalism. Within this approach, short and long length scale physics is addressed by the DCA cluster calculation, while intermediate length scale physics is addressed diagrammatically using dual fermions. The bare and dressed dual fermionic Green functions scale as O(1/Lc), so perturbation theory on the dual lattice converges very quickly. Lastly, we study the responses to dynamical modulation of the optical lattice potential by analyzing properties of the repulsive fermionic Hubbard model in an optical lattice. We provide numerical evidence showing the modulations by on-site local interaction cannot be ignored, and instead can even strongly contribute to the dynamical behaviors of the system in highly-doped cases.
In a recent publication [Chen et al., Phys. Rev. B 86, 165136 (2012)], we identified a line of Lifshitz transition points separating the Fermi liquid and pseudogap regions in the hole-doped two-dimensional Hubbard model. Here, we extend the study to further determine the superconducting transition temperature in the phase diagram. By means of large-scale dynamical cluster quantum Monte Carlo simulations, we are able to identify the evolution of the ddd-wave superconducting dome in the hole-dope side of the phase diagram, with next-nearest-neighbor hopping (${t}^{\ensuremath{'}}$), chemical potential, and temperature as control parameters. To obtain the superconducting transition temperature Tc{T}_{c}Tc​, we employ two-particle measurements of the pairing susceptibilities. As ${t}^{\ensuremath{'}}$ goes from positive to negative values, we find the ddd-wave projected irreducible pairing vertex function is enhanced, and the curvature of its doping dependence changes from convex to concave, which fixes the position of the maximum superconducting temperature at the same filling ($n\ensuremath{\approx}0.85$) and constraints the dome from precisely following the Lifshitz line. We furthermore decompose the irreducible vertex function into fully irreducible, charge and spin components via the parquet equations, and consistently find that the spin component dominates the pairing vertex function in the doping range where the dome is located. Our investigations deepen the understanding of the phase diagram of the two-dimensional Hubbard model and, more importantly, pose new questions to the field. For example, we found as ${t}^{\ensuremath{'}}$ goes from positive to negative values, the curvature of the pairing strength as a function of doping changes from convex to concave, and the nature of the dominant fluctuations changes from charge degree of freedom to spin degree of freedom. The study of these issues will lead to further understanding of the phase diagram of the two-dimensional Hubbard model and also the physics of the hole-doped cuprate high-temperature superconductors.
The competition between ddd-wave superconductivity (SC) and antiferromagnetism (AF) in the high-Tc{T}_{c}Tc​ cuprates is investigated by studying the hole- and electron-doped two-dimensional Hubbard model with a recently proposed variational quantum-cluster theory. The approach is shown to provide a thermodynamically consistent determination of the particle number, provided that an overall shift of the on-site energies is treated as a variational parameter. The consequences for the single-particle excitation spectra and for the phase diagram are explored. By comparing the single-particle spectra with quantum Monte Carlo and experimental data, we verify that the low-energy excitations in a strongly correlated electronic system are described appropriately. The cluster calculations also reproduce the overall ground-state phase diagram of the high-temperature superconductors. In particular, they include salient features such as the enhanced robustness of the antiferromagnetic state as a function of electron doping and the tendency towards phase separation into a mixed antiferromagnetic-superconducting phase at low doping and a pure superconducting phase at high (both hole and electron) doping.