The QM/MM AIMD simulations were initiated from 40 structures taken from the equilibrium MD trajectory. The QM region included all waters within 6.5 Å radius from the central water (this selection criterion resulted in the 35-40 QM water molecules, depending on a snapshot). The QM part was treated by ωB97X-D/6-31G* and the MM part by TIP3P.