Note that here we only ran 5 trajectories for NVE -eqp and 15 for NVE -nms methods, but 30 trajectories for the NVT-NVE scheme. This is because the spectra calculated from the direct NVE schemes need fewer trajectories to converge. Figure shows the convergence test results for power spectra as a function of number of trajectories used in different methods.