  We present in detail the recently developed multi-configurational symmetrized
projector quantum Monte Carlo (MSPQMC) method for excited states of the Hubbard
model. We describe the implementation of the Monte Carlo method for a
multi-configurational trial wavefunction. We give a detailed discussion of
issues related to the symmetry of the projection procedure which validates our
Monte Carlo procedure for excited states and leads naturally to the idea of
symmetrized sampling for correlation functions, developed earlier in the
context of ground state simulations. It also leads to three possible averaging
schemes. We have analyzed the errors incurred in these various averaging
procedures and discuss and detail the preferred averaging procedure for
correlations that do not have the full symmetry of the Hamiltonian. We study
the energies and correlation functions of the low-lying excited states of the
half-filled Hubbard model in 1-D. We have used this technique to study the
pair-binding energies of two holes in $4n$ and $4n+2$ systems, which compare
well the Bethe ansatz data of Fye, Martins and Scalettar. We have also studied
small clusters amenable to exact diagonalization studies in 2-D and have
reproduced their energies and correlation functions by the MSPQMC method. We
identify two ways in which a multiconfigurational trial wavefunction can lead
to a negative sign problem. We observe that this effect is not severe in 1-D
and tends to vanish with increasing system size. We also note that this does
not enhance the severity of the sign problem in two dimensions.
