The Raman spectrum of diopside has been calculated by using three m>purelym> Density Functional Theory (DFT) Hamiltonians (PBE, WCPBE, LDA), the Hartree-Fock Hamiltonian (HF) and three hybrid HF/DFT ones (B3LYP, WC1LYP, PBE0). A comparison has been done between the calculated frequencies with those measured by Raman spectroscopy on a natural sample, along with several different orientations and beam polarizations, or retrieved from literature; such a comparison demonstrated the excellent performances of the hybrid Hamiltonians in reproducing the vibrational spectrum of the mineral, in line with what it is generally observed in literature concerning other mineral phases. In particular, the mean average absolute discrepancies of the calculated frequencies with respect to the experimental data were: 3.2 (WC1LYP), 4.7 (B3LYP), 6.5 (PBE0), 18.0 (PBE), 9.7 (WCPBE), 7.3 (LDA), and 40.6 cm−1 (HF). The very good quality of the WC1LYP results allowed for a reliable assignment of all of the experimentally observed Raman signals, and the corresponding assignments to specific patterns of atomic vibrational motion (normal modes).