These excited state CMESs can be used in an ad hoc manner to perform MD simulations. The underlying assumption for these simulations is that the nuclear wave functions adiabatically keep their excited state character during the dynamics and do not relax to the lower vibrational states. In the next section we will show that we can use them to obtain vibrational excited state absorption frequencies. For these frequency calculations, a natural question is that since we need to construct excited state CMESs, why not we directly take the energy difference between CMESs to obtain the frequencies? We note that because of the requirement on the orthogonality with the lower states of the same expectation position, the "excited states" for CMESs do not rigorously correspond to the vibrational excited states solved from the Schrödinger equation, therefore, the energy gaps between CMESs do not necessarily match well with the reference excited state absorption values. Additionally, intensity information will be more directly available from dynamics simulations than pure static energy calculation. Furthermore, the success of vibrational frequency calculations will demonstrate the good quality, at least near the equilibrium region, of the vibrational excited state CMESs, which we can use in the future to perform excited state dynamics simulations.