We use determinant quantum Monte Carlo to study the single-particle properties of quasiparticles and phonons in a variant of the two-dimensional Holstein model that includes an additional nonlinear electron-phonon (e-ph) interaction. We find that a small positive nonlinear interaction reduces the effective coupling between the electrons and the lattice, suppresses charge-density-wave (CDW) correlations, and hardens the effective phonon frequency. Conversely, a small negative nonlinear interaction can enhance the e-ph coupling resulting in heavier quasiparticles, an increased tendency towards a CDW phase at all fillings, and a softened phonon frequency. An effective linear model with a renormalized interaction strength and phonon frequency can qualitatively capture this physics; however, the quantitative effects of the nonlinearity on both the electronic and phononic degrees of freedom cannot be captured by such a model. These results are significant for typical nonlinear coupling strengths found in real materials, indicating that nonlinearity can have an important influence on the physics of many e-ph coupled systems.
Abstract Determining the range of validity of Migdal’s approximation for electron-phonon (e-ph) coupled systems is a long-standing problem. Many attempts to answer this question employ the Holstein Hamiltonian, where the electron density couples linearly to local lattice displacements. When these displacements are large, however, nonlinear corrections to the interaction must also be included, which can significantly alter the physical picture obtained from this model. Using determinant quantum Monte Carlo and the self-consistent Migdal approximation, we compared superconducting and charge-density-wave correlations in the Holstein model with and without second-order nonlinear interactions. We find a disagreement between the two cases, even for relatively small values of the e-ph coupling strength, and, importantly, that this can occur in the same parameter regions where Migdal’s approximation holds. Our results demonstrate that questions regarding the validity of Migdal’s approximation go hand in hand with questions of the validity of a linear e-ph interaction.
The electron-phonon $(e\text{\ensuremath{-}}\mathrm{ph})$ interaction remains of great interest in condensed matter physics and plays a vital role in realizing superconductors, charge density waves (CDW), and polarons. We study the two-dimensional Holstein model for $e\text{\ensuremath{-}}\mathrm{ph}$ coupling using determinant quantum Monte Carlo across a wide range of its phase diagram as a function of temperature, electron density, dimensionless $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength, and the adiabatic ratio of the phonon frequency to the Fermi energy. We describe the behavior of the CDW correlations, the competition between superconducting and CDW orders and polaron formation, the optimal conditions for superconductivity, and the transition from the weak-coupling regime to the strong-coupling regime. Superconductivity is optimized at intermediate $e\text{\ensuremath{-}}\mathrm{ph}$ coupling strength and intermediate electron density, and the superconducting correlations increase monotonically with phonon frequency. The global maximum for superconductivity in the Holstein model occurs at large phonon frequency, the limit where an attractive Hubbard model effectively describes the physics.
We study the effects of anharmonicity on the physics of the Holstein model, which describes the coupling of itinerant fermions and localized quantum phonons, by introducing a quartic term in the phonon potential energy. We find that the presence of this anharmonic term reduces the extent of the charge density wave (CDW) phase at half-filling as well as the transition temperature to this phase. Doping away from half-filling, we observe a first-order phase transition between the CDW and a homogeneous phase which is also present in the harmonic model. In addition, we study the evolution of the superconducting susceptibility in the doped region and show that anharmonicity can enhance the superconducting response.
Abstract Several state-of-the-art numerical methods have observed static or fluctuating spin and charge stripes in doped two-dimensional Hubbard models, suggesting that these orders play a significant role in shaping the cuprate phase diagram. Many experiments, however, also indicate that the cuprates have strong electron-phonon (e-ph) coupling, and it is unclear how this interaction influences stripe correlations. We study static and fluctuating stripe orders in the doped single-band Hubbard-Holstein model using zero temperature variational Monte Carlo and finite temperature determinant quantum Monte Carlo. We find that the lattice couples more strongly with the charge component of the stripes, leading to an enhancement or suppression of stripe correlations, depending on model parameters like the next-nearest-neighbor hopping
t′ t′ or phonon energy Ω. Our results help elucidate how thee-ph interaction can tip the delicate balance between stripe and superconducting correlations in the Hubbard-Holstein model with implications for our understanding of the high-Tccuprates.
Studies of Hamiltonians modeling electron-electron (e-e) and electron-phonon (e-ph) coupling have been fundamental in capturing the novel ordering seen in many quasi-one-dimensional condensed matter systems. Extending such studies to quasi-two-dimensional (2D) systems is of great current interest, as e-ph couplings are predicted to play a major role in the stabilization or enhancement of novel phases in 2D material systems. In this work, we study model systems that describe the interplay between the Hubbard coupling and the phonon modes in the Holstein and Su-Schrieffer-Heeger (SSH) Hamiltonians using the functional renormalization group. For both types of e-ph couplings, we find the predicted charge density wave phases in competition with antiferromagnetic ordering. As the system is doped, the transition shifts, with both orders showing incommensurate peaks. We compare the evolution of the quasiparticle weight for the Holstein model with that of the SSH model as the systems transition from antiferromagnetic to charge-ordered ground states. Finally, we calculate the self-energy of the phonon and capture the impact of charge ordering on the phonon modes.
The Holstein Hamiltonian describes fermions hopping on a lattice and interacting locally with dispersionless phonon degrees of freedom. In the low-density limit, dressed quasiparticles, polarons and bipolarons, propagate with an effective mass. At higher densities, pairs can condense into a low-temperature superconducting phase and, at or near commensurate filling on a bipartite lattice, to charge density wave (CDW) order. CDW formation breaks a discrete symmetry and hence occurs via a second-order (Ising) transition and therefore at a finite
Tcdw in two dimensions. Quantum Monte Carlo calculations have determined
Tcdw for a variety of geometries, including square, honeycomb, and Lieb lattices. The superconducting transition, on the other hand, in
d=2 is in the Kosterlitz-Thouless universality class and is much less well characterized. In this paper we determine
Tsc for the square lattice for several values of the density $\ensuremath{\rho}$ and phonon frequency ${\ensuremath{\omega}}_{0}$. We find that quasilong-range order sets in at ${T}_{\mathrm{sc}}\ensuremath{\lesssim}t/20$, where
t is the near-neighbor hopping amplitude, consistent with previous rough estimates from simulations which extrapolated to only the temperatures we reach from considerably higher
T. We also show evidence of a discontinuous evolution of the density as the CDW transition is approached at half filling.