Water network-mediated, electron-induced proton transfer in [C5H5N ⋠(H2O)n](-) clusters.

The role of proton-assisted charge accommodation in electron capture by a heterocyclic electron scavenger is investigated through theoretical analysis of the vibrational spectra of cold, gas phase [Py ⋅ (H2O)n=3-5](-) clusters. These radical anions are formed when an excess electron is attached to water clusters containing a single pyridine (Py) molecule in a supersonic jet ion source. Under these conditions, the cluster ion distribution starts promptly at n = 3, and the photoelectron spectra, combined with vibrational predissociation spectra of the Ar-tagged anions, establish that for n > 3, these species are best described as hydrated hydroxide ions with the neutral pyridinium radical, PyH((0)), occupying one of the primary solvation sites of the OH(-). The n = 3 cluster appears to be a special case where charge localization on Py and hydroxide is nearly isoenergetic, and the nature of this species is explored with ab initio molecular dynamics calculations of the trajectories that start from metastable arrangements of the anion based on a diffuse, essentially dipole-bound electron. These calculations indicate that the reaction proceeds via a relatively slow rearrangement of the water network to create a favorable hydration configuration around the water molecule that eventually donates a proton to the Py nitrogen atom to yield the product hydroxide ion. The correlation between the degree of excess charge localization and the evolving shape of the water network revealed by this approach thus provides a microscopic picture of the "solvent coordinate" at the heart of a prototypical proton-coupled electron transfer reaction.

Quantum mechanical study of sulfuric acid hydration: atmospheric implications.

The role of the binary nucleation of sulfuric acid in aerosol formation and its implications for global warming is one of the fundamental unsettled questions in atmospheric chemistry. We have investigated the thermodynamics of sulfuric acid hydration using ab initio quantum mechanical methods. For H(2)SO(4)(H(2)O)(n) where n = 1-6, we used a scheme combining molecular dynamics configurational sampling with high-level ab initio calculations to locate the global and many low lying local minima for each cluster size. For each isomer, we extrapolated the Møller-Plesset perturbation theory (MP2) energies to their complete basis set (CBS) limit and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled harmonic vibrational frequencies. We found that ionic pair (HSO(4)(-)·H(3)O(+))(H(2)O)(n-1) clusters are competitive with the neutral (H(2)SO(4))(H(2)O)(n) clusters for n ≥ 3 and are more stable than neutral clusters for n ≥ 4 depending on the temperature. The Boltzmann averaged Gibbs free energies for the formation of H(2)SO(4)(H(2)O)(n) clusters are favorable in colder regions of the troposphere (T = 216.65-273.15 K) for n = 1-6, but the formation of clusters with n ≥ 5 is not favorable at higher (T > 273.15 K) temperatures. Our results suggest the critical cluster of a binary H(2)SO(4)-H(2)O system must contain more than one H(2)SO(4) and are in concert with recent findings (1) that the role of binary nucleation is small at ambient conditions, but significant at colder regions of the troposphere. Overall, the results support the idea that binary nucleation of sulfuric acid and water cannot account for nucleation of sulfuric acid in the lower troposphere.

Benchmark structures and binding energies of small water clusters with anharmonicity corrections.

For (H(2)O)(n) where n = 1-10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 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 (H(2)O)(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.