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

The Nature of Hydrogen Bonds: A Delineation of the Role of Different Energy Components on Hydrogen Bond Strengths and Lengths

Hydrogen bonds are a complex interplay between different energy components, and their nature is still subject of an ongoing debate. In this minireview, we therefore provide an overview of the different perspectives on hydrogen bonding. This will be done by discussing the following individual energy components: 1) electrostatic interactions, 2) charge-transfer interactions, 3) π-resonance assistance, 4) steric repulsion, 5) cooperative effects, 6) dispersion interactions and 7) secondary electrostatic interactions. We demonstrate how these energetic factors are essential in a correct description of the hydrogen bond, and discuss several examples of systems whose energetic and geometrical features are not captured by easy-to-use predictive models.

Nine questions on energy decomposition analysis

The paper collects the answers of the authors to the following questions: Is the lack of precision in the definition of many chemical concepts one of the reasons for the coexistence of many partition schemes? Does the adoption of a given partition scheme imply a set of more precise definitions of the underlying chemical concepts? How can one use the results of a partition scheme to improve the clarity of definitions of concepts? Are partition schemes subject to scientific Darwinism? If so, what is the influence of a community's sociological pressure in the "natural selection" process? To what extent does/can/should investigated systems influence the choice of a particular partition scheme? Do we need more focused chemical validation of Energy Decomposition Analysis (EDA) methodology and descriptors/terms in general? Is there any interest in developing common benchmarks and test sets for cross-validation of methods? Is it possible to contemplate a unified partition scheme (let us call it the "standard model" of partitioning), that is proper for all applications in chemistry, in the foreseeable future or even in principle? In the end, science is about experiments and the real world. Can one, therefore, use any experiment or experimental data be used to favor one partition scheme over another? © 2019 Wiley Periodicals, Inc.

NBO 7.0: New vistas in localized and delocalized chemical bonding theory

We briefly outline some leading features of the newest version, NBO 7.0, of the natural bond orbital (NBO) wavefunction analysis program. Major extensions include: (1) a new NPEPA module implementing Karafiloglou's "polyelectron population analysis" in the NBO framework; (2) new RDM2 program infrastructure for describing electron correlation effects based on full evaluation of the second-order reduced density matrix; (3) improved convex-solver implementation of natural resonance theory (NRT), allowing a greatly expanded range of applications and associated "resonance NBO" (RNBO) visualization of chemical reactivity; (4) a variety of other improvements in well-established NBO algorithms. We also provide brief introduction to the new NBOPro@Jmol utility program, a plugin to the Jmol chemical structure viewer that serves as a convenient tool to provide on-demand NBO descriptors or orbital visualizations for a broad variety of chemical inquiries in research or classroom applications. © 2019 Wiley Periodicals, Inc.

What Is the Nature of Supramolecular Bonding? Comprehensive NBO/NRT Picture of Halogen and Pnicogen Bonding in RPH2···IF/FI Complexes (R = CH3, OH, CF3, CN, NO2)

We employ a variety of natural bond orbital (NBO) and natural resonance theory (NRT) tools to comprehensively investigate the nature of halogen and pnicogen bonding interactions in RPH₂···IF/FI binary complexes (R = CH₃, OH, CF₃, CN, and NO₂) and the tuning effects of R-substituents. Though such interactions are commonly attributed to “sigma-hole”-type electrostatic effects, we show that they exhibit profound similarities and analogies to the resonance-type 3-center, 4-electron (3c/4e) donor-acceptor interactions of hydrogen bonding, where classical-type “electrostatics” are known to play only a secondary modulating role. The general 3c/4e resonance perspective corresponds to a continuous range of interatomic A···B bond orders (b_AB), spanning both the stronger “covalent” interactions of the molecular domain (say, b_AB ≥ ½) and the weaker interactions (b_AB ˂ ½, often misleadingly termed “noncovalent”) that underlie supramolecular complexation phenomena. We show how a unified NBO/NRT-based description of hydrogen, halogen, pnicogen, and related bonding yields an improved predictive utility and intuitive understanding of empirical trends in binding energies, structural geometry, and other measurable properties that are expected to be manifested in all such supramolecular interaction phenomena.

An ab initio Study of the Structure and Energetics of Hydrogen Bonding in Ionic Liquids

Unlike typical hydrogen-bonded networks such as water, hydrogen bonded ionic liquids display some unusual characteristics due to the complex interplay of electrostatics, polarization, and dispersion forces in the bulk. Protic ionic liquids in particular contain close-to traditional linear hydrogen bonds that define their physicochemical properties. This work investigates whether hydrogen bonded ionic liquids (HBILs) can be differentiated from aprotic ionic liquids with no linear hydrogen bonds using state-of-the-art ab initio calculations. This is achieved through geometry optimizations of a series of single ion pairs of HBILs in the gas phase and an implicit solvent. Using benchmark CCSD(T)/CBS calculations, the electrostatic and dispersion components of the interaction energy of these systems are compared with those of aprotic ionic liquids. The inclusion of the implicit solvent significantly influenced geometries of single ion pairs, with the gas phase shortening the hydrogen bond to reduce electrostatic interactions. HBILs were found to have stronger interactions by at least 10EtMeNH0 kJ mol⁻¹ over aprotic ILs, clearly highlighting the electrostatic nature of hydrogen bonding. Geometric and energetic parameters were found to complement each other in determining the extent of hydrogen bonding present in these ionic liquids.

Correlations of NHN hydrogen bond energy with geometry and 1H NMR chemical shift difference of NH protons for aniline complexes

In this computational work, we propose to use the NMR chemical shift difference of NH₂ protons for 1:1 complexes formed by aniline and nitrogen-containing proton acceptors for the estimation of the hydrogen bond energy and geometry (N⋯H and N⋯N distances). The proposed correlations could be applied to other aromatic amines as well, in a gas phase, a solution, or a solid state, for both inter- and intramolecular hydrogen bonds. We considered a set of 21 complexes with the NHN hydrogen bond without proton transfer, including hydrogen bonds from weak to medium strong ones (2-21 kcal/mol), with neutral or anionic bases and with sp³ and sp² hybridized nitrogen proton acceptors. For each complex apart from direct hydrogen bond energy calculation, we have tested several other ways to estimate the energy: (a) using a correlation between NH stretching band intensity and hydrogen bond energy and (b) using correlations between electron density properties at (3, -1) bond critical point (quantum theory of atoms in molecules analysis) and hydrogen bond energy. Besides for the studied type of complexes, we obtained refined linear correlations linking the local electron kinetic (G) and potential (V) energy densities with the hydrogen bond energy.

Resonance Theory Reboot

What is now called "resonance theory" has a long and conflicted history. We first sketch the early roots of resonance theory, its heritage of diverse physics and chemistry conceptions, and its subsequent rise to reigning chemical bonding paradigm of the mid-20th century. We then outline the alternative "natural" pathway to localized Lewis- and resonance-structural conceptions that was initiated in the 1950s, given semi-empirical formulation in the 1970s, recast in ab initio form in the 1980s, and successfully generalized to multi-structural "natural resonance theory" (NRT) form in the 1990s. Although earlier numerical applications were often frustrated by the ineptness of then-available numerical solvers, the NRT variational problem was recently shown to be amenable to highly efficient convex programming methods that yield provably optimal resonance weightings at a small fraction of previous computational costs. Such convexity-based algorithms now allow a full "reboot" of NRT methodology for tackling a broad range of chemical applications, including the many familiar resonance phenomena of organic and biochemistry as well as the still broader range of resonance attraction effects in the inorganic domain. We illustrate these advances for prototype chemical applications, including (i) stable near-equilibrium species, where resonance mixing typically provides only small corrections to a dominant Lewis-structural picture, (ii) reactive transition-state species, where strong resonance mixing of reactant and product bonding patterns is inherent, (iii) coordinative and related supramolecular interactions of the inorganic domain, where sub-integer resonance bond orders are the essential origin of intermolecular attraction, and (iv) exotic long-bonding and metallic delocalization phenomena, where no single "parent" Lewis-structural pattern gains pre-eminent weighting in the overall resonance hybrid.

A-H…σ Hydrogen Bonds: Dihydrogen and Cycloalkanes as Proton Acceptors

ωB97XD/aug-cc-pVTZ calculations were performed for complexes of dihydrogen, cyclopropane, cyclobutane and cyclopentane, with simple proton donating species such as hydrogen fluoride, hydrogen chloride, water, hydrogen cyanide and acetylene. Numerous dependencies between geometrical, energetic and topological parameters of complexes considered were found, since various theoretical approaches were applied: Quantum Theory of 'Atoms in Molecules' (QTAIM), Natural Bond Orbital (NBO) method and energy decomposition analysis (EDA). It was confirmed that complexes of dihydrogen and cyclopropane are linked through the A-H…σ interactions that may be classified as hydrogen bonds. In the case of complexes of cyclobutane such hydrogen bonds are rather weak. Other type and also weak A-H…C hydrogen bonds are formed for complexes with cyclopentane.

A theoretical investigation on Cu/Ag/Au bonding in XH2P⋯MY(X = H, CH3, F, CN, NO2; M = Cu, Ag, Au; Y = F, Cl, Br, I) complexes

Intermolecular interaction of XH₂P···MY (X = H, CH₃, F, CN, NO₂; M = Cu, Ag, Au; Y = F, Cl, Br, I) complexes was investigated by means of an ab initio method. The molecular interaction energies are in the order Ag < Cu < Au and increased with the decrease of R_P···M. Interaction energies are strengthened when electron-donating substituents X connected to XH₂P, while electron-withdrawing substituents produce the opposite effect. The strongest P···M bond was found in CH₃H₂P···AuF with -70.95 kcal/mol, while the weakest one was found in NO₂H₂P···AgI with -20.45 kcal/mol. The three-center/four-electron (3c/4e) resonance-type of P:-M-:Y hyperbond was recognized by the natural resonance theory and the natural bond orbital analysis. The competition of P:M-Y ↔ P-M:Y resonance structures mainly arises from hyperconjugation interactions; the bond order of b_P-M and b_M-Y is in line with the conservation of the idealized relationship b_P-M + b_M-Y ≈ 1. In all MF-containing complexes, P-M:F resonance accounted for a larger proportion which leads to the covalent characters for partial ionicity of MF. The interaction energies of these Cu/Ag/Au complexes are basically above the characteristic values of the halogen-bond complexes and close to the observed strong hydrogen bonds in ionic hydrogen-bonded species.

Noncovalent Interactions between Molecular Hydrogen and the Alkali Fluorides: H-H···F-M (M = Li, Na, K, Rb, Cs). High Level Theoretical Predictions and SAPT Analysis

Various types of hydrogen bonds have been recognized during the past century. In this research, a new type of noncovalent interaction, the dipole-induced hydrogen bond formed between a hydrogen molecule and an alkali halide, H-H···F-M, is studied. Proposed by Zhang and co-workers ( Phys. Chem. Chem. Phys. 2015, 17, 20361), these systems are extensively investigated initially using the "gold standard" CCSD(T) method in conjunction with augmented correlation-consistent polarized core-valence basis sets up to quadruple-ζ. The full triple excitations CCSDT method has been used to further refine the energies. Several properties including geometries, bond energies, vibrarional frequencies, charge distributions, and dipole moments have been reported. The earlier Zhang research considered only the linear H-H···F-M structures. However, we find these linear stationary points to be separated by very small barriers from the much lower lying bent C _s structures. The CCSDT/aug-cc-pCVQZ(-PP) method predicts the dissociation energies for bent H-H···F-M (M = Li, Na, K, Rb, Cs) are 2.76, 2.96, 3.00, 2.89, and 2.49 kcal mol^-1, respectively, suggesting that the H···F hydrogen bond becomes gradually stronger when alkali metal M goes down the periodic table from Li to K but becomes slightly weaker for Rb and even more for Cs. This Li < Na < K > Rb > Cs order is consistent with that for the dipole moments for the isolated MF (M = Li, Na, K, Rb, Cs) diatomics. Symmetry adapted perturbation theory (SAPT) is used to understand these unusual noncovalent interactions.

Using the localized pair model to describe the nature of the hydrogen bond

Herein, we present a novel study of hydrogen bonding using the localized pair model. Using localized molecular orbitals, we examine how the two localized electron pairs involved in hydrogen bonding change upon the formation of the interaction, as well as while the distance between the two species within the complex change. This is achieved through the use of extracule densities, which describe the probability associated with the centre of mass of an electron pair. Quantitative assessments are performed by analyzing the extracule density demonstrating that this method is on par with other density-based tools for the analysis of hydrogen bonds and their strengths.

A new strategy to efficiently cleave and form C-H bonds using proton-coupled electron transfer

Oxidative activation and reductive formation of C-H bonds are crucial in many chemical, industrial, and biological processes. Reported here is a new strategy for these transformations, using a form of proton-coupled electron transfer (PCET): intermolecular electron transfer coupled to intramolecular proton transfer with an appropriately placed cofactor. In a fluorenyl-benzoate, the positioned carboxylate facilitates rapid cleavage of a benzylic C-H bond upon reaction with even weak 1e- oxidants, for example, decamethylferrocenium. Mechanistic studies establish that the proton and electron transfer to disparate sites in a single concerted kinetic step, via multi-site concerted proton-electron transfer. This work represents a new elementary reaction step available to C-H bonds. This strategy is extended to reductive formation of C-H bonds in two systems. Molecular design considerations and possible utility in synthetic and enzymatic systems are discussed.

Understanding and modulating the high-energy properties of noble-gas hydrides from their long-bonding: an NBO/NRT investigation on HNgCO+/CS+/OSi+ and HNgCN/NC (Ng = He, Ar, Kr, Xe, Rn) molecules

The noble-gas hydrides, HNgX (X is an electronegative atom or fragment), represent potential high-energy materials because their two-body decomposition process, HNgX → Ng + HX, is strongly exoergic. Our previous studies have shown that each member of the HNgX (X = halogen atom or CN/NC fragment) molecules is composed of three leading resonance structures: two ω-bonding structures (H-Ng+ :X- and H:- Ng+-X) and one long-bonding structure (H∧X). The last one paints a novel [small sigma, Greek, circumflex]-type long-bonding picture. The present study focuses on the relationship between this novel bonding motif and the unusual energetic properties. We chose HNgCO+/CS+/OSi+/CN/NC, with the formula HNgAB (Ng = He, Ar, Kr, Xe, Rn; AB = CO+/CS+/OSi+/CN/NC) as the research system. We first investigated the bonding of HNgCO+ and its analogous HNgCS+/OSi+ species using NBO/NRT methods, and quantitatively compared the bonding with that in HNgCN/NC molecules. NBO/NRT results showed that each of the HNgCO+/CS+/OSi+ molecules could be better represented as a resonance hybrid of ω-bonding and long-bonding structures, but the long-bonding is much weaker than that in HNgCN/NC molecules. Furthermore, we introduced the long-bonding concept into the rationalization of the high-energy properties, and found a good correlation between the highly exothermic two-body dissociation channel and the long-bond order, bH-A. We also found that the long-bond order is highly tunable for these noble-gas hydrides due to its dependence on the nature of the electronegative AB fragments or the central noble-gas atoms, Ng. On the basis of these results, we could optimize the energetic properties by changing the long-bonding motif of our studied molecules. Overall, this study shows that the long-bonding model provides an easy way to rationalize and modulate the unusual energy properties of noble-gas hydrides, and that it is helpful to predict some noble-gas hydrides as potential energetic materials.

On the Computational Characterization of Charge-Transfer Effects in Noncovalently Bound Molecular Complexes

Charge-transfer (CT) is an important binding force in the formation of intermolecular complexes, and there have been a variety of theoretical models proposed to quantify this effect. These approaches, which typically rely on a definition of a "CT-free" state based on a partition of the system, sometimes yield significantly different results for a given intermolecular complex. Two widely used definitions of the "CT-free" state, the absolutely localized molecular orbitals (ALMO) method (where only on-fragment orbital mixings are permitted) and the constrained density functional theory (CDFT) approach (where fragment electron populations are fixed), are carefully examined in this work. Natural bond orbital (NBO) and the regularized symmetry-adapted perturbation theory (SAPT) are also briefly considered. Results for the ALMO and CDFT definitions of CT are compared on a broad range of model systems, including hydrogen-bonding systems, borane complexes, metal-carbonyl complexes, and complexes formed by water and metal cations. For most of these systems, CDFT yields a much smaller equilibrium CT energy compared to that given by the ALMO-based definition. This is mainly because the CDFT population constraint does not fully inhibit CT, which means that the CDFT "CT-free" state is in fact CT-contaminated. Examples of this contamination include (i) matching forward and backward donation (e.g., formic acid dimer) and (ii) unidirectional CT without changing fragment populations. The magnitude of the latter effect is quantified in systems such as the water dimer by employing a 3-space density constraint in addition to the orbital constraint. Furthermore, by means of the adiabatic EDA, it is shown that several observable effects of CT, such as the "pyramidalization" of the planar BH₃ molecule upon the complexation with Lewis bases, already appear on the "CT-free" CDFT surface. These results reveal the essential distinctions between the ALMO and CDFT definitions of CT and suggest that the former is more consistent with accepted understanding of the role of CT in intermolecular binding.

Refinement of labile hydrogen positions based on DFT calculations of 1H NMR chemical shifts: comparison with X-ray and neutron diffraction methods

Numerous gas phase electron diffraction, ultra-fast electron diffraction, X-ray and neutron diffraction experiments on β-dicarbonyl compounds exhibiting enol-enol tautomeric equilibrium, with emphasis on acetylacetone and dibenzoylmethane, have so far been reported with conflicting results on the structural details of the O-HO intramolecular hydrogen bond and resulted in alternative hypotheses on the intramolecular hydrogen bond potential function either a double minimum potential corresponding to two tautomeric forms in equilibrium or a single symmetrical one. We demonstrate herein, firstly, that the DFT calculated OH ¹H NMR chemical shifts of acetylacetone and dibenzoylmethane exhibit a strong linear dependence on the computed OO hydrogen bond length of ∼-50 ppm Å^-1 and as a function of the O-HO bond angle of ∼1 ppm per degree, upon the transfer of the hydrogen atom from the ground state toward the transition state. Secondly, the refinement of labile hydrogen atomic positions in intramolecular hydrogen bonds based on the root-mean-square deviation between experimentally determined and DFT calculated ¹H NMR chemical shifts in solution can provide high resolution structures of O-H and O(H)O bond lengths and O-HO bond angles with an accuracy of ∼10^-2 Å and ∼0.5°, respectively. Thirdly, the calculated ¹H NMR chemical shifts in solution of the two ground state tautomers in equilibrium of acetylacetone and dibenzoylmethane are in excellent agreement with the experimental value, even for moderate basis sets for energy minimization. In contrast, the single symmetrical structure in a strongly delocalized system is a transition state with calculated ¹H NMR chemical shifts which strongly deviate from the experimental value. Fourth, the DFT calculated ground state O-H bond lengths of acetylacetone and dibenzoylmethane are in quantitative agreement with the literature data which take into account the effect of quantum nuclear motion. The DFT structural results are critically discussed with respect to the state-of-the-art variable temperature X-ray and neutron diffraction methods.

Energy decomposition analysis in an adiabatic picture

Energy decomposition analysis (EDA) of electronic structure calculations has facilitated quantitative understanding of diverse intermolecular interactions. Nevertheless, such analyses are usually performed at a single geometry and thus decompose a "single-point" interaction energy. As a result, the influence of the physically meaningful EDA components on the molecular structure and other properties are not directly obtained. To address this gap, the absolutely localized molecular orbital (ALMO)-EDA is reformulated in an adiabatic picture, where the frozen, polarization, and charge transfer energy contributions are defined as energy differences between the stationary points on different potential energy surfaces (PESs), which are accessed by geometry optimizations at the frozen, polarized and fully relaxed levels of density functional theory (DFT). Other molecular properties such as vibrational frequencies can thus be obtained at the stationary points on each PES. We apply the adiabatic ALMO-EDA to different configurations of the water dimer, the water-Cl^- and water-Mg²⁺/Ca²⁺ complexes, metallocenes (Fe²⁺, Ni²⁺, Cu²⁺, Zn²⁺), and the ammonia-borane complex. This method appears to be very useful for unraveling how physical effects such as polarization and charge transfer modulate changes in molecular properties induced by intermolecular interactions. As an example of the insight obtained, we find that a linear hydrogen bond geometry for the water dimer is preferred even without the presence of polarization and charge transfer, while the red shift in the OH stretch frequency is primarily a charge transfer effect; by contrast, a near-linear geometry for the water-chloride hydrogen bond is achieved only when charge transfer is allowed.

Theoretical Prediction of Robust Second-Row Oxyanion Clusters in the Metastable Domain of Antielectrostatic Hydrogen Bonding

We provide ab initio and density functional theory evidence for a family of surprisingly robust like-charged clusters of common HSO₄^- and H₂PO₄^- oxyanions, ranging up to tetramers of net charge 4-. Our results support other recent theoretical and experimental evidence for "antielectrostatic" hydrogen-bonded (AEHB) species that challenge conventional electrostatic conceptions and force-field modeling of closed-shell ion interactions. We provide structural and energetic descriptors of the predicted kinetic well-depths (in the range 3-10 kcal/mol) and barrier widths (in the range 2-4 Å) for simple AEHB dimers, including evidence of extremely strong hydrogen bonding in the fluoride-bisulfate dianion. For more complex polyanionic species, we employ natural-bond-orbital-based descriptors to characterize the electronic features of the cooperative hydrogen-bonding network that are able to successfully defy Coulomb explosion. The computational results suggest a variety of kinetically stable AEHB species that may be suitable for experimental detection as long-lived gas-phase species or structural units of condensed phases, despite the imposing electrostatic barriers that oppose their formation under ambient conditions.

Revealing Electron-Electron Interactions within Lewis Pairs in Chemical Systems

The so-called "Lewis pair" is a ubiquitous phenomenon in chemistry and is often used as an intuitive construct to predict and rationalize chemical structure and behavior. Concepts from the very general Valence Shell Electron Pair Repulsion (VSEPR) model to the most esoteric reaction mechanism routinely rely on the notion that electrons tend to exist in pairs and that these pairs can be thought of as being localized to a particular region of space. It is precisely this localization that allows one to intuit how these pairs might behave, generally speaking, so that reasonable predictions may be made regarding molecular structure, intermolecular interactions, property trends, and reaction mechanisms, etc. Of course, it is rather unfortunate that the Lewis model is entirely qualitative and yields no information regarding how any specific electron pair is distributed. Here we demonstrate a novel electronic structure analysis technique that predicts and analyzes precise quantitative details about the relative and absolute distribution of individual electron pairs. This Single Electron Pair Distribution Analysis (SEPDA) reveals quantitative details about the distribution of the well-known Lewis pairs, such as how they are distributed in space and how their relative velocities change in various chemical contexts. We show that these distributions allow one to image the explicitly pairwise electronic behavior of bonds and lone pairs. We further demonstrate how this electronic behavior changes with several conditions to explore the nature of the covalent chemical bond, non-covalent interactions, bond formation, and exotic 3-center-2-electron species. It is shown that indications of the strength of bonded and non-bonded interactions may also be gleaned from such distributions and SEPDA can be used as a tool to differentiate between interaction types. We anticipate that SEPDA will be of broad utility in a wide variety of chemical contexts because it affords a very detailed, visual and intuitive analysis technique that is generally applicable.

Hydrogen bonds and other interactions as a response to protect doublet/octet electron structure

MP2/aug-cc-pVTZ calculations were performed for complexes linked by hydrogen bonds. Three types of proton donating species were taken into account: H₂O, CCl₃H, and H₃O⁺. These calculations are supported by the natural bond orbital (NBO) method and the quantum theory of atoms in molecules (QTAIM) approach. Numerous correlations between parameters of H-bonded systems were found. The most important are those which show the response of the system on the H-bond formation; for example, the increase of polarization of the A-H bond correlates with the strength of the hydrogen bond. Similar relationships were found for the σ-hole bonds while the π-hole bonds do not follow the trends known for the hydrogen bonds. Graphical abstract Hydrogen bonds and other interactions as a response to protect doublet/octet electron structureᅟ.

Characterizing the interplay of Pauli repulsion, electrostatics, dispersion and charge transfer in halogen bonding with energy decomposition analysis

Variational energy decomposition analysis establishes charge-transfer as the origin of halogen bond strength differences that go against electrostatics.

Polyion Covalency: Exotic Species from the Unexplored World of Electrostatically Shielded Molecular Ion Chemistry

Standard quantum chemical methods have been employed to describe a variety of kinetically stable polyionic molecular species that are trapped in appreciable potential wells by chemical bonding forces, despite powerful electrostatic opposition that challenges conventional chemical detection and characterization. The studied species are covalent or dative analogs of "anti-electrostatic" hydrogen-bonded (AEHB) species, all illustrating how short-range quantum covalency can overcome the powerful "shielding" opposition of long-range electrostatic forces to form highly charged molecular species, analogous to known neutral or singly ionic counterparts. Computational predictions of representative structural, spectroscopic, and NBO-based electronic signatures of multiply charged analogs of common neutral species (CH₃ CH₃ , CO₂ , FeCO) are provided to suggest the unique material properties characteristic of this shielded domain of polyionic chemical phenomena.

Predicting the Ionic Product of Water

We present a first-principles calculation and mechanistic characterization of the ion product of liquid water (K _W ), based on Quantum Cluster Equilibrium (QCE) theory with a variety of ab initio and density functional methods. The QCE method is based on T-dependent Boltzmann weighting of different-sized clusters and consequently enables the observation of thermodynamically less favored and therefore low populated species such as hydronium and hydroxide ions in water. We find that common quantum chemical methods achieve semi-quantitative accuracy in predicting K _W and its T-dependence. Dominant ion-pair water clusters of the QCE equilibrium distribution are found to exhibit stable 2-coordinate buttress-type motifs, all with maximally Grotthus-ordered H-bond patterns that successfully prevent recombination of hydronium and hydroxide ions at 3-coordinate bridgehead sites. We employ standard quantum chemistry techniques to describe kinetic and mechanistic aspects of ion-pair formation, and we obtain NBO-based bonding indices to characterize other electronic, structural, spectroscopic, and reactive properties of cluster-mediated ionic dissociation.

Hydrogen-Bond Strength of CC and GG Pairs Determined by Steric Repulsion: Electrostatics and Charge Transfer Overruled

Theoretical and experimental studies have elucidated the bonding mechanism in hydrogen bonds as an electrostatic interaction, which also exhibits considerable stabilization by charge transfer, polarization, and dispersion interactions. Therefore, these components have been used to rationalize the differences in strength of hydrogen‐bonded systems. A completely new viewpoint is presented, in which the Pauli (steric) repulsion controls the mechanism of hydrogen bonding. Quantum chemical computations on the mismatched DNA base pairs CC and GG (C=cytosine, G=guanine) show that the enhanced stabilization and shorter distance of GG is determined entirely by the difference in the Pauli repulsion, which is significantly less repulsive for GG than for CC. This is the first time that evidence is presented for the Pauli repulsion as decisive factor in relative hydrogen‐bond strengths and lengths.

NMR and IR Investigations of Strong Intramolecular Hydrogen Bonds

For the purpose of this review, strong hydrogen bonds have been defined on the basis of experimental data, such as OH stretching wavenumbers, νOH, and OH chemical shifts, δOH (in the latter case, after correction for ring current effects). Limits for O-H···Y systems are taken as 2800 > νOH > 1800 cm-1, and 19 ppm > δOH > 15 ppm. Recent results as well as an account of theoretical advances are presented for a series of important classes of compounds such as β-diketone enols, β-thioxoketone enols, Mannich bases, proton sponges, quinoline N-oxides and diacid anions. The O···O distance has long been used as a parameter for hydrogen bond strength in O-H···O systems. On a broad scale, a correlation between OH stretching wavenumbers and O···O distances is observed, as demonstrated experimentally as well as theoretically, but for substituted β-diketone enols this correlation is relatively weak.

Hydrogen Atomic Positions of O-H···O Hydrogen Bonds in Solution and in the Solid State: The Synergy of Quantum Chemical Calculations with ¹H-NMR Chemical Shifts and X-ray Diffraction Methods

The exact knowledge of hydrogen atomic positions of O-H···O hydrogen bonds in solution and in the solid state has been a major challenge in structural and physical organic chemistry. The objective of this review article is to summarize recent developments in the refinement of labile hydrogen positions with the use of: (i) density functional theory (DFT) calculations after a structure has been determined by X-ray from single crystals or from powders; (ii) ¹H-NMR chemical shifts as constraints in DFT calculations, and (iii) use of root-mean-square deviation between experimentally determined and DFT calculated ¹H-NMR chemical shifts considering the great sensitivity of ¹H-NMR shielding to hydrogen bonding properties.

Quantitative probing of subtle interactions among H-bonds in alpha hydroxy carboxylic acid complexes

The alpha OH stretching frequency may be affected upon complexing with water and formic acid.

Dissection of H-bonding interactions in a glycolic acid–water dimer

The binding strength and collective effects of multiple H-bonds in the glycolic acid–water dimer were studied in comparison to the aromatic analog, 9-hydroxy-9-fluorene carboxylic acid (9HFCA).

Unravelling hydrogen bonding interactions of tryptamine–water dimer from neutral to cation

The physical origin for the three intermolecular hydrogen bonds in the neutral and cationic forms of the tryptamine–water dimer is explored.

Toward a uniform description of hydrogen bonds and halogen bonds: correlations of interaction energies with various geometric, electronic and topological parameters

Correlations between interaction energies and various structural parameters were established to reveal the differences between hydrogen bonds and halogen bonds.

Some Brief Notes on Theoretical and Experimental Investigations of Intramolecular Hydrogen Bonding

A review of selected literature data related to intramolecular hydrogen bonding in ortho-hydroxyaryl Schiff bases, ortho-hydroxyaryl ketones, ortho-hydroxyaryl amides, proton sponges and ortho-hydroxyaryl Mannich bases is presented. The paper reports on the application of experimental spectroscopic measurements (IR and NMR) and quantum-mechanical calculations for investigations of the proton transfer processes, the potential energy curves, tautomeric equilibrium, aromaticity etc. Finally, the equilibrium between the intra- and inter-molecular hydrogen bonds in amides is discussed.

Assessment of the Presence and Strength of H-Bonds by Means of Corrected NMR

The downfield shift of the NMR signal of the bridging proton in a H-bond (HB) is composed of two elements. The formation of the HB causes charge transfer and polarization that lead to a deshielding. A second factor is the mere presence of the proton-accepting group, whose electron density and response to an external magnetic field induce effects at the position of the bridging proton, exclusive of any H-bonding phenomenon. This second positional shielding must be subtracted from the full observed shift in order to assess the deshielding of the proton caused purely by HB formation. This concept is applied to a number of H-bonded systems, both intramolecular and intermolecular. When the positional shielding is removed, the remaining chemical shift is in much better coincidence with other measures of HB strength.

Hematite(001)-liquid water interface from hybrid density functional-based molecular dynamics

The atom-scale characterisation of interfaces between transition metal oxides and liquid water is fundamental to our mechanistic understanding of diverse phenomena ranging from crystal growth to biogeochemical transformations to solar fuel production. Here we report on the results of large-scale hybrid density functional theory-based molecular dynamics simulations for the hematite(001)-liquid water interface. A specific focus is placed on understanding how different terminations of the same surface influence surface solvation. We find that the two dominant terminations for the hematite(001) surface exhibit strong differences both in terms of the active species formed on the surface and the strength of surface solvation. According to present simulations, we find that charged oxyanions (-O(-)) and doubly protonated oxygens (-OH[Formula: see text]) can be formed on the iron terminated layer via autoionization of neutral -OH groups. No such charged species are found for the oxygen terminated surface. In addition, the missing iron sublayer in the iron terminated surface strongly influences the solvation structure, which becomes less well ordered in the vicinity of the interface. These pronounced differences are likely to affect the reactivity of the two surface terminations, and in particular the energetics of excess charge carriers at the surface.

Lone-pair-π interactions: analysis of the physical origin and biological implications

Lone-pair-π (lp-π) interactions have been suggested to stabilize DNA and protein structures, and to participate in the formation of DNA-protein complexes. To elucidate their physical origin, we have carried out a theoretical multi-approach analysis of two biologically relevant model systems, water-indole and water-uracil complexes, which we compared with the structurally similar chloride-tetracyanobenzene (TCB) complex previously shown to contain a strong charge-transfer (CT) binding component. We demonstrate that the CT component in lp-π interactions between water and indole/uracil is significantly smaller than that stabilizing the Cl(-)-TCB reference system. The strong lp(Cl(-))-π(TCB) orbital interaction is characterized by a small energy gap and an efficient lp-π* overlap. In contrast, in lp-π interactions between water and indole or uracil, the corresponding energy gap is larger and the overlap less efficient. As a result, water-uracil and water-indole interactions are weak forces composed by smaller contributions from all energy components: electrostatics, polarization, dispersion, and charge transfer. In addition, indole exhibits a negative electrostatic potential at its π-face, making lp-π interactions less favorable than O-Hπ hydrogen bonding. Consequently, some of the water-tryptophan contacts observed in X-ray structures of proteins and previously interpreted as lp-π interactions [Luisi, et al., Proteins, 2004, 57, 1-8], might in fact arise from O-Hπ hydrogen bonding.