What is a hydrogen bond? Mutually consistent theoretical and experimental criteria for characterizing H-bonding interactions

Hydrogen bond types which do not fit accepted definitions

There are various interactions that either partially fit or do not fit the currently accepted definitions of the hydrogen bond. However, they possess characteristics of this interaction. It seems that it is partly connected to the fact that these definitions are not precise. The typical 3c-4e (three centres - four electrons) A-H⋯B hydrogen bond is characterized by the single-atom A and B centres that are highly electronegative. On the other hand, non-typical interactions that do not fit the hydrogen bond definitions well are characterised by uncommon proton donors and/or proton acceptors. The cases of multi-centre proton acceptors, π-electron or σ-electron systems are well known - such interactions are designated as A-H⋯π and A-H⋯σ hydrogen bonds, respectively. However, the cases of interactions with the multi-centre proton donors and proton acceptors do not fit the majority of definitions of hydrogen bond. The π⋯H+⋯π system in the proton-bound homodimer of acetylene is an example. This system can be classified as a hydrogen bond according to the two-sites hydrogen bond, 2sHB, definition. There are various types of interactions discussed in this review; among them, those that are undoubtedly unclassified as hydrogen bonds, i.e., hydride bonds, and charge inverted hydrogen bonds, CIHBs. Special emphasis is also put here on the proton sponges and other systems such as the [FHF]- anion or [NgHNg]+ cation (Ng is the noble gas centre).

Blueshift in Trifurcated Hydrogen Bonds: A Tradeoff between Tetrel Bonding and Steric Repulsion

We have quantum chemically investigated the origin of the atypical blueshift of the H-C bond stretching frequency in the hydrogen-bonded complex X- •••H3 C-Y (X, Y=F, Cl, Br, I), as compared to the corresponding redshift occurring in Cl- •••H3 N and Cl- •••H3 C-H, using relativistic density functional theory (DFT) at ZORA-BLYP-D3(BJ)/QZ4P. Previously, this blueshift was attributed, among others, to the contraction of the H-C bonds as the H3 C moiety becomes less pyramidal. Herein, we provide quantitative evidence that, instead, the blueshift arises from a direct and strong X- •••C interaction of the HOMO of A- with the backside lobe on carbon of the low-lying C-Y antibonding σ* LUMO of the H3 C-Y fragment. This X- •••C bond, in essence a tetrel bond, pushes the H atoms towards a shorter H-C distance and makes the H3 C moiety more planar. The blueshift may, therefore, serve as a diagnostic for tetrel bonding.

Effects of Cell-Adhesive Ligand Presentation on Pentapeptide Supramolecular Assembly and Gelation: Simulations and Experiments

The extracellular matrix (ECM) is a complex, hierarchical material containing structural and bioactive components. This complexity makes decoupling the effects of biomechanical properties and cell-matrix interactions difficult, especially when studying cellular processes in a 3D environment. Matrix mechanics and cell adhesion are both known regulators of specific cellular processes such as stem cell proliferation and differentiation. However, more information is required about how such variables impact various neural lineages that could, upon transplantation, therapeutically improve neural function after a central nervous system injury or disease. Rapidly Assembling Pentapeptides for Injectable Delivery (RAPID) hydrogels are one biomaterial approach to meet these goals, consisting of a family of peptide sequences that assemble into physical hydrogels in physiological media. In this study, we studied our previously reported supramolecularly-assembling RAPID hydrogels functionalized with the ECM-derived cell-adhesive peptide ligands RGD, IKVAV, and YIGSR. Using molecular dynamics simulations and experimental rheology, we demonstrated that these integrin-binding ligands at physiological concentrations (3-12 mm) did not impact the assembly of the KYFIL peptide system. In simulations, molecular measures of assembly such as hydrogen bonding and pi-pi interactions appeared unaffected by cell-adhesion sequence or concentration. Visualizations of clustering and analysis of solvent-accessible surface area indicated that the integrin-binding domains remained exposed. KYFIL or AYFIL hydrogels containing 3 mm of integrin-binding domains resulted in mechanical properties consistent with their non-functionalized equivalents. This strategy of doping RAPID gels with cell-adhesion sequences allows for the precise tuning of peptide ligand concentration, independent of the rheological properties. The controllability of the RAPID hydrogel system provides an opportunity to investigate the effect of integrin-binding interactions on encapsulated neural cells to discern how hydrogel microenvironment impacts growth, maturation, or differentiation.

Valence Bond Theory Allows a Generalized Description of Hydrogen Bonding

This paper describes the nature of the hydrogen bond (HB), B:---H-A, using valence bond theory (VBT). Our analysis shows that the most important HB interactions are polarization and charge transfer, and their corresponding sum displays a pattern that is identical for a variety of energy decomposition analysis (EDA) methods. Furthermore, the sum terms obtained with the different EDA methods correlate linearly with the corresponding VB quantities. The VBT analysis demonstrates that the total covalent-ionic resonance energy (RECS) of the HB portion (B---H in B:---H-A) correlates linearly with the dissociation energy of the HB, ΔEdiss. In principle, therefore, RECS(HB) can be determined by experiment. The VBT wavefunction reveals that the contributions of ionic structures to the HB increase the positive charge on the hydrogen of the corresponding external/free O-H bonds in, for example, the water dimer compared with a free water molecule. This increases the electric field of the external O-H bonds of water clusters and contributes to bringing about catalysis of reactions by water droplets and in water-hydrophobic interfaces.

C-H Groups as Donors in Hydrogen Bonds: A Historical Overview and Occurrence in Proteins and Nucleic Acids

Hydrogen bonds constitute a unique type of non-covalent interaction, with a critical role in biology. Until fairly recently, the canonical view held that these bonds occur between electronegative atoms, typically O and N, and that they are mostly electrostatic in nature. However, it is now understood that polarized C-H groups may also act as hydrogen bond donors in many systems, including biological macromolecules. First recognized from physical chemistry studies, C-H…X bonds were visualized with X-ray crystallography sixty years ago, although their true significance has only been recognized in the last few decades. This review traces the origins of the field and describes the occurrence and significance of the most important C-H…O bonds in proteins and nucleic acids.

"Noncovalent Interaction": A Chemical Misnomer That Inhibits Proper Understanding of Hydrogen Bonding, Rotation Barriers, and Other Topics

We discuss the problematic terminology of "noncovalent interactions" as commonly applied to hydrogen bonds, rotation barriers, steric repulsions, and other stereoelectronic phenomena. Although categorization as "noncovalent" seems to justify classical-type pedagogical rationalizations, we show that these phenomena are irreducible corollaries of the same orbital-level conceptions of electronic covalency and resonance that govern all chemical bonding phenomena. Retention of such nomenclature is pedagogically misleading in supporting superficial dipole-dipole and related "simple, neat, and wrong" conceptions as well as perpetuating inappropriate bifurcation of the introductory chemistry curriculum into distinct "covalent" vs. "noncovalent" modules. If retained at all, the line of dichotomization between "covalent" and "noncovalent" interaction should be re-drawn beyond the range of quantal exchange effects (roughly, at the contact boundary of empirical van der Waals radii) to better unify the pedagogy of molecular and supramolecular bonding phenomena.

The Role of H-Bonds in the Excited-State Properties of Multichromophoric Systems: Static and Dynamic Aspects

Given their importance, hydrogen bonds (H-bonds) have been the subject of intense investigation since their discovery. Indeed, H-bonds play a fundamental role in determining the structure, the electronic properties, and the dynamics of complex systems, including biologically relevant materials such as DNA and proteins. While H-bonds have been largely investigated for systems in their electronic ground state, fewer studies have focused on how the presence of H-bonds could affect the static and dynamic properties of electronic excited states. This review presents an overview of the more relevant progress in studying the role of H-bond interactions in modulating excited-state features in multichromophoric biomimetic complex systems. The most promising spectroscopic techniques that can be used for investigating the H-bond effects in excited states and for characterizing the ultrafast processes associated with their dynamics are briefly summarized. Then, experimental insights into the modulation of the electronic properties resulting from the presence of H-bond interactions are provided, and the role of the H-bond in tuning the excited-state dynamics and the related photophysical processes is discussed.

Exploring the Non-Covalent Bonding in Water Clusters

QTAIM and source function analysis were used to explore the non-covalent bonding in twelve different water clusters (H2O)n obtained by considering n = 2–7 and various geometrical arrangements. A total of seventy-seven O−H⋯O hydrogen bonds (HBs) were identified in the systems under consideration, and the examination of the electron density at the bond critical point (BCP) of these HBs revealed the existence of a great diversity of O−H⋯O interactions. Furthermore, the analysis of quantities, such as |V(r)|/G(r) and H(r), allowed a further description of the nature of analogous O−H⋯O interactions within each cluster. In the case of 2-D cyclic clusters, the HBs are nearly equivalent between them. However, significant differences among the O−H⋯O interactions were observed in 3-D clusters. The assessment of the source function (SF) confirmed these findings. Finally, the ability of SF to decompose the electron density (ρ) into atomic contributions allowed the evaluation of the localized or delocalized character of these contributions to ρ at the BCP associated to the different HBs, revealing that weak O−H⋯O interactions have a significant spread of the atomic contributions, whereas strong interactions have more localized atomic contributions. These observations suggest that the nature of the O−H⋯O hydrogen bond in water clusters is determined by the inductive effects originated by the different spatial arrangements of the water molecules in the studied clusters.

Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions

Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.

Molecular Motion in the Interconverting σ-H2, Di-, and Tri-hydride Regimes: Mo(PH3)5H2

The transition-metal complex Mo(PH3)5H2 is known to exist in three possible isomeric forms, including a nonclassical, σ-bound dihydrogen complex and two classical dihydride isomers. As such, it has served as a model complex for the energies of conversion between these limiting structural regimes. In the present study, ab initio molecular dynamics computer simulations, combined with enhanced sampling techniques, were utilized to directly assess the degree of motion and isomerization of the dihydrogen/dihydride moieties in this complex. Ligand rotations (for both the H₂ unit and the phosphine units) were found to be dominant in the low-temperature (298 K) regime, and the classical thermodynamic distribution showed no probability of thermally accessing dihydride forms, although unrestrained molecular dynamics trajectories showed fleeting configurations outside of the σ-H₂ configuration. Simulations at higher temperatures surprisingly revealed new tri-hydride isomers that are energetically competitive with the σ-H₂ and cis-/trans-dihydride isomers. Low-energy pathways to hydrogen/hydride transfer and phosphine dissociation were readily accessible, which considerably expands the known isomeric flexibility of this complex.

Sandwich, Triple-Decker and Other Sandwich-like Complexes of Cyclopentadienyl Anions with Lithium or Sodium Cations

Density functional theory, DFT, calculations were carried out on complexes containing cyclopentadienyl anions and lithium or sodium cations; half-sandwich, sandwich and sandwich-like complexes (among them triple-decker ones) are analyzed. Searches performed through the Cambridge Structural Database revealed that crystal structures containing these motifs exist, mostly structures with sodium cations. The DFT calculations performed here include geometry optimization and frequency calculations of the complexes at the ωB97XD/aug-cc-pVTZ level, followed by the partitioning of the energy of interaction via the Energy Decomposition Analysis scheme, EDA, at the BP86-D3/TZ2P level. Additional calculations and analyses were performed using both the Quantum Theory of Atoms in Molecules, QTAIM, and the Natural Bond Orbital analyses, NBO. The results of this work show that the electrostatic interaction energy is the most important attractive contribution to the total interaction energy of each of the complex systems analyzed here, and that complexation itself leads to minor electron charge shifts.

Estimations of OH·N hydrogen bond length from positions and intensities of IR bands

In this computational work applicability of IR spectral parameters for evaluations of OH···N hydrogen bond length is discussed. For a set of 124 complexes with OH···N hydrogen bond formed by combinations of methanol/acetic acid and pyridine (and their fluorine substituted versions) geometries, energies and IR parameters were calculated at MP2/def2-TZVP level of theory. For a number of IR parameters (the shift of proton donor group stretching vibration Δν_s, increase of its intensity I, the low-frequency hydrogen bond stretching vibration ν_σ, bending in-plane δ and out-of-plane γ vibrations) equations linking them with interatomic distances are proposed, the robustness and accuracy of such equations are discussed. The enthalpy of OH···N hydrogen bond formation ΔH was also linked with electron density parameters in (3; -1) critical point.

High-Density “Windowpane” Coordination Patterns of Water Clusters and Their NBO/NRT Characterization

Cluster mixture models for liquid water at higher pressures suggest the need for water clusters of higher coordination and density than those commonly based on tetrahedral H-bonding motifs. We show here how proton-ordered water clusters of increased coordination and density can assemble from a starting cyclic tetramer or twisted bicyclic (Möbius-like) heptamer to form extended Aufbau sequences of stable two-, three-, and four-coordinate “windowpane” motifs. Such windowpane clusters exhibit sharply reduced (~90°) bond angles that differ appreciably from the tetrahedral angles of idealized crystalline ice I_h. Computed free energy and natural resonance theory (NRT) bond orders provide quantitative descriptors for the relative stabilities of clusters and strengths of individual coordinative linkages. The unity and consistency of NRT description is demonstrated to extend from familiar supra-integer bonds of the molecular regime to the near-zero bond orders of the weakest linkages in the present H-bond clusters. Our results serve to confirm that H-bonding exemplifies resonance–covalent (fractional) bonding in the sub-integer range and to further discount the dichotomous conceptions of “electrostatics” for intermolecular bonding vs. “covalency” for intramolecular bonding that still pervade much of freshman-level pedagogy and force-field methodology.

The proton and the lithium cation linked with π-electron and σ-electron systems: are such interactions beyond or within the definition of hydrogen/lithium bond?

MP2/aug-cc-pVTZ calculations were performed on systems containing the proton or the lithium cation located between two π-electron systems or between π-electron and σ-electron units. The proton or the lithium cation attached to the acetylene or its derivative may be treated as the Lewis acid unit while the remaining part of the complex, the π-electron species or the dihydrogen, act as the Lewis base through their π-electrons or σ-electrons, respectively. The complexes analysed here are linked by the π∙∙∙H + /Li + ∙∙∙π and π∙∙∙H + /Li + ∙∙∙σ interactions. It is discussed whether these interactions are classified as hydrogen and lithium bonds. Therefore, different definitions of the latter interactions are presented. The Electron Localization Function (ELF) and the Natural Bond Orbital (NBO) approaches were applied to analyse the above-mentioned complexes. The unique properties of interactions with the proton and with the lithium cation that occur in complexes analysed here are described.

The chemical bond in solids-revisited

This article complements an earlier topical review of the chemical bond (Jones 2018J. Phys.: Condens. Matter30153001), starting in the mid-19th century and seen from the perspective of a condensed matter physicist. The discussion of applications focused on the structure and properties of phase change materials. We review here additional aspects of chemistry, particularly some that have raised interest recently in this context. Concepts such as 'electron-rich', 'electron-deficient (excess orbital)', 'hypervalent', 'three-centre', and 'metavalent' bonds, and 'multicentre hyperbonding' are now found in the condensed matter literature. They are surveyed here, as well as the bond in metals and the 'Peierls' distortion. What are these concepts, are they related, and are they sometimes new labels for established, but unfamiliar ideas? 'Half bonds' and 'fractional valencies' play a central role in this discussion. It is remarkable that they were introduced 100 years ago, but ignored or forgotten, and have needed to be rediscovered more than once.

Influence of an Oriented External Electric Field on the Mechanism of Double Proton Transfer between Pyrazole and Guanidine: from an Asynchronous Plateau Transition State to a Synchronous or Stepwise Mechanism

The double proton transfer (DPT) reaction between pyrazole and guanidine, a concerted reaction but strongly asynchronous and presenting a "plateau transition region", has been theoretically reinvestigated in the presence of an external uniform electric field. First, we computed the reaction path by DFT and proposed a very detailed description of the constitutive electronic events, based on the ELF topology and the bond evolution theory. Then, we studied the effect of an oriented external electric field (OEEF) on the reaction mechanism, for an OEEF oriented along the proton transfer axis. We observe that in one direction, the DPT reaction can be transformed into a stepwise reaction, going through a stabilized single proton transferred intermediate. Contrarily, the two proton transfers occur simultaneously when the electric field is applied in the opposite direction. In the latter case, the order in which the two protons are transferred in the same elementary step can even be reversed if the OEEF is intense enough. Finally, it has been shown that the evolution of the double proton transfer reaction in the presence of an electric field could be quantitatively anticipated by analyzing the ELF value at the bifurcation point between V(A, H) proton donor and V(B) proton acceptor of the double hydrogen bonded complex in the entrance channel.

Chlorine Dioxide: An Exception that Proves the Rules of Localized Chemical Bonding

We employ natural bond orbital (NBO) and natural resonance theory (NRT) tools to analyze the enigmatic properties of the C2v-symmetric isomer of chlorine dioxide radical (ClO2), whose many challenges to Pauling-type localized bonding concepts were recognized by Linus Pauling himself. Although spin-contamination is minimal in this species, ClO2 exhibits an unusually strong form of "different Lewis structures for different spins" bonding pattern, intrinsically outside the framework of "maximal pairing" concepts. We show how the novel spin-unpaired donor-acceptor interactions lead to weakened bonding in the supramolecular domain of polyradical (ClO2)n homoclusters and aqueous ClO2(H2O)n heteroclusters. Despite feeble binding energies and large inter-radical separations, the polyradical clusters are found to maintain coherent spin patterns in each cluster component, attesting to the quantal donor-acceptor nature of their interactions and the cooperative and anticooperative couplings that govern intra- and intermolecular spin distributions in such spin-clusters.

Accelerating therapeutic protein design

Protein structures provide for defined microenvironments that can support complex pharmacological functions, otherwise unachievable by small molecules. The advent of therapeutic proteins has thus greatly broadened the range of manageable disorders. Leveraging the knowledge and recent advances in de novo protein design methods has the prospect of revolutionizing how protein drugs are discovered and developed. This review lays out the main challenges facing therapeutic proteins discovery and development, and how present and future advancements of protein design can accelerate the protein drug pipelines.

Anti-Electrostatic Pi-Hole Bonding: How Covalency Conquers Coulombics

Intermolecular bonding attraction at π-bonded centers is often described as "electrostatically driven" and given quasi-classical rationalization in terms of a "pi hole" depletion region in the electrostatic potential. However, we demonstrate here that such bonding attraction also occurs between closed-shell ions of like charge, thereby yielding locally stable complexes that sharply violate classical electrostatic expectations. Standard DFT and MP2 computational methods are employed to investigate complexation of simple pi-bonded diatomic anions (BO-, CN-) with simple atomic anions (H-, F-) or with one another. Such "anti-electrostatic" anion-anion attractions are shown to lead to robust metastable binding wells (ranging up to 20-30 kcal/mol at DFT level, or still deeper at dynamically correlated MP2 level) that are shielded by broad predissociation barriers (ranging up to 1.5 Å width) from long-range ionic dissociation. Like-charge attraction at pi-centers thereby provides additional evidence for the dominance of 3-center/4-electron (3c/4e) nD-π*AX interactions that are fully analogous to the nD-σ*AH interactions of H-bonding. Using standard keyword options of natural bond orbital (NBO) analysis, we demonstrate that both n-σ* (sigma hole) and n-π* (pi hole) interactions represent simple variants of the essential resonance-type donor-acceptor (Bürgi-Dunitz-type) attraction that apparently underlies all intermolecular association phenomena of chemical interest. We further demonstrate that "deletion" of such π*-based donor-acceptor interaction obliterates the characteristic Bürgi-Dunitz signatures of pi-hole interactions, thereby establishing the unique cause/effect relationship to short-range covalency ("charge transfer") rather than envisioned Coulombic properties of unperturbed monomers.

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.

Modeling the peroxyl radical scavenging behavior of Carnosic acid: Mechanism, kinetics, and effects of physiological environments

Carnosic acid (CA), a phenolic diterpene and abietane-type compound, is a potent natural antioxidant with medical benefits. The present paper elucidates, for the first time, the kinetics and the exact mechanism of the peroxyl radical scavenging activity of CA in the gas phase and under physiological conditions. According to the obtained results, the reaction of CA with HOO^• is significantly faster in aqueous solution than in the gas phase and nonpolar environments. The abstraction of the hydrogen atom from 2-OH is the decisive mechanism in the gas phase and nonpolar media, while both hydrogen abstraction (15%) and electron transfer (85%) mechanisms can take place in aqueous solution. The overall rate coefficient in water (4.73 × 10⁶ M^-1 s^-1) is about 36 times higher than that of the reference antioxidant Trolox (1.30 × 10⁵ M^-1 s^-1), suggesting that CA is a potent scavenger of peroxyl radicals in polar media.

Non-covalent interactions in molecular systems: thermodynamic evaluation of the hydrogen bond strength in aminoalcohols

In molecules with two functional groups that form hydrogen bonds, the structure-property relationship can depend significantly on the strength of intra-molecular hydrogen bonding. This bonding can cause a substantial conformational change that is accompanied by a frequency shift in the infrared spectrum, which provides the basis for experimental studies. Despite its great importance in biological systems, the available literature data for the strength of this bonding are scarce and not in agreement. In this work, we present the results of four thermodynamic methods for the determination of the strength of intramolecular hydrogen bonds. Comprehensive thermochemical analysis of 1-amino-2-alcohols and 2-amino-1-alcohols was performed with Fourier-transform infrared spectroscopy, high-level G4 quantum-chemical calculations, the homomorph scheme with enthalpies of vaporization and a group contribution method. With the combination of these four thermodynamic methods, the strength of intramolecular hydrogen bonding in 1,2-aminoalcohols and 2,1-aminoalcohols was evaluated quantitatively. The results were correlated with NBO parameters to find an explanation for the different strengths of intramolecular hydrogen bonds in total charge transfer and second order stabilization energies.

1H NMR Studies of Intramolecular OH/OH Hydrogen Bonds via Titratable Isotope Shifts

Methanol titrations of partially deuterated 1,4- and 1,3-diols dissolved in nonpolar solvents such as CD₂Cl₂ and benzene-d₆ have provided ¹H NMR measurements of OH/OD isotope shifts, diagnostic for intact intramolecular hydrogen bonds, under conditions of increasing protic solvent concentration. 1,4- and 1,3-diols with conformationally favored intramolecular OH/OH hydrogen bonds can be titrated to constant isotope shift values, albeit with variable sign, in the presence of excess methanol equivalents, providing evidence for intact intramolecular hydrogen bonds under these conditions. Conversely, the isotope shift in a 1,3-diol with a conformationally labile intramolecular hydrogen bond titrated to zero when in the presence of excess equivalents methanol, consistent with intramolecular hydrogen bond rupture under these conditions. Additionally, the titration behavior of hydroxyl chemical shifts in diols and protected derivatives has revealed significant OH/OD isotope shifts in the absence of chemical shift differences (δOH_in = δOH_out) that are necessary for an equilibrium isotope effect, lending evidence for an intrinsic contribution to the isotope effect. OH/OD isotope shift titration thus provides a means for understanding the origins of these isotope effects and for probing the intact or nonintact nature of intramolecular OH/OH hydrogen bonds in response to intermolecular hydrogen bonds provided by a protic solvent.

Three in One: The Versatility of Hydrogen Bonding Interaction in Halide Salts with Hydroxy-Functionalized Pyridinium Cations

The paradigm of supramolecular chemistry relies on the delicate balance of noncovalent forces. Here we present a systematic approach for controlling the structural versatility of halide salts by the nature of hydrogen bonding interactions. We synthesized halide salts with hydroxy-functionalized pyridinium cations [HOCn Py]+ (n=2, 3, 4) and chloride, bromide and iodide anions, which are typically used as precursor material for synthesizing ionic liquids by anion metathesis reaction. The X-ray structures of these omnium halides show two types of hydrogen bonding: 'intra-ionic' H-bonds, wherein the anion interacts with the hydroxy group and the positively charged ring at the same cation, and 'inter-ionic' H-bonds, wherein the anion also interacts with the hydroxy group and the ring system but of different cations. We show that hydrogen bonding is controllable by the length of the hydroxyalkyl chain and the interaction strength of the anion. Some molten halide salts exhibit a third type of hydrogen bonding. IR spectra reveal elusive H-bonds between the OH groups of cations, showing interaction between ions of like charge. They are formed despite the repulsive interaction between the like-charged ions and compete with the favored cation-anion H-bonds. All types of H-bonding are analyzed by quantum chemical methods and the natural bond orbital approach, emphasizing the importance of charge transfer in these interactions. For simple omnium salts, we evidenced three distinct types of hydrogen bonds: Three in one!

Neat, Simple, and Wrong: Debunking Electrostatic Fallacies Regarding Noncovalent Interactions

Multipole moments such as charge, dipole, and quadrupole are often invoked to rationalize intermolecular phenomena, but a low-order multipole expansion is rarely a valid description of electrostatics at the length scales that characterize nonbonded interactions. This is illustrated by examining several common misunderstandings rooted in erroneous electrostatic arguments. First, the notion that steric repulsion originates in Coulomb interactions is easily disproved by dissecting the interaction potential for Ar₂. Second, the Hunter-Sanders model of π-π interactions, which is based on quadrupolar electrostatics, is shown to have no basis in accurate calculations. Third, curved "buckybowls" exhibit unusually large dipole moments, but these are ancillary to the forces that control their intermolecular interactions, as illustrated by two examples involving corannulene. Finally, the assumption that interactions between water and small anions are dictated by the dipole moment of H₂O is shown to be false in the case of binary halide-water complexes. These examples present a compelling case that electrostatic explanations based on low-order multipole moments are very often counterfactual for nonbonded interactions at close range and should not be taken seriously in the absence of additional justification.

Highs and Lows of Bond Lengths: Is There Any Limit?

Two distinct points on the potential energy curve (PEC) of a pairwise interaction, the zero‐energy crossing point and the point where the stretching force constant vanishes, allow us to anticipate the range of possible distances between two atoms in diatomic, molecular moieties and crystalline systems. We show that these bond‐stability boundaries are unambiguously defined and correlate with topological descriptors of electron‐density‐based scalar fields, and can be calculated using generic PECs. Chemical databases and quantum‐mechanical calculations are used to analyze a full set of diatomic bonds of atoms from the s‐p main block. Emphasis is placed on the effect of substituents in C−C covalent bonds, concluding that distances shorter than 1.14 Å or longer than 2.0 Å are unlikely to be achieved, in agreement with ultra‐high‐pressure data and transition‐state distances, respectively. Presumed exceptions are used to place our model in the correct framework and to formulate a conjecture for chained interactions, which offers an explanation for the multimodal histogram of O−H distances reported for hundreds of chemical systems.

Comment on "Superposition of Waves or Densities: Which Is the Nature of Chemical Resonance?" [J. Comput. Chem. 2021, 42, 412-417]

We reply to specific criticisms and misrepresentations of natural resonance theory (NRT) in a recent article [Y. Wang, J. Comput. Chem. 2021, 42, 412-417] and argue that it presents a false dichotomy with respect to theoretical efforts to comprehend the nature of resonance-type phenomena.

Interfacial interactions between protective, surface-engineered shells and encapsulated bacteria with different cell surface composition

Surface-engineered encapsulation is a non-genetic method to protect living organisms against harsh environmental conditions. Different cell encapsulation methods exist, yielding shells with different interfacial-interactions with encapsulated, bacterial surfaces. However, the impact of interfacial-interactions on the protection offered by different shells is unclear and can vary for bacteria with different surface composition. Probiotic bacteria require protection against gastro-intestinal fluids and antibiotics. Here, we encapsulated two probiotic strains using ZIF-8 (zeolitic imidazolate framework) biomineralization (strong-interaction by coordinate-covalent bonding), alginate gelation (intermediate-interaction by hydrogen bonding) or protamine-assisted packing of SiO₂ nanoparticles yielding a yolk-shell (weak-interaction across a void between shells and bacterial surfaces). The surface of probiotic Lactobacillus acidophilus was rich in protein, yielding a hydrophilic, positively-charged surface below and a negatively-charged one above pH 4.0. Probiotic Bifidobacterium infantis had a hydrophilic, uncharged surface, rich in polysaccharides with little proteins. Although amino groups are required for coordinate-covalent bonding of zinc and hydrogen bonding of alginate, both L. acidophilus and B. infantis could be encapsulated using ZIF-8 biomineralization and alginate gelation. Weakly, intermediately and strongly interacting shells all yielded porous shells. The strongly interacting ZIF-8 biomineralized shell made encapsulated bacteria more susceptible to antibiotics, presumably due to the cell wall damage already inflicted during Zif-8 biomineralization. Overall, weakly interacting yolk-shells and intermediately interacting alginate gels protected best and maintained probiotic activity of encapsulated bacteria. The impact of interfacial-interactions between shells and encapsulated bacteria on different aspect of protection described here, contributes to the further development of effective surface-engineered shells and its application for protecting bacteria.

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N₂, and BF bound to a [Ru(II)(NH₃)₅]^{2 +} moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

Accurate location of hydrogen atoms in hydrogen bonds of tizoxanide from the combination of experimental and theoretical models

To obtain detailed information about the position of hydrogen atoms in hydrogen bonds, HBs, of crystalline organic molecular compounds is not an easy task. In this work we propose a combination of ssNMR experimental data with theoretical procedures to get such information. Furthermore, the combination of experimental and theoretical models provides us with well-defined grounds to analyse the strength of π-stacking interactions between layers of hydrogen bonded molecules. Two different theoretical models were considered, both approaches being quite different. The first one is a solid-state model, so that the periodicity of a crystalline system underlies calculations of the electronic energy, the electronic density and NMR parameters. The other one is a molecular model in which molecules are taken as isolated monomers, dimers and tetramers. These two models were applied to the tizoxanide, TIZ, molecular crystal though it can widely be applied to any other molecular crystal. By the application of the quantum molecular model it was possible to learn about the way the intermolecular HBs affect the position of hydrogen atoms that belong to HBs in TIZ. This molecule has two intermolecular HBs that stabilize the structure of a basic dimer, but it also has an intramolecular HB in each monomer whose position should be optimized together with the other ones. We found that by doing this it is possible to obtain reliable results of calculations of NMR spectroscopic parameters. Working with the solid-state model we found that any local variation of the TIZ crystalline structure is correlated with the variation of the values of the NMR parameters of each nucleus. The excellent agreement between experimental and calculated chemical shifts leads to the conclusion that the N10-H10 bond distance should be (1.00 ± 0.02) Å.

A reliable and efficient resonance theory based on analysis of DFT wave functions

Due to methodological difficulties and limitations of applicability, a quantitative bonding analysis based on the theory of resonance is presently not as convenient and popular as that based on the molecular orbital (MO) methods. Here, we propose an efficient quantitative resonance theory by expanding the DFT wave function in terms of a complete set of Lewis structures. By rigorously separating the resonance subsystem represented by a set of localized MOs, this approach is able to treat large molecules, nonplanar π-conjugate systems, and bonding systems mixing both σ and π electrons. Assessment in 2c-2e systems suggests a new projection-weighted symmetric orthogonalization method to evaluate the weights of resonance contributors, which overcomes the drawbacks of other weighting schemes. Applications to benzene, naphthalene and chlorobenzene show that the present method is insensitive to the basis set employed in the DFT calculations, and to the choices of the independent Lewis set determined by Rumer's rule. Advanced applications to diverse chemical problems provide unique and valuable insights into the understanding of hydrogen bonding, the π substituent effect on benzene, and the mechanism of Diels-Alder reactions.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

Clusters of Hydroxyl-Functionalized Cations Stabilized by Cooperative Hydrogen Bonds: The Role of Polarizability and Alkyl Chain Length

We explore quantum chemical calculations for studying clusters of hydroxyl-functionalized cations kinetically stabilized by hydrogen bonding despite strongly repulsive electrostatic forces. In a comprehensive study, we calculate clusters of ammonium, piperidinium, pyrrolidinium, imidazolium, pyridinium, and imidazolium cations, which are prominent constituents of ionic liquids. All cations are decorated with hydroxy-alkyl chains allowing H-bond formation between ions of like charge. The cluster topologies comprise linear and cyclic clusters up to the size of hexamers. The ring structures exhibit cooperative hydrogen bonds opposing the repulsive Coulomb forces and leading to kinetic stability of the clusters. We discuss the importance of hydrogen bonding and dispersion forces for the stability of the differently sized clusters. We find the largest clusters when hydrogen bonding is maximized in cyclic topologies and dispersion interaction is properly taken into account. The kinetic stability of the clusters with short-chained cations is studied for the different types of cations ranging from hard to polarizable or exhibiting additional functional groups such as the acidic C(2)-H position in the imidazolium-based cation. Increasing the alkyl chain length, the cation effect diminishes and the kinetic stability is exclusively governed by the alkyl chain tether increasing the distance between the positively charged rings of the cations. With adding the counterion tetrafluoroborate (BF₄⁻) to the cationic clusters, the binding energies immediately switch from strongly positive to strongly negative. In the neutral clusters, the OH functional groups of the cations can interact either with other cations or with the anions. The hexamer cluster with the cyclic H-bond motive and “released” anions is almost as stable as the hexamer built by H-bonded ion pairs exclusively, which is in accord with recent IR spectra of similar ionic liquids detecting both types of hydrogen bonding. For the cationic and neutral clusters, we discuss geometric and spectroscopic properties as sensitive probes of opposite- and like-charge interaction. Finally, we show that NMR proton chemical shifts and deuteron quadrupole coupling constants can be related to each other, allowing to predict properties which are not easily accessible by experiment.

Non-Bonded Radii of the Atoms Under Compression

We present quantum mechanical estimates for non-bonded, van der Waals-like, radii of 93 atoms in a pressure range from 0 to 300 gigapascal. Trends in radii are largely maintained under pressure, but atoms also change place in their relative size ordering. Multiple isobaric contractions of radii are predicted and are explained by pressure-induced changes to the electronic ground state configurations of the atoms. The presented radii are predictive of drastically different chemistry under high pressure and permit an extension of chemical thinking to different thermodynamic regimes. For example, they can aid in assignment of bonded and non-bonded contacts, for distinguishing molecular entities, and for estimating available space inside compressed materials. All data has been made available in an interactive web application.

Cooperativity in a cycloalkane-1,2/1,3-polyol corona: Topological hydrogen bonding in 1,2-diol motifs

A corona, consisting of 18 carbon atoms bearing 12 hydroxy groups in a continuous hydrogen-bonded chain, is built up by alternating degenerate conformations of alternating alkane-1,2-diol and 1,3-diol motifs. Geometries, proton nuclear magnetic resonance shifts and interaction energies for the dodecahydroxycyclo-octadecane and selected fragments are determined by density functional calculations at the B3LYP/6-311+G(d,p) level. Cooperative effects of O-H⋯O-H bonding are evident from the simple juxtaposition of these two motifs with a common OH group in butane-1,2,4-triol conformers. Bracketing a 1,2-diol motif with two 1,3-diol motifs in hexane-1,3,4,6-tetrol leads to a structure in which the 1,2-diol motif displays a bond critical point for hydrogen bonding. This is associated with enhancement of the shift of the hydrogen-bonded OH proton and of the corresponding H⋯O interaction energy. The full corona has a complete outer ring of O-H⋯O-H bond paths, and an inner ring of bond paths, due to C-H⋯H-C hydrogen-hydrogen bonding, which result in a central ring critical point. The topological O-H⋯O-H hydrogen bond, never seen in simple alkane-1,2-diols, is associated with cooperative enhancement of the H⋯O interaction energy, but this is not a necessary condition for a bond path: values for topological C-H⋯H-C hydrogen-hydrogen bonds can be as low as -0.4 kcal mol^-1 .

The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide

The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala₂₅ ). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala₂₅ . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala₂₅ . This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala₂₅ . Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.

Do weak interactions affect the biological behavior of DNA? A DFT study of CpG island-like chains

The origin, stability, and contribution to the formation of noncovalent interactions, such as hydrogen bonds and π - π stacking, have been already widely discussed. However, there are few discussions about the relevance of these weak interactions in DNA performance. In this work, we seek to shed light on the effect of hydrogen bonds and π - π stacking interactions on the biological behavior of DNA through the description of these intermolecular forces in CpG island-like (GC-rich) chains. Furthermore, we made some comparisons with TATA box-like (TA-rich) chains in order to describe hydrogen bond and π - π stacking interactions as a function of the DNA sequence. For hydrogen bonds, we found that there is not a significant effect related to the number of base pairs. Whereas for π - π stacking interactions, the energy tended to decrease as the number of base pairs increased. We observed anticooperative effects for both hydrogen bonds and π - π stacking interactions. These results are in contrast with those of TATA box-like chains since cooperative and additive effects were found for both hydrogen bonds and π - π stacking, respectively. Based on the chemical hardness and density of states, we can conclude that proteins may interact easier with GC-rich chains. We conclude that regardless of the chain length, a protein could interact more easily with these genomics regions because the π - π stacking energies did not increase as a function of the number of base pairs, making, for the first time, a first approximation of the influence of noncovalent interaction on DNA behavior. We did all this work by means of DFT framework included in the DMol³ code (M06-L/DNP). Graphical Abstract Cartoon representation of how nocovalent interactions affect the interaction of DNA with a protein, i.e., how hydrogen bond and π - π stacking interactions influence the biological behavior of DNA.

B3LYP, M06 and B3PW91 DFT assignment of nd8 metal-bis-(N-heterocyclic carbene) complexes

This paper is focused on the examination of the bonding properties of a series of [M(NHC)₂X₂] (M = nd⁸ transition metal; X = Cl, Br and I) complexes in normal, abnormal and mixed C^∩C coordination modes. Structures have been optimised in gas phase using B3LYP, M06 and P3BW91 functionals. Two basis sets have been used: the LanL2DZ and a mixed basis set (LanL2DZ for nd⁸ transition metals as well as halogen atoms and 6-311+G(d,p) for other atoms). Results obtained indicate that the B3PW91 bond distances are closer to experimental data. The complexation energies obtained for each binding mode increase in the order: Ni²⁺ < Pd²⁺ < Pt²⁺, independently of the halogen atom adopted. From the Quantum Theory of Atoms in Molecule (QTAIM) approach, the instability has been found to follow this trend: M - X < M - C. The analysis of metal-ligand interactions using the natural bond orbital (NBO) revealed that the strongest metal-ligand interactions are observed in the normal binding mode. The NCH → MX₂ donation terms were found to be interestingly predominant compared with back donation ones in the complexes studied, except in Pt chloride ones. The contribution of electrostatic interaction energy between the above fragments (∆E_elstat term) is in the range 57.48-63.95% traducing the fact that the interactions are mostly electrostatic. Graphical abstract.

Molecular Hydrogen as a Lewis Base in Hydrogen Bonds and Other Interactions

The second-order Møller–Plesset perturbation theory calculations with the aug-cc-pVTZ basis set were performed for complexes of molecular hydrogen. These complexes are connected by various types of interactions, the hydrogen bonds and halogen bonds are most often represented in the sample of species analysed; most interactions can be classified as σ-hole and π-hole bonds. Different theoretical approaches were applied to describe these interactions: Quantum Theory of ‘Atoms in Molecules’, Natural Bond Orbital method, or the decomposition of the energy of interaction. The energetic, geometrical, and topological parameters are analysed and spectroscopic properties are discussed. The stretching frequency of the H-H bond of molecular hydrogen involved in intermolecular interactions is considered as a parameter expressing the strength of interaction.

Microscopic origin of breakdown of Stokes-Einstein relation in binary mixtures: Inherent structure analysis

Aqueous binary mixtures often exhibit dramatic departure from the predicted hydrodynamic behavior when transport properties are plotted against composition. We show by inherent structure (IS) analysis that this sharp composition dependent breakdown of the Stokes-Einstein relation can be attributed to the non-monotonic variation in the average inherent structure energy of these mixtures. Further IS analysis reveals the existence of a unique ground state, stabilized by both the formation of an optimum number of H-bonds and a favorable hydrophobic interaction at this composition. The surprisingly sharp turnaround behavior observed in the effective hydrodynamic radius also owes its origin to the same combination of these two factors. Interestingly, the temperature dependence of isothermal compressibility shows a minimum at the particular composition. Extensive studies on water-dimethyl sulfoxide and water-ethanol mixtures using two different force-fields of water reveal many features that are nearly universal. A justification of this quasi-universal behavior is provided in terms of a mode-coupling theory (MCT) of viscosity, which can serve as the starting point of a remarkable correlation observed with the nearest neighbor structure, as captured by the first peaks of the radial distribution function, and the slowdown in the intermediate scattering function at intermediate wavenumbers. Therefore, the formation of the local structure captured through IS analysis can be correlated with the MCT.

Charge-Shift Bonding in Xenon Hydrides: An NBO/NRT Investigation on HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CCH, CN) via H-Xe Blue-Shift Phenomena

Noble-gas bonding represents curiosity. Some xenon hydrides, such as HXeY (Y = Cl, Br, I) and their hydrogen-bonded complexes HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH), have been identified in matrixes by observing H-Xe frequencies or its monomer-to-complex blue shifts. However, the H-Xe bonding in HXeY is not yet completely understood. Previous theoretical studies provide two answers. The first one holds that it is a classical covalent bond, based on a single ionic structure H-Xe+ Y-. The second one holds that it is resonance bonding between H-Xe+ Y- and H- Xe+-Y. This study investigates the H-Xe bonding, via unusual blue-shifted phenomena, combined with some NBO/NRT calculations for chosen hydrogen-bonded complexes HXeY···HX (Y = Cl, Br, I; X = OH, Cl, Br, I, CN, CCH). This study provides new insights into the H-Xe bonding in HXeY. The H-Xe bond in HXeY is not a classical covalent bond. It is a charge-shift (CS) bond, a new class of electron-pair bonds, which is proposed by Shaik and Hiberty et al. The unusual blue shift in studied hydrogen-bonded complexes is its H-Xe CS bonding character in IR spectroscopy. It is expected that these studies on the H-Xe bonding and its IR spectroscopic property might assist the chemical community in accepting this new-class electron-pair bond concept.

H/D Isotope Effects on 1H-NMR Chemical Shifts in Cyclic Heterodimers and Heterotrimers of Phosphinic and Phosphoric Acids

Hydrogen-bonded heterocomplexes formed by POOH-containing acids (diphenylphosphoric 1, dimethylphosphoric 2, diphenylphosphinic 3, and dimethylphosphinic 4) are studied by the low-temperature (100 K) ¹H-NMR and ³¹P-NMR using liquefied gases CDF₃/CDF₂Cl as a solvent. Formation of cyclic dimers and cyclic trimers consisting of molecules of two different acids is confirmed by the analysis of vicinal H/D isotope effects (changes in the bridging proton chemical shift, δH, after the deuteration of a neighboring H-bond). Acids 1 and 4 (or 1 and 3) form heterotrimers with very strong (short) H-bonds (δH ca. 17 ppm). While in the case of all heterotrimers the H-bonds are cyclically arranged head-to-tail, ···O=P-O-H···O=P-O-H···, and thus their cooperative coupling is expected, the signs of vicinal H/D isotope effects indicate an effective anticooperativity, presumably due to steric factors: when one of the H-bonds is elongated upon deuteration, the structure of the heterotrimer adjusts by shortening the neighboring hydrogen bonds. We also demonstrate the formation of cyclic tetramers: in the case of acids 1 and 4 the structure has alternating molecules of 1 and 4 in the cycle, while in case of acids 1 and 3 the cycle has two molecules of 1 followed by two molecules of 3.

Hydrogen Bonding in the Liquid State of Linear Alcohols: Molecular Dynamics and Thermodynamics

Linear monohydroxy alcohols are strongly hydrogen-bonded liquids that are considered to be homologues of water. Here, we report ab initio molecular dynamics simulations of the liquid alcohols, methanol to pentanol, and from the combined radial-angular probability distribution of the intermolecular O···O distances and HO···O angles determine the geometrical parameters that define the hydrogen bonds in these systems. The key feature of hydrogen bonds in the liquid alcohols, irrespective of the size of the alkyl group, is the strong orientation dependence with the donor-acceptor HO···O angle being close to zero, similar to that observed in liquid water. Hydrogen bond formation is consequently considered to be the passage from a state where donor-acceptor pairs show no preferred orientation to one where they are almost linear. The potential of mean force, the reversible work associated with this process, is computed from the pair probability density distributions obtained from the simulations and that for a hypothetical state where donor-acceptor pairs are randomly oriented. We find that the magnitude of the free energy for hydrogen bond formation is maximum for ethanol and show that this arises from a larger electrostatic contribution to hydrogen bond formation in ethanol as compared to the other alcohols.

Magnetic descriptors of hydrogen bonds in malonaldehyde and its derivatives

The nature of the hydrogen bond, HB, as such is still unknown, though a few of its most fundamental features has been uncovered during the last few decades. At the moment, it is possible to obtain reliable results for only a few of its broadest properties, like magnetic properties. They could give new insights into the physics underlying the strength and features of HBs. In this article we analyze the electronic origin of the NMR spectroscopic parameters of malonaldehyde, MA, and some substituted MAs. These substituted MAs are such that the H-bonds are assisted by one of two phenomena: resonance, RAHB, or charge, CAHB. We have studied the dependences of these parameters on two of the main factors which contribute the most to both phenomena, the geometrical and electronic factors, and found out how they can be used to characterize RAHB or CAHB by means of reliable theoretical calculations. We show that in the set of compounds analyzed here (i) the shielding of the proton of the H-bond can be used as a measure of the strength of the HB and (ii) the relation between the contact and non-contact mechanisms of J-couplings between donor and acceptor atoms is a reliable descriptor of whether the H-bond is resonance assisted or charge assisted.

Can 2-X-Ethanols Form Intramolecular Hydrogen Bonds?

For 2-X-ethanols, where X = F, OH, or NH₂, the gauche conformer is favored over the trans conformer by at least 2 kcal/mol. Initially, this preference, ΔE, was attributed to an intramolecular hydrogen bond, IMHB, between the OH and X groups. Over the years, this conclusion has been challenged by two major arguments. One claim is that the entirety of ΔE can be accounted for by the gauche effect. Against this, calculations using five different methods show that the maximum contribution of the gauche effect to ΔE is less than 1 kcal/mol. A second argument employs the quantum theory of atoms in molecules to contend that the absence of a bond critical point (BCP) between the OH and X groups in 2-X-ethanols denotes the lack of an IMHB. By looking at the 2-X-ethanols at fixed XCCO torsional angles ranging from 0° to 60°, it is shown that the BCP criterion is inconsistent with other properties such as energy, bond lengths, and stretching frequencies. These inconsistencies are removed when the theory of noncovalent interactions is used. The IMHBs in 2-X-ethanols are found to be similar in form but smaller in magnitude than their intermolecular counterparts. This work concludes that 2-X-ethanols form IMHBs.