Anharmonicity and hydrogen bonding

Vibrational Spectroscopy of the Water Dimer at Jet-Cooled and Atmospheric Temperatures

The vibrational spectroscopy of the water dimer provides an understanding of basic hydrogen bonding in water clusters, and with about one water dimer for every 1,000 water molecules, it plays a critical role in atmospheric science. Here, we review how the experimental and theoretical progress of the past decades has improved our understanding of water dimer vibrational spectroscopy under both cold and warm conditions. We focus on the intramolecular OH-stretching transitions of the donor unit, because these are the ones mostly affected by dimer formation and because their assignment has proven a challenge. We review cold experimental results from early matrix isolation to recent mass-selected jet expansion techniques and, in parallel, the improvements in the theoretical anharmonic models. We discuss and illustrate changes in the vibrational spectra of complexes upon increasing temperature, and the difficulties in recording and calculating these spectra. In the atmosphere, water dimer spectra at ambient temperature are crucial. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

On the synergy of matrix-isolation infrared spectroscopy and vibrational configuration interaction computations

The key feature of matrix-isolation infrared (MI-IR) spectroscopy is the isolation of single guest molecules in a host system at cryogenic conditions. The matrix mostly hinders rotation of the guest molecule, providing access to pure vibrational features. Vibrational self-consistent field (VSCF) and configuration interaction computations (VCI) on ab initio multimode potential energy surfaces (PES) give rise to anharmonic vibrational spectra. In a single-sourced combination of these experimental and computational approaches, we have established an iterative spectroscopic characterization procedure. The present article reviews the scope of this procedure by highlighting the strengths and limitations based on the examples of water, carbon dioxide, methane, methanol, and fluoroethane. An assessment of setups for the construction of the multimode PES on the example of methanol demonstrates that CCSD(T)-F12 level of theory is preferable to compute (a) accurate vibrational frequencies and (b) equilibrium or vibrationally averaged structural parameters. Our procedure has allowed us to uniquely assign unknown or disputed bands and enabled us to clarify problematic spectral regions that are crowded with combination bands and overtones. Besides spectroscopic assignment, the excellent agreement between theory and experiment paves the way to tackle questions of rather fundamental nature as to whether or not matrix effects are systematic, and it shows the limits of conventional notations used by spectroscopists.

Quantum aspects of torsional vibrations in the HO3H, DO3H and DO3D molecules

The hydrogen trioxide (HT) molecule HOOOH is 1) a prototype for a class of molecules of the form X(ZY)₂ with two equivalent internal tops (ZY) and 2) the second representative of the polyoxides series of the form HO_nH (n ≥ 2). Due to this, it is the subject of close attention of researchers. In this paper, we performed a group theory analysis of the torsional and spin states of the HOOOH, DOOOH, and DOOOD molecules. The relationships have been established between the symmetry species of the C_2V(M) molecular symmetry group to which the HOOOH and DOOOD molecules belong, and the symmetry species of the C₂ and C_S point groups to which the equilibrium configurations of trans- and cis-conformers of the above molecules belong, respectively. 2D PES and 2D surfaces of kinematic coefficients related to torsional vibrations of hydroxyl groups were calculated at the complete basis set (CBS) limit by extrapolating the results of calculations at the MP2/cc-pVTZ and MP2/cc-pVQZ levels of theory. For all the three molecules, the energies of the stationary torsional states were computed using the Fourier method for a numerical solution of the 2D vibrational Schrödinger equation. Symmetry species of torsional states and the values of quantum numbers defining the types of torsional vibrations were found by analyzing the torsional wave functions. The selection rules for transitions between torsional states in the dipole approximation were also formulated.

Near-Infrared Spectroscopy and Anharmonic Theory of Protonated Water Clusters: Higher Elevations in the Hydrogen Bonding Landscape

Near-infrared spectroscopy measurements are presented for protonated water clusters, H⁺(H₂O) _n, in the size range of n = 1-8. Clusters are produced in a pulsed-discharge supersonic expansion, mass selected, and studied with infrared laser photodissociation spectroscopy in the regions of 3600-4550 and 4850-7350 cm^-1. Although there is some variation with cluster size, the main features of these spectra are a broad absorption near 5300 cm^-1, a sharp doublet near 7200 cm^-1, as well as a structured absorption near 4100 cm^-1 for n ≥ 2. The vibrational patterns measured for the hydronium, Zundel, and Eigen ions are compared to those predicted by different forms of anharmonic theory. Second-order vibrational perturbation theory (VPT2) and a local mode treatment of the OH stretches both capture key aspects of the spectra but suffer understandable deficiencies in the quantitative description of band positions and intensities.

Temperature dependence of the intensity of the vibration-rotational absorption band ν2 of H2O trapped in an argon matrix

Using two sets of effective rotational constants for the ground (000) and the excited bending (010) vibrational states the calculation of frequencies and intensities of vibration-rotational transitions for J^″=0-2; and J^'=0-3; was carried out in frame of the model of a rigid asymmetric top for temperatures from 0 to 40K. The calculation of the intensities of vibration-rotational absorption bands of H₂O in an Ar matrix was carried out both for thermodynamic equilibrium and for the case of non-equilibrium population of para- and ortho-states. For the analysis of possible interaction of vibration-rotational and translational motions of a water molecule in an Ar matrix by 3D Schrödinger equation solving using discrete variable representation (DVR) method, calculations of translational frequencies of Н₂О in a cage formed after one argon atom deleting were carried out. The results of theoretical calculations were compared to experimental data taken from literature.

Acetic acid-water complex: The first observation of structures containing the higher-energy acetic acid conformer

Non-covalent interaction of acetic acid (AA) and water is studied experimentally by IR spectroscopy in a nitrogen matrix and theoretically at the MP2 and coupled-cluster with single and double and perturbative triple excitations [CCSD(T)]/6-311++G(2d,2p) levels of theory. This work is focused on the first preparation and characterization of complexes of higher-energy (cis) conformer of AA with water. The calculations show three 1:1 structures for the trans-AA⋯H2O complexes and three 1:1 structures for the cis-AA⋯H2O complexes. Two trans-AA⋯H2O and two cis-AA⋯H2O complexes are found and structurally assigned in the experiments. The two cis-AA⋯ ⋅ H2O complexes are obtained by annealing of a matrix containing water and cis-AA molecules prepared by selective vibrational excitation of the ground-state trans form. The less stable trans-AA⋯H2O complex is obtained by vibrational excitation of the less stable cis-AA⋯H2O complex. In addition, the 1:2 complexes of trans-AA and cis-AA with water molecules are studied computationally and the most stable forms of the 1:2 complexes are experimentally identified.

Phonon coupling of water monomers in a solid nitrogen matrix

The infrared absorption spectra of the H2O, HDO, and D2O monomers isolated in solid N2 have been recorded at various temperatures between 4 and 30 K. A study of the absorption features of the ν1, ν2, and ν3 vibrational modes for each monomer shows their optical line shapes to be strongly temperature dependent. For all three modes, a decrease in the absorption amplitude and a proportional broadening of the linewidth was observed with increasing temperature, while the integrated absorbance remained constant. These observations were explained in terms of phonon coupling, by which high frequency intramolecular modes decay by exciting matrix phonons. Fits of the linewidth for the lowest frequency ν2 vibrational mode to the predicted vibrational relaxation rate in a solid medium gave average phonon mode frequencies consistent with the Debye frequency for solid N2.

Communication: determination of the bond dissociation energy (D0) of the water dimer, (H2O)2, by velocity map imaging

The bond dissociation energy (D(0)) of the water dimer is determined by using state-to-state vibrational predissociation measurements following excitation of the bound OH stretch fundamental of the donor unit of the dimer. Velocity map imaging and resonance-enhanced multiphoton ionization (REMPI) are used to determine pair-correlated product velocity and translational energy distributions. H(2)O fragments are detected in the ground vibrational (000) and the first excited bending (010) states by 2 + 1 REMPI via the C̃ (1)B(1) (000) ← X̃ (1)A(1) (000 and 010) transitions. The fragments' velocity and center-of-mass translational energy distributions are determined from images of selected rovibrational levels of H(2)O. An accurate value for D(0) is obtained by fitting both the structure in the images and the maximum velocity of the fragments. This value, D(0) = 1105 ± 10 cm(-1) (13.2 ± 0.12 kJ/mol), is in excellent agreement with the recent theoretical value of D(0) = 1103 ± 4 cm(-1) (13.2 ± 0.05 kJ∕mol) suggested as a benchmark by Shank et al. [J. Chem. Phys. 130, 144314 (2009)].

Interaction energy and the shift in OH stretch frequency on hydrogen bonding for the H2O --> H2O, CH3OH --> H2O, and H2O --> CH3OH dimers

The ability to use calculated OH frequencies to assign experimentally observed peaks in hydrogen bonded systems hinges on the accuracy of the calculation. Here we test the ability of several commonly employed model chemistries--HF, MP2, and several density functionals paired with the 6-31+G(d) and 6-311++G(d,p) basis sets--to calculate the interaction energy (D(e)) and shift in OH stretch fundamental frequency on dimerization (delta(nu)) for the H(2)O --> H(2)O, CH(3)OH --> H(2)O, and H(2)O --> CH(3)OH dimers (where for X --> Y, X is the hydrogen bond donor and Y the acceptor). We quantify the error in D(e) and delta(nu) by comparison to experiment and high level calculation and, using a simple model, evaluate how error in D(e) propagates to delta(nu). We find that B3LYP and MPWB1K perform best of the density functional methods studied, that their accuracy in calculating delta(nu) is approximately 30-50 cm(-1) and that correcting for error in D(e) does little to heighten agreement between the calculated and experimental delta(nu). Accuracy of calculated delta(nu) is also shown to vary as a function of hydrogen bond donor: while the PBE and TPSS functionals perform best in the calculation of delta(nu) for the CH(3)OH --> H(2)O dimer their performance is relatively poor in describing H(2)O --> H(2)O and H(2)O --> CH(3)OH.

Hydrogen bonded OH-stretching vibration in the water dimer

We have calculated the frequencies and intensities of the hydrogen-bonded OH-stretching transitions in the water dimer complex. The potential-energy curve and dipole-moment function are calculated ab initio at the coupled cluster with singles, doubles, and perturbative triples level of theory with correlation-consistent Dunning basis sets. The vibrational frequencies and wavefunctions are found from a numerical solution to a one-dimensional Schrödinger equation. The corresponding transition intensities are found from numerical integration of these vibrational wavefunctions with the ab initio calculated dipole moment function. We investigate the effect of counterpoise correcting both the potential-energy surface and dipole-moment function. We find that the effect of using a numeric potential is significant for higher overtones and that inclusion of a counterpoise correction for basis set superposition error is important.

Geometrical behavior of hydrogen bonding patterns in the alpha-dodecyl-omega-hydroxy-tris(oxyethylene)-water system monitored by near infrared spectroscopy

Changes in the geometry of hydrogen bonding patterns in the alpha-dodecyl-omega-hydroxy-tris(oxyethylene) (C(12)E(3))-water system have been investigated by near infrared (NIR) spectroscopy. In the 5,300-4,600 cm(-1) region, the characteristic bands for C(12)E(3) and water can be separately investigated, since the combination bands of the OH stretching and its COH bending of alcohols are observed at 5,000-4,650 cm(-1), whereas the combination bands of the OH stretching and its HOH bending of water, at 5,300-5,000 cm(-1). The NIR result has revealed that the addition of water to C(12)E(3) promotes the formation of the OHcdots, three dots, centeredOHcdots, three dots, centeredO hydrogen bonds.

Overtone spectroscopy of H2O clusters in the vOH=2 manifold: Infrared-ultraviolet vibrationally mediated dissociation studies

Spectroscopy and predissociation dynamics of (H2O)2 and Ar-H2O are investigated with vibrationally mediated dissociation (VMD) techniques, wherein upsilon(OH) = 2 overtones of the complexes are selectively prepared with direct infrared pumping, followed by 193 nm photolysis of the excited H2O molecules. As a function of relative laser timing, the photolysis breaks H2O into OH and H fragments either (i) directly inside the complex or (ii) after the complex undergoes vibrational predissociation, with the nascent quantum state distribution of the OH photofragment probed via laser-induced fluorescence. This capability provides the first rotationally resolved spectroscopic analysis of (H2O)2 in the first overtone region and vibrational predissociation dynamics of water dimer and Ar-water clusters. The sensitivity of the VMD approach permits several upsilon(OH) = 2 overtone bands to be observed, the spectroscopic assignment of which is discussed in the context of recent anharmonic theoretical calculations.

Effect of Temperature and Concentration on the Structure of N-Methylacetamide−Water Complexes: Near-Infrared Spectroscopic Study

Generalized two-dimensional (2D) FT-NIR correlation spectroscopy and chemometric methods have been used to study temperature-dependent spectral changes in pure N-methylacetamide (NMA) and NMA-water mixtures. We also examined the effect of varying water content on the structure of the mixture. It has been found that the extent of self-association of NMA in CCl4 is very high; the association occurs even at concentration of 0.001 M. In the pure liquid NMA, the population of the monomers is negligible and the structure is dominated by the linear associates. An increase in temperature reduces the number of hydrogen bonds, but in contrast to alcohols their strength remains nearly the same. This reflects a difference in the mechanism of thermal breaking of the associates of NMA and alcohols. The present results reveal that the interaction between NMA and water in the NMA-rich region (X(H2O) < 0.1) does not have a significant effect on the intrinsic structure of NMA. The structure of NMA is dominant, and the molecules of water do not form separate clusters but are dispersed and incorporated into the structure of NMA. We did not observe the presence of the free OH groups in the mixture. This led to the suggestion that each molecule of water forms two hydrogen bonds to two different molecules of NMA. An analysis of the asynchronous spectra reveals that most of the peaks observed in the asynchronous spectra, constructed from the temperature-dependent data, simply result from the frequency shift. This assumption is supported by the simulation studies.