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 ΔH° and ΔG° 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 ΔH° and ΔG° 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.

Comparison of CBS-QB3, CBS-APNO, G2, and G3 thermochemical predictions with experiment for formation of ionic clusters of hydronium and hydroxide ions complexed with water

The GAUSSIAN 2, GAUSSIAN 3, complete basis set-QB3, and complete basis set-APNO methods have been used to calculate ΔH∘ and ΔG∘ values for ionic clusters of hydronium and hydroxide ions complexed with water. Results for the clusters H3O+(H2O)n andOH−(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.

Atmospheric Implications for Formation of Clusters of Ammonium and 1−10 Water Molecules

A mixed molecular dynamics/quantum mechanics model has been applied to the ammonium/water clustering system. The use of the high level MP2 calculation method and correlated basis sets, such as aug-cc-pVDZ and aug-cc-pVTZ, lends confidence in the accuracy of the extrapolated energies. These calculations provide electronic and free energies for the formation of clusters of ammonium and 1−10 water molecules at two different temperatures. Structures and thermodynamic values are in good agreement with previous experimental and theoretical results. The estimated concentration of these clusters in the troposphere was calculated using atmospheric amounts of ammonium and water. Results show the favorability of forming these clusters and implications for ion-induced nucleation in the atmosphere.

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 ΔH° and ΔG° values for neutral clusters of water, (H2O)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 × 1014 molecules/cm3, water trimer concentrations of 2.6 × 1012 molecules/cm3, tetramer concentrations of approximately 5.8 × 1011 molecules/cm3, and pentamer concentrations of approximately 3.5 × 1010 molecules/cm3 in saturated air at 298 K. These results have important implications for understanding the gas-phase chemistry of the lower troposphere.

The ability of the Gaussian-2, Gaussian-3, Complete Basis Set–QB3, and Complete Basis Set–APNO model chemistries to model the geometries of small water clusters

The Gaussian-2, Gaussian-3, Complete Basis Set-QB3, and Complete Basis Set-APNO methods have been used to calculate geometries of neutral clusters of water, (H2O)n, where n = 2–6. The structures are in excellent agreement with those determined from experiment and those predicted from previous high-level calculations. These methods also provide excellent thermochemical predictions for water clusters, and thus can be used with confidence in evaluating the structures and thermochemistry of water clusters.

Pople's Gaussian-3 model chemistry applied to an investigation of (H2O)8 water clusters

The Gaussian-3 method developed by Pople and coworkers has been used to calculate the free energy of neutral octamer clusters of water, (H2O)8. The most energetically stable structures are in excellent agreement with those determined from experiment and those predicted from previous high-level calculations. Cubic structures are favored over noncubic structures over all temperature ranges studied. The D2d cubic structure is the lowest free energy structure and dominates the potential energy and free energy hypersurfaces from 0 K to 298 K.

Accurate Predictions of Water Cluster Formation, (H2O)n=2-10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Article Accurate Experimental Values for the Free Energies of Hydration of H+, OH-, and H3O+

Accurate experimental values for the free energies of hydration, or the free energies of solvation, of the H+, OH-, and H3O+ ions are of fundamental importance. By use of the most accurate value for the free energy of solvation of H+, the known value for the free energy of solvation of water, and the known values for the gas phase and aqueous phase deprotonation of water, the corresponding experimental free energy of solvation for OH- is −106.4 ± 0.5 kcal/mol. Similarly, by use of the known values for ΔGf 0 for H3O, H2O+, and OH-, the known values for ΔGs for H+ and OH-, and the known value for the aqueous phase autoionization of water, we obtain an experimental free energy of solvation value for H3O+ of −103.4 ± 0.5 kcal/mol. These values are in excellent agreement with the commonly accepted values and with the value for ΔGs(OH-) obtained from embedding clusters of OH-(H2O)n in a dielectric continuum.

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 ΔG° 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.

Hydration of OCS with One to Four Water Molecules in Atmospheric and Laboratory Conditions

Carbonyl sulfide is the most abundant sulfur gas in the atmosphere. We have used MP2 and CCSD(T) theory to study the structures and thermochemistries of carbonyl sulfide interacting with one to four water molecules. We have completed an extensive search for clusters of OCS(H2O)n, where n = 1−4. We located three dimers, two trimers, five tetramers, and four pentamers with the MP2/aug-cc-pVDZ method. In each of the complexes with two or more waters, OCS preferentially interacts with low-energy water clusters. Our results match current theoretical and experimental literature, showing correlation with available geometries and frequencies for the OCS(H2O) species. The CCSD(T)/aug-cc-pVTZ thermochemical values combined with the average amount of OCS and the saturated concentration of H2O in the troposphere, lead to the prediction of 106 OCS(H2O) clusters·cm−3 and 102 OCS(H2O)2 clusters·cm−3 at 298 K. We predict the structures of OCS(H2O)n, n = 1−4 that should predominate in a low-temperature molecular beam and identify specific infrared vibrations that can be used to identify these different clusters.

Accurate Predictions of Water Cluster Formation, (H2O)n=2−10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.