To sum up, using an inner-sphere radius of 16 Å and a sample size of 250 snapshots with a 2 ps step size, our benchmark calculations indicate that the additional error from these parameters is minor compared to the basis-set incompleteness error. We estimate that the overall average error for the frequencies is approximately 4 cm -1 , with a maximum error of 16 cm -1 .