Accurate Calculations of OH-Stretching Intensities with a Reduced-Dimensional Local Mode Model Including Eckart Axis Embedding

Effect of Temperature on the OH-Stretching Bands of the Methanol Dimer

We present a conceptually simple model for understanding the significant spectral changes that occur with the temperature in the infrared spectra of hydrogen-bound complexes. We have measured room-temperature spectra of the methanol dimer and two deuterated isotopologues in the OH(D)-stretching region. We correctly predict spectral changes observed in the gas phase for the bound OH stretch in the methanol dimer from jet-cooled to room temperature and corroborate this with experimental and theoretical results for deuterated isotopologues. The origin of the observed spectral features is explained based on a reduced-dimensional vibrational model, which includes the two high-frequency OH stretches, the two methyl torsions, and the six intermolecular low-frequency vibrations. Key to the success of the model is a new coordinate definition to describe the intrinsic large-amplitude curvilinear motion of low-frequency vibrations. Despite the deceivingly simple appearance of the room temperature bound OH-stretching fundamental band, it consists of ∼107 vibrational transitions.

Quantum Chemical Investigation of the Cold Water Dimer Spectrum in the First OH-Stretching Overtone Region Provides a New Interpretation

Intramolecular vibrational transition wavenumbers and intensities were calculated in the fundamental HOH-bending, fundamental OH-stretching, first OH-stretching-HOH-bending combination, and first OH-stretching overtone (ΔvOH = 2) regions of the water dimer's spectrum. Furthermore, the rotational-vibrational spectrum was calculated in the ΔvOH = 2 region at 10 K, corresponding to the temperature of the existing jet-expansion experiments. The calculated spectrum was obtained by combining results from a full-dimensional (12D) vibrational and a reduced-dimensional vibrational-rotational-tunneling model. The ΔvOH = 2 spectral region is rich in features due to contributions from multiple vibrational-rotational-tunneling sub-bands. Origins of the experimental vibrational bands depend on the assignment of the observed sub-bands. Based on our calculations, we assign the observed sub-bands, and our reassignment leads to new values for the vibrational band origins of the free donor and antisymmetric acceptor OH-stretching first overtones of ∼7227 and ∼7238 cm-1, respectively. The observed bands with origins at 7192.34 and ∼7366 cm-1 are assigned to the symmetric acceptor OH-stretching first overtone and the OH-stretching combination of the donor, respectively.

Setting up the HyDRA blind challenge for the microhydration of organic molecules

The procedure leading to the first HyDRA blind challenge for the prediction of water donor stretching vibrations in monohydrates of organic molecules is described. A training set of 10 monohydrates...

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.

Estimations of FH···X hydrogen bond energies from IR intensities: Iogansen's rule revisited

In this work the possibility of using the IR intensity of the stretching vibration ν_s of proton donor group for estimation of hydrogen bond strength was investigated. For a set of complexes with FH···X (X = F, N, O) hydrogen bonds in the wide range of energies (0.1-49.2 kcal/mol) vibrational frequencies ν_s and their intensities A were calculated (CCSD at complete basis set limit). The validity of the previously proposed linear proportionality between the intensification of the stretching vibration ν_s in IR spectra and hydrogen bond enthalpy -ΔH = 12.2 ∆ A (A. V. Iogansen, Spectrochimica Acta A 1999) was examined. It is shown that for a range of similar hydrogen bond types with complexation energies ∆E <15 kcal/mol the ∆E( ∆ A ) function remains similar to that proposed in the Iogansen's work, while upon strengthening this dependency becomes significantly nonlinear. We examined two other parameters ( ∆ A ν s and ∆ A ∙ m R ) related to IR intensity as descriptors of hydrogen bond strength which are proportional to transition dipole moment matrix element and mass-independent dipole moment derivative. It was found that the dependency ∆E( ∆ A ν s ) stays linear in the whole studied range of complexation energies and it can be used for evaluation of ∆E from infrared spectral data with the accuracy about 2 kcal/mol. The mass-independent product ∆ A ∙ m R is an appropriate descriptor for sets of complexes with various hydrogen bond types. Simple equations proposed in this work can be used for estimations of hydrogen bond strength in various systems, where experimental thermodynamic methods or direct calculations are difficult or even impossible.

Graphical Transition Moment Decomposition and Conceptual Density Functional Theory Approaches to Study the Fundamental and Lower-Level Overtone Absorption Intensities of Some OH Stretching Vibrations

The investigation of electron density migrations caused by molecular structure changes is of central importance in various fields of chemistry. To address this topic in general and to study absorption intensities of vibrations, we analyze sensitive dipole moment functions (DMFs) of a molecule by combining the linear response function of conceptual DFT and bond dipoles separated by the quantum theory of atoms in molecules with a graphical transition moment decomposition scheme. The fundamental intensities of OH stretching vibrations depend strongly on the substituents but only weakly on the molecular conformations. Interestingly, in some alcohols, completely opposite trends have been observed for the lower-level overtone intensities: a weak substituent dependence but a stronger conformation dependence. It is well known that the formation of a hydrogen-bonded complex increases the OH stretching fundamental intensity, but less well known is the decrease in their overtone intensities. To investigate these characteristics comprehensively, we calculated their intensities (Δv = 1, 2, and 3) for conformers of ethanol and trifluoroethanol (TFE) and hydrogen-bonded phenol (PhOH) systems via the DFT method in the local mode model for the OH stretching coordinate ΔR. Their first and second derivatives of the electron density with respect to ΔR were calculated and interpreted using their bond moments. For ethanol and TFE, the OH, CC, and CH bond moments were found to make an important contribution to the molecular DMF derivatives parallel to the OH bond. The OH bond contributes only to the first derivative of DMF, and its conformational dependence is determined by the magnitude of the charge polarization of each structure. The electron density derivatives in the CC bond region were largely maintained during the internal rotation; thus, their conformation-dependent contributions were expressed by a geometrical factor of the CC bond direction. The CH bond at the antiperiplanar position of the OH bond was found to make a remarkably large contribution to the second derivative of DMF in the gauche conformer. The importance of electron density migration on substituents was also identified in the hydrogen-bonded phenol, in which the π-electron density change on the aromatic ring was clearly shown. This migration creates the DMF derivatives both perpendicular and parallel to the OH bond and strongly affects the absorption intensities. In all the cases, some bond moments on the substituents contribute to the first and second DMF derivatives in a structure-dependent manner, thus explaining their stereoelectronic effects.

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