The electron-phonon interaction is one of the fundamental interactions in almost all condensed matter materials. In conventional superconductors, the electron-phonon interaction is the glue that attracts two electrons to one another to form a pair. A strong electron-phonon interaction leads to the concept of a polaron, which is an electron with lattice distortions around it. The small polaron is a polaron with spatial extent comparable to an interatomic dimension of the solid. Evidence for polarons has been identified in many experiments in superconductors and semiconductors. In this thesis we present exact calculations of the polaron. Specifically we have refined Trugman's method to solve the ground state of an electron-phonon coupled system in the whole parameter regime, and we also generalized this method to treat spin-orbit coupled systems. The most difficult regimes, which is the strong-coupling regime and the small phonon frequency limit, have been solved by these refinements. There are three representative kinds of electron-phonon interaction, the Holstein model, the Fr"ohlich model and the BLF-SSH model. In this thesis we have addressed the first and the third one. The second one, the Fr"ohlich model, is very similar to the Holstein model but the interaction is nonlocal. For the Holstein model we have observed the expected smooth crossover from free electrons to small polarons, while for the BLF-SSH model, we have studied the weak coupling regime with perturbation theory and derived a new analytical result for the one-dimensional problem.
Theoretical studies suggest that Su-Schrieffer-Heeger-like electron-phonon (
e-ph) interactions can mediate high-temperature bipolaronic superconductivity that is robust against repulsive electron-electron interactions. Here we present a comparative analysis of the pairing and competing charge/bond correlations in the two-dimensional Holstein and optical Su-Schrieffer-Heeger (SSH) models using numerically exact determinant quantum Monte Carlo. We find that the SSH interactions support light bipolarons and strong superconducting correlations out to relatively large values of the
e-ph coupling
λ and densities near half-filling, while the Holstein interaction does not due to the formation of heavy bipolarons and competing charge-density-wave order. We further find that the Holstein and SSH models have comparable pairing correlations in the weak coupling limit for carrier concentrations
⟨n⟩≪1, where competing orders and polaronic effects are absent. These results support the proposal that SSH (bi)polarons can support superconductivity to larger values of
λ in comparison to the Holstein polaron, but that the resulting
Tc gains are small in the weak coupling limit. We also find that the SSH model's pairing correlations are suppressed after including a weak on-site Hubbard repulsion. These results have important implications for identifying and engineering bipolaronic superconductivity.
The bound states of two electrons in the adiabatic Holstein-Hubbard model are studied numerically in one and two dimensions from the anticontinuous limit. This model involves a competition between a local electron-phonon coupling (with a classical lattice) which tends to form pairs of electrons and the repulsive Hubbard interaction
U≥0 which tends to break them. In 1D, the ground-state always consists in a pair of localized polarons in a singlet state. They are located at the same site for U=0. Increasing U, there is a first order transition at which the bipolaron becomes a spin singlet pair of two polarons bounded by a magnetic interaction. The pinning mode of the bipolaron soften in the vicinity of this transition leading to a higher mobility of the bipolaron which is tested numerically. In 2D, and for any
U, the electron-phonon coupling needs to be large enough in order to form small polarons or bipolarons instead of extended electrons. We calculate the phase diagram of the bipolaron involving first order transitions lines with a triple point. A pair of polarons can form three types of bipolarons: a) on a single site at small
U, b) a spin singlet state on two nearest neighbor sites for larger
U as in 1D and c) a new intermediate state obtained as the resonant combination of four 2-sites singlet states sharing a central site, called quadrisinglet. The breathing and pinning internal modes of bipolarons in 2D generally only weakly soften and thus, they are practically not mobile. On the opposite, in the vicinity of the triple point involving the quadrisinglet, both modes exhibit a significant softening. However, it was not sufficient for allowing the existence of a classical mobile bipolaron (at least in that model).