Application of the density matrix method to spectroscopy and dynamics of photosynthetic reaction centers