Clusters of Ammonium Cation—Hydrogen Bond versus σ‐Hole Bond

MP2/6‐311++G(d,p) calculations were performed on the NH₄⁺ ??? (HCN)_n and NH₄⁺ ??? (N₂)_n clusters (n=1–8), and interactions within them were analyzed. It was found that for molecules of N₂ and HCN, the N centers play the role of the Lewis bases, whereas the ammonium cation acts as the Lewis acid, as it is characterized by sites of positive electrostatic potential, that is, H atoms and the sites located at the N atom in the extension of the H?N bonds. Hence, the coordination number for the ammonium cation is eight, and two types of interactions of this cation with the Lewis base centers are possible: N?H ??? N hydrogen bonds and H?N ??? N interactions that are classified as σ‐hole bonds. Redistribution of the electronic charge resulting from complexation of the ammonium cation was analyzed. On the one hand, the interactions are similar, as they lead to electronic charge transfer from the Lewis base (HCN or N₂ in this study) to NH₄⁺. On the other hand, the hydrogen bond results in the accumulation of electronic charge on the N atom of the NH₄⁺ ion, whereas the σ‐hole bond results in the depletion of the electronic charge on this atom. Quantum theory of “atoms in molecules” and the natural bond orbital method were applied to deepen the understanding of the nature of the interactions analyzed. Density functional theory/natural energy decomposition analysis was used to analyze the interactions of the ammonium ion with various types of Lewis bases. Different correlations between the geometrical, energetic, and topological parameters were found and discussed.     

π‐Hole Bonds: Boron and Aluminum Lewis Acid Centers

MP2/aug‐cc‐pVTZ calculations were performed on complexes of boron and aluminum trihydrides and trihalides with hydrogen cyanide (ZH₃‐NCH and ZX₃‐NCH; Z=B, Al; X=F, Cl). The complexes are linked through the B???N and Al???N interactions, which are named as triel bonds and which are classified as π‐hole bonds. It was found that they possess numerous characteristics of typical covalent bonds, since they are ruled mainly by processes of the electron charge shift from the Lewis base to the Lewis acid unit. Other configurations of the ZH₃‐NCH and ZX₃‐NCH complexes linked by the dihydrogen, hydrogen, and halogen bonds were found. However, these interactions are much weaker than the corresponding π‐hole bonds. The quantum theory of atoms in molecules and the natural bond orbital approaches were applied to characterize the complexes and interactions analyzed. The crystal structures of triel trihydrides and triel trihalides were also analyzed for comparison with the results of calculations.     

Creating σ‐Holes through the Formation of Beryllium Bonds

Through the use of ab initio theoretical models based on MP2/aug‐cc‐pVDZ‐optimized geometries and CCSD(T)/aug‐cc‐pVTZ and CCSD(T)/aug‐c‐pVDZ total energies, it has been shown that the significant electron density rearrangements that follow the formation of a beryllium bond may lead to the appearance of a σ‐hole in systems that previously do not exhibit this feature, such as CH₃OF, NO₂F, NO₃F, and other fluorine‐containing systems. The creation of the σ‐hole is another manifestation of the bond activation–reinforcement (BAR) rule. The appearance of a σ‐hole on the F atoms of CH₃OF is due to the enhancement of the electronegativity of the O atom that participates in the beryllium bond. This atom recovers part of the charge transferred to Be by polarizing the valence density of the F into the bonding region. An analysis of the electron density shows that indeed this bond becomes reinforced, but the F atom becomes more electron deficient with the appearance of the σ‐hole. Importantly, similar effects are also observed even when the atom participating in the beryllium bond is not directly attached to the F atom, as in NO₂F, NO₃F, or NCF. Hence, whereas the isolated CH₃OF, NO₂F, and NO₃F are unable to yield F ??? Base halogen bonds, their complexes with BeX₂ derivatives are able to yield such bonds. Significant cooperative effects between the new halogen bond and the beryllium bond reinforce the strength of both noncovalent interactions.     

Homonuclear chalcogen–chalcogen bond interactions in complexes pairing YO3 and YHX molecules (Y = S,Se; X = H,Cl, Br,CCH, NC,OH, OCH3): Influence of substitution and cooperativity

MP2/aug‐cc‐pVTZ calculations are performed on complexes of YO₃ (Y = S, Se) with a series of electron‐donating chalcogen bases YHX (X = H, Cl, Br, CCH, NC, OH, OCH₃). These complexes are formed through the interaction of a positive electrostatic potential region (π‐hole) on the YO₃ molecule with the negative region in YHX. Interaction energies of the binary O₃Y???YHX complexes are in the range of ?4.37 to ?12.09 kcal/mol. The quantum theory of atoms in molecules and the natural bond orbital analysis were applied to characterize the nature of interactions. It was found that the formation and stability of these binary complexes are ruled mainly by electrostatic effects, although the electron charge transfer from YHX to YO₃ unit also seems to play an important role. In addition, mutual influence between the Y???N and Y???Y interactions is studied in the ternary HCN???O₃Y???YHX complexes. The results indicate that the formation of a Y???N interaction tends to weaken Y???Y bond in the ternary systems. Although the Y???Y interaction is weaker than the Y???N one, however, both types of interactions seem to compete with each other in the HCN???O₃Y???YHX complexes. © 2016 Wiley Periodicals, Inc.     

Comparative studies on group III σ‐hole and π‐hole interactions

The σ‐hole of M₂H₆ (M = Al, Ga, In) and π‐hole of MH₃ (M = Al, Ga, In) were discovered and analyzed, the bimolecular complexes M₂H₆···NH₃ and MH₃···N₂P₂F₄ (M = Al, Ga, In) were constructed to carry out comparative studies on the group III σ‐hole interactions and π‐hole interactions. The two types of interactions are all partial‐covalent interactions; the π‐hole interactions are stronger than σ‐hole interactions. The electrostatic energy is the largest contribution for forming the σ‐hole and π‐hole interaction, the polarization energy is also an important factor to form the M···N interaction. The electrostatic energy contributions to the interaction energy of the σ‐hole interactions are somewhat greater than those of the π‐hole interactions. However, the polarization contributions for the π‐hole interactions are somewhat greater than those for the σ‐hole interactions. © 2016 Wiley Periodicals, Inc.     

The competition of Y⋯o and X⋯n halogen bonds to enhance the group V σ‐hole interaction in the NCY⋯oPH3⋯NCX and OPH3⋯NCX⋯NCY (X,YF,Cl, and Br) complexes

The positive electrostatic potentials (ESP) outside the σ‐hole along the extension of O? P bond in O?PH₃ and the negative ESP outside the nitrogen atom along the extension of the C? N bond in NCX could form the Group V σ‐hole interaction O?PH₃?NCX. In this work, the complexes NCY?O?PH₃?NCX and O?PH₃?NCX?NCY (X, Y?F, Cl, Br) were designed to investigate the enhancing effects of Y?O and X?N halogen bonds on the P?N Group V σ‐hole interaction. With the addition of Y?O halogen bond, the V _{S, max} values outside the σ‐hole region of O?PH₃ becomes increasingly positive resulting in a stronger and more polarizable P?N interaction. With the addition of X?N halogen bond, the V _{S, min} values outside the nitrogen atom of NCX becomes increasingly negative, also resulting in a stronger and more polarizable P?N interaction. The Y?O halogen bonds affect the σ‐hole region (decreased density region) outside the phosphorus atom more than the P?N internuclear region (increased density region outside the nitrogen atom), while it is contrary for the X?N halogen bonds. © 2015 Wiley Periodicals, Inc.     

Cooperativity between the Halogen Bond and the Hydrogen Bond in H3N⋅⋅⋅XY⋅⋅⋅HF Complexes (X,Y=F,Cl, Br)

Ab initio calculations are used to provide information on H₃N???XY???HF triads (X, Y=F, Cl, Br) each having a halogen bond and a hydrogen bond. The investigated triads include H₃N???Br₂‐HF, H₃N???Cl₂???HF, H₃N???BrCI???HF, H₃N???BrF???HF, and H₃N???ClF???HF. To understand the properties of the systems better, the corresponding dyads are also investigated. Molecular geometries, binding energies, and infrared spectra of monomers, dyads, and triads are studied at the MP2 level of theory with the 6‐311++G(d,p) basis set. Because the primary aim of this study is to examine cooperative effects, particular attention is given to parameters such as cooperative energies, many‐body interaction energies, and cooperativity factors. The cooperative energy ranges from ?1.45 to ?4.64 kcal mol^?1, the three‐body interaction energy from ?2.17 to ?6.71 kcal mol^?1, and the cooperativity factor from 1.27 to 4.35. These results indicate significant cooperativity between the halogen and hydrogen bonds in these complexes. This cooperativity is much greater than that between hydrogen bonds. The effect of a halogen bond on a hydrogen bond is more pronounced than that of a hydrogen bond on a halogen bond.     

Strong N−X⋅⋅⋅O−N Halogen Bonds: A Comprehensive Study on N‐Halosaccharin Pyridine N‐Oxide Complexes

A study of the strong N?X???^?O?N⁺ (X=I, Br) halogen bonding interactions reports 2×27 donor×acceptor complexes of N‐halosaccharins and pyridine N‐oxides (PyNO). DFT calculations were used to investigate the X???O halogen bond (XB) interaction energies in 54 complexes. A simplified computationally fast electrostatic model was developed for predicting the X???O XBs. The XB interaction energies vary from ?47.5 to ?120.3 kJ mol^?1; the strongest N?I???^?O?N⁺ XBs approaching those of 3‐center‐4‐electron [N?I?N]⁺ halogen‐bonded systems (ca. 160 kJ mol^?1). ¹H NMR association constants (K_XB) determined in CDCl₃ and [D₆]acetone vary from 2.0×10⁰ to >10⁸ m ^?1 and correlate well with the calculated donor×acceptor complexation enthalpies found between ?38.4 and ?77.5 kJ mol^?1. In X‐ray crystal structures, the N‐iodosaccharin‐PyNO complexes manifest short interaction ratios (R_XB) between 0.65–0.67 for the N?I???^?O?N⁺ halogen bond.     

Ab Initio and DFT Studies on CO2 Interacting with Znq+–Imidazole (q=0, 1, 2) Complexes: Prediction of Charge Transfer through σ‐ or π‐Type Models

Using first‐principles methodologies, the equilibrium structures and the relative stability of CO₂@[Zn^q+Im] (where q=0, 1, 2; Im=imidazole) complexes are studied to understand the nature of the interactions between the CO₂ and Zn^q+–imidazole entities. These complexes are considered as prototype models mimicking the interactions of CO₂ with these subunits of zeolitic imidazolate frameworks or Zn enzymes. These computations are performed using both ab initio calculations and density functional theory. Dispersion effects accounting for long‐range interactions are considered. Solvent (water) effects were also considered using a polarizable continuum model approach. Natural bond orbital, charge, frontier orbital and vibrational analyses clearly reveal the occurrence of charge transfer through covalent and noncovalent interactions. Moreover, it is found that CO₂ can adsorb through more favorable π‐type stacking as well as σ‐type hydrogen‐bonding interactions. The inter‐monomer interaction potentials show a significant anisotropy that might induce CO₂ orientation and site‐selectivity effects in porous materials and in active sites of Zn enzymes. Hence, this study provides valuable information about how CO₂ adsorption takes place at the microscopic level within zeolitic imidazolate frameworks and biomolecules. These findings might help in understanding the role of such complexes in chemistry, biology and material science for further development of new materials and industrial applications.     

Dimers of N‐Heterocyclic Carbene Copper,Silver, and Gold Halides: Probing Metallophilic Interactions through Electron Density Based Concepts

Homobimetallic metallophilic interactions between copper, silver, and gold‐based [(NHC)MX]‐type complexes (NHC=N‐heterocyclic carbene, i.e, 1,3,4‐trimethyl‐4,5‐dihydro‐1H‐1,2,4‐triazol‐5‐ylidene; X=F, Cl, Br, I) were investigated by means of ab initio interaction energies, Ziegler–Rauk‐type energy‐decomposition analysis, the natural orbital for chemical valence (NOCV) framework, and the noncovalent interaction (NCI) index. It was found that the dimers of these complexes predominantly adopt a head‐to‐tail arrangement with typical M ??? M distance of 3.04–3.64 Å, in good agreement with the experimental X‐ray structure determined for [{(NHC)AuCl}₂], which has an Au ??? Au distance of 3.33 Å. The interaction energies between silver‐ and gold‐based monomers are calculated to be about ?25 kcal mol^?1, whereas that for the Cu congener is significantly lower (?19.7 kcal mol^?1). With the inclusion of thermal and solvent contributions, both of which are destabilizing, by about 15 and 8 kcal mol^?1, respectively, an equilibrium process is predicted for the formation of dimer complexes. Energy‐decomposition analysis revealed a dominant electrostatic contribution to the interaction energy, besides significantly stabilizing dispersion and orbital interactions. This electrostatic contribution is rationalized by NHC(δ⁺) ??? halogen(δ^?) interactions between monomers, as demonstrated by electrostatic potentials and derived charges. The dominant NOCV orbital indicates weakening of the π backdonation in the monomers on dimer formation, whereas the second most dominant NOCV represents an electron‐density deformation according to the formation of a very weak M ??? M bond. One of the characteristic signals found in the reduced density gradient versus electron density diagram corresponds to the noncovalent interactions between the metal centers of the monomers in the NCI plots, which is the manifestation of metallophilic interaction.     

Hydrogen bond and σ‐hole interaction in M2C=S···HCN (M = H,F, Cl,Br, HO,H3C,H2N) complex: Dual roles of C=S group and substitution effect

An ab initio computational study of the dual functions of C?S group in the M₂C?S ··· HCN (M = H, F, Cl, Br, HO, H₃C, H₂N) complex has been performed at the MP2(Full)/aug‐cc‐pVTZ level. The C?S group can act as both the electron donor and acceptor, thus two minima complexes were found for each molecular pairs. The interaction energy of hydrogen bond in the F, Cl, or Br substituted complexes is less negative than that in the corresponding H₂CS one, while the interaction energy of the σ‐hole interaction is more negative. The OH substitution weakens the hydrogen bond, whereas the H₃C and H₂N substitution strengthens it. The σ‐hole interaction in the HO, H₃C, and H₂N complexes is very weak. The substitution effect has been understood with electrostatic induction and conjugation effects. The energy decomposition analysis has been performed for the halogen‐substituted complexes. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012.     

CH⋅⋅⋅N Hydrogen‐Bonding Interaction in 7‐Azaindole:CHX3 (X=F,Cl) Complexes

The C?H???Y (Y=hydrogen‐bond acceptor) interactions are somewhat unconventional in the context of hydrogen‐bonding interactions. Typical C?H stretching frequency shifts in the hydrogen‐bond donor C?H group are not only small, that is, of the order of a few tens of cm^?1, but also bidirectional, that is, they can be red or blue shifted depending on the hydrogen‐bond acceptor. In this work we examine the C?H???N interaction in complexes of 7‐azaindole with CHCl₃ and CHF₃ that are prepared in the gas phase through supersonic jet expansion using the fluorescence depletion by infra‐red (FDIR) method. Although the hydrogen‐bond acceptor, 7‐azaindole, has multiple sites of interaction, it is found that the C?H???N hydrogen‐bonding interaction prevails over the others. The electronic excitation spectra suggest that both complexes are more stabilized in the S₁ state than in the S₀ state. The C?H stretching frequency is found to be red shifted by 82 cm^?1 in the CHCl₃ complex, which is the largest redshift reported so far in gas‐phase investigations of 1:1 haloform complexes with various substrates. In the CHF₃ complex the observed C?H frequency is blue shifted by 4 cm^?1. This is at variance with the frequency shifts that are predicted using several computational methods; these predict at best a redshift of 8.5 cm^?1. This discrepancy is analogous to that reported for the pyridine‐CHF₃ complex [W. A. Herrebout, S. M. Melikova, S. N. Delanoye, K. S. Rutkowski, D. N. Shchepkin, B. J. van der Veken, J. Phys. Chem. A­ 2005 , 109, 3038], in which the blueshift is termed a pseudo blueshift and is shown to be due to the shifting of levels caused by Fermi resonance between the overtones of the C?H bending and stretching modes. The dissociation energies, (D₀), of the CHCl₃ and CHF₃ complexes are computed (MP2/aug‐cc‐pVDZ level) as 6.46 and 5.06 kcal mol^?1, respectively.     

Structures of Trichloromethyl Thiocyanate,CCl3SCN,in Gaseous and Crystalline State

Trichloromethyl thiocyanate, CCl₃SCN, was structurally studied in both the gas and crystal phases by means of gas electron diffraction (GED) and single‐crystal X‐ray diffraction (XRD), respectively. Both experimental studies and quantum chemical calculations indicate a staggered orientation of the CCl₃ group relative to the SCN group. This conclusion is supported by the similarity of the C?SCN bond length to that of the anti‐structure of CH₂ClSCN (Berrueta Martínez et al. Phys. Chem. Chem. Phys. 2015, 17, 15805–15812). 1 Bond lengths and angles are similar for gas and crystal CCl₃SCN structures; however, the crystal structure presents different intermolecular interactions. These include halogen and chalcogen type interactions, the geometry of which was studied. Characteristic C‐Y???N angles (Y=Cl or S) close to 180° provide evidence for typical σ‐hole interactions along the halogen/chalcogen?carbon bond in N???Cl and N???S, intermolecular units.     

The Prominent Enhancing Effect of the Cation–π Interaction on the Halogen–Hydride Halogen Bond in M1⋅⋅⋅C6H5X⋅⋅⋅HM2

We designed M¹???C₆H₅X???HM² (M¹=Li⁺, Na⁺; X=Cl, Br; M²=Li, Na, BeH, MgH) complexes to enhance halogen–hydride halogen bonding with a cation–π interaction. The interaction strength has been estimated mainly in terms of the binding distance and the interaction energy. The results show that halogen–hydride halogen bonding is strengthened greatly by a cation–π interaction. The interaction energy in the triads is two to six times as much as that in the dyads. The largest interaction energy is ?8.31 kcal mol^?1 for the halogen bond in the Li⁺???C₆H₅Br???HNa complex. The nature of the cation, the halogen donor, and the metal hydride influence the nature of the halogen bond. The enhancement effect of Li⁺ on the halogen bond is larger than that of Na⁺. The halogen bond in the Cl donor has a greater enhancement than that in the Br one. The metal hydride imposes its effect in the order HBeH<HMgH<HNa<HLi for the Cl complex and HBeH<HMgH<HLi<HNa for the Br complex. The large cooperative energy indicates that there is a strong interplay between the halogen–hydride halogen bonding and the cation–π interaction. Natural bond orbital and energy decomposition analyses indicate that the electrostatic interaction plays a dominate role in enhancing halogen bonding by a cation–π interaction.     

Influence of the Li⋅⋅⋅π Interaction on the H/X⋅⋅⋅π Interactions in HOLi⋅⋅⋅C6H6⋅⋅⋅HOX/XOH (X=F,Cl, Br,I) Complexes

The influences of the Li???π interaction of C₆H₆???LiOH on the H???π interaction of C₆H₆???HOX (X=F, Cl, Br, I) and the X???π interaction of C₆H₆???XOH (X=Cl, Br, I) are investigated by means of full electronic second‐order Møller–Plesset perturbation theory calculations and “quantum theory of atoms in molecules” (QTAIM) studies. The binding energies, binding distances, infrared vibrational frequencies, and electron densities at the bond critical points (BCPs) of the hydrogen bonds and halogen bonds prove that the addition of the Li???π interaction to benzene weakens the H???π and X???π interactions. The influences of the Li???π interaction on H???π interactions are greater than those on X???π interactions; the influences of the H???π interactions on the Li???π interaction are greater than X???π interactions on Li???π interaction. The greater the influence of Li???π interaction on H/X???π interactions, the greater the influences of H/X???π interactions on Li???π interaction. QTAIM studies show that the intermolecular interactions of C₆H₆???HOX and C₆H₆???XOH are mainly of the π type. The electron densities at the BCPs of hydrogen bonds and halogen bonds decrease on going from bimolecular complexes to termolecular complexes, and the π‐electron densities at the BCPs show the same pattern. Natural bond orbital analyses show that the Li???π interaction reduces electron transfer from C₆H₆ to HOX and XOH.     

The Role of Molecular Electrostatic Potentials in the Formation of a Halogen Bond in Furan⋅⋅⋅XY and Thiophene⋅⋅⋅XY Complexes

The halogen bonding of furan???XY and thiophene???XY (X=Cl, Br; Y=F, Cl, Br), involving σ‐ and π‐type interactions, was studied by using MP2 calculations and quantum theory of “atoms in molecules” (QTAIM) studies. The negative electrostatic potentials of furan and thiophene, as well as the most positive electrostatic potential (V_S,max) on the surface of the interacting X atom determined the geometries of the complexes. Linear relationships were found between interaction energy and V_S,max of the X atom, indicating that electrostatic interactions play an important role in these halogen‐bonding interactions. The halogen‐bonding interactions in furan???XY and thiophene???XY are weak, “closed‐shell” noncovalent interactions. The linear relationship of topological properties, energy properties, and the integration of interatomic surfaces versus V_S,max of atom X demonstrate the importance of the positive σ hole, as reflected by the computed V_S,max of atom X, in determining the topological properties of the halogen bonds.     

Heterolytic activation of C?H bond in methane with (HNCHCHNH)M(CH3) (M = Pd+, Pt+, Rh+, Ir+, Rh,Ir): Comparative density functional study of activation mechanisms

Activation of methane by oxidative addition and σ‐bond metathesis has been investigated for (N‐N)M(CH₃) (M = Pd⁺, Pt⁺, Rh⁺, Ir⁺, Rh, Ir; N‐N = (HN?CH? CH?NH) using different density functional approaches. The pathway of oxidative addition is in general favored, the exceptions being Pd⁺ and Rh⁺. Oxidative addition is clearly more favorable for the third‐row metal complexes than those of the second row. The third‐row metal complexes also tend to have a lower activation barrier for σ‐bond metathesis than those of the second row. In each case, the oxidative addition is preceded by formation of a sigma complex. The bonding energies of these complexes are significantly stronger for the cationic systems. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Does Fluorine Participate in Halogen Bonding?

When R is sufficiently electron withdrawing, the fluorine in the R?F molecules could interact with electron donors (e.g., ammonia) and form a noncovalent bond (F ??? N). Although these interactions are usually categorized as halogen bonding, our studies show that there are fundamental differences between these interactions and halogen bonds. Although the anisotropic distribution of electronic charge around a halogen is responsible for halogen bond formations, the electronic charge around the fluorine in these molecules is spherical. According to source function analysis, F is the sink of electron density at the F ??? N BCP, whereas other halogens are the source. In contrast to halogen bonds, the F ??? N interactions cannot be regarded as lump–hole interactions; there is no hole in the valence shell charge concentration (VSCC) of fluorine. Although the quadruple moment of Cl and Br is mainly responsible for the existence of σ‐holes, it is negligibly small in the fluorine. Here, the atomic dipole moment of F plays a stabilizing role in the formation of F ??? N bonds. Interacting quantum atoms (IQA) analysis indicates that the interaction between halogen and nitrogen in the halogen bonds is attractive, whereas it is repulsive in the F ??? N interactions. Virial‐based atomic energies show that the fluorine, in contrast to Cl and Br, stabilize upon complex formation. According to these differences, it seems that the F ??? N interactions should be referred to as “fluorine bond” instead of halogen bond.     

A Computational Study on 3‐Azonia‐, 3‐Phosphonia‐, and 3‐Arsoniaspiro[2.2]pentanes and Related Three‐Membered Heterocycles

A theoretical study at the ab initio MP2/6‐311++G(d,p) level of theory is carried out to characterize several heterocyclic spiro[2.2]pentane cations with N, P, and As as spiro atoms. The strain and relative stability of the spiropentanes are obtained through isodesmic reactions. Nucleus‐independent chemical shifts (NICS) and 3D NICS isosurfaces show σ‐aromatic characteristics, similar to those found in cyclopropane. The interaction with the Cl^? anion, which results in four different stationary structures, is studied and characterized by means of the atoms in molecules methodology, and Cl ??? pnicogen, Cl ??? H, and Cl ??? C interactions are found. The most stable structure in all cases corresponds to opening of one of the three‐membered rings, due to the attack of the Cl atom, and C?Cl bond formation. Furthermore, the reaction with the 3‐boranuidaspiro[2.2]pentane anion results in the formation of a new compound through cleavage of one ring of both reactants.