The workflow covers steps essential for generating structures, calculating properties, and post-processing the properties to the final IR and Raman spectra. This includes configurational sampling through an MD simulation, calculating analytically the partial Hessian and the dipole and polarizability gradients, which are then post-processed to obtain frequencies and IR and Raman intensities for the different configurations.