We employ the DPMD model trained with the SCAN functional that was introduced to construct the water phase diagram . It produces desirable properties in a wide range of thermodynamic conditions. The DW model implemented in DeePMD-kit estimates the dipole moment of water molecules of the DPMD model. Although the potential energy surface of the DPMD model is labeled with the meta-GGA SCAN functional, the GGA functional (BLYP) produces dipole moments nearly identical to those produced by the SCAN functional. To reduce the computational costs, we utilize the BLYP functional to generate data set for the DW model. Our decision to use the BLYP functional instead of the SCAN functional for calculating maximally localized Wannier function center positions is informed by the validation described in the Supporting Information. The ab initio calculations are performed with Vienna Ab initio Simulation Package (VASP) . More details about the quantum calculations are given in sec. 2.2.2.