Ab initio MP2/aug'-cc-pVTZ calculations have been carried out to determine the structures and binding energies of complexes of BeF 2 , BeF 3 -1 , and BeF 4 -2 with n H 2 O molecules, with n = 1-6. For each series of complexes with BeF 2 , BeF 3 -1 , and BeF 4 -2 , the binding energies increase as the number of water molecules increases, but the binding energies per intermolecular bond for complexes formed between BeF 2 and BeF 4 -2 and nH 2 O decrease as the number of water molecules increases. The binding energies per bond of the BeF 3 -1 complexes show little dependence on the number of H 2 O molecules present. Intermolecular O-H ... F hydrogen bonds (HB) stabilise all complexes of H 2 O with the beryllium bases except for BeF 2 :1H 2 O and BeF 2 :2H 2 O. Complexes of BeF 2 are also stabilised by O ... Be beryllium bonds (BeB) and O-H ... O hydrogen bonds. EOM-CCSD calculations have also been performed to determine the spin-spin coupling constants. The one-bond coupling constants 1 J(Be-F) increase as the Be-F distance decreases, and exhibit an excellent second-order correlation with that distance. 2h J(O-F) coupling constants across O-H ... F hydrogen bonds also exhibit a secondorder dependence on distance. Coupling constants 1be J(O-Be) and 2h J(O-O) are found in complexes with BeF 2 and exhibit a linear dependence on the Be-O and O-O distances, respectively.