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:cond-mat/9809321 · Physical review. B, Condensed matter

Numerical renormalization-group approach to fluctuation exchange in the presence of electron-phonon coupling: Pairing in the Holstein-Hubbard model

C.-H. Pao, H.-B. Schüttler

Published Jul 1, 1999• 5 citations
arXivHTMLPDF

Abstract

The fluctuation exchange (FLEX) approximation is applied to study the Holstein-Hubbard model. Due to the retarded nature of the phonon-mediated electron-electron interaction, neither fast Fourier transform nor previously developed numerical renormalization-group (NRG) methods for Hubbard-type purely electronic models are applicable, while brute force solutions are limited by the demands on computational time and storage which increase rapidly at low temperature T. Here, we describe a NRG technique to solve the FLEX equations efficiently. Several orders of magnitude of CPU time and storage can be saved at low T $(\ensuremath{\sim}80\mathrm{K})$. To test our approach, we compare our NRG results to brute force calculations on small lattices at elevated temperatures. Both s-wave and d-wave superconducting phase diagrams are then obtained by applying the NRG approach at low T. The isotope effect for s-wave pairing is BCS-like in a realistic phonon frequency range, but vanishes at unphysically large phonon frequency $(\ensuremath{\sim}\mathrm{b}\mathrm{a}\mathrm{n}\mathrm{d}\mathrm{w}\mathrm{i}\mathrm{d}\mathrm{t}\mathrm{h})$. For d-wave pairing, the isotope exponent is negative and small compared to the typical observed values in nonoptimally doped cuprates.

DOI
10.1103/physrevb.60.1283

Related papers

Renormalization Group Approach to Self-Consistent-Field Approximation for the Interacting Electrons

C.-H. Pao

Due to the retarded nature of the phonon-mediated electron-electron interaction, neither fast Fourier transform (FFT) nor previously developed NRG methods for Hubbard-type purely electronic models are applicable, while brute-force solutions are limited by the demands on computational time and storage which increase rapidly at low temperature T. A new numerical renormalization group (NRG) approach to the solution of self-consistent-field equations for correlated lattice electrons is illustrated for the two-dimensional Hubbard-type model. This approach allows us to handle the fermion frequency matrix inversions in the FLEX equations efficiently. Several orders of magnitude of CPU time and storage can be saved at the temperature scales relevant to electronic phase transitions.

Chinese Journal of Physics•Apr 1, 2000•0 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

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

Intertwined fluctuations and isotope effects in the Hubbard-Holstein model on the square lattice from functional renormalization

Aiman Al-Eryani, Sabine Andergassen, Michael M. Scherer

Electron-electron and electron-phonon interactions are responsible for the formation of spin, charge, and superconducting correlations in layered quantum materials. A paradigmatic model for such materials that captures both kinds of interactions is the two-dimensional Hubbard-Holstein model with a dispersionless Einstein phonon. In this work, we provide a detailed analysis of the magnetic, density, and superconducting fluctuations at and away from half filling. To that end, we employ the functional renormalization group using the recently introduced extension of the single-boson exchange formulation. More precisely, we go beyond previous approaches to the model by resolving the full frequency dependence of the two-particle vertex and taking into account the feedback from the electronic self-energy. We perform broad parameter scans in the space of Hubbard repulsion, electron-phonon coupling strength, and phonon frequency to explore the leading magnetic, density, and superconducting susceptibilities from the adiabatic to the anti-adiabatic regime. Our numerical data reveal that self-energy effects lead to an enhancement of the d -wave superconducting susceptibility toward larger phonon frequencies, in contrast to earlier isotope-effect studies. At small phonon frequencies, large density contributions to the s -wave superconducting susceptibility change sign and eventually lead to a reduction of s -wave superconductivity with increasing electron-phonon coupling, signaling the breakdown of Migdal-Eliashberg theory. We analyze our findings systematically, employing detailed diagnostics of the intertwined fluctuations and pinning down the various positive and negative isotope effects of the physical susceptibilities. Our results pave the way for a more comprehensive understanding of the interplay of electron-electron and electron-phonon interactions in correlated quantum materials.

arXiv:2504.10863•Physical Review Research•Aug 26, 2025•3 citations

Efficient fluctuation-exchange approach to low-temperature spin fluctuations and superconductivity: From the Hubbard model to Na x CoO 2 · y H 2 O

Niklas Witt, Erik G. C. P. van Loon, Takuya Nomoto +2 more

Superconductivity arises mostly at energy and temperature scales that are much smaller than the typical bare electronic energies. Since the computational effort of diagrammatic many-body techniques increases with the number of required Matsubara frequencies and thus with the inverse temperature, phase transitions that occur at low temperatures are typically hard to address numerically. In this work, we implement a fluctuation exchange (FLEX) approach to spin fluctuations and superconductivity using the ``intermediate representation basis'' (IR) [Shinaoka et al., Phys. Rev. B 96, 035147 (2017)] for Matsubara Green functions. This $\mathrm{FLEX}+\mathrm{IR}$ approach is numerically very efficient and enables us to reach temperatures on the order of ${10}^{\ensuremath{-}4}$ in units of the electronic bandwidth in multiorbital systems. After benchmarking the method in the doped repulsive Hubbard model on the square lattice, we study the possibility of spin-fluctuation-mediated superconductivity in the hydrated sodium cobalt material ${\mathrm{Na}}_{x}{\mathrm{CoO}}_{2}\ifmmode\cdot\else\textperiodcentered\fi{}{y\mathrm{H}}_{2}\mathrm{O}$ reaching the scale of the experimental transition temperature ${T}_{\mathrm{c}}=4.5\phantom{\rule{0.16em}{0ex}}\mathrm{K}$ and below.

arXiv:2012.04562•Physical review. B./Physical review. B•May 26, 2021•26 citations

Wilson's Numerical Renormalization Group Calculation of the Anderson–Holstein Model: Phonon Propagator

Han Yong Choi, Tae-Ho Park, Gun Sang Jeon

We study the symmetric Anderson–Holstein (AH) model at zero temperature with Wilson's numerical renormalization group (NRG) technique to study the interplay between the electron-electron and electron-phonon interactions. An improved method for calculating the phonon propagator using the NRG technique is presented and applied to the AH model as the onsite Coulomb repulsion U and electron-phonon coupling constant g are varied. As g is increased, the phonon propagator is successively renormalized, and for g≳gco crosses over to the regime where the mode splits into two components, one of which approaches back to the bare frequency and the other develops into a soft mode. The initial renormalization of the phonon mode, as g is increased from 0, can be either positive or negative depending on U and the hybridization Δ.

International Journal of Modern Physics B•Aug 10, 2003•1 citation
The fluctuation exchange (FLEX) approximation is applied to study the Holstein-Hubbard model. Due to the retarded nature of the phonon-mediated electron-electron interaction, neither fast Fourier transform nor previously developed numerical renormalization-group (NRG) methods for Hubbard-type purely electronic models are applicable, while brute force solutions are limited by the demands on computational time and storage which increase rapidly at low temperature T. Here, we describe a NRG technique to solve the FLEX equations efficiently. Several orders of magnitude of CPU time and storage can be saved at low T (\ensuremath∼80K)(\ensuremath{\sim}80\mathrm{K})(\ensuremath∼80K). To test our approach, we compare our NRG results to brute force calculations on small lattices at elevated temperatures. Both s-wave and d-wave superconducting phase diagrams are then obtained by applying the NRG approach at low T. The isotope effect for s-wave pairing is BCS-like in a realistic phonon frequency range, but vanishes at unphysically large phonon frequency (\ensuremath∼bandwidth)(\ensuremath{\sim}\mathrm{b}\mathrm{a}\mathrm{n}\mathrm{d}\mathrm{w}\mathrm{i}\mathrm{d}\mathrm{t}\mathrm{h})(\ensuremath∼bandwidth). For d-wave pairing, the isotope exponent is negative and small compared to the typical observed values in nonoptimally doped cuprates.
Due to the retarded nature of the phonon-mediated electron-electron interaction, neither fast Fourier transform (FFT) nor previously developed NRG methods for Hubbard-type purely electronic models are applicable, while brute-force solutions are limited by the demands on computational time and storage which increase rapidly at low temperature T. A new numerical renormalization group (NRG) approach to the solution of self-consistent-field equations for correlated lattice electrons is illustrated for the two-dimensional Hubbard-type model. This approach allows us to handle the fermion frequency matrix inversions in the FLEX equations efficiently. Several orders of magnitude of CPU time and storage can be saved at the temperature scales relevant to electronic phase transitions.
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 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.
Electron-electron and electron-phonon interactions are responsible for the formation of spin, charge, and superconducting correlations in layered quantum materials. A paradigmatic model for such materials that captures both kinds of interactions is the two-dimensional Hubbard-Holstein model with a dispersionless Einstein phonon. In this work, we provide a detailed analysis of the magnetic, density, and superconducting fluctuations at and away from half filling. To that end, we employ the functional renormalization group using the recently introduced extension of the single-boson exchange formulation. More precisely, we go beyond previous approaches to the model by resolving the full frequency dependence of the two-particle vertex and taking into account the feedback from the electronic self-energy. We perform broad parameter scans in the space of Hubbard repulsion, electron-phonon coupling strength, and phonon frequency to explore the leading magnetic, density, and superconducting susceptibilities from the adiabatic to the anti-adiabatic regime. Our numerical data reveal that self-energy effects lead to an enhancement of the d -wave superconducting susceptibility toward larger phonon frequencies, in contrast to earlier isotope-effect studies. At small phonon frequencies, large density contributions to the s -wave superconducting susceptibility change sign and eventually lead to a reduction of s -wave superconductivity with increasing electron-phonon coupling, signaling the breakdown of Migdal-Eliashberg theory. We analyze our findings systematically, employing detailed diagnostics of the intertwined fluctuations and pinning down the various positive and negative isotope effects of the physical susceptibilities. Our results pave the way for a more comprehensive understanding of the interplay of electron-electron and electron-phonon interactions in correlated quantum materials.
Superconductivity arises mostly at energy and temperature scales that are much smaller than the typical bare electronic energies. Since the computational effort of diagrammatic many-body techniques increases with the number of required Matsubara frequencies and thus with the inverse temperature, phase transitions that occur at low temperatures are typically hard to address numerically. In this work, we implement a fluctuation exchange (FLEX) approach to spin fluctuations and superconductivity using the ``intermediate representation basis'' (IR) [Shinaoka et al., Phys. Rev. B 96, 035147 (2017)] for Matsubara Green functions. This FLEX+IR\mathrm{FLEX}+\mathrm{IR}FLEX+IR approach is numerically very efficient and enables us to reach temperatures on the order of ${10}^{\ensuremath{-}4}$ in units of the electronic bandwidth in multiorbital systems. After benchmarking the method in the doped repulsive Hubbard model on the square lattice, we study the possibility of spin-fluctuation-mediated superconductivity in the hydrated sodium cobalt material ${\mathrm{Na}}_{x}{\mathrm{CoO}}_{2}\ifmmode\cdot\else\textperiodcentered\fi{}{y\mathrm{H}}_{2}\mathrm{O}$ reaching the scale of the experimental transition temperature Tc=4.5K{T}_{\mathrm{c}}=4.5\phantom{\rule{0.16em}{0ex}}\mathrm{K}Tc​=4.5K and below.
We study the symmetric Anderson–Holstein (AH) model at zero temperature with Wilson's numerical renormalization group (NRG) technique to study the interplay between the electron-electron and electron-phonon interactions. An improved method for calculating the phonon propagator using the NRG technique is presented and applied to the AH model as the onsite Coulomb repulsion U and electron-phonon coupling constant g are varied. As g is increased, the phonon propagator is successively renormalized, and for g≳gco crosses over to the regime where the mode splits into two components, one of which approaches back to the bare frequency and the other develops into a soft mode. The initial renormalization of the phonon mode, as g is increased from 0, can be either positive or negative depending on U and the hybridization Δ.