Interplay between Metal⋅⋅⋅π Interactions and Hydrogen Bonds: Some Unusual Synergetic Effects of Coinage Metals and Substituents

The ternary systems of C₂H₄ (C₂H₂ or C₆H₆)‐MCN‐HF (M=Cu, Ag, Au) and the respective binary systems were investigated to study the interplay between metal???π interactions and hydrogen bonds. The metal???π interactions in C₂H₄‐MCN become stronger with the irregular order Ag<Cu<Au, while the hydrogen bonds in MCN‐HF become weaker following the same order. The metal???π interactions are weakened as the H atoms in the π system are replaced with electron‐withdrawing groups and enhanced by electron‐donating groups. Type 1 of these ternary systems, in which MCN acts as Lewis base and acid simultaneously, is more stable than type 2, in which C₂H₄ acts as a double Lewis base. Negative cooperativity is present in type 2 ternary systems with a weakening of the metal???π interactions and the hydrogen bonds. Positive cooperativity is found in type 1 ternary systems with an enhancement of the metal???π interactions and the hydrogen bonds, except for C₂(CN)₄‐AuCN‐HF‐1. The weaker metal???π interaction in C₆H₆‐AuCN has a greater enhancing effect on the hydrogen bond in AuCN‐HF than those in C₂H₄‐AuCN and C₂H₂‐AuCN. These synergetic effects were analyzed with the natural bond orbital and energy decomposition.     

Behavior of Halogen Bonds of the Y−X⋅⋅⋅π Type (X,Y=F,Cl, Br,I) in the Benzene π System,Elucidated by Using a Quantum Theory of Atoms in Molecules Dual‐Functional Analysis

The nature of halogen bonds of the Y?X‐?‐π(C₆H₆) type (X, Y=F, Cl, Br, and I) have been elucidated by using the quantum theory of atoms in molecules (QTAIM) dual‐functional analysis (QTAIM‐DFA), which we proposed recently. Asterisks (?) emphasize the presence of bond‐critical points (BCPs) in the interactions in question. Total electron energy densities, H_b( r _c), are plotted versus H_b( r _c)?V_b( r _c)/2 [=(?²/8m)?²ρ_b( r _c)] for the interactions in QTAIM‐DFA, in which V_b( r _c) are potential energy densities at the BCPs. Data for perturbed structures around fully optimized structures were used for the plots, in addition to those of the fully optimized ones. The plots were analyzed by using the polar (R, θ) coordinate for the data of fully optimized structures with (θ_p, κ_p) for those that contained the perturbed structures; θ_p corresponds to the tangent line of the plot and κ_p is the curvature. Whereas (R, θ) corresponds to the static nature, (θ_p, κ_p) represents the dynamic nature of the interactions. All interactions in Y?X‐?‐π(C₆H₆) are classified by pure closed‐shell interactions and characterized to have vdW nature, except for Y?I‐?‐π(C₆H₆) (Y=F, Cl, Br) and F?Br‐?‐π(C₆H₆), which have typical hydrogen‐bond nature without covalency. I?I‐?‐π(C₆H₆) has a borderline nature between the two. Y?F‐?‐π(C₆H₆) (Y=Br, I) were optimized as bent forms, in which Y‐?‐π interactions were detected. The Y‐?‐π interactions in the bent forms are predicted to be substantially weaker than those in the linear F?Y‐?‐π(C₆H₆) forms.     

The nido‐Cage⋅⋅⋅π Bond: A Non‐covalent Interaction between Boron Clusters and Aromatic Rings and Its Applications

Non‐covalent interactions involving multicenter multielectron skeletons such as boron clusters are rare. Now, a non‐covalent interaction, the nido‐cage???π bond, is discovered based on the boron cluster C₂B₉H₁₂^? and an aromatic π system. The X‐ray diffraction studies indicate that the nido‐cage???π bonding presents parallel‐displaced or T‐shaped geometries. The contacting distance between cage and π ring varies with the type and the substituent of the aromatic ring. Theoretical calculations reveal that this nido‐cage???π bond shares a similar nature to the conventional anion???π or π???π bonds found in classical aromatic ring systems. This nido‐cage???π interaction induces variable photophysical properties such as aggregation‐induced emission and aggregation‐caused quenching in one molecule. This work offers an overall understanding towards the boron cluster‐based non‐covalent bond and opens a door to investigate its properties.     

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.     

Exploring the (Very Flat) Potential Energy Landscape of R−Br⋅⋅⋅π Interactions with Accurate CCSD(T) and SAPT Techniques

Halogen bonds involving an aromatic moiety as an acceptor, otherwise known as R?X???π interactions, have increasingly been recognized as being important in materials and in protein–ligand complexes. These types of interactions have been the subject of many recent investigations, but little is known about the ways in which the strengths of R?X???π interactions vary as a function of the relative geometries of the interacting pairs. Here we use the accurate CCSD(T) and SAPT2+3δMP2 methods to investigate the potential energy landscapes for systems of HBr, HCCBr, and NCBr complexed with benzene. It is found that only the separation between the complexed molecules have a strong effect on interaction strength while other geometric parameters, such as tilting and shifting R?Br???π donor relative to the benzene plane, affect these interactions only mildly. Importantly, it is found that the C_6v (T‐shaped) configuration is not the global minimum for any of the dimers investigated.     

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.     

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.     

The role of π electrons in the formation of benzene-containing lithium-bonded complexes

The intermolecular interactions in C₆H₆???LiX (X=OH, NH₂, F, Cl, Br, NC, CN) complexes are investigated by using second‐order Møller–Plesset perturbation theory (MP2) calculations and quantum theory of “atoms in molecules” (QTAIM) studies, and the role of π electrons is studied in the formation of these benzene‐containing lithium‐bonded complexes. The molecular electrostatic potentials of benzene and LiX determine the geometries of the lithium‐bonded complexes. The electron densities at the lithium bond critical points in the πC₆H₆???LiX complexes are obviously stronger than those in the σC₆H₆???LiX complexes, which indicates that the intermolecular interactions in the C₆H₆???LiX complexes are mainly attributable to π‐type interaction. The topological and energy properties at the lithium bond critical points in both the C₆H₆???LiX and πC₆H₆???LiX complexes are linear with the interaction energies, thereby showing the crucial role of the π electrons in the formation of these complexes. Electron localization function (ELF) analysis indicates that the formation of the lithium bonds leads to the reduction of the ELF π‐electron density and volume, and the reduction of the π‐electron volume is linear with the interaction energies with the correction coefficient 0.9949.     

Luminescent Benzoquinolate‐Isocyanide Platinum(II) Complexes: Effect of Pt⋅⋅⋅Pt and π⋅⋅⋅π Interactions on their Photophysical Properties

The neutral compounds [Pt(bzq)(CN)(CNR)] (R=tBu ( 1 ), Xyl ( 2 ), 2‐Np ( 3 ); bzq= benzoquinolate, Xyl=2,6‐dimethylphenyl, 2‐Np=2‐napthyl) were isolated as the pure isomers with a trans‐C_bzq,CNR configuration, as confirmed by ¹³C{¹H} NMR spectroscopy in the isotopically marked [Pt(bzq)(¹³CN)(CNR)] (R=tBu ( 1′ ), Xyl ( 2′ ), 2‐Np ( 3′ )) derivatives (δ¹³C_CN≈110 ppm; ¹J(Pt,¹³C)≈1425 Hz]. By contrast, complex [Pt(bzq)(C≡CPh)(CNXyl)] ( 4 ) with a trans‐N_bzq,CNR configuration, has been selectively isolated from [Pt(bzq)Cl(CNXyl)] (trans‐N_bzq,CNR) using Sonogashira conditions. X‐ray diffraction studies reveal that while 1 adopts a columnar‐stacked chain structure with Pt–Pt distances of 3.371(1) Å and significant π???π interactions (3.262 Å), complex 2 forms dimers supported only by short Pt???Pt (3.370(1) Å) interactions. In complex 4 the packing is directed by weak bzq???Xyl and bzq???C≡E (C, N) interactions. In solid state at room temperature, compounds 1 and 2 both show a bright red emission (?=42.1 % 1 , 57.6 % 2 ). Luminescence properties in the solid state at 77 K and concentration‐dependent emission studies in CH₂Cl₂ at 298 K and at 77 K are also reported for 1 , 1·CHCl3 , 2 , 2' , 2·CHCl3 , 3 , 4 .     

Theoretical Investigation on the Effect of Different Nitrogen Donors on Intramolecular Se⋅⋅⋅N Interactions

The effect of different donor nitrogen atoms on the strength and nature of intramolecular Se ??? N interactions is evaluated for organoselenium compounds having N,N‐dimethylaminomethyl (dime), oxazoline (oxa) and pyridyl (py) substituents. Quantum chemical calculations on three series of compounds [2‐(dime)C₆H₄SeX ( 1 a – g ), 2‐(oxa)C₆H₄SeX ( 2 a – g ), 2‐(py)C₆H₄SeX ( 3 a – g ); X=Cl, Br, OH, CN, SPh, SePh, CH₃] at the B3LYP/6‐31G(d) level show that the stability of different conformers depends on the strength of intramolecular nonbonded Se ??? N interactions. Natural bond orbital (NBO), NBO deletion and atoms in molecules (AIM) analyses suggest that the nature of the Se ??? N interaction is predominantly covalent and involves nN→σ*_{Se? X} orbital interaction. In the three series of compounds, the strength of the Se ??? N interaction decreases in the order 3 > 2 > 1 for a particular X, and it decreases in the order Cl>Br>OH>SPh≈CN≈SePh>CH₃ for all the three series 1 – 3 . However, further analyses suggest that the differences in strength of Se ??? N interaction in 1 – 3 is predominantly determined by the distance between the Se and N atoms, which in turn is an outcome of specific structures of 1 , 2 and 3 , and the nature of the donor nitrogen atoms involved has very little effect on the strength of Se ??? N interaction. It is also observed that Se ??? N interaction becomes stronger in polar solvents such as CHCl₃, as indicated by the shorter r_{Se ??? N} and higher E_{Se ??? N} values in CHCl₃ compared to those observed in the gas phase.     

Short but Weak: The Z‐DNA Lone‐Pair⋅⋅⋅π Conundrum Challenges Standard Carbon Van der Waals Radii

Current interest in lone‐pair???π (lp???π) interactions is gaining momentum in biochemistry and (supramolecular) chemistry. However, the physicochemical origin of the exceptionally short (ca. 2.8 Å) oxygen‐to‐nucleobase plane distances observed in prototypical Z‐DNA CpG steps remains unclear. High‐level quantum mechanical calculations, including SAPT2+3 interaction energy decompositions, demonstrate that lp???π contacts do not result from n→π* orbital overlaps but from weak dispersion and electrostatic interactions combined with stereochemical effects imposed by the locally strained structural context. They also suggest that the carbon van der Waals (vdW) radii, originally derived for sp³ carbons, should not be used for smaller sp² carbons attached to electron‐withdrawing groups. Using a more adapted carbon vdW radius results in these lp???π contacts being no longer of the sub‐vdW type. These findings challenge the whole lp???π concept that refers to elusive orbital interactions that fail to explain short interatomic contact distances.     

Resonance Character of Copper/Silver/Gold Bonding in Small Molecule⋅⋅⋅MX (X=F,Cl, Br,CH3, CF3) Complexes

The resonance character of Cu/Ag/Au bonding is investigated in B???M?X (M=Cu, Ag, Au; X=F, Cl, Br, CH₃, CF₃; B=CO, H₂O, H₂S, C₂H₂, C₂H₄) complexes. The natural bond orbital/natural resonance theory results strongly support the general resonance‐type three‐center/four‐electron (3c/4e) picture of Cu/Ag/Au bonding, B:M?X?B⁺?M:X^?, which mainly arises from hyperconjugation interactions. On the basis of such resonance‐type bonding mechanisms, the ligand effects in the more strongly bound OC???M?X series are analyzed, and distinct competition between CO and the axial ligand X is observed. This competitive bonding picture directly explains why CO in OC???Au?CF₃ can be readily replaced by a number of other ligands. Additionally, conservation of the bond order indicates that the idealized relationship b_B???M+b_MX=1 should be suitably generalized for intermolecular bonding, especially if there is additional partial multiple bonding at one end of the 3c/4e hyperbonded triad.     

Influence of the Homopolar Dihydrogen Bonding CH⋅⋅⋅HC on Coordination Geometry: Experimental and Theoretical Studies

The reaction of the N‐thiophosphorylated thiourea (HOCH₂)(Me)₂CNHC(S)NHP(S)(OiPr)₂ (HL), deprotonated by the thiophosphorylamide group, with NiCl₂ leads to green needles of the pseudotetrahedral complex [Ni(L‐1,5‐S,S′)₂] ? 0.5 (n‐C₆H₁₄) or pale green blocks of the trans square‐planar complex trans‐[Ni(L‐1,5‐S,S′)₂]. The former complex is stabilized by homopolar dihydrogen C?H???H?C interactions formed by n‐hexane solvent molecules with the [Ni(L‐1,5‐S,S′)₂] unit. Furthermore, the dispersion‐dominated C?H??? H?C interactions are, together with other noncovalent interactions (C?H???N, C?H???Ni, C?H???S), responsible for pseudotetrahedral coordination around the Ni^II center in [Ni(L ‐1,5‐S,S′)₂] ? 0.5 (n‐C₆H₁₄).     

C—H⋯π and π–π interactions in the supramolecular structure of 3‐methyl‐1,4‐diphenyl‐1H‐pyrazolo[3,4‐b]pyridine

The supramolecular structure of the title compound, C₁₉H₁₅N₃, is defined by π–π‐stacking and C—H?π interactions. There are no conventional hydrogen bonds in the structure.     

(η5‐Cyclopentadienyl)(3‐hydroxy‐3‐methylbut‐1‐ynyl‐κC1)(triphenylphosphine‐κP)nickel(II): novel O—H...π bonding in an organometallic molecule

In the title compound, [Ni(C₅H₅)(C₅H₇O)(C₁₈H₁₅P)], the molecule adopts the expected half‐sandwich structure with no unusual metal–ligand distances. No classical hydrogen bonds are found in the structure; instead, the OH group of the butynol unit is involved in an unusual O—H...π interaction with the C[triple‐bond]C group of an adjacent molecule. The crystal structure is further stabilized by C—H...O and C—H...π interactions, leading to an extensive network of spiral columns.     

Insight into the π‐hole···π‐electrons tetrel bonds between F2ZO (Z = C,Si, Ge) and unsaturated hydrocarbons

The intermolecular π‐hole···π‐electrons interactions between F₂ZO (Z = C, Si, Ge) molecules and unsaturated hydrocarbons including acetylene, ethylene, 1,3‐butadiene and benzene were constructed to reveal the differences of tetrel bonds forming by carbon and heavier tetrel atoms. The ab initio computation in association with topological analysis of electron density, natural bond orbital, and energy decomposition analysis demonstrate that the strength of Si···π and Ge···π tetrel bonds is much stronger than that of C···π tetrel bonds. The Si···π and Ge···π tetrel bonds exhibit covalent or partially covalent interaction nature, while the weak C···π tetrel bonds display the hallmarks of noncovalent interaction, the electrostatic interaction is the primary influencing factor. The Si···π and Ge···π interactions are determined by both the σ‐ and π‐electron densities, while the C···π interactions are dominated mainly by the π‐electron densities. The π‐hole···π‐electrons tetrel bonds are dominated by electrostatic interaction, and polarization has a comparable contribution in the Si···π and Ge···π tetrel bonds.     

Probing the CH⋅⋅⋅π Weak Hydrogen Bond in Anesthetic Binding: The Sevoflurane–Benzene Cluster

Cooperativity between weak hydrogen bonds can be revealed in molecular clusters isolated in the gas phase. Here we examine the structure, internal dynamics, and origin of the weak intermolecular forces between sevoflurane and a benzene molecule, using multi‐isotopic broadband rotational spectra. This heterodimer is held together by a primary C? H???π hydrogen bond, assisted by multiple weak C? H???F interactions. The multiple nonbonding forces hinder the internal rotation of benzene around the isopropyl C? H bond in sevoflurane, producing detectable quantum tunneling effects in the rotational spectrum.