We find that in these systems, CNEO-MD can further improve over CNEO-DFT harmonic analysis in vibrational spectra calculations due to the incorporation of strong mode coupling effects on top of quantum nuclear delocalization effects. In contrast, conventional DFT-based AIMD may not necessarily perform better than DFT harmonic analysis because of the lack of quantum nuclear delocalization effects.