Vibrational tug-of-war: The pKAdependence of the broad vibrational features of strongly hydrogen-bonded carboxylic acids

Strong H-bonding from Zeolite Bro̷nsted Acid Site to Water: Origin of the Broad IR Doublet

Hydrogen bonding between water molecules and zeolite Bro̷nsted acid sites (BAS) has received much attention due to the significant influence of water on the adsorption and catalytic properties of these widely used porous materials. When a single water molecule is adsorbed at the BAS, the zeolite O-H stretch vibration decreases in frequency and splits into two extraordinarily broad bands peaked near 2500 and 2900 cm-1 in the infrared (IR) spectrum. This broad doublet feature is the predominant IR signature used to identify and interpret water-BAS H-bonding at low hydration levels, but the origin of the band splitting is not well understood. In this study, we used broadband two-dimensional infrared (2D IR) spectroscopy to investigate zeolite HZSM-5 prepared with a single water molecule per BAS. We find that the 2D IR spectrum is not explained by the most common interpretation of Fermi resonance coupling between the stretch and the bend of the BAS OH group, which predicts intense excited-state transitions that are absent from the experimental results. We present an alternative model of a double-well proton stretch potential, where the band splitting is caused by excited-state tunneling through the proton-transfer barrier. This one-dimensional model reproduces the basic experimental pattern of transition frequencies and amplitudes, suggesting that the doublet bands may originate from a highly anharmonic potential in which the excited state proton wave functions are delocalized over the H-bond between zeolite BAS and adsorbed H2O. Additional details about molecular orientation and coordination of the adsorbed water molecule are also resolved in the 2D IR spectroscopy.

Molecular Dynamics Simulation of Association Processes in Aqueous Solutions of Maleate Salts of Drug-like Compounds: The Role of Counterion

The study of the formation of microstructures during the interaction of a protonated drug-like compound (API) with a maleic acid monoanion sheds light on the assembly processes in an aqueous solution at the molecular level. Molecular dynamics (MD) simulations coupled with density functional theory (DFT) calculations made it possible to find initial hydrogen bonding motifs during the assembly process, leading to the formation of heterodimers and trimers. The process of trimer formation [protonated API—maleic acid monoanion—protonated API] proceeds through the formation of three intermolecular H-bonds by the CO₂⁻ group of the maleic acid monoanion in both systems. The total enthalpy/energy of these H-bonds is more than 70 kJ/mol. Thus, the maleic acid monoanion plays a key role in the processes of association in aqueous solution, and the interaction of the maleic acid monoanion with API is more preferable than the interaction of API molecules with each other. DFT computations in the discrete continuum approximation reveal the spectral features of heterodimers and trimers, and the ATR-IR spectra confirmed these findings. MD simulations followed by DFT calculations made it possible to describe the initial stages of the formation of pharmaceutical cocrystals in an aqueous solution.

Vibrational Stark fields in carboxylic acid dimers

Carboxylic acids form exceptionally stable dimers and have been used to model proton and double proton transfer processes. The stabilization energies of the carboxylic acid dimers are very weakly dependent on the nature of substitution. However, the electric field experienced by the OH group of a particular carboxylic acid is dependent more on the nature of the substitution on the dimer partner. In general, the electric field was higher when the partner was substituted with an electron-donating group and lower with an electron-withdrawing substituent on the partner. The Stark tuning rate (Δ) of the O-H stretching vibrations calculated at the MP2/aug-cc-pVDZ level was found to be weakly dependent on the nature of substitution on the carboxylic acid. The average Stark tuning rate of O-H stretching vibrations of a particular carboxylic acid when paired with other acids was 5.7 cm^-1 (MV cm^-1)^-1, while the corresponding average Stark tuning rate of the partner acids due to a particular carboxylic acid was 21.9 cm^-1 (MV cm^-1)^-1. The difference in the Stark tuning rate is attributed to the primary and secondary effects of substitution on the carboxylic acid. The average Stark tuning rate for the anharmonic O-D frequency shifts is about 40-50% higher than the corresponding harmonic O-D frequency shifts calculated at the B3LYP/aug-cc-pVDZ level, much greater than the typical scaling factors used, indicating the strong anharmonicity of O-H/O-D oscillators in carboxylic acid dimers. Finally, the linear correlation observed between pK_a and the electric field was used to estimate the pK_a of fluoroformic acid to be around 0.9.

Quantum mechanical effects in acid–base chemistry

Acid–base chemistry has immense importance for explaining and predicting the chemical products formed by an acid and a base when mixed together. However, the traditional chemistry theories used to describe acid–base reactions do not take into account the effect arising from the quantum mechanical nature of the acidic hydrogen shuttling potential and its dependence on the acid base distance. Here, infrared and NMR spectroscopies, in combination with first principles simulations, are performed to demonstrate that quantum mechanical effects, including electronic and nuclear quantum effects, play an essential role in defining the acid–base chemistry when 1-methylimidazole and acetic acid are mixed together. In particular, it is observed that the acid and the base interact to form a complex containing a strong hydrogen bond, in which the acidic hydrogen atom is neither close to the acid nor to the base, but delocalized between them. In addition, the delocalization of the acidic hydrogen atom in the complex leads to characteristic IR and NMR signatures. The presence of a hydrogen delocalized state in this simple system challenges the conventional knowledge of acid–base chemistry and opens up new avenues for designing materials in which specific properties produced by the hydrogen delocalized state can be harvested.

Acid-based theories do not consider the quantum mechanical nature of the acidic hydrogen shuttling potential. Here, it is demonstrated that this particularity is needed to explain the formation acid-base complex with a delocalized acidic hydrogen.

Probing Hydrogen-Bonding Interactions within Phenol-Benzimidazole Proton-Coupled Electron Transfer Model Complexes with Cryogenic Ion Vibrational Spectroscopy

Hydrogen-bonding interactions within a series of phenol-benzimidazole model proton-coupled electron transfer (PCET) dyad complexes are characterized using cryogenic ion vibrational spectroscopy. A highly red-shifted and surprisingly broad (>1000 cm^-1) transition is observed in one of the models and assigned to the phenolic OH stretch strongly H-bonded to the N₍₃₎ benzimidazole atom. The breadth is attributed to a combination of anharmonic Fermi-resonance coupling between the OH stretch and background doorway states involving OH bending modes and strong coupling of the OH stretch frequency to structural deformations along the proton-transfer coordinate accessible at the vibrational zero-point level. The other models show unexpected protonation of the benzimidazole group upon electrospray ionization instead of at more basic remote amine/amide groups. This leads to the formation of HO-⁺HN₍₃₎ H-bond motifs that are much weaker than the OH-N₍₃₎ H-bond arrangement. H-bonding between the N₍₁₎H⁺ benzimidazole group and the carbonyl on the tyrosine backbone is the stronger and preferred interaction in these complexes. The results show that conjugation effects, secondary H-bond interactions, and H-bond soft modes strongly influence the OH-N₍₃₎ interaction and highlight the importance of the direct monitoring of proton stretch transitions in characterizing the proton-transfer reaction coordinate in PCET systems.

Composition-Dependent Hydrogen-Bonding Motifs and Dynamics in Brønsted Acid-Base Mixtures

In recent years the interaction of organophosphates and imines, which is at the core of Brønsted acid organocatalysis, has been established to be based on strong ionic hydrogen bonds. Yet, besides the formation of homodimers consisting of two acid molecules and heterodimers consisting of one acid and one base, also multimeric molecular aggregates are formed in solution. These multimeric aggregates consist of one base and several acid molecules. The details of the intermolecular bonding in such aggregates, however, have remained elusive. To characterize composition-dependent bonding and bonding dynamics in these aggregates, we use linear and nonlinear infrared (IR) spectroscopy at varying molar ratios of diphenyl phosphoric acid and quinaldine. We identify the individual aggregate species, giving rise to the structured, strong, and very broad infrared absorptions, which span more than 1000 cm^–1. Linear infrared spectra and density functional theory calculations of the proton transfer potential show that doubly ionic intermolecular hydrogen bonds between the acid and the base lead to absorptions which peak at ∼2040 cm^–1. The contribution of singly ionic hydrogen bonds between an acid anion and an acid molecule is observed at higher frequencies. As common to such strong hydrogen bonds, ultrafast IR spectroscopy reveals rapid, ∼ 100 fs, dissipation of energy from the proton transfer coordinate. Yet, the full dissipation of the excess energy occurs on a ∼0.8–1.1 ps time scale, which becomes longer when multimers dominate. Our results thus demonstrate the coupling and collectivity of the hydrogen bonds within these complexes, which enable efficient energy transfer.

IR Spectra of Hydrogen-Bonded Complexes of Trifluoroacetic Acid with Acetone and Diethyl Ether in the Gas Phase. Interaction between CH and OH Stretching Vibrations

The spectra of complexes of trifluoroacetic acid (TFA) with acetone and diethyl ether (DEE) and their perdeuterated isotopologues were extracted from the spectra of the mixture of the compounds recorded at room temperature. The ν(OH) bands of the complexes with protiated and deuterated acetone notably differ from each other, whereas these ν(OH) bands are practically not affected by the deuteration of DEE. An assumption about the interaction of CH and OH groups in the (CH₃)-C═O···HO fragment is made. According to density functional theory calculations, complexes of TFA with both acetone and DEE have a cyclic structure with one strong ═O···HO hydrogen bond and one weak CH···O═ bond. The structural, spectroscopic, and electronic properties indicate an essential role of weak bonds in the total complexation energy of the systems studied.

Calculation of vibrational spectroscopic and geometrical characteristics of the [F(HF)2]- and [F(DF)2]- complexes using the second-order vibrational perturbation theory and a 6D variational method

Vibrational spectroscopic and average geometrical parameters of the strong H-bonded complexes [F(HF)₂]^- and [F(DF)₂]^- are determined for the first time from nine-dimensional (9D) perturbative and 6D variational calculations. The frequencies and intensities for all fundamental and some combination and overtone transitions obtained by the method of second-order vibrational perturbation theory (VPT2) are reported. A two-fold decrease in the H-F (D-F) stretching band frequency and a more than ten-fold increase in the intensity of this band upon complexation are predicted. The theoretical frequencies for both isolated isotopologues are in satisfactory agreement (to better than 70 cm^-1) with the scarce experimental data obtained in condensed phases. The main purpose of variational calculations is to analyze the intermode anharmonic coupling and the changes in the geometrical parameters upon vibrational excitation and H/D isotopic substitution. The equilibrium nuclear configuration and the 2D potential energy surface (PES) of [F(HF)₂]^- for H-F stretches are calculated in the MP2/6-311++G(3df,3pd), CCSD(T)/6-311++G(3df,3pd), CCSD(T)/aug-cc-pVTZ, and CCSD(T)/d-aug-cc-pVTZ approximations with the basis set superposition error taken into account. Anharmonic vibrational problems are solved by the variational method for 2D, 4D, and 6D systems of H-bond and H-F (D-F) stretches and in-plane bends. The VPT2 calculations and calculations of the PESs for 4D and 6D systems are performed in the MP2/6-311++G(3df,3pd) approximation. Comparison of variational anharmonic solutions for different vibrational subsystems demonstrates the influence of intermode anharmonic coupling on the mixing of wave functions and spectroscopic and geometrical characteristics. The inverse Ubbelohde effect is predicted and substantiated.