In this work, we derived and implemented the QM/DIM method and its high-order energy derivatives for the calculation of noble MNP-mediated molecular vibrational spectra, including IR and Raman scattering. Taking the complexes composed of a pyridine molecule with MNPs consisting of 18 or 32 gold atoms as testing examples, we assessed the accuracy of the current QM/DIM scheme in the calculations of molecular vibrational frequencies and IR/Raman intensities. The QM/DIM method demonstrates the ability to reproduce the molecular vibrational frequencies and IR spectral intensities obtained using the full QM approach. The Raman profiles from QM/DIM calculations behave well for large intermolecular distances such as 5 Å for the tested complexes, while it shows artificial enhancements at shorter molecule-MNP distances for some of the vibrational modes. Even though, with large MNPs such as Au 2057 and Ag 2057 , the simulated normal SERS of 4,4 -BPy by QM/DIM method are comparable with experimental spectra.