In the following, we summarize the key equations for the MD-based IR and Raman intensities, together with the different methodologies we tested for the calculation of dipole moments and polarizability tensors. The expressions for IR and Raman intensities for static calculations based on conventional normal mode analysis are given in Section S1 in the Supplementary Information (SI).