Symmetrized DMRG Method for Excited States of Hubbard Models

We extend the density matrix renormalization group method to exploit Parity, $C_2$ (rotation by $\pi$) and electron-hole symmtries of model Hamiltonians. We demonstrate the power of this method by obtaining the lowest energy states in all the eight symmetry subspaces of Hubbard chains with upto 50 sites. The ground-state energy, optical gap and spin gap of regular $U=4t$ and $U=6t$ Hubbard chains agree very well with exact results. This development extends the scope of the DMRG method and allows future applications to study of optical properties of low-dimensional conjugated polymeric systems.