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.     

Source of Cooperativity in Halogen‐Bonded Haloamine Tetramers

Inspired by the isostructural motif in α‐bromoacetophenone oxime crystals, we investigated halogen–halogen bonding in haloamine quartets. Our Kohn–Sham molecular orbital and energy decomposition analysis reveal a synergy that can be traced to a charge‐transfer interaction in the halogen‐bonded tetramers. The halogen lone‐pair orbital on one monomer donates electrons into the unoccupied σ*_N?X orbital on the perpendicular N?X bond of the neighboring monomer. This interaction has local σ symmetry. Interestingly, we discovered a second, somewhat weaker donor–acceptor interaction of local π symmetry, which partially counteracts the aforementioned regular σ‐symmetric halogen‐bonding orbital interaction. The halogen–halogen interaction in haloamines is the first known example of a halogen bond in which back donation takes place. We also find that this cooperativity in halogen bonds results from the reduction of the donor–acceptor orbital‐energy gap that occurs every time a monomer is added to the aggregate.     

The many flavours of halogen bonds – message from experimental electron density and Raman spectroscopy

Experimental electron‐density studies based on high‐resolution diffraction experiments allow halogen bonds between heavy halogens to be classified. The topological properties of the electron density in Cl…Cl contacts vary smoothly as a function of the interaction distance. The situation is less straightforward for halogen bonds between iodine and small electronegative nucleophiles, such as nitrogen or oxygen, where the electron density in the bond critical point does not simply increase for shorter distances. The number of successful charge–density studies involving iodine is small, but at least individual examples for three cases have been observed. (a) Very short halogen bonds between electron‐rich nucleophiles and heavy halogen atoms resemble three‐centre–four‐electron bonds, with a rather symmetric heavy halogen and without an appreciable σ hole. (b) For a narrow intermediate range of halogen bonds, the asymmetric electronic situation for the heavy halogen with a pronounced σ hole leads to rather low electron density in the (3,?1) critical point of the halogen bond; the properties of this bond critical point cannot fully describe the nature of the associated interaction. (c) For longer and presumably weaker contacts, the electron density in the halogen bond critical point is only to a minor extent reduced by the presence of the σ hole and hence may be higher than in the aforementioned case. In addition to the electron density and its derived properties, the halogen–carbon bond distance opposite to the σ hole and the Raman frequency for the associated vibration emerge as alternative criteria to gauge the halogen‐bond strength. We find exceptionally long C—I distances for tetrafluorodiiodobenzene molecules in cocrystals with short halogen bonds and a significant red shift for their Raman vibrations.     

The ability of Ex2Box4+ to interact with guests containing π‐electron‐rich and π‐electron‐poor moieties

The ability of Ex 2 Box⁴⁺ as a host, able to trap guests containing both π‐electron rich (polycyclic aromatic hydrocarbons‐PAHs) and π‐electron poor (quinoid‐ and nitro‐PAHs) moieties was investigated to shed light on the main factors that control the host–guest (HG) interaction. The nature of the HG interactions was elucidated by energy decomposition (EDA‐NOCV), noncovalent interaction (NCI), and magnetic response analyses. EDA‐NOCV reveals that dispersion contributions are the most significant to sustain the HG interaction, while electrostatic and orbital contributions are very tiny. In fact, no significant covalent character in the HG interactions was observed. The obtained results point strictly to NCIs, modulated by dispersion contributions. Regardless of whether the guests contain π‐electron‐rich or π‐electron‐poor moieties, and no significant charge‐transfer was observed. All in all, HG interactions between guests 3‐14 and host 2 are predominantly modulated by π‐π stacking.     

Theoretical study on the interactions of halogen‐bonds and pnicogen‐bonds in phosphine derivatives with Br2, BrCl,and BrF

The MP2 ab initio quantum chemistry methods were utilized to study the halogen‐bond and pnicogen‐bond system formed between PH₂X (X = Br, CH₃, OH, CN, NO₂, CF₃) and BrY (Y = Br, Cl, F). Calculated results show that all substituent can form halogen‐bond complexes while part substituent can form pnicogen‐bond complexes. Traditional, chlorine‐shared and ion‐pair halogen‐bonds complexes have been found with the different substituent X and Y. The halogen‐bonds are stronger than the related pnicogen‐bonds. For halogen‐bonds, strongly electronegative substituents which are connected to the Lewis acid can strengthen the bonds and significantly influenced the structures and properties of the compounds. In contrast, the substituents which connected to the Lewis bases can produce opposite effects. The interaction energies of halogen‐bonds are 2.56 to 32.06 kcal·mol^?1; The strongest halogen‐bond was found in the complex of PH₂OH???BrF. The interaction energies of pnicogen‐bonds are in the range 1.20 to 2.28 kcal·mol^?1; the strongest pnicogen‐bond was found in PH₂Br???Br₂ complex. The charge transfer of lp(P) ? σ*(Br? Y), lp(F) ? σ*(Br? P), and lp(Br) ? σ*(X? P) play important roles in the formation of the halogen‐bonds and pnicogen‐bonds, which lead to polarization of the monomers. The polarization caused by the halogen‐bond is more obvious than that by the pnicogen‐bond, resulting in that some halogen‐bonds having little covalent character. The symmetry adapted perturbation theory (SAPT) energy decomposition analysis showes that the halogen‐bond and pnicogen‐bond interactions are predominantly electrostatic and dispersion, respectively.     

Halogen Bonding: An Interim Discussion

Halogen bonding is a noncovalent interaction that is receiving rapidly increasing attention because of its significance in biological systems and its importance in the design of new materials in a variety of areas, for example, electronics, nonlinear optical activity, and pharmaceuticals. The interactions can be understood in terms of electrostatics/polarization and dispersion; they involve a region of positive electrostatic potential on a covalently bonded halogen and a negative site, such as the lone pair of a Lewis base. The positive potential, labeled a σ hole, is on the extension of the covalent bond to the halogen, which accounts for the characteristic near‐linearity of halogen bonding. In many instances, the lateral sides of the halogen have negative electrostatic potentials, allowing it to also interact favorably with positive sites. In this discussion, after looking at some of the experimental observations of halogen bonding, we address the origins of σ holes, the factors that govern the magnitudes of their electrostatic potentials, and the properties of the resulting complexes with negative sites. The relationship of halogen and hydrogen bonding is examined. We also point out that σ‐hole interactions are not limited to halogens, but can also involve covalently bonded atoms of Groups IV–VI. Examples of applications in biological/medicinal chemistry and in crystal engineering are mentioned, taking note that halogen bonding can be “tuned” to fit various requirements, that is, strength of interaction, steric factors, and so forth.     

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.     

Nature of MoH···I bonds in Cp2Mo(L)H···I‐C≡C‐R Complexes (L=H,CN, PPh2, C(CH3)3; R=NO2, Cl,Br, H,OH, CH3, NH2)

The nature of the MoH···I bond in Cp₂Mo(L)H···I‐C≡C‐R (L= H, CN, PPh₂, C(CH₃)₃; R=NO₂, Cl, Br, H, OH, CH₃, NH₂) was investigated using electrostatic potential analysis, topological analysis of the electron density, energy decomposition analysis and natural bond orbital analysis. The calculated results show that MoH···I interactions in the title complexes belong to halogen‐hydride bond, which is similar to halogen bonds, not hydrogen bonds. Different to the classical halogen bonds, the directionality of MoH···I bond is low; Although electrostatic interaction is dorminant, the orbital interactions also play important roles in this kind of halogen bond, and steric interactions are weak; the strength of H···I bond can tuned by the most positive electrostatic potential of the I atom. As the electron‐withdrawing ability of the R substituent in the alkyne increases, the electrostatic potential maximum of the I atom increases, which enhances the strength of the H···I halogen bond, as well as the electron transfer.     

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.     

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.     

Is a MH (M = Be and Mg) radical a better electron donor in halogen‐hydride interaction?: A theoretical comparison with HMH

Quantum chemical calculations have been performed to investigate the complexes of HMH? XCCH, HMH? XCF₃, MH? XCCH, and MH? XCF₃ (M = Be and Mg; X = Cl, Br, and I) at the MP2/aug‐cc‐pVTZ level. The geometrical, energetic, and spectroscopic parameters were analyzed for these complexes. The results show that the MH is a better electron donor in the halogen‐hydride interaction than HMH. The enhancement of halogen‐hydride interaction increases in the order of Cl < Br < I, Be < Mg, and XCCH < XCF₃. The halogen‐hydride interaction was understood with natural bond orbital, atoms in molecules, and electrostatic potentials. © 2012 Wiley Periodicals, Inc.     

A Crystallographic Charge Density Study of the Partial Covalent Nature of Strong N⋅⋅⋅Br Halogen Bonds

The covalent nature of strong N?Br???N halogen bonds in a cocrystal ( 2 ) of N‐bromosuccinimide ( NBS ) with 3,5‐dimethylpyridine ( lut ) was determined from X‐ray charge density studies and compared to a weak N?Br???O halogen bond in pure crystalline NBS ( 1 ) and a covalent bond in bis(3‐methylpyridine)bromonium cation (in its perchlorate salt ( 3 ). In 2 , the donor N?Br bond is elongated by 0.0954 Å, while the Br???acceptor distance of 2.3194(4) is 1.08 Å shorter than the sum of the van der Waals radii. A maximum electron density of 0.38 e Å^?3 along the Br???N halogen bond indicates a considerable covalent contribution to the total interaction. This value is intermediate to 0.067 e Å^?3 for the Br???O contact in 1 , and approximately 0.7 e Å^?3 in both N?Br bonds of the bromonium cation in 3 . A calculation of the natural bond order charges of the contact atoms, and the σ*(N1?Br) population of NBS as a function of distance between NBS and lut , have shown that charge transfer becomes significant at a Br???N distance below about 3 Å.     

A theoretical approach to substituent effects. Structural consequences of electrostatic and orbital interactions in model mono- and disubstituted methanes

Ab initio molecular orbital theory is used to examine the effect of substituents on bond lengths in mono- and disubstituted methanes. The relative importance of electrostatic and orbital interaction terms are assessed. The results suggest that for substituents (X) which show powerful σ effects and weak π interactions (e.g., F), the changes in bond length are due primarily to the electrostatic component except in some disubstituted methanes in which case the change in the hyperconjugative ability of the C—X bond is also important. On the other hand, substituents X which show weak σ effects but powerful π interactions (e.g., NH₂) affect bond lengths primarily through hyperconjugative interaction of a filled or vacant π-type orbital on X with the adjacent bonds.     

Why are the Homoleptic Diyl Complexes M(InR)4 with M = Ni and Pt Stable Compounds while only Ni(CO)4 but not Pt(CO)4 can become Isolated? A Theoretical Study of M(EMe)44 and M(CO)4 (M = Ni,Pd, Pt; E = B,Al, Ga,In, Tl)

The equilibrium geometries and first bond dissociation energies of the homoleptic complexes M(EMe)₄ and M(CO)₄ with M = Ni, Pd, Pt and E = B, Al, Ga, In, Tl have been calculated at the gradient corrected DFT level using the BP86 functionals. The electronic structure of the metal‐ligand bonds has been examined with the topologial analysis of the electron density distribution. The nature of the bonding is revealed by partitioning the metal‐ligand interaction energies into contributions by electrostatic attraction, covalent bonding and Pauli repulsion. The calculated data show that the M‐CO and M‐EMe bonding is very similar. However, the M‐EMe bonds of the lighter elements E are much stronger than the M‐CO bonds. The bond energies of the latter are as low or even lower than the M‐TlMe bonds. The main reason why Pd(CO)₄ and Pt(CO)₄ are unstable at room temperature in a condensed phase can be traced back to the already rather weak bond energy of the Ni‐CO bond. The Pd‐L bond energies of the complexes with L = CO and L = EMe are always 10 — 20 kcal/mol lower than the Ni‐L bond energies. The calculated bond energy of Ni(CO)₄ is only D_o = 27 kcal/mol. Thus, the bond energy of Pd(CO)₄ is only D_o = 12 kcal/mol. The first bond dissociation energy of Pt(CO)₄ is low because the relaxation energy of the Pt(CO)₃ fragment is rather high. The low bond energies of the M‐CO bonds are mainly caused by the relatively weak electrostatic attraction and by the comparatively large Pauli repulsion. The σ and π contributions to the covalent M‐CO interactions have about the same strength. The π bonding in the M‐EMe bonds is less than in the M‐CO bonds but it remains an important part of the bond energy. The trends of the electrostatic and covalent contributions to the bond energies and the σ and π bonding in the metal‐ligand bonds are discussed.     

Substituent effects in halogen bonding complexes between aromatic donors and acceptors: a comprehensive ab initio study

Substituent effects in halogen bonding complexes involving aromatic rings are investigated. We have analyzed how the interaction energy (the RI-MP2/aug-cc-pVDZ level of theory) is affected by the substitution in both halogen bond donor and acceptor aromatic moieties. In addition, we have used two different aromatic electron donor molecules pyridine and cyanobenzene, which allow us to study the effect of having the electron donor nitrogen atom forming part of the ring or outside the ring (-CN). Interestingly, the effect of the substituents on the interaction energies is similar in both cases. We have obtained the Hammett's plots for four combinations of aromatic donors and acceptors and in all cases we have obtained good regression plots (interaction energies vs. Hammett's σ parameter). We have also studied and compared bifurcated halogen bonds using both possible combinations, that is two donors and one acceptor and vice versa. In addition, we have analyzed the effect of the solvent on the interaction energies using COSMO. Finally, we have used Bader's theory of "atoms-in-molecules" to demonstrate that the electron density computed at the bond critical point that emerges upon complexation can be used as a measure of bond order in this noncovalent interaction.     

Nature of the Hydrogen Bond Enhanced Halogen Bond

The factors responsible for the enhancement of the halogen bond by an adjacent hydrogen bond have been quantitatively explored by means of state-of-the-art computational methods. It is found that the strength of a halogen bond is enhanced by ca. 3 kcal/mol when the halogen donor simultaneously operates as a halogen bond donor and a hydrogen bond acceptor. This enhancement is the result of both stronger electrostatic and orbital interactions between the XB donor and the XB acceptor, which indicates a significant degree of covalency in these halogen bonds. In addition, the halogen bond strength can be easily tuned by modifying the electron density of the aryl group of the XB donor as well as the acidity of the hydrogen atoms responsible for the hydrogen bond.     

Theoretical ab Initio Study on Cooperativity Effects between Nitro π-hole and Halogen Bonding Interactions

This article analyzes the interplay between nitro's π-hole and halogen–bonding (XB) interactions in nitroarenes. Remarkable cooperativity effects are observed when π–hole and XB interactions coexist in the same complex. The nitroarene presents two π-holes, one approximately over the N atom of the nitro group and the other over the aromatic ring, being the former more positive. The interplay between both interactions has been analyzed in terms of energetic and geometric features of the complexes, which are computed at the RI-MP2/def2-TZVPD level of theory. Molecular electrostatic potential (MEP) surface calculations have been used to explore the variation of the MEP values at the π-hole upon the formation of halogen bonding interactions between the nitroarene and CF₃X (X=Cl, Br and I) molecules. In addition, the Bader's theory of atoms in molecules” (AIM) is used to characterize the interactions by means of the distribution of bond critical points and bond paths and to analyze their strengthening or weakening depending upon the variation of charge density at critical points. The aforementioned computational methods are adequate to examine how these interactions mutually influence each other. Natural bond orbital (NBO) and noncovalent interaction plot (NCIPlot) computational tools have been also used in some representative complexes to further analyze cooperativity effects. Finally, the Cambridge Structural Database (CSD) is used to provide some experimental evidence.     

Competition between chalcogen bond and halogen bond interactions in YOX<Subscript>4</Subscript>:NH<Subscript>3</Subscript> (Y = S,Se; X = F,Cl, Br) complexes: An ab initio investigation

Using ab initio calculations, the geometries, interaction energies and bonding properties of chalcogen bond and halogen bond interactions between YOX₄ (Y = S, Se; X = F, Cl, Br) and NH₃ molecules are studied. These binary complexes are formed through the interaction of a positive electrostatic potential region (σ-hole) on the YOX₄ with the negative region in the NH₃. The ab initio calculations are carried out at the MP2/aug-cc-pVTZ level, through analysis of molecular electrostatic potentials, quantum theory of atoms in molecules and natural bond orbital methods. Our results indicate that even though the chalcogen and halogen bonds are mainly dominated by electrostatic effects, but the polarization and dispersion effects also make important contributions to the total interaction energy of these complexes. The examination of interaction energies suggests that the chalcogen bond is always favored over the halogen bond for all of the binary YOX₄:NH₃ complexes.     

A new unconventional halogen bond C?X···H?M between HCCX (X = Cl and Br) and HMH (M = Be and Mg): An ab initio study

In this article, a new type of halogen‐bonded complex YCCX···HMY (X = Cl, Br; M = Be, Mg; Y = H, F, CH₃) has been predicted and characterized at the MP2/aug‐cc‐pVTZ level. We named it as halogen‐hydride halogen bonding. In each YCCX···HMY complex, a halogen bond is formed between the positively charged X atom and the negatively charged H atom. This new kind of halogen bond has similar characteristics to the conventional halogen bond, such as the elongation of the C? X bond and the red shift of the C? X stretch frequency upon complexation. The interaction strength of this type of halogen bond is in a range of 3.34–10.52 kJ/mol, which is smaller than that of dihydrogen bond and conventional halogen bond. The nature of the electrostatic interaction in this type of halogen bond has also been unveiled by means of the natural bond orbital, atoms in molecules, and energy decomposition analyses. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Halogen bond versus hydrogen bond: The many‐body interactions approach

The analysis of interrelation between halogen bond and hydrogen bond in the (RX)(HNC)(HCN) complexes (R = CH₃, CF₃ and X = Cl, Br, I) was performed on the basis of DFT calculations. Both two‐body additive contributions and three‐body nonadditive contributions to the total interaction energy were discussed. QTAIM was used for topological analysis of electron density. Additionally, QTAIM analysis of electron density was performed for both two‐ and three‐body complexes. The electron charge transfer in trimers showed the dual character of the fragment with halogen atom involved into the investigated interactions—it acts as Lewis acid and Lewis base, depending on the type of interaction considered. The effect of cooperativity of X‐ and H‐bonding was assessed on the basis of many‐body interaction energy and electron density analysis. Additionally, an alternative two‐body model with the same situation (in the context of intermolecular interactions) is investigated. The anti‐cooperative effect was found also for this model.