Zero-temperature phases of the two-dimensional Hubbard-Holstein model: A non-Gaussian exact diagonalization study

We propose a numerical method which embeds the variational non-Gaussian wavefunction approach within exact diagonalization, allowing for efficient treatment of correlated systems with both electron-electron and electron-phonon interactions. Using a generalized polaron transformation, we construct a variational wavefunction that absorbs entanglement between electrons and phonons into a variational non-Gaussian transformation; exact...