Hydrated Complexes: Relevance to Atmospheric Chemistry and Climate

Global Search for Minimum Energy (H2O)n Clusters, n = 3−5

The Gaussian-3 (G3) model chemistry method has been used to calculate the relative deltaG(o) values for all possible conformers of neutral clusters of water, (H2O)n, where n = 3-5. A complete 12-fold conformational search around each hydrogen bond produced 144, 1728, and 20,736 initial starting structures of the water trimer, tetramer, and pentamer. These structures were optimized with PM3, followed by HF/6-31G* optimization, and then with the G3 model chemistry. Only two trimers are present on the G3 potential energy hypersurface. We identified 5 tetramers and 10 pentamers on the potential energy and free-energy hypersurfaces at 298 K. None of these 17 structures were linear; all linear starting models folded into cyclic or three-dimensional structures. The cyclic pentamer is the most stable isomer at 298 K. On the basis of this and previous studies, we expect the cyclic tetramers and pentamers to be the most significant cyclic water clusters in the atmosphere.

Comparison of Model Chemistry and Density Functional Theory Thermochemical Predictions with Experiment for Formation of Ionic Clusters of the Ammonium Cation Complexed with Water and Ammonia; Atmospheric Implications

The G2, G3, CBS-QB3, and CBS-APNO model chemistry methods and the B3LYP, B3P86, mPW1PW, and PBE1PBE density functional theory (DFT) methods have been used to calculate deltaH(o) and deltaG(o) values for ionic clusters of the ammonium ion complexed with water and ammonia. Results for the clusters NH4(+) (NH3)n and NH4(+) (H2O)n, where n = 1-4, are reported in this paper and compared against experimental values. Agreement with the experimental values for deltaH(o) and deltaG(o) for formation of NH4(+) (NH3)n clusters is excellent. Comparison between experiment and theory for formation of the NH4(+) (H2O)n clusters is quite good considering the uncertainty in the experimental values. The four DFT methods yield excellent agreement with experiment and the model chemistry methods when the aug-cc-pVTZ basis set is used for energetic calculations and the 6-31G* basis set is used for geometries and frequencies. On the basis of these results, we predict that all ions in the lower troposphere will be saturated with at least one complete first hydration shell of water molecules.

Thermodynamics of Forming Water Clusters at Various Temperatures and Pressures by Gaussian-2, Gaussian-3, Complete Basis Set-QB3, and Complete Basis Set-APNO Model Chemistries; Implications for Atmospheric Chemistry

The Gaussian-2, Gaussian-3, complete basis set- (CBS-) QB3, and CBS-APNO methods have been used to calculate Delta H degrees and Delta G degrees values for neutral clusters of water, (H(2)O)(n), where n = 2-6. The structures are similar to those determined from experiment and from previous high-level calculations. The thermodynamic calculations by the G2, G3, and CBS-APNO methods compare well against the estimated MP2(CBS) limit. The cyclic pentamer and hexamer structures release the most heat per hydrogen bond formed of any of the clusters. While the cage and prism forms of the hexamer are the lowest energy structures at very low temperatures, as temperature is increased the cyclic structure is favored. The free energies of cluster formation at different temperatures reveal interesting insights, the most striking being that the cyclic trimer, cyclic tetramer, and cyclic pentamer, like the dimer, should be detectable in the lower troposphere. We predict water dimer concentrations of 9 x 10(14) molecules/cm(3), water trimer concentrations of 2.6 x 10(12) molecules/cm(3), tetramer concentrations of approximately 5.8 x 10(11) molecules/cm(3), and pentamer concentrations of approximately 3.5 x 10(10) molecules/cm(3) in saturated air at 298 K. These results have important implications for understanding the gas-phase chemistry of the lower troposphere.