Computational solution of chemistry problems
AB initio quantum chemical techniques have been used to investigate weakly bound complexes of H2O and SO2. An energy gradient program was used to locate stable structures for the H2O, SO2 complexes, and SCF calculations were carried out to determine the binding energies of complexes with multiple water molecules. A 4-31G basis set was used for most potential energy searches. More accurate basis sets including a generally contracted basis set with d orbitals on the sulfur were used for geometry and binding energy verification. For single water complexes, five different stable geometries were located with binding energies between 4 and 11 Kcal mol(-1), suggesting a binding shell for H2O around SO2 and a mechanism for the formation of an SO2-containing water droplet. Calculations on one of the complexes utilizing a larger double zeta basis and d functions on the sulfur atom lead to adjusted binding energies in the range 3 to 8 Kcal mol(-1). Very little charge transfer between SO2 and H2O was present. Addition of more than one H2O was found to be energetically favorable although the addition of the fourth water in certain geometries did not increase the stability of the complex. An alternative mechanism for the tropospheric gas phase production of acid rain is suggested.