It is worth noting that these two properties (the dipole moment and the enthalpy of vapor-ization) are evaluated a posteriori and are not explicitly used to train the models during the optimization process. Moreover, both such observables are related to the extent of the intermolecular interactions, as the enthalpy of vaporization is proportional to the potential energy in the system, and the interaction between dipoles of the molecules plays a key role in it. Figure However, implementing such an approach encounters challenges due to the substantial differences between the geometric and electrostatic descriptions at the QM level compared to the all-atom models, as demonstrated by recent research on OPC3. An interesting outcome of these considerations is that while our approach can achieve optimal solutions at least as good as the state-of-the-art models in a very efficient way, it also underlines how such all atom water models are de facto a coarse grained description of the real water molecule features. Like other coarse-grained models, they encounter a degeneracy due to simplified representations of system degrees of freedom, resulting in a certain level of precision combined with inherent indeterminacy. This indeterminacy implies that different parametrizations can lead to similar behaviors. Notably, in our case, this degeneracy yields a set of non-identical solutions that belong to a minimum that is identified by the PCA.