
==== Front
J Phys Chem A
J Phys Chem A
jx
jpcafh
The Journal of Physical Chemistry. a
1089-5639
1520-5215
American Chemical Society

39155731
10.1021/acs.jpca.4c03468
Article
Regularities and Anomalies in Neon Matrix Shifts of Hydrogen-Bonded O–H Stretching Fundamentals
https://orcid.org/0000-0002-5883-3801
Bödecker Margarethe †§
https://orcid.org/0000-0002-2594-7028
Mihrin Dmytro ‡§
https://orcid.org/0000-0001-8841-7705
Suhm Martin A. *†
https://orcid.org/0000-0003-2983-6795
Wugt Larsen René *‡
† Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
‡ Department of Chemistry, Technical University of Denmark, Kemitorvet 206, 2800 Kgs. Lyngby, Denmark
* Email: msuhm@gwdg.de.
* Email: rewl@kemi.dtu.dk.
19 08 2024
29 08 2024
128 34 71247136
26 05 2024
06 08 2024
05 08 2024
© 2024 The Authors. Published by American Chemical Society
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/).

O–H bond stretching vibrations in hydrogen-bonded complexes embedded into cryogenic neon matrices are subtly downshifted from cold gas phase reference wavenumbers. To the extent that this shift is systematic, it enables neon matrices as more universally applicable spectroscopic benchmark environments for quantum chemical predictions. Outliers are indicative of either an assignment problem in one of the two cryogenic experiments or they reveal interesting dynamics or structural effects on the complexes as a function of the environment. We compile 6 literature-known pairs of experimental data in jet and neon matrix expansions and realize a 6-fold expansion of that number through targeted matrix isolation and/or slit jet expansion spectroscopy presented in this work. In many cases, the neon matrix shift is less than the uncertainty of the currently best-performing blind quantum chemical predictions for the gas phase, but in specific cases, it may exceed the currently achievable theoretical accuracy. Some evidence for a positive correlation of the matrix shift with the hydrogen bond shift is found, similar to observations for helium nanodroplets. Outliers in particular for water acting as a donor are discussed, and in a few cases they call for a future reinvestigation. Substantial improvement in the correlation of the matrix shift with the hydrogen bond shift is achieved for ketone monohydrates by removing a vibrational resonance. New insights into nitrile hydration isomerism are obtained, and the linear OH stretching spectrum of the jet-cooled ammonia–water complex is presented for the first time. Vibrational spectroscopy in weakly perturbing solid rare gas quantum matrices for the benchmarking of gas phase theory and future explicit theoretical treatments of the quantum matrix environment to better understand the outliers are both encouraged.

Deutsche Forschungsgemeinschaft 10.13039/501100001659 389479699 / GRK2455 document-id-old-9jp4c03468
document-id-new-14jp4c03468
ccc-price
Special Issue

Published as part of The Journal of Physical Chemistry Aspecial issue “Richard J. Saykally Festschrift”.
==== Body
pmc1 Introduction

For the empirical performance test of state-of-the-art electronic structure and quantum nuclear dynamics predictions about hydrogen-bonded complexes,1 experimental gas phase studies in cold supersonic jets2 and traps3 are the gold standard. However, they can suffer from restrictions such as sensitivity, indirect spectral probing, or metastable conformations.4 Control or even replacement experiments that embed the complex in an inert matrix are highly desirable because they can reduce some of those restrictions by longer observation times, by strictly linear absorption features, and by annealing experiments.5,6 The associated spectral matrix interaction shift is not easy to model7−14 but often minimized in quantum matrices built from light species such as H2, He, or Ne. In contrast, Ar matrices are known for inducing substantial conformational switches, such as in the case of ethanol dimers.15,16 Although there are more than a dozen research groups in the field of neon matrix isolation spectroscopy of hydrogen bonding and even more working on corresponding supersonic jet characterizations, there is a surprisingly low number of published OH stretching data pairs which would allow us to judge the influence of the matrix. Here, we extend previous coordinated jet plus Ne-matrix experiments17−20 and increase the number of available data pairs from 6 to 36 in three ways: by experimentally complementing single literature values with their missing matrix or cold gas phase counterpart (5 cases); by making available some previously unpublished reference spectra obtained in our groups (22 cases); and by targeting a few new systems with both experimental techniques (3 cases). The size of the expected Ne matrix downshifts is on the order of 10 cm–1 for the particularly environment-sensitive O–H stretches, perhaps slightly increasing with increasing hydrogen bond strength. This is comparable to rare residual site splittings in this spectral range, still well below 1% of the absolute transition wavenumber, and it is roughly within explicit anharmonic contributions to O–H stretching fundamentals which cannot be accounted for by simple computational harmonic scaling approaches.2,21

New jet/Ne-matrix data pairs and critical analysis of the published literature can lead to a win-win situation: In regular cases, the data pairs allow to identify any systematic trends in such Ne-matrix shifts and thus widen the benchmark potential of spectroscopic experiments toward quantum theory in the field of hydrogen bonding with its mix of intermolecular interactions. In anomalous cases, they may hint at potential spectral interpretation issues on the gas phase or matrix isolation side.22,23 Otherwise, they reveal interesting matrix embedding effects on the complex structure and vibrational dynamics, such as conformational switches, tuned anharmonic resonances,24 possibly also constrained large amplitude motion or weak cooperative hydrogen bonding to matrix atoms. Indeed, studies of complexes with more classical matrix species M such as N2 have previously revealed sizable XH spectral shifts,25,26 whenever an XH···XH···M motif is possible. The popular cryogenic tagging technique, which brings action spectroscopy of hydrogen-bonded complexes closer to direct absorption measurements, profits from sufficiently small shifts induced by the tag, which is often chosen to be a typical (quantum) matrix gas.27 In rotational spectroscopy, neon complexes are regularly observed and characterized as side products, but definitely worth a closer investigation, whenever they show interesting structures or dynamics.28 More often, complexes with heavier rare gas atoms are observed and discussed, and systematic trends across group 18 of the periodic table can be particularly informative.29 Therefore, a quantitative understanding of spectral shifts induced by single, multiple, or infinite numbers of rare gas atoms attached to a molecular system has multiple benefits.

After a brief description of the employed experimental techniques and uncertainty assessments, we present an initial table of dominant peak wavenumbers in the gas and Ne matrix phases, categorized according to the spectral technique and resulting in formal band position shifts due to the matrix environment. In the next step, several of these shifts are critically assessed with possible resonance partners or other spectral complications in mind, thus leading to a cautious adjustment and uncertainty estimate of the underlying O–H oscillator position and its shift in the neon matrix. We discuss selected outliers as a function of hydrogen bond strength and type, suggest a reinvestigation in some cases, and end with a bold toy model proposal to trigger theoretical and experimental extensions of this work.

Given that O–H stretches are particularly sensitive to tagging30 and matrix effects,26 the present study probably comes close to the worst-case scenario for the use of vibrational matrix data in comparison with gas phase computations. Only acidic X–H stretches are expected to be more sensitive to the matrix environment.31−33 We thus conclude that Ne matrix isolation spectroscopy has a substantial benchmark potential34 for quantum chemistry and dynamics beyond the hydride stretching region, in particular for the class of large-amplitude intermolecular modes generated by hydrogen bonding, which are observed in the experimentally challenging far-infrared spectral region.35−37

2 Experimental Section

Direct infrared absorption by molecular complexes in the cold gas phase is a powerful technique to detect the position of hydride stretching fundamentals, often in combination with long slit jet expansions.38,39 Here, we employ the most recent variant of an FTIR-probing approach40 which uses gas recycling to make expensive compounds and carrier gases accessible despite a high gas throughput. Relevant parameters for the present work include a spectral resolution of 2 cm–1 (in explicitly denoted cases 1 cm–1 or less), the use of an InSb detector together with an optical filter to address OH stretching fundamentals, the application of a gas stagnation pressure of less than 1 bar, and the use of typical dilutions of hydrogen bond acceptor and donor molecules in a lower than 1:100 ratio to the carrier gas (for specific values see the specific spectral figures, also in the Supporting Information). The pulsed admission of such gas mixtures to a vacuum container through the Bruker VERTEX 70v FTIR-probed 700 mm slit nozzle far exceeds the available pumping capacity of 2000 m3/h. This is compensated for by a 4 m3 vacuum buffer and a sufficiently long waiting time between pulses. Such a brute-force gas-phase approach is still orders of magnitude less sensitive than matrix isolation, which is again combined with FTIR spectroscopy in the present work. This experimental setup employs a 4K liquid helium-free closed-cycle cryocooler from Advanced Research Systems interfaced with a Bruker VERTEX 80v FTIR spectrometer.37 In brief, doped matrices of neon are prepared by mass flow-controlled deposition of LN2 precooled neon gas onto a wedged diamond window housed in an oxygen-free high-conductivity copper window holder attached to the cold head. The outer rotatable vacuum shroud of the cryocooler is equipped with two additional wedged diamond windows, providing an optical port for combined infrared and terahertz investigations. An independent dual inlet copper tubing system enables both the codeposition of selectively isotopically enriched samples and spatial “spot-to-spot” investigations of solute-rich and water-rich regions of the doped neon matrices.41,42 A combination of resistive heaters, Si diode temperature sensors, and feedback electronics furthermore promotes the option to anneal the doped neon matrices up to 9.5 K. This annealing procedure both circumvents rare site effects and triggers further complexation events in the soft neon environment. A spectral resolution of 0.6 cm–1 is often chosen as the best compromise between the signal-to-noise ratio and sufficient spectral resolution to resolve observed band structures. For the neon matrix isolation spectra, the mixing ratios are less important than the spatial distribution of the dopants within the matrix, which is monitored by the 3.5 mm IR probe beam. This is due to the separate deposition inlet tubing approach, which has been designed for partial D-enrichment investigations. The distributions of dimeric and trimeric species differ significantly spatially within the approximately 10 × 10 mm2 circular matrices even for a specific mixing ratio. This allows for the easy distinction between monomeric, dimeric, and trimeric species within the same experiment. Furthermore, a series of spectra have been collected after annealing to 9.5 K along with independent complementary reference spectra for each of the subunits. The combined data sets are often sufficient to confidently establish the dimer assignments. In both experiments, decadic absorbance spectra are collected by reference to a background spectrum without gas flow or deposited matrix. These absorbances are in the range of 1 for the neon matrix spectra but only 10–4–10–5 in the jet spectra, thus underscoring the desire to replace jet studies with matrix isolation studies for systematic benchmarking work. Used chemicals are listed in Table S4.

Normally, it is sufficient to determine the band maximum of a vibrational transition in either the jet or the matrix spectrum for a neon matrix shift calculation, given the modest asymmetry of the observed bands with a typical width (FWHM) between 2 and 10 cm–1 and the integer wavenumber precision targeted in this work. One would then assign a base uncertainty of ±1 cm–1 to the observed shifts, which may already arise from limiting integer rounding effects. The spectral intensity, which is often harder to quantify, would not be relevant in the shift determination. However, in several cases involving water complexes, there are anharmonic resonances43−45 which may be tuned or detuned by the matrix environment. The nature of these resonances is that a dark state with essentially no intrinsic IR intensity (typically less than 5% of the total intensity involved) steals the intensity from the bright OH stretching state by wave function mixing. In this dark state model limit, the observed intensity is spread over two or more signals, but the position of the original OH bright state can still be obtained in a good approximation by working out the center of the intensity spread over the resonance partners. This conservation of the center of spectral intensity follows from the trace invariance upon diagonalization of an effective normal mode coupling matrix. It may not hold for vibrational Franck–Condon progressions.46,47 Nevertheless, such an effective low-dimensional Hamiltonian model may be the best one can do in the absence of a rigorous anharmonic treatment, as it would be possible for water-dimer-sized complexes.48,49 Beyond simply locating the OH stretching state with the strongest transition observed in the jet or in the matrix (which can jump discontinuously when the interacting states switch their spectral sequence), we thus investigate the shift of the intensity centers of any resonance multiplets wherever these can be clearly identified. For this, we need both the positions and relative intensities of the observed multiplet partners. Our direct absorption spectroscopy techniques are favorable for such an approach. However, there is still considerable uncertainty in the analysis of a real spectral trace due to noise, baseline effects, spectral wings, spectral overlap, monomer and cluster impurities, matrix site splittings, annealing effects, or hot band contributions in the jet due to insufficient vibrational cooling. To minimize the arbitrariness of the proposed resonance multiplet analysis, we used a set of pragmatic rules: 1. Ne matrix and jet spectra for a given system are treated on an equal footing when choosing the width of integration windows

2. Classes of compounds such as ketones are treated on an equal footing

3. Uncertainties from noise, baseline effects, and integration window size are propagated from two limiting bandwidth choices

Details on the signal intensity analysis procedure are described in the Supporting Information (Figure S1, Tables S1 and S2).

In situations of limited signal-to-noise ratio (typical for jet spectra), where resonances may be expected but are not detected, we include the potential effect on the uncertainty of a one-sided (and in this sense worst-case) resonance barely hidden in the noise. Specifically, we assume that this hidden resonance is separated from the main signal by twice the typical resonance coupling element W (W ≈ 10 cm–1 for b2lib resonances,44W ≈ 30 cm–1 for b2ON resonances45 and W ≈ 50 cm–1 for OH bend overtone (b2) resonances45,50). The threshold wavenumbers below which such water-specific resonances may become relevant are chosen as 3550 cm–1 (b2lib), 3450 cm–1 (b2ON) and 3350 cm–1 (b2), respectively. This increases the uncertainty of the intensity center for water complexes compared to the uncertainty of a single absorption band, even in the absence of any other spectral feature. For details, see Table S3 in the Supporting Information. In matrix spectra, the signal-to-noise ratio is typically less limiting than residual site splittings or impurity bands from other aggregates. For literature Ne matrix values of water complexes, we typically assume a rough uncertainty of 10 cm–1 for the intensity center. For Ne matrix spectra obtained in this work, more specific analyses can be carried out based on the stochastic integration procedure.51

An error source that is difficult to quantify involves possible CH combination bands of the acceptor molecule underneath the signals assigned to the OH stretching fundamental and its resonance partners. These could distort the relative intensities and thus the intensity center obtained by deperturbation, in particular, for large residues and high acceptor concentrations. Variation of the donor/acceptor concentration ratio (if possible down to pure acceptor spectra) can help to rule out major distortions from such non–OH signals in the jet and matrix spectra.

3 Results

We start with a compilation of hydrogen-bonded OH stretching fundamentals for the available pairs of (He or Ne) jet-cooled and (Ne) matrix-cooled hydrogen-bonded 1:1 complexes, including both literature data and results reported in this work (TW) for the first time. Where the intensity of the OH stretch is believed to be redistributed among several bands due to resonances, we list only the dominant transition (raw). Table 1 summarizes hydrogen bonds to chalcogen atoms (OH–O(S)) and Table 2 lists OH–N contacts.

Table 1 Raw Wavenumber Downshifts Δ of Dominant O–H Stretching Vibrations in the Ne Matrix Compared to the Gas Phase Values for Different OH–O(S) Donor–acceptor Pairs Using IR-BD (Beam Depletion), CRD (Cavity Ring-down), fs-IR (Femtosecond Ionization Coupled with IR), IR-UV (Double Resonance, Including IR-REMPI Variants), and FTIR Techniquesa

OH–O(S) species	jet experiment	ν̃jet	ν̃mat	matrix experiment	Δ	
w–2,2,2-trifluoroacetophenone	FTIR slit2	3611(1)	3608(1)	Ne matrixTW	3(2)	
w–w	IR-BD52	3601(1)	3591(1)	Ne matrix53	10(2)	
CRD54	3600	 	 	 	
Raman55	3602	 	 	 	
fs-IR56	3601	 	 	 	
FTIR slit44	3602	 	 	 	
w–diacetyl	FTIR slit20	3599(1)	3594(1)	Ne matrix20	5(2)	
w–diacetyl (*)	 	3575(1)	3568(1)	 	7(2)	
w–formaldehyde	FTIR slit2	3591(1)	3585(1)	Ne matrixTW	6(2)	
methanol–methanol	CRD57	3574(1)	3567(1)	Ne matrix58	7(2)	
FTIR slit59	3575	 	 	 	
w–methanol	FTIR slit60	3567(1)	3556(1)	Ne matrixTW	11(2)	
w–cyclobutanone	FTIR slit44	3548(1)	3540(1)	Ne matrixTW	8(2)	
w–ethanol	FTIR slit60	3548(1)	3537(1)	Ne matrixTW	11(2)	
w–oxirane	FTIR slit18	3542(1)	3536(1)	Ne matrix18	6(2)	
w–acetone	FTIR slit44	3538(1)	3515(1)	Ne matrixTW	23(2)	
w–2-propanol	FTIR slitTW	3537(1)	3529(1)	Ne matrixTW	8(2)	
w–acetophenone	FTIR slit44	3536(1)	3513(1)	Ne matrixTW	23(2)	
w–acetophenone (*)	FTIR slit44	3567(1)	3557(1)	Ne matrixTW	10(2)	
w–cyclohexanone	FTIR slit44	3532(1)	3510(1)	Ne matrixTW	22(2)	
ethanol (g)–ethanol (g)(★)	FTIR slit16	3531(1)	3519(1)	Ne matrixTW	12(2)	
Raman61	3532	 	 	 	
w–t-butanol	FTIR slit62	3530(1)	3523(1)	Ne matrixTW	7(2)	
phenol–phenol	IR-UV63	3530(1)	3515(1)	Ne matrixTW	15(2)	
 	 	3507	Ne matrix64	 	
methanol–t-butanol	FTIR slit65	3529(1)	3513(1)	Ne matrixTW	16(2)	
phenol–w	IR-UV66,67	3522(1)	3499(1)	Ne matrixTW	23(2)	
Raman68	3523	 	 	 	
IR-UV69,70	3524	 	 	 	
w–cycloheptanone (★)	FTIR slit44	3512(1)	3504(1)	Ne matrixTW	8(2)	
w–tetrahydrothiophene	FTIR slit2	3507(1)	3499(1)	Ne matrixTW	8(2)	
w–cyclooctanone (★)	FTIR slit2	3503(1)	3493(1)	Ne matrixTW	10(2)	
w–cyclooctanone (*)(★)	FTIR slit2	3525(1)	3505(1)	Ne matrixTW	20(2)	
t-butanol–t-butanol	FTIR slit65	3497(1)	3491(1)	Ne matrixTW	6(2)	
w–tetrahydrofuran	FTIR slit2	3491(1)	3482(1)	Ne matrixTW	9(2)	
a All numbers are in cm–1. TW = This work. w = H2O. Largely sorted in the sequence of decreasing wavenumber without separating metastable isomers (*). Systems, where a conformational switch upon matrix embedding could potentially be happening, are marked with (★). If TW data are not discussed in the main document, the experimental data can be found in Figures S4 and S6–S17 in the Supporting Information.

Table 2 Raw Wavenumber Downshifts Δ of Dominant O–H Stretching Vibrations in the Ne Matrix Compared to the Jet-Isolated Values for Different OH–N Donor–acceptor Pairs Using IR-UV (Double Resonance, Including IR-REMPI Variants) and FTIR Techniquesa

OH–N species	jet experiment	ν̃jet	ν̃mat	matrix experiment	Δ	
w–acrylonitrile	FTIR slitTW	3622(1)	3621(1)	Ne matrix71	1(2)	
w–acrylonitrile (*)	FTIR slitTW	3609(1)	3599(1)	Ne matrix71	10(2)	
w–acetonitrile	FTIR slitTW	3603(1)	3592(1)	Ne matrix53	11(2)	
2-aminoethanol (M)	FTIR slit72	3568(1)	3554(1)	Ne matrixTW	14(2)	
FTIR slit73	3569	 	 	 	
w–ammonia	FTIR slitTW	3486(1)	3456(1)	Ne matrix74	30(2)	
IR-UV75	≈3485	 	 	 	
gas phase76	3515	 	 	 	
w–pyridine	FTIR slit2	3454(1)	3419(1)	Ne matrix77	35(2)	
gas phase76	3480	 	 	 	
w–methylamine	FTIR slitTW	3417(1)	3380(1)	Ne matrix43	37(2)	
methanol–methylamine	FTIR slitTW	3396(1)	3377(1)	Ne matrixTW	19(2)	
methanol–dimethylamine	FTIR slitTW	3339(1)	3323(2)	Ne matrixTW	16(3)	
IR-VUV78	3334	 	 	 	
2-aminoethanol–2-aminoethanol	FTIR slit72	3304(1)	3287(1)	Ne matrixTW	17(2)	
FTIR slit73	≈3300	 	 	 	
a All numbers in cm–1. TW = This work. (M) stands for internally hydrogen-bonded monomers. w = H2O. Sorted in the sequence of decreasing wavenumber without separating metastable isomers (*). The experimental data on methylamine monohydrate and 2-aminoethanol are available from Figures S3 and S5 in the Supporting Information.

The tables include 8 new jet measurements and 24 new neon matrix data to maximize the number and diversity of available pairs. In the final columns of Tables 1 and 2, the resulting matrix downshift in cm–1 is listed together with its estimated error bar. We use integer numbers throughout because of the limited spectral resolution and unresolved band contours in most gas phase data and the dependence on annealing conditions and other heterogeneities in the matrix. We are not aware of any known cases of conformational switching induced by neon matrix embedding. Still, for a few systems, conformational switching upon matrix embedding in the sense of one conformation disappearing and another one appearing upon matrix embedding (both within a monomer and with respect to the contact point between monomers) may be happening (see (★) in Table 1). These cases are discussed in more detail in the Supporting Information (Figures S7 and S10–S12). Future positive evidence for such a neon-matrix-induced conformational switch in a hydrogen-bonded complex would provide a welcome testing ground for any models describing regular matrix effects.

Wherever the infrared OH stretching intensity of a water complex is found or even expected to spread to resonance partners (in 16 cases), we apply an appropriate deperturbation strategy or at least deperturbation uncertainty based on assigned or potential resonance partners44,45 and report the results in Table 3. The deperturbation relies on a negligible intrinsic intensity of the resonance partners, usually overtone or combination bands, and will be exemplified later. It also relies on the assumption that the redistributed intensity arises from OH stretching wave function mixing into dark states, rather than from Franck–Condon-like progressions which do not involve wave function mixing.

Table 3 Deperturbed Wavenumber Downshifts Δ of Water (w) Hydrogen-Bonded O–H Stretching Vibrations in the Ne Matrix Compared to the Cold Gas Phase Values for Different OH–O(S) and OH–N Donor–acceptor Pairs Using FTIR Techniquesa

 	ν̃jet	ν̃mat	Δ	
OH–O(S) species	
w–cyclobutanone	3548(2)44	3540(1)TW	8(3)	
w–ethanol	3548(2)60	3537(1)TW	11(3)	
w–oxirane	3542(1)18	3536(1)18	6(2)	
w–2-propanol	3537(2)TW	3529(1)TW	8(3)	
w–acetone	3532(1)44	3519(1)TW	13(2)	
w–acetophenone	3531(2)44	3522(1)TW	9(3)	
w–t-butanol	3530(2)62	3523(1)TW	7(3)	
w–cyclohexanone	3523(3)44	3514(1)TW	9(4)	
w–cycloheptanone	3515(1)44	3504(1)TW	11(2)	
w–tetrahydrothiophene	3507(3)2	3499(1)TW	8(4)	
w–cyclooctanone	3503(1)2	3493(1)TW	10(2)	
w–cyclooctanone (*)	3525(3)2	3505(1)TW	20(4)	
w–tetrahydrofuran	3491(3)2	3482(1)TW	9(4)	
OH–N species	
w–ammonia	3486(4)TW	3456(1)74	30(5)	
w–pyridine	3454(2)2	3419(10)77	35(12)	
w–methylamine	3417(8)TW	3380(10)43	37(18)	
a All numbers in cm–1. TW = This work. w = H2O. Sorted by decreasing wavenumber but keeping isomers (*) together. A list of all data pairs with water as the donor is given in Table S6 in the Supporting Information.

All of the Ne matrix shifts in Table 1 (OH–O(S)) are downshifts, and most are within 15 cm–1, including their uncertainty. The only outliers are phenol when donating a hydrogen bond to water and a few ketone monohydrates with OH-stretching wavenumbers near 3500 cm–1, which are explained below.

A rather different situation is found for nitrogen as hydrogen bond acceptor, as shown in Table 2 (OH–N). The most prominent matrix downshifts are caused by the monohydrates of pyridine and amines, including ammonia. In contrast, the studied methanol complexes appear to shift significantly less.

The smallest neon downshift of any hydrogen-bonded system, in fact, a negligible one is found for acrylonitrile monohydrate but only for the isomer in which the water coordinates the nitrogen side-on, exposing the nitrogen lone pair to the matrix. Such small or even negative neon matrix downshifts are quite common for the studied OH donor molecules without a hydrogen bond partner. Available values are summarized in the Supporting Information (Table S5).

Figure 1 summarizes the key empirical findings of this work by plotting the neon matrix shifts as a function of the observed jet wavenumber (see Table S6 for the entries). In this figure, deperturbation of the jet or matrix spectra with respect to identified or suspected anharmonic resonances is already taken into account as far as possible. The corresponding figure prior to deperturbation, a deperturbed variant that plots matrix shifts as a function of hydrogen bond-induced downshift rather than absolute jet wavenumber on the abscissa (for intermolecular hydrogen bonds) and one variant only containing literature-known data pairs can be found in the Supporting Information (Figures S18–S20). An alternative would be to plot the matrix shift against the enthalpy of proton transfer.79 In black, some relevant non-hydrogen-bonded monomer results are given for comparison and show small shifts of varying sign. The depicted matrix downshifts for hydrogen-bonded systems (red from Table 1 and blue from Table 2) are uniformly positive and most data points scatter within a 10 cm–1 wide guiding band with a small positive slope of 1/30. This guiding band basically assumes that a hydrogen bond downshift of 30 cm–1 induces an extra matrix downshift of roughly 1 cm–1. Such an empirical, approximately linear relationship has been discussed before for hydrogen-bonded complexes in helium nanodroplets.80 For neon matrices, the relationship needs to be further explored and refined in the future by adding more examples. However, as mentioned when the corresponding data tables are introduced, there are some pronounced outliers from it, mostly involving water. Some of these outliers will be discussed in the following, and several of them are resolved by the deperturbation strategy.

Figure 1 Deperturbed Ne matrix shifts plotted against the hydrogen-bonded OH stretching jet wavenumber for 36 dimers or internally hydrogen bonded monomers (and in black for 4 free monomers that act as donors in the dimers). See Table S6 for the entries.

3.1 Discussion of Outliers

3.1.1 Phenol Monohydration

One of the largest matrix shifts in the table of OH–O(S) (Table 1) bonded systems (red symbols in Figure 1) involves the monohydrate of phenol. In the gas phase, it has been studied several times66−70 and it is consensus that phenol preferentially acts as a hydrogen bond donor toward water. To ensure that this is also the case in the Ne matrix, Figure 2 summarizes a series of matrix experiments, including individual deuteration of the water component only, which would be difficult to realize in a jet experiment. A more generally applicable but expensive alternative is to use oxygen isotopes, which is feasible for both jet65 and cryogenic matrix isolation experiments.81 The isotope-edited matrix isolation experiments for phenol monohydrate show that the assigned OHb signal is indeed due to a phenol vibration. Because this vibration depends in a subtle way on the deuteration level of the accepting water (0, 1, or 2 D atoms), the band shifts slightly and broadens somewhat, but deuterated water as a donor would absorb in the entirely different OD stretching region. The figure shows the absorptions of both hetero- and homodimers and traces of larger clusters together with ro-vibrational transitions of water molecules.82 Compared to the jet spectra, the matrix bands are only somewhat broader and, at best, slightly split. The peak positions and the band integrals are sufficiently well-defined for the purpose of this work, and the signal-to-noise ratio is much better than in typical jet experiments (vide infra). Therefore, any uncertainty in the matrix spectra of the heterodimers is likely to arise from possible spectral overlap with homodimers and larger clusters.

Figure 2 Series of Ne-matrix experiments identifying the phenol OH stretching vibration in its monohydrate among homocluster signals and proving that water is the hydrogen bond acceptor by isotope substitution.

It remains to be shown by appropriate theoretical calculations why the Ne matrix shift of the phenol monohydrate is more than twice as large as that of cycloheptanone monohydrate or t-butyl alcohol monohydrate, where water acts as a donor and the OHb fundamental vibration is similar in wavenumber. However, the intermediate value for phenol dimer (see Table 1, we use our own matrix value rather than the one reported elsewhere64 which deviates systematically, also for monomeric phenol) indicates that phenolic OH groups are more susceptible to neon perturbation than water OH groups. This might be related to their planar nature and needs to be explored more systematically.

3.1.2 Ketone Hydrates

The recently discovered44 systematic resonance between the water OHb fundamental and a three-quantum state involving the water bending overtone (b2) and a librational fundamental (lib) in ketone monohydrates with an effective coupling matrix element W of about 10 cm–1 offers an ideal playground for matrix isolation effects. Depending on the substitution pattern of the ketone, the OHb fundamental will be more or less downshifted, whereas the three-quantum state remains remarkably stationary slightly above 3500 cm–1 in the investigated jet spectra.44 For ketones with an OHb fundamental close to that position, subtle matrix shifts can modify or even invert the resonance situation, formally swapping the order of the zero-order states before engaging in the resonance. Figure 3 gives three examples. In the case of acetophenone, there are actually two hydrogen bond isomers depending on whether water docks on the methyl or on the phenyl side. Phenyl docking is somewhat less stable and also features a less downshifted OHb signal (OHb′). In the jet, there is no appreciable resonance with the three-quantum b2lib state in this case, and the neon matrix downshift is regular. Methyl docking of the water molecule leads to a stronger OHb downshift, still above the dark three-quantum state in the gas phase. Therefore, only part of the OHb intensity is transferred to the dark state, and the higher frequency transition is more intense. This changes in the Ne matrix, which further downshifts the OHb transition, whereas it slightly blue-shifts the three-quantum state (which is plausible for a mode including bending and librational character). That happens to such an extent that the two zero-order states come much closer, intensifying the resonance. They actually cross, and as a consequence, the dominant OHb character (and infrared intensity) is observed at the lower frequency. If one only looks at the dominant transitions, this would suggest an unusual Ne matrix shift of more than 20 cm–1. A simple two-state deperturbation brings the matrix shift of the zero-order OHb positions down below 10 cm–1, into a regular range for OH–O(S) bonded systems. To experimentally locate the zero-order OHb positions, it is sufficient to calculate the center or weighted average of the intensity of the participating transitions, assuming that the three-quantum state brings no intrinsic intensity with it. A related case is found for the monohydrate of cyclohexanone (Figure 3, bottom). Here, the OHb-b2lib resonance is almost perfect, with equal intensity in both transitions within the signal-to-noise ratio of the jet experiment. Each of the transitions has about 50% OHb stretching character. Now, the Ne matrix shifts necessarily detune the resonance and they do so in the same direction as in acetophenone. The three-quantum state is somewhat upshifted, the OH stretching vibration downshifted and therefore more intensity accumulates in the lower wavenumber transition. The higher sensitivity in the matrix experiment makes it still possible to capture the relative intensity with some degree of confidence. Without deperturbation, the calculation of the matrix shift is somewhat arbitrary, because the jet reference is unclear. It is either anomalously small or anomalously large (Table 1). With deperturbation, once more, one finds a fairly regular Ne-matrix shift of 10 cm–1. The most pronounced swap of the resonance partners with matrix isolation is found for acetone monohydrate (Figure 3, top). The strong-weak pattern of the jet doublet is mirrored to a weak-strong pattern, and the distance between the deperturbed OHb positions (13 cm–1) is significantly more regular than the distance between the dominant peaks in the jet and in the matrix (23 cm–1).

Figure 3 Three cases where the OHb-b2lib resonance is swapped (acetophenone, acetone) or detuned (cyclohexanone) by Ne matrix embedding and a simple resonance treatment regularizes the matrix shift of the OHb mode. The arrows represent the deperturbed band position. W = coupling matrix element required to explain the intensity ratio, listed on the left side of the spectra is 2W, which corresponds to the distance between the two interacting states when perfectly in resonance (e.g., in the second-lowest trace). See text for further explanations.

As a side product and consistency check of the b2lib deperturbation, one may compare the coupling matrix elements W required to explain the intensity pattern between the two environments (Figure 3, which lists 2W as the closest possible approach to the perturbed states). Interestingly, it appears to decrease somewhat from the jet to the matrix for the aliphatic ketones, whereas the coupling quite visibly increases for the aromatic ketone when moving from the jet to the matrix. This could be a real effect, given the modification of the carbonyl group in conjugation with the π system of the aromatic ring. However, one should also consider an alternative explanation that involves the metastable conformer observed for the aromatic ketone (acetophenone, OHb′). While there is no resonance in the jet, the downshift of OHb′ in the matrix compared to the jet wavenumber may light up an OHb′-b2lib′ resonance, whose resonance partner b2lib′ might then coincide with b2lib in the perturbed spectrum. This would lead to an overestimation of b2lib intensity sharing for the global minimum transition and hence an overestimation of its coupling matrix element. However, to bring down the global minimum acetophenone coupling matrix element to the same size as that of the other ketones in the neon matrix (2W ≈ 15 cm–1) one would have to invoke a far too large coupling matrix element for OHb′, the metastable conformer signature. Therefore, the most consistent interpretation still involves an unusually strong resonance for OHb in acetophenone monohydrate and no resonance for OHb′, but any partial contribution of b2lib′ intensity to the b2lib signal might slightly increase the deperturbed matrix shift of 9 cm–1 for OHb and at the same time also slightly increase the listed 10 cm–1 matrix shift of OHb′.

The similar size of regular Ne-matrix shifts and the OHb-b2lib coupling matrix elements (both about 10 cm–1) in ketone hydrates thus creates a perfect situation to achieve two things: A deeper clarification of the resonance situation observed in the gas phase and a regularization of the matrix shifts for the purpose of understanding their magnitude. A future theoretical model of ketone monohydrates in neon matrices might be successful even in scaled harmonic approximation if it is compared to deperturbed experimental data. It will be more challenging to model matrix-isolated monohydrates of ketones including higher-order anharmonic effects responsible for the b2lib resonance.

An example where the b2lib resonance presumably also plays a role but cannot be unambiguously disentangled with the available spectroscopic material is the monohydrate of cyclooctanone. It was shown to coexist in two conformations83 and was part of the HyDRA blind challenge,2 where the focus was on the global minimum conformation, but a metastable conformation was observed in the vibrational jet spectrum as well. In the neon matrix, there is also evidence for both conformations, but now the spectral splitting between the two OH transitions is only 12 instead of 22 cm–1 in the jet. Whereas the 10(2) cm–1 neon matrix shift of the global minimum structure is regular, the shift of the metastable isomer appears too large. Indeed, there is some intensity in the neon spectrum at higher wavenumber, which might be responsible for the downshift of the main signal of the metastable isomer via a b2lib resonance mechanism. While the situation is too congested for a final clarification, the combination of neon matrix and jet spectra may still serve as a resonance alert in such a case, where the jet spectrum alone2 had a regular appearance.

3.1.3 H2O–NH3

The lightest strongly hydrogen-bonded monohydrate in Table 2 is the complex of water with ammonia. Despite its prototype character which makes it suitable for high-level theory treatments and the expected large amplitude zero-point motion of the 5 hydrogen atoms around the pseudodiatomic O–N frame, there is little known about the OH stretching fundamental. In the context of an important experimental dissociation energy study,75 which found preferential energy flow out of the OH stretch excitation into the symmetric (umbrella) NH bend, but not into the OH bend or asymmetric NH bend, a first IR action spectrum was reported. It involved two rotationally unresolved bands, with the higher frequency one remaining somewhat unexplained. Our linear absorption spectra (Figure 4) do not find the higher-frequency band and support the previous assignment of the lower-frequency band to the OH stretching mode. We tentatively resolve a PQR-like structure that is consistent with a ≈10 K rotational temperature and an exponential decay time of >5 ps (corresponding to a Lorentzian FWHM profile of <1 cm–1). For significantly shorter lifetimes, the contour would not gain structure when switching from our standard 2 cm–1 resolution to 1 cm–1. The previously reported band center75 (3485 cm–1) is nicely confirmed within the joint uncertainties of the two experiments, at 3486 cm–1.

Figure 4 Cold gas phase spectra of ammonia and water (w) with He as the carrier gas at different spectral resolutions and stagnation pressures. The suspected R-branch of ammonia monohydrate peaks at 3488 cm–1 and the suspected P-branch at 3484 cm–1. By comparing the mixed spectra with an ammonia-only spectrum (yellow, inverted) in the region of the heterodimer signal, the ammonia monomer signals are identified. The positions of water oligomer bands are labeled with wn (n ≤ 4), in a water-only spectrum (green, inverted, expanded in Ne, and enhanced clustering) to show the low concentration of water trimers and tetramers in the mixed expansions. The three lines around 3200 cm–1 are caused by the water monomer bending overtone.

The corresponding feature in the Ne matrix spectra74 is remarkably downshifted from the gas phase value (Table 2), without any evident nearby resonance partner in the jet and no reported resonance partner in the matrix. The lack of evidence for a resonance is consistent with the experimentally observed energy redistribution,75 which does not involve water bending motion. However, it is still conceivable that a resonance partner is hidden underneath the other absorptions in the relatively dense spectra. This could also be true in neon matrix isolation, where only the peak wavenumber is listed.74 It may be worthwhile to further investigate this particularly simple species in the neon matrix and to explore isotope substitution. Nevertheless, ammonia serves as a first indicator that NH groups coordinated by OH groups might lead to unusually large neon matrix effects, even for relatively modest hydrogen bond shifts.

3.1.4 Other Amine Hydrates

Adding methyl groups to ammonia has a major impact on the strength of the hydrogen bond. Therefore, it is of interest to see how the unusually large matrix shift of the monohydrate of ammonia (vide supra) evolves with increasing methylation. We leave out trimethylamine because the literature is controversial on both sides (jet spectroscopy84 and matrix isolation85) and will be the subject of a separate study, as will the monohydrates of dimethylamine and of tertiary amines.45

The jet spectrum of methylamine monohydrate has a limited signal-to-noise ratio but is surprisingly simple (Figure S3). Besides slight baseline undulations in the region where one might suspect a weak resonance partner (see Section 2), it exhibits a single strong band at 3417 cm–1, indicative of stronger hydrogen bonding than to ammonia. Due to the poor signal-to-noise ratio which may hide a resonance peak, we apply an uncertainty of ±8 cm–1 to the center of the OH intensity. 18O substitution makes sure the band is due to OH and not NH stretching because it is downshifted by the expected amount. The neon matrix spectrum of methylamine monohydrate has been published before43 and shows a significant downshift to 3380 cm–1 as well as significant b2 intensity stealing (quantified by the authors at about 20%). Rather than attempting a deperturbation based on published intensities or an experimental reinvestigation targeting the resonances, which we reserve for a later study, we apply a tentative uncertainty of 10 cm–1 to the literature neon matrix value. Although the effect of the reported b2 intensity would be to increase the matrix shift of the intensity center, there might be a compensating effect from a b2ON state, like in tertiary amines.45 In this way, very much like ammonia, methylamine monohydrate is indicative of an unusually large neon shift of amine hydrates but requires further investigation to narrow down the error bar.

A similar statement can be made for pyridine monohydrate, if one takes the clear-cut single jet transition2 at 3454 cm–1 and compares it to a neon matrix assignment77 of 3419 cm–1, with an estimated uncertainty of ±10 cm–1 due to the congested nature of the matrix spectrum. Though uncertain, it is the third nitrogen acceptor that appears to produce a relatively large neon matrix shift in the water donor OH position, which is worth being reinvestigated.

3.1.5 Replacing Water by Methanol

An important clue for the unusually large water neon matrix shifts indicated by literature values in combination with the current jet results comes from the replacement of water by methanol. This removes some resonance opportunities because the OH bending mode and its overtone (b2) now shift to a lower energy. Indeed, the spectra of methanol–methylamine and methanol–dimethylamine (Figure 5) do not show evidence for resonances and their neon matrix shifts are now down to more normal values (Table 2).

Figure 5 Jet and Ne matrix infrared spectra of methylamine (MA) and dimethylamine (DMA) complexed by a single methanol molecule. The dominance of a single transition demonstrates the absence of anharmonic resonances in the methanol complexes. This facilitates the shift analysis compared to monohydrates of amines, which will be the subject of a future study.

This indicates that the anomalous shift found for methylamine monohydrate may be specific for the light water molecule, potentially due to the large amplitude librational motion of the hydrogen bond, which may be more affected by the neon environment than in the more bulky and heavy methanol case. In the future, this calls for a more systematic investigation of different amine monosolvates.

An important role in such a reanalysis could be played by DOH, which binds to the amines almost in the same way as H2O does. However, monodeuteration removes the b2 resonance option, so that one expects single bands near the HOH features, which might be easily compared between jet and matrix spectra. If the matrix shift of DOH toward amines is close to that of methanol, the special resonance possibilities of HOH might be responsible for the shift anomaly. If it is close to that of HOH, then the accessibility of neon atoms to strong OH···N hydrogen bonds may be a relevant factor. However, due to the expected metastability86 of DOH···N arrangements in comparison to HOD···N coordination (in particular during the long observation times in a matrix) and the accompanying deuteration of nontertiary amines, it is not clear whether monodeuteration will be a viable strategy.

Any deeper analysis of amine monohydrate shifts in neon matrices will need to encompass the more bulky tertiary amines, which have revealed strong resonances in the cold gas phase45 but also offer strong downshifts which will help identify a systematic correlation between hydrogen bond shifts and matrix shifts.

3.1.6 Nitrile Hydrates

The microhydration of nitriles is more subtle because the water interaction is weaker and there is a range of possible coordination geometries involving either the sp-hybridized lone pair of the nitrogen or the π-cloud of the triple bond. The latter allows for some secondary interaction of the water oxygen with an adjacent C–H bond and thus depends on the residue attached to the cyano group. Spectroscopically, they are easy to distinguish because the linear lone pair coordination gives rise to a much larger intensity enhancement and stronger downshift than the bent π-cloud arrangement.71 For acetonitrile, the Ne matrix shift of the observed OH stretching band is in the expected range, and it is plausible to assign the theoretically favored lone pair coordination in both the gas phase and the matrix environment (Figure 6). By spectral analogy, the more downshifted out of two acrylonitrile monohydrate signals can be attributed to the same lone pair coordination. However, now, this coordination type is clearly metastable, as suggested by the addition of neon to the carrier gas and by a microwave structural investigation, which did not observe it.88 The persistent higher wavenumber transition corresponds to a stronger overall interaction and has a negligible neon matrix shift. It must correspond to one of the two nonequivalent bent in-plane π-cloud coordinations. In the jet, it is clear that the less strained global minimum structure (shown in Figure 6) is formed. The experimentally observed, anomalously small neon matrix shift for this isomer (Table 2) might be a consequence of exposure of the negatively polarized nitrogen atom to the matrix. Further investigation extending to benzonitrile89,90 may clarify the assignment.

Figure 6 FTIR jet spectra of the OH stretching region of water–nitrile clusters. While for water–acetonitrile (AceCN, top) only one conformer was observed, two conformers are observed for water–acrylonitrile (AcrCN, bottom) in He as a carrier gas. The admixture of Ne results in depopulation of the metastable conformer. The bars below the spectra indicate relative harmonic wavenumber and intensity predictions (B3LYP-D3(BJ,abc)/def2-TZVP level of theory, using Orca,87 more details in Table S9), wavenumber-scaled by 0.9702 to match the single AceCN band. A more detailed figure is available from the Supporting Information (Figure S2).

3.1.7 Toy Model Capturing O–H···N Matrix Shift Anomalies

It is beyond the scope of this contribution to discuss in detail the possible origins of unusual matrix shifts besides those that can be blamed on tuned anharmonic resonances. Mechanical hydrogen bond compression along the bond direction,14 restrictions of the hydrogen bond libration,80 mechanical matrix cavity size and shape effects,7,13 dielectric effects of the matrix environment,8 accessibility for stabilizing neon atoms in close vicinity of the polarized OH bonds, coupling between the IR photons and matrix phonons, matrix preparation protocols and matrix granularity, or specific host–guest interactions such as weak but cooperative hydrogen bonding26 are among the conceivable origins. Several of these effects will contribute to the 10 cm–1 scatter which is highlighted by the gray ribbon in Figure 1.

We would just like to introduce a simple toy model that might provide a starting point for the future discussion of severe outliers. It assumes that the atoms and groups attached to the central O–H···N unit (which seems to cause particularly strong anomalies) partially block the matrix atoms from laterally perturbing the hydrogen bond. By using dispersion-corrected hybrid-DFT calculations to obtain a plausible structure of the hydrogen-bonded complex (for computational details, see Table S9), one can work out distances from the groups and atoms to the central H atom (di) and sum over the ratio between the van der Waals radius91 of the group or atom ri and this distance di, after raising it to the power of 3. The resulting simple-minded and pairwise additive crowdedness index C = ∑ iN (ri/di)3 can be used as a crude proxy for the accessibility of the hydrogen bond to the matrix environment. By plotting the ratio of the observed Ne matrix shift to the hydrogen bond complexation downshift as a function of C, one obtains a surprisingly monotonous correlation (see Figure 7), within error bars. The less accessible the hydrogen bond circumference is to the matrix atoms (assuming that they do not lead to a major distortion from the predicted gas phase structure), the less its OH stretching mode appears to become affected by the matrix. This might help to rationalize why methanol is a more regular hydrogen bond donor than water in the neon matrix, among other observations. The task ahead is to fill this correlation or variants thereof with many more experimental data points and to check whether the relationship is more coincidental or more systematic and is worth being refined.

Figure 7 Dimensionless ratio of the observed Ne matrix downshift to the hydrogen bond complexation downshift (donor OH wavenumber – complex OH wavenumber, therefore excluding intramolecular hydrogen bonds) plotted against the computed dimensionless crowdedness index C of the complex. More details on ri and di are given in Tables S7 and S8 and all xyz coordinates are shown in Tables S10–S16 in the Supporting Information.

Further experimental investigations should obviously include bulky groups attached to the OH group and the N acceptor because, in the limit of a fully substituent-embedded strong OH–N hydrogen bond, it appears plausible that the neon matrix shift may again vanish due to the inaccessibility of the matrix host to the interaction center. Such new data pairs will thus be particularly valuable in checking any crowdedness index or alternative concepts.

4 Conclusions

Neon matrix shifts are frequently considered to be minor in comparison to uncertainties associated with quantum-chemical calculations, at least in neutral systems. For example, a value typically less than 4 cm–1 has recently been quoted for molecules.92 For diatomics, a systematic study93 revealed small shifts for neon, particularly for diatomic hydrides. In this work, we confirm that many hydrogen-bonded systems stay within about 15 cm–1 of the OH stretching gas phase value after neon matrix embedding. In several cases, this regular behavior is achieved only by systematic deperturbation from resonance partners, which would otherwise distort the picture because they act differently in the matrix and vacuum environment. We also report a few preliminary Ne matrix downshifts for OH–N hydrogen bonds which exceed the typical range and are even larger than the binding energy of a Ne dimer,94 whether or not its zero-point vibrational energy is included. While such large shifts are not unprecedented, they are usually associated with strong acids such as HF, e.g., in combination with reasonably strong bases like ammonia.31 The reported Ne matrix shift of about 109 cm–1 for HF···NH3 actually refers to a thermal gas phase experiment and would likely shrink upon jet cooling. Our tentative evidence for large neon matrix shifts involving water as a very weak acid strongly bound to certain nitrogen-containing acceptors is subject to further investigation. It remains to be seen whether the replacement of water by R–OH always brings back regular matrix shifts, as indicated by a few examples that we have presented.

Highlights of the present extension of the neon shift database include: 1. The first analysis of complexes as elementary as phenol–water (in the matrix, see Figure 2) or water–ammonia (see Figure 4) and water–acetonitrile (see Figure 6, in direct absorption under jet cooling). These can serve as further training systems in future issues of the HyDRA blind challenge.2

2. The consistent removal of matrix shift outliers for ketone hydrates by combining matrix and jet studies with a simple resonance deperturbation model based on intensity centers (see Figure 3). This multi-experimental approach to vibrational resonances promises to become a valuable tool beyond the actual matrix shift topic.

3. The regularization of large matrix shifts through the replacement of microsolvating water by methanol and the associated suppression of stretch-bend Fermi resonance opportunities in strongly hydrogen-bonded amine complexes (see Figure 5). It remains to be seen whether this will also be found for more bulky amine hydrates.45

An extension of the present OH stretching database to NH,95 FH,96 ClH,32,33 and other hydrides97 would be desirable, but for the time being, OH provides the broadest coverage and variation of complexes and our next goal is to systematically extend the OH–N database into the water Fermi resonance regime45 while trying to reduce the uncertainties. This will put the proposed simple toy model to a more rigorous test.

The present data coverage already appears large enough to invite more rigorous attempts to model the matrix shifts, e.g., by polarizable embedding models98 or DFT cluster and bulk studies.99,100 Because one may expect substantial cancellation between repulsive and attractive interaction components in a matrix environment, it appears imperative to have both systematic behavior and outliers included in the growing data set. For this purpose, it would be valuable to build a living database of OH stretching spectra, which are known under both jet cooling and neon matrix isolation conditions. If more such spectral pairs become available in the literature (or if we have missed a literature entry), we would appreciate being contacted. The web site https://qmbench.net/ offers a place where such critically curated living databases can be made accessible in a convenient form.

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpca.4c03468.Experimental details and infrared spectra; additional shift correlations and tables; details on the crowdedness index; and computational details (PDF)

Supplementary Material

jp4c03468_si_001.pdf

Author Contributions

§ M.B. and D.M. contributed equally.

The authors declare no competing financial interest.

Acknowledgments

We thank Jonas Andersen, Taija Fischer, Beppo Hartwig, Nils Lüttschwager, and Eaindra Lwin for valuable discussions on and help with the experiments. MAS acknowledges funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation – 389479699/GRK2455).
==== Refs
References

Del Bene J. E. ; Jordan M. J. T. Vibrational Spectroscopy of the Hydrogen Bond: An Ab Initio Quantum-Chemical Perspective. Int. Rev. Phys. Chem. 1999, 18 , 119–162. 10.1080/014423599230026.
Fischer T. L. ; Bödecker M. ; Schweer S. M. ; Dupont J. ; Lepère V. ; Zehnacker-Rentien A. ; Suhm M. A. ; Schröder B. ; Henkes T. ; Andrada D. M. ; et al. The First HyDRA Challenge for Computational Vibrational Spectroscopy. Phys. Chem. Chem. Phys. 2023, 25 , 22089–22102. 10.1039/D3CP01216F.37610422
Heine N. ; Asmis K. R. Cryogenic Ion Trap Vibrational Spectroscopy of Hydrogen-Bonded Clusters Relevant to Atmospheric Chemistry. Int. Rev. Phys. Chem. 2015, 34 , 1–34. 10.1080/0144235X.2014.979659.
Ruoff R. S. ; Klots T. D. ; Emilsson T. ; Gutowsky H. S. Relaxation of Conformers and Isomers in Seeded Supersonic Jets of Inert Gases. J. Chem. Phys. 1990, 93 , 3142–3150. 10.1063/1.458848.
Kollipost F. ; Andersen J. ; Mahler D. W. ; Heimdal J. ; Heger M. ; Suhm M. A. ; Wugt Larsen R. The Effect of Hydrogen Bonding on Torsional Dynamics: a Combined Far-Infrared Jet and Matrix Isolation Study of Methanol Dimer. J. Chem. Phys. 2014, 141 , 174314 10.1063/1.4900922.25381521
Heger M. ; Andersen J. ; Suhm M. A. ; Wugt Larsen R. The Donor OH Stretching-Libration Dynamics of Hydrogen-Bonded Methanol Dimers in Cryogenic Matrices. Phys. Chem. Chem. Phys. 2016, 18 , 3739–3745. 10.1039/C5CP07387A.26763101
Knözinger E. ; Schuller W. ; Langel W. Structure and Dynamics in Pure and Doped Rare-Gas Matrices. Faraday Discuss. Chem. Soc. 1988, 86 , 285–293. 10.1039/DC9888600285.
Jordan M. J. T. ; Del Bene J. E. Unraveling Environmental Effects on Hydrogen-Bonded Complexes: Matrix Effects on the Structures and Proton-Stretching Frequencies of Hydrogen–Halide Complexes with Ammonia and Trimethylamine. J. Am. Chem. Soc. 2000, 122 , 2101–2115. 10.1021/ja993981s.
Bochenkova A. V. ; Suhm M. A. ; Granovsky A. A. ; Nemukhin A. V. Hybrid Diatomics-In-Molecules-Based Quantum Mechanical/Molecular Mechanical Approach Applied to the Modeling of Structures and Spectra of Mixed Molecular Clusters Arn(HCl)m and Arn(HF)m. J. Chem. Phys. 2004, 120 , 3732–3743. 10.1063/1.1642596.15268536
Bochenkova A. V. ; Bochenkov V. E. ; Khriachtchev L. HArF in Solid Argon Revisited: Transition from Unstable to Stable Configuration. J. Phys. Chem. A 2009, 113 , 7654–7659. 10.1021/jp810457h.19243121
Moudens A. ; Georges R. ; Goubet M. ; Makarewicz J. ; Lokshtanov S. E. ; Vigasin A. A. Direct Absorption Spectroscopy of Water Clusters Formed in a Continuous Slit Nozzle Expansion. J. Chem. Phys. 2009, 131 , 204312 10.1063/1.3264576.19947685
Simon A. ; Iftner C. ; Mascetti J. ; Spiegelman F. Water Clusters in an Argon Matrix: Infrared Spectra from Molecular Dynamics Simulations with a Self-Consistent Charge Density Functional-Based Tight Binding/Force-Field Potential. J. Phys. Chem. A 2015, 119 , 2449–2467. 10.1021/jp508533k.25650885
Ito F. Modeling and Spectral Simulation of Formic Acid Dimer in Ar Matrix Using ONIOM Calculations. Comput. Theor. Chem. 2019, 1161 , 18–25. 10.1016/j.comptc.2019.06.005.
Bader F. ; Tremblay J. C. ; Paulus B. A Pair Potential Modeling Study of F3– in Neon Matrices. Phys. Chem. Chem. Phys. 2021, 23 , 886–899. 10.1039/D0CP05031H.33346758
Coussan S. ; Alikhani M. E. ; Perchard J. P. ; Zheng W. Q. Infrared-Induced Isomerization of Ethanol Dimers Trapped in Argon and Nitrogen Matrices: Monochromatic Irradiation Experiments and DFT Calculations. J. Phys. Chem. A 2000, 104 , 5475–5483. 10.1021/jp0001283.
Emmeluth C. ; Dyczmons V. ; Kinzel T. ; Botschwina P. ; Suhm M. A. ; Yáñez M. Combined Jet Relaxation and Quantum Chemical Study of the Pairing Preferences of Ethanol. Phys. Chem. Chem. Phys. 2005, 7 , 991–997. 10.1039/B417870J.19791391
Asselin P. ; Goubet M. ; Lewerenz M. ; Soulard P. ; Perchard J. P. Rovibrational and Dynamical Properties of the Hydrogen Bonded Complex (CH2)2S-HF: A Combined Free Jet, Cell, and Neon Matrix-Fourier Transform Infrared Study. J. Chem. Phys. 2004, 121 , 5241–5252. 10.1063/1.1780167.15352817
Cirtog M. ; Asselin P. ; Soulard P. ; Tremblay B. ; Madebène B. ; Alikhani M. E. ; Georges R. ; Moudens A. ; Goubet M. ; Huet T. R. ; et al. The (CH2)2O-H2O Hydrogen Bonded Complex. Ab Initio calculations and Fourier Transform Infrared Spectroscopy from Neon Matrix and a New Supersonic Jet Experiment Coupled to the Infrared AILES Beamline of Synchrotron SOLEIL. J. Phys. Chem. A 2011, 115 , 2523–2532. 10.1021/jp111507z.21381647
Ceponkus J. ; Engdahl A. ; Uvdal P. ; Nelander B. Structure and Dynamics of Small Water Clusters, Trapped in Inert Matrices. Chem. Phys. Lett. 2013, 581 , 1–9. 10.1016/j.cplett.2013.06.046.
Dargent D. ; Madebène B. ; Soulard P. ; Tremblay B. ; Zins E. L. ; Alikhani M. E. ; Asselin P. Conformational Landscape of the 1/1 Diacetyl/Water Complex Investigated by Infrared Spectroscopy and Ab Initio Calculations. J. Phys. Chem. A 2017, 121 , 88–97. 10.1021/acs.jpca.6b10492.27959532
Heger M. ; Suhm M. A. ; Mata R. A. Communication: Towards The Binding Energy and Vibrational Red shift of the Simplest Organic Hydrogen Bond: Harmonic Constraints for Methanol Dimer. J. Chem. Phys. 2014, 141 , 101105 10.1063/1.4895728.25217897
Nelander B. The Intramolecular Fundamentals of the Water Dimer. J. Chem. Phys. 1988, 88 , 5254–5256. 10.1063/1.454584.
Wuelfert S. ; Herren D. ; Leutwyler S. Reply to “The Intramolecular Fundamentals of the Water Dimer. J. Chem. Phys. 1988, 88 , 5256–5257. 10.1063/1.454585.
Sakpal S. ; Chakrabarty S. ; Reddy K. D. ; Deshmukh S. H. ; Biswas R. ; Bagchi S. ; Ghosh A. Perturbation of Fermi Resonance on Hydrogen-Bonded CO: 2D IR Studies of Small Ester Probes. J. Phys. Chem. B 2024, 128 , 4440–4447. 10.1021/acs.jpcb.3c06698.38686937
Oswald S. ; Suhm M. A. Experimental Reference Data for Hexafluorinated Propanol by Exploring an Unusual Intermolecular Torsional Balance. Angew. Chem., Int. Ed. 2017, 56 , 12672–12676. 10.1002/anie.201705301.
Oswald S. ; Wallrabe M. ; Suhm M. A. Cooperativity in Alcohol-Nitrogen Complexes: Understanding Cryomatrices through Slit Jet Expansions. J. Phys. Chem. A 2017, 121 , 3411–3422. 10.1021/acs.jpca.7b01265.28443670
Carrascosa E. ; Pellegrinelli R. P. ; Rizzo T. R. ; Muyskens M. A. Cryogenic Infrared Action Spectroscopy Fingerprints the Hydrogen Bonding Network in Gas-Phase Coumarin Cations. J. Phys. Chem. A 2020, 124 , 9942–9950. 10.1021/acs.jpca.0c06430.33201700
Xu Y. ; van Wijngaarden J. ; Jäger W. Microwave Spectroscopy of Ternary and Quaternary van der Waals Clusters. Int. Rev. Phys. Chem. 2005, 24 , 301–338. 10.1080/01442350500252039.
Lu T. ; Obenchain D. A. ; Zhang J. ; Grabow J.-U. ; Feng G. Van der Waals Interactions of the Disulfide Bond Revealed: A Microwave Spectroscopic Study of the Diethyl Disulfide-Argon Complex. J. Chem. Phys. 2021, 154 , 124306 10.1063/5.0043615.33810705
Rock C. A. ; Tschumper G. S. Insight into the Binding of Argon to Cyclic Water Clusters from Symmetry-Adapted Perturbation Theory. Int. J. Mol. Sci. 2023, 24 , 17480 10.3390/ijms242417480.38139311
Andrews L. ; Wang X. ; Mielke Z. Infrared Spectrum of the H3N–HCl Complex in Solid Ne, Ne/Ar, Ar, and Kr. Matrix Effects on a Strong Hydrogen-Bonded Complex. J. Phys. Chem. A 2001, 105 , 6054–6064. 10.1021/jp010323x.
Weiss N. M. ; Waller A. W. ; Phillips J. A. Infrared Spectrum of CH3CN–HCl in Solid Neon, and Modeling Matrix Effects in CH3CN–HCl and H3N–HCl. J. Mol. Struct. 2016, 1105 , 341–349. 10.1016/j.molstruc.2015.10.014.
Soares C. ; Ley A. R. ; Zehner B. C. ; Treacy P. W. ; Phillips J. A. Matrix Effects on Hydrogen Bonding and Proton Transfer in Fluoropyridine - HCl Complexes. Phys. Chem. Chem. Phys. 2022, 24 , 2371–2386. 10.1039/D1CP04110J.35019906
Mata R. A. ; Suhm M. A. Benchmarking Quantum Chemical Methods: Are We Heading in the Right Direction?. Angew. Chem., Int. Ed. 2017, 56 , 11011–11018. 10.1002/anie.201611308.
Keutsch F. N. ; Saykally R. J. Water Clusters: Untangling the Mysteries of the Liquid, one Molecule at a Time. Proc. Natl. Acad. Sci. U. S. A. 2001, 98 , 10533–10540. 10.1073/pnas.191266498.11535820
Mihrin D. ; Andersen J. ; Jakobsen P. W. ; Wugt Larsen R. Highly Localized H2O Librational Motion as a Far-Infrared Spectroscopic Probe for Microsolvation of Organic Molecules. Phys. Chem. Chem. Phys. 2019, 21 , 1717–1723. 10.1039/C8CP05985C.30623967
Mihrin D. ; Voute A. ; Jakobsen P. W. ; Feilberg K. L. ; Wugt Larsen R. The Effect of Alkylation on the Micro-Solvation of Ethers Revealed by Highly Localized Water Librational Motion. J. Chem. Phys. 2022, 156 , 084305 10.1063/5.0081161.35232174
Lovejoy C. M. ; Nesbitt D. J. Slit Pulsed Valve for Generation of Long-Path-Length Supersonic Expansions. Rev. Sci. Instrum. 1987, 58 , 807–811. 10.1063/1.1139636.
Liu K. ; Fellers R. S. ; Viant M. R. ; McLaughlin R. P. ; Brown M. G. ; Saykally R. J. A Long Path Length Pulsed Slit Valve Appropriate for High Temperature Operation: Infrared Spectroscopy of Jet-Cooled Large Water Clusters and Nucleotide Bases. Rev. Sci. Instrum. 1996, 67 , 410–416. 10.1063/1.1146605.
Gottschalk H. C. ; Fischer T. L. ; Meyer V. ; Hildebrandt R. ; Schmitt U. ; Suhm M. A. A Sustainable Slit Jet FTIR Spectrometer for Hydrate Complexes and Beyond. Instruments 2021, 5 , 12 10.3390/instruments5010012.
Andersen J. ; Heimdal J. ; Wugt Larsen R. The Influence of Large-Amplitude Librational Motion on the Hydrogen Bond Energy for Alcohol-Water Complexes. Phys. Chem. Chem. Phys. 2015, 17 , 23761–23769. 10.1039/C5CP04321B.26304774
Andersen J. ; Heimdal J. ; Wugt Larsen R. Spectroscopic Identification of Ethanol-Water Conformers by Large-Amplitude Hydrogen Bond Librational Modes. J. Chem. Phys. 2015, 143 , 224315 10.1063/1.4937482.26671383
Soulard P. ; Tremblay B. Vibrational Study of Methylamine Dimer and Hydrated Methylamine Complexes in Solid Neon Supported by Ab Initio Calculations. J. Mol. Struct. 2021, 1236 , 130308 10.1016/j.molstruc.2021.130308.
Fischer T. L. ; Wagner T. ; Gottschalk H. C. ; Nejad A. ; Suhm M. A. A Rather Universal Vibrational Resonance in 1:1 Hydrates of Carbonyl Compounds. J. Phys. Chem. Lett. 2021, 12 , 138–144. 10.1021/acs.jpclett.0c03197.33315407
Lwin E. ; Fischer T. L. ; Suhm M. A. Microhydration of Tertiary Amines: Robust Resonances in Red-Shifted Water. J. Phys. Chem. Lett. 2023, 14 , 10194–10199. 10.1021/acs.jpclett.3c02517.37930195
Zabuga A. V. ; Kamrath M. Z. ; Rizzo T. R. Franck-Condon-like Progressions in Infrared Spectra of Biological Molecules. J. Phys. Chem. A 2015, 119 , 10494–10501. 10.1021/acs.jpca.5b08801.26436667
Nejad A. ; Pérez Mellor A. F. ; Lange M. ; Alata I. ; Zehnacker A. ; Suhm M. A. Subtle Hydrogen Bond Preference and Dual Franck-Condon Activity - The Interesting Pairing of 2-Naphthol with Anisole. Phys. Chem. Chem. Phys. 2023, 25 , 10427–10439. 10.1039/D3CP00160A.37000584
Liu Y. ; Li J. ; Felker P. M. ; Bačić Z. HCl-H2O Dimer: An Accurate Full-Dimensional Potential Energy Surface and Fully Coupled Quantum Calculations of Intra- and Intermolecular Vibrational States and Frequency Shifts. Phys. Chem. Chem. Phys. 2021, 23 , 7101–7114. 10.1039/D1CP00865J.33876076
Wang X.-G. ; Carrington T. Computing Excited OH Stretch States of Water Dimer in 12D Using Contracted Intermolecular and Intramolecular Basis Functions. J. Chem. Phys. 2023, 158 , 084107 10.1063/5.0139586.36859104
Lüttschwager N. O. B. The Strength of the OH-Bend/OH-Stretch Fermi Resonance in Small Water Clusters. Phys. Chem. Chem. Phys. 2024, 26 , 10120–10135. 10.1039/D3CP06255D.38487881
Lüttschwager N. NoisySignalIntegration.jl: A Julia Package for Uncertainty Evaluation of Numeric Integrals. JOSS 2021, 6 , 3526 10.21105/joss.03526.
Huisken F. ; Kaloudis M. ; Kulcke A. Infrared Spectroscopy of Small Size-Selected Water Clusters. J. Chem. Phys. 1996, 104 , 17–25. 10.1063/1.470871.
Bouteiller Y. ; Perchard J. P. The Vibrational Spectrum of (H2O)2: Comparison between Anharmonic Ab Initio Calculations and Neon Matrix Infrared Data Between 9000 and 90 cm–1. Chem. Phys. 2004, 305 , 1–12. 10.1016/j.chemphys.2004.06.028.
Paul J. B. ; Collier C. P. ; Saykally R. J. ; Scherer J. J. ; O’Keefe A. Direct Measurement of Water Cluster Concentrations by Infrared Cavity Ringdown Laser Absorption Spectroscopy. J. Phys. Chem. A 1997, 101 , 5211–5214. 10.1021/jp971216z.
Otto K. E. ; Xue Z. ; Zielke P. ; Suhm M. A. The Raman Spectrum of Isolated Water Clusters. Phys. Chem. Chem. Phys. 2014, 16 , 9849–9858. 10.1039/c3cp54272f.24398903
León I. ; Montero R. ; Castaño F. ; Longarte A. ; Fernández J. A. Mass-Resolved Infrared Spectroscopy of Complexes without Chromophore by Nonresonant Femtosecond Ionization Detection. J. Phys. Chem. A 2012, 116 , 6798–6803. 10.1021/jp303937h.22639969
Provencal R. A. ; Paul J. B. ; Roth K. ; Chapo C. ; Casaes R. N. ; Saykally R. J. ; Tschumper G. S. ; Schaefer H. F. Infrared Cavity Ringdown Spectroscopy of Methanol Clusters: Single Donor Hydrogen Bonding. J. Chem. Phys. 1999, 110 , 4258–4267. 10.1063/1.478309.
Kollipost F. ; Papendorf K. ; Lee Y.-F. ; Lee Y.-P. ; Suhm M. A. Alcohol Dimers - How Much Diagonal OH Anharmonicity?. Phys. Chem. Chem. Phys. 2014, 16 , 15948–15956. 10.1039/C4CP01418A.24963724
Larsen R. W. ; Zielke P. ; Suhm M. A. Hydrogen-Bonded OH Stretching Modes of Methanol Clusters: A Combined IR and Raman Isotopomer Study. J. Chem. Phys. 2007, 126 , 194307 10.1063/1.2732745.17523804
Nedić M. ; Wassermann T. N. ; Larsen R. W. ; Suhm M. A. A Combined Raman- and Infrared Jet Study of Mixed Methanol-Water and Ethanol-Water Clusters. Phys. Chem. Chem. Phys. 2011, 13 , 14050–14063. 10.1039/c1cp20182d.21491035
Wassermann T. N. ; Suhm M. A. Ethanol Monomers and Dimers Revisited: A Raman Study of Conformational Preferences and Argon Nanocoating Effects. J. Phys. Chem. A 2010, 114 , 8223–8233. 10.1021/jp104861q.20701332
Fischer T. L. Of Resonances and Radicals – Hydrate Studies for Benchmarking, 10.53846/goediss-10364, PhD thesis. 2024 (accessed Apr 15, 2024).
Ebata T. ; Watanabe T. ; Mikami N. Evidence for the Cyclic Form of Phenol Trimer: Vibrational Spectroscopy of the OH Stretching Vibrations of Jet-Cooled Phenol Dimer and Trimer. J. Phys. Chem. 1995, 99 , 5761–5764. 10.1021/j100016a006.
Plokhotnichenko A. M. ; Radchenko E. D. ; Blagoi Y. P. ; Karachevtsev V. A. Dimers of Phenol in Argon and Neon Matrices. Low Temp. Phys. 2001, 27 , 666–675. 10.1063/1.1399207.
Emmeluth C. ; Dyczmons V. ; Suhm M. A. Tuning the Hydrogen Bond Donor/Acceptor Isomerism in Jet-Cooled Mixed Dimers of Aliphatic Alcohols. J. Phys. Chem. A 2006, 110 , 2906–2915. 10.1021/jp0540407.16509612
Doi A. ; Mikami N. Dynamics of Hydrogen-Bonded OH Stretches as Revealed by Single-Mode Infrared-Ultraviolet Laser Double Resonance Spectroscopy on Supersonically Cooled Clusters of Phenol. J. Chem. Phys. 2008, 129 , 154308 10.1063/1.2988494.19045194
Kwasniewski D. ; Butler M. ; Reisler H. Vibrational Predissociation of the Phenol-Water Dimer: A View from the Water. Phys. Chem. Chem. Phys. 2019, 21 , 13968–13976. 10.1039/C8CP06581K.30511053
Hartland G. V. ; Henson B. F. ; Venturo V. A. ; Felker P. M. Ionization-Loss Stimulated Raman Spectroscopy of Jet-Cooled Hydrogen-Bonded Complexes Containing Phenols. J. Phys. Chem. 1992, 96 , 1164–1173. 10.1021/j100182a028.
Tanabe S. ; Ebata T. ; Fujii M. ; Mikami N. OH Stretching Vibrations of Phenol–(H2O)n (n = 1–3) Complexes Observed by IR-UV Double-Resonance Spectroscopy. Chem. Phys. Lett. 1993, 215 , 347–352. 10.1016/0009-2614(93)85726-5.
Watanabe T. ; Ebata T. ; Tanabe S. ; Mikami N. Size-Selected Vibrational Spectra of Phenol-(H2O)n (n = 1–4) Clusters Observed by IR–UV Double Resonance and Stimulated Raman-UV Double Resonance Spectroscopies. J. Chem. Phys. 1996, 105 , 408–419. 10.1063/1.471917.
Soulard P. ; Tremblay B. Vibrational Study of Acrylonitrile Dimer and Acrylonitrile-Water Hydrogen-Bonded Complexes in Solid Neon Supported by Ab Initio Calculations. J. Mol. Struct. 2022, 1266 , 133503 10.1016/j.molstruc.2022.133503.
Asselin P. ; Madebène B. ; Soulard P. ; Georges R. ; Goubet M. ; Huet T. R. ; Pirali O. ; Zehnacker-Rentien A. Competition Between Inter- and Intra-Molecular Hydrogen Bonding: An Infrared Spectroscopic Study of Jet-Cooled Amino-Ethanol and its Dimer. J. Chem. Phys. 2016, 145 , 224313 10.1063/1.4972016.27984904
Liu Y. ; Rice C. A. ; Suhm M. A. Torsional Isomers in Methylated Aminoethanols - A Jet-FTIR Study. Can. J. Chem. 2004, 82 , 1006–1012. 10.1139/v04-046.
Engdahl A. ; Nelander B. The Intramolecular Vibrations of the Ammonia Water Complex. A Matrix Isolation Study. J. Chem. Phys. 1989, 91 , 6604–6612. 10.1063/1.457379.
Mollner A. K. ; Casterline B. E. ; Ch’ng L. C. ; Reisler H. Imaging the State-Specific Vibrational Predissociation of the Ammonia-Water Hydrogen-Bonded Dimer. J. Phys. Chem. A 2009, 113 , 10174–10183. 10.1021/jp904566w.19715290
Millen D. J. ; Mines G. W. Hydrogen Bonding in the Gas Phase. Part 5.—Infrared Spectroscopic Investigation of O–H ··· N Complexes Formed by Water: Ammonia Monohydrate and Amine and Pyridine Monohydrates. J. Chem. Soc., Faraday Trans. 2 1977, 73 (3 ), 369–377. 10.1039/F29777300369.
Esteves-López N. ; Coussan S. UV Photochemistry of Pyridine-Water and Pyridine-Ammonia Complexes Trapped in Cryogenic Matrices. J. Mol. Struct. 2018, 1172 , 65–73. 10.1016/j.molstruc.2017.11.037.
Jiang S. ; Kong X. ; Wang C. ; Zang X. ; Su M. ; Zheng H. ; Zhang B. ; Li G. ; Xie H. ; Yang X. ; et al. Infrared Spectroscopy of Hydrogen-Bonding Interactions in Neutral Dimethylamine-Methanol Complexes. J. Phys. Chem. A 2019, 123 , 10109–10115. 10.1021/acs.jpca.9b08630.31660735
Ali O. Y. ; Jewer E. ; Fridgen T. D. Infrared Spectroscopic Characterization of Hydrogen-Bonded Propylene Oxide–Ethanol and Propylene Oxide–2-Fluoroethanol Complexes Isolated in Solid Neon Matrices. Can. J. Chem. 2013, 91 , 1292–1302. 10.1139/cjc-2013-0355.
Choi M. Y. ; Douberly G. E. ; Falconer T. M. ; Lewis W. K. ; Lindsay C. M. ; Merritt J. M. ; Stiles P. L. ; Miller R. E. Infrared Spectroscopy of Helium Nanodroplets: Novel Methods for Physics and Chemistry. Int. Rev. Phys. Chem. 2006, 25 , 15–75. 10.1080/01442350600625092.
Andersen J. ; Heimdal J. ; Mahler D. W. ; Nelander B. ; Larsen R. W. Communication: THz Absorption Spectrum of the CO2–H2O Complex: Observation and Assignment of Intermolecular van der Waals Vibrations. J. Chem. Phys. 2014, 140 , 091103 10.1063/1.4867901.24606346
Dinu D. F. ; Podewitz M. ; Grothe H. ; Liedl K. R. ; Loerting T. Toward Elimination of Discrepancies between Theory and Experiment: Anharmonic Rotational-Vibrational Spectrum of Water in Solid Noble Gas Matrices. J. Phys. Chem. A 2019, 123 , 8234–8242. 10.1021/acs.jpca.9b07221.31433184
Burevschi E. ; Peña I. ; Sanz M. E. Geminal Diol Formation from the Interaction of a Ketone with Water in the Gas Phase: Structure and Reactivity of Cyclooctanone- (H2O)1,2 Clusters. J. Phys. Chem. Lett. 2021, 12 , 12419–12425. 10.1021/acs.jpclett.1c03493.34939809
Jiang S. ; Su M. ; Yang S. ; Wang C. ; Huang Q.-R. ; Li G. ; Xie H. ; Yang J. ; Wu G. ; Zhang W. ; et al. Vibrational Signature of Dynamic Coupling of a Strong Hydrogen Bond. J. Phys. Chem. Lett. 2021, 12 , 2259–2265. 10.1021/acs.jpclett.1c00168.33636082
Rozenberg M. ; Loewenschuss A. ; Nielsen C. J. H-Bonded Clusters in the Trimethylamine/Water System: A Matrix Isolation and Computational Study. J. Phys. Chem. A 2012, 116 , 4089–4096. 10.1021/jp3020035.22489577
Némethy G. ; Scheraga H. A. Structure of Water and Hydrophobic Bonding in Proteins. IV. The Thermodynamic Properties of Liquid Deuterium Oxide. J. Chem. Phys. 1964, 41 , 680–689. 10.1063/1.1725946.
Neese F. Software Update: The ORCA Program System—Version 5.0. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2022, 12 , e1606 10.1002/wcms.1606.
Calabrese C. ; Vigorito A. ; Maris A. ; Mariotti S. ; Fathi P. ; Geppert W. D. ; Melandri S. Millimeter Wave Spectrum of the Weakly Bound Complex CH2CHCN–H2O: Structure, Dynamics, and Implications for Astronomical Search. J. Phys. Chem. A 2015, 119 , 11674–11682. 10.1021/acs.jpca.5b08426.26524193
Ishikawa S. ; Ebata T. ; Mikami N. Structures and the Vibrational Relaxations of Size-Selected Benzonitrile–(H2O)n=1–3 and–(CH3OH)n=1–3 Clusters Studied by Fluorescence Detected Raman and Infrared Spectroscopies. J. Chem. Phys. 1999, 110 , 9504–9515. 10.1063/1.478915.
Khatri J. ; Roy T. K. ; Chatterjee K. ; Schwaab G. ; Havenith M. Vibrational Spectroscopy of Benzonitrile-(Water)1–2 Clusters in Helium Droplets. J. Phys. Chem. A 2021, 125 , 6954–6963. 10.1021/acs.jpca.1c04553.34355893
Bondi A. Van der Waals Volumes and Radii. J. Phys. Chem. 1964, 68 , 441–451. 10.1021/j100785a001.
Dinu D. F. ; Podewitz M. ; Grothe H. ; Loerting T. ; Liedl K. R. On the Synergy of Matrix-Isolation Infrared Spectroscopy and Vibrational Configuration Interaction Computations. Theor. Chem. Acc. 2020, 139 , 174 10.1007/s00214-020-02682-0.33192169
Jacox M. E. Comparison of the Ground State Vibrational Fundamentals of Diatomic Molecules in the Gas Phase and in Inert Solid Matrices. J. Mol. Spectrosc. 1985, 113 , 286–301. 10.1016/0022-2852(85)90268-1.
Wüest A. ; Merkt F. Determination of the Interaction Potential of the Ground Electronic State of Ne2 by High-Resolution Vacuum Ultraviolet Laser Spectroscopy. J. Chem. Phys. 2003, 118 , 8807–8812. 10.1063/1.1566944.
Oswald S. ; Suhm M. A. ; Coussan S. Incremental NH Stretching Downshift through Stepwise Nitrogen Complexation of Pyrrole: A Combined Jet Expansion and Matrix Isolation Study. Phys. Chem. Chem. Phys. 2019, 21 , 1277–1284. 10.1039/C8CP07053A.30569059
Goubet M. ; Asselin P. ; Soulard P. ; Perchard J. P. Neon Matrix Induced Perturbations of Weak Hydrogen-Bonded Complexes. Phys. Chem. Chem. Phys. 2003, 5 , 5365 10.1039/b311252g.
Graneri M. H. V. ; Spagnoli D. ; Wild D. A. ; McKinley A. J. Twin Peaks: Matrix Isolation Studies of H2S· Amine Complexes Shedding Light on Fundamental S-H··· N Bonding. J. Chem. Phys. 2024, 160 , 124312 10.1063/5.0191308.38545947
van den Heuvel W. ; Reinholdt P. ; Kongsted J. Embedding Beyond Electrostatics: The Extended Polarizable Density Embedding Model. J. Phys. Chem. B 2023, 127 , 3248–3256. 10.1021/acs.jpcb.2c08721.37002869
Mahjoubi K. ; Benoit D. M. ; Jaidane N.-E. ; Al-Mogren M. M. ; Hochlaf M. Understanding of Matrix Embedding: A Theoretical Spectroscopic Study of CO Interacting with Ar Clusters, Surfaces and Matrices. Phys. Chem. Chem. Phys. 2015, 17 , 17159–17168. 10.1039/C5CP01672J.26067278
Bader F. ; Lindic T. ; Paulus B. A Validation of Cluster Modeling in the Description of Matrix Isolation Spectroscopy. J. Comput. Chem. 2020, 41 , 751–758. 10.1002/jcc.26123.31804712
