Proton tunnelling in the hydrogen bonds of halogen-substituted derivatives of benzoic acid studied by NMR relaxometry: the case of large energy asymmetry

1H/17O Chemical Shift Waves in Carboxyl-Bridged Hydrogen Bond Networks in Organic Solids

We report solid-state 1H and 17O NMR results for four 17O-labeled organic compounds each containing an extensive carboxyl-bridged hydrogen bond (CBHB) network in the crystal lattice: tetrabutylammonium hydrogen di-[17O2]salicylate (1), [17O4]quinolinic acid (2), [17O4]dinicotinic acid (3), and [17O2]Gly/[17O2]Gly·HCl cocrystal (4). The 1H isotropic chemical shifts found for protons involved in different CBHB networks are between 8.2 and 20.5 ppm, which reflect very different hydrogen-bonding environments. Similarly, the 17O isotropic chemical shifts found for the carboxylate oxygen atoms in CBHB networks, spanning a large range between 166 and 341 ppm, are also remarkably sensitive to the hydrogen-bonding environments. We introduced a simple graphical representation in which 1H and 17O chemical shifts are displayed along the H and O atomic chains that form the CBHB network. In such a depiction, because wavy patterns are often observed, we refer to these wavy patterns as 1H/17O chemical shift waves. Typical patterns of 1H/17O chemical shift waves in CBHB networks are discussed. The reported 1H and 17O NMR parameters for the CBHB network models examined in this study can serve as benchmarks to aid in spectral interpretation for CBHB networks in proteins.

Hydrogen Bonds: Raman Spectroscopic Study

The work outlines general ideas on how the frequency and the intensity of proton vibrations of X–H···Y hydrogen bonding are formed as the bond evolves from weak to maximally strong bonding. For this purpose, the Raman spectra of different chemical compounds with moderate, strong, and extremely strong hydrogen bonds were obtained in the temperature region of 5 K–300 K. The dependence of the proton vibrational frequency is schematically presented as a function of the rigidity of O-H···O bonding. The problems of proton dynamics on tautomeric O–H···O bonds are considered. A brief description of the N–H···O and C–H···Y hydrogen bonds is given.

Tautomeric Hydrogen Bond in Dimers of Ibuprofen

The Raman spectra of polycrystalline samples of ( RS)-2-(4-isobutylphenyl)-propionic acid of the common name ibuprofen have been measured in the temperature range 5-300 K. In the low-frequency spectrum of the normal C₁₂H₁₇(COOH) and deuterated C₁₂H₁₇(COOD) species, modes with ∼103 and ∼95 cm^-1 wavenumbers were detected, which corresponded to translational vibrations of O-H(D)···O hydrogen bonds of two different tautomers: left L and right R , respectively. At temperatures below 150 K, only the L-tautomer is found, and at T ≥ 150 K, both tautomers are observed. The energy difference Δ E of the ground vibrational state of potential minima for L- and R-tautomers is ∼80 meV for COOH and ∼70 meV for COOD. At T ≥ 150 K, the vibrational frequency of the C═O bond in the COOH moiety exhibits an unusual temperature dependence.

Unusual behavior of benzoic acid at low temperature: Raman spectroscopic study

The Raman spectra of benzoic acid single crystals have been measured in the temperature range of 5-300K. At T<60K the spectra show at least two anomalous features, one of which is of direct relevance to intensity changes of the lattice modes in the low-wavenumber region. The intensity of modes at ∼86 and ∼146cm(-1) tends to zero at T→0K. It is associated with appearance of two H-bonds of different length in the same l-tautomer, and with the loss of the inversion center in the dimer. The modes at ∼86 and ∼146cm(-1) are assigned to symmetric stretching intra-dimer vibrations of the OH⋯O hydrogen bonds of the first and second order, respectively. The assignment is based on the measurements of spectral parameters as function of temperature. The other anomaly is that the series of weak and narrow bands arises in the high-wavenumber region of 2500-3700cm(-1). The bands are assigned to combination tones of O-H hydrogen bonded stretching vibration and intramolecular modes. This effect results from a low-temperature transition of a conventional two wells potential of short H-bond in the l-tautomer to asymmetrical single well potential, and is due to a strong coupling of intramolecular vibrations to O-H stretching.

Proton tunnelling in the hydrogen bonds of the benzoic acid dimer: (18)O substitution and isotope effects of the heavy atom framework

Field-cycling (1)H NMR relaxometry has been used to measure the rate of concerted double proton transfer in the hydrogen bonds of (16)O and (18)O isotopologues of benzoic acid dimers. The experiments have been conducted in the solid state at low temperature 13.3 ≤ T ≤ 80 K where the dynamics are dominated by incoherent proton tunnelling. The low temperature tunnelling rate in the (16)O isotopologue is observed to be approximately 15% faster than in the (18)O isotopologue. The difference is attributed to an isotope effect of the heavy atom framework of the benzoic acid dimer resulting from displacements of the oxygen atoms that accompany the proton transfer. Sources of systematic uncertainty have been minimized in the design of the experimental protocols and the experiments are critically appraised in formally assigning the measured differences to an effect of mass on the tunnelling dynamics.

Complex methyl group and hydrogen-bonded proton motions in terms of the Arrhenius and Schrödinger equations

Equations for the temperature dependence of the spectral densities J(is)(m)(momega(I) +/-omega(T)), where m=1, 2, omega(I) and omega(T) are the resonance and tunnel splitting angular frequencies, in the presence of a complex motion, have been derived. The spin pairs of the protons or deuterons of the methyl group perform a complex motion consisting of three component motions. Two of them involve mass transportation over the barrier and through the barrier. They are characterized by k((H)) (Arrhenius) and k((T)) (Schrödinger) rate constants, respectively. The third motion causes fluctuations of the frequencies (nomega(I)+/-omega(T)) and it is related to the lifetime of the methyl spin at the energy level influenced by the rotor-bath interactions. These interactions induce rapid transitions, changing the symmetry of the torsional sublevels either from A to E or from E(a) to E(b). The correlation function for this third motion (k((omega)) rate constant) has been proposed by Müller-Warmuth et al. The spectral densities of the methyl group hindered rotation (k((H)), k((T)) and k((omega)) rate constants) differ from the spectral densities of the proton transfer (k((H)) and k((T)) rate constants) because three compound motions contribute to the complex motion of the methyl group. The recently derived equation [Formula: see text] , where [Formula: see text] and [Formula: see text] are the fraction and energy of particles with energies from zero to E(H), is taken into account in the calculations of the spectral densities. This equation follows from Maxwell's distribution of thermal energy. The spectral densities derived are applied to analyse the experimental temperature dependencies of proton and deuteron spin-lattice relaxation rate in solids containing the methyl group. A wide range of temperatures from zero Kelvin up to the melting point is considered. It has been established that the motion characterized by k((omega)) influences the spin-lattice relaxation up to the temperature T(tun) only. This temperature is directly determined by the equation C(p)T=E(H) (thermal energy=activation energy), where C(p) is the molar heat capacity. Probably the cessation of the third motion is a result of the de Broglie wavelength related to this motion becoming too short. As shown recently, the potential barrier can be an obstacle for the de Broglie wave. The theoretical equations derived in this paper are compared to those known in the literature.

Application of Schrödinger Equation to Study the Tunnelling Dynamics of Proton Transfer in the Hydrogen Bond of 2,5-Dinitrobenzoic Acid: Proton T1, T1ρ, and Deuteron T1 Relaxation Methods

Temperature measurements of proton T1 (24.7 MHz), deuteron (deuterated hydroxyl group) T1 (55.2 MHz), and proton T1(rho) (B1 = 9 G) spin-lattice relaxation times of 2,5-dinitrobenzoic acid have been performed. An analysis of present experimental data together with previously published proton T1 (55.2 MHz) data has revealed the following molecular motions: proton/deuteron transfer in the hydrogen bond and two-site hopping of the whole dimer. It is shown that the proton-transfer dynamics are characterized by two correlation times tau(ov) and tau(tu), describing two fundamentally different motional processes, namely, thermally activated jumps over the barrier and tunneling through the barrier. The temperature dependence of 1/tau(tu) is the solution of Schrödinger's equation, which also yields the temperature T(tun), where begins the tunnel pathway for proton transfer. A new equation for the spectral density function of complex motion consisting of the three motions is derived. The third motion (two-site hopping of the whole dimer characterized by tau(lib) correlation time) is responsible for a proton T1(rho) minimum in high temperatures, just below the melting point. Such a minimum is not reached by T1 temperature dependencies. The minimum of T1(rho) assigned to the classical hopping of a hydrogen-bonded proton occurs in the same low-temperature regime in which the flattening of the temperature dependencies of T1 points to the dominance of incoherent tunneling. This experimental fact denies the known theories predicting the intermediate temperature regime where a smooth transition between classical and quantum tunneling dynamics is expected. The fit of the derived theoretical equations to the experimental data T1(rho) and T1 is satisfactory. The correlation times obtained for deuterons indicate deuteron-transfer dynamics much slower than proton-transfer dynamics. It is concluded that the classical proton transfer takes place over the whole temperature regime, while the incoherent tunneling occurs below 46.5 (hydrogen) or 87.2 K (deuterium) only.

Isotope effects associated with tunneling and double proton transfer in the hydrogen bonds of benzoic acid

The isotope effects associated with double proton transfer in the hydrogen bonds of benzoic acid (BA) dimers have been measured using field-cycling (1)H NMR relaxometry and quasielastic neutron scattering. By studying mixed isotope (hydrogen and deuterium) samples, the dynamics of three isotopologues, BA-HH, BA-HD, and BA-DD, have been investigated. Low temperature measurements provide accurate measurements of the incoherent tunneling rate, k(0). This parameter scales accurately with the mass number, m, according to the formula k(0)=(E/m)e(-Fm) providing conclusive evidence that the proton transfer process is a strongly correlated motion of two hydrons. Furthermore, we conclude that the tunneling pathway is the same for the three isotopologue species. Measurements at higher temperatures illuminate the through barrier processes that are mediated via intermediate or excited vibrational states. In parallel with the investigation of proton transfer dynamics, the theoretical and experimental aspects of studying spin-lattice relaxation in single crystals of mixed isotope samples are investigated in depth. Heteronuclear dipolar interactions between (1)H and (2)H isotopes contribute significantly to the overall proton spin-lattice relaxation and it is shown that these must be modeled correctly to obtain accurate values for the proton transfer rates. Since the sample used in the NMR measurements was a single crystal, full account of the orientation dependence of the spin-lattice relaxation with respect to the applied B field was incorporated into the data analysis.

A 13C field-cycling NMR relaxometry investigation of proton tunnelling in the hydrogen bond: dynamic isotope effects, the influence of heteronuclear interactions and coupled relaxation

Concerted double proton transfer in the hydrogen bonds of a carboxylic acid dimer has been studied using 13C field-cycling NMR relaxometry. Heteronuclear 13C-1H dipolar interactions dominate the 13C spin-lattice relaxation which is significantly influenced by the polarisation state of the 1H Zeeman reservoir. The methodology of field-cycling experiments for such heteronuclear spin-coupled systems is studied experimentally and theoretically, including an investigation of various saturation-recovery and polarisation-recovery pulse sequence schemes. A theoretical model of the spin-lattice relaxation of this coupled system is presented which is corroborated by experiment. Spectral density components with frequencies omega(C), omega(C) + omega(H), and omega(C) - omega(H) are mapped out experimentally from the magnetic field dependence of the 13C and 1H spin-lattice relaxation and the proton transfer rate at low temperature is determined from their widths. Any dynamic isotope effect on the proton tunnelling in the hydrogen bond arising from 13C enrichment in the skeletal framework of the dimer is found to be smaller than experimental uncertainties (approximately 5%).

The ground-state tunneling splitting of various carboxylic acid dimers

Carboxylic acid dimers in gas phase reveal ground-state tunneling splittings due to a double proton transfer between the two subunits. In this study we apply a recently developed accurate semiclassical method to determine the ground-state tunneling splittings of eight different carboxylic acid derivative dimers (formic acid, benzoic acid, carbamic acid, fluoro formic acid, carbonic acid, glyoxylic acid, acrylic acid, and N,N-dimethyl carbamic acid) and their fully deuterated analogs. The calculated splittings range from 5.3e-4 to 0.13 cm(-1) (for the deuterated species from 2.8e-7 to 3.3e-4 cm(-1)), thus indicating a strong substituent dependence of the splitting, which varies by more than two orders of magnitude. One reason for differences in the splittings could be addressed to different barriers heights, which vary from 6.3 to 8.8 kcal/mol, due to different mesomeric stabilization of the various transition states. The calculated splittings were compared to available experimental data and good agreement was found. A correlation could be found between the tunneling splitting and the energy barrier of the double proton transfer, as the splitting increases with increased strength of the hydrogen bonds. From this correlation an empirical formula was derived, which allows the prediction of the ground-state tunneling splitting of carboxylic acid dimers at a very low cost and the tunneling splittings for parahalogen substituted benzoic acid dimers is predicted.