The inverse scattering method is applied to a coupled electron-phonon system of one dimension.Adiabatic potentials are calculated exactly for the case of any number of independent electrons coupled to deformable continuum.Polarons, bipolarons and their interaction are discussed.
Characterizing bipolaron binding, and understanding how it depends on electron-phonon interaction, is crucial to unraveling the nature of emergent many-body states in strongly interacting electron-phonon systems. So far, most studies of bipolarons have been limited to the Holstein model, in which the coupling constant is momentum-independent. The paradigmatic example of momentum-dependent electron-phonon interaction comes from the system in which phonon distortions modify electron hopping, the SSH model. Already individual polarons in the SSH model are richer than the Holstein model counterparts, and feature a phase transition into the finite momentum ground state with increasing electron-phonon interaction. In this paper, we use a variational approach to study bipolarons in the one-dimensional SSH model and discuss their ground state, dispersion, and excitation spectra. We explore the full parameter range of the system, including the adiabatic regime of slow phonons, which was inaccessible to previous theoretical studies. In agreement with earlier studies, we find that in the anti-adiabatic strongly interacting regime, bipolarons have low effective mass. By contrast, in the adiabatic case, we find that increasing electron-phonon interactions results in an exponential increase of the bipolaron mass. We establish the existence of multiple branches of bound excited states of SSH bipolaron and discuss the signatures of these bound states in dynamics. We show that in the anti-adiabatic regime, response functions obey a parity selection rule, that imposes symmetry constraints on the excitation spectra and provides a clear signature of SSH bipolarons.
The Holstein Hubbard and Holstein t-J models are studied for a wide range of phonon frequencies, electron-electron, and electron-phonon interaction strengths on finite lattices with up to ten sites by means of direct Lanczos diagonalization. Previously the necessary truncation of the phononic Hilbert space caused serious limitations to either very small systems (four or even two sites) or to weak electron-phonon coupling, in particular in the adiabatic regime. Using parallel computers we were able to investigate the transition from ``large'' to ``small'' polarons in detail. By resolving the low-lying eigenstates of the Hamiltonian and by calculating the spectral function, we can identify a polaron band in the strong-coupling case, whose dispersion deviates from the free-particle dispersion at low and intermediate phonon frequencies. For two electrons (holes) we establish the existence of bipolaronic states and discuss the formation of a bipolaron band. For the two-dimensional Holstein t-J model, we demonstrate that the formation of hole polarons is favored by strong Coulomb correlations. Analyzing hole-hole correlation functions, we find that hole binding is enhanced as a dynamical effect of the electron-phonon interaction. \textcopyright{} 1996 The American Physical Society.
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.