Lammps is used to perform MD simaultions of classical pairwise potentials(SPC/Fw and TIP4P/ε) and DPMD models. All systems each containing 556 water molecules are equilibrated during short NVT simulations followed by a 5 ns NPT simulation, using the Nosè-Hoover barostat with a damping parameter of τ P = 1.0 ps. Every simulation used the Nosè-Hoover thermostat with a damping parameter of τ T = 0.1 ps and a timestep of 1.0 fs. After the equilibration run, we perform 30 ns long production run using the same setup. During the run, the polarization vectors are sampled every 10 fs. Subsequently, to calculate the vibrational density of states (vDOS), additional simulations lasting 100 ps are conducted after the main production run. The trajectory data are collected every 1.0 fs for these simulations.