We investigate the quantum dynamics of two defect centers in solids,which are coupled by vacuum-induced dipole-dipole interactions.When the interaction between defects and phonons is taken into account,the two coupled electron-phonon systems make up two equivalent multilevel atoms.By making Born-Markov and rotating wave approximations,we derive a master equation describing the dynamics of the coupled multilevel atoms.The results indicate the concepts of subradiant and superradiant states can be applied to these systems and the population transfer process presents different behaviors from those of the two dipolar-coupled two-level atoms due to the participation of phonons.