  Combining density-matrix and Lanczos algorithms we propose a new optimized
phonon approach for finite-cluster diagonalizations of interacting
electron-phonon systems. To illustrate the efficiency and reliability of our
method, we investigate the problem of bipolaron band formation in the extended
Holstein Hubbard model.
