Regularities and Anomalies in Neon Matrix Shifts of Hydrogen-Bonded O-H Stretching Fundamentals

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.

Exploring structures of small anionic nickel-ethanol clusters with infrared spectroscopy

Small anionic nickel clusters with ethanol are investigated with a combination of mass-selective infrared photodissociation spectroscopy in a molecular beam and density functional theory simulations at the BLYP/6-311g(d,p) and TPSSh/def2-TZVPP level. In this context, the O-H stretching vibration of the ethanol is analyzed to obtain information about the structural motif, the geometry of the metal core, and the spin state of the clusters. For the [Ni2(EtOH)]- and [Ni3(EtOH)]- clusters, we assign quartet states of motifs with a hydrogen bond from the ethanol to the linear nickel core. The aggregation of a further ethanol molecule, yielding the [Ni3(EtOH)2]- cluster, results in the formation of a cooperative hydrogen bond network between the nickel core and the two ethanol molecules.

The impact of solvation on the structure and electric field strength in Li+GlyGly complexes

To scrutinise the impact of electric fields on the structure and vibrations of biomolecules in the presence of water, we study the sequential solvation of lithium diglycine up to three water molecules with cryogenic infrared action spectroscopy. Conformer-specific IR-IR spectroscopy and H2O/D2O isotopic substitution experiments provide most of the information required to decipher the structure of the observed conformers. Additional confirmation is provided by scaled harmonic vibrational frequency calculations using MP2 and DFT. The first water molecule is shown to bind to the Li+ ion, which weakens the electric field experienced by the peptide and as a consequence, also the strength of an internal NH⋯NH2 hydrogen bond in the diglycine backbone. The strength of this hydrogen bond decreases approximately linearly with the number of water molecules as a result of the decreasing electric field strength and coincides with an increase in the number of conformers observed in our spectra. The addition of two water molecules is already sufficient to change the preferred conformation of the peptide backbone, allowing for Li+ coordination to the lone pair of the terminal amine group.

Intermolecular Hydrogen Bonding in Associated Fluids: The Case of Isopentyl Alcohol Dissolved in Carbon Tetrachloride

Fourier-transform infrared (FTIR) spectra of isopentyl-alcohol dissolved in carbon tetrachloride (CCl4) were recorded as a function of concentration and temperature. Dilute isopentyl alcohol/CCl4 solutions were prepared in alcohol at concentrations of 1, 0.5, 0.3, 0.2, 0.1, 0.05, 0.02, 0.01, 0.005, 0.001 and 0.0005 M. Infrared absorption measurements were taken within a temperature range of 17-67 °C below the boiling point of the solutions. Decomposition of the spectral features corresponding to associated and unassociated species was performed to quantitatively follow the effect of temperature and concentration on intermolecular hydrogen bonding (HB) in isopentyl alcohol. The spectral feature in the 3600-3650 cm-1 frequency range attributed to the free OH stretching band was studied in detail to determine changes based on concentration and temperature variations. Computational methodologies were applied to evaluate the energetics and vibrational properties of the species involved in the structure in the gaseous state where no interactions are present. The results are discussed in view of relevant structural models to gain quantitative information concerning the effect of concentration and temperature on intermolecular hydrogen bonding.

The first HyDRA challenge for computational vibrational spectroscopy

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

Raising the benchmark potential of a simple alcohol-ketone intermolecular balance

2-Butanone offers two hydrogen bond docking variants to a solvating methanol which are cleanly separated by supersonic jet infrared absorption spectroscopy in the OH-stretching range, resolving earlier action spectroscopy indeterminacies for this elementary case of an intermolecular alcohol-ketone balance. The solvent preference for the shorter chain side is unambiguously derived from the spectra of homologous compounds. It is analysed in terms of competing steric and dispersion interactions and the resulting energy differences across a low interconversion barrier. Fortuitous cancellations are discussed and quantitative energy deficiencies of the employed DFT approaches are suggested. Some benchmarkable experimental observations: at low temperature, a single methanol molecule prefers the methyl-sided oxygen lone pair of 2-butanone over the ethyl-sided lone pair by 1-2 kJ mol^-1, the trans butane backbone is conserved in both low-lying isomers, the OH-stretching fundamentals differ by 47(2) cm^-1.

Chromone-methanol clusters in the electronic ground and lowest triplet state: a delicate interplay of non-covalent interactions

Chromone offers two energetically almost equivalent docking sites for alcohol molecules, in which the hydroxyl group is hydrogen bonded to one of the free electron pairs of the carbonyl O atom. Here, the delicate balance between these two competing arrangements is studied by combining IR/R2PI and UV/IR/UV spectroscopy in a molecular beam supported by quantum-chemical calculations. Most interestingly, chromone undergoes an efficient intersystem crossing into the triplet manifold upon electronic excitation, so that the studies on aromatic molecule-solvent complexes are for the first time extended to such a cluster in a triplet state. As the lowest triplet state (T₁) is of ground state character, powerful energy decomposition approaches such as symmetry-adapted perturbation theory (SAPT) and local energy decomposition using the domain-based local pair natural orbital coupled-cluster method (DLPNO-CCSD(T)/LED) are applied. From the theoretical analysis we infer for the T₁ state a loss of planarity (puckering) of the 4-pyrone ring of the chromone unit, which considerably affects the interplay between different types of non-covalent interactions at the two possible binding sites.

Hydrogen Delocalization in an Asymmetric Biomolecule: The Curious Case of Alpha-Fenchol

Rotational microwave jet spectroscopy studies of the monoterpenol α-fenchol have so far failed to identify its second most stable torsional conformer, despite computational predictions that it is only very slightly higher in energy than the global minimum. Vibrational FTIR and Raman jet spectroscopy investigations reveal unusually complex OH and OD stretching spectra compared to other alcohols. Via modeling of the torsional states, observed spectral splittings are explained by delocalization of the hydroxy hydrogen atom through quantum tunneling between the two non-equivalent but accidentally near-degenerate conformers separated by a low and narrow barrier. The energy differences between the torsional states are determined to be only 16(1) and 7(1) cm-1hc for the protiated and deuterated alcohol, respectively, which further shrink to 9(1) and 3(1) cm-1hc upon OH or OD stretch excitation. Comparisons are made with the more strongly asymmetric monoterpenols borneol and isopinocampheol as well as with the symmetric, rapidly tunneling propargyl alcohol. In addition, the third-in contrast localized-torsional conformer and the most stable dimer are assigned for α-fenchol, as well as the two most stable dimers for propargyl alcohol.

Subtle hydrogen bonds: benchmarking with OH stretching fundamentals of vicinal diols in the gas phase

The theoretical description of spectral signatures for weakly bound hydrogen contacts between alcohol groups is challenging and remains poorly characterised. By combining Raman jet spectroscopy with appropriately scaled harmonic DFT predictions and relaxation path analyses for 16 vicinal diols (ethylene glycol (ethane-1,2-diol), propane-1,2-diol, 3,3,3-trifluoro-propane-1,2-diol, rac-butane-2,3-diol, 2-methyl-propane-1,2-diol, 2-methyl-butane-2,3-diol, pinacol (2,3-dimethyl-butane-2,3-diol), 3-butene-1,2-diol, 1-phenyl-ethane-1,2-diol, trans-cyclobutane-1,2-diol, trans-cyclopentane-1,2-diol, trans-cyclohexane-1,2-diol, trans-cycloheptane-1,2-diol, cis-cyclohexane-1,2-diol, 1-(1-hydroxy-1-methylethyl)-cyclopentanol and [1,1'-bicyclopentyl]-1,1'-diol), 69 conformational assignments become possible in a two-tier approach with a 5 diol training and an 11 diol test set. The latter reveals systematic deviations for ring strain and secondary π interactions, but otherwise a remarkably robust correction and correlation model based on hybrid DFT with a minimally augmented triple-zeta basis set is obtained, whereas GGA functionals perform significantly worse. Raw experimental data in the 3560-3700 cm^-1 wavenumber range as well as computed geometries of all conformations invite further vibrational and structural benchmarking at the onset of hydrogen bonding. Beyond this diol-probed threshold, the accurate prediction of hydrogen bond induced shifts of different magnitudes remains one of the challenges for DFT functionals.

Predicting OH stretching fundamental wavenumbers of alcohols for conformational assignment: different correction patterns for density functional and wave-function-based methods

A model is presented for the prediction of OH stretching fundamental wavenumbers of alcohol conformers in the gas phase by application of a small set of empirical anharmonicity corrections to calculations in the harmonic approximation. In contrast to the popular application of a uniform scaling factor, the local chemical structure of the alcohol is taken into account to greatly improve accuracy. Interestingly, different correction patterns emerge for results of hybrid density functional (B3LYP-D3 and PBE0-D3) and wave-function-based methods (SCS-LMP2, LCCSD(T*)-F12a and CCSD(T)-F12a 1D). This raises questions about electronic structure deficiencies in these methods and differences in anharmonicity between alcohols. After its initial construction on the basis of literature assignments the model is tested with Raman jet spectroscopy of propargyl alcohol, cyclohexanol, borneol, isopinocampheol and 2-methylbutan-2-ol. For propargyl alcohol a spectral splitting attributed to tunneling is resolved. PBE0-D3 is identified as a well performing and broadly affordable electronic structure method for this model. A mean absolute error of 1.3 cm-1 and a maximum absolute error of 3 cm-1 result for 46 conformers of 24 alcohols in a 60 cm-1 range, when a single parameter is adjusted separately for each alcohol substitution class (methanol, primary, secondary, tertiary).

Effect of solvents on the chiral recognition mechanisms of immobilized cellulose-based chiral stationary phase

The effect of solvents on the enantioselectivities of four structurally similar chiral solutes with a cellulose derivative-based chiral stationary phase, Chiralpak IB, were studied using acetone (AC), 2-propanol (IPA), and tert-butanol (TBA) separately as polar modifiers. The enantioselectivities α of benzoin and methyl mandelate decrease with an increase in modifier concentration C_M, whereas the enantioselectivity of pantolactone increased with increasing AC concentration. These results were attributed to the heterogeneous adsorption mechanisms of enantiomers. To interpret the dependence of enantioselectivity on modifier content, an enantioselectivity model based on a two-site adsorption model was proposed. The dependence of α on C_M was inferred to be mainly due to the distinct modulating effects of modifier concentration on the two adsorption sites: the nonselective type-I site and enantioselective type-II site. The model fitted the benzoin data satisfactorily over a wide TBA concentration range. The retention factors as a function of TBA concentration were successfully deconvoluted for each site. With the use of the proposed model, it was inferred that the chiral recognitions of benzoin and methyl mandelate were mainly achieved by the presence of an aromatic group adjacent to the hydroxyl group. When using IPA and TBA separately as modifiers, the presence of an aromatic group adjacent to the ketone group mainly contributed to the nonselective π interactions and enantioselective steric interactions, respectively. These results, along with those of the modifier adsorption isotherms, determined using the perturbation method, as well as the retention behaviors of various achiral solutes, indicate that the molecular recognition mechanism of IB sorbent is highly sensitive to the adsorbate's molecular geometry. The molecular environment of the sorbent can be controlled using different modifiers, leading to distinct adsorption and retention mechanisms.

Activity coefficients of binary methanol alcohol mixtures from cluster weighting

The hydrogen bond network of different small alcohols is investigated via cluster analysis. Methanol/alcohol mixtures are studied with increasing chain length and branching of the molecule. Those changes can play an important role in different fields, including solvent and metal extraction. The extended tight binding method GFN2-xTB allows the evaluation and geometry optimization of thousands of clusters built via a genetic algorithm. Interaction energies and geometries are evaluated and discussed for the neat systems. Thermodynamic properties, such as vaporization enthalpies and activity coefficients, are calculated with the binary quantum cluster equilibrium (bQCE) approach using our in-house code peacemaker 2.8. Combined distribution functions of the distances against the angles of the hydrogen bonds are evaluated for neat and mixed clusters and weighted by the equilibrium populations achieved from bQCE calculations.

Shifting formic acid dimers into perspective: vibrational scrutiny in helium nanodroplets

A metastable dimer of formic acid has been prepared inside superfluid helium nanodroplets and examined using IR spectroscopy and quantum chemical calculations.

Exploring the hydrogen bond kinetics of methanol–water solutions using Raman scattering

Stimulated Raman scattering was used to clearly show the hydrogen bond kinetics of water–methanol mixed solutions.

Solvation effects on wavenumbers and absorption intensities of the OH-stretch vibration in phenolic compounds – electrical- and mechanical anharmonicity via a combined DFT/Numerov approach

A previously measured oscillatory intensity pattern in phenolic compounds between different solvents was successfully reproduced for the first time, employing modern grid-based methods to solve the time-independent vibrational Schrödinger equation.

Conformational isomerizations triggered by vibrational excitation of second stretching overtones

Vibrational excitation using frequency-tunable IR laser light has been developed as a powerful tool for selective manipulation of molecular conformations. In this methodology, vibrational excitation has been typically applied to the first stretching overtones (∼80 kJ mol^-1) but also to the fundamental modes (∼40 kJ mol^-1). Here, we demonstrate that selective conformational isomerizations are also achieved using excitation to second stretching overtones (∼120 kJ mol^-1). The extremely weak absorptions of the second stretching overtones of molecules isolated in low-temperature matrices were measured for the first time; here using three prototype molecules: hydroxyacetone (HA), glycolic acid (GAc) and glycolamide (GAm). Benchmarking of computed anharmonic IR spectra showed that the B3LYP/SNSD method provides the best agreement with experimental frequencies of the ν(OH), 2ν(OH) and 3ν(OH) modes for the studied molecules in argon matrices. Selective irradiation at the 3ν(OH) frequencies (9850-10 500 cm^-1) of HA, GAc and GAm monomers in argon matrices at 15 K successfully triggers their conformational isomerization. These results open the door to extend control over conformations separated by higher barriers and to induce other transformations not energetically accessible by excitation to the fundamental or first stretching overtone modes.

C-D Vibration at C2 Position of Imidazolium Cation as a Probe of the Ionic Liquid Microenvironment

Unlike molecular solvents, imidazolium-based ionic liquids are entirely made of ions with spatial heterogeneity. There is a need for spectroscopic probes that can assess the microenvironment near the cations of these complex liquids. In this manuscript, we describe simple chemical procedures to label the C2 position of imidazolium cation with a C-D vibrational probe and show, through linear and nonlinear vibrational spectroscopies, that this C-D stretching mode can be a useful analytical tool to assess both the solvent microenvironment and solute-solvent interactions in imidazolium-based ionic liquids from the cation point of view. It is expected that this C-D vibration probe on the cation will lead to the development of innovative experimental strategies that can provide a better understanding of such ionic liquids.

Cooperativity in Alcohol-Nitrogen Complexes: Understanding Cryomatrices through Slit Jet Expansions

FTIR spectroscopy of supersonic expansions is used to characterize alcohol dimers with one, two, and several nitrogen molecules attached to them. The nitrogen coating causes progressive spectral downshifts of the OH stretching fundamentals which are related to and explain matrix isolation shifts. Comparison of methanol, tert-butyl alcohol and ethanol as well as deuteration of methanol assist in the assignment. Alcohol monomers and trimers are significantly more resistant to nitrogen coating due to a lack of cooperativity and dangling bonds, respectively. In the case of ethanol, the role of conformational isomerism and combination bands is further elucidated. The experimental findings help rationalize the anomalously small OH stretching dimerization shift of methanol in the gas phase, in comparison to that of tert-butyl alcohol.

Infrared absorption of methanol-water clusters (CH3OH)n(H2O), n = 1-4, recorded with the VUV-ionization/IR-depletion technique

We recorded infrared (IR) spectra in the CH- and OH-stretching regions of size-selected clusters of methanol (M) with one water molecule (W), represented as M_nW, n = 1-4, in a pulsed supersonic jet using the photoionization/IR-depletion technique. Vacuum ultraviolet emission at 118 nm served as the source of ionization in a time-of-flight mass spectrometer to detect clusters M_nW as protonated forms M_n-1WH⁺. The variations in intensities of M_n-1WH⁺ were monitored as the wavelength of the IR laser light was tuned across the range 2700-3800 cm^-1. IR spectra of size-selected clusters were obtained on processing of the observed action spectra of the related cluster-ions according to a mechanism that takes into account the production and loss of each cluster due to IR photodissociation. Spectra of methanol-water clusters in the OH region show significant variations as the number of methanol molecules increases, whereas those in the CH region are similar for all clusters. Scaled harmonic vibrational wavenumbers and relative IR intensities predicted with the M06-2X/aug-cc-pVTZ method for the methanol-water clusters are consistent with our experimental results. For dimers, absorption bands of a structure WM with H₂O as a hydrogen-bond donor were observed at 3570, 3682, and 3722 cm^-1, whereas weak bands of MW with methanol as a hydrogen-bond donor were observed at 3611 and 3753 cm^-1. For M₂W, the free OH band of H₂O was observed at 3721 cm^-1, whereas a broad feature was deconvoluted to three bands near 3425, 3472, and 3536 cm^-1, corresponding to the three hydrogen-bonded OH-stretching modes in a cyclic structure. For M₃W, the free OH shifted to 3715 cm^-1, and the hydrogen-bonded OH-stretching bands became much broader, with a weak feature near 3179 cm^-1 corresponding to the symmetric OH-stretching mode of a cyclic structure. For M₄W, the observed spectrum agrees unsatisfactorily with predictions for the most stable cyclic structure, indicating significant contributions from branched isomers, which is distinctly different from M₅ of which the cyclic form dominates.

Matrix isolation FTIR study of hydrogen-bonded complexes of methanol with heterocyclic organic compounds

Hydrogen bonded complexes of heterocyclic compounds with methanol were studied using matrix isolation FTIR spectroscopy and theoretical calculations.

Hydrogen bonding in alcohol–ethylene oxide and alcohol–ethylene sulfide complexes

The O–H⋯O and O–H⋯S hydrogen bonds are of similar strength in the corresponding alcohol–EO and alcohol–ES complexes.

The donor OH stretching–libration dynamics of hydrogen-bonded methanol dimers in cryogenic matrices

FTIR spectra of the methanol dimer trapped in neon matrices are presented.

Non-covalent interaction of benzene with methanol and diethyl ether solid surfaces

We have investigated the interactions involved at the interface of binary, layered ices (benzene on methanol and on diethyl ether) by means of laboratory experiments and ab initio calculations on model clusters.

OH-stretch overtone of methanol: empirical assignment using a two temperature technique in a supersonic jet

This paper describes a novel approach for empirical lower state assignments in complex high resolution ro-vibrational overtone spectra of molecules with low rotational constants and complex intramolecular dynamics. Methanol, CH3OH, was chosen as a representative of such molecules - it is an asymmetric top with two non-hydrogen nuclei and hindered internal rotation leading to dense and disordered rotational structure of vibrational overtone bands. We report the first rotationally resolved methanol spectra of the OH-stretch overtone 2ν1 band using sub-Doppler diode laser spectroscopy in a supersonic jet, and describe how the combination of two temperature analysis (TTA) and analysis by ground state combination differences (GSCDs) is used to reliably identify spectral lines that originate from lowest rotational states. In the first step of the analysis, the TTA was utilized to obtain a set of possible rotational assignments for each spectral line using the line intensity variation between two different temperatures in the supersonic jet (13, and 56 K respectively). Thereafter, the GSCDs were used to confirm specific lower state assignment for those spectral lines that have been identified to have low rotational ground states by the TTA. We show that the TTA pre-selection leads to fast and reliable confirmation by GSCDs and avoids false assignments due to accidental GSCD matches. The procedure yields an important subset of reliably assigned spectral lines in the complex ro-vibrational structure that provides a convenient starting point for subsequent application of traditional spectral analysis techniques.

Soft hydrogen bonds to alkenes: the methanol-ethene prototype under experimental and theoretical scrutiny

An FTIR spectroscopic study of the elusive hydrogen-bonded methanol-ethene complex, the most elementary example for weak intermolecular alcohol hydrogen bonding to a π cloud, is presented. By isolating the complex in a supersonic jet, the rigorous comparability to high-level quantum chemical calculations is ensured. In stark contrast to classical hydrogen bonds, experimental overtone analysis reveals the harmonic oscillator approximation for the OH red shift to be accurate. Harmonic calculations up to explicitly correlated local coupled-cluster level are thus found to agree very well with experiment. The experimental OH values for the red shift (45 cm^-1), the small change in diagonal anharmonicity (-3 cm^-1) and the overtone intensity attenuation (2 × 10²-fold) together with theoretical predictions for the preferred structural arrangement and the zero-point-corrected dissociation energy (8 kJ mol^-1) may thus be regarded as definitive reference values for related systems and for more approximate computational methods. In particular, MP2 calculations are shown to fail for this kind of weak intermolecular interaction.

To π or not to π – how does methanol dock onto anisole?

Supersonic jet spectroscopy reveals that methanol opts against aromatic docking onto anisole despite a very close competition.