In order to assign the spectra, the average positions of the atoms over the last 500 fs of the NVT trajectory were used as a starting point for geometry optimization calculations using static-DFT simulations with Cp2k. Upon complete minimization (within the constraint of a fixed unit cell), frequency calculations were performed within the harmonic approximation.