Currently, the extended Hubbard model, which is a good approximation only for narrow-band system, suggests that the width of polaron in polymer is heavily reduced by the electron-electron interaction. However, the polymers possess wide band, so the off-diagonal correlation cannot be neglected. We considered the effects of full electron-electron correlation on the width and the creation energy of polaron in polymer. In the case of strong screening, our results are greatly different from those of Hubbard model. This suggests that in polymer, the Hubbard model should be gone beyond.
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
CuO4 cluster results in the construction of local polaronic eigenstates $\left|p\right\ensuremath{\rangle}$ with hole numbers
nh=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
f. The Fermi-type single-electron quasiparticle dispersion and spectral weight are calculated for the undoped antiferromagnetic parent insulator like
La2CuO4. 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.