Self-Association and Microhydration of Phenol: Identification of Large-Amplitude Hydrogen Bond Librational Modes

The self-association mechanisms of phenol have represented long-standing challenges to quantum chemical methodologies owing to the competition between strongly directional intermolecular hydrogen bonding, weaker non-directional London dispersion forces and C–H⋯π interactions between the aromatic rings. The present work explores these subtle self-association mechanisms of relevance for biological molecular recognition processes via spectroscopic observations of large-amplitude hydrogen bond librational modes of phenol cluster molecules embedded in inert neon “quantum” matrices complemented by domain-based local pair natural orbital-coupled cluster DLPNO-CCSD(T) theory. The spectral signatures confirm a primarily intermolecular O-H⋯H hydrogen-bonded structure of the phenol dimer strengthened further by cooperative contributions from inter-ring London dispersion forces as supported by DLPNO-based local energy decomposition (LED) predictions. In the same way, the hydrogen bond librational bands observed for the trimeric cluster molecule confirm a pseudo-C3 symmetric cyclic cooperative hydrogen-bonded barrel-like potential energy minimum structure. This structure is vastly different from the sterically favored “chair” conformations observed for aliphatic alcohol cluster molecules of the same size owing to the additional stabilizing London dispersion forces and C–H⋯π interactions between the aromatic rings. The hydrogen bond librational transition observed for the phenol monohydrate finally confirms that phenol acts as a hydrogen bond donor to water in contrast to the hydrogen bond acceptor role observed for aliphatic alcohols.


Introduction
Hydrogen bonding stands as a cornerstone in both the structure and function of supramolecular chemical interactions, governing the majority of biological recognition mechanisms and the properties of solvents, and guiding chemical reactivity.In intricate environments, highly directional hydrogen bonds represent just one among several other classes of non-covalent intermolecular interactions contributing to molecular association mechanisms.Non-directional attractive London dispersion forces constitute an integral component of every non-covalent interaction, contributing a substantial portion of the total association energy within any given molecular aggregate growing with the size of the molecular system.London dispersion forces become particularly important when the primary bond is weaker than typical conventional hydrogen bonds, such as those between O-H• • •N and O-H• • •O motifs [1].In addition, secondary interactions such as non-conventional O-H• • •π hydrogen bonds and even weaker C-H• • •π contacts often compete effectively with conventional hydrogen bonds in aromatic molecular systems [2,3].The complex interplay between London dispersion, electrostatic forces, and weak donor-acceptor interactions represents a significant theoretical challenge for predictive ab initio methods targeting realistic supramolecular assemblies.Currently, only high-level post-Hartree-Fock ab initio quantum chemistry approaches such as coupled cluster theory CCSD(T) [4][5][6] consistently provide accurate predictions of non-covalent interactions, particularly in the presence of extensive van der Waals forces.However, the computational cost of these advanced ab initio methods scales rapidly, rendering them impractical for larger molecular systems, whereas lower-tier methods such as density functional theory (DFT), despite reliance on empirical corrections, often fail to achieve the necessary chemical accuracy [7].The prototypical homocluster molecules of phenol (PhOH) have been investigated extensively and have earned a place in benchmark databases for biomolecular systems [8], although the computational costs for theoretical CCSD(T) predictions become prohibitive already at the (PhOH) 2 level due to the size of the aromatic fragments [9,10].The introduction of these aromatic structures facilitates the formation of intermolecular contacts involving the π-electron clouds.The self-association mechanisms of phenol thus involve an interplay of both conventional primary O-H• • •O hydrogen bonds, substantial dispersive contributions and secondary C-H• • •π contacts constituting a perfect playground for the validation of lower-tier theoretical methodologies with accurate spectroscopic findings.
The currently available experimental spectroscopic data primarily consist of structural and vibrational observations for phenol dimer (PhOH) 2 and phenol trimer (PhOH) 3 as well as various phenol microhydrate cluster molecules PhOH•(H 2 O) 1-3 [11][12][13][14][15][16][17][18][19][20].While both microwave spectroscopy guided substitution structures and mid-IR transitions associated with strongly active intramolecular vibrational OH-stretching modes have been reported in the literature for these homo-and microhydrate cluster molecules, there are currently no reported observations of the low-energy transitions associated with the large-amplitude and highly anharmonic OH librational (hindered torsional) modes linked directly to these (cooperative) hydrogen bond networks.This class of large-amplitude hydrogen bond librational motion has been shown to be an accurate probe of the strength and directionality of the intermolecular hydrogen bonds [21][22][23].For aliphatic alcohols, the alkyl groups in cyclic cluster molecules have previously been shown experimentally to avoid each other due to steric hindrance resulting in chair-like potential energy minima for alcohol trimers and "up-down-up-down" (relative to the hydrogen bond plane) minima for alcohol tetramers [24].Furthermore, the spectroscopic detection of these large-amplitude hydrogen bond librational modes reveals the most dominating contributions to the change in vibrational zero-point energy upon complexation (∆ZPE) [25].These observations help to translate high-level ab initio quantum chemical predictions of electronic equilibrium dissociation energies (D e ) into accurate semi-empirical ground-state dissociation energies (D 0 ), which is notoriously challenging from first principles quantum chemistry alone due to the highly anharmonic nature of this class of large-amplitude vibrational motion [26][27][28][29].
In the present work, we provide for the first time spectroscopic assignments for the experimentally less-accessible vibrational transitions associated with large-amplitude hydrogen bond librational motions for phenol dimer, phenol trimer and isotopic variants of the phenol monohydrate (PhOH•H 2 O/PhOH•D 2 O) embedded in inert cryogenic "quantum" neon matrices at 4 K.These experimental findings are complemented by systematic theoretical conformation searches and harmonic force field predictions employing both dispersion-corrected DFT and ab initio methodologies of the most stable potential energy minima.In addition, detailed analyses of the total interaction energies in these phenol cluster molecules are carried out by means of a local energy decomposition (LED) scheme to account for the competition between the different classes of intermolecular non-covalent forces.

Phenol Cluster Molecules
Figure 1 shows the mid-infrared absorption spectra of doped neon matrices embedded with phenol at two different mixing ratios (by thermostating the phenol sublimation vessel at 0 • C and at −5 • C, respectively).In the former, most concentrated neon matrix experiment, a so-called pre-annealing spectrum (blue trace) was recorded immediately after the deposition of the sample and a post-annealing spectrum (red trace) was collected after the annealing procedure.The primary distinct bands observed in the spectra above 3400 cm −1 are associated with the strongly IR-active intramolecular OH stretching modes of phenol monomer and the intensity-enhanced donor OH-stretching modes of phenol cluster molecules abundant in sub-% concentrations relative to the monomer.The bands observed at 3655 cm −1 , 3515 cm −1 and 3431 cm −1 belong to the free OH-stretching mode of phenol monomer [11,13], the hydrogenbonded donor OH-stretching modes of (PhOH) 2 [11,13] and (PhOH) 3 [11], respectively, in agreement with the observed concentration dependency, the effects of annealing and reported band origins from previous neon [11] and argon matrix isolation [14] and jet spectroscopy investigations [15,16].The band origins for phenol monomer and phenol dimer embedded in neon have recently been revised with the latest observations at 3655 and 3515 cm −1 , respectively [13].The donor OH-stretching mode of (PhOH) 2 has previously been observed at 3530 cm −1 in a jet REMPI investigation [15,16], whereas the dangling OH-stretching mode of (PhOH) 2 and the free OH-stretching mode of the monomer have previously been assigned at 3654 cm −1 and 3657 cm −1 in jet expansions, respectively.Three transitions associated with the OH-stretching modes of (PhOH) 3 have previously been assigned upon jet-cooling at 3394, 3441 and 3449 cm −1 suggesting a cyclic structure of (PhOH) 3 [15], which has subsequently been confirmed by jet microwave spectroscopy [17][18][19].

P h O H P h O H P h O H P h O H P h O H P h O H P h O H
In the 400-800 cm −1 spectral range normally associated with the class of largeamplitude (hindered) torsional motions, several additional broad but still distinct bands, with intensities depending on both sample concentration and matrix annealing are observed.In order to distinguish the size of the cluster molecules, the concentration dependence and growth rate upon annealing of these bands were explored in additional experiments employing the phenol sublimation temperature in the 0-10 • C range during deposition.Figure 2a shows the spectral sampling of the doped neon matrix below 800 cm −1 at several different spatial locations limited to a diameter of approximately 3.5 mm of the infrared probe beam.This provides a variation of the local mixing ratio between neon and phenol molecules within the same matrix.The individual spectra have been normalized using the phenol monomer transition intensities in the 650-575 cm −1 range, which are not overlapped by the cluster molecule features.Figure 2b shows spectra below 800 cm −1 recorded before (black trace) and during annealing of the matrix at different steps (colored traces), with the latter calculated by subtracting the cold absorbance spectrum from the annealing spectra.The experimental findings from the combined spatial spot probing and annealing of the doped neon matrices allow the identification of two different sets of bands belonging to phenol cluster molecules.The most distinct absorption feature below 800 cm −1 , behaving as a cluster molecule band even for very low phenol/neon mixing ratios and always appearing in the spectra prior to annealing, is located at 605 cm −1 .Its band intensity correlates well with the OH-stretching band at 3515 cm −1 and this transition is therefore assigned to the only large-amplitude hydrogen bond librational (strongly hindered OH torsional) fundamental expected for (PhOH) 2 .The two other distinct bands observed at 557 cm −1 and 781 cm −1 show higher stoichiometry compared to the assigned (PhOH) 2 features.The most intense far-infrared band located at 557 cm −1 reveals the same concentration and annealing dependencies as the mid-infrared OH-stretching band at 3431 cm −1 assigned for (PhOH) 3 .The intensity of the band at 781 cm −1 similarly shows strong correlation with the 557 cm −1 and 3431 cm −1 bands with respect to the phenol mixing ratio and annealing procedures; however, in the high concentration experiments the spectral overlap with some broad satellite bands around 762 cm −1 impedes the exact correlations, tentatively assigned to CH out-of-plane bending motions of (PhOH) 2 /(PhOH) 3 .These exclusively experimental considerations suggest that the linked 557 cm −1 and 781 cm −1 bands should be assigned to transitions associated with large-amplitude hydrogen bond librational motion of the most stable conformation of (PhOH) 3 , which will be supported by complementary theoretical predictions in the following.
Figure 3 shows the optimized molecular geometries of phenol and its hydrogenbonded cluster molecules (PhOH) 2 and (PhOH) 3 using various theoretical methodologies, and provides the relative zero-point energy corrected ground-state dissociation energies for the four most stable conformations predicted for the (PhOH) 3 system.The computational predictions of the molecular geometry for the (PhOH) 2 system remain challenging even for modern computer resources.The MP2 quantum chemical method, which in general tends to provide consistent ab initio predictions for conventional hydrogen-bonded cluster molecules, runs into an issue in the treatment of the intermolecular interactions between the aromatic rings of the phenol fragments.Previously, a qualitatively correct structure has only been obtained with this methodology in conjunction with counterpoise corrections to account for the small basis set [10].Considering the large basis set required to achieve accurate vibrational frequency predictions, this theoretical approach becomes prohibitively expensive.The optimized geometry of (PhOH) 2 has been calculated using the MP2/AVQZ methodology (Figure 3a), which has shown significant underestimation of the hydrogen bonding in the complex, compared to what is observed in the experiment.Even with this relatively large basis set, the conventional MP2 method suffers from a strong overbinding effect between the aromatic rings, which in this case directly competes with the hydrogen bond.In order to achieve the desired accuracy for all cluster molecules investigated, the empirically scaled SCS-MP2 method and the dispersion-corrected hybrid PW6B95-D4 DFT functional have been selected, both known to perform comparatively well in the description of intermolecular non-covalent interactions [7].These methods have reproduced experimental data with reasonable consistency; however, there is variance in the predicted dissociation energies and harmonic vibrational frequencies between the two different approaches.
The predicted harmonic frequencies of vibrational modes for the conformations of (PhOH) 2 , (PhOH) 3 and PhOH•H 2 O are given in Table S1 in the Electronic Supplementary Information (ESI).For the (PhOH) 2 system, for which only one conformation has been optimized (Figure 3b), both methods predict one highly IR-active transition where the band associated with the hydrogen bond librational mode is experimentally observed in the present work.For the even more challenging (PhOH) 3 system, four different conformations have been optimized.Both methodologies show a large energy gap between the most stable cyclic pseudo-C 3 symmetric barrel-like potential energy minimum structure involving three cooperative OH• • • O hydrogen bonds and the other three predicted energy minima (Figure 3c).The less stable conformations 2-4 differ in the arrangements of the aromatic rings, or substitute the conventional intermolecular OH• • • O hydrogen bonds for the non-conventional OH• • • π bond or the stacking of aromatic rings.The differences in the optimized geometries result in different rankings of the dissociation energy calculated for the less stable conformations (Figure 3d,f).The simulated vibrational spectra for the four different conformations of (PhOH) 3 are shown in Figure 4.In this simulation, the calculated harmonic frequencies of the normal modes have been scaled using a separate factor for the modes involving large-amplitude hindered torsional motion of the OH-groups, which are estimated from the experimental spectra.The predicted pseudo-C 3 -symmetric global potential energy minimum (conf.1) has two near-degenerate highly IR-active librational modes (Figure 5b) estimated around 550 cm −1 .Owing to the pseudo-C 3 symmetry, there are no other predicted intense transitions associated with hydrogen bond librational modes in the 500-700 cm −1 range, unlike the asymmetric conformations 2-4, which all have at least two distinct IR-active hydrogen bond librational transitions in this region (indicated by filled colored band areas).In the case of conformation 2, one of the predicted hydrogen bond librational transitions has severe spectral overlaps with other less perturbed intramolecular transitions (DFT simulation); however, the two remaining highly IR-active hydrogen bond librational transitions should still be observed clearly between 500 and 600 cm −1 for this conformation, which is not supported in the present experiments.The second distinct transition assigned experimentally for (PhOH) 3 in the present experimental work at 781 cm −1 is also supported by the DFT and SCS-MP2 quantum chemical predictions.While the positions of the intramolecular transitions seem better predicted by the DFT approach, both methods agree on the position of the third IR-active transition associated with the concerted hydrogen bond librational motion of the most stable C 3 conformation of (PhOH) 3 (Figure 5c), involving all three O-H• • • O motifs, which is close to the experimental value.The relative transition intensity ratio between the two hydrogen bond librational transitions at 557 and 781 cm −1 of (PhOH) 3 is also qualitatively reproduced, although the overlap with the transitions attributed with perturbed intramolecular out-of-plane CH bending modes of (PhOH) 2 and (PhOH) 3 at higher concentrations makes more accurate analyses impossible.Together with the near-degenerate donor OH-stretching transitions observed at 3431 cm −1 , the two transitions associated with hydrogen bond librational motion allow us to unambiguously conclude that only the most stable pseudo-C 3 barrel-like conformation is observed in the present neon matrix experiments, both in the lower concentration regime and after being formed via annealing under cryogenic conditions.This is in strong contrast to previous findings for cyclic trimers of aliphatic alcohols, which have been shown to form "chair" structures due to steric hindrance [24,30].

The Phenol Monohydrate
Figure 6 shows the spectra collected for neon matrices doped with pure phenol, pure H 2 O and pure D 2 O together with spectra of phenol/H 2 O and phenol/D 2 O mixtures in the OH-stretching range (3400-3700 cm −1 ) and the relevant range for the transitions associated with large-amplitude hydrogen bond librational motion (500-670 cm −1 ).In between the hydrogen-bonded OH-stretching transitions for (PhOH) 2 , (H 2 O) 2 and (H 2 O) 3 , it is evident that the simultaneous deposition of phenol and H 2 O reveal a new strongly IR-active distinct band at 3499 cm −1 , which has recently been assigned to the phenol monohydrate PhOH•H 2 O [13].The corresponding neon matrix isolation experiments for phenol/D 2 O show that the assigned signal is indeed due to most stable conformation of this mixed cluster molecule, as the band shifts only slightly (and broadens somewhat due to spectral overlap with traces of the phenol-HDO complex), whereas the less stable conformation with water as hydrogen bond donor to phenol would absorb in the entirely different OD stretching region.This is supported by the quantum chemical results showing that the predicted zero-point energy corrected dissociation energy D 0 of the global potential energy minimum (Figure 7a) is 9 kJ•mol −1 higher than predicted for the local minimum (Figure 7b).The donor OH-stretching transitions of the phenol hydrates PhOH•(H 2 O) 1-3 have previously been observed in jet expansions at 3524 cm −1 , 3388 cm −1 and 3345 cm −1 [20], respectively, and the OH-stretching transition of the phenol monohydrate embedded in argon matrices has been assigned at 3461 cm −1 [14].In the spectral region below 700 cm −1 , a highly IR-active distinct band is observed at 642 cm −1 in the neon matrices doped simultaneously with phenol and H 2 O. Furthermore, the band intensity responds significantly to annealing, indicating complex formation events in the neon matrix.The band shifts slightly down to 638 cm −1 in the experiment where phenol is deposited together with D 2 O.The small isotopic shift of this transition confirms the proposed structure of the phenol monohydrate PhOH•H 2 O where the phenol molecule acts as the hydrogen bond donor and the observed transition is straightforwardly assigned to the large-amplitude donor librational motion of PhOH•H 2 O as visualized in Figure 5d.We do not observe any features in the low-energy part of the spectrum which could be clearly attributed to the less stable conformation.

Supporting Theoretical Analysis
The results of the LED analyses for the phenol cluster molecules PhOH•H 2 O, (PhOH) 2 and (PhOH) 3 are summarized in Table 1.The optimized potential energy minima structures obtained from the PW6B95-D4/ma-def2-QZVP and SCS-MP2/AVTZ methodologies were used for the subsequent LED analyses at the DLPNO-CCSD(T)/AVQZ level of theory.The electronic equilibrium dissociation energies D e and the zero-point corrected dissociation energies D 0 obtained from the higher DLPNO-CCSD(T) level on the optimized geometries are denoted D e (CC) and D 0 (CC), respectively.In the case of the smaller (PhOH) 2 and PhOH•H 2 O systems, geometry optimizations were also performed using the highest feasible SCS-MP2/AVQZ level to access the potential influence of the basis set size.The basis set does not have a significant influence on the dissociation energies for (PhOH) 2 as seen from the minor 0.3 kJ•mol −1 (0.1 %) increase in D e (CC), when calculated on the optimized geometry employing the larger AVQZ basis set.However, the molecular properties obtained directly from the method, specifically the value of change in zero-point energy upon complexation ∆ZPE, is much more dependent on the basis set.In that respect, the cheaper DFT approach performs better for (PhOH) 2 , likely due to lower basis set requirements.In general, when approaching the basis set limit, the SCS-MP2 method tends to underestimate non-covalent interactions, as can be seen from the results for the PhOH•H 2 O system.Still, the SCS-MP2 wave function approach provides more consistent results when dealing with the larger (PhOH) 3 conformations, when evaluated against the high-level DLPNO-CCSD(T) energies based on the same respective optimized molecular geometries probably due the larger contributions from London dispersion forces.
Table 1.The electronic equilibrium dissociation energies (D e ) for PhOH•H 2 O, (PhOH) 2 and (PhOH) 3 calculated both directly using the theoretical approaches PW6B95-D4/ma-def2-QZVP (denoted as DFT) and SCS-MP2 and the higher DLPNO-CCSD(T)/AVQZ level (D e (CC)) on the respective optimized geometries, the change of zero-point energy (∆ZPE), the zero-point energy corrected dissociation energies (D 0 , D 0 (CC)), the interfragment dispersion energies (∆E disp, frag ) and the sums of the dispersive (∑E disp ) and non-dispersive parts of D e (∑E N/D , given as D e (CC)− ∑E disp [31,32]).The SCS-MP2 results for PhOH•H 2 O are obtained using the AVQZ basis set.All values are given in units of kJ•mol −1 .The interfragment dispersion energies, obtained from the LED analysis, highlight the primary differences between the intermolecular forces dominating in the PhOH•H 2 O and the (PhOH) 2 systems.Both for the DFT and the SCS-MP2 approach, the resulting electronic equilibrium and zero-point energy corrected dissociation energies are almost identical for the global potential energy minima conformations.However, whereas the interfragment dispersion energy of (PhOH) 2 constitutes as much as 72% of the D e -value, this contribution constitutes only 29% for the less electron-rich PhOH•H 2 O system.

DFT SCS-MP2
While the increase in the basis set size in the calculated SCS-MP2 model of (PhOH) 2 leads to only a minor change in the cumulative value of D e , it causes over 2 kJ•mol −1 difference in the dispersion energy, suggesting a noticeable change in the predicted balance of intermolecular forces and the strong competition with the hydrogen bond.When comparing the predictions for the (PhOH) 2 to (PhOH) 3 systems with the AVTZ basis, the interfragment dispersion energy of the (PhOH) 3 system is lower by 1.3 kJ•mol −1 , and the proportion of E disp in the total dissociation energy is also lowered to 67%.The distance between the aromatic rings increases, e.g., the calculated interatomic distance between C-atoms in the para-positions increases from 6 to 6.25 Å.Despite the cooperative effect, the hydrogen bond length increases from 1.93 to 1.98 Å, while at the same time due to mutual rotation of the rings, the hydrogen in the meta-position becomes oriented towards the π-cloud of the adjacent fragment.The same can be observed for the DFT geometries.The DFT structures overall show significantly lower weighting of the dispersive interactions for the (PhOH) 2 and (PhOH) 3 structures.
The effect of hydrogen bond cooperativity is common for cyclic cluster molecules and can be probed experimentally by the relative average hydrogen bond librational energy compared to the single hydrogen bond librational transition energy of the respective dimer [24,30,33].In the case of (PhOH) 3 , we observe experimental evidence of cooperativity in the larger average observed hydrogen bond librational energy of 631 cm −1 (assuming the assigned 557 cm −1 transition is degenerate) relative to the observed librational band at 605 cm −1 for (PhOH) 2 .It is worth noting that only the high-level wavefunction-based methods properly capture this effect.The DFT approach predicts an almost exactly three-times larger dissociation energy for the (PhOH) 3 system with three cooperative intermolecular hydrogen bonds relative to the dissociation energy of (PhOH) 2 .The high-level DLPNO-CCSD(T) predictions of the DFT optimized molecular structures capture the experimental findings that D 0 ((PhOH) 3 ) > 3×D 0 ((PhOH) 2 ) and so does the SCS-MP2 approach even before the DLPNO-CCSD(T) refinement.
While we cannot quantitatively isolate the CH• • • π interactions using the LED approach, we can roughly estimate the relative roles of this contact between the phenol cluster molecules via selected structural parameters.In the (PhOH) 2 system, the O-H• • • O angle is 127 • , and the planes of the aromatic rings are oriented near perpendicular around the hydrogen bond axis.The bond angle in the (PhOH) 3 system is 90 • .The structure of (PhOH) 4 was calculated using the DFT approach to extend this series (see ESI), and the bond angle is 103 • for this system.The CH• • • π interaction is estimated as a function of the distance from the hydrogen in the ortho-position to the center of the adjacent phenol fragment.This distance is 3.75 Å, 3.30 Å and 3.37 Å for (PhOH) 2-4 , respectively.The relatively unobstructed rotation of the rings in the dimer results in optimal orientation that does not benefit a strong CH• • • π contact.In the (PhOH) 3 and (PhOH) 4 systems, the cooperative hydrogen bond networks constrict the mutual orientation of the rings, and therefore the hydrogen bond angles become smaller, but this orientation favors the CH• • • π interaction.The CH• • • π bond distance is shorter in (PhOH) 3 , but then increases again when transitioning to (PhOH) 4 .This appears to be due to strengthening of the cooperative hydrogen bond network resulting in an opening of the phenyl ring "barrel"-like structure.
The experimental setup consists of a Bruker Vertex 80V Fourier-transform spectrometer used with a DE-204 4 K (Advanced Research Systems, Inc., Macungie, PA, USA) closedcycle cryocooler.A transmission sample holder with a wedged diamond cold window (Diamond Materials, GmbH, Freiburg im Breisgau, Germany) which can be cooled to 4.5-5.0K is mounted onto the cold plate of the cooler.The cryocooler is surrounded by a rotatable vacuum shroud equipped with two wedged diamond windows and a gate valve allowing the insertion of the inlet system into the vacuum space of the cryostat.The wedged diamond windows minimize spectral interference fringes from internal reflections in the windows.The optical setup of the spectrometer consisted of a globar radiation source, a germanium-coated KBr beam splitter together with a broad-band HgCdTe detector for the complete infrared spectral range (4000-450 cm −1 ).A spectral resolution of 0.6 cm −1 was used throughout the measurements.
The inert gas matrix is obtained by simultaneous deposition of neon gas (99.999%,Air Liquide Danmark A/S, Taastrup, Denmark) with vapors of the samples through separate inlet tubes, which are brought to within 5 mm of the cold sample window using a motorized stage.The neon is supplied from an MKS G-series mass flow controller, passing a 1 mlong LN 2 -cooled coiled trap for impurities, at a flow rate of 6-8 sccm depending on the desired mixing ratio.Water vapor was supplied from a variable-pressure (30-80 Pa) supply volume through a small-flow Swagelok metering valve, and phenol vapor was deposited directly from the sublimation vessel thermostated in the temperature range −5-10 ± 0.2 • C throughout the deposition process.
Following acquisition of the initial ("pre-annealing") matrix spectra, the samples were annealed by raising the temperature to 9 K for 60 min using a LakeShore controller operating a resistive heater and a Si-diode temperature sensor attached to the cold head.Subsequently, the heater was deactivated, allowing the matrix to cool again to 4-5 K before recording the "post-annealing" spectra.This annealing procedure softens the neon matrix and allows the diffusion of phenol molecules within the matrix and triggers the self-association of phenol into hydrogen-bonded cluster molecules.
The deposition of the samples via separate inlet tubes creates intentional spatial inhomogeneity in the phenol/water/neon mixing ratio, which is exploited to enable spatial "spot-to-spot" matrix probing and differentiate between cluster sizes in the same experiment.The effective diameter of the probed region of the matrix is comparable to the employed 3.5 mm aperture size.The background spectra were all collected of the evacuated cryostat at room temperature.
The electronic equilibrium dissociation energies (D e ) of the cluster molecules were obtained by performing domain-based local pair natural orbital coupled cluster (DLPNO-CCSD(T)) [45,46] single-point calculations on the respective DFT and SCS-MP2 geometries.The cut-offs for the DLPNO methodology were set by the TightPNO keyword [47].The zero-point corrected ground-state dissociation energies (D 0 ) of the cluster molecules were calculated as the sum of the DLPNO-CCSD(T) electronic dissociation energies (D e ) and the ∆ZPE values from the respective base calculations.
The intermolecular London dispersion energies of the cluster molecules were extracted using the DLPNO-based local energy decomposition (LED) scheme [31].The procedure is described in the literature [31,32].
The resolution of the identity approximation was used in all of the calculations with the appropriate fitting basis sets [48][49][50][51].The numerical "chain-of-spheres" approximation (COSX [52,53]) was used with the MP2 level calculations.The calculations were all carried out using the DTU Computing Center (DCC) cluster [54].

Conclusions
The self-association and microhydration molecular recognition mechanisms of phenol have been explored experimentally by means of the less accessible vibrational transitions associated with large-amplitude hydrogen bond librational motion.The present work provides for the first time unambiguous spectroscopic assignments for the large-amplitude and highly anharmonic hydrogen bond librational modes of the phenol dimer (605 cm −1 ), the phenol trimer (557 and 781 cm −1 ) and two isotopic variants of the phenol monohydrate (642 cm −1 for PhOH•H 2 O and 638 for PhOH•D 2 O) embedded in inert cryogenic "quantum" neon matrices at 4 K.While accurate microwave spectroscopic information concerned with the molecular structures of these systems has been reported in the literature previously, these new experimental findings provide important constraints for future high-level quantum chemical models of the corresponding intermolecular potential energy surfaces.
The experimental findings have been complemented by extensive conformation search sampling and further computational optimizations employing both dispersion-corrected DFT and wavefunction ab initio methodologies of the most stable potential energy minima.In addition, the local energy decompositions of the total interaction energies for these

0
°C P o s t -a n n e a l i n g 0 °C A s -d e p o s i t e d − 5 °C s a m p

Figure 1 .
Figure 1.The infrared absorption spectra of phenol embedded in solid neon at 4 K.In the two different experiments shown, the phenol source sublimation vessel was thermostated at 0 • C (bottom) and at −5 • C (top and bottom, respectively) during deposition to achieve different mixing ratios with neon.The spectra of the 0 • C experiment were recorded before (blue trace) and after annealing of the matrix at 9.5 K (red trace).The spectra have been normalized to the phenol monomer transitions in the 650-575 cm −1 range.The experimentally observed transitions associated with the intramolecular OH-stretching and the large-amplitude hydrogen bond librational modes of (PhOH) 2 and (PhOH) 3 are indicated with the respective band origins.

Figure 2 .
Figure 2. (a)The spectral dependence of the phenol/neon mixing ratio by spatial spot probing.The three shown spectra collected for spots with increasing phenol/ratio (the blue and red traces collected for the lowest and highest phenol/neon mixing ratios, respectively) have been normalized to the monomer transitions.The proposed dimer and trimer absorption features are differentiated based on their different growth rates relative to the monomer bands.(b) The evolution of the neon matrix spectra over time during annealing at 9.5 K.A series of difference spectra (annealing spectrum subtracted the pre-annealing spectrum) collected at 3 min intervals during the annealing procedure is shown.The colored difference spectra (the blue trace collected after 3 min and the red trace collected after 15 min, respectively) show excess absorption relative to the pre-annealing spectrum (black trace) demonstrating the progressive formation of phenol cluster molecules.

Figure 6 .
Figure 6.The infrared absorption spectra of phenol (black trace, small traces of H 2 O), H 2 O (purple trace) and D 2 O (brown trace) together with spectra of phenol/H 2 O (blue trace) and phenol/D 2 O (green trace) mixtures.For the two mixtures the spectra after annealing of the matrix to 9.5 K (red and orange traces, respectively) are shown below 650 cm −1 .The new assigned transitions associated with the large-amplitude hydrogen bond librational modes of the PhOH•H 2 O and PhOH•D 2 O monohydrates are indicated in the spectra with their respective band positions.

2 Figure 7 .
Figure 7.The optimized global (a) and local (b) potential energy minima structures of the phenol monohydrate (PhOH•H 2 O) employing the SCS-MP2/AVQZ level of theory.