The GAUSSIAN 2, GAUSSIAN 3, complete basis set-QB3, and complete basis set-APNO methods have been used to calculate DeltaH( composite function) and DeltaG( composite function) values for ionic clusters of hydronium and hydroxide ions complexed with water. Results for the clusters H3O+(H2O)n and OH-(H2O)n, where n=1-4 are reported in this paper, and compared against experimental values contained in the National Institutes of Standards and Technology (NIST) database. Agreement with experiment is excellent for the three ab initio methods for formation of these clusters. The high accuracy of these methods makes them reliable for calculating energetics for the formation of ionic clusters containing water. In addition this allows them to serve as a valuable check on the accuracy of experimental data reported in the NIST database, and makes them useful tools for addressing unresolved issues in atmospheric chemistry.