Unlike the primary relaxation mode, the characteristic time (τ 2 ) and the spectral amplitude (∆ϵ 2 ) of the secondary relaxation are significantly different between the models. Indeed, the location and physical interpretation of the secondary relaxation mode has been the subject of ongoing debate . The characteristic time τ 2 reported by earlier works ranges from 0.18 to 1.43 ps, and the amplitude is not negligible (1.60 < ∆ϵ 2 < 2.50) . Our simulation results suggest that the DPMD model again outperforms the pairwise potentials for predicting both the spectral amplitude and the characteristic time.