Divalent carbon atom as the proton acceptor in hydrogen bonding

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

Halogen Bond to Experimentally Significant N-Heterocyclic Carbenes (I, IMe2, IiPr2, ItBu2, IPh2, IMes2, IDipp2, IAd2; I = Imidazol-2-ylidene)

The subjects of the article are halogen bonds between either XCN or XCCH (X = Cl, Br, I) and the carbene carbon atom in imidazol-2-ylidene (I) or its derivatives (IR2) with experimentally significant and systematically increased R substituents at both nitrogen atoms: methyl = Me, iso-propyl = iPr, tert-butyl = tBu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad. It is shown that the halogen bond strength increases in the order Cl < Br < I and the XCN molecule forms stronger complexes than XCCH. Of all the carbenes considered, IMes2 forms the strongest and also the shortest halogen bonds with an apogee for complex IMes2⋯ICN for which D0 = 18.71 kcal/mol and dC⋯I = 2.541 Å. In many cases, IDipp2 forms as strong halogen bonds as IMes2. Quite the opposite, although characterized by the greatest nucleophilicity, ItBu2 forms the weakest complexes (and the longest halogen bonds) if X ≠ Cl. While this finding can easily be attributed to the steric hindrance exerted by the highly branched tert-butyl groups, it appears that the presence of the four C-H⋯X hydrogen bonds may also be of importance here. Similar situation occurs in the case of complexes with IAd2.

Halogen bonds with carbenes acting as Lewis base units: complexes of imidazol-2-ylidene: theoretical analysis and experimental evidence

ωB97XD/aug-cc-pVDZ and ωB97XD/aug-cc-pVTZ calculations were performed on complexes of imidazol-2-ylidene that are linked by halogen bonds. This singlet carbene acts as the Lewis base through a lone electron pair located at the carbon centre. The XCCH, XCN and X₂ units were chosen here as those that interact through the X Lewis acid halogen centre (X = Cl, Br and I); if X = F the complexes are linked by interactions which are not classified as halogen bonds. The properties of interactions that occur in complexes are analyzed using the results of DFT calculations which are supported by parameters derived from the Quantum Theory of Atoms in Molecules, QTAIM, and the Natural Bond Orbital, NBO, approaches. The energy decomposition analysis, EDA, applied here provided additional characteristics of interactions linking complexes analyzed. The majority of complexes are linked by the medium in strength and strong halogen bonds which often possess characteristics typical for covalent bonds. Searches through the Cambridge Structural Database were also performed and structures analogues to complexes analyzed theoretically were found, and these structures are also discussed in this study.

Nature of Beryllium, Magnesium, and Zinc Bonds in Carbene⋯MX2 (M = Be, Mg, Zn; X = H, Br) Dimers Revealed by the IQA, ETS-NOCV and LED Methods

The nature of beryllium−, magnesium− and zinc−carbene bonds in the cyclopropenylidene⋯MX2 (M = Be, Mg, Zn; X = H, Br) and imidazol-2-ylidene⋯MBr2 dimers is investigated by the joint use of the topological QTAIM-based IQA decomposition scheme, the molecular orbital-based ETS-NOCV charge and energy decomposition method, and the LED energy decomposition approach based on the state-of-the-art DLPNO-CCSD(T) method. All these methods show that the C⋯M bond strengthens according to the following order: Zn < Mg << Be. Electrostatics is proved to be the dominant bond component, whereas the orbital component is far less important. It is shown that QTAIM/IQA underestimates electrostatic contribution for zinc bonds with respect to both ETS-NOCV and LED schemes. The σ carbene→MX2 donation appears to be much more important than the MX2→ carbene back-donation of π symmetry. The substitution of hydrogen atoms by bromine (X in MX2) strengthens the metal−carbene bond in all cases. The physical origin of rotational barriers has been unveiled by the ETS-NOCV approach.

On the Coexistence of the Carbene⋯H-D Hydrogen Bond and Other Accompanying Interactions in Forty Dimers of N-Heterocyclic-Carbenes (I, IMe2, IiPr2, ItBu2, IMes2, IDipp2, IAd2; I = imidazol-2-ylidene) and Some Fundamental Proton Donors (HF, HCN, H2O, MeOH, NH3)

The subject of research is forty dimers formed by imidazol-2-ylidene (I) or its derivative (IR2) obtained by replacing the hydrogen atoms in both N-H bonds with larger important and popular substituents of increasing complexity (methyl = Me, iso-propyl = iPr, tert-butyl = tBu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad) and fundamental proton donor (HD) molecules (HF, HCN, H2O, MeOH, NH3). While the main goal is to characterize the generally dominant C⋯H-D hydrogen bond engaging a carbene carbon atom, an equally important issue is the often omitted analysis of the role of accompanying secondary interactions. Despite the often completely different binding possibilities of the considered carbenes, and especially HD molecules, several general trends are found. Namely, for a given carbene, the dissociation energy values of the IR2⋯HD dimers increase in the following order: NH3< H2O < HCN ≤ MeOH ≪ HF. Importantly, it is found that, for a given HD molecule, IDipp2 forms the strongest dimers. This is attributed to the multiplicity of various interactions accompanying the dominant C⋯H-D hydrogen bond. It is shown that substitution of hydrogen atoms in both N-H bonds of the imidazol-2-ylidene molecule by the investigated groups leads to stronger dimers with HF, HCN, H2O or MeOH. The presented results should contribute to increasing the knowledge about the carbene chemistry and the role of intermolecular interactions, including secondary ones.

Carbon-Centered Hydrogen Bonds in Proteins

Hydrogen bonding (H-bonding) without lone pair(s) of electrons and π-electrons is a concept developed 2-3 years ago. H-bonds involving less electronegative tetrahedral carbon are beyond the classical concept of H-bonds. Herein, we present the first report on H-bonds with tetravalent carbons in proteins. A special bonding arrangement is needed to increase the negative charge density around the sp³-hybridized carbon atom. Therefore, less electronegative elements such as As and Mg, when bonded to sp³-C, enable the C-atoms as H-bond acceptors. Careful protein structure analysis aided by several quantum chemical calculations suggests that these H-bonds are weak to moderate in strength. We developed an empirical equation to estimate the C-H···C H-bond energy in proteins from the distances between the C- and H-atoms. In proteins, the binding energies range from -5.4 to -14.0 kJ/mol. The C-H···C H-bonds assist the substrate binding in proteins. We also explored the potential role of these carbon-centered H-bonds in C-H bond activation through σ-bond metathesis. To our surprise, contribution from these H-bonds is almost of similar magnitude as that from C-H···π H-bonds for C-H bond activation.

π-Hole Tetrel Bonds—Lewis Acid Properties of Metallylenes

The MP2/aug-cc-pVTZ calculations were performed on the dihalometallylenes to indicate their Lewis acid and Lewis base sites. The results of the Cambridge Structural Database search show corresponding and related crystal structures where the tetrel center often possesses the configuration of a trigonal bipyramid or octahedron. The calculations were also carried out on dimers of dichlorogermylene and dibromogermylene and on complexes of these germylenes with one and two 1,4-dioxide molecules. The Ge⋯Cl, Ge⋯Br, and Ge⋯O interactions are analyzed. The Ge⋯O interactions in the above mentioned germylene complexes may be classified as the π-hole tetrel bonds. The MP2 calculations are supported by the results of the Quantum Theory of Atoms in Molecules (QTAIM) and the Natural Bond Orbital (NBO) approaches.

Theoretical Study of N-Heterocyclic-Carbene-ZnX2 (X = H, Me, Et) Complexes

This article discusses the properties of as many as 30 carbene–ZnX2 (X = H, Me, Et) complexes featuring a zinc bond C⋯Zn. The group of carbenes is represented by imidazol-2-ylidene and its nine derivatives (labeled as IR), in which both hydrogen atoms of N-H bonds have been substituted by R groups with various spatial hindrances, from the smallest Me, iPr, tBu through Ph, Tol, and Xyl to the bulkiest Mes, Dipp, and Ad. The main goal is to study the relationship between type and size of R and X and both the strength of C⋯Zn and the torsional angle of the ZnX2 plane with respect to the plane of the imidazol-2-ylidene ring. Despite the considerable diversity of R and X, the range of dC⋯Zn is quite narrow: 2.12–2.20 Å. On the contrary, D0 is characterized by a fairly wide range of 18.5–27.4 kcal/mol. For the smallest carbenes, the ZnX2 molecule is either in the plane of the carbene or is only slightly twisted with respect to it. The twist angle becomes larger and more varied with the bulkier R. However, the value of this angle is not easy to predict because it results not only from the presence of steric effects but also from the possible presence of various interatomic interactions, such as dihydrogen bonds, tetrel bonds, agostic bonds, and hydrogen bonds. It has been shown that at least some of these interactions may have a non-negligible influence on the structure of the IR–ZnX2 complex. This fact should be taken into account in addition to the commonly discussed R⋯X steric repulsion.

Methyl groups as widespread Lewis bases in noncovalent interactions

It is well known that, under certain conditions, C(sp³) atoms behave, via their σ-hole, as Lewis acids in tetrel bonding. Here, we show that methyl groups, when bound to atoms less electronegative than carbon, can counterintuitively participate in noncovalent interactions as electron density donors. Thousands of experimental structures are found in which methyl groups behave as Lewis bases to establish alkaline, alkaline earth, triel, tetrel, pnictogen, chalcogen and halogen bonds. Theoretical calculations confirm the high directionality and significant strength of the interactions that arise from a common pattern based on the electron density holes model. Moreover, despite the absence of lone pairs, methyl groups are able to transfer charge from σ bonding orbitals into empty orbitals of the electrophile to reinforce the attractive interaction.

Study of Beryllium, Magnesium, and Spodium Bonds to Carbenes and Carbodiphosphoranes

The aim of this article is to present results of theoretical study on the properties of C⋯M bonds, where C is either a carbene or carbodiphosphorane carbon atom and M is an acidic center of MX2 (M = Be, Mg, Zn). Due to the rarity of theoretical data regarding the C⋯Zn bond (i.e., the zinc bond), the main focus is placed on comparing the characteristics of this interaction with C⋯Be (beryllium bond) and C⋯Mg (magnesium bond). For this purpose, theoretical studies (ωB97X-D/6-311++G(2df,2p)) have been performed for a large group of dimers formed by MX2 (X = H, F, Cl, Br, Me) and either a carbene ((NH2)2C, imidazol-2-ylidene, imidazolidin-2-ylidene, tetrahydropyrymid-2-ylidene, cyclopropenylidene) or carbodiphosphorane ((PH3)2C, (NH3)2C) molecule. The investigated dimers are characterized by a very strong charge transfer effect from either the carbene or carbodiphosphorane molecule to the MX2 one. This may even be over six times as strong as in the water dimer. According to the QTAIM and NCI method, the zinc bond is not very different than the beryllium bond, with both featuring a significant covalent contribution. However, the zinc bond should be definitely stronger if delocalization index is considered.

Theoretical study of hydrogen bonding interactions in substituted nitroxide radicals

Interaction energy (E_int) of hydrogen bonded complexes of nitroxide radicals can be assessed in terms of the deepest minimum of molecular electrostatic potential (V_min).

Toxic target of trans-crotonaldehyde in mitochondria altered by diallyl disulfides for anti-myocardial ischemia

The purposes of this study were to probe spectral behaviors of the toxic targets in trans-crotonaldehyde (TCA) in mitochondria altered by diallyl disulfides (DADS) derived from garlic. Ultraviolet absorption spectra showed that when ethanol as a solvent, the DADS blue shifted the peak of TCA from 318 nm to 312 nm. In mitochondria, the DADS further blue displaced the peak of TCA from 312 nm to 308 nm. Raman spectra displayed that the SS of DADS directly interacted with the CC toxic target of TCA, then the CCC of DADS interacted with the CHO toxic target of TCA. When DADS to TCA was 1:2, the DADS was the most powerful for the removal of the CC and CHO toxic targets of TCA. Study suggested that the SS of DADS altered the CC toxic target of TCA, while the CCC of DADS eliminated the CHO toxic target of TCA via local electron delocalization. The above two together clearly depicted the spectral behaviors of the toxic targets of TCA in mitochondria altered by DADS. These results are of great significance and value to elucidate the effects of garlic organic polysulfide on myocardial ischemia for the extensive development and use of garlic extracts.

The electronic donation and frequency shifts on the YCCH⋯BH₄⁻ boron-bonded complexes (Y=H, CH₃, CF₃ and CCl₃)

In this theoretical work, the tetrahydroborate ion (BH4(-)) was used as proton acceptor in the formation of the YCC-H⋯BH4(-) complexes (Y=H, CH3, CCl3 and CF3). Using B3LYP/6-311++G(d,p) level of theory, the results of structure corroborate with the analyses of infrared spectra showing that the changes in the bond lengths are in good agreement with the frequency shifts of the HCC-H, H3CCC-H, Cl3CCC-H and F3CCC-H proton donors. Based on the calculations carried out by the Quantum Theory of Atoms in Molecules (QTAIM), the reductions of electronic density corroborate with the red shifts in the frequencies of the C-H bonds. In addition to that, the C-H bonds are polarized because the contributions of s orbital diminish whereas of p increase. In line with this, the variations on the atomic radii computed via QTAIM calculations show that carbon outweigh hydrogen as follows (ΔrC>ΔrH). This scenario is indirectly supported by the Bent's rule of the chemical bonding. Although the interaction energies (corrected with BSSE and ZPE) vary between -19 and -67 kJ mol(-1), these complexes interact without covalent character.

The nature of NO-bonding in N-oxide group

The nature of the NO-bond in the N-oxide group was investigated by means of combined theoretical calculations (including QTAIM and NBO approaches) and statistical analyses of the contents of crystal structure databases. The N-O bond in the N-oxide group should be classified as the NO donating bond with an important contribution of ON back-donation (of the π-electron type, when available). The visualization of the Laplacian of electron density in the region of an oxygen valence sphere suggests the presence of two lone pairs for the imine-N-oxide group (characterized by effective ON back-donation). A detailed bonding analysis performed by means of natural resonance theory indicates that the N→O bond is of an order of magnitude clearly greater than 1. In addition, the stability of the N→O bond in various N-oxides was estimated. The analyses of the hydrogen- and halogen-bonded complexes of the N-oxides reveal strong Lewis basicity of the N-oxide group. The formation of H- and X-bonding leads to N→O bond elongation due to its structural, topological and spectroscopic characteristics. Moreover, in pyridine-N-oxide, the electron-withdrawing -NO2 group additionally stabilizes the N→O bond, whereas the opposite effect can be observed for the electron-donating-NH2 substituent. This is due to a substituent effect on the π-type ON back-donation. As a result, the oxygen atom in pyridine-N-oxide may change its availability during intermolecular interaction formation, as revealed in the interaction energy, which changes by about half of the estimated total interaction energy.

Nature of a hydride-halogen bond. A SAPT-, QTAIM-, and NBO-based study

The nature of a hydride-halogen bond is investigated by means of high-level quantum mechanical calculations expended with symmetry-adapted perturbation theory (SAPT), quantum theory of atoms in molecules (QTAIM), and natural bond orbital (NBO) methods. As model hydride-halogen bonded systems complexes between either LiH or HBeH and either XCF(3) or XCCH (X = F, Cl, Br, I) are used. It is shown that the formation of a hydride-halogen bond leads to the elongation of the R(δ+)-H(δ-) hydride bond, which is accompanied by the blue shift of the ν(R-H) stretching vibration frequency and the increase of the IR intensity of this mode. All these effects, although untypical in the case of, e.g., hydrogen bonds, can be considered as rather typical for hydride-halogen bonded systems. The decomposition of the interaction energy based on the SAPT method clearly indicates the dominant role of the induction term, thus the inductive nature of a hydride-halogen bond in opposition to previous findings. NBO-based analysis indicates the charge transfer from the hydride molecule to the more remote parts of the halogen donor and that the elongation of the R-H bond is caused by the charge outflow from the σ(RH) bonding orbital.

Divalent N(I) compounds with two lone pairs on nitrogen

Carbon with the C(0) state has been reported recently, examples of which were known for the past decades. Silicon in the Si(0) state and phosphorus in the P(I) state are also known experimentally. This prompted us to search for divalent N(I) compounds, which resulted in the identification of ::N(←L)(2)(⊕) systems with bicoordinated nitrogen in the N(I) formal oxidation state. It was found that several biguanide derivatives (especially in their protonated state) belong to this class. Quantum chemical analysis provided the structural details, molecular orbitals, charge localization (vs delocalization) trends, etc. This class of compounds has been found to be characterized by two lone pairs on the central nitrogen, very similar to the central carbon in divalent C(0) compounds (::C(←L)(2)). The new bonding environment for nitrogen reported in this article, divalent nitrogen N(I), is clearly different from the nitrenium ions NR(2)(⊕). The electronic structure and reactivity of representative examples of this novel class of divalent nitrogen N(I) systems (::N(←L)(2)(⊕)) have been analyzed in detail, in terms of molecular orbitals, atomic charges, protonation energies, complexation energies with Lewis acids like BH(3), AlCl(3), and AuCl and compared with those of divalent C(0) systems.

Anhydrous TEMPO-H: reactions of a good hydrogen atom donor with low-valent carbon centres

In this paper, we report a novel synthesis of anhydrous 1-hydroxy-2,2,6,6-tetramethyl-piperidine (TEMPO-H). An X-ray crystal structure and full characterization of the compound are included. Compared to hydrated TEMPO-H, its anhydrous form exhibits improved stability and a differing chemical reactivity. The reactions of anhydrous TEMPO-H with a variety of low-valent carbon centres are described. For example, anhydrous TEMPO-H was reacted with 1,3-bis(2,4,6-trimethylphenyl)imidazol-2-ylidene (IMes), an unsaturated NHC. Crystals of [CHNC(6)H(2)(CH(3))(3)](2)C···HO(NC(5)H(6)(CH(3))(4)), IMes···TEMPO-H, were isolated and a crystal structure determined. The experimental structure is compared to the results of theoretical calculations on the hydrogen-bonded dimer. Anhydrous TEMPO-H was also reacted with the saturated NHC, 1,3-bis(2,6-diisopropylphenyl)imidazolidin-2-ylidene (SIPr), giving the product [CH(2)Ni-Pr(2)C(6)H(3)](2)CH···O(NC(5)H(6)(CH(3))(4)). In contrast, the reaction of hydrated TEMPO-H with 1,3-bis(2,6-diisopropylphenyl)imidazolidin-2-ylidene gave small amounts of the hydrolysis product, N-(2,6-diisopropylphenyl)-N-[2-(2,6-diisopropylphenylamino)ethyl]formamide. Finally, anhydrous TEMPO-H was reacted with (triphenylphosphoranylidene)ketene to generate Ph(3)PC(H)C(=O)O(NC(5)H(6)(CH(3))(4)). A full characterization of the product, including an X-ray crystal structure, is described.

When does carbonylation of carbenes yield ketenes? A theoretical study with implications for synthesis

Quantum-chemical calculations using DFT and ab initio methods have been carried out for 32 carbenes RR'C which comprise different classes of compounds and the associated ketenes RR'C═C═O. The calculated singlet-triplet gaps ΔE(S-T) of the carbenes exhibit a very high correlation with the bond dissociation energies (BDEs) of the ketenes. An energy decomposition analysis of the RR'C-CO bond using the triplet states of the carbene and CO as interacting fragments supports the assignment of ΔE(S-T) as the dominant factor for the BDE but also shows that the specific interactions of the carbene may sometimes compensate for the S/T gap. The trend of the interaction energy ΔE(int) values is mainly determined by the Pauli repulsion between the carbene and CO. The stability of amino-substituted ketenes strongly depends on the destabilizing conjugation between the nitrogen lone-pair orbital and the ketene double bonds. There is a ketene structure of the unsaturated N-heterocyclic carbene parent compound NHC1 with CO as a local energy minimum on the potential-energy surface. However, the compound NHC1-CO is thermodynamically unstable toward dissociation. The saturated homologue NHC2-CO has only a very small bond dissociation energy of D(e) = 3.2 kcal/mol. The [3]ferrocenophane-type compound FeNHC-CO has a BDE of D(e) = 16.0 kcal/mol.