No Arabic abstract
The Hubbard-Holstein model describes fermions on a discrete lattice, with on-site repulsion between fermions and a coupling to phonons that are localized on sites. Generally, at half-filling, increasing the coupling $g$ to the phonons drives the system towards a Peierls charge density wave state whereas increasing the electron-electron interaction $U$ drives the fermions into a Mott antiferromagnet. At low $g$ and $U$, or when doped, the system is metallic. In one-dimension, using quantum Monte Carlo simulations, we study the case where fermions have a long range coupling to phonons, with characteristic range $xi$, interpolating between the Holstein and Frohlich limits. Without electron-electron interaction, the fermions adopt a Peierls state when the coupling to the phonons is strong enough. This state is destabilized by a small coupling range $xi$, and leads to a collapse of the fermions, and, consequently, phase separation. Increasing interaction $U$ will drive any of these three phases (metallic, Peierls, phase separation) into a Mott insulator phase. The phase separation region is once again present in the $U e 0$ case, even for small values of the coupling range.
We study the interplay between the electron-phonon (e-ph) and on-site electron-electron (e-e) interactions in a three-orbital Hubbard-Holstein model on an extended one-dimensional lattice using determinant quantum Monte Carlo. For weak e-e and e-ph interactions, we observe a competition between an orbital-selective Mott phase (OSMP) and a (multicomponent) charge-density-wave (CDW) insulating phase, with an intermediate metallic phase located between them. For large e-e and e-ph couplings, the OSMP and CDW phases persist, while the metallic phase develops short-range orbital correlations and becomes insulating when both the e-e and e-ph interactions are large but comparable. Many of our conclusions are in line with those drawn from a prior dynamical mean field theory study of the two-orbital Hubbard-Holstein model [Phys. Rev. B 95, 12112(R) (2017)] in infinite dimension, suggesting that the competition between the e-ph and e-e interactions in multiorbital Hubbard-Holstein models leads to rich physics, regardless of the dimension of the system.
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 (PDW), a charge 4e (and even a charge 6e) superconductor, regimes of phase separation, and a variety of distinct charge-density-wave (CDW), spin-density-wave (SDW) 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.
We study the charge conductivity of the one-dimensional repulsive Hubbard model at finite temperature using the method of dynamical quantum typicality, focusing at half filling. This numerical approach allows us to obtain current autocorrelation functions from systems with as many as 18 sites, way beyond the range of standard exact diagonalization. Our data clearly suggest that the charge Drude weight vanishes with a power law as a function of system size. The low-frequency dependence of the conductivity is consistent with a finite dc value and thus with diffusion, despite large finite-size effects. Furthermore, we consider the mass-imbalanced Hubbard model for which the charge Drude weight decays exponentially with system size, as expected for a non-integrable model. We analyze the conductivity and diffusion constant as a function of the mass imbalance and we observe that the conductivity of the lighter component decreases exponentially fast with the mass-imbalance ratio. While in the extreme limit of immobile heavy particles, the Falicov-Kimball model, there is an effective Anderson-localization mechanism leading to a vanishing conductivity of the lighter species, we resolve finite conductivities for an inverse mass ratio of $eta gtrsim 0.25$.
Using the momentum average approximation we study the importance of adding higher-than-linear terms in the electron-phonon coupling on the properties of single polarons described by a generalized Holstein model. For medium and strong linear coupling, even small quadratic electron-phonon coupling terms are found to lead to very significant quantitative changes in the properties of the polaron, which cannot be captured by a linear Holstein Hamiltonian with renormalized parameters. We argue that the bi-polaron phase diagram is equally sensitive to addition of quadratic coupling terms if the linear coupling is large. These results suggest that the linear approximation is likely to be inappropriate to model systems with strong electron-phonon coupling, at least for low carrier concentrations.
We study finite-temperature transport properties of the one-dimensional Hubbard model using the density matrix renormalization group. Our aim is two-fold: First, we compute both the charge and the spin current correlation function of the integrable model at half filling. The former decays rapidly, implying that the corresponding Drude weight is either zero or very small. Second, we calculate the optical charge conductivity sigma(omega) in presence of small integrability-breaking next-nearest neighbor interactions (the extended Hubbard model). The DC conductivity is finite and diverges as the temperature is decreased below the gap. Our results thus suggest that the half-filled, gapped Hubbard model is a normal charge conductor at finite temperatures. As a testbed for our numerics, we compute sigma(omega) for the integrable XXZ spin chain in its gapped phase.