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.     

Enhancement of Iodine–Hydride Interaction by Substitution and Cooperative Effects in NCX–NCI–HMY Complexes

The NCX‐NCI‐HMY (X=H, Cl, Br, I, Li; M=Be, Mg; Y=H, Li, Na) trimers are investigated to find ways to enhance the iodine–hydride interaction. The interaction energy in the NCI–HMH dimer is ?2.87 and ?5.87 kcal mol^?1 for M=Be and Mg, respectively. When the free H atom in the NCI–HMH dimer is replaced with an alkali atom, the interaction energy is enhanced greatly. When NCX is added into this dimer, the interaction energy of the iodine–hydride interaction is increased by 9–45 % and its increased percentage follows the order X=Cl<Br<H<I<Li and M=Be<Mg. The combination of the alkali substitution and the cooperativity results in a more prominent enhancing effect. The largest interaction energy is found for the NCLi–NCI–HMgLi trimer (?7.03 kcal mol^?1). The influence of the I???H interaction on the X???N interaction is also studied in the trimers. Both types of interactions are analyzed with NBO, AIM, and MEP. The interaction energy in the trimer is also unveiled by a many‐body analysis.     

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.     

Density functional theory study of calix[4]arene‐N‐azacrown‐5, calix[4]arene‐N‐phenyl‐azacrown‐5, and their complexes with alkali‐metal cations: Na+, K+, and Rb+

Theoretical studies of 1,3‐alternate‐25,27‐bis(1‐methoxyethyl)calix[4]arene‐azacrown‐5 ( L₁ ), 1,3‐alternate‐25,27‐bis(1‐methoxyethyl)calix[4]arene‐N‐phenyl‐azacrown‐5 ( L₂ ), and the corresponding complexes M⁺/ L of L₁ and L₂ with the alkali‐metal cations: Na⁺, K⁺, and Rb⁺ have been performed using density functional theory (DFT) at B3LYP/6‐31G* level. The optimized geometric structures obtained from DFT calculations are used to perform natural bond orbital (NBO) analysis. The two main types of driving force metal–ligand and cation–π interactions are investigated. The results indicate that intermolecular electrostatic interactions are dominant and the electron‐donating oxygen offer lone pair electrons to the contacting RY* (1‐center Rydberg) or LP* (1‐center valence antibond lone pair) orbitals of M⁺ (Na⁺, K⁺, and Rb⁺). What's more, the cation–π interactions between the metal ion and π‐orbitals of the two rotated benzene rings play a minor role. For all the structures, the most pronounced changes in geometric parameters upon interaction are observed in the calix[4]arene molecule. In addition, an extra pendant phenyl group attached to nitrogen can promote metal complexation by 3D encapsulation greatly. In addition, the enthalpies of complexation reaction and hydrated cation exchange reaction had been studied by the calculated thermodynamic data. The calculated results of hydrated cation exchange reaction are in a good agreement with the experimental data for the complexes. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

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.     

Synthesis,X‐ray characterization and density functional theory studies of N6‐benzyl‐N6‐methyladenine–M(II) complexes (M = Zn,Cd): The prominent role of π–π, C–H···π and anion–π interactions

We report the synthesis and X‐ray characterization of the N⁶‐benzyl‐N⁶‐methyladenine ligand (L) and three metal complexes, namely [Zn(HL)Cl₃]·H₂O ( 1 ), [Cd(HL)₂Cl₄] ( 2 ) and [H₂L]₂[Cd₃(μ‐L)₂(μ‐Cl)₄Cl₆]·3H₂O ( 3 ). Complex 1 consists of the 7H‐adenine tautomer protonated at N3 and coordinated to a tetrahedral Zn(II) metal centre through N9. The octahedral Cd(II) in complex 2 is N⁹‐coordinated to two N⁶‐benzyl‐N⁶‐methyladeninium ligands (7H‐tautomer protonated at N3) that occupy apical positions and four chlorido ligands form the basal plane. Compound 3 corresponds to a trinuclear Cd(II) complex, where the central Cd atom is six‐coordinated to two bridging μ‐L and four bridging μ‐Cl ligands. The other two Cd atoms are six‐coordinated to three terminal chlorido ligands, to two bridging μ‐Cl ligands and to the bridging μ‐L through N3. Essentially, the coordination patterns, degree of protonation and tautomeric forms of the nucleobase dominate the solid‐state architectures of 1 – 3 . Additionally, the hydrogen‐bonding interactions produced by the endocyclic N atoms and NH groups stabilize high‐dimensional‐order supramolecular assemblies. Moreover, energetically strong anion–π and lone pair (lp)–π interactions are important in constructing the final solid‐state architectures in 1 – 3 . We have studied the non‐covalent interactions energetically using density functional theory calculations and rationalized the interactions using molecular electrostatic potential surfaces and Bader's theory of atoms in molecules. We have particularly analysed cooperative lp–π and anion–π interactions in 1 and π⁺–π⁺ interactions in 3 .     

Hydrogen‐Bond Cooperativity in Formamide2–Water: A Model for Water‐Mediated Interactions

The rotational spectrum of formamide₂–H₂O formed in a supersonic jet has been characterized by Fourier‐transform microwave spectroscopy. This adduct provides a simple model of water‐mediated interaction involving the amide linkages, as occur in protein folding or amide‐association processes, showing the interplay between self‐association and solvation. Mono‐substituted ¹³C, ¹⁵N, ¹⁸O, and ²H isotopologues have been observed and their data used to investigate the structure. The adduct forms an almost planar three‐body sequential cycle. The two formamide molecules link on one side through an N?H???O hydrogen bond and on the other side through a water‐mediated interaction with the formation of C=O???H?O and O???H?N hydrogen bonds. The analysis of the quadrupole coupling effects of two ¹⁴N‐nuclei reveals the subtle inductive forces associated to cooperative hydrogen bonding. These forces are involved in the changes in the C=O and C?N bond lengths with respect to pure formamide.     

Energy decomposition analysis of cation–π, metal ion–lone pair,hydrogen bonded,charge‐assisted hydrogen bonded,and π–π interactions

This study probes the nature of noncovalent interactions, such as cation–π, metal ion–lone pair (M–LP), hydrogen bonding (HB), charge‐assisted hydrogen bonding (CAHB), and π–π interactions, using energy decomposition schemes—density functional theory (DFT)–symmetry‐adapted perturbation theory and reduced variational space. Among cation–π complexes, the polarization and electrostatic components are the major contributors to the interaction energy (IE) for metal ion–π complexes, while for onium ion–π complexes ( , , , and ) the dispersion component is prominent. For M–LP complexes, the electrostatic component contributes more to the IE except the dicationic metal ion complexes with H₂S and PH₃ where the polarization component dominates. Although electrostatic component dominates for the HB and CAHB complexes, dispersion is predominant in π–π complexes.Copyright © 2015 Wiley Periodicals, Inc.     

Intermolecular Interactions in Li+‐glyme and Li+‐glyme–TFSA− Complexes: Relationship with Physicochemical Properties of [Li(glyme)][TFSA] Ionic Liquids

The stabilization energies (ΔE_form) calculated for the formation of the Li⁺ complexes with mono‐, di‐ tri‐ and tetra‐glyme (G1, G2, G3 and G4) at the MP2/6‐311G** level were ?61.0, ?79.5, ?95.6 and ?107.7 kcal mol^?1, respectively. The electrostatic and induction interactions are the major sources of the attraction in the complexes. Although the ΔE_form increases by the increase of the number of the O???Li contact, the ΔE_form per oxygen atom decreases. The negative charge on the oxygen atom that has contact with the Li⁺ weakens the attractive electrostatic and induction interactions of other oxygen atoms with the Li⁺. The binding energies calculated for the [Li(glyme)]⁺ complexes with TFSA^? anion (glyme=G1, G2, G3, and G4) were ?106.5, ?93.7, ?82.8, and ?70.0 kcal mol^?1, respectively. The binding energies for the complexes are significantly smaller than that for the Li⁺ with the TFSA^? anion. The binding energy decreases by the increase of the glyme chain length. The weak attraction between the [Li(glyme)]⁺ complex (glyme=G3 and G4) and TFSA^? anion is one of the causes of the fast diffusion of the [Li(glyme)]⁺ complex in the mixture of the glyme and the Li salt in spite of the large size of the [Li(glyme)]⁺ complex. The HOMO energy level of glyme in the [Li(glyme)]⁺ complex is significantly lower than that of isolated glyme, which shows that the interaction of the Li⁺ with the oxygen atoms of glyme increases the oxidative stability of the glyme.     

Enhanced Aerogen–π Interaction by a Cation–π Force

The interaction between a noble gas atom and an aromatic π‐electron system, which mainly originates from the London dispersion force, is very weak and has not attracted enough attention yet. Herein, we reported a type of notably enhanced aerogen–π interaction between cation–π systems and noble gas atoms. The binding strength of a divalent cation–π system with a xenon atom is comparable to a moderate hydrogen bond (up to ca. 7 kcal mol^?1), whereas krypton and argon atoms produce slightly weaker interactions. Energy‐decomposition analysis reveals that the induction interaction is responsible for the stabilization of divalent cation–π?Xe species besides the dispersion interaction. Our results might be helpful to increase the understanding of some unsolved mysteries of aerogens.     

A Carbohydrate–Anion Recognition System in Aprotic Solvents

A carbohydrate–anion recognition system in nonpolar solvents is reported, in which complexes form at the B‐faces of β‐D ‐pyranosides with H1‐, H3‐, and H5‐cis patterns similar to carbohydrate–π interactions. The complexation effect was evaluated for a range of carbohydrate structures; it resulted in either 1:1 carbohydrate–anion complexes, or 1:2 complex formation depending on the protection pattern of the carbohydrate. The interaction was also evaluated with different anions and solvents. In both cases it resulted in significant binding differences. The results indicate that complexation originates from van der Waals interactions or weak CH ??? A^? hydrogen bonds between the binding partners and is related to electron‐withdrawing groups of the carbohydrates as well as increased hydrogen‐bond‐accepting capability of the anions.     

Hydrogen‐Bonding Interactions and Properties of Energetic Nitroamino[1,3,5]triazine‐Based Guanidinium Salts: DFT‐D and QTAIM Studies

The intramolecular hydrogen‐bonding interactions and properties of a series of nitroamino[1,3,5]triazine‐based guanidinium salts were studied by using the dispersion‐corrected density functional theory method (DFT‐D). Results show that there are evident LP(N or O; LP=lone pair)→σ*(N? H) orbital interactions related to O???H? N or N???H? N hydrogen bonds. Quantum theory of atoms in molecules (QTAIM) was applied to characterize the intramolecular hydrogen bonds. For the guanidinium salts studied, the intramolecular hydrogen bonds are associated with a seven‐ or eight‐membered pseudo‐ring. The guanylurea cation is more helpful for improving the thermal stabilities of the ionic salts than other guanidinium cations. The contributions of different substituents on the triazine ring to the thermal stability increase in the order of ? NO₂23 (? ONO₂)2. Energy decomposition analysis shows that the salts are stable owing to electrostatic and orbital interactions between the ions, whereas the dispersion energy has very small contributions. Moreover, the salts exhibit relatively high densities in the range of 1.62–1.89 g cm^?3. The detonation velocities and pressures lie in the range of 6.49–8.85 km s^?1 and 17.79–35.59 GPa, respectively, which makes most of them promising explosives.     

PI by NMR: Probing CH–π Interactions in Protein–Ligand Complexes by NMR Spectroscopy

While CH–π interactions with target proteins are crucial determinants for the affinity of arguably every drug molecule, no method exists to directly measure the strength of individual CH–π interactions in drug–protein complexes. Herein, we present a fast and reliable methodology called PI (π interactions) by NMR, which can differentiate the strength of protein–ligand CH–π interactions in solution. By combining selective amino‐acid side‐chain labeling with ¹H‐¹³C NMR, we are able to identify specific protein protons of side‐chains engaged in CH–π interactions with aromatic ring systems of a ligand, based solely on ¹H chemical‐shift values of the interacting protein aromatic ring protons. The information encoded in the chemical shifts induced by such interactions serves as a proxy for the strength of each individual CH–π interaction. PI by NMR changes the paradigm by which chemists can optimize the potency of drug candidates: direct determination of individual π interactions rather than averaged measures of all interactions.     

σ‐Hole Bond Versus Hydrogen Bond: From Tetravalent to Pentavalent N,P, and As Atoms

Ab initio calculations were performed on complexes of ZH₄⁺ (Z=N, P, As) and their fluoro derivatives, ZFH₃⁺ and ZF₄⁺, with a HCN (or LiCN) molecule acting as the Lewis base through the nitrogen electronegative center. It was found that the complexes are linked by the Z? H???N hydrogen bond or another type of noncovalent interaction in which the tetravalent heavy atom of the cation acts as the Lewis acid center, that is, when the Z???N link exists, which may be classified as the σ‐hole bond. The formation of the latter interaction is usually preferable to the formation of the corresponding hydrogen bond. The Z???N interaction may be also considered as the preliminary stage of the S_N2 reaction. This is supported by the observation that for a short Z???N contact, the corresponding complex geometry coincides with the trigonal‐bipyramidal geometry typical for the transition state of the S_N2 reaction. The Z???N interaction for some of complexes analyzed here possesses characteristics typical for covalent bonds. Numerous interrelations between geometrical, topological and energetic parameters are discussed. The natural bond orbital method as well as the Quantum Theory of “Atoms in Molecules” is applied to characterize interactions in the analyzed complexes. The experimental evidences of the existence of these interactions, based on the Cambridge Structure Database search, are also presented. In addition, it is justified that mechanisms of the formation of the Z???N interactions are similar to the processes occurring for the other noncovalent links. The formation of Z???N interaction as well as of other interactions may be explained with the use of the σ‐hole concept.     

Rational Design of a Nonbasic Molecular Receptor for Selective NH4+/K+ Complexation in the Gas Phase

A new molecular receptor ( 1 ) for ammonium recognition has been designed and constructed by using only carbon atoms. This molecular receptor can co‐exist in two different isoenergetic conformations but, upon complexation, the conformers are no longer isoenergetic, and a basket‐shaped conformation becomes clearly more stable. The pre‐organised tetrahedral structure of this basket‐shaped molecule favours the complexation of ammonium ions by N? H???π interactions with the four phenyl groups of the host. A similar behaviour is not observed in a similar, but less pre‐organised, reference molecule. ESI‐MS competition experiments show that 1 is able to bind NH₄⁺ over K⁺ selectively. This is the first example of a neutral molecular receptor that shows a remarkable NH₄⁺/K⁺ selectivity. DFT‐calculations provide insight into the nature of host–guest interactions of both 1? NH₄⁺ and 1? K⁺ complexes as well as in the mechanism involved in multiple cation–π interactions and the influence of these interactions on the conformational stability and the selective binding of the host.     

Influence of metal ion on chelate–aryl stacking interactions

CCSD(T)/CBS and DFT methods are employed to study the stacking interactions of acetylacetonate‐type (acac‐type) chelates of nickel, palladium, and platinum with benzene. The strongest chelate–aryl stacking interactions are formed by nickel and palladium chelate, with interaction energies of −5.75 kcal mol⁻¹ and −5.73 kcal mol⁻¹, while the interaction of platinum chelate is weaker, with interaction energy of −5.36 kcal mol⁻¹. These interaction energies are significantly stronger than stacking of two benzenes, −2.73 kcal mol⁻¹. The strongest nickel and palladium chelate–aryl interactions are with benzene center above the metal area, while the strongest platinum chelate–aryl interaction is with the benzene center above the C2 atom of the acac‐type chelate ring. These preferences arise from very different electrostatic potentials above the metal ions, ranging from very positive above nickel to slightly negative above platinum. While the differences in electrostatic potentials above metal atoms cause different geometries with the most stable interaction among the three metals, the dispersion (correlation energy) component is the largest contribution to the total interaction energy for all three metals.     

Exclusive π Encapsulation of Light Alkali Metal Cations by a Neutral Molecule

Cation–π interactions are one of the most important classes of noncovalent bonding, and are seen throughout biology, chemistry, and materials science. However, in almost every documented case, these interactions play only a supporting role to much stronger covalent or dative bonds, thus making examples of exclusive cation–π bonding exceedingly rare. In this study, a neutral diboryne molecule is found to encapsulate the light alkali metal cations Li⁺ and Na⁺ in the absence of a net charge, covalent bonds, or lone‐pair donor groups. The resulting encapsulation complexes are, to our knowledge, the first structurally authenticated species in which a neutral molecule binds the light alkali metals exclusively through cation–π interactions.     

The Cation–π Interaction in Small‐Molecule Catalysis

Catalysis by small molecules (≤1000 Da, 10^?9 m) that are capable of binding and activating substrates through attractive, noncovalent interactions has emerged as an important approach in organic and organometallic chemistry. While the canonical noncovalent interactions, including hydrogen bonding, ion pairing, and π stacking, have become mainstays of catalyst design, the cation–π interaction has been comparatively underutilized in this context since its discovery in the 1980s. However, like a hydrogen bond, the cation–π interaction exhibits a typical binding affinity of several kcal mol^?1 with substantial directionality. These properties render it attractive as a design element for the development of small‐molecule catalysts, and in recent years, the catalysis community has begun to take advantage of these features, drawing inspiration from pioneering research in molecular recognition and structural biology. This Review surveys the burgeoning application of the cation–π interaction in catalysis.     

CH–π and CF–π Interactions Lead to Structural Changes of N‐Heterocyclic Carbene Palladium Complexes

The role of CH–π and CF–π interactions in determining the structure of N‐heterocyclic carbene (NHC) palladium complexes were studied using ¹H NMR spectroscopy, X‐ray crystallography, and DFT calculations. The CH–π interactions led to the formation of the cis‐anti isomers in 1‐aryl‐3‐isopropylimidazol‐2‐ylidene‐based [(NHC)₂PdX₂] complexes, while CF–π interactions led to the exclusive formation of the cis‐syn isomer of diiodobis(3‐isopropyl‐1‐pentafluorophenylimidazol‐2‐ylidene) palladium(II).