Polaronic approach to strongly correlated electron systems with strong electron-phonon interaction

The three-band $p\text{\ensuremath{-}}d$ model of strongly correlated electrons interacting with optical phonons via diagonal and off-diagonal electron-phonon interactions is considered within the cluster perturbation theory. In the beginning, the exact diagonalization of the Hamiltonian of a ${\mathrm{CuO}}_{4}$ cluster results in the construction of local polaronic eigenstates $\left|p\right\ensuremath{\rangle}$ with hole numbers ${n}_{h}=0,1,2$ per unit cell. The intercluster hoppings and interactions are exactly written in terms of Hubbard operators ${X}_{f}^{pq}=\left|p\right\ensuremath{\rangle}\ensuremath{\langle}q|$ determined via the multielectron polaronic eigenstates $\left|p\right\ensuremath{\rangle}$ at site $\mathbf{f}$. The Fermi-type single-electron quasiparticle dispersion and spectral weight are calculated for the undoped antiferromagnetic parent insulator like ${\mathrm{La}}_{2}{\mathrm{CuO}}_{4}$. The quasiparticle dispersion of Hubbard polarons is determined by a hybridization of the Hubbard fermion subbands with local Franck-Condon resonances so the main polaronic effect of the quasiparticle band structure is a splitting of the Hubbard bands on the number of Hubbard polaron subbands. Increasing of the EPI constant results in an increase of splitting, decrease of the subband width, transfer of the spectral weight to high-energy multiphonon excitations, and subsequent localization of the charge carriers. Herewith, the effect of such renormalization for the conduction band and the valence one differs depending on the ratio of the diagonal and off-diagonal EPI. In the framework of the GTB method, the Franck-Condon broadening of the spectral function of polaronic excitations is also reproduced for strongly correlated systems with strong electron-phonon interaction.