subsequently to the desired level of 6-311++G(d,p)[36]/6-31+G(d). The choice of the thirty-six atoms that made use of the higher 6-311++G(d,p) basis sets was done from defining a reaction center -the set of molecules with the anticipated bond-breaking and bond formation and the immediately neighboring water molecules. This is the starting configuration of the pure-water model for dimerization.