Interpretation of Hydrogen/Deuterium Isotope Effects on NMR Chemical Shifts of [FHF]− Ion Based on Calculations of Nuclear Magnetic Shielding Tensor Surface

Weak, Broken, but Working-Intramolecular Hydrogen Bond in 2,2'-bipyridine

From an academic and practical point of view, it is desirable to be able to assess the possibility of the proton exchange of a given molecular system just by knowing the positions of the proton acceptor and the proton donor. This study addresses the difference between intramolecular hydrogen bonds in 2,2'-bipyridinium and 1,10-phenanthrolinium. Solid-state ¹⁵N NMR and model calculations show that these hydrogen bonds are weak; their energies are 25 kJ/mol and 15 kJ/mol, respectively. Neither these hydrogen bonds nor N-H stretches can be responsible for the fast reversible proton transfer observed for 2,2'-bipyridinium in a polar solvent down to 115 K. This process must have been caused by an external force, which was a fluctuating electric field present in the solution. However, these hydrogen bonds are the grain that tips the scales precisely because they are an integral part of a large system of interactions, including both intramolecular interactions and environmental influence.

Modeling of the Response of Hydrogen Bond Properties on an External Electric Field: Geometry, NMR Chemical Shift, Spin-Spin Scalar Coupling

The response of the geometric and NMR properties of molecular systems to an external electric field has been studied theoretically in a wide field range. It has been shown that this adduct under field approach can be used to model the geometric and spectral changes experienced by molecular systems in polar media if the system in question has one and only one bond, the polarizability of which significantly exceeds the polarizability of other bonds. If this requirement is met, then it becomes possible to model even extreme cases, for example, proton dissociation in hydrogen halides. This requirement is fulfilled for many complexes with one hydrogen bond. For such complexes, this approach can be used to facilitate a detailed analysis of spectral changes associated with geometric changes in the hydrogen bond. For example, in hydrogen-bonded complexes of isocyanide C≡¹⁵N-¹H⋯X, ¹J(¹⁵N¹H) depends exclusively on the N-H distance, while δ(¹⁵N) is also slightly influenced by the nature of X.

NMR Properties of the Cyanide Anion, a Quasisymmetric Two-Faced Hydrogen Bonding Acceptor

The isotopically enriched cyanide anion, (13C≡15N)−, has a great potential as the NMR probe of non-covalent interactions. However, hydrogen cyanide is highly toxic and can decompose explosively. It is therefore desirable to be able to theoretically estimate any valuable results of certain experiments in advance in order to carry out experimental studies only for the most suitable molecular systems. We report the effect of hydrogen bonding on NMR properties of 15N≡13CH···X and 13C≡15NH···X hydrogen bonding complexes in solution, where X = 19F, 15N, and O=31P, calculated at the ωB97XD/def2tzvp and the polarizable continuum model (PCM) approximations. In many cases, the isotropic 13C and 15N chemical shieldings of the cyanide anion are not the most informative NMR properties of such complexes. Instead, the anisotropy of these chemical shieldings and the values of scalar coupling constants, including those across hydrogen bonds, can be used to characterize the geometry of such complexes in solids and solutions. 1J(15N13C) strongly correlates with the length of the N≡C bond.

Actual Symmetry of Symmetric Molecular Adducts in the Gas Phase, Solution and in the Solid State

This review discusses molecular adducts, whose composition allows a symmetric structure. Such adducts are popular model systems, as they are useful for analyzing the effect of structure on the property selected for study since they allow one to reduce the number of parameters. The main objectives of this discussion are to evaluate the influence of the surroundings on the symmetry of these adducts, steric hindrances within the adducts, competition between different noncovalent interactions responsible for stabilizing the adducts, and experimental methods that can be used to study the symmetry at different time scales. This review considers the following central binding units: hydrogen (proton), halogen (anion), metal (cation), water (hydrogen peroxide).

Modeling of Solute-Solvent Interactions Using an External Electric Field-From Tautomeric Equilibrium in Nonpolar Solvents to the Dissociation of Alkali Metal Halides

An implicit account of the solvent effect can be carried out using traditional static quantum chemistry calculations by applying an external electric field to the studied molecular system. This approach allows one to distinguish between the effects of the macroscopic reaction field of the solvent and specific solute-solvent interactions. In this study, we report on the dependence of the simulation results on the use of the polarizable continuum approximation and on the importance of the solvent effect in nonpolar solvents. The latter was demonstrated using experimental data on tautomeric equilibria between the pyridone and hydroxypyridine forms of 2,6-di-tert-butyl-4-hydroxy-pyridine in cyclohexane and chloroform.

Editorial to the Special Issue "Gulliver in the Country of Lilliput: An Interplay of Noncovalent Interactions"

Noncovalent interactions allow our world to exist [...].

Electric field effect on 31P NMR magnetic shielding

Magnetic shielding depends on molecular structure and noncovalent interactions. This study shows that it is also measurably dependent on the electric field generated by surrounding molecules. This effect has been observed explicitly for ³¹P nucleus using the adduct under field approach. The results obtained indicate that the field strength experienced by molecules in crystals consisting of molecules with large dipole moments is similar to that in polar solvents. Therefore, magnetic shielding should explicitly depend on solvent polarity. It is important to note that this effect cannot be reproduced correctly within the polarizable continuum model approach.

Inter- vs. Intramolecular Hydrogen Bond Patterns and Proton Dynamics in Nitrophthalic Acid Associates

Noncovalent interactions are among the main tools of molecular engineering. Rational molecular design requires knowledge about a result of interplay between given structural moieties within a given phase state. We herein report a study of intra- and intermolecular interactions of 3-nitrophthalic and 4-nitrophthalic acids in the gas, liquid, and solid phases. A combination of the Infrared, Raman, Nuclear Magnetic Resonance, and Incoherent Inelastic Neutron Scattering spectroscopies and the Car-Parrinello Molecular Dynamics and Density Functional Theory calculations was used. This integrated approach made it possible to assess the balance of repulsive and attractive intramolecular interactions between adjacent carboxyl groups as well as to study the dependence of this balance on steric confinement and the effect of this balance on intermolecular interactions of the carboxyl groups.

Adduct under Field-A Qualitative Approach to Account for Solvent Effect on Hydrogen Bonding

The location of a mobile proton in acid-base complexes in aprotic solvents can be predicted using a simplified Adduct under Field (AuF) approach, where solute–solvent effects on the geometry of hydrogen bond are simulated using a fictitious external electric field. The parameters of the field have been estimated using experimental data on acid-base complexes in CDF₃/CDClF₂. With some limitations, they can be applied to the chemically similar CHCl₃ and CH₂Cl₂. The obtained data indicate that the solute–solvent effects are critically important regardless of the type of complexes. The temperature dependences of the strength and fluctuation rate of the field explain the behavior of experimentally measured parameters.

Solvent effects on acid-base complexes. What is more important: A macroscopic reaction field or solute-solvent interactions?

Can the geometry of an acid-base complex in solution be reproduced in calculations using an implicit accounting for the solvent effect in the form of a macroscopic reaction field? The answer is, "Yes, it can." Is this field equal to the real electric field experienced by the complex in solution? The answer is, "No, it is not." How can the geometry be correct under wrong conditions? This question is answered using density functional theory modeling of geometric and NMR parameters of pyridine⋯HF⋯(HCF₃)_n adducts in the absence and presence of an external electric field. This adduct under field approach shows that the N⋯H distance is a function of the H-F distance whatever method is used to change the geometry of the latter. An explicit account for solute-solvent interactions is required to get a realistic value of the solvent reaction field. Besides that, this approach reveals how certain NMR parameters depend on the solvent reaction field, the solute-solvent interactions, and the geometry of the N⋯H-F hydrogen bond. For some of them, the obtained dependences are far from self-evident.

Local-structure effects on 31P NMR chemical shift tensors in solid state

The effect of the local structure on the ³¹P NMR chemical shift tensor (CST) has been studied experimentally and simulated theoretically using the density functional theory gauge-independent-atomic-orbital approach. It has been shown that the dominating impact comes from a small number of noncovalent interactions between the phosphorus-containing group under question and the atoms of adjacent molecules. These interactions can be unambiguously identified using the Bader analysis of the electronic density. A robust and computationally effective approach designed to attribute a given experimental ³¹P CST to a certain local morphology has been elaborated. This approach can be useful in studies of surfaces, complex molecular systems, and amorphous materials.

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.

Symmetry and dynamics of FHF− anion in vacuum, in CD2Cl2 and in CCl4. Ab initio MD study of fluctuating solvent–solute hydrogen and halogen bonds

Asymmetric solvation of FHF⁻ by halogen- and hydrogen-bonding solvents breaks the symmetry of the anion.

Theoretical study of the H/D isotope effect on phase transition of hydrogen-bonded organic conductor κ-H3(Cat-EDT-TTF)2

κ-H₃(Cat-EDT-TTF)₂ (H-TTF) is a hydrogen-bonded π-electron system. Only its isotopologue, D-TTF, shows the phase transition. We obtained a symmetric single-well effective-PEC for H-TTF and low-barrier effective-PEC for D-TTF.

Theoretical analysis of geometry and NMR isotope shift in hydrogen-bonding center of photoactive yellow protein by combination of multicomponent quantum mechanics and ONIOM scheme

Multicomponent quantum mechanical (MC_QM) calculation has been extended with ONIOM (our own N-layered integrated molecular orbital + molecular mechanics) scheme [ONIOM(MC_QM:MM)] to take account of both the nuclear quantum effect and the surrounding environment effect. The authors have demonstrated the first implementation and application of ONIOM(MC_QM:MM) method for the analysis of the geometry and the isotope shift in hydrogen-bonding center of photoactive yellow protein. ONIOM(MC_QM:MM) calculation for a model with deprotonated Arg52 reproduced the elongation of O-H bond of Glu46 observed by neutron diffraction crystallography. Among the unique isotope shifts in different conditions, the model with protonated Arg52 with solvent effect reasonably provided the best agreement with the corresponding experimental values from liquid NMR measurement. Our results implied the availability of ONIOM(MC_QM:MM) to distinguish the local environment around hydrogen bonds in a biomolecule.

Energy Analysis of Competing Non-Covalent Interaction in 1:1 and 1:2 Adducts of Collidine with Benzoic Acids by Means of X-Ray Diffraction

The hydrogen bond pattern and the types of non-covalent interactions in the crystals of the 1:1 and 1:2 adducts of 2,4,6-trimethylpyridine and benzoic acids are studied using high-resolution X-ray diffraction. The geometries of the hydrogen bonds are estimated using a combined XRD/DFT approach that provides the geometrical parameters within the margin of error of neutron diffraction studies. The energies of the non-covalent interactions are estimated on the base of the experimental electron density distribution function. It is shown that the structures of the adducts are governed by the NOH and OHO hydrogen bonds. In turn, C-H...O contacts and stacking interactions define the packing of the adducts in the crystal. On the other hand, it is important to note that the latter interactions affect the competition of the former hydrogen bonds in some 1:2 adducts.

NMR study of conformational exchange and double-well proton potential in intramolecular hydrogen bonds in monoanions of succinic acid and derivatives

We present a (1)H, (2)H, and (13)C NMR study of the monoanions of succinic (1), meso- and rac-dimethylsuccinic (2, 3), and methylsuccinic (4) acids (with tetraalkylammonium as the counterion) dissolved in CDF(3)/CDF(2)Cl at 300-120 K. In all four monoanions, the carboxylic groups are linked by a short intramolecular OHO hydrogen bond revealed by the bridging-proton chemical shift of about 20 ppm. We show that the flexibility of the carbon skeleton allows for two gauche isomers in monoanions 1, 2, and 4, interconverting through experimental energy barriers of 10-15 kcal/mol (the process itself and the energy barrier are also reproduced in MP2/6-311++G** calculations). In 3, one of the gauche forms is absent because of the steric repulsion of the methyl groups. In all four monoanions, the bridging proton is located in a double-well potential and subject, at least to some extent, to proton tautomerism, for which we estimate the two proton positions to be separated by ca. 0.2 Å. In 1 and 3, the proton potential is symmetric. In 2, slowing the conformational interconversion introduces an asymmetry to the proton potential, an effect that might be strong enough even to synchronize the proton tautomerism with the interconversion of the two gauche forms. In 4, the asymmetry of the proton potential is due to the asymmetric substitution. The intramolecular H-bond is likely to remain intact during the interconversion of the gauche forms in 1, 3, and 4, whereas the situation in 2 is less clear.

NMR Study of Blue-Shifting Hydrogen Bonds Formed by Fluoroform in Solution

Low-temperature (193 K) 1H, 13C and 15N NMR spectra of blue- and red-shifting H-bonded complexes formed by fluoroform with various proton acceptors were measured. Experimental NMR parameters were plotted versus the ab initio calculated H-bond strength (MP2/6–31+G(d, p); interaction energy varies from ~5 up to 25 kJ · mol−1 in the series). We show that experimental 1H and 15N shieldings, as well as the H/D isotope effect on 13C shielding change monotonously with the calculated H-bond strengthening. The 13C chemical shift and the CH scalar coupling change non-monotonously and the extremum points are situated approximately in the region of transformation from blue- to red-shifting H-bonds. The most informative NMR feature is the H/D isotope effect on 15N shielding which changes its sign upon transformation from blue- to red-shifting H-bonds.

To rationalize these observations, ab initio calculations of 13C and 15N shieldings as functions of C···H and C···N distances were performed for complexes of CHF3 with acetonitrile (blue-shifting) and pyridine (red-shifting). The coupling of the vibrations of the covalent and hydrogen bonds has been accounted for by direct computation of the distance C···N = f(C···H) dependence. We demonstrate that the unusual sign of the H/D isotope effect on 15N chemical shift across a blue-shifting H-bond can be explained as a result of the inversion of the dynamic coupling of two vibrations.

Geometric H/D Isotope Effects and Cooperativity of the Hydrogen Bonds in Porphycene

We investigate the primary, secondary, and vicinal hydrogen/deuterium (H/D) isotope effects on the geometry of the two intramolecular hydrogen bonds in porphycene. Multidimensional potential energy surfaces describing the anharmonic motion in the vicinity of the trans isomer are calculated for the different symmetric (HH/DD) and asymmetric (HD) isotopomers. From the solution of the nuclear Schrödinger equation the ground-state wavefunction is obtained, which is further used to determine the quantum corrections to the classical equilibrium geometries of the hydrogen bonds and thus the geometric isotope effects. In particular, it is found that the hydrogen bonds are cooperative, that is, both expand simultaneously even in the case of an asymmetric isotopic substitution. The theoretical predictions compare favorably with NMR chemical-shift data.

Geometry and Cooperativity Effects in Adenosine−Carboxylic Acid Complexes

NMR experiments and theoretical investigations were performed on hydrogen bonded complexes of specifically 1- and 7-15N-labeled adenine nucleosides with carboxylic acids. By employing a freonic solvent of CDClF2 and CDF3, NMR spectra were acquired at temperatures as low as 123 K, where the regime of slow hydrogen bond exchange is reached and several higher-order complexes were found to coexist in solution. Unlike acetic acid, chloroacetic acid forms Watson-Crick complexes with the proton largely displaced from oxygen to the nitrogen acceptor in an ion pairing structure. Calculated geometries and chemical shifts of the proton in the hydrogen bridge favorably agree with experimentally determined values if vibrational averaging and solvent effects are taken into account. The results indicate that binding a second acidic ligand at the adenine Hoogsteen site in a ternary complex weakens the hydrogen bond to the Watson-Crick bound carboxylic acid. However, substituting a second adenine nucleobase for a carboxylic acid in the trimolecular complex leads to cooperative binding at Watson-Crick and Hoogsteen faces of adenosine.