By using variational wave functions and quantum Monte Carlo techniques, we investigate the interplay between electron-electron and electron-phonon interactions in the two-dimensional Hubbard-Holstein model. Here, the ground-state phase diagram is triggered by several energy scales, i.e., the electron hopping
t, the on-site electron-electron interaction
U, the phonon energy ${\ensuremath{\omega}}_{0}$, and the electron-phonon coupling
g. At half filling, the ground state is an antiferromagnetic insulator for $U\ensuremath{\gtrsim}2{g}^{2}/{\ensuremath{\omega}}_{0}$, while it is a charge-density-wave (or bipolaronic) insulator for $U\ensuremath{\lesssim}2{g}^{2}/{\ensuremath{\omega}}_{0}$. In addition to these phases, we find a superconducting phase that intrudes between them. For ${\ensuremath{\omega}}_{0}/t=1$, superconductivity emerges when both
U/t and $2{g}^{2}/t{\ensuremath{\omega}}_{0}$ are small; then, by increasing the value of the phonon energy ${\ensuremath{\omega}}_{0}$, it extends along the transition line between antiferromagnetic and charge-density-wave insulators. Away from half filling, phase separation occurs when doping the charge-density-wave insulator, while a uniform (superconducting) ground state is found when doping the superconducting phase. In the analysis of finite-size effects, it is extremely important to average over twisted boundary conditions, especially in the weak-coupling limit and in the doped case.
Monte Carlo simulations are performed to examine superconductivity and charge-density-wave fluctuations in the infinite-dimensional electron-phonon Holstein model. The electron-phonon system interpolates between an attractive, static, Falicov-Kimball model that always exhibits charge-density-wave order and an attractive, instantaneous, Hubbard model that always superconducts as a function of phonon frequency. The maximum charge-density-wave transition temperature at half-filling is an order of magnitude smaller than the effective electronic bandwidth and is virtually independent of the phonon frequency. The maximum superconducting transition temperature depends strongly on phonon frequency and is bounded from above by the maximum charge-density-wave transition temperature. The exact solution is compared to both weak-coupling expansions and strong-coupling expansions. The effective phonon potential becomes anharmonic and develops a double-well structure that deepens as the electron-phonon interaction increases.
We study the Holstein-Hubbard model at half filling to explore ordered phases including superconductivity (SC), antiferromagnetism (AF), and charge order (CO) in situations where the electron-electron and electron-phonon interactions are strong (comparable to the electronic bandwidth). The model is solved in the dynamical mean-field approximation with a continuous-time quantum Monte Carlo impurity solver. We determine the superconducting transition temperature
Tc and the SC order parameter and show that the phonon-induced retardation or the strong Coulomb interaction leads to a significant reduction and shift of the
Tc dome against the effective electron-electron interaction
Ueff given by the Hubbard
U reduced by the phonon-mediated attraction in the static limit. This behavior is analyzed by comparison to an effective static model in the polaron representation with a renormalized bandwidth. In addition, we discuss the superconducting gap $\ensuremath{\Delta}$ and $2\ensuremath{\Delta}/{T}_{c}$ to reveal the effect of the retardation and the Coulomb interaction. We also determine the finite-temperature phase diagram including AF and CO. In the moderate-coupling regime, there is a hysteretic region of AF and CO around
Ueff=0, while the two phases are separated by a paramagnetic metal in the weak-coupling regime and a paramagnetic insulator in the strong-coupling regime.
We analyze the quantum phase diagram of the Holstein-Hubbard model using an asymptotically exact strong coupling expansion. We find all sorts of interesting phases including a pair-density wave, a charge 4e (and even a charge 6e) superconductor, regimes of phase separation, and a variety of distinct charge-density-wave, spin-density-wave, and superconducting regimes. We chart the crossovers that occur as a function of the degree of retardation, i.e., the ratio of characteristic phonon frequency to the strength of interactions.
The Holstein–Hubbard model is considered to study the competition between the electron-phonon interaction and the electron-electron interaction in a one-dimensional correlated Fermi systemat half-filling. First, a series of canonical transformations that incorporate phonon coherence and correlation effectsis performed and then the transformed Hamiltonian is averaged with respect to a zero-phonon state. The resulting effective electronic Hamiltonian is finally solved exactly by the Bethe-ansatz technique. It is shown that the present method suggests the existence of a wider intermediate metallic phase at the SDW-CDW crossover region than the one predicted by the simple variational method of Takada and Chatterjee. The presence of this metallic phase is indeed a favorable situation from the point of view of high temperature superconductivity.
We reconsider the strong-coupling expansion for the Hubbard model recently introduced by Sarker and Pairault et al. By introducing slave particles that act as projection operators onto the empty, singly occupied and doubly occupied atomic states, the perturbation theory around the atomic limit distinguishes between processes that do conserve or do not conserve the total number of doubly occupied sites. This allows for a systematic t/U expansion that does not break down at low temperature (t being the intersite hopping amplitude and U the local Coulomb repulsion). The fermionic field becomes a two-component field, which reflects the presence of the two Hubbard bands. The single-particle propagator is naturally expressed as a function of a 2×2 matrix self-energy. Furthermore, by introducing a time- and space-fluctuating spin-quantization axis in the functional integral, we can expand around a "non-degenerate" ground-state where each singly occupied site has a well defined spin direction (which may fluctuate in time). This formalism is used to derive the effective action of charge carriers in the lower Hubbard band to first order in t/U. We recover the action of the t–J model in the spin-hole coherent-state path integral. We also compare our results with those previously obtained by studying fluctuations around the large-U Hartree–Fock saddle point.