The harmonic oscillator model was used to justify the use of classical MD to calculate vibrational spectra since its classical autocorrelation functions match exactly with the quantum autocorrelation functions in peak positions. Our group also showed that in the harmonic oscillator model, the ground state CMES differs from the underlying PES only by a universal shift of 1 2 ℏω. Therefore, the classical dynamics picture remains the same, and CMES-MD is also exact for the harmonic oscillator. Here, to extend out CMES theory to excited states, we will analytically solve for the excited state CMESs for a harmonic oscillator and investigate their properties.