Intermolecular Non-Covalent Carbon-Bonding Interactions with Methyl Groups: A CSD, PDB and DFT Study

Attraction versus Repulsion between Methyl and Related Groups: (CH3NHCH3)2 and (CH3SeBr2CH3)2

The starting point for this work was a set of crystal structures containing the motif of interaction between methyl groups in homodimers. Two structures were selected for which QTAIM, NCI and NBO analyses suggested an attractive interaction. However, the calculated interaction energy was negative for only one of these systems. The ability of methyl groups to interact with one another is then examined by DFT calculations. A series of (CH3PnHCH3)2 homodimers were allowed to interact with each other for a range of Pn atoms N, P, As, and Sb. Interaction energies of these C⋅⋅⋅C tetrel-bonded species were below 1 kcal/mol, but could be raised to nearly 3 kcal/mol if the C atom was changed to a heavier tetrel. A strengthening of the C⋅⋅⋅C intermethyl bonds can also be achieved by introducing an asymmetry via an electron-withdrawing substituent on one unit and a donor on the other. The attractions between the methyl and related groups occur in spite of a coulombic repulsion between σ-holes on the two groups. NBO, AIM, and NCI tools must be interpreted with caution as they can falsely suggest bonding when the potentials are repulsive.

Metal-Induced Enhancement of Tetrel Bonding. The Case of C⋅⋅⋅X-IrIII (X=Cl, Br) Tetrel Bond Involving a Methyl Group

In X-ray structures of the isomorphic mer-[IrX3(THT)(CNXyl)2] (X=Cl 1, Br 2; THT=tetrahydrothiophene; Xyl=2,6-Me2C6H3-) complexes, we revealed short intermolecular contacts between the C-atom of an isocyanide methyl group and halide ligands of another molecule. Geometrical consideration of the X-ray data and analysis of appropriate DFT studies allowed the attribution of these contacts to CMe⋅⋅⋅X-IrIII (X=Cl, Br) tetrel bond. Specifically, through the application of DFT calculations and various theoretical models, the presence of tetrel bonding interactions was validated, and the contribution of the CMe⋅⋅⋅X-IrIII interaction was assessed. The reinforcement of the tetrel bond upon the isocyanide coordination to iridium(III) is substantiated by molecular electrostatic potential (MEP) surface calculations. To distinguish the tetrel bonding characteristics of CMe⋅⋅⋅X-IrIII (X=Cl, Br) interactions from conventional hydrogen bonding, we employed multiple computational methodologies, including Natural Bond Orbital (NBO) analysis and Electron Localization Function (ELF) analysis. Additionally, Energy Decomposition Analysis (EDA) was applied to selected model systems to explore the underlying physical nature of these interactions.

Cooperativity between Intermolecular Hydrogen and Carbon Bonds in ZY···CH3CN/CH3NC···HX Trimers (ZY = H2O, H2S, HF, HCl, HBr, NH3, and H2CO; HX = HF, HCl, and HBr)

Hydrogen-bonding and carbon-bonding interactions are widespread in nature. We studied the cooperativity between these interactions in 42 trimeric complexes ZY···CH3CN/CH3NC···HX, where ZY molecules are H2O, H2S, HF, HCl, HBr, NH3, and H2CO, and HX molecules are HF, HCl, and HBr. Acetonitrile (CH3CN) and isoacetonitrile (CH3NC) act as hydrogen bond acceptors as well as carbon bond donors in these trimers. Various theoretical methods, such as electronic structure calculations, quantum theory of atoms in molecule (QTAIM), natural bond orbital (NBO), and reduced density gradient analysis, are employed to study these trimers, and the results are compared with the corresponding ZY···CH3CN/CH3NC and CH3CN/CH3NC···HX dimers. Electronic structure calculations are performed at the second-order Mo̷ller-Plesset perturbation theory using the 6-311++G(2d,2p) basis set. We show that both the interactions act synergistically in these trimers leading to an increase in their bond strength as compared to the strength in the individual dimers. The cooperative energies for these trimers are in the range of 0.69 to 3.22 kJ/mol. It is seen that the carbon bonds benefit more from the cooperativity than the hydrogen bonds. The trends of cooperativity and correlations of interaction energies and cooperative energies with relevant QTAIM and NBO parameters are reported.

Spodium bonding with noble gas atoms

The nature of the bonding between a neutral group 12 member (Zn3, Cd3 and Hg3) ring and a noble gas atom was explored using quantum chemical simulations. Natural bond orbital, quantum theory of atoms in molecules, symmetry-adapted perturbation theory, and molecular electrostatic potential surface analysis were also used to investigate the type of interaction between the noble gas atom and the metal rings (Zn3, Cd3 and Hg3). The Zn3, Cd3 and Hg3 rings are bonded to the noble gas through non-covalent interactions, which was revealed by the non-covalent interaction index. Additionally, energy decomposition analysis reveals that dispersion energy is the key factor in stabilizing these systems.

Topology of electrostatic potential and electron density reveals a covalent to non-covalent carbon-carbon bond continuum

The covalent and non-covalent nature of carbon-carbon (CC) interactions in a wide range of molecular systems can be characterized using various methods, including the analysis of molecular electrostatic potential (MESP), represented as V(r), and the molecular electron density (MED), represented as ρ(r). These techniques provide valuable insights into the bonding between carbon atoms in different molecular environments. By uncovering a fundamental exponential relationship between the distance of the CC bond and the highest eigenvalue (λv1) of V(r) at the bond critical point (BCP), this study establishes the continuum model for all types of CC interactions, including transition states. The continuum model is further delineated into three distinct regions, namely covalent, borderline cases, and non-covalent, based on the gradient, , with the bond distance of the CC interaction. For covalent interactions, this parameter exhibits a more negative value than -5.0 a.u. Å-1, while for non-covalent interactions, it is less negative than -1.0 a.u. Å-1. Borderline cases, which encompass transition state structures, fall within the range of -1.0 to -5.0 a.u. Å-1. Furthermore, this study expands upon Popelier's analysis of the Laplacian of the MED, denoted as ∇2ρ, to encompass the entire spectrum of covalent, non-covalent, and borderline cases of CC interactions. Therefore, the present study presents compelling evidence supporting the concept of a continuum model for CC bonds in chemistry. Additionally, this continuum model is further explored within the context of C-N, C-O, C-S, N-N, O-O, and S-S interactions, albeit with a limited dataset.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Beryllium bonding with noble gas atoms

Quantum chemical calculations were carried out to investigate the nature of the bonding between a neutral Be₃ ring and noble gas atom. Electronic structure calculation for these complexes was carried out at different computational levels in association with natural bond orbital, quantum theory of atoms in molecules, electron localization function, symmetry adapted perturbation theory, and molecular electrostatic potential surface analysis of Be₃ complexes. The Be atoms in the Be₃ moiety are chemically bonded to one another, with the BeBe bond dissociation energy being ~125 kJ mol^-1 . The Be₃ ring interacts with the noble gases through non-covalent interactions. The binding energies of the noble gas atoms with the Be₃ ring increases with increase in their atomic number. The non-covalent interaction index, density overlap region indicator and independent gradient model analyses reveal the presence of non-covalent inter-fragment interactions in the complexes. Energy decomposition analysis reveals that dispersion plays the major role towards stabilizing these systems.

Intermolecular noncovalent interactions with carbon in solution

One of the most familiar carbon-centered noncovalent interactions (NCIs) involving an antibonding π*-orbital situated at the Bürgi-Dunitz angle from the electron donor, mostly lone pairs of electrons, is known as n → π* interactions, and if it involves a σ* orbital in a linear fashion, then it is known as the carbon bond. These NCIs can be intra- or inter-molecular and are usually weak in strength but have a paramount effect on the structure and function of small-molecular crystals and proteins. Surprisingly, the experimental evidence of such interactions in the solution phase is scarce. It is even difficult to determine the interaction energy in the solution. Using NMR spectroscopy aided with molecular dynamics (MD) simulation and high-level quantum mechanical calculations, herein we provide the experimental evidence of intermolecular carbon-centered NCIs in solution. The challenge was to find appropriate heterodimers that could sustain room temperature thermal energy and collisions from the solvent molecules. However, after several trial model compounds, the pyridine-N-oxide:dimethyltetracyanocyclopropane (PNO-DMTCCP) complex was found to be a good candidate for the investigation. NBO analyses show that the PNO:DMTCCP complex is stabilized mainly by intermolecular n → π* interaction when a weaker carbon bond gives extra stability to the complex. From the NMR study, it is observed that the NCIs between DMTCCP and PNO are enthalpy driven with an enthalpy change of -28.12 kJ mol-1 and dimerization energy of ∼-38 kJ mol-1 is comparable to the binding energies of a conventional hydrogen-bonded dimer. This study opens up a new strategy to investigate weak intermolecular interactions such as n → π* interaction and carbon bonds in the solution phase.

Noncovalent interactions in proteins and nucleic acids: beyond hydrogen bonding and π-stacking

Understanding the noncovalent interactions (NCIs) among the residues of proteins and nucleic acids, and between drugs and proteins/nucleic acids, etc., has extraordinary relevance in biomolecular structure and function. It helps in interpreting the dynamics of complex biological systems and enzymatic activity, which is esential for new drug design and efficient drug delivery. NCIs like hydrogen bonding (H-bonding) and π-stacking have been researchers' delight for a long time. Prominent among the recently discovered NCIs are halogen, chalcogen, pnictogen, tetrel, carbo-hydrogen, and spodium bonding, and n → π* interaction. These NCIs have caught the imaginations of various research groups in recent years while explaining several chemical and biological processes. At this stage, a holistic view of these new ideas and findings lying scattered can undoubtedly trigger our minds to explore more. The present review attempts to address NCIs beyond H-bonding and π-stacking, which are mainly n → σ*, n → π* and σ → σ* type interactions. Five of the seven NCIs mentioned earlier are linked to five non-inert end groups of the modern periodic table. Halogen (group-17) bonding is one of the oldest and most explored NCIs, which finds its relevance in biomolecules due to the phase correction and inhibitory properties of halogens. Chalcogen (group 16) bonding serves as a redox-active functional group of different active sites of enzymes and acts as a nucleophile in proteases and phosphates. Pnictogen (group 15), tetrel (group 14), triel (group 13) and spodium (group 12) bonding does exist in biomolecules. The n → π* interactions are linked to backbone carbonyl groups and protein side chains. Thus, they are crucial in determining the conformational stability of the secondary structures in proteins. In addition, a more recently discovered to and fro σ → σ* type interaction, namely carbo-hydrogen bonding, is also present in protein-ligand systems. This review summarizes these grand epiphanies routinely used to elucidate the structure and dynamics of biomolecules, their enzymatic activities, and their application in drug discovery. It also briefs about the future perspectives and challenges posed to the spectroscopists and theoreticians.

Pushing the Limits of the Hydrogen Bond Enhanced Halogen Bond —The Case of the C–H Hydrogen Bond

C–H hydrogen bonds have remarkable impacts on various chemical systems. Here we consider the influence of C–H hydrogen bonds to iodine atoms. Positioning a methyl group between two iodine halogen bond donors of the receptor engendered intramolecular C–H hydrogen bonding (HBing) to the electron-rich belt of both halogen bond donors. When coupled with control molecules, the role of the C–H hydrogen bond was evaluated. Gas-phase density functional theory studies indicated that methyl C–H hydrogen bonds help bias a bidentate binding conformation. Interaction energy analysis suggested that the charged C–H donors augment the halogen bond interaction—producing a >10 kcal mol⁻¹ enhancement over a control lacking the C–H⋯I–C interaction. X-ray crystallographic analysis demonstrated C–H hydrogen bonds and bidentate conformations with triflate and iodide anions, yet the steric bulk of the central functional group seems to impact the expected trends in halogen bond distance. In solution, anion titration data indicated elevated performance from the receptors that utilize C–H Hydrogen Bond enhanced Halogen Bonds (HBeXBs). Collectively, the results suggest that even modest hydrogen bonds between C–H donors and iodine acceptors can influence molecular structure and improve receptor performance.

C–H hydrogen bonds to iodine halogen bond donors are shown to improve halogen bonding and molecular preorganization.

Enhancement of the Tetrel Bond by the Effects of Substituents, Cooperativity, and Electric Field: Transition from Noncovalent to Covalent Bond

The T⋅⋅⋅N tetrel bond (TB) formed between TX₃ OH (T=C, Si, Ge; X=H, F) and the Lewis base N≡CM (M=H, Li, Na) is studied by ab initio calculations at the MP2/aug-cc-pVTZ level. Complexes involving TH₃ OH contain a conventional TB with interaction energy less than 10 kcal/mol. This bond is substantially strengthened, approaching 35 kcal/mol and covalent character, when fluorosubstituted TF₃ OH is combined with NCLi or NCNa. Along with this enhanced binding comes a near equalization of the TB T⋅⋅⋅N and the internal T-O bond lengths, and the associated structure acquires a trigonal bipyramidal shape, despite a high internal deformation energy. This structural transformation becomes more complete, and the TB is further strengthened upon adding an electron acceptor BeCl₂ to the Lewis acid and a base to the NCM unit. This same TB strengthening can be accomplished also by imposition of an external electric field.

Computational Evaluation of Me2 TCCP as Lewis Acid

Supramolecular adducts between dimethyl-2,2,3,3-tetracyanocyclopropane (Me2 TCCP) with 21 small (polar) molecules and 10 anions were computed with DFT (B3LYP-D3/def2-TZVP). Their optimized geometries were used to obtain interaction energies, and perform energy decomposition and 'atoms-in-molecules' analyses. A set of 38 other adducts were also evaluated for comparison purposes. Selected examples were further scrutinized by inspection of the molecular electrostatic potential maps, Noncovalent Interaction index plots, the Laplacian, the orbital interactions, and by estimating the Gibbs free energy of complexation in hexane solution. These calculations divulge the thermodynamic feasibility of Me2 TCCP adducts and show that complexation is typically driven by dispersion with less polarized partners, but by orbital interactions when more polarized or anionic guests are deployed. Most Me2 TCCP adducts are more stable than simple hydrogen bonding with water, but less stable than traditional Lewis adducts involving Me3 B, or a strong halogen bond such as with Br2 . Several bonding analyses showed that the locus of interaction is found near the electron poor sp3 -hydridized (NC)2 C-C(CN)2 carbon atoms. An empty hybrid σ*/π* orbital on Me2 TCCP was identified that can be held responsible for the stability of the most stable adducts due to donor-acceptor interactions.

Methyl groups as widespread Lewis bases in noncovalent interactions

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

The Effects of Electronegativity of X and Hybridization of C on the X-C⋅⋅⋅O Interactions: A Statistical Analysis on Tetrel Bonding

Cone and distance-cone corrected statistical analyses have been performed on X-C⋅⋅⋅O (X=H, B, C, N, O and F; the C atom is sp² and sp³ hybridized) tetrel bonds. The sp³ -C and sp² -C prefer to form the interactions through σ-hole (∠XCO≈180°) and π-hole (∠XCO≈90°), respectively. With the increase in electronegativity of X, the preference for the particular angles of the respective geometries increases and the C⋅⋅⋅O distance becomes shorter. The angular preference is found to be more prominent in the cases of π-hole interactions than that in the σ-hole interactions. A similar distance-cone corrected statistical analysis on O=C⋅⋅⋅O interaction also suggests that the preferred ∠OCO angle is ∼90° and the preferred C⋅⋅⋅O distance is around the sum of van der Waals radii (3.22 Å) of the C and O atoms. However, a cone-corrected statistical analysis on X-Si⋅⋅⋅O interactions suggests that the preference for linearity in this case is much higher than that for the X-C⋅⋅⋅O σ-hole interactions.

Relative Strengths of a Pnicogen and a Tetrel Bond and Their Mutual Effects upon One Another

The ability of the T and Z atoms of TR₃ZR₂ to engage in a noncovalent interaction with NH₃ is assessed by DFT calculations, where the T atom refers to C, Si, and Ge; Z = As, Sb, and P; and substituents R = H and F. In most instances, the tetrel bond (TB) is both stronger and shorter than the pnicogen bond (ZB). These two bond strengths can be equalized, or preference shifted to the ZB, if F substituents are placed on the Z and H on the T atoms. Employing C as the T atom results in a very weak TB, with the ZB clearly favored energetically. The simultaneous formation of both TB and ZB weakens both, particularly the latter, but both bonds survive intact. Geometric and spectroscopic perturbations of the subunits reflect the two types of noncovalent bonds.

Intermolecular π-π Stacking Interactions Made Visible

Mixing the liquids hexafluorobenzene (1) and 1,3,5-trimethylbenzene (mesitylene, 2) results in a crystalline solid with a melting point of 34 °C. The solid consists of alternating π-π stacked pillars of both aromatics. This simple experiment can be used to visually demonstrate the existence and the effect of noncovalent intermolecular π-π stacking interactions. Both benzene derivatives are relatively benign and widely available, and the experiment can be performed within minutes for less than $15 when done on a 22 mL scale (total volume). The demonstration is very robust, as 1:2 mixtures in volume ratios between 2/3 and 3/2 all give a visually similar result (molar ratios of 1.8-0.8). Substituting 2 with the liquid aromatics o-xylene, p-xylene, and aniline also resulted in the formation of a crystalline solid, while using many other liquid aromatics did not.

Polycentric binding in complexes of trimethylamine-N-oxide with dihalogens

Dihalogens readily interact with trimethylamine-N-oxide under ambient conditions. Accordingly, herein, stable 1 : 1 adducts were obtained in the case of iodine chloride and iodine bromide. The crystal and molecular structure of the trimethylamine-N-oxide-iodine chloride adduct was solved. Furthermore, the geometry and electronic structure of the trimethylamine-N-oxide-dihalogen complexes were studied computationally. Only molecular ensembles were found in the global minimum for the 1 : 1 stoichiometry. The O⋯X-Y halogen bond is the main factor for the thermodynamic stability of these complexes. Arguments for electrostatic interactions as the driving force for this noncovalent interaction were discussed. Also, the equilibrium structures are additionally stabilised by weak C-H⋯X hydrogen bonds. Consequently, formally monodentate ligands are bound in a polycentric manner.

DFT and IsoStar Analyses to Assess the Utility of σ- and π-Hole Interactions for Crystal Engineering

The interpretation of 36 charge neutral 'contact pairs' from the IsoStar database was supported by DFT calculations of model molecules 1-12, and bimolecular adducts thereof. The 'central groups' are σ-hole donors (H2 O and aromatic C-I), π-hole donors (R-C(O)Me, R-NO2 and R-C6 F5 ) and for comparison R-C6 H5 (R=any group or atom). The 'contact groups' are hydrogen bond donors X-H (X=N, O, S, or R2 C, or R3 C) and lone-pair containing fragments (R3 C-F, R-C≡N and R2 C=O). Nearly all the IsoStar distributions follow expectations based on the electrostatic potential of the 'central-' and 'contact group'. Interaction energies (ΔEBSSE ) are dominated by electrostatics (particularly between two polarized molecules) or dispersion (especially in case of large contact area). Orbital interactions never dominate, but could be significant (∼30 %) and of the n/π→σ*/π* kind. The largest degree of directionality in the IsoStar plots was typically observed for adducts more stable than ΔEBSSE ≈-4 kcal⋅mol-1 , which can be seen as a benchmark-value for the utility of an interaction in crystal engineering. This benchmark could be met with all the σ- and π-hole donors studied.

On the importance of RH3C⋯N tetrel bonding interactions in the solid state of a dinuclear zinc complex with a tetradentate Schiff base ligand

The tetrel bonding and π-stacking interactions in a new dinuclear zinc complex using a tetradentate N₂O₂ donor Schiff base have been analysed energetically using DFT calculations and several computational tools.

Origins and properties of the tetrel bond

The tetrel bond (TB) recruits an element drawn from the C, Si, Ge, Sn, Pb family as electron acceptor in an interaction with a partner Lewis base. The underlying principles that explain this attractive interaction are described in terms of occupied and vacant orbitals, total electron density, and electrostatic potential. These principles facilitate a delineation of the factors that feed into a strong TB. The geometric deformation that occurs within the tetrel-bearing Lewis acid monomer is a particularly important issue, with both primary and secondary effects. As a first-row atom of low polarizability, C is a reluctant participant in TBs, but its preponderance in organic and biochemistry make it extremely important that its potential in this regard be thoroughly understood. The IR and NMR manifestations of tetrel bonding are explored as spectroscopy offers a bridge to experimental examination of this phenomenon. In addition to the most common σ-hole type TBs, discussion is provided of π-hole interactions which are a result of a common alternate covalent bonding pattern of tetrel atoms.

Understanding noncovalent bonds and their controlling forces

The fundamental underpinnings of noncovalent bonds are presented, focusing on the σ-hole interactions that are closely related to the H-bond. Different means of assessing their strength and the factors that control it are discussed. The establishment of a noncovalent bond is monitored as the two subunits are brought together, allowing the electrostatic, charge redistribution, and other effects to slowly take hold. Methods are discussed that permit prediction as to which site an approaching nucleophile will be drawn, and the maximum number of bonds around a central atom in its normal or hypervalent states is assessed. The manner in which a pair of anions can be held together despite an overall Coulombic repulsion is explained. The possibility that first-row atoms can participate in such bonds is discussed, along with the introduction of a tetrel analog of the dihydrogen bond.

Tetrel Bonding Interactions Involving Carbon at Work: Recent Advances in Crystal Engineering and Catalysis

The σ- and π-hole interactions are used to define attractive forces involving elements of groups 12–18 of the periodic table acting as Lewis acids and any electron rich site (Lewis base, anion, and π-system). When the electrophilic atom belongs to group 14, the resulting interaction is termed a tetrel bond. In the first part of this feature paper, tetrel bonds formed in crystalline solids involving sp3-hybridized carbon atom are described and discussed by using selected structures retrieved from the Cambridge Structural Database. The interaction is characterized by a strong directionality (close to linearity) due to the small size of the σ-hole in the C-atom opposite the covalently bonded electron withdrawing group. The second part describes the utilization of two allotropic forms of carbon (C60 and carbon nanotubes) as supramolecular catalysts based on anion–π interactions (π-hole tetrel bonding). This part emphasizes that the π-hole, which is considerably more accessible by nucleophiles than the σ-hole, can be conveniently used in supramolecular catalysis.

Engineering Crystals Using sp3 -C Centred Tetrel Bonding Interactions

1,1,2,2-Tetracyanocyclopropane derivatives 1 and 2 were designed and synthesized to probe the utility of sp3 -C centred tetrel bonding interactions in crystal engineering. The crystal packing of 1 and 2 and their 1,4-dioxane cocrystals is dominated by sp3 -C(CN)2 ⋅⋅⋅O interactions, has significant C⋅⋅⋅O van der Waals overlap (≤0.266 Å) and DFT calculations indicate interaction energies of up to -11.0 kcal mol-1 . A cocrystal of 2 with 1,4-thioxane reveals that the cyclopropane synthon prefers interacting with O over S. Computational analyses revealed that the electropositive C2 (CN)4 pocket in 1 and 2 can be seen as a strongly directional 'tetrel-bond donor', similar to halogen bond or hydrogen bond donors. This disclosure is expected to have implications for the utility of such 'tetrel bond donors' in molecular disciplines such as crystal engineering, supramolecular chemistry, molecular recognition and medicinal chemistry.

Observations of tetrel bonding between sp3-carbon and THF

We report the direct observation of tetrel bonding interactions between sp3-carbons of the supramolecular synthon 3,3-dimethyl-tetracyanocyclopropane (1) and tetrahydrofuran in the gas and crystalline phase. The intermolecular contact is established via σ-holes and is driven mainly by electrostatic forces. The complex manifests distinct binding geometries when captured in the crystalline phase and in the gas phase. We elucidate these binding trends using complementary gas phase quantum chemical calculations and find a total binding energy of -11.2 kcal mol-1 for the adduct. Our observations pave the way for novel strategies to engineer sp3-C centred non-covalent bonding schemes for supramolecular chemistry.

The ditetrel bond: noncovalent bond between neutral tetrel atoms

The ability of a tetrel atom to serve in the capacity of electron donor in a σ-hole noncovalent bond is tested by quantum calculations.