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).

Complexes of phosphine oxides with substituted phenols: hydrogen bond characterization based on shifts of PO stretching bands

In this work IR spectral characteristics of PO groups are used to evaluate the strength of OHO hydrogen bonds. Three phosphine oxides: triphenylphosphine oxide, tributylphosphine oxide and hexamethylphosphoramide are investigated as proton acceptors. The results of the experimental IR study and DFT calculation of 30 complexes formed by phosphine oxides with various substituted phenols or CF3CH2OH in CCl4 solution at room temperature are reported. We show that the PO vibrational frequency changes non-linearly upon hydrogen bond formation and strengthening and that the shift of the PO band could be used for the estimation of hydrogen bond strength in complexes with phosphine oxides. The accuracy of these estimations and the influence of solvation effects on the main characteristics of complexes are discussed.

Coexistence of Cationic and Anionic Phosphate Moieties in Solids: Unusual but Not Impossible

Phosphoric acid is commonly known either as a neutral molecule or as an anion (phosphate). We theoretically confirm by ab initio molecular dynamics simulations (AIMD) that a cationic form H4PO4+ coexists with the anionic form H2PO4- in the same salt. This paradoxical situation is achieved by partial substitution of Cs+ by H4PO4+ in CsH2PO4. Thus, HnPO4 acts simultaneously as both the positive and the negative ion of the salt. We analyze the dynamical protonation pattern within the unusual hydrogen bond network that is established between the ions. Our AIMD simulations show that a conventional assignment of protonation states of the phosphate groups is not meaningful. Instead, a better description of the protonation situation is achieved by an efficiently fractional assignment of the strongly hydrogen-bonded protons to both its nearest and next-nearest oxygen neighbors.

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

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

Acid Dissociation in (HX)n (H2 O)n Clusters (X=F, Cl, Br, I; n=2, 3)

In this work, we analyze the interactions between two or three hydrogen halide molecules and the same number of water moieties through a systematic exploration of their potential energy surfaces. Our results indicate that the most stable HF and HCl aggregates do not experience dissociation of any of the acid fragments, even with three water molecules. In contrast, in the HBr and HI clusters, one of the acid fragments does dissociate. While the global minimum of (HBr)₃ (H₂ O)₃ is a hydrogen-bridged bihalide anion (BrHBr^- ), which is persistent at temperatures up to 203 K, the lowest energy structure of (HI)₃ (H₂ O)₃ has a separated ion pair, but the motif with a bihalide anion (IHI^- ) is only 0.2 kcal mol^-1 above the global minimum. Among the more stable structures is a broad spectrum of contacts, including water⋯water, HX⋯water, and HX⋯HX hydrogen bonds, halogen bonds, ionic and long-range X⋯H contacts.

Unusually Large Effects of Charge-assisted C-H⋅⋅⋅F Hydrogen Bonds to Anionic Fluorine in Organic Solvents: Computational Study of 19 F NMR Shifts versus Thermochemistry

A comparison of computed ¹⁹ F NMR chemical shifts and experiment provides evidence for large specific solvent effects for fluoride-type anions interacting with the σ*(C-H) orbitals in organic solvents like MeCN or CH₂ Cl₂ . We show this for systems ranging from the fluoride ion and the bifluoride ion [FHF]^- to polyhalogen anions [ClF_x ]^- . Discrepancies between computed and experimental shifts when using continuum solvent models like COSMO or force-field-based descriptions like the 3D-RISM-SCF model show specific orbital interactions that require a quantum-mechanical treatment of the solvent molecules. This is confirmed by orbital analyses of the shielding constants, while less negatively charged fluorine atoms (e. g., in [EF₄ ]^- ) do not require such quantum-mechanical treatments to achieve reasonable accuracy. The larger ¹⁹ F solvent shift of fluoride in MeCN compared to water is due to the larger coordination number in the former. These observations are due to unusually strong charge-assisted C-H⋅⋅⋅F^- hydrogen bonds, which manifest beyond some threshold negative natural charge on fluorine of ca. < -0.6 e. The interactions are accompanied by sizable free energies of solvation, in the order F^- ≫[FHF]^- >[ClF₂ ]^- >[ClF₄ ]^- . COSMO-RS solvation free energies tend to moderately underestimate those from the micro-solvated cluster treatment. Red-shifted and intense vibrational C-H stretching bands, potentially accessible in bulk solution, are further spectroscopic finger prints.

Hydrogen and Lithium Bonds-Lewis Acid Units Possessing Multi-Center Covalent Bonds

MP2/aug-cc-pVTZ calculations were carried out on complexes wherein the proton or the lithium cation is located between π-electron systems, or between π-electron and σ-electron units. The acetylene or its fluorine and lithium derivatives act as the Lewis base π-electron species similarly to molecular hydrogen, which acts as the electron donor via its σ-electrons. These complexes may be classified as linked by π-H∙∙∙π/σ hydrogen bonds and π-Li∙∙∙π/σ lithium bonds. The properties of these interactions are discussed, and particularly the Lewis acid units are analyzed, because multi-center π-H or π-Li covalent bonds may occur in these systems. Various theoretical approaches were applied here to analyze the above-mentioned interactions—the Quantum Theory of Atoms in Molecules (QTAIM), the Symmetry-Adapted Perturbation Theory (SAPT) and the Non-Covalent Interaction (NCI) method.

Classification of So-Called Non-Covalent Interactions Based on VSEPR Model

The variety of interactions have been analyzed in numerous studies. They are often compared with the hydrogen bond that is crucial in numerous chemical and biological processes. One can mention such interactions as the halogen bond, pnicogen bond, and others that may be classified as σ-hole bonds. However, not only σ-holes may act as Lewis acid centers. Numerous species are characterized by the occurrence of π-holes, which also may play a role of the electron acceptor. The situation is complicated since numerous interactions, such as the pnicogen bond or the chalcogen bond, for example, may be classified as a σ-hole bond or π-hole bond; it ultimately depends on the configuration at the Lewis acid centre. The disadvantage of classifications of interactions is also connected with their names, derived from the names of groups such as halogen and tetrel bonds or from single elements such as hydrogen and carbon bonds. The chaos is aggravated by the properties of elements. For example, a hydrogen atom can act as the Lewis acid or as the Lewis base site if it is positively or negatively charged, respectively. Hence names of the corresponding interactions occur in literature, namely hydrogen bonds and hydride bonds. There are other numerous disadvantages connected with classifications and names of interactions; these are discussed in this study. Several studies show that the majority of interactions are ruled by the same mechanisms related to the electron charge shifts, and that the occurrence of numerous interactions leads to specific changes in geometries of interacting species. These changes follow the rules of the valence-shell electron-pair repulsion model (VSEPR). That is why the simple classification of interactions based on VSEPR is proposed here. This classification is still open since numerous processes and interactions not discussed in this study may be included within it.

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

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

Pseudo Jahn-Teller Origin of the Proton-transfer Energy Barrier in the Hydrogen-bonded [FHF]-System

The results of ab initio calculations of the adiabatic potential energy surfaces for the proton-bound [FHF]- system at different F-F distances have been rationalized in the framework of the vibronic theory. It is shown that the instability of the symmetric D∞h structure at increased F∙∙∙F distances and the proton displacement to one of the fluorine atoms is due to the pseudo Jahn–Teller mixing of the ground 1Σg electronic state with the lowest excited state of 1Σu symmetry through the asymmetric σu vibrational mode.

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

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

Intramolecular Hydrogen Bond Energy and Its Decomposition—O–H∙∙∙O Interactions

The method to calculate the energy of intramolecular hydrogen bond is proposed and tested for a sample of malonaldehyde and its fluorine derivatives; the corresponding calculations were performed at the ωB97XD/aug-cc-pVTZ level. This method based on relationships found for related intermolecular hydrogen bonds is compared with other approaches which may be applied to estimate the intramolecular hydrogen bond energy. Particularly, methods based on the comparison of the system that contains the intramolecular hydrogen bond compared with corresponding conformations where such interaction does not occur are discussed. The function-based energy decomposition analysis, FB-EDA, of the intramolecular hydrogen bonds is also proposed here.

CP2K: An electronic structure and molecular dynamics software package - Quickstep: Efficient and accurate electronic structure calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially aimed at massively parallel and linear-scaling electronic structure methods and state-of-the-art ab initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.

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

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

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

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

Elucidating Solvent Effects on Strong Intramolecular Hydrogen Bond: DFT-MD Study of Dibenzoylmethane in Methanol Solution

The dynamic aspect of solvation plays a crucial role in determining properties of strong intramolecular hydrogen bonds since solvent fluctuations modify instantaneous hydrogen-bonded proton transfer barriers. Previous studies pointed out that solvent-solute interactions in the first solvation shell govern the position of the proton but the ability of the electric field due to other solvent molecules to localize the proton remains an important issue. In this work, we examine the structure of the O-H⋅⋅⋅O intramolecular hydrogen bond of dibenzoylmethane in methanol solution by employing density functional theory-based molecular dynamics and quantum chemical calculations. Our computations showed that homogeneous electric fields with intensities corresponding to those found in polar solvents are able to considerably alter the proton transfer barrier height in the gas phase. In methanol solution, the proton position is correlated with the difference in electrostatic potentials on the oxygen atoms of dibenzoylmethane even when dibenzoylmethane-methanol hydrogen bonding is lacking. On a timescale of our simulation, the hydrogen bonding and solvent electrostatics tend to localize the proton on different oxygen atoms. These findings provide an insight into the importance of the solvent electric field on the structure of a strong intramolecular hydrogen bond.

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

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

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

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