The two-dimensional Holstein-Hubbard model is studied by means of continuous-time quantum Monte Carlo simulations. Using renormalization-group-invariant correlation ratios and finite-size extrapolation, the critical temperature of the charge-density-wave transition is determined as a function of coupling strength, phonon frequency, and Hubbard repulsion. The phase transition is demonstrated to be in the universality class of the two-dimensional Ising model and detectable via the fidelity susceptibility. The structure of the ground-state phase diagram and the possibility of a bipolaronic metal with a single-particle gap above ${T}_{c}$ are explored.
The two-dimensional Holstein-Hubbard model is studied by means of continuous-time quantum Monte Carlo simulations. Using renormalization-group-invariant correlation ratios and finite-size extrapolation, the critical temperature of the charge-density-wave transition is determined as a function of coupling strength, phonon frequency, and Hubbard repulsion. The phase transition is demonstrated to be in the universality class of the two-dimensional Ising model and detectable via the fidelity susceptibility. The structure of the ground-state phase diagram and the possibility of a bipolaronic metal with a single-particle gap above Tc are explored.
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.
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 dynamical properties of the two-dimensional square lattice Holstein Hamiltonian using numerically exact determinant quantum Monte Carlo simulations. In particular, we report systematic calculations of the model's single-particle spectral function and optical conductivity over a range of phonon energies, electron-phonon coupling strengths, carrier concentrations, and temperatures. In doing so, we map the evolution of the system from a dressed metallic phase to a (bi)polaron liquid/insulator to a charge-density-wave insulator as the carrier concentration is tuned from dilute values to half-filling. We corroborate these results using several equal-time correlation measurements and related proxy quantities reconstructed from time-displaced correlation measurements. This paper is also accompanied by an extensive open data set, covering over 16,640 unique parameter values.
We present a diagrammatic Monte Carlo study of the properties of the Hubbard-Holstein bipolaron on a two-dimensional square lattice. With a small Coulomb repulsion U and with increasing electron-phonon interaction, and when reaching a value about two times smaller than the one corresponding to the transition of light polaron to heavy polaron, the system suffers a sharp transition from a state formed by two weakly bound light polarons to a heavy, strongly bound on-site bipolaron. Aside from this rather conventional bipolaron a new bipolaron state is found for large U at intermediate and large electron-phonon coupling, corresponding to two polarons bound on nearest-neighbor sites. We discuss both the properties of the different bipolaron states and the transition from one state to another. We present a phase diagram in parameter space defined by the electron-phonon coupling and U. Our numerical method does not use any artificial approximation and can be easily modified to other bipolaron models with longer range electron-phonon and∕or electron-electron interaction.
We use determinant quantum Monte Carlo to study the half-filled `bond-Holstein' model on a square lattice. We find that the model exhibits a charge-density-wave (CDW) phase transition with a critical temperature Tcdw considerably higher than that of the canonical `site-Holstein' model. Using a finite-size scaling analysis of the charge structure factor Scdw, we obtain Tcdw to greater than one percent accuracy. At the same time, local observables also show clear signatures consistent with the transition temperatures inferred from our scaling analysis. We attribute the enhanced CDW tendencies to a phonon-mediated nearest-neighbor electron repulsion that is directly proportional to the dimensionless electron-phonon coupling λ in the atomic (t→0) limit. This behavior contrasts with the site-Holstein case, where the same limit yields only an on-site attraction. We supplement our analysis with results from several unsupervised machine learning methods, which not only confirm our estimates of Tcdw but also provide insight into the high-temperature crossover between a metallic and bipolaron liquid regime.
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.