Effect of cooperativity in lithium bonding on the strength of halogen bonding and tetrel bonding: (LiCN)n···ClYF3 and (LiCN)n···YF3Cl (Y = C, Si and n = 1-5) complexes as a working model

The Tetrel Bonds of Hypervalent Halogen Compounds

The tetrel bond between PhXF2Y(TF3) (T = C and Si; X = Cl, Br, and I; Y = F and Cl) and the electron donor MCN (M = Li and Na) was investigated at the M06-2X/aug-cc-pVDZ level of theory. As the electronegativity of the halogen atom X increases, the strength of the tetrel bond also increases, but as the electronegativity of the halogen atom Y increases, the strength of the tetrel bond decreases. The magnitude of the interaction energy in most -CF3 complexes was found to be less than 10 kcal/mol, but to exceed 11 kcal/mol for PhClF2Cl(CF3)⋯NCNa. The tetrel bond is greatly enhanced when the -SiF3 group interacts with LiCN or NaCN, with the largest interaction energy approaching 100 kcal/mol and displaying a covalent Si⋯N interaction. Along with this enhancement, the Si⋯N distance was found to be less than the X-Si bond length, the -SiF3 group to be closer to the N atom, and in most -SiF3 systems, the X-Si-F angle to be less than 90°; the -SiF3 group therefore undergoes inversion and complete transfer in some systems.

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.

Promotion of TH3 (T = Si and Ge) group transfer within a tetrel bond by a cation-π interaction

The possibility of the transfer of the TH₃ group across a tetrel bond is considered by ab initio calculations. The TB is constructed by pairing PhTH₃ (Ph = phenyl; T = Si and Ge) with bases NH₃, NHCH₂, and the C₃N₂H₄ carbene. The TH₃ moves toward the base but only by a small amount in these dimers. However, when a Be²⁺ or Mg²⁺ dication is placed above the phenyl ring, the tetrel bond strength is greatly magnified reaching up to nearly 100 kcal mol^-1. This dication also induces a much higher degree of transfer which can be best categorized as half-transfer for the two N-bases and a near complete transfer for the carbene.

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.

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.

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.

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.

Strengthening of halogen bond in XCl∙∙∙FH∙∙∙F- through cooperativity with a strong hydrogen bond and proton transfer

A theoretical calculation has been performed for the ternary complexes XCl∙∙∙FH∙∙∙F^- (X = CCH, CN, OH, NC, and F) and the corresponding binary complexes. The halogen bond in the dyad is very weak with the interaction energy less than 2.5 kcal/mol. Interestingly, the halogen bond gets a big enhancement when it combines with a very strong hydrogen bond in FH∙∙∙F^-, and the largest interaction energy is up to ∼25.6 kcal/mol in FCl∙∙∙FH∙∙∙F^-. The enhancement of halogen bond not only results in a larger elongation of X-Cl bond and a bigger redshift of the bond stretch vibration but also makes the blue-shifting halogen bond in NCCl∙∙∙FH be a red-shifting one in NCCl∙∙∙FH∙∙∙F^-. The halogen bond belongs to a purely close-shell interaction in the dyad, while it becomes a partially covalent interaction in XCl∙∙∙FH∙∙∙F^- (X = OH, NC, and F) with negative energy density. In FH∙∙∙F^-, the proton is shared between the two F atoms, however, this proton transfers towards the F^- end in XCl∙∙∙FH∙∙∙F^-.

Calculation of VS,max and Its Use as a Descriptor for the Theoretical Calculation of pKa Values for Carboxylic Acids

The theoretical calculation of pKa values for Brønsted acids is a challenging task that involves sophisticated and time-consuming methods. Therefore, heuristic approaches are efficient and appealing methodologies to approximate these values. Herein, we used the maximum surface electrostatic potential (V_S,max) on the acidic hydrogen atoms of carboxylic acids to describe the H-bond interaction with water (the same descriptor that is used to characterize σ-bonded complexes) and correlate the results with experimental pKa values to obtain a predictive model for other carboxylic acids. We benchmarked six different methods, all including an implicit solvation model (water): Five density functionals and the Møller–Plesset second order perturbation theory in combination with six different basis sets for a total of thirty-six levels of theory. The ωB97X-D/cc-pVDZ level of theory stood out as the best one for consistently reproducing the reported pKa values, with a predictive power of 98% correlation in a test set of ten other carboxylic acids.

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.

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.

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.