Computing excited OH stretch states of water dimer in 12D using contracted intermolecular and intramolecular basis functions

Infrared spectra of partially deuterated water dimers in the fundamental O-D stretch region

Spectra of mixed H/D water dimers are studied in the mid-infrared region of the O-D stretch fundamental (2630-2800 cm-1) using a pulsed supersonic slit jet and a tunable optical parametric oscillator infrared source. Over 30 bands, belonging to nine of the ten possible isotopologues (only H2O-HOD is missed), are observed and analyzed. The implications for excited state tunneling splittings, lifetime effects, and vibrational shifts are discussed. These are the first significant new experimental results on (gas phase) mixed water dimers in over 25 years, and they are valuable for testing water dimer theoretical calculations, a field which continues to be of lively current interest.

Tunneling splittings in the vibrationally excited states of water trimer

Tunneling splitting (TS) patterns in vibrationally excited states of the water trimer are calculated, taking into account six tunneling pathways that describe the flips of free OH bonds and five bifurcation mechanisms that break and reform hydrogen bonds in the trimer ring. A tunneling matrix (TM) model is used to derive the energy shifts due to tunneling in terms of the six distinct TM elements in symbolic form. TM elements are calculated using the recently-developed modified WKB (Wentzel-Kramers-Brillouin) method in full dimensionality. Convergence was achieved for the lowest six excited vibrational modes. Bifurcation widths of the pseudorotational quartets are shown to be of comparable size to the ground-state widths, obtained using instanton theory, or increased for some particular modes of vibration. The largest increase is obtained for the excited out-of-phase flip of two adjacent water monomers with free OH bonds pointing in opposite directions relative to the ring plane. Bifurcation widths in (D2O)3 are found to be two orders of magnitude smaller than in (H2O)3. Geometrical arguments were used to explain the order of states in some TS multiplets in vibrationally excited water trimers.

HCl trimer: HCl-stretch excited intramolecular and intermolecular vibrational states from 12D fully coupled quantum calculations employing contracted intra- and inter-molecular bases

We present fully coupled, full-dimensional quantum calculations of the inter- and intra-molecular vibrational states of HCl trimer, a paradigmatic hydrogen-bonded molecular trimer. They are performed utilizing the recently developed methodology for the rigorous 12D quantum treatment of the vibrations of the noncovalently bound trimers of flexible diatomic molecules [Felker and Bačić, J. Chem. Phys. 158, 234109 (2023)], which was previously applied to the HF trimer by us. In this work, the many-body 12D potential energy surface (PES) of (HCl)3 [Mancini and Bowman, J. Phys. Chem. A 118, 7367 (2014)] is employed. The calculations extend to the intramolecular HCl-stretch excited vibrational states of the trimer with one- and two-quanta, together with the low-energy intermolecular vibrational states in the two excited v = 1 intramolecular vibrational manifolds. They reveal significant coupling between the intra- and inter-molecular vibrational modes. The 12D calculations also show that the frequencies of the v = 1 HCl stretching states of the HCl trimer are significantly redshifted relative to those of the isolated HCl monomer. Detailed comparison is made between the results of the 12D calculations on the two-body PES, obtained by removing the three-body term from the original 2 + 3-body PES, and those computed on the 2 + 3-body PES. It demonstrates that the three-body interactions have a strong effect on the trimer binding energy as well as on its intra- and inter-molecular vibrational energy levels. Comparison with the available spectroscopic data shows that good agreement with the experiment is achieved only if the three-body interactions are included. Some low-energy vibrational states localized in a secondary minimum of the PES are characterized as well.

The strength of the OH-bend/OH-stretch Fermi resonance in small water clusters

A novel Raman jet-spectrometer is used to study the Fermi resonance between the OH bending overtone and OH stretching fundamental in small cyclic water clusters (H2O)n with n = 3, 4, 5. The new setup features a recirculating vacuum system which reduces the gas consumption by 2 to 3 orders of magnitude and enables long-term measurements of very weak Raman signals. Raman spectra measured from highly diluted expansions with unprecedented signal-to-noise ratio are presented and cluster-specific intensity ratios and effective coupling constants are derived using Markov-Chain Monte-Carlo methods, yielding a high probability for an almost "perfect" resonance for the tetramer and pentamer, i.e. a close frequency match of bend overtone and stretch fundamental with intensity ratios close to 1, but a larger coupling constant for the trimer, with best estimates close to W5 ≲ 50 cm-1 < W4 ≲ 60 cm-1 < W3 ≈ 65 cm-1.

Spectra of the D2O dimer in the O-D fundamental stretch region: The acceptor symmetric stretch fundamental and new combination bands

The O-D stretch fundamental region of the deuterated water dimer, (D2O)2, is further studied using a pulsed supersonic slit jet and a tunable optical parametric oscillator infrared source. The previously unobserved acceptor symmetric O-D stretch fundamental vibration is detected, with Ka = 0 ← 0 and 1 ← 0 sub-bands at about 2669 and 2674 cm-1, respectively. The analysis indicates that the various water dimer tunneling splittings generally decrease in the excited vibrational state, similar to the three other previously observed O-D stretch fundamentals. Two new (D2O)2 combination bands are observed, giving information on intermolecular vibrations in the excited O-D stretch states. The likely vibrational assignments for these and a previously observed combination band are discussed.

Effect of Temperature on the OH-Stretching Bands of the Methanol Dimer

We present a conceptually simple model for understanding the significant spectral changes that occur with the temperature in the infrared spectra of hydrogen-bound complexes. We have measured room-temperature spectra of the methanol dimer and two deuterated isotopologues in the OH(D)-stretching region. We correctly predict spectral changes observed in the gas phase for the bound OH stretch in the methanol dimer from jet-cooled to room temperature and corroborate this with experimental and theoretical results for deuterated isotopologues. The origin of the observed spectral features is explained based on a reduced-dimensional vibrational model, which includes the two high-frequency OH stretches, the two methyl torsions, and the six intermolecular low-frequency vibrations. Key to the success of the model is a new coordinate definition to describe the intrinsic large-amplitude curvilinear motion of low-frequency vibrations. Despite the deceivingly simple appearance of the room temperature bound OH-stretching fundamental band, it consists of ∼107 vibrational transitions.

Rovibrational states calculations of the H2O-HCN heterodimer with the multiconfiguration time dependent Hartree method

Water and hydrogen cyanide are two of the most common species in space and the atmosphere with the ability of binding to form dimers such as H2O-HCN. In the literature, while calculations characterizing various properties of the H2O-HCN cluster (equilibrium distance, vibrational frequencies and rotational constants) have been done in the past, extensive calculations of the rovibrational states of this system using a reliable quantum dynamical approach have yet to be reported. In this work, we intend to mend that by performing the first calculation of the rovibrational states of the H2O-HCN van der Waals complex on a recently developed potential energy surface. We use the block improved relaxation procedure implemented in the Heidelberg MultiConfiguration Time-Dependent Hartree (MCTDH) package to compute the states of the H2O-HCN isomer, from which we extract the transition frequencies and rotational constants of the complex. We further adapt an approach first suggested by Wang and Carrington-and supported here by analysis routines of the Heidelberg MCTDH package-to properly characterize the computed rovibrational states. The subsequent assignment of rovibrational states was done by theoretical analysis and visual inspection of the wavefunctions. Our simulations provide a Zero Point Energy (ZPE) and intermolecular vibrational frequencies in good agreement with past ab initio calculations. The transition frequencies and rotational constants obtained from our simulations match well with the available experimental data. This work has the broad aim to propose the MCTDH approach as a reliable option to compute and characterize rovibrational states of van der Waals complexes such as the current one.

Microhydration of Tertiary Amines: Robust Resonances in Red-Shifted Water

Tertiary amines are strong hydrogen bond acceptors. When a water molecule donates one of the OH groups, its in-phase stretching vibration wavenumber is decreased to such an extent that it comes close to the water bending overtone. This gives rise to anharmonic phenomena such as classical Fermi resonance, resonance with multiple-quantum dark states, or combination transitions with low-frequency intermolecular modes. These effects, which contribute to the spectral breadth of room-temperature hydrogen-bonded amine complexes, are disentangled by Fourier transform infrared spectroscopy in pulsed supersonic slit jet expansions. Monohydrates of the amines quinuclidine, N-methylpyrrolidine, N-methylpiperidine, and dimethylcyclohexylamine exhibit systematic mode coupling signatures. These suggest relatively fast energy flow out of the excited OH stretching fundamental into intra- and intermolecular degrees of freedom of the hydrogen-bonded water molecule. Trimeric complexes are spectroscopically separated from the amine monohydrates.

Quantum Chemical Investigation of the Cold Water Dimer Spectrum in the First OH-Stretching Overtone Region Provides a New Interpretation

Intramolecular vibrational transition wavenumbers and intensities were calculated in the fundamental HOH-bending, fundamental OH-stretching, first OH-stretching-HOH-bending combination, and first OH-stretching overtone (ΔvOH = 2) regions of the water dimer's spectrum. Furthermore, the rotational-vibrational spectrum was calculated in the ΔvOH = 2 region at 10 K, corresponding to the temperature of the existing jet-expansion experiments. The calculated spectrum was obtained by combining results from a full-dimensional (12D) vibrational and a reduced-dimensional vibrational-rotational-tunneling model. The ΔvOH = 2 spectral region is rich in features due to contributions from multiple vibrational-rotational-tunneling sub-bands. Origins of the experimental vibrational bands depend on the assignment of the observed sub-bands. Based on our calculations, we assign the observed sub-bands, and our reassignment leads to new values for the vibrational band origins of the free donor and antisymmetric acceptor OH-stretching first overtones of ∼7227 and ∼7238 cm-1, respectively. The observed bands with origins at 7192.34 and ∼7366 cm-1 are assigned to the symmetric acceptor OH-stretching first overtone and the OH-stretching combination of the donor, respectively.

H2O-HCN complex: A new potential energy surface and intermolecular rovibrational states from rigorous quantum calculations

In this work the H2O-HCN complex is quantitatively characterized in two ways. First, we report a new rigid-monomer 5D intermolecular potential energy surface (PES) for this complex, calculated using the symmetry-adapted perturbation theory based on density functional theory method. The PES is based on 2833 ab initio points computed employing the aug-cc-pVQZ basis set, utilizing the autoPES code, which provides a site-site analytical fit with the long-range region given by perturbation theory. Next, we present the results of the quantum 5D calculations of the fully coupled intermolecular rovibrational states of the H2O-HCN complex for the total angular momentum J values of 0, 1, and 2, performed on the new PES. These calculations rely on the quantum bound-state methodology developed by us recently and applied to a variety of noncovalently bound binary molecular complexes. The vibrationally averaged ground-state geometry of H2O-HCN determined from the quantum 5D calculations agrees very well with that from the microwave spectroscopic measurements. In addition, the computed ground-state rotational transition frequencies, as well as the B and C rotational constants calculated for the ground state of the complex, are in excellent agreement with the experimental values. The assignment of the calculated intermolecular vibrational states of the H2O-HCN complex is surprisingly challenging. It turns out that only the excitations of the intermolecular stretch mode can be assigned with confidence. The coupling among the angular degrees of freedom (DOFs) of the complex is unusually strong, and as a result most of the excited intermolecular states are unassigned. On the other hand, the coupling of the radial, intermolecular stretch mode and the angular DOFs is weak, allowing straightforward assignment of the excitation of the former.

The first HyDRA challenge for computational vibrational spectroscopy

Vibrational spectroscopy in supersonic jet expansions is a powerful tool to assess molecular aggregates in close to ideal conditions for the benchmarking of quantum chemical approaches. The low temperatures achieved as well as the absence of environment effects allow for a direct comparison between computed and experimental spectra. This provides potential benchmarking data which can be revisited to hone different computational techniques, and it allows for the critical analysis of procedures under the setting of a blind challenge. In the latter case, the final result is unknown to modellers, providing an unbiased testing opportunity for quantum chemical models. In this work, we present the spectroscopic and computational results for the first HyDRA blind challenge. The latter deals with the prediction of water donor stretching vibrations in monohydrates of organic molecules. This edition features a test set of 10 systems. Experimental water donor OH vibrational wavenumbers for the vacuum-isolated monohydrates of formaldehyde, tetrahydrofuran, pyridine, tetrahydrothiophene, trifluoroethanol, methyl lactate, dimethylimidazolidinone, cyclooctanone, trifluoroacetophenone and 1-phenylcyclohexane-cis-1,2-diol are provided. The results of the challenge show promising predictive properties in both purely quantum mechanical approaches as well as regression and other machine learning strategies.

Understanding the high-resolution spectral signature of the N2-H2O van der Waals complex in the 2OH stretch region

We present the observation of the N2-H2O van der Waals complex in the 2OH stretch overtone region. The high-resolution jet cooled spectra were measured using a sensitive continuous wave cavity ringdown spectrometer. Several bands were observed and vibrationally assigned in terms of ν1, ν2, and ν3, the vibrational quantum numbers of the isolated H2O molecule, as (ν1'ν2'ν3')←(ν1″ν2″ν3″)=(200)←(000) and (101) ← (000). A combination band involving the excitation of the in-plane bending motion of N2 and the (101) vibration of water is also reported. The spectra were analyzed using a set of four asymmetric top rotors, each associated with a nuclear spin isomer. Several local perturbations of the (101) vibrational state were observed. These perturbations were assigned to the presence of the nearby (200) vibrational state and to the combination of (200) with intermolecular modes.