Toward Elimination of Discrepancies between Theory and Experiment: Anharmonic Rotational-Vibrational Spectrum of Water in Solid Noble Gas Matrices.

Rotational-vibrational spectroscopy of water in solid noble gas matrices has been studied for many decades. Despite that, discrepancies persist in the literature about the assignment of specific bands. We tackle the involved rotational-vibrational spectrum of the water isotopologues H216O, HD16O, and D216O with an unprecedented combination of experimental high-resolution matrix isolation infrared (MI-IR) spectroscopy and computational anharmonic vibrational spectroscopy by vibrational configuration interaction (VCI) on high-level ab initio potential energy surfaces. With VCI, the average deviation to gas-phase experiments is reduced from >100 to ≈1 cm-1 when compared to harmonic vibrational spectra. Discrepancies between MI-IR and VCI spectra are identified as matrix effects rather than missing anharmonicity in the theoretical approach. Matrix effects are small in Ne (≈1.5 cm-1) and a bit larger in Ar (≈10 cm-1). Controversial assignments in Ne MI-IR spectra are resolved, for example, concerning the ν3 triad in HDO. We identify new transitions, for example, the ν2 101 ← 110 transition in D2O and H2O or the ν3 000 ← 101 transition in D2O, and reassign bands, for example, the band at 3718.9 cm-1 that is newly assigned as the 110 ← 111 transition. The identification and solution of discrepancies for a well-studied benchmark system such as water prove the importance of an iterative and one-hand combination of theory and experiment in the field of high-resolution infrared spectroscopy of single molecules. As the computational costs involved in the VCI approach are reasonably low, such combined experimental/theoretical studies can be extended to molecules larger than triatomics.

Introduction

Water in the gas phase is of key importance in several fields, including engineering, meteorology, astrochemistry, and physical chemistry. Water molecules are abundant in Earth’s atmosphere and play a crucial role in cloud formation, the hydrologic cycle, and climate. Water is an important greenhouse gas, closing large parts of the IR window of the atmosphere to space. It efficiently absorbs thermal radiation from Earth’s surface but also from incoming sunlight. Without atmospheric water, Earth would be in a permanent ice age. Its identification in the atmosphere and in space is possible by its spectroscopic properties, based on rotational–vibrational spectroscopy and on microwave spectroscopy.[1,2] Pure rotational spectra[3,4] and IR spectra[5,6] are also used for remote detection. Water is a simple, triatomic, asymmetric rotator with C2 symmetry, yet its spectra are highly complex. In the gas-phase, symmetric stretch (ν1), asymmetric stretch (ν3), and bending (ν2) of the covalent bonds are coupled to rotations. Combination bands and overtones complicate the picture even more. The spectroscopic complexity of water in the gas-phase is reflected in a comprehensive summary by the IUPAC.[6−9] In 1957, Pimentel et al. isolated water in an inert N2 matrix at 20 K,[10] providing the first matrix isolation infrared (MI-IR) spectrum. After this pioneering study, single molecule spectroscopy of water has seen a bloom[11,12,21,13−20] where MI has emerged as the standard for the study of the water monomer, the water dimer, and higher oligomers. In the atmosphere, the identification of the dimer has been a challenge as there is a strong overlap of monomer and dimer bands. It was only possible in 2003 through IR spectroscopy[22] (for a review, see Mukhopadhyay et al.[23,24]). In contrast to gas-phase IR, MI-IR spectroscopy allows for a dissection of the water monomer and its oligomers. However, the rotation of the water molecule in a cage of noble gas atoms and matrix splitting of bands caused by different cage geometries have been a continuous source of debate and reassignments. Also, complexes of water, for example, the H2O*N2 cluster, have been considered to assign all observed bands.[25,26] Water is one of the most studied molecules in MI-IR, with experimental insights available for a wide range of different conditions, for example, different matrix materials, different temperatures, different isotopes, or in the presence of added dopants.[11,12,25−34,13−15,17−21] A recent historical overview is given by Ceponkus et al.[32] For the Ar matrix, the rotational–vibrational assignments for H216O, HD16O, and D216O are summed up by Engdahl and Nelander,[33] and a more recent study of H216O is given by Perchard.[27] For the Ne matrix, a comprehensive assignment is reported by Forney et al.[28]

The rotational–vibrational spectrum is in theory derived by solving the nuclear Schrödinger equation. For a simple triatomic system such as water, sophisticated theoretical predictions are available,[35−38] with global potential energy surfaces (PES) calculated at a very high level of electronic structure theory and empirical refinement, while the exact solution of the nuclear Schrödinger equation provides a full line list of rotational–vibrational transitions that agree with gas-phase spectra within hundredths of a wavenumber. Such approaches yield the most accurate computational results possible and are strongly required, for example, in the field of astrochemistry. Considering water in noble gas matrices, tailor-made exact models were reported,[32] too. However, exact approaches are tedious to implement, and thus, they are not necessarily flexible, that is, not universally applicable but specific for the molecule under study. In order to be able to understand single molecule spectra in the matrix, a compromise needs to be found that provides sufficient accuracy on the theory side to allow for a sound assignment of experimentally observed bands. Ideally, the computation is efficient enough so that the method is suitable for a variety of molecules, also larger than triatomics. That is, computational approaches are needed that combine efficiency and accuracy yet also rigor and flexibility. When it comes to molecular identification by infrared spectroscopy, it is common to rely on computed vibrational spectra that are routinely predicted within the harmonic approximation (HA). Its inherent theoretical shortcomings are well known.[39] The most prominent drawback of the HA is its inaccuracy-vibrational frequencies are generally overestimated, with an error bar larger than the spectral range of the vibrational transitions. In order to compare HA spectra with experimental spectra, the easiest yet least rigorous practice is to shift the entire calculated HA spectrum with respect to the most intense band in the experimental spectrum by a specific scaling factor taken from the literature. As anharmonicity is specific for each vibration, such a linear scaling fails to correct for the whole calculated spectrum. The approach of scaled quantum mechanics[40] provides a theoretically more rigorous route to improve the accuracy, while flexibility is somewhat reduced, as the improvement is based on empirical refinement and, thus, biased toward an underlying experimental data set. In theoretical spectroscopy, two rigorous models that go beyond the HA received particular attention: the vibrational perturbation theory (VPT2), as used within MI studies of water by Tremblay et al.,[30,31] and the vibrational self-consistent field (VSCF)[41−43] in combination with vibrational configuration interaction (VCI) based on multimode representations of the PES. A correlation-corrected VSCF approach was demonstrated, for example, by Jung and Gerber for water oligomers.[44] Both approaches come with improved accuracy and retain flexibility and efficiency. They are comprehensively reviewed.[45,46] In the present work, we employ the VSCF/VCI approach by Rauhut et al.,[47−50] which is concisely described in the next section.

To this end, we not only perform VCI calculations ourselves but also report new MI-IR experiments of H216O, HD16O, and D216O. As matrix materials, we use argon, which is commonly used in the literature,[11,12,25,27,13−15,17−21] and neon, which is less frequently used,[28−31] as its low freezing point is technically challenging to handle. Our cutting-edge experimental MI-IR setup, with a vacuumed spectrometer, a high-quality light source and interferometer, and a liquid nitrogen cooled detector, allows for a very good signal/noise ratio that helps to detect even the weakest bands (see next section for more details). It should be stressed that water is a well-studied molecule. As mentioned above, state-of-the-art calculations provide extensive and highly accurate rotational–vibrational line lists. Those can be employed to investigate the MI-IR spectrum of water. However, with the one-hand combination of experiment and calculation in this work, we are in a position to perform an iterative assignment. We do not need to rely merely on literature data. The merit of such a combined approach becomes clear under the consideration that it can be extended to molecules where computationally predicted spectroscopic properties are not as accurately accessible as for the water molecule.

Methodology

Experimental Spectrum by Matrix Isolation Infrared Spectroscopy

The matrix isolation (MI) technique used in this study has been previously described by Grothe et al.[51] A detailed depiction of the apparatus is found in the Supporting Information of ref (52). The highlights of the setup are a 180° rotatable high-vacuum cryostat that can be cooled to 5.8 K, allowing for stable solid Ne matrices, and the FTIR spectrometer Vertex 80v (Bruker, Karlsruhe), which is operating under 2 mbar with a liquid-nitrogen-cooled MCT detector to eliminate all atmospheric noise in the infrared region. This is especially necessary when studying the infrared spectrum of the omnipresent water. The spectrometer itself is vibration-cushioned and uncoupled from the matrix isolation apparatus to avoid vibration-induction. The arising gap is constantly flushed with pure N2 gas. In all experiments, the atmospheric water is not observed with this setup. To record the spectra of water vapor trapped in the noble gas matrix, 512 scans are accumulated at a resolution of 0.3 cm–1.

We use a liquid sample of D2O (99.96%, delivered by Eurisotop, D215ES, lot = Q2141, batch = BR575) in all experiments. H2O and HDO species are generated spontaneously by isotope exchange with protons on the surface of the mixing chamber steel tubes. About 2.5 mL of the sample is added by a sterile syringe into a glass flask under argon flux to reduce contamination with atmospherically gases such as H2O, CO2, O2, and N2. The sample is mounted to the mixing chamber of the MI apparatus and herein degassed by multiple freeze and thaw procedures. Before mixing of the sample with a rare gas, the vapor above the liquid is evacuated in order to provide “freshly” evaporated chemicals. The gas mixture is then prepared by barometric monitoring within a volume of about 2 L that was previously evacuated at 10–5 mbar. A mixture of 2 mbar sample gas in dilution with about 990–1010 mbar rare gas is denoted as 1:500. We measure the spectra of 1:250, 1:500, and 1:1000 dilutions both in Ar and in Ne. The matrix is deposited onto a gold plate in the high-vacuum chamber at 10–7 mbar and 5.8 K, generated by a Gifford-McMahon cooler.[53] The gas mixture is introduced at a constant flow of 4 mbar/min from a volume of about 200 mL and a pressure of 900–980 mbar, depending on the previously prepared mixture. The deposition time is about 45 min per matrix layer. As water adsorbs on the surface of the mixing chamber steel tubes, the initial ratio of dilution is somewhat altered. Hence, it is necessary to deposit up to three layers to produce a matrix containing enough water species for an intensive infrared spectrum. All IR spectra are recorded in the mid-infrared region of 4000–500 cm–1. Background spectra of the gold plate at 5.8 K are taken before each sample measurement. The spectra shown in Figure are corrected for this background.

Figure 2

MI-IR spectrum of a D2O sample containing H2O and HDO, diluted 1:250, 1:500, and 1:1000 in (a) Ar and (b) Ne, both at 6 K. Pure vibrational transitions are labeled as ν along the vertical dotted lines. Rotational–vibrational transitions are labeled as J′ ← J″, with J being the rotational quantum number and Ka and Kc being the projection quantum numbers of the asymmetric top. *N2 are supposedly due to a H2O*N2 complex. Water dimers are labeled as proton donor (PD) or proton acceptor (PA) vibrations. New or tentative assignments are marked in green. Expected bands from literature data are marked in gray.

Calculated Spectrum by Vibrational Self-Consistent Field and Configuration Interaction

All calculations are performed with the MOLPRO software package.[54,55] In all that follows, the electronic structure is calculated with the coupled cluster ansatz including single, double, and perturbative triple excitations and the explicit correlation of interelectronic distances, that is, the CCSD(T)-F12 ansatz.[56,57] This ansatz allows for a faster basis-set convergence than conventional coupled cluster approaches with little additional computational cost. The accuracy of this method with the here used triple-ζ (VTZ-F12) basis set[58,59] corresponds roughly to a conventional calculation with a quintuple-ζ basis set, which is sufficiently close to the complete basis-set limit. The Born–Oppenheimer equilibrium geometry of the water monomer, numerical second derivatives of the energy near the equilibrium, and the harmonic frequencies are obtained by using the MOLPRO default procedures optg and freq. This yields also the commonly known normal modes of the water molecule: asymmetric stretch q3, symmetric stretch q1, and in-plane bend q2. These normal-modes q are used as a coordinate system for the subsequent calculation of a multimode representation of the potential energy surface (PES), a sum of one-mode (1D) potentials V(q), two-mode (2D) potentials V(q,q), and so forth, as depicted in Figure . The strongly anharmonic nature, especially of V(q1), is evident. The MOLPRO surf algorithm constructs the PES on a grid where each grid point is generated by elongation along the normal modes (1D), along two combined normal modes (2D), and so on. Symmetry considerations, iterative interpolation, and prescreening techniques significantly lower the amount of ab initio single point energy calculations, with negligible loss of accuracy in the final PES, as pioneered by Rauhut.[47] After testing various schemes for the PES construction, we find in accordance with the literature[60] that the use of the CCSD(T)-F12 level of theory in 1D, 2D, and 3D sub-potentials is preferable in terms of accuracy in the subsequently calculated vibrational spectrum. As the water molecule has a rather simple electronic structure, the CCSD(T)-F12 single point energy calculations are also not too expensive. Based on the multimode representation of the PES, the nuclear Schrödinger equation using the Watson operator is solved variationally within the VSCF[41−43] where the wave function is approximated as a Hartree product of one-mode functions. This approximated solution is subsequently correlation-corrected by the configuration-selective VCI.[48,50] This yields the VCI wave function as a linear combination of various excited Hartree products, with the anharmonic vibrational energies corresponding to specific VCI states that inherit the nomenclature of the normal-mode picture. The nuclear wave function allows for the calculation of IR intensities by employing a multimode representation of the dipole moment surface. The latter is calculated at the HF/VTZ-F12 level of theory. The rotational fine structure is calculated within the adiabatic rotation approximation (ARA).[61] All these calculations are achieved by the vscf and vci algorithms that are provided in MOLPRO. In a fairly rough approximation, the R-Branch is considered as the mirror image of the P-Branch in close vicinity of the Q-Branch. In order to assign rotational–vibrational transitions from the initial state J″ to the excited state J′, we use the spectroscopic nomenclature J′ ← J″, with J being the rotational quantum number and Ka and Kc being the projection quantum numbers of the asymmetric top.

Figure 1

Potential energy surface of H2O. Grid representation of the multimode PES, calculated with CCSD(T)-F12/VTZ-F12. (left) One-mode sub-potentials V(q) with vector depiction of the normal-modes q of water. (right) Two-mode sub-potentials V(q,q). Three-mode sub-potentials are not depicted yet considered in the calculation.

Results and Discussion

The results of our assignment in the MI-IR spectra are depicted in Figure a considering the Ne matrix and in Figure b considering the Ar matrix. The assignment is summed up in Tables and . With a pure D216O sample, we simultaneously observe H216O, HD16O, and D216O in the spectrum due to proton exchange in the mixing chamber. In Ar, the H2O shows the weakest bands. D2O and HDO show higher intensities, inversely proportional to the sample dilution. The same holds for Ne, yet here, the intensities are lower overall. In both Ar and Ne MI-IR, we observe roughly 10 bands that can be accounted to water dimers (see Table S1 in the Supporting Information). All spectra exhibit multiple rotational–vibrational transitions (rovibs) with high intensities in the P-branch and a few less intense rovibs in the R-branch. In Ar, most of the pure vibrational transitions (vibs) are observed, whereas in Ne only two are observed, presumably because of the low overall intensity. Table presents a comparison of the vibs between the HA and VCI calculations and experimental data (comprising the measured Ar and Ne MI-IR data here, some literature MI-IR data,[17,28,32,33] and literature gas-phase data.[62−66]) Table presents a similar comparison, considering the rovibs. When compared to the gas-phase, it is immediately evident that the HA is far from being sufficient to reproduce the experiment. Considering the vibs, the mean absolute deviation (MAD) between HA and experiment amounts to 107.9 cm–1, with a maximum deviation (MAX) of 187.9 cm–1. In contrast, the VCI results perfectly agree with the gas-phase, with an MAD of 1.1 cm–1 (MAX = 1.6 cm–1). That is, the VCI approach results in a roughly 100-fold reduction of the discrepancy between experiment and theory. This has also been found for other molecules in the literature.[60] Also, when considering the 42 rovibs as listed in Table , the VCI results perfectly agree with the gas-phase, with an MAD of 1.1 cm–1 (MAX = 5.7 cm–1). In view of the assignment of MI-IR spectra, we compare the VCI results with MI-IR data. For the vibs listed in Table , we observe with respect to the VCI results in Ar matrices an MAD of 12.9 cm–1 (MAX of 20.4 cm–1) and in Ne matrices an MAD of 2.1 cm–1 (MAX = 5.2 cm–1). The gas-phase and Ne matrix show similar MADs, whereas the MAD of the Ar matrix increases by a factor of 6. Similar conclusions can be drawn for the rovibs listed in Table . With respect to the VCI results, in Ar matrices, an MAD of 13.1 cm–1 (MAX = 24.4 cm–1) and in Ne matrices an MAD of 3.2 cm–1 (MAX = 7.9 cm–1) are found. In this case, an improvement of a factor of 4 is reached when using Ne instead of Ar. We can confirm that Ne matrices have little effect on the trapped water molecules, therefore being an excellent approximation to the gas-phase situation, whereas Ar does exert a vivid effect and is only a good approximation to the gas-phase situation. The mean matrix-shift incurred by changing from gas-phase (i.e., MAD wrt gas-phase) to Ar matrix of 13.2 cm–1 (rovibs: 13.6 cm–1) and to Ne matrix of 1.5 cm–1 (rovibs: 3.2 cm–1) is much smaller than the improvement gained by anharmonic VCI calculations compared to the HA, which is in mean roughly at 100 cm–1. This means that matrix-shifts, especially in the case of Ne matrices, are not problematic when the assignment is supported by VCI calculations. In other words, the correction for anharmonicity by the VCI approach, with an MAD of 1.1 cm–1 with respect to the gas-phase, perfectly allows identification of matrix-shifts. In all cases, the largest discrepancies between theory and experiments (i.e., the MAX values) arise from the calculated rovibs within the adiabatic rotation approximation (ARA). The H2O molecule as an asymmetric top is known to be an unfavorable case for this method.[61] This holds especially with increasing quantum number J. However, as the here observed rovibs are J = 1 and J = 2, the results are still adequate. In general, this discrepancy is a minor issue for the investigation in MI-IR spectra, as this technique usually suppresses rotation completely. In fact, the water molecule is an exceptional case.

Table 1

Fundamental Vibrational Transitions of Water Monomers as Observed in Ar and Ne MI-IR at 6 K and Compared to Gas-phase Reference and the here Calculated Transitionsa

Values in cm–1. Bold face: obtained in the present study. Italics: Reference data.

VCI(SDTQ) on a 3-mode PES at CCSD(T)-F12/VTZ-F12.

Harmonic Approximation at CCSD(T)-F12/VTZ-F12.

Ref (17).

Ref (28).

Ref (32).

Ref (33).

Ref (62).

Ref (63).

Ref (64).

Ref (65).

Ref (66).

Table 2

Rotational-Vibrational Transitions of Water Monomers as Observed in Ar and Ne MI-IR at 6 K and Compared to Gas-phase Reference and the here Calculated Transitionsa

Bold face: obtained in the present study. Italics: Reference data. In green: tentative or new assignments.

VCI(SDTQ) on a 3-mode PES at CCSD(T)-F12/VTZ-F12.

Harmonic Approximation at CCSD(T)-F12/VTZ-F12.

Ref (17).

Ref (28).

Ref (32).

Ref (33).

Ref (62).

Ref (63).

Ref (64).

Ref (65).

Ref (66).

With the use of the VCI approach, we can confirm most of the transitions assigned in the literature and additionally identify bands that were previously not assigned. In Ne MI-IR, a band at 1580.0 cm–1 can be assigned as H2O ν2 101 ← 110 (VCI: 1578.9 cm–1) and another band at 1170.2 cm–1 can be assigned as D2O ν2 101 ← 110 (VCI: 1169.1 cm–1). As the underlying normal-mode for both species have the same A1 irrep, this assignment is also plausible from symmetry considerations. In previous Ne experiments,[28] the H2O ν3 000 ← 101 transition is observed, yet the corresponding transition in D2O is not. We find a weak signal at 2778.7 cm–1 that can be assigned as the D2O ν3 000 ← 101 (VCI: 2773.9 cm–1). Furthermore, an inconsistency in the literature considering the HDO ν3, ν1, and ν2 can be resolved. The HDO normal-modes have the same A′ irrep, and hence, we expect a similar band shape and assignment for all three spectral regions. In the HDO ν1 and ν2 region, the assignments in Ar MI-IR by Engdahl and Nelander[33] and in Ne MI-IR by Forney et al.[28] are consistent and also plausible from the symmetry considerations. In the HDO ν3 region, we observe a triad with an intensity of 1:5:1 in both Ar and Ne matrices. The assignment of these transitions by Engdahl and Nelander in Ar as 111 ← 000, 101 ← 000, and HDO*N2 differs from the assignment by Forney et al. in Ne as 212 ← 101, 111 ← 000, and 110 ← 101. Our VCI calculations confirm the assignment given by Engdahl and Nelander, which is consistent within all three regions also from symmetry consideration. We adopt the assignment of the HDO ν3 region by Engdahl and Nelander also for the Ne MI-IR. One intriguing point still remains: The band at 3690.6 cm–1 in Ar is given by Engdahl and Nelander[33] either as the HDO*N2 band or as the 110 ← 111 transition. In the Ne matrix, this band occurs at 3718.9 cm–1 where it is assigned as 111 ← 000 by Forney et al. Based on our calculations, the 110 ← 111 transition is at the origin of this band and it is not necessary to invoke a N2-induced band. Other discrepancies in the assignment of the HDO species are noticeable between Forney et al.[28] and a more recent publication by Ceponkus et al.[32] The HDO ν3 is assigned at 3699.0 cm–1 in ref (28) and at 3713.8 cm–1 in ref (32). We observe both bands with almost identical but very weak intensities. Considering the VCI calculation, the band should in fact be located at neither of these two positions but in the middle between them (VCI: 3706.1 cm–1). Based on the good agreement between VCI calculation and the residual Ne MI-IR spectrum, we suggest that HDO ν3 should be at roughly 3707 cm–1. This is not observed in our experiments and not reported in the experiments in refs (32) and (28). The band either red- or blueshifts by 7 cm–1. In our experiment, it seems to basically split by 14 cm–1. VCI-calculated (and gas-phase) rovibs cannot explain the bands at 3699.0 and 3713.8 cm–1. Further experiments may verify whether the HDO ν3 fundamental vibration can be found at 3707 cm–1. The HDO ν1 is assigned at 2722.9 cm–1 in ref (28), which is in excellent agreement with our VCI calculations. Ceponkus et al.[32] claim it to be at 2727.5 cm–1. This band, however, is astonishingly close to a transition observed in ref (28) for HD18O species (2727.7 cm–1). In our experiment, HD18O can be clearly ruled out. If naturally occurring 18O was the origin, then we should observe a 1000-fold more intense HD16O band, which is not the case. Also, rovibs can be excluded based on our VCI calculations and gas-phase reference. As shown in Figure , there are some bands in the 2729–2727 cm–1 region and a weak band at 2718.9 cm–1, which is indirectly proportional to the sample concentration. The observed bands must have a different origin and do not represent fundamental vibs or rovibs of the water monomer. The origin of this remains unclear—it might be explained based on oligomers or effects not considered in our calculation, for example, a splitting caused by different kinds of matrix cages.

Besides the improved spectroscopic assignment, the VCI wave function allows for the computation of vibrationally averaged structural parameters,[67,68] here collectively abbreviated as rAv and not to be confused with equilibrium structural parameters, abbreviated as rE, from geometry optimization within the Born–Oppenheimer approximation. In Table , rAv and rE from the present multimode PES-based VSCF/VCI calculation and from the extensive calculations by Czakó et al.[68] on an adiabatic semiglobal PES[69] are listed together with structural parameters deduced from microwave (MW) spectroscopy[70] and gas-phase electron diffraction (GED) experiments.[71] The experimental structural parameters from MW, for example, the rotational constants A, B, and C for all three isotopomers, are in a much better agreement with the computed rAv than with the computed rE, for both the calculations by Czakó et al.[68] and the present VCI calculations. The rE bond length and angle are the same for all three isotopomers, as the Born–Oppenheimer PES is mass-independent. In experiments, however, bond lengths and angles are not observed as mass-independent. Both in MW and in GED experiments, a difference of roughly 0.06 Å between r(OH) and r(OD) is observed. This experimentally observed mass dependency is also reproduced in the calculated rAv bond lengths. Although a comparison of ab initio computed with spectroscopically deduced structural parameters is involved—we refer to Császár et al.[72] for a discussion on the pitfalls of such comparison—the good agreement substantiates the physical validity of the present VCI approach and it is clear evidence that anharmonicity significantly impacts the observed molecular structure.

Table 3

Calculated and Experimental Gas-phase Structural Parameters of Watera

Bold face: obtained in the present study. Italics: Reference data.

Geometry optimization at CCSD(T)-F12/VTZ-F12.

Vibrational Ground State from VCI(SDTQ) on a 3-mode PES at CCSD(T)-F12/VTZ-F12.

Calculations by Czakó et al.

eMicrowave experiments in Ref. (70)

Gas-phase electron diffraction experiments in Ref. (71).

Conclusions

Our results demonstrate that deviations of computed vibrational spectra to experimental spectra are minimized when anharmonicity and mode coupling are incorporated with a multimode PES representation and the time-independent nuclear Schrödinger equation is solved with the variational VSCF/VCI approach. This approach allows for a prediction of pure vibrational transitions with an error of just 1 cm–1 on average and is also suitable for describing the challenging situation of a molecule rotating in the matrix cage when the rotation is not too hindered as in the case of water. Only this accuracy allows us to confidently assign every single observed band, identify inconsistencies in the state-of-the-art assignment, and to reassign some bands in the Ne MI-IR spectrum of the HDO species. From our MI-IR experiments, we gain profound understanding of matrix-shifts where Ne shifts are rather small with a difference to the gas phase of only 1.5 cm–1 on average. With Ne MI-IR being close to the gas-phase situation and the capability of the VSCF/VCI approach to accurately predict gas-phase transitions, we encounter a powerful combination of experiment and theory. Even for water that has been studied for more than 60 years by MI-IR spectroscopy and theoretical calculations, our combined approach has allowed us to improve our understanding of the rotational–vibrational transitions. The theoretical framework and the way to include anharmonicity as employed here will allow for sound assignments of many small molecules studied by MI-IR spectroscopy. The benefit will become even more evident for molecules where exact solutions are not yet available in the literature.