Fit a general data set (or specially the phonon spectral energy density) as a Lorentzian function to obtain the peak position (or phonon frequency) and full width at half maximum (or relaxation time).