Communication: Towards the binding energy and vibrational red shift of the simplest organic hydrogen bond: Harmonic constraints for methanol dimer

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.

Dispersion Control over Molecule Cohesion: Exploiting and Dissecting the Tipping Power of Aromatic Rings

ConspectusWe have learned over the past years how London dispersion forces can be effectively used to influence or even qualitatively tip the structure of aggregates and the conformation of single molecules. This happens despite the fact that single dispersion contacts are much weaker than competing polar forces. It is a classical case of strength by numbers, with the importance of London dispersion forces scaling with the system size. Knowledge about the tipping points, however difficult to attain, is necessary for a rational design of intermolecular forces. One requires a careful assessment of the competing interactions, either by sensitive spectroscopic techniques for the study of the isolated molecules and aggregates or by theoretical approaches. Of particular interest are the systems close to the tipping point, when dispersion interactions barely outweigh or approach the strength of the other interactions. Such subtle cases are important milestones for a scale-up to realistic multi-interaction situations encountered in the fields of life and materials science. In searching for examples that provide ideal competing interactions in complexes and small clusters, aromatic systems can offer a diverse set of molecules with a variation of dispersion and electrostatic forces that control the dominant and peripheral interactions. Our combined spectroscopic and theoretical investigations provide valuable insights into the balance of intermolecular forces because they typically allow us to switch the aromatic substituent on and off. High-resolution rotational spectroscopy serves as a benchmark for molecular structures, as correct calculations should be based on correct geometries. When discussing the competition with other noncovalent interactions, obvious competitors are directional hydrogen bonds. As a second counterweight to aryl interactions, we will discuss aurophilic/metallophilic interactions, which also have a strong stabilization with a small number of atoms involved. Vibrational spectroscopy is most sensitive to interactions of light atoms, and the competition of OH hydrogen bonds with dispersion forces in a molecular aggregate can be judged well by the OH stretching frequency. Experiments in the gas phase are ideal for gauging the accuracy of quantum chemical predictions free of solvent forces. A tight collaboration utilizing these three methods allows experiment vs experiment vs theory benchmarking of the overall influence of dispersion in molecular structures and energetics.

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.

A plethora of isomerization processes and hydrogen scrambling in the fragmentation of the methanol dimer cation: a PEPICO study

The valence photoionization of light and deuterated methanol dimers was studied by imaging photoelectron photoion coincidence spectroscopy in the 10.00-10.35 eV photon energy range. Methanol clusters were generated in a rich methanol beam in nitrogen after expansion into vacuum. They generally photoionize dissociatively to protonated methanol cluster cations, (CH₃OH)_nH⁺. However, the stable dimer parent ion (CH₃OH)₂⁺ is readily detected below the dissociation threshold to yield the dominant CH₃OH₂⁺ fragment ion. In addition to protonated methanol, we could also detect the water- and methyl-loss fragment ions of the methanol dimer cation for the first time. These newly revealed fragmentation channels are slow and cannot compete with protonated methanol cation formation at higher internal energies. In fact, the water- and methyl-loss fragment ions appear together and disappear at a ca. 150 meV higher energy in the breakdown diagram. Experiments with selectively deuterated methanol samples showed H scrambling involving two hydroxyl and one methyl hydrogens prior to protonated methanol formation. These insights guided the potential energy surface exploration to rationalize the dissociative photoionization mechanism. The potential energy surface was further validated by a statistical model including isotope effects to fit the experiment for the light and the perdeuterated methanol dimers simultaneously. The (CH₃OH)₂⁺ parent ion dissociates via five parallel channels at low internal energies. The loss of both CH₂OH and CH₃O neutral fragments leads to protonated methanol. However, the latter, direct dissociation channel is energetically forbidden at low energies. Instead, an isomerization transition state is followed by proton transfer from a methyl group, which leads to the CH₃(H)OH⁺⋯CH₂OH ion, the precursor to the CH₂OH-, H₂O-, and CH₃-loss fragments after further isomerization steps, in part by a roaming mechanism. Water loss yields the ethanol cation, and two paths are proposed to account for m/z 49 fragment ions after CH₃ loss. The roaming pathways are quickly outcompeted by hydrogen bond breaking to yield CH₃OH₂⁺, which explains the dominance of the protonated methanol fragment ion in the mass spectrum.

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.

Halogens in Acetophenones Direct the Hydrogen Bond Docking Preference of Phenol via Stacking Interactions

Phenol is added to acetophenone (methyl phenyl ketone) and to six of its halogenated derivatives in a supersonic jet expansion to determine the hydrogen bonding preference of the cold and isolated 1:1 complexes by linear infrared spectroscopy. Halogenation is found to have a pronounced effect on the docking site in this intermolecular ketone balance experiment. The spectra unambiguously decide between competing variants of phenyl group stacking due to their differences in hydrogen bond strength. Structures where the phenyl group interaction strongly distorts the hydrogen bond are more difficult to quantify in the experiment. For unsubstituted acetophenone, phenol clearly prefers the methyl side despite a predicted sub-kJ/mol advantage that is nearly independent of zero-point vibrational energy, turning this complex into a challenging benchmark system for electronic structure methods, which include long range dispersion interactions in some way.

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).

Understanding benzyl alcohol aggregation by chiral modification: the pairing step

A combination of linear infrared and Raman spectroscopy in supersonic slit jet expansions is used to clarify the conformational preferences in the dimer of the transiently chiral benzyl alcohol (phenylmethanol) under vacuum isolation. By experimentally exploring close analogies with the permanently chiral 1-phenylethanol, which is conformationally locked in the jet through intramolecular chirality induction, and by computational analysis of their conformational energy landscapes, several conclusions are drawn. The lowest energy dimer is confirmed to be cooperatively OHOHπ-bonded and shown to be homochiral. Its heterochiral counterpart is slightly higher in energy and can be spectrally assigned as a minor constituent. A metastable heterochiral OHπ/OHπ structure with weakly coupled hydrogen bonds is efficiently trapped behind a Ci symmetry-enhanced barrier and can be assigned by IR/Raman mutual exclusion. Its homochiral counterpart is kinetically less stable but might be addressed by rotational spectroscopy. Ratings of standard density functionals with a standard basis set in terms of reproducing these experimental chirality synchronization benchmarks are presented.

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.

Concerted Pair Motion Due to Double Hydrogen Bonding: The Formic Acid Dimer Case

Formic acid dimer as the prototypical doubly hydrogen-bonded gas-phase species is discussed from the perspective of the three translational and the three rotational degrees of freedom which are lost when two formic acid molecules form a stable complex. The experimental characterisation of these strongly hindered translations and rotations is reviewed, as are attempts to describe the associated fundamental vibrations, their combinations, and their thermal shifts by different electronic structure calculations and vibrational models. A remarkable match is confirmed for the combination of a CCSD(T)-level harmonic treatment and an MP2-level anharmonic VPT2 correction. Qualitatively correct thermal shifts of the vibrational spectra can be obtained from classical molecular dynamics in CCSD(T)-quality force fields. A detailed analysis suggests that this agreement between experiment and composite theoretical treatment is not strongly affected by fortuitous error cancellation but fully converged variational treatments of the six pair or intermolecular modes and their overtones and combinations in this model system would be welcome.

Soft experimental constraints for soft interactions: a spectroscopic benchmark data set for weak and strong hydrogen bonds

An experimental benchmark data base on rotational constants, vibrational properties and energy differences for weakly and more strongly hydrogen-bonded complexes and their constituents from the spectroscopic literature is assembled. It is characterized in detail and finally contracted to a more compact, discriminatory set (ENCH-51, for Experimental Non-Covalent Harmonic with 51 entries). The meeting points between theory and experiment consist of equilibrium rotational constants and harmonic frequencies and energies, which are back-corrected from experimental observables and are very easily accessible by quantum chemical calculations. The relative performance of B3LYP-D3, PBE0-D3 and M06-2X density functional theory predictions with a quadruple-zeta basis set is used to illustrate systematic errors, error compensation and selective performance for structural, vibrational and energetical observables. The current focus is on perspectives and different benchmarking methodologies, rather than on a specific theoretical method or a specific class of compounds. Extension of the data base in chemical, observable and quantum chemical method space is encouraged.

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.

Strained hydrogen bonding in imidazole trimer: a combined infrared, Raman, and theory study

This is not how three imidazole molecules prefer to arrange, as a combined IR, Raman and computational analysis unambiguously shows.

Reply to the 'Comment on "Methanol dimer formation drastically enhances hydrogen abstraction from methanol by OH at low temperature"' by D. Heard, R. Shannon, J. Gomez Martin, R. Caravan, M. Blitz, J. Plane, M. Antiñolo, M. Agundez, E. Jimenez, B. Ballesteros, A. Canosa, G. El Dib, J. Albaladejo and J. Cernicharo, Phys. Chem. Chem. Phys., 2018, 20, DOI: 10.1039/C7CP04561A

In this Reply we answer the two main arguments raised in the Comment. The first argument is related to the binding energy of the methanol dimer and its influence on the dimerization rate constant. We show that the dimerization rate constants calculated in the Comment are unphysically low. We report values that are about two orders of magnitude higher than the values of the Comment, which confirm the conclusions of the original article that dimers can be present in a small amount. The second argument based on the dependence of the pseudo-first order rates on the methanol concentration was already explained in detail in the Supporting Information of the original article.

Comment on “Methanol dimer formation drastically enhances hydrogen abstraction from methanol by OH at low temperature” by W. Siebrand, Z. Smedarchina, E. Martínez-Núñez and A. Fernández-Ramos, Phys. Chem. Chem. Phys., 2016, 18, 22712

A low fraction of methanol dimers is formed in the Laval nozzle used for kinetics studies at very low temperatures.

Benchmarking Quantum Chemical Methods: Are We Heading in the Right Direction?

Theoreticians and experimentalists should work together more closely to establish reliable rankings and benchmarks for quantum chemical methods. Comparison to carefully designed experimental benchmark data should be a priority. Guidelines to improve the situation for experiments and calculations are proposed.

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.

Spectroscopy and dynamics of double proton transfer in formic acid dimer

We present the isolated gas phase infrared spectra of formic acid dimer, (HCOOH)2, and its deuterated counterpart formic-d acid, (DCOOH)2, at room temperature. The formic acid dimer spectrum was obtained by spectral subtraction of a spectrum of formic acid vapor recorded at low pressure from that recorded at a higher pressure. The spectra of formic acid vapor contain features from both formic acid monomer and formic acid dimer, but at low and high pressures of formic acid, the equilibrium is pushed towards the monomer and dimer, respectively. A similar approach was used for the formic-d acid dimer. Building on the previous development of the Molecular Mechanics with Proton Transfer (MMPT) force field for simulating proton transfer reactions, molecular dynamics (MD) simulations were carried out to interpret the experimental spectra in the OH-stretching region. Within the framework of MMPT, a combination of symmetric single and double minimum potential energy surfaces (PESs) provides a good description of the double proton transfer PES. In a next step, potential morphing together with electronic structure calculations at the B3LYP and MP2 level of theory was used to align the computed and experimentally observed spectral features in the OH-stretching region. From this analysis, a barrier for double proton transfer between 5 and 7 kcal mol(-1) was derived, which compares with a CCSD(T)/aug-cc-pVTZ calculated barrier of 7.9 kcal mol(-1). Such a combination of experimental and computational techniques for estimating barriers for proton transfer in gas phase systems is generic and holds promise for further improved PESs and energetics of these important systems. Additional MD simulations at the semi-empirical DFTB level of theory agree quite well for the center band position but underestimate the width of the OH-stretching band.

A Symmetric Recognition Motif between Vicinal Diols: The Fourfold Grip in Ethylene Glycol Dimer

Ethylene glycol has a transiently chiral, asymmetric global minimum structure, but it favors a highly symmetric, achiral dimer arrangement which has not been considered or found in previous quantum‐chemical studies. Complementary FTIR and Raman spectroscopy in supersonic jets allows for the detection and straightforward assignment of this four‐fold hydrogen‐bonded dimer, which introduces an interesting supramolecular binding motif for vicinal diols and provides a strong case for transient chirality synchronization.

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.

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.