Rate equation analysis of three phonon-photon-phonon terahertz quantum cascade lasers