where θ t +t is the angle between bond vector at the times t and t . To calculate the spectral density, we fitted a sum of exponential functions with the large number, N , of pre-fixed timescales, τ i , to the correlation functions from simulations using the Python scipy optimize.nnls solver: