Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules

The force constants and equilibrium geometry of water were calculated ab initio from Hartree-Fock self-consistent field wavefunctions using the force method, and a 53/31, a 73/3 + 1 and a 95/41 + 2 gaussian lobe basis set. The force method proved to be very economic and numerically accurate.

The calculated values agree well with the experimental ones. Especially good is the agreement for the interaction constant F_rα . Calculated force constants for the 95/41 + 2 basis are: F_r = 9·158 mdyn/å, F_α = 0·8513 mdyn/å, F_rα = +0·3007, F_rr = -0·1724; the experimental values from isotope frequencies are: F_r = 8·456, F_α = 0·762, F_rα = +0·246, F_rr = -0·100 mdyn/å. The agreement in the geometry is also good for the 95/41 + 2 basis r _e = 0·9518 å, α = 107·49°, and the experimental r _e = 0·9572 å, α = 104·52°.

Force constants were found to be quite insensitive to variation in the basis functions, all three sets giving almost the same result. On the other hand, the dipole moment derivative is very sensitive.