Approaching the full set of energy levels of water

Thermal rate constants and kinetic isotope effects of the H + H2O2 reactions: barrier height and reaction energy from single- and multireference methods

CONTEXT

One of the more significant sub-mechanisms of H2/O2 combustion involves the reaction of hydrogen peroxide with hydrogen atoms (H + H2O2), resulting in the production of OH + H2O (R1) and H2 + HO2 (R2) paths. Previous experimental and ab initio calculations reveal some variations in the barrier height for (R1). To improve the energetics of both (R1) and (R2), single reference and multireference ab initio methods are employed, and the rate constants and H/D kinetic isotope effects (KIEs) are calculated as a function of temperature. For (R1), the best results for the barrier height and reaction energies computed with the CASPT2(15,11)/aug-cc-pV6Z are 5.2 and - 70.3 kcal.mol-1, respectively. CCSD(T)/aug-cc-pV5Z + CV (core-valence) calculations for (R2) give 9.7 and - 15.6 kcal.mol-1 to those parameters. The CVT/SCT rate constants of both paths agree well with the fitted rate constants from uncertainty-weighted statistical analysis of the 14-mechanism of H2/O2. The kinetic isotopic effect (kH/kD) for the reaction D + H2O2 → DH + HO2 was found to be 0.47, which is in excellent agreement with the experimental value of 0.43.

METHODS

The structures of reactants, transition state, and products of (R1) and (R2) are calculated with the aug-cc-pVTZ basis set and M062X DFT, CCSD(T), and CASSCF methods. The barrier heights and reaction energies of (R1) and (R2) are computed using the M06-2X, CCSD(T), MRCI, and CASPT2 methods and various basis sets. The rate constants are calculated with the variational transition state theory including multidimensional tunneling corrections (VTST-MT), with potential energy surfaces built by the M06-2X/aug-cc-pVTZ approach.

High Accuracy ab Initio Calculations of Rotational-Vibrational Levels of the HCN/HNC System

Highly accurate ab initio calculations of vibrational and rotational-vibrational energy levels of the HCN/HNC (hydrogen cyanide/hydrogen isocyanide) isomerising system are presented for several isotopologues. All-electron multireference configuration interaction (MRCI) electronic structure calculations were performed using basis sets up to aug-cc-pCV6Z on a grid of 1541 geometries. The ab initio energies were used to produce an analytical potential energy surface (PES) describing the two minima simultaneously. An adiabatic Born-Oppenheimer diagonal correction (BODC) correction surface as well as a relativistic correction surface were also calculated. These surfaces were used to compute vibrational and rotational-vibrational energy levels up to 25 000 cm^-1 which reproduce the extensive set of experimentally known HCN/HNC levels with a root-mean-square deviation σ = 1.5 cm^-1. We studied the effect of nonadiabatic effects by introducing opportune radial and angular corrections to the nuclear kinetic energy operator. Empirical determination of two nonadiabatic parameters results in observed energies up to 7000 cm^-1 for four HCN isotopologues (HCN, DCN, H¹³CN, and HC¹⁵N) being reproduced with σ = 0.37 cm^-1. The height of the isomerization barrier, the isomerization energy and the dissociation energy were computed using a number of models; our best results are 16 809.4, 5312.8, and 43 729 cm^-1, respectively.

A database of water transitions from experiment and theory (IUPAC Technical Report)

The report and results of an IUPAC Task Group (TG) formed in 2004 on “A Database of Water Transitions from Experiment and Theory” (Project No. 2004-035-1-100) are presented. Energy levels and recommended labels involving exact and approximate quantum numbers for the main isotopologues of water in the gas phase, H216O, H218O, H217O, HD16O, HD18O, HD17O, D216O, D218O, and D217O, are determined from measured transition frequencies. The transition frequencies and energy levels are validated using first-principles nuclear motion computations and the MARVEL (measured active rotational–vibrational energy levels) approach. The extensive data including lines and levels are required for analysis and synthesis of spectra, thermochemical applications, the construction of theoretical models, and the removal of spectral contamination by ubiquitous water lines. These datasets can also be used to assess where measurements are lacking for each isotopologue and to provide accurate frequencies for many yet-to-be measured transitions. The lack of high-quality frequency calibration standards in the near infrared is identified as an issue that has hindered the determination of high-accuracy energy levels at higher frequencies. The generation of spectra using the MARVEL energy levels combined with transition intensities computed using high accuracy ab initio dipole moment surfaces are discussed. A recommendation of the TG is for further work to identify a single, suitable model to represent pressure- (and temperature-) dependent line profiles more accurately than Voigt profiles.

Absorption cross sections of surface-adsorbed H2O in the 295-370 nm region and heterogeneous nucleation of H2O on fused silica surfaces

We have determined absorption cross sections of a monolayer of H2O adsorbed on the fused silica surfaces in the 295-370 nm region at 293 ± 1 K by using Brewster angle cavity ring-down spectroscopy. Absorption cross sections of surface-adsorbed H2O vary between (4.66 ± 0.83) × 10(-20) and (1.73 ± 0.52) × 10(-21) cm(2)/molecule over this wavelength range, where errors quoted represent experimental scatter (1σ). Our experimental study provides direct evidence that surface-adsorbed H2O is an absorber of the near UV solar radiation. We also varied the H2O pressure in the surface study cell over the 0.01-17 Torr range and obtained probe laser absorptions at 295, 340, and 350 nm by multilayer of adsorbed H2O molecules until the heterogeneous nucleation of water occurred on fused silica surfaces. The average absorption cross sections of multilayer adsorbed H2O are (2.17 ± 0.53) × 10(-20), (2.48 ± 0.67) × 10(-21), and (2.34 ± 0.59) × 10(-21) cm(2)/molecule at 295, 340, and 350 nm. The average absorption cross sections of transitional H2O layer are (6.06 ± 2.73) × 10(-20), (6.48 ± 3.85) × 10(-21), and (8.04 ± 4.92) × 10(-21) cm(2)/molecule at 295, 340, and 350 nm. The average thin water film absorption cross sections are (2.39 ± 0.50) × 10(-19), (3.21 ± 0.81) × 10(-20), and (3.37 ± 0.94) × 10(-20) cm(2)/molecule at 295 nm, 340 nm, and 350 nm. Atmospheric implications of the results are discussed.

Global spectroscopy of the water monomer

Given the large energy required for its electronic excitation, the most important properties of the water molecule are governed by its ground potential energy surface (PES). Novel experiments are now able to probe this surface over a very extended energy range, requiring new theoretical procedures for their interpretation. As part of this study, a new, accurate, global spectroscopic-quality PES and a new, accurate, global dipole moment surface are developed. They are used for the computation of the high-resolution spectrum of water up to the first dissociation limit and beyond as well as for the determination of Stark coefficients for high-lying states. The water PES has been determined by combined ab initio and semi-empirical studies. As a first step, a very accurate, global, ab initio PES was determined using the all-electron, internally contracted multi-reference configuration interaction technique together with a large Gaussian basis set. Scalar relativistic energy corrections are also determined in order to move the energy determinations close to the relativistic complete basis set full configuration interaction limit. The electronic energies were computed for a set of about 2500 geometries, covering carefully selected configurations from equilibrium up to dissociation. Nuclear motion computations using this PES reproduce the observed energy levels up to 39 000 cm(-1) with an accuracy of better than 10 cm(-1). Line positions and widths of resonant states above dissociation show an agreement with experiment of about 50 cm(-1). An improved semi-empirical PES is produced by fitting the ab initio PES to accurate experimental data, resulting in greatly improved accuracy, with a maximum deviation of about 1 cm(-1) for all vibrational band origins. Theoretical results based on this semi-empirical surface are compared with experimental data for energies starting at 27 000 cm(-1), going all the way up to dissociation at about 41 000 cm(-1) and a few hundred wavenumbers beyond it.

State-resolved spectroscopy of high vibrational levels of water up to the dissociative continuum

We summarize here our experimental studies of the high rovibrational energy levels of water. The use of double-resonance vibrational overtone excitation followed by energy-selective photofragmentation and laser-induced fluorescence detection of OH fragments allowed us to measure previously inaccessible rovibrational energies above the seventh OH-stretch overtone. Extension of the experimental approach to triple-resonance excitation provides access to rovibrational levels via transitions with significant transition dipole moments (mainly OH-stretch overtones) up to the dissociation threshold of the O-H bond. A collisionally assisted excitation scheme enables us to probe vibrations that are not readily accessible via pure laser excitation. Observation of the continuous absorption onset yields a precise value for the O-H bond dissociation threshold, 41 145.94 ± 0.15 cm(-1). Finally, we detect long-lived resonances as sharp peaks in spectra above the dissociation threshold.

Communication: Feshbach resonances in the water molecule revealed by state-selective spectroscopy

We employ triple-resonance vibrational overtone excitation to access quasibound states of water from several fully characterized bound states of the molecule. Comparison of the measured dissociation spectra allows a rigorous assignment of rotational quantum numbers J, nuclear spin and parity, and a tentative vibrational characterization of the observed resonances. Their asymmetrical shapes (Fano profiles) reflect interference of dipole moments for transitions to these resonances with that to the dissociative continuum. The assignments and Fano profile parameters of the resonances stand as a benchmark for the extension of accurate quantum-mechanical calculations to activated complexes of water. The narrow widths of some of these resonances indicate that water molecules may survive for as long as up to 60 ps in states above the dissociation threshold. We consider the possible implication of such long-lived states for the kinetics of water dissociation and the OH+H association reaction.

State-selective spectroscopy of water up to its first dissociation limit

A joint experimental and first-principles quantum chemical study of the vibration-rotation states of the water molecule up to its first dissociation limit is presented. Triple-resonance, quantum state-selective spectroscopy is used to probe the entire ladder of water's stretching vibrations up to 19 quanta of OH stretch, the last stretching state below dissociation. A new ground state potential energy surface of water is calculated using a large basis set and an all-electron, multireference configuration interaction procedure, which is augmented by relativistic corrections and fitted to a flexible functional form appropriate for a dissociating system. Variational nuclear motion calculations on this surface are used to give vibrational assignments. A total of 44 new vibrational states and 366 rotation-vibration energy levels are characterized; these span the region from 35,508 to 41,126 cm(-1) above the vibrational ground state.