Isotopic Quantum Effects in Liquid Methanol

Thermodynamically Stable Cationic Dimers in Carboxyl-Functionalized Ionic Liquids: The Paradoxical Case of "Anti-Electrostatic" Hydrogen Bonding

We show that carboxyl-functionalized ionic liquids (ILs) form doubly hydrogen-bonded cationic dimers (c⁺=c⁺) despite the repulsive forces between ions of like charge and competing hydrogen bonds between cation and anion (c⁺–a⁻). This structural motif as known for formic acid, the archetype of double hydrogen bridges, is present in the solid state of the IL 1−(carboxymethyl)pyridinium bis(trifluoromethylsulfonyl)imide [HOOC−CH₂−py][NTf₂]. By means of quantum chemical calculations, we explored different hydrogen-bonded isomers of neutral (HOOC–(CH₂)_n–py⁺)₂(NTf₂⁻)₂, single-charged (HOOC–(CH₂)_n–py⁺)₂(NTf₂⁻), and double-charged (HOOC– (CH₂)_n−py⁺)₂ complexes for demonstrating the paradoxical case of “anti-electrostatic” hydrogen bonding (AEHB) between ions of like charge. For the pure doubly hydrogen-bonded cationic dimers (HOOC– (CH₂)_n−py⁺)₂, we report robust kinetic stability for n = 1–4. At n = 5, hydrogen bonding and dispersion fully compensate for the repulsive Coulomb forces between the cations, allowing for the quantification of the two equivalent hydrogen bonds and dispersion interaction in the order of 58.5 and 11 kJmol⁻¹, respectively. For n = 6–8, we calculated negative free energies for temperatures below 47, 80, and 114 K, respectively. Quantum cluster equilibrium (QCE) theory predicts the equilibria between cationic monomers and dimers by considering the intermolecular interaction between the species, leading to thermodynamic stability at even higher temperatures. We rationalize the H-bond characteristics of the cationic dimers by the natural bond orbital (NBO) approach, emphasizing the strong correlation between NBO-based and spectroscopic descriptors, such as NMR chemical shifts and vibrational frequencies.

Size of the hydrogen bond network in liquid methanol: a quantum cluster equilibrium model with extensive structure search

Studies have debated what is a favorable cluster size in liquid methanol. Applications of the quantum cluster equilibrium (QCE) model on a limited set of cluster structures have demonstrated the dominance of cyclic hexamers in liquid methanol. In this study, we examined the aforementioned question by integrating our implementation of QCE with a molecular-dynamics-based structural searching scheme. QCE simulations were performed using a database comprising extensively searched stable conformers of (MeOH)n for n = 2-14, which were optimized by B3LYP/6-31+G(d,p) with and without the dispersion correction. Our analysis indicated that an octamer structure can contribute significantly to cluster probability. By reoptimizing selected conformers with high probability at the MP2 level, we found that the aforementioned octamer became the dominant species due to favorable vibrational free energy, which was attributed to modes of intermolecular vibration.

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.

Competition between hydrogen bonds and van der Waals forces in intermolecular structure formation of protonated branched-chain alcohol clusters

To investigate the influence of bulky alkyl groups on hydrogen-bonded (H-bonded) network structures of alcohols, infrared (IR) spectra of protonated clusters of 2-propanol (2-PrOH) and tert-butyl alcohol (t-BuOH) were observed in the OH and CH stretch regions. In addition, by varying the tag species, the temperature dependence profile of the isomer population of H+(t-BuOH)n was revealed. An extensive search for stable isomers was performed using dispersion-corrected density functional theory methods, and temperature-dependent IR spectral simulations were done on the basis of the harmonic superposition approximation. The computational results qualitatively agreed with the observed size and temperature dependence of the H-bonded network structures of these protonated bulky alcohol clusters. However, the difficulty in the quantitative evaluation of dispersion was also demonstrated. It was shown that H+(2-PrOH)n (n = 4-7) have essentially the same network structures as the protonated normal alcohol clusters studied so far. On the other hand, H+(t-BuOH)n (n = 4-8) showed a clear preference for the smaller-membered ring structures, that is very different from the preference of the protonated normal alcohol clusters. The origin of the different structure preferences was discussed in terms of the steric effect and dispersion.

Microheterogeneity in CH3OH/CD3OH mixture

Recently, we demonstrated the presence of microheterogeneity in binary mixtures of unlike alcohols. [RSC Adv. 2016, 6, 37195-37202] The aim of this work was examination if this phenomenon occurs also in the mixture of very similar alcohols like CH₃OH and CD₃OH. Theoretical calculations suggest that the isotopic substitution in methyl group influences properties of the OH group. Hence, one can expect that this effect may lead to partial separation of CH₃OH and CD₃OH at a molecular level and it contributes to deviation from the ideal mixture. This work evidences that CH₃OH/CD₃OH mixture also deviates from the ideal one, but the extent of this deviation is much smaller as compared with the mixtures of other alcohols. It is of particular note that this deviation results mainly from the difference between the CH₃ and CD₃ groups, while the contribution from the OH groups is small. The structure of CH₃OH/CD₃OH mixture at a molecular level is similar to the structure of binary mixtures of other alcohols. The mixture is composed of the homoclusters of both alcohols and the mixed clusters. The homoclusters existing in the mixture are similar to those present in bulk alcohols. The highest population of the heteroclusters and the largest deviation from the ideal mixture were observed at equimolar mixture. Both the experimental and theoretical results reveal that in CH₃OH/CD₃OH mixture dominate the cyclic tetramers and larger clusters, while the population of the linear clusters is negligible. Though the extent and strength of hydrogen bonding in both alcohols are the same, the position and intensity of the 2ν(OH) band for CH₃OH and CD₃OH are different. We propose possible explanation of this observation.

Hydrogen bond network structures of protonated short-chain alcohol clusters

Protonated alcohol clusters enable extraction of the physical essence of the nature of hydrogen bond networks.

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.

Anharmonic effects in the quantum cluster equilibrium method

The well-established quantum cluster equilibrium (QCE) model provides a statistical thermodynamic framework to apply high-level ab initio calculations of finite cluster structures to macroscopic liquid phases using the partition function. So far, the harmonic approximation has been applied throughout the calculations. In this article, we apply an important correction in the evaluation of the one-particle partition function and account for anharmonicity. Therefore, we implemented an analytical approximation to the Morse partition function and the derivatives of its logarithm with respect to temperature, which are required for the evaluation of thermodynamic quantities. This anharmonic QCE approach has been applied to liquid hydrogen chloride and cluster distributions, and the molar volume, the volumetric thermal expansion coefficient, and the isobaric heat capacity have been calculated. An improved description for all properties is observed if anharmonic effects are considered.

Quantum Chemical Methods for the Prediction of Energetic, Physical, and Spectroscopic Properties of Ionic Liquids

The accurate prediction of physicochemical properties of condensed systems is a longstanding goal of theoretical (quantum) chemistry. Ionic liquids comprising entirely of ions provide a unique challenge in this respect due to the diverse chemical nature of available ions and the complex interplay of intermolecular interactions among them, thus resulting in the wide variability of physicochemical properties, such as thermodynamic, transport, and spectroscopic properties. It is well understood that intermolecular forces are directly linked to physicochemical properties of condensed systems, and therefore, an understanding of this relationship would greatly aid in the design and synthesis of functionalized materials with tailored properties for an application at hand. This review aims to give an overview of how electronic structure properties obtained from quantum chemical methods such as interaction/binding energy and its fundamental components, dipole moment, polarizability, and orbital energies, can help shed light on the energetic, physical, and spectroscopic properties of semi-Coulomb systems such as ionic liquids. Particular emphasis is given to the prediction of their thermodynamic, transport, spectroscopic, and solubilizing properties.

Cation-cation clusters in ionic liquids: Cooperative hydrogen bonding overcomes like-charge repulsion

Direct spectroscopic evidence for H-bonding between like-charged ions is reported for the ionic liquid, 1-(2-hydroxyethyl)-3-methylimidazolium tetrafluoroborate. New infrared bands in the OH frequency range appear at low temperatures indicating the formation of H-bonded cation-cation clusters similar to those known for water and alcohols. Supported by DFT calculations, these vibrational bands can be assigned to attractive interaction between the hydroxyl groups of the cations. The repulsive Coulomb interaction is overcome by cooperative hydrogen bonding between ions of like charge. The transition energy from purely cation-anion interacting configurations to those including cation-cation H-bonds is determined to be 3-4 kJmol(-1). The experimental findings and DFT calculations strongly support the concept of anti-electrostatic hydrogen bonds (AEHBs) as recently suggested by Weinhold and Klein. The like-charge configurations are kinetically stabilized with decreasing temperatures.

Structural properties of methanol–water binary mixtures within the quantum cluster equilibrium model

The Quantum Cluster Equilibrium (QCE) method computes cluster distributions and thermodynamic properties of binary methanol–water mixtures in agreement with experiments.

The microscopic structure of liquid methanol from Raman spectroscopy

The microscopic structure of liquid methanol has been systematically investigated with Raman spectroscopy in a temperature range of 15-55 degrees C. The unbonded free -OH stretching vibrational band has been observed at approximately 3660 cm(-1) in pure liquid. With the aid of depolarization measurements and theoretical calculations, four featured spectral components have been unambiguously identified and assigned to four well-defined vibrational modes of clusters in chain or ring forms. Furthermore, the cluster size distribution and its temperature dependence have been derived from the spectral fittings for the first time, which lead to the conclusion that the trimer, tetramer, and pentamer are the dominant clusters in liquid methanol, taking up more than 50% of total clusters.

Cooperative versus dispersion effects: What is more important in an associated liquid such as water?

We implemented the quantum cluster equilibrium theory in our postprocessing program PEACEMAKER. This program may be run in conjunction with the very efficient vibrational frequency analysis code SNF and can therefore provide access to all electronic structure programs combined with this program. We applied the quantum cluster equilibrium theory in order to investigate the influence of a wide range of electronic structure models on the description of the liquid state. This investigation revealed much about the relevance of approximations in modern simulations of associated liquids such as water. While it is often claimed that the use of density-functional theory in condensed matter is leading to gravely erroneous results, we found that, contrary to these assertions, the exact exchange functional B3LYP and the gradient-corrected functional BP perform very well in combination with sizable basis sets as compared to second-order Moller-Plesset perturbation theory employing the same basis set. The use of density-functional theory with smaller basis sets does, in fact, lead to better results in the liquid state than the use of second-order Moller-Plesset perturbation theory in combination with these small basis sets. Most importantly, the neglect of cooperative effects disturbs a good description much more evenly if we apply second-order Moller-Plesset perturbation theory in combination with large basis sets than density-functional theory including cooperativity with smaller basis sets or Hartree-Fock using a very small basis set.