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:2610.37746

A $B_{1g}$ phonon unlocks $d$-wave superconductivity in the repulsive Hubbard model

Wen O. Wang, Tianshu Li, Xiaoyue Ma, Rong Zhang, Lei Zhang, Shruti Agarwal, Benjamin Cohen-Stead, Emily Z. Zhang, Brian Moritz, Steven Johnston, Zhi-Xun Shen, T. P. Devereaux

Published Oct 27, 2026• 0 citations
arXivHTMLPDF

Abstract

Using numerically exact finite-temperature quantum Monte Carlo, we demonstrate that a $B_{1g}$ phonon coupled to the doped $t$-$t'$ repulsive Hubbard model produces genuine $d$-wave superconductivity at a substantial transition temperature, in excess of either interaction alone. The observed trends suggest that stronger Hubbard repulsion may cooperate with the phonon to further raise the pairing scale. The enhancement persists at stronger repulsion within the accessible window, limited by the fermion sign problem, and follows a common trend across phonon frequencies $Ω$ when organized by a retardation scale $2πT/Ω$. We also find a well-defined superconducting dome in doping, echoing a characteristic feature of correlated superconductors. Our results identify the $B_{1g}$ phonon as a plausible microscopic ingredient for robust $d$-wave superconductivity in the Hubbard model and suggest a concrete, physically motivated route toward materials with higher superconducting transition temperatures.

DOI
10.48550/arxiv.2610.37746

Related papers

Finite-temperature signatures of underlying superconductivity in the electron-doped Hubbard model

Wen O. Wang, Thomas Devereaux

We perform numerically exact determinant quantum Monte Carlo simulations of the Hubbard model and analyze pairing tendencies by evaluating correlation functions at the imaginary-time midpoint ($τ=β/2$), which suppresses high-frequency weight and emphasizes low-energy physics. Using this diagnostic, we identify clear finite-temperature signatures of underlying $d$-wave superconductivity for electron doping, while finding no clear indication upon cooling for hole doping. Our analysis enables direct comparison with ground-state DMRG, revealing consistent real-space pairing patterns. These results provide a practical route to bridge the gap between finite-temperature and ground-state numerically exact simulations of the Hubbard model despite the fermion sign problem.

arXiv:2510.16616•Oct 18, 2025•0 citations

Finite-temperature superconducting correlations of the Hubbard model

Ehsan Khatami, Richard T. Scalettar, Rajiv R. P. Singh

We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square-lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model for the $d$-wave and extended $s$-wave pairing symmetries. The pairing structure factor shows a strong tendency to $d$-wave pairing and peaks at an interaction strength comparable to the bandwidth. The extended $s$-wave structure factor and correlation length are larger at higher temperatures but clearly saturate as temperature is lowered, whereas the $d$-wave counterparts, which start off lower at high temperatures, continue to rise near half filling. This rise is even more dramatic in the $d$-wave susceptibility. The convergence of NLCEs breaks down as the susceptibilities and correlation lengths become large, so we are unable to determine the onset of long-range order. However, our results extend the conclusion, previously restricted to only magnetic and charge correlations, that NLCEs offer a unique window into pairing in the Hubbard model at strong coupling.

arXiv:1503.06213•Physical Review B•Jun 18, 2015•24 citations

Sign change of the extended s -wave pairing vertex in the dynamic Hubbard model: A quantum Monte Carlo study

Karim Bouadim, Matthew Enjalran, F. Hébert +2 more

The ``dynamic'' Hubbard Hamiltonian describes interacting fermions on a lattice whose on-site repulsion is modulated by a coupling to a fluctuating bosonic field. We investigate one such model, introduced by Hirsch, using the determinant quantum Monte Carlo method. Our key result is that the extended $s$-wave pairing vertex, repulsive in the usual static Hubbard model, becomes attractive as the coupling to the fluctuating Bose field increases. The sign problem prevents us from exploring a low enough temperature to see if a superconducting transition occurs. We also observe a stabilization of antiferromagnetic correlations and the Mott gap near half-filling, and a near linear behavior of the energy as a function of particle density which indicates a tendency toward phase separation.

arXiv:0710.1356•Physical Review B•Jan 22, 2008•9 citations

Phonon enhancement of electronic order and negative isotope effect in the Hubbard-Holstein model on a square lattice

Da Wang, Wansheng Wang, Qiang-Hua Wang

In phonon mediated conventional $s$-wave superconductors, higher-frequency phonon (or smaller atomic mass) leads to a higher superconducting transition temperature, known as the isotope effect. However, in correlated systems, various competing electronic order (such as spin-density-wave, charge-density-wave, and unconventional superconductivity) arises and the effect of electron-phonon coupling on these orders is a long-standing problem. Using the functional renormalization group, here we investigated the interplay between the electron correlation and electron-phonon coupling in the Hubbard-Holstein model on a square lattice. At half-filling, we found spin-density-wave and charge-density-wave phases and the transition between them, while no superconducting phase arises. Upon finite doping, $d$-wave/$s$-wave superconductivity emerges in proximity to the spin-density-wave/charge-density-wave phase. Surprisingly, lower-frequency Holstein phonons are either less destructive or even beneficial to the various phases, resulting in a negative isotope effect. For the superconducting phases, such an effect is apparently beyond the Bardeen-Cooper-Schrieffer theory.

arXiv:1503.00471•Physical Review B•Nov 2, 2015•23 citations

Superconductivity in a Hubbard-Fröhlich model and in cuprates

Thomas M. Hardy, J. P. Hague, John H. Samson +1 more

Using the variational Monte Carlo method, we find that a relatively weak long-range electron-phonon interaction induces a $d$-wave superconducting state in doped Mott-Hubbard insulators and/or strongly correlated metals with a condensation energy significantly larger than can be obtained with Coulomb repulsion only. Moreover, the superconductivity is shown to exist for infinite on-site Coulomb repulsion without the need for additional mechanisms such as spin fluctuations to mediate $d$-wave superconductivity. We argue that our superconducting state is robust with respect to a more intricate choice of the trial-wave function and that a possible origin of high-temperature superconductivity lies in a proper combination of strong electron-electron correlations with poorly screened Fr\"ohlich electron-phonon interaction.

arXiv:0806.2810•Physical Review B•Jun 2, 2009•51 citations

Suppression of d -wave superconductivity in the checkerboard Hubbard model

D. G. S. P. Doluweera, Alexandru Macridin, Thomas Maier +2 more

Using a dynamical cluster quantum Monte Carlo approximation, we investigate the $d$-wave superconducting transition temperature ${T}_{c}$ in the doped two-dimensional repulsive Hubbard model with a weak inhomogeneity. The inhomogeneity is introduced in the hoppings ${t}^{\ensuremath{'}}$ and $t$ in the form of a checkerboard pattern where $t$ is the hopping within a $2\ifmmode\times\else\texttimes\fi{}2$ plaquette, and ${t}^{\ensuremath{'}}$ is the hopping between the plaquettes. We find inhomogeneity suppresses ${T}_{c}$. The characteristic spin excitation energy and the strength of $d$-wave pairing interaction decrease with decreasing ${T}_{c}$, suggesting a strong correlation between these quantities.

arXiv:0806.1504•Physical Review B•Jul 9, 2008•19 citations
A B1gB_{1g}B1g​ phonon unlocks ddd-wave superconductivity in the repulsive Hubbard model
Using numerically exact finite-temperature quantum Monte Carlo, we demonstrate that a B1gB_{1g}B1g​ phonon coupled to the doped ttt-t′t't′ repulsive Hubbard model produces genuine ddd-wave superconductivity at a substantial transition temperature, in excess of either interaction alone. The observed trends suggest that stronger Hubbard repulsion may cooperate with the phonon to further raise the pairing scale. The enhancement persists at stronger repulsion within the accessible window, limited by the fermion sign problem, and follows a common trend across phonon frequencies ΩΩΩ when organized by a retardation scale 2πT/Ω2πT/Ω2πT/Ω. We also find a well-defined superconducting dome in doping, echoing a characteristic feature of correlated superconductors. Our results identify the B1gB_{1g}B1g​ phonon as a plausible microscopic ingredient for robust ddd-wave superconductivity in the Hubbard model and suggest a concrete, physically motivated route toward materials with higher superconducting transition temperatures.
We perform numerically exact determinant quantum Monte Carlo simulations of the Hubbard model and analyze pairing tendencies by evaluating correlation functions at the imaginary-time midpoint (τ=β/2τ=β/2τ=β/2), which suppresses high-frequency weight and emphasizes low-energy physics. Using this diagnostic, we identify clear finite-temperature signatures of underlying ddd-wave superconductivity for electron doping, while finding no clear indication upon cooling for hole doping. Our analysis enables direct comparison with ground-state DMRG, revealing consistent real-space pairing patterns. These results provide a practical route to bridge the gap between finite-temperature and ground-state numerically exact simulations of the Hubbard model despite the fermion sign problem.
We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square-lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model for the ddd-wave and extended sss-wave pairing symmetries. The pairing structure factor shows a strong tendency to ddd-wave pairing and peaks at an interaction strength comparable to the bandwidth. The extended sss-wave structure factor and correlation length are larger at higher temperatures but clearly saturate as temperature is lowered, whereas the ddd-wave counterparts, which start off lower at high temperatures, continue to rise near half filling. This rise is even more dramatic in the ddd-wave susceptibility. The convergence of NLCEs breaks down as the susceptibilities and correlation lengths become large, so we are unable to determine the onset of long-range order. However, our results extend the conclusion, previously restricted to only magnetic and charge correlations, that NLCEs offer a unique window into pairing in the Hubbard model at strong coupling.
The ``dynamic'' Hubbard Hamiltonian describes interacting fermions on a lattice whose on-site repulsion is modulated by a coupling to a fluctuating bosonic field. We investigate one such model, introduced by Hirsch, using the determinant quantum Monte Carlo method. Our key result is that the extended sss-wave pairing vertex, repulsive in the usual static Hubbard model, becomes attractive as the coupling to the fluctuating Bose field increases. The sign problem prevents us from exploring a low enough temperature to see if a superconducting transition occurs. We also observe a stabilization of antiferromagnetic correlations and the Mott gap near half-filling, and a near linear behavior of the energy as a function of particle density which indicates a tendency toward phase separation.
In phonon mediated conventional sss-wave superconductors, higher-frequency phonon (or smaller atomic mass) leads to a higher superconducting transition temperature, known as the isotope effect. However, in correlated systems, various competing electronic order (such as spin-density-wave, charge-density-wave, and unconventional superconductivity) arises and the effect of electron-phonon coupling on these orders is a long-standing problem. Using the functional renormalization group, here we investigated the interplay between the electron correlation and electron-phonon coupling in the Hubbard-Holstein model on a square lattice. At half-filling, we found spin-density-wave and charge-density-wave phases and the transition between them, while no superconducting phase arises. Upon finite doping, ddd-wave/sss-wave superconductivity emerges in proximity to the spin-density-wave/charge-density-wave phase. Surprisingly, lower-frequency Holstein phonons are either less destructive or even beneficial to the various phases, resulting in a negative isotope effect. For the superconducting phases, such an effect is apparently beyond the Bardeen-Cooper-Schrieffer theory.
Using the variational Monte Carlo method, we find that a relatively weak long-range electron-phonon interaction induces a ddd-wave superconducting state in doped Mott-Hubbard insulators and/or strongly correlated metals with a condensation energy significantly larger than can be obtained with Coulomb repulsion only. Moreover, the superconductivity is shown to exist for infinite on-site Coulomb repulsion without the need for additional mechanisms such as spin fluctuations to mediate ddd-wave superconductivity. We argue that our superconducting state is robust with respect to a more intricate choice of the trial-wave function and that a possible origin of high-temperature superconductivity lies in a proper combination of strong electron-electron correlations with poorly screened Fr\"ohlich electron-phonon interaction.
Using a dynamical cluster quantum Monte Carlo approximation, we investigate the ddd-wave superconducting transition temperature Tc{T}_{c}Tc​ in the doped two-dimensional repulsive Hubbard model with a weak inhomogeneity. The inhomogeneity is introduced in the hoppings ${t}^{\ensuremath{'}}$ and ttt in the form of a checkerboard pattern where ttt is the hopping within a $2\ifmmode\times\else\texttimes\fi{}2$ plaquette, and ${t}^{\ensuremath{'}}$ is the hopping between the plaquettes. We find inhomogeneity suppresses Tc{T}_{c}Tc​. The characteristic spin excitation energy and the strength of ddd-wave pairing interaction decrease with decreasing Tc{T}_{c}Tc​, suggesting a strong correlation between these quantities.