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

Correlations and renormalization of the electron-phonon coupling in the honeycomb Hubbard ladder and superconductivity in polyacene

George Karakonstantakis, L. Liu, Ronny Thomale, Steven Allan Kivelson

Published Dec 27, 2013• 18 citations
arXivHTMLPDF

Abstract

We have performed extensive density matrix renormalization group (DMRG) studies of the Hubbard model on a honeycomb ladder. The band structure (with Hubbard $U=0$) exhibits an unusual quadratic band touching at half-filling, which is associated with a quantum Lifshitz transition from a band insulator to a metal. For one electron per site, nonzero $U$ drives the system into an insulating state in which there is no pair-binding between added electrons; this implies that superconductivity driven directly by the repulsive electron-electron interactions is unlikely in the regime of small doping, $x\ensuremath{\ll}1$. However, the divergent density of states as $x\ensuremath{\rightarrow}0$, the large values of the phonon frequencies, and an unusual correlation induced enhancement of the electron-phonon coupling imply that lightly doped polyacenes, which approximately realize this structure, are good candidates for high-temperature electron-phonon driven superconductivity.

DOI
10.1103/physrevb.88.224512

Related papers

Renormalization of the electron-phonon interaction by strong electronic correlations in high- Tc superconductors

R. Zeyher, Miodrag L. Kulić

The renormalization of the electron-phonon interaction by strong electronic correlations is studied using a one-band Hubbard model with infinite repulsion and nearest (t model) and nearest and second-nearest (tt\ensuremath{'} model) neighbor hopping terms and an on-site electron-phonon coupling. Using Hubbard's X operators and an extension from 2 to N degrees of freedom for the electrons the leading contributions for the electron self-energy and the vertex function in 1/N and the electron-phonon coupling constant are given and numerically evaluated for a square lattice. We find that the momentum dependence of the vertex function depends strongly on doping: For large dopings it is rather weak, with decreasing doping it becomes more and more pronounced leading to a strong reduction of the vertex at larger momentum transfers until, for very small dopings, the vertex essentially consists of a forward scattering peak with a width proportional to the doping. This behavior occurs both in the t and the tt\ensuremath{'} models and also in one dimension where analytic expressions are derived. Correlation effects also change ${\mathrm{\ensuremath{\alpha}}}^{2}$F in general: The full-symmetric component ${\mathrm{\ensuremath{\alpha}}}^{2}$${\mathit{F}}_{1}$ is in the t model somewhat, in the tt\ensuremath{'} model rather strongly suppressed, especially near half-filling; the other symmetry components ${\mathrm{\ensuremath{\alpha}}}^{2}$${\mathit{F}}_{\mathit{i}}$ with i=2,...,5 increase strongly with decreasing doping and are of similar magnitude as ${\mathrm{\ensuremath{\alpha}}}^{2}$${\mathit{F}}_{1}$ near half-filling.Including also direct Coulomb repulsion nontrivial symmetries such as d wave become more stable than the s-wave order parameter below a critical value for the doping even for the considered phonon-mediated superconductivity. Most dramatic, however, is the quenching of the resistivity due to electron-phonon scattering both in the t and the tt\ensuremath{'} models at intermediate and small dopings. This result may explain the absence of phonon features in the experimental transport coefficients of high-${\mathit{T}}_{\mathit{c}}$ compounds. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Feb 1, 1996•147 citations

Phase diagram for the one-dimensional Hubbard-Holstein model: A density-matrix renormalization group study

Masaki Tezuka, Ryotaro Arita, Hideo Aoki

Phase diagram of the Hubbard-Holstein model in the coexistence of electron-electron and electron-phonon interactions has been theoretically obtained with the density-matrix renormalization group method for one-dimensional systems, where an improved warm-up (the recursive sweep) procedure has enabled us to calculate various correlation functions. We have examined the cases of (i) the systems half-filled by electrons for the full parameter space spanned by the electron-electron and electron-phonon coupling constants and the phonon frequency, (ii) non-half-filled system, and (iii) trestle lattice. For (i), we have detected a region where both the charge and on-site pairing correlations decay with power laws in real space, which suggests a metallic behavior. While pairing correlations are not dominant in (i), we have found that they become dominant as the system is doped in (ii) or as the electronic band structure is modified (with a broken electron-hole symmetry) in (iii) in certain parameter regions.

arXiv:0707.3197•Physical Review B•Oct 19, 2007•78 citations

Phase diagram and isotope effects of the quasi-one-dimensional electron gas coupled to phonons

Ian P. Bindloss

Using a multistep renormalization group method, we study the low-temperature phases of the interacting one-dimensional (1D) electron gas coupled to phonons. We obtain analytic expressions for the weak-coupling quantum phase boundaries of the 1D extended Holstein-Hubbard model and the 1D extended Peierls-Hubbard model for general band-filling and phonon frequency. Away from half-filling, the phase diagrams are characterized by a delicate competition between spin density wave, charge density wave, and superconducting orders. We study the dependence of the ground state on the electron-phonon (el-ph) and electron-electron (el-el) coupling strengths, the screening length, electron bandwidth, phonon frequency, doping, and type of phonon. Unlike the case in Fermi liquids, in 1D the el-ph coupling is strongly renormalized, often to stronger values. Even when the bare phonon-induced attraction is weak compared to the bare el-el repulsion, a small amount of retardation can cause the renormalized el-ph interaction to dominate the problem. We find cases in which a repulsive el-el interaction enhances the superconducting susceptibility in the presence of a retarded el-ph interaction. The spin gap and superconducting susceptibility are found to be strongly dependent on the deviation from half-filling (doping). In some cases, the superconducting susceptibility varies nonmonotonically with doping and exhibits a maximum at a particular doping. For a quasi-1D array of weakly coupled, fluctuating 1D chains, the superconducting transition temperature ${T}_{c}$ also exhibits a maximum as a function of doping. The effect of changing the ion mass (isotope effect) on ${T}_{c}$ is found to be largest near half-filling and to decrease rapidly with doping.

arXiv:cond-mat/0404154•Physical Review B•May 27, 2005•13 citations

Density matrix renormalization group study of superconductivity in the triangular lattice Hubbard model

Jordan Venderley, Eun-Ah Kim

With the discovery of strong-coupling physics and superconductivity in moir\'e superlattices, it is essential to have an understanding of strong-coupling driven superconductivity in systems with trigonal symmetry. The simplest lattice model with trigonal symmetry is the triangular lattice Hubbard model. Although the triangular lattice spin model is a heavily studied model in the context of frustration, studies of the hole-doped triangular lattice Hubbard model are rare. Here, we use density matrix renormalization group to investigate the dominant superconducting channels in the hole-doped triangular lattice Hubbard model over a range of repulsive interaction strengths. We find a clear transition from $p$-wave superconductivity at a moderate on-site repulsion strength ($U/t=2$) at filling above 1/4 ($n\ensuremath{\sim}0.65$) to $d$-wave superconductivity at a strong on-site repulsion strength ($U/t=10$) at filling below 1/4 ($n\ensuremath{\sim}0.4$). The unusual tunability that moir\'e superlattices offer in controlling $U/t$ would open up the opportunity to realize this transition between $d$-wave and $p$-wave superconductivity.

arXiv:1901.11034•Physical review. B./Physical review. B•Aug 26, 2019•38 citations

Stripes versus superconductivity in the doped Hubbard model on the honeycomb lattice

Mingpu Qin

We study the ground state of the doped Hubbard model on the honeycomb lattice in the small-doping and strongly interacting region. The nature of the ground state after doping holes into antiferromagnetic Mott insulating states on the honeycomb lattice remains a long-standing unsolved issue, even though tremendous effort has been spent investigating this challenging problem. In this work, we employ two complementary, state-of-the-art, many-body computational methods: constrained path (CP) auxiliary-field quantum Monte Carlo (AFQMC) with self-consistent constraint and density matrix renormalization group (DMRG) methods. Systematic and detailed cross validations are performed between these two methods for narrow systems for which DMRG can produce reliable results. AFQMC is then utilized to study wider systems to investigate the thermodynamic limit properties. The ground state is found to be a half-filled stripe state in the small-doping and strongly interacting region. The pairing correlation shows $d$-wave symmetry locally but decays exponentially with the distance between two pairs.

arXiv:2104.14160•Physical review. B./Physical review. B•Jan 7, 2022•11 citations

Doping asymmetry in the three-band Hamiltonian for cuprate ladders: Failure of the standard model of superconductivity in cuprates

Jeong-Pil Song, S. Mazumdar, Rudolf Torsten Clay

The relevance of the single-band two-dimensional Hubbard model to superconductivity in the doped cuprates has recently been questioned, based on density matrix renormalization group (DMRG) computations that found superconductivity over an unrealistically broad doping region upon electron-doping, yet a complete absence of superconductivity for hole-doping. We report very similar results from DMRG calculations on a ${\mathrm{Cu}}_{2}{\mathrm{O}}_{3}$ two-leg ladder within the parent three-band correlated-electron Hamiltonian. The strong asymmetry found in our calculations are in contradiction to the deep and profound symmetry in the experimental phase diagrams of electron- and hole-doped cuprate superconductors, as seen from the occurrence of quantum critical points within the superconducting domes in both cases that are characterized by Fermi surface reconstruction, large jumps in carrier density, and strange metal behavior.

arXiv:2302.08356•Physical review. B./Physical review. B•Jun 27, 2023•5 citations
We have performed extensive density matrix renormalization group (DMRG) studies of the Hubbard model on a honeycomb ladder. The band structure (with Hubbard U=0U=0U=0) exhibits an unusual quadratic band touching at half-filling, which is associated with a quantum Lifshitz transition from a band insulator to a metal. For one electron per site, nonzero UUU drives the system into an insulating state in which there is no pair-binding between added electrons; this implies that superconductivity driven directly by the repulsive electron-electron interactions is unlikely in the regime of small doping, x\ensuremath≪1x\ensuremath{\ll}1x\ensuremath≪1. However, the divergent density of states as x\ensuremath→0x\ensuremath{\rightarrow}0x\ensuremath→0, the large values of the phonon frequencies, and an unusual correlation induced enhancement of the electron-phonon coupling imply that lightly doped polyacenes, which approximately realize this structure, are good candidates for high-temperature electron-phonon driven superconductivity.
The renormalization of the electron-phonon interaction by strong electronic correlations is studied using a one-band Hubbard model with infinite repulsion and nearest (t model) and nearest and second-nearest (tt\ensuremath{'} model) neighbor hopping terms and an on-site electron-phonon coupling. Using Hubbard's X operators and an extension from 2 to N degrees of freedom for the electrons the leading contributions for the electron self-energy and the vertex function in 1/N and the electron-phonon coupling constant are given and numerically evaluated for a square lattice. We find that the momentum dependence of the vertex function depends strongly on doping: For large dopings it is rather weak, with decreasing doping it becomes more and more pronounced leading to a strong reduction of the vertex at larger momentum transfers until, for very small dopings, the vertex essentially consists of a forward scattering peak with a width proportional to the doping. This behavior occurs both in the t and the tt\ensuremath{'} models and also in one dimension where analytic expressions are derived. Correlation effects also change ${\mathrm{\ensuremath{\alpha}}}^{2}$F in general: The full-symmetric component ${\mathrm{\ensuremath{\alpha}}}^{2}$F1{\mathit{F}}_{1}F1​ is in the t model somewhat, in the tt\ensuremath{'} model rather strongly suppressed, especially near half-filling; the other symmetry components ${\mathrm{\ensuremath{\alpha}}}^{2}$Fi{\mathit{F}}_{\mathit{i}}Fi​ with i=2,...,5 increase strongly with decreasing doping and are of similar magnitude as ${\mathrm{\ensuremath{\alpha}}}^{2}$F1{\mathit{F}}_{1}F1​ near half-filling.Including also direct Coulomb repulsion nontrivial symmetries such as d wave become more stable than the s-wave order parameter below a critical value for the doping even for the considered phonon-mediated superconductivity. Most dramatic, however, is the quenching of the resistivity due to electron-phonon scattering both in the t and the tt\ensuremath{'} models at intermediate and small dopings. This result may explain the absence of phonon features in the experimental transport coefficients of high-Tc{\mathit{T}}_{\mathit{c}}Tc​ compounds. \textcopyright{} 1996 The American Physical Society.
Phase diagram of the Hubbard-Holstein model in the coexistence of electron-electron and electron-phonon interactions has been theoretically obtained with the density-matrix renormalization group method for one-dimensional systems, where an improved warm-up (the recursive sweep) procedure has enabled us to calculate various correlation functions. We have examined the cases of (i) the systems half-filled by electrons for the full parameter space spanned by the electron-electron and electron-phonon coupling constants and the phonon frequency, (ii) non-half-filled system, and (iii) trestle lattice. For (i), we have detected a region where both the charge and on-site pairing correlations decay with power laws in real space, which suggests a metallic behavior. While pairing correlations are not dominant in (i), we have found that they become dominant as the system is doped in (ii) or as the electronic band structure is modified (with a broken electron-hole symmetry) in (iii) in certain parameter regions.
Using a multistep renormalization group method, we study the low-temperature phases of the interacting one-dimensional (1D) electron gas coupled to phonons. We obtain analytic expressions for the weak-coupling quantum phase boundaries of the 1D extended Holstein-Hubbard model and the 1D extended Peierls-Hubbard model for general band-filling and phonon frequency. Away from half-filling, the phase diagrams are characterized by a delicate competition between spin density wave, charge density wave, and superconducting orders. We study the dependence of the ground state on the electron-phonon (el-ph) and electron-electron (el-el) coupling strengths, the screening length, electron bandwidth, phonon frequency, doping, and type of phonon. Unlike the case in Fermi liquids, in 1D the el-ph coupling is strongly renormalized, often to stronger values. Even when the bare phonon-induced attraction is weak compared to the bare el-el repulsion, a small amount of retardation can cause the renormalized el-ph interaction to dominate the problem. We find cases in which a repulsive el-el interaction enhances the superconducting susceptibility in the presence of a retarded el-ph interaction. The spin gap and superconducting susceptibility are found to be strongly dependent on the deviation from half-filling (doping). In some cases, the superconducting susceptibility varies nonmonotonically with doping and exhibits a maximum at a particular doping. For a quasi-1D array of weakly coupled, fluctuating 1D chains, the superconducting transition temperature Tc{T}_{c}Tc​ also exhibits a maximum as a function of doping. The effect of changing the ion mass (isotope effect) on Tc{T}_{c}Tc​ is found to be largest near half-filling and to decrease rapidly with doping.
With the discovery of strong-coupling physics and superconductivity in moir\'e superlattices, it is essential to have an understanding of strong-coupling driven superconductivity in systems with trigonal symmetry. The simplest lattice model with trigonal symmetry is the triangular lattice Hubbard model. Although the triangular lattice spin model is a heavily studied model in the context of frustration, studies of the hole-doped triangular lattice Hubbard model are rare. Here, we use density matrix renormalization group to investigate the dominant superconducting channels in the hole-doped triangular lattice Hubbard model over a range of repulsive interaction strengths. We find a clear transition from ppp-wave superconductivity at a moderate on-site repulsion strength (U/t=2U/t=2U/t=2) at filling above 1/4 ($n\ensuremath{\sim}0.65$) to ddd-wave superconductivity at a strong on-site repulsion strength (U/t=10U/t=10U/t=10) at filling below 1/4 ($n\ensuremath{\sim}0.4$). The unusual tunability that moir\'e superlattices offer in controlling U/tU/tU/t would open up the opportunity to realize this transition between ddd-wave and ppp-wave superconductivity.
We study the ground state of the doped Hubbard model on the honeycomb lattice in the small-doping and strongly interacting region. The nature of the ground state after doping holes into antiferromagnetic Mott insulating states on the honeycomb lattice remains a long-standing unsolved issue, even though tremendous effort has been spent investigating this challenging problem. In this work, we employ two complementary, state-of-the-art, many-body computational methods: constrained path (CP) auxiliary-field quantum Monte Carlo (AFQMC) with self-consistent constraint and density matrix renormalization group (DMRG) methods. Systematic and detailed cross validations are performed between these two methods for narrow systems for which DMRG can produce reliable results. AFQMC is then utilized to study wider systems to investigate the thermodynamic limit properties. The ground state is found to be a half-filled stripe state in the small-doping and strongly interacting region. The pairing correlation shows ddd-wave symmetry locally but decays exponentially with the distance between two pairs.
The relevance of the single-band two-dimensional Hubbard model to superconductivity in the doped cuprates has recently been questioned, based on density matrix renormalization group (DMRG) computations that found superconductivity over an unrealistically broad doping region upon electron-doping, yet a complete absence of superconductivity for hole-doping. We report very similar results from DMRG calculations on a Cu2O3{\mathrm{Cu}}_{2}{\mathrm{O}}_{3}Cu2​O3​ two-leg ladder within the parent three-band correlated-electron Hamiltonian. The strong asymmetry found in our calculations are in contradiction to the deep and profound symmetry in the experimental phase diagrams of electron- and hole-doped cuprate superconductors, as seen from the occurrence of quantum critical points within the superconducting domes in both cases that are characterized by Fermi surface reconstruction, large jumps in carrier density, and strange metal behavior.