Symmetry of Hydrogen Bonds: Application of NMR Method of Isotopic Perturbation and Relevance of Solvatomers

Short, strong, symmetric, low-barrier hydrogen bonds (H-bonds) are thought to be of special significance. We have been searching for symmetric H-bonds by using the NMR technique of isotopic perturbation. Various dicarboxylate monoanions, aldehyde enols, diamines, enamines, acid-base complexes, and two sterically encumbered enols have been investigated. Among all of these, we have found only one example of a symmetric H-bond, in nitromalonamide enol, and all of the others are equilibrating mixtures of tautomers. The nearly universal lack of symmetry is attributed to the presence of these H-bonded species as a mixture of solvatomers, meaning isomers (or stereoisomers or tautomers) that differ in their solvation environment. The disorder of solvation renders the two donor atoms instantaneously inequivalent, whereupon the hydrogen attaches to the less well solvated donor. We therefore conclude that there is no special significance to short, strong, symmetric, low-barrier H-bonds. Moreover, they have no heightened stability or else they would have been more prevalent.

Isotopic-Perturbation NMR Study of Hydrogen-Bond Symmetry in Solution: Temperature Dependence and Comparison of OHO and ODO Hydrogen Bonds

Is a hydrogen bond symmetric, with the hydrogen centered between two donor atoms, or is it asymmetric, with the hydrogen closer to one but jumping to the other? The NMR method of isotopic perturbation has been used to distinguish these. Previous evidence from isotope shifts implies that a wide variety of dicarboxylate monanions are asymmetric, present as a rapidly equilibrating mixture of tautomers. However, calculations of hydrogen trajectories across an anharmonic potential-energy surface could reproduce the observed isotope shifts in a phthalate monoanion. Therefore, it was concluded that those isotope shifts are instead consistent with isotope-induced desymmetrization on a symmetric potential-energy surface. To distinguish between these two interpretations, the ¹⁸O-induced isotope effects on the ¹³C NMR chemical shifts of cyclohexene-1,2-dicarboxylate monoanion in chloroform-d and on the ¹⁹F NMR chemical shifts of difluoromaleate monoanion in D₂O have been investigated. In both cases the isotope effects are larger at lower temperature and also with deuterium in the hydrogen bond. It is concluded that these behaviors are consistent with the perturbation of an equilibrium between asymmetric tautomers and inconsistent with isotope-induced desymmetrization on a symmetric potential-energy surface.

Simplified calculation approaches designed to reproduce the geometry of hydrogen bonds in molecular complexes in aprotic solvents

The impact of the environment onto the geometry of hydrogen bonds can be critically important for the properties of the questioned molecular system. The paper reports on the design of calculation approaches capable to simulate the effect of aprotic polar solvents on the geometric and NMR parameters of intermolecular hydrogen bonds. A hydrogen fluoride and pyridine complex has been used as the main model system because the experimental estimates of these parameters are available for it. Specifically, F-H, F⋯N, and H-N distances, the values of ¹⁵N NMR shift, and spin-spin coupling constants ¹J(¹⁹F¹H), ^1hJ(¹H¹⁵N), and ^2hJ(¹⁹F¹⁵N) have been analyzed. Calculation approaches based on the gas-phase and the Polarizable Continuum Model (PCM) approximations and their combinations with geometric constraints and additional noncovalent interactions have been probed. The main result of this work is that the effect of an aprotic polar solvent on the geometry of a proton-donor⋯H⋯proton-acceptor complex cannot be reproduced under the PCM approximation if no correction for solvent-solute interactions is made. These interactions can be implicitly accounted for using a simple computational protocol.

NMR Study of Solvation Effect on the Geometry of Proton-Bound Homodimers of Increasing Size

Hydrogen bond geometries in the proton-bound homodimers of quinoline and acridine derivatives in an aprotic polar solution have been experimentally studied using ¹H NMR at 120 K. The reported results show that an increase of the dielectric permittivity of the medium results in contraction of the N···N distance. The degree of contraction depends on the homodimer's size and its substituent-specific solvation features. Neither of these effects can be reproduced using conventional implicit solvent models employed in computational studies. In general, the N···N distance in the homodimers of pyridine, quinoline, and acridine derivatives decreases in the sequence gas phase > solid state > polar solvent.

How short is the strongest hydrogen bond in the proton-bound homodimers of pyridine derivatives?

Hydrogen bond geometries in the proton-bound homodimers of ortho-unsubstituted and ortho-methylsubstituted pyridine derivatives in aprotic polar solution were estimated using experimental NMR data. Within the series of homodimers studied the hydrogen bond lengths depend on the proton affinity of pyridines and--at least for the ortho-methylsubstituted pyridines--on the pKa of the conjugate acids in an approximately quadratic manner. The shortest possible hydrogen bond in the homodimers of ortho-unsubstituted pyridines is characterized by the N···N distance of 2.613 Å. Steric repulsion between the methyl groups of the ortho-methylsubstituted pyridines becomes operative at an N···N distance of ∼2.7 Å and limits the closest approach to 2.665 Å.

Variable-temperature study of hydrogen-bond symmetry in cyclohexene-1,2-dicarboxylate monoanion in chloroform-d

The symmetry of the hydrogen bond in hydrogen cyclohexene-1,2-dicarboxylate monoanion was determined in chloroform using the NMR method of isotopic perturbation. As the temperature decreases, the (18)O-induced (13)C chemical-shift separations increase not only at carboxyl carbons but also at ipso (alkene) carbons. The magnitude of the ipso increase is consistent with an (18)O isotope effect on carboxylic acid acidity. Therefore it is concluded that this monoanion is a mixture of tautomers in rapid equilibrium, rather than a single symmetric structure in which a chemical-shift separation arises from coupling between a desymmetrizing vibration and anharmonic isotope-dependent vibrations, which is expected to show the opposite temperature dependence.

Hydrogen-bond symmetry in difluoromaleate monoanion

The symmetry of the hydrogen bond in hydrogen difluoromaleate monoanion is probed by X-ray crystallography and by the NMR method of isotopic perturbation in water, in two aprotic organic solvents, and in an isotropic liquid crystal. The X-ray crystal structure of potassium hydrogen difluoromaleate shows a remarkably short O-O distance of 2.41 Å and equal O-H distances of 1.206 Å, consistent with a strong and symmetric hydrogen bond. Incorporation of (18)O into one carboxyl group allows investigation of the symmetry of the H-bond in solution by the method of isotopic perturbation. The (19)F NMR spectra of the mono-(18)O-substituted monoanion in water, CD(2)Cl(2), and CD(3)CN show an AB spin system, corresponding to fluorines in different environments. The difference is attributed to the perturbation of the acidity of a carboxylic acid by (18)O, not to the mere presence of the (18)O, because the mono-(18)O dianion shows equivalent fluorines. Therefore, it is concluded that the monoanion exists as an equilibrating pair of interconverting tautomers and not as a single symmetric structure not only in water but also in organic solvents. However, in the isotropic liquid crystal phase of 4-cyanophenyl 4-heptylbenzoate, tetrabutylammonium hydrogen difluoromaleate-(18)O shows equivalent fluorines, consistent with a single symmetric structure. These results support earlier studies, which suggested that the symmetry of hydrogen bonds can be determined by the local environment.

On the Intramolecular Hydrogen Bond in Solution: Car-Parrinello and Path Integral Molecular Dynamics Perspective

The issue of the symmetry of short, low-barrier hydrogen bonds in solution is addressed here with advanced ab initio simulations of a hydrogen maleate anion in different environments, starting with the isolated anion, going through two crystal structures (sodium and potassium salts), then to an aqueous solution, and finally in the presence of counterions. By Car-Parrinello and path integral molecular dynamics simulations, it is demonstrated that the position of the proton in the intramolecular hydrogen bond of an aqueous hydrogen maleate anion is entirely related to the solvation pattern around the oxygen atoms of the intramolecular hydrogen bond. In particular, this anion has an asymmetric hydrogen bond, with the proton always located on the oxygen atom that is less solvated, owing to the instantaneous solvation environment. Simulations of water solutions of hydrogen maleate ion with two different counterions, K(+) and Na(+), surprisingly show that the intramolecular hydrogen-bond potential in the case of the Na(+) salt is always asymmetric, regardless of the hydrogen bonds to water, whereas for the K(+) salt, the potential for H motion depends on the location of the K(+). It is proposed that repulsion by the larger and more hydrated K(+) is weaker than that by Na(+) and competitive with solvation by water.

Geometry and spectral properties of the protonated homodimer of pyridine in the liquid and solid states. A combined NMR, X-ray diffraction and inelastic neutron scattering study

The structure and spectral signatures of the protonated homodimer of pyridine in its complex with a poorly coordinating anion have been studied in solution in CDF(3)/CDClF(2) down to 120 K and in a single crystal. In both phases, the hydrogen bond is asymmetric. In the solution, the proton is involved in a fast reversible transfer that determines the multiplicity of NMR signals and the sign of the primary H/D isotope effect of --0.95 ppm. The proton resonates at 21.73 ppm that is above any value reported in the past and is indicative of a very short hydrogen bond. By combining X-ray diffraction analysis with model computations, the position of the proton in the crystal has been defined as d(N-H) = 1.123 Å and d(H···N) = 1.532 Å. The same distances have been estimated using a (15)N NMR correlation. The frequency of the protonic out-of-plane bending mode is 822 cm(-1) in agreement with Novak's correlation.

Are short, low-barrier hydrogen bonds unusually strong?

In a symmetric hydrogen bond (H-bond), the hydrogen atom is perfectly centered between the two donor atoms. The energy diagram for hydrogen motion is thus a single-well potential, rather than the double-well potential of a more typical H-bond, in which the hydrogen is covalently bonded to one atom and H-bonded to the other. Examples of symmetric H-bonds are often found in crystal structures, and they exhibit the distinctive feature of unusually short length: for example, the O-O distance in symmetric OHO H-bonds is found to be less than 2.5 Å. In comparison, the O-O distance in a typical asymmetric H-bond, such as ROH···OR(2), ranges from about 2.7 to 3.0 Å. In this Account, we briefly review and update our use of the method of isotopic perturbation to search for a symmetric, centered, or single-well-potential H-bond in solution. Such low-barrier H-bonds are thought to be unusually strong, owing perhaps to the resonance stabilization of two identical resonance forms [A-H···B ↔ A···H-B]. This presumptive bond strength has been invoked to explain some enzyme-catalyzed reactions. Yet in solution, a wide variety of OHO, OHN, and NHN H-bonds have all been found to be asymmetric, in double-well potentials. Examples include the monoanion of (±)-2,3-di-tert-butylsuccinic acid and a protonated tetramethylnaphthalenediamine, even though these two ions are often considered prototypes of species with strong H-bonds. In fact, all of the purported examples of strong, symmetric H-bonds have been found to exist in solution as pairs of asymmetric tautomers, in contrast to their symmetry in some crystals. The asymmetry can be attributed to the disorder of the local solvation environment, which leads to an equilibrium among solvatomers (that is, isomers that differ in solvation). If the disorder of the local environment is sufficient to break symmetry, then symmetry itself is not sufficient to stabilize the H-bond, and symmetric H-bonds do not have an enhanced stability or an unusual strength. Nor are short H-bonds unusually strong. We discuss previous evidence for "short, strong, low-barrier" H-bonds and show it to be based on ambiguous comparisons. The role of such H-bonds in enzyme-catalyzed reactions is then ascribed not to any unusual strength of the H-bond itself but to relief of "strain."

Counteranion-dependent mechanisms of intramolecular proton transfer in aprotic solution

Using the freon mixture CDF(3)/CDClF(2) as solvent we have been able to measure the (1)H and (15)N NMR spectra of the doubly (15)N labeled 2,2'-bipyridinium cation (BpyH(+)) at temperatures down to 115 K. The obtained NMR parameters strongly depend on the type of counteranions indicating the formation of ion pairs. In the case of the bulky poorly coordinating tetrakis[3,5-bis(trifluoromethyl)phenyl]-borate as the counteranion a strong intramolecular NHN hydrogen bond was observed in BpyH(+) exhibiting a degenerate intramolecular proton transfer which is of the order of 10(6) s(-1) even at 120 K. By contrast, the weak hydrogen bond acceptor tetrafluoroborate favors a weak intermolecular FHN interaction and quenches the intramolecular proton transfer. The intramolecular proton transfer requires in this case a dissociation of the ion pair which is hindered by the Coulomb interaction. A slow intramolecular proton transfer was observed in the case of dichloroacetate which forms a strong intermolecular OHN hydrogen bond to BpyH(+). The mechanism of this transfer presumably involves a preliminary intermolecular proton transfer from nitrogen towards oxygen followed by a hydrogen bond switch to the neighboring nitrogen to which the proton is then transferred.