Rate equation description of multi-phonon creation of color centers and simultaneous one-phonon annihilation