Hydrogen-bonded clusters of hydrogen fluoride (HF)(n)(n = 2-15), have been studied at Hartree Fock level using the 6-31G basis set. Generated wavefunctions were used to calculate the electrostatic potential V(r) on the cluster surfaces. Calculated physical properties are expressed quantitatively in terms of molecular surface area. volume and three statistically based quantities obtained from the surface electrostatic potential. The maxima and minima in the potentials V-S.min(r), V-S.max(r) have predicted the sites and directionality of the H bonding in the cluster, The energy of the HOMO of the clusters correlates linearly with its V-S.min. Energy gaps between frontier orbitals, epsilon(HOMO) (LUMO) have been discussed as a measure of cluster reactivity. (C) 2002 Elsevier Science B.V. All rights reserved.