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.

Structure and thermodynamics of H3O+(H2O)8 clusters: A combined molecular dynamics and quantum mechanics approach

We have studied the structure and stability of (H3O+)(H2O)8 clusters using a combination of molecular dynamics sampling and high-level ab initio calculations. 20 distinct oxygen frameworks are found within 2 kcal/mol of the electronic or standard Gibbs free energy minimum. The impact of quantum zero-point vibrational corrections on the relative stability of these isomers is quite significant. The box-like isomers are favored in terms of electronic energy, but with the inclusion of zero-point vibrational corrections and entropic effects tree-like isomers are favored at higher temperatures. Under conditions from 0 to 298.15 K, the global minimum is predicted to be a tree-like structure with one dangling singly coordinated water molecule. Above 298.15 K, higher entropy tree-like isomers with two or more singly coordinated water molecules are favored. These assignments are generally consistent with experimental IR spectra of (H3O+)(H2O)8 obtained at 150 K.

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.

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.

Investigation of the Potential Energy Surface for the First Step in the Alkaline Hydrolysis of Methyl Acetate

The potential energy surface for the first step of the alkaline hydrolysis of methyl acetate was explored by a variety of methods. The conformational search routine within SPARTAN was used to determine the lowest energy am1 and pm3 structures for the anionic tetrahedral intermediate. Ab initio single point and geometry optimization calculations were performed to determine the lowest energy conformer, and the linear synchronous transition (lst) method was used to provide an initial structure for transition state optimization. Transition states were obtained at the am1, pm3, 3-21G, and 3-21 + G levels of theory. These transition states were compared with the anionic tetrahedral intermediates to examine the assumption that the intermediate is a good model for the transition state. In addition, the Cramer/Truhlar sm3 solvation model was used at the semiempirical level to compare gas phase and aqueous alkaline hydrolysis of methyl acetate.

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.

The Effect of Grade Point Average on Undergraduates' Understanding of Heat Transfer

Engineering students continue to develop and show misconceptions due to prior knowledge and experiences (Miller, Streveler, Olds, Chi, Nelson, & Geist, 2007). Misconceptions have been documented in students’ understanding of heat transfer(Krause, Decker, Niska, Alford, & Griffin, 2003) by concept inventories (e.g., Jacobi,Martin, Mitchell, & Newell, 2003; Nottis, Prince, Vigeant, Nelson, & Hartsock, 2009). Students’ conceptual understanding has also been shown to vary by grade point average (Nottis et al., 2009). Inquiry-based activities (Nottis, Prince, & Vigeant, 2010) haveshown some success over traditional instructional methods (Tasoglu & Bakac, 2010) in altering misconceptions. The purpose of the current study was to determine whether undergraduate engineering students’ understanding of heat transfer concepts significantly changed after instruction with eight inquiry-based activities (Prince & Felder, 2007) supplementing instruction and whether students’ self reported GPA and prior knowledge, as measured by completion of specific engineering courses, affected these changes. The Heat and Energy Concept Inventory (Prince, Vigeant, & Nottis, 2010) was used to assess conceptual understanding. It was found that conceptual understanding significantly increased from pre- to post-test. It was also found that GPA had an effect on conceptual understanding of heat transfer; significant differences were found in post-test scores onthe concept inventory between GPA groups. However, there were mixed results when courses previously taken were analyzed. Future research should strive to analyze how prior knowledge effects conceptual understanding and aim to reduce the limitations of the current study such as, sampling method and methods of measuring GPA and priorknowledge.

Assessing the Validity of Brain Glucose Concentration Measurement Using Microdialysis

Brain functions, such as learning, orchestrating locomotion, memory recall, and processing information, all require glucose as a source of energy. During these functions, the glucose concentration decreases as the glucose is being consumed by brain cells. By measuring this drop in concentration, it is possible to determine which parts of the brain are used during specific functions and consequently, how much energy the brain requires to complete the function. One way to measure in vivo brain glucose levels is with a microdialysis probe. The drawback of this analytical procedure, as with many steadystate fluid flow systems, is that the probe fluid will not reach equilibrium with the brain fluid. Therefore, brain concentration is inferred by taking samples at multiple inlet glucose concentrations and finding a point of convergence. The goal of this thesis is to create a three-dimensional, time-dependent, finite element representation of the brainprobe system in COMSOL 4.2 that describes the diffusion and convection of glucose. Once validated with experimental results, this model can then be used to test parameters that experiments cannot access. When simulations were run using published values for physical constants (i.e. diffusivities, density and viscosity), the resulting glucose model concentrations were within the error of the experimental data. This verifies that the model is an accurate representation of the physical system. In addition to accurately describing the experimental brain-probe system, the model I created is able to show the validity of zero-net-flux for a given experiment. A useful discovery is that the slope of the zero-net-flux line is dependent on perfusate flow rate and diffusion coefficients, but it is independent of brain glucose concentrations. The model was simplified with the realization that the perfusate is at thermal equilibrium with the brain throughout the active region of the probe. This allowed for the assumption that all model parameters are temperature independent. The time to steady-state for the probe is approximately one minute. However, the signal degrades in the exit tubing due to Taylor dispersion, on the order of two minutes for two meters of tubing. Given an analytical instrument requiring a five μL aliquot, the smallest brain process measurable for this system is 13 minutes.

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.