Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)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.

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 H2SO4(H2O)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 (HSO4–·H3O+)(H2O)n−1clusters are competitive with the neutral (H2SO4)(H2O)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 H2SO4(H2O)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 H2SO4–H2O system must contain more than one H2SO4 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.

Characterization of Halyomorpha halys (brown marmorated stink bug) biogenic volatile organic compound emissions and their role in secondary organic aerosol formation.

The formation of aerosols is a key component in understanding cloud formation in the context of radiative forcings and global climate modeling. Biogenic volatile organic compounds (BVOCs) are a significant source of aerosols, yet there is still much to be learned about their structures, sources, and interactions. The aims of this project were to identify the BVOCs found in the defense chemicals of the brown marmorated stink bug Halymorpha halys and quantify them using gas chromatography-mass spectrometry (GC/MS) and test whether oxidation of these compounds by ozone-promoted aerosol and cloud seed formation. The bugs were tested under two conditions: agitation by asphyxiation and direct glandular exposure. Tridecane, 2(5H)-furanone 5-ethyl, and (E)-2-decenal were identified as the three most abundant compounds. H. halys were also tested in the agitated condition in a smog chamber. It was found that in the presence of 100-180 ppm ozone, secondary aerosols do form. A scanning mobility particle sizer (SMPS) and a cloud condensation nuclei counter (CCNC) were used to characterize the secondary aerosols that formed. This reaction resulted in 0.23 microg/ bug of particulate mass. It was also found that these secondary organic aerosol particles could act as cloud condensation nuclei. At a supersaturation of 1%, we found a kappa value of 0.09. Once regional populations of these stink bugs stabilize and the populations estimates can be made, the additional impacts of their contribution to regional air quality can be calculated.