Prediction of Accurate Anharmonic Experimental Vibrational Frequencies for Water Clusters, (H2O)n, n = 2−5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2−5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.

The Role of Anharmonicity in Hydrogen-Bonded Systems: The Case of Water Clusters

The nature of vibrational anharmonicity has been examined for the case of small water clusters using second-order vibrational perturbation theory (VPT2) applied on second-order Møller–Plesset perturbation theory (MP2) potential energy surfaces. Using a training set of 16 water clusters (H2O)n=2–6,8,9 with a total of 723 vibrational modes, we determined scaling factors that map the harmonic frequencies onto anharmonic ones. The intermolecular modes were found to be substantially more anharmonic than intramolecular bending and stretching modes. Due to the varying levels of anharmonicity of the intermolecular and intramolecular modes, different frequency scaling factors for each region were necessary to achieve the highest accuracy. Furthermore, new scaling factors for zero-point vibrational energies (ZPVE) and vibrational corrections to the enthalpy (ΔHvib) and the entropy (Svib) have been determined. All the scaling factors reported in this study are different from previous works in that they are intended for hydrogen-bonded systems, while others were built using experimental frequencies of covalently bonded systems. An application of our scaling factors to the vibrational frequencies of water dimer and thermodynamic functions of 11 larger water clusters highlights the importance of anharmonic effects in hydrogen-bonded systems.

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.

The Role of Anharmonicity in Hydrogen-Bonded Systems: The Case of Water Clusters

The nature of vibrational anharmonicity has been examined for the case of small water clusters using second-order vibrational perturbation theory (VPT2) applied on second-order Møller–Plesset perturbation theory (MP2) potential energy surfaces. Using a training set of 16 water clusters (H2O)n=2–6,8,9 with a total of 723 vibrational modes, we determined scaling factors that map the harmonic frequencies onto anharmonic ones. The intermolecular modes were found to be substantially more anharmonic than intramolecular bending and stretching modes. Due to the varying levels of anharmonicity of the intermolecular and intramolecular modes, different frequency scaling factors for each region were necessary to achieve the highest accuracy. Furthermore, new scaling factors for zero-point vibrational energies (ZPVE) and vibrational corrections to the enthalpy (ΔHvib) and the entropy (Svib) have been determined. All the scaling factors reported in this study are different from previous works in that they are intended for hydrogen-bonded systems, while others were built using experimental frequencies of covalently bonded systems. An application of our scaling factors to the vibrational frequencies of water dimer and thermodynamic functions of 11 larger water clusters highlights the importance of anharmonic effects in hydrogen-bonded systems.

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−1 clusters 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.

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 (

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.

Exploration of the Potential Energy Surfaces, Prediction of Atmospheric Concentrations, and Prediction of Vibrational Spectra for the HO2···(H2O)n (n = 1−2) Hydrogen Bonded Complexes

The hydroperoxy radical (HO2) plays a critical role in Earth's atmospheric chemistry as a component of many important reactions. The self-reaction of hydroperoxy radicals in the gas phase is strongly affected by the presence of water vapor. In this work, we explore the potential energy surfaces of hydroperoxy radicals hydrogen bonded to one or two water molecules, and predict atmospheric concentrations and vibrational spectra of these complexes. We predict that when the HO2 concentration is on the order of 108molecules·cm-3 at 298 K, that the number of HO2···H2O complexes is on the order of 107molecules·cm-3 and the number of HO2···(H2O)2 complexes is on the order of 106molecules·cm-3. Using the computed abundance of HO2···H2O, we predict that, at 298 K, the bimolecular rate constant for HO2···H2O + HO2 is about 10 times that for HO2 + HO2.

Experimental and Theoretical Study of the OH Vibrational Spectra and Overtone Chemistry of Gas-Phase Vinylacetic Acid

In this study we present the gas-phase vibrational spectrum of vinylacetic acid with a focus on the ν = 1−5 vibrational states of the OH stretching transitions. Cross sections for ν = 1, 2, 4 and 5 of the OH stretching vibrational transitions are derived on the basis of the vapor pressure data obtained for vinylacetic acid. Ab initio calculations are used to assist in the band assignments of the experimental spectra, and to determine the threshold for the decarboxylation of vinylacetic acid. When compared to the theoretical energy barrier to decarboxylation, it is found that the νOH = 4 transition with thermal excitation of low frequency modes or rotational motion and νOH = 5 transitions have sufficient energy for the reaction to proceed following overtone excitation.

Vibrational spectroscopic monitoring of CO2-O energy transfer: cooling processes in atmospheres of Venus & Mars

The vibrational excitation of CO2 by a fast-moving O atom followed by infrared emission from the vibrationally excited CO2 has been shown to be an important cooling mechanism in the upper atmospheresof Venus, Earth and Mars. We are trying to determine more precisely the efficiency (rate coefficient) of the CO2-O vibrational energy transfer. For experimental ease the reverse reaction is used, i.e. collision of a vibrationally excited CO2 with atomic O, where we are able to convert to the atmospherically relevant reaction via a known equilibrium constant. The goal of this experiment was to measure the magnitudes of rate coefficients for vibrational energy states above the first excited state, a bending mode in CO2. An isotope of CO2, 13CO2, was used for experimental ease. The rate coefficients for given vibrational energy transfers in 13CO2 are not significantly different from 12CO2 at this level of precision. A slow-flowing gas mixture was flowed through a reaction cell: 13CO2 (vibrational specie of interest), O3(atomic O source), and Ar (bath gas). Transient diode laser absorption spectroscopy was used to monitor thechanging absorption of certain vibrational modes of 13CO2 after a UV pulse from a Nd:YAG laser was fired. Ozone absorbed the UV pulse in a process which vibrationally excited 13CO2 and liberated atomic O.Transient absorption signals were obtained by tuning the diode laser frequency to an appropriate ν3 transition and monitoring the population as a function of time following the Nd:YAG pulse. Transient absorption curves were obtained for various O atom concentrations to determine the rate coefficient of interest. Therotational states of the transitions used for detection were difficult to identify, though their short reequilibration timescale made the identification irrelevant for vibrational energy transfer measurements. The rate coefficient for quenching of the (1000) state was found to be (4 ± 8) x 10-12 cm3 s-1 which is the same order of magnitude as the lowest-energy bend-excited mode: (1.8 ± 0.3) x 10-12 cm3 s-1. More data is necessary before it can be certain that the numerical difference between the two is real.

Vibrational Relaxation of N2O (v2) by Atomic Oxygen Measured by Tunable Diode Laser Absorption Spectroscopy

Through the use of Transient Diode Laser Absorption Spectroscopy (TDLAS), the rate coefficient for the vibrational relaxation of N2O (ν2) by O(3P) at room temperature (32 ºC)) was determined to be (1.51 ± 0.11)x10-12 cm3molecule-1sec-1. A Q-switched, frequency quadrupled (266 nm) Nd:YAG laser pulse was used as the pump for this experiment. This pulse caused the photodissociation of O3 into O2 and O atoms.Excited oxygen (O(1D)) was collisionally quenched to ground state (O(3P)) by Ar and/or Xe. Photodissociation also caused a temperature jump within the system, exciting the ν2 state of N2O molecules. Population in the ν2 state was monitored through a TDLASobservation of a ν3 transition. Data were fit using a Visual Fortran 6.0 Global Fitting program. Analysis of room temperature data taken using only Ar to quench O atoms to the ground state gave the same rate coefficient as analysis of data taken using an Ar/Xe mixture, suggesting Ar alone is a sufficient bath gas. Experimentation was alsoperformed at -27 ºC and -82 ºC for a temperature dependence analysis. A linear regression analysis gave a rate coefficient dependence on temperature of ... for the rate coefficient of the vibrational relaxation of N2O (ν2) by atomic oxygen.

Vibrational Relaxation of 13CO2(v2) by Atomic Oxygen

Carbon dioxide (CO2) has been of recent interest due to the issue of greenhouse cooling in the upper atmosphere by species such as CO2 and NO. In the Earth’s upper atmosphere, between altitudes of 75 and 110 km, a collisional energy exchange occurs between CO2 and atomic oxygen, which promotes a population of ground state CO2 to the bend excited state. The relaxation of CO2 following this excitation is characterized by spontaneous emission of 15-μm. Most of this energy is emitted away from Earth. Due to the low density in the upper atmosphere, most of this energy is not reabsorbed and thus escapes into space, leading to a local cooling effect in the upper atmosphere. To determine the efficiency of the CO2- O atom collisional energy exchange, transient diode laser absorption spectroscopy was used to monitor the population of the first vibrationally excited state, 13CO2(0110) or ν2, as a function of time. The rate coefficient, kO(ν2), for the vibrational relaxation 13CO2 (ν2)-O was determined by fitting laboratory measurements using a home-written linear least squares algorithm. The rate coefficient, kO(ν2), of the vibrational relaxation of 13CO2(ν2), by atomic oxygen at room temperature was determined to be (1.6 ± 0.3 x 10-12 cm3 s-1), which is within the uncertainty of the rate coefficient previously found in this group for 12CO2(ν2) relaxation. The cold temperature kO(ν2) values were determined to be: (2.1 ± 0.8) x 10-12 cm3 s-1 at Tfinal = 274 K, (1.8 ± 0.3) x 10-12 cm3 s-1 at Tfinal = 239 K, (2 ± 1) x 10-12 cm3 s-1 at Tfinal = 208 K, and (1.7 ± 0.3) x 10-12 cm3 s-1 at Tfinal = 186 K. These data did not show a definitive negative temperature dependence comparable to that found for 12CO2 previously.

Vibrational Relaxation of O3(nu2) by O((3)P)

Laboratory measurements of the rate coefficient for quenching of O3(nu2) by ground-state atomic oxygen, kO(nu2), at room temperature are presented. kO(nu2) is currently not well known and is necessary for appropriate nonlocal thermodynamic equilibrium modeling of the upper mesosphere and lower thermosphere. In this work, a 266 nm laser pulse photolyzes a small amount of O3 in a slow-flowing gas mixture of O3, Xe, and Ar. This process simultaneously produces atomic oxygen and increases the temperature of the gas mixture slightly, thereby increasing the population in the O3(nu2) state. Transient diode laser absorption spectroscopy is used to monitor the populations of the O3(nu2) and ground vibrational states as the system re-equilibrates. Relaxation rates are measured over a range of quencher concentrations to extract the rate coefficient of interest. The value of kO(nu2) was determined to be (2.2 0.5) * 10(-12) cm(3) s(-1).