A theoretical study of water clusters: the relation between hydrogen-bond topology and interaction energy from quantum-chemical computations for clusters with up to 22 molecules

Quantum-chemical calculations of a variety of water clusters with eight, ten and twelve molecules were performed, as well as for selected clusters with up to 22 water molecules. Geometry optimizations were carried out at the B3LYP/cc-pVDZ level and single-point energies were calculated at the B3LYP/aug-cc-pVDZ level for selected clusters. The electronic energies were studied with respect to the geometry of the oxygen arrangement and six different characteristics of the hydrogen-bond arrangement in the cluster. Especially the effect of the placement of the non-hydrogen bonding hydrogens on the interaction energy was studied. Models for the interaction energy with respect to different characteristics of the hydrogen-bond arrangement were derived through least-square fits. The results from the study of the clusters with eight, ten and twelve molecules are used to predict possible low-energy structures for various shapes of clusters with up to 22 molecules.