Tetrel bond cooperativity in open-chain (CH3CN)n and (CH3NC)n clusters (n=2–7): An ab initio study

From 'halogen' to 'tetrel' bonds: matrix isolation IR spectroscopic and quantum mechanical studies of the effect of central atom substitution in donor tetrahalogens on binary complex formation with formic acid

Binary complex formation between silicon tetrachloride (SiCl4) and formic acid (FA) has been observed in an argon matrix environment. Such complex formation manifests as spectral shifts in signature vibrations of the latter, namely the νCO, νC-O and νO-H vibrations. Quantum chemical calculations reveal that the most stable conformers of the complex involve predominantly the tetrel bond, which has been defined in existing literature as a variant of the "σ-hole" interactions. Here, regions of positive electrostatic potential on the tetrahedral face of SiCl4 act as electrophilic centers (σ-hole) to which the nucleophilic carbonyl group of FA is able to bind. Atoms-in-molecules analysis predicts a bond critical point along the non-covalent contact between the tetrel atom Si and the carbonyl oxygen on FA, corroborating the presence of the tetrel bond. The hyperconjugative interaction parameters at the binding interface obtained from Natural Bond Orbital (NBO) analysis are also consistent with such observations. Although apparently similar to SiCl4, there are noticeable differences in the binding preferences of the lower homologue carbon tetrachloride (CCl4). The binary complexes of the latter with the same FA acceptor molecule have been previously shown to involve halogen bonded, rather than tetrel bonded interactions (Banerjee and Bhattacharya, Spectrochim. Acta Mol. and Biomol. Spectrosc., 2021, 250, 119355). Such variations in the nature of non-covalent interactions of these tetrahalogens are attributed to differences in the distribution of electronic charge density surrounding the central tetrel atom, as obtained from mappings of their electrostatic potential surfaces. Our combined experimental and theoretical findings therefore provide direct evidence of the growing propensity of tetrel atoms to engage in tetrel bonding as we move lower down Group 14, and re-assert the reluctance of the smaller and more electronegative carbon atom to serve as a tetrel bond participant.

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.

Unraveling pnicogen bonding cooperativity: Insights from molecular electrostatic potential analysis

A theoretical investigation on the cooperativity of a series of binary, ternary, and quaternary complexes interconnected by pnicogen bonds has been conducted using calculations at the M06-2X/aug-cc-pVTZ level of density functional theory. By measuring changes in the molecular electrostatic potential (MESP) at the nucleus of interacting atoms in all of the complexes, it is possible to quantify the substantial reorganization of the electron density triggered by the formation of pnicogen bonds. The positive change in MESP, indicating a loss of electron density from the donor molecule in a dimer, facilitates the acceptance of electron density from a third molecule, resulting in the formation of a ternary complex with a stronger pnicogen bond compared to the one present in the binary complex. Similarly, the acceptor molecule in a dimer with a negative change in MESP showed an enhanced tendency to donate electron density to an electron-deficient third molecule. The MESP analysis provided valuable insights into the donor/acceptor characteristics of pnicogen bonds within the quaternary complexes. The proposed MESP hypotheses are consistent with the positive cooperativity observed in the pnicogen-bonded clusters. To quantify the changes in MESP, both at the donor atom (ΔVdonor ) and the acceptor atom (ΔVacceptor ), for all pnicogen bonds in the cluster, the total change in MESP (ΔΔVn ) was measured as ΔΔVn  = ∑(ΔVdonor )-∑(ΔVacceptor ). Remarkably, ΔΔVn exhibited a strong linear relationship with the sum of the bond energies of the pnicogen bonds in the cluster. This establishes the MESP analysis as a robust approach for understanding the strength and cooperative behavior of pnicogen-bonded clusters. Additionally, the MESP features provided clear evidence of pnicogen bond formation, further supporting the reliability of this approach.

Intermolecular interactions between the heavy-atom analogues of acetylene T2H2 (T = Si, Ge, Sn, Pb) and HCN

METHODS

The intermolecular interactions between the heavy-atom analogues of acetylene T₂H₂ (T = Si, Ge, Sn, Pb) and HCN have been investigated by theoretical calculations at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVDZ level.

RESULTS

The global energy minimum of T₂H₂ is the butterfly structure A, and another energy minimum is the planar structure B. Both structures A and B exhibit the dual behavior when binding with HCN. The various hydrogen bond (HB), dihydrogen bond (DB) and tetrel bond (TB) complexes can be found according to the MEP maps of T₂H₂. One TB and three HB complexes formed between structure A and HCN can be located for Si₂H₂ and Ge₂H₂. One TB, two HB and one DB complexes formed between structure A and HCN can be located for Sn₂H₂ and Pb₂H₂. Four TB and one HB complexes formed between structure B and HCN can be located for all the T₂H₂. The geometries and binding strengths of the complexes are compared and analyzed.

CONCLUSIONS

The interactions in these complexes are generally weak, and the interaction energies of these complexes range from -0.53 to -8.23 kcal/mol. The interaction energies of the TB complexes are larger than those of the corresponding HB and DB complexes for structure A···HCN systems. The relative binding strength of the four TB complexes exhibits different order for different structure B···HCN systems, which is consistent with the MEP maps of the isolated monomers.

Properties and Stabilities of Cyclic and Open Chains of Halogen Bonds

Open and cyclic chains from two to eight units of ICl and IF are constructed and examined by density functional theory (DFT) calculations. These chains contain either I···I or I···X halogen bonds (XBs) where X refers to Cl or F. The closed rings are more stable than the open chains due to the presence of an additional XB and enhanced cooperativity. This pattern is true even for most trimers where there is sizable geometric distortion in the rings. I···F rings are generally more stable than the corresponding I···I cycles as the I···F bond is stronger than I···I even in the simple dimer. However, I···I rings are comparable in energy to I···Cl. It is possible to construct I···I rings of at least as large as eight units, which are held together exclusively by XBs. On the other hand, the maximum possible size of I···X rings is 6. Red shifts are observed in the I-X stretching frequency bands, which magnify as the chain, both cyclic and open, grows longer. The NMR chemical shielding of the I atoms increases for I···I chains but diminishes when I···Cl bonds are present.

Tetrel Bonds between Phenyltrifluorosilane and Dimethyl Sulfoxide: Influence of Basis Sets, Substitution and Competition

Ab initio calculations have been performed for the complexes of DMSO and phenyltrifluorosilane (PTS) and its derivatives with a substituent of NH₃, OCH₃, CH₃, OH, F, CHO, CN, NO₂, and SO₃H. It is necessary to use sufficiently flexible basis sets, such as aug’-cc-pVTZ, to get reliable results for the Si···O tetrel bonds. The tetrel bond in these complexes has been characterized in views of geometries, interaction energies, orbital interactions and topological parameters. The electron-donating group in PTS weakens this interaction and the electron-withdrawing group prominently strengthens it to the point where it exceeds that of the majority of hydrogen bonds. The largest interaction energy occurs in the p-HO₃S-PhSiF₃···DMSO complex, amounting to −122 kJ/mol. The strong Si···O tetrel bond depends to a large extent on the charge transfer from the O lone pair into the empty p orbital of Si, although it has a dominant electrostatic character. For the PTS derivatives of NH₂, OH, CHO and NO₂, the hydrogen bonded complex is favorable to the tetrel bonded complex for the NH₂ and OH derivatives, while the σ-hole interaction prefers the π-hole interaction for the CHO and NO₂ derivatives.

Comparison for Electron Donor Capability of Carbon-Bound Halogens in Tetrel Bonds

The tetrel bond formed by HC≡CX, H2C=CHX, and H3CCH2X (X=F, Cl, Br, I) as an electron donor and TH3F (T=C, Si, Ge) was explored by ab initio calculations. The tetrel bond formed by H3CCH2X is the strongest, as high as -3.45 kcal/mol for the H3CCH2F···GeH3F dimer, followed by H2C=CHX, and the weakest bond is from HC≡CX, where the tetrel bond can be as small as -0.8 kcal/mol. The strength of the tetrel bond increases in the order of C < Si < Ge. For the H3CCH2X and HC≡CX complexes, the tetrel bond strength shows a similar increasing tendency with the decrease of the electronegativity of the halogen atom. Electrostatic interaction plays the largest role in the stronger tetrel bonds, while dispersion interaction makes an important contribution to the H2C=CHX complexes.

Noncovalent bond between tetrel π-hole and hydride

The π-hole above the plane of the X2T'Y molecule (T' = Si, Ge, Sn; X = F, Cl, H; Y = O, S) was allowed to interact with the TH hydride of TH(CH3)3 (T = Si, Ge, Sn). The resulting THT' tetrel bond is quite strong, with interaction energies exceeding 30 kcal mol-1. F2T'O engages in the strongest such bonds, as compared to F2T'S, Cl2T'O, or Cl2T'S. The bond weakens as T' grows larger as in Si > Ge > Sn, despite the opposite trend in the depth of the π-hole. The reverse pattern of stronger tetrel bond with larger T is observed for the Lewis base TH(CH3)3, even though the minimum in the electrostatic potential around the H is nearly independent of T. The THT' arrangement is nonlinear which can be understood on the basis of the positions of the extrema in the molecular electrostatic potentials of the monomers. The tetrel bond is weakened when H2O forms an OT' tetrel bond with the second π-hole of F2T'O, and strengthened if H2O participates in an OHO H-bond.

Intermolecular interactions between the heavy alkenes H2Si = TH2 (T = C, Si, Ge, Sn, Pb) and acetylene

The intermolecular interactions between the heavy alkenes H₂Si = TH₂ (T = C, Si, Ge, Sn, Pb) and C₂H₂ have been calculated at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVDZ level, and the nature of these complexes has been investigated by natural bond orbital. The four types (type-A, type-B, type-C and type-D) of complexes can be located for H₂Si = TH₂···C₂H₂ system. The complexes involving H₂Si = TH₂···C₂F₂ and H₂Si = TH₂···C₂(CN)₂ have also been examined to explore the substituent effects. Some complexes which are stable for H₂Si = TH₂···C₂H₂ system become unstable for H₂Si = TH₂···C₂F₂ or H₂Si = TH₂···C₂(CN)₂ system, while other complexes which are unstable for H₂Si = TH₂···C₂H₂ system become stable for H₂Si = TH₂···C₂F₂ or H₂Si = TH₂···C₂(CN)₂ system.

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.

Intermolecular Interactions Involving Heavy Alkenes H2Si=TH2 (T = C, Si, Ge, Sn, Pb) with H2O and HCl: Tetrel Bond and Hydrogen Bond

The intermolecular interactions between the heavy alkenes H₂Si=TH₂ (T = C, Si, Ge, Sn, Pb) and H₂O or HCl have been explored at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVDZ level. The various hydrogen bond (HB) and tetrel bond (TB) complexes can be located on the basis of molecular electrostatic potential maps of the isolated monomers. The competition between TB and HB interactions has been investigated through the relaxed potential energy surface scan. The results indicate that the HB complexes become more and more unstable relative to the TB complexes with the increase of the T atomic number, and cannot even retain as a minimum in some cases, for H₂Si=TH2···H₂O systems. In contrast, the HB complexes are generally more stable than TB complexes, and the TB complexes exhibit rather weak binding strength, for H₂Si=TH₂···HCl systems. The majority of the TB complexes formed between H₂Si=TH₂ and H₂O possesses very strong binding strength with covalent characteristics. The noncovalent TB complexes can be divided into two types on the basis of the orbital interactions: π-hole complexes, with binding angles ranging from 91 to 111°, and hybrid σ/π-hole complexes, with binding angles ranging from 130 to 165°. The interplay between different molecular interactions has been explored, and an interesting result is that the covalent TB interaction is significantly abated and becomes noncovalent because of the competitive effect.

Weak Interactions Get Strong: Synergy between Tetrel and Alkaline-Earth Bonds

Weak and strong noncovalent interactions such as tetrel bonds and alkaline-earth bonds, respectively, cooperate and get reinforced when acting together in ternary complexes of general formula RN··· SiH₃F···MY, where MY is a Be or Mg derivative and RN is a N-containing Lewis base with different hybridization patterns. Cooperativity has been studied in the optimized MP2/aug'-cc-pVTZ ternary complexes by looking at changes on geometries, binding energies, ²⁹Si NMR chemical shifts, and topological features according to the atoms in molecules theoretical framework. Our study shows that cooperativity in terms of energy is in general significant: more than 40 kJ/mol, and up to 83.6 kJ/mol in the most favorable case. The weakest the isolated interaction, the strongest the reinforcement in the ternary complex; in this sense, the tetrel bond is shortened enormously, between 0.3 and 0.6 Å. This dramatic reinforcement of the tetrel bond is also nicely reflected in the positive variations of the ²⁹Si chemical shifts in all the ternary complexes. At the same time the ternary complexes are characterized by the presence of totally planar silyl group, due to the pentacoordination of the Si atom. Both the hybridization of the N base and the geometry imposed by the alkaline-earth ligands have a strong influence on the binding energies, as they modify the donor ability of N and the Lewis acid character of the alkaline-earth metal.

Quantitative Assessment of Tetrel Bonding Utilizing Vibrational Spectroscopy

A set of 35 representative neutral and charged tetrel complexes was investigated with the objective of finding the factors that influence the strength of tetrel bonding involving single bonded C, Si, and Ge donors and double bonded C or Si donors. For the first time, we introduced an intrinsic bond strength measure for tetrel bonding, derived from calculated vibrational spectroscopy data obtained at the CCSD(T)/aug-cc-pVTZ level of theory and used this measure to rationalize and order the tetrel bonds. Our study revealed that the strength of tetrel bonds is affected by several factors, such as the magnitude of the σ-hole in the tetrel atom, the negative electrostatic potential at the lone pair of the tetrel-acceptor, the positive charge at the peripheral hydrogen of the tetrel-donor, the exchange-repulsion between the lone pair orbitals of the peripheral atoms of the tetrel-donor and the heteroatom of the tetrel-acceptor, and the stabilization brought about by electron delocalization. Thus, focusing on just one or two of these factors, in particular, the σ-hole description can only lead to an incomplete picture. Tetrel bonding covers a range of -1.4 to -26 kcal/mol, which can be strengthened by substituting the peripheral ligands with electron-withdrawing substituents and by positively charged tetrel-donors or negatively charged tetrel-acceptors.

Strong Tetrel Bonds: Theoretical Aspects and Experimental Evidence

In recent years, noncovalent interactions involving group-14 elements of the periodic table acting as a Lewis acid center (or tetrel-bonding interactions) have attracted considerable attention due to their potential applications in supramolecular chemistry, material science and so on. The aim of the present study is to characterize the geometry, strength and bonding properties of strong tetrel-bond interactions in some charge-assisted tetrel-bonded complexes. Ab initio calculations are performed, and the results are supported by the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) approaches. The interaction energies of the anionic tetrel-bonded complexes formed between XF₃M molecule (X=F, CN; M=Si, Ge and Sn) and A- anions (A-=F-, Cl-, Br-, CN-, NC- and N₃-) vary between -16.35 and -96.30 kcal/mol. The M atom in these complexes is generally characterized by pentavalency, i.e., is hypervalent. Moreover, the QTAIM analysis confirms that the anionic tetrel-bonding interaction in these systems could be classified as a strong interaction with some covalent character. On the other hand, it is found that the tetrel-bond interactions in cationic tetrel-bonded [p-NH₃(C₆H₄)MH₃]⁺···Z and [p-NH₃(C₆F₄)MH₃]⁺···Z complexes (M=Si, Ge, Sn and Z=NH₃, NH₂CH₃, NH₂OH and NH₂NH₂) are characterized by a strong orbital interaction between the filled lone-pair orbital of the Lewis base and empty BD*M-C orbital of the Lewis base. The substitution of the F atoms in the benzene ring provides a strong orbital interaction, and hence improved tetrel-bond interaction. For all charge-assisted tetrel-bonded complexes, it is seen that the formation of tetrel-bond interaction is accompanied bysignificant electron density redistribution over the interacting subunits. Finally, we provide some experimental evidence for the existence of such charge-assisted tetrel-bond interactions in crystalline phase.

Complexes of O=C=S with Nitrogen Bases: Chalcogen Bonds, Tetrel Bonds, and Other Secondary Interactions

Ab initio MP2/aug'-cc-pVTZ calculations have been carried out to investigate chalcogen-bond formation through the σ-hole at S and tetrel-bond formation through the π-hole at C in complexes of OCS with a series of nitrogen bases. The binding energies of chalcogen- and tetrel-bonded complexes with the sp-hybridized bases correlate exponentially with the N-S and N-C distances, respectively. The presence of secondary interactions between an N-H or C-H group of an sp2 -hybridized base and OCS in chalcogen-bonded complexes decreases the correlation between binding energies and the N-S distance. These secondary interactions are stronger in the tetrel-bonded complexes with the sp2 bases, particularly in the isomers of OCS:imidazole and OCS : N2 H2 , where they may be described as distorted N-H⋅⋅⋅O or N-H⋅⋅⋅S hydrogen bonds. Charge-transfer interactions are consistent with the nature of the primary and secondary interactions in these complexes. The in-plane OCS bending frequencies are blue-shift in the chalcogen-bonded complexes, and red-shifted in the tetrel-bonded complexes. EOM-CCSD spin-spin coupling constants 1c J(N4-S) across chalcogen bonds have absolute values less than 9.0 Hz, while the two-bond coupling constants 2c J(N4-C) do not exceed 4.0 Hz. These are greater in absolute value that the one-bond coupling constants 1t J(N4-C) across tetrel bonds that are less than 0.5 Hz at much shorter N-C distances.

The effect of anions on noncovalent interactions in model clusters of chalcogen-containing (CH3)2X (X = O, S, Se) molecules

A computational study of F⁻⋯(CH₃)₂O⋯CH₃F with F⁻ bound to the protons of the two methyl groups, found significant enhancement of the O⋯C interaction relative to the neutral (CH₃)₂O⋯CH₃F dyad.

Cooperativity between the hydrogen bonding and σ-hole interaction in linear NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 complexes (X = Cl, Br; Z = Ar, Kr): a comparative study

Quantum chemical calculations are performed to investigate the cooperativity of hydrogen bonding with halogen or aerogen bonding interactions in linear NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 clusters, where X = Cl, Br and Z = Ar, Kr. To understand the cooperativity mechanism in these systems, the corresponding binary NCX···NCH and O3Z···NCH complexes are also considered. The binding distances, interaction energies, and bonding properties of the NCX···(NCH)n=2–5 and O3Z···(NCH)n=2–5 clusters are analyzed in detail. It is found that the cooperative effects in the hydrogen bonding tend to strengthen X···N and Z···N interactions. For both NCX···(NCH)n and O3Z···(NCH)n clusters, a small bond shrinkage is observed from n = 4 to n = 5, which suggests that the cooperativity effects are almost saturated in the larger clusters (n > 5). As the size of the X or Z atom is increased, the magnitude of the cooperative energy in these systems is also increased, which is mainly ascribed to changes in electrostatic potentials and orbital interactions. Our results indicate that the cooperative effects lead to a substantial change in the 14N nuclear quadrupole coupling constants of the NCH molecule.

Regulation of coin metal substituents and cooperativity on the strength and nature of tetrel bonds

Ab initio calculations have been performed for the tetrel-bonded dyad MCN⋯TF₄ (M = Cu, Ag, and Au; T = C, Si, Ge, and Sn) and C₂H₄…MCN…TF₄ and C₂(CN)₄…MCN…TF₄.

Comparison of tetrel bonds in neutral and protonated complexes of pyridineTF3and furanTF3(T = C, Si, and Ge) with NH3

Protonation not only changes the primary interaction mode between α/β-furanCF₃/p-PyCF₃and NH₃but also prominently enhances the strength of the Si/Ge⋯N tetrel bond.

Modulating the strength of tetrel bonding through beryllium bonding

Quantum chemical calculations were performed to investigate the stability of the ternary complexes BeH2···XMH3···NH3 (X = F, Cl, and Br; M = C, Si, and Ge) and the corresponding binary complexes at the atomic level. Our results reveal that the stability of the XMH3···BeH2 complexes is mainly due to both a strong beryllium bond and a weak tetrel-hydride interaction, while the XMH3···NH3 complexes are stabilized by a tetrel bond. The beryllium bond with a halogen atom as the electron donor has many features in common with a beryllium bond with an O or N atom as the electron donor, although they do exhibit some different characteristics. The stability of the XMH3···NH3 complex is dominated by the electrostatic interaction, while the orbital interaction also makes an important contribution. Interestingly, as the identities of the X and M atoms are varied, the strength of the tetrel bond fluctuates in an irregular manner, which can explained by changes in electrostatic potentials and orbital interactions. In the ternary systems, both the beryllium bond and the tetrel bond are enhanced, which is mainly ascribed to increased electrostatic potentials on the corresponding atoms and charge transfer. In particular, when compared to the strengths of the tetrel and beryllium bonds in the binary systems, in the ternary systems the tetrel bond is enhanced to a greater degree than the beryllium bond. Graphical Abstract A tetrel bond can be strengthened greatly by a beryllium bond.

On the Importance of Noncovalent Carbon-Bonding Interactions in the Stabilization of a 1D Co(II) Polymeric Chain as a Precursor of a Novel 2D Coordination Polymer

A new cobalt(II) coordination polymer 2 with μ1,5 dicyanamide (dca) and a bidentate ligand 3,5-dimethyl-1-(2'-pyridyl)pyrazole (pypz) is prepared in a stepwise manner using the newly synthesized one-dimensional linear Co(II) coordination polymer 1 as a precursor. The structural and thermal characterizations elucidate that the more stable complex 2 shows a two-dimensional layer structural feature. Here, Co(II) atoms with μ1,5 dicyanamido bridges are linked by the ligand pypz forming a macrocyclic chain that runs along the crystallographic 'c' axis having 'sql' (Shubnikov notation) net topology with a 4-connected uninodal node having point symbol {4(4).6(2)}. The remarkable noncovalent carbon-bonding contacts detected in the X-ray structure of compound 1 are analyzed and characterized by density functional theory calculations and the analysis of electron charge density (atoms in molecules).