Protonated water clusters described by an empirical valence bond potential.

The properties of low-lying stationary points on the potential energy surfaces of singly protonated water clusters (H(2)O)(n)H(+), are investigated using an empirical valence bond potential. Candidate global minima are reported for n=2-4, 8, and 20-22. For n=8, the variation in the energies and structures of low-lying minima with the number of valence bond states included in the model is studied. For n=4 and 8, disconnectivity graphs are also reported and are compared to results for the equivalent neutral water clusters as described by the rigid TIP3P potential. For the larger clusters, n=20-22, the structural properties of the low energy minima are compared with recently published spectroscopic data on these systems. The observed differences between the n=20 and n=21 systems are qualitatively reproduced by the model potential, but the similarities between the n=21 and n=22 systems are not.