On the Power of Geometry over Tetrel Bonds

Sigma-Hole and Lone-Pair-Hole Site-Based Interactions of Seesaw Tetravalent Chalcogen-Bearing Molecules with Lewis Bases

For the first time, sigma (σ)- and lone-pair (lp)-hole site-based interactions of SF4 and SeF4 molecules in seesaw geometry with NH3 and FH Lewis bases were herein comparatively investigated. The obtained findings from the electrostatic potential analysis outlined the emergence of sundry holes on the molecular entity of the SF4 and SeF4 molecules, dubbed the σ- and lp-holes. The energetic viewpoint announced splendid negative binding energy values for σ-hole site-based interactions succeeded by lp-hole analogues, which were found to be -9.21 and -0.50 kcal/mol, respectively, for SeF4···NH3 complex as a case study. Conspicuously, a proper concurrence between the strength of chalcogen σ-hole site-based interactions and the chalcogen's atomic size was obtained, whereas a reverse pattern was proclaimed for the lp-hole counterparts. Further, a higher preference for the YF4···NH3 complexes with elevated negative binding energy was promulgated over the YF4···FH ones, indicating the eminent role of Lewis basicity. The indications of the quantum theory of atoms in molecules generally asserted the closed-shell nature of all the considered interactions. The observation of symmetry-adapted perturbation theory revealed the substantial contributing role of the electrostatic forces beyond the occurrence of σ-hole site-based interactions. In comparison, the dispersion forces were specified to govern the lp-hole counterparts. Such emerging findings would be a gate for the fruitful forthcoming applications of chalcogen bonding interactions in crystal engineering and biological systems.

Static and Dynamical Quantum Studies of CX3-AlX2 and CSiX3-BX2 (X = F, Cl, Br) Complexes with Hydrocyanic Acid: Unusual Behavior of Strong π-Hole at Triel Center

The set of TX₃-TrX₂ (T = C, Si, Ge; Tr = B, Al, Ga; X = F, Cl, Br) molecules offers a rather unique opportunity to study both σ-hole and π-hole dimerization on the tetrel and triel ends, respectively. According to the molecular electrostatic potential (MEP) distribution, the π-hole extrema (acidic sites) were more intense than their σ-hole counterparts. The molecules owning the most (CX₃-AlX₂) and least (SiX₃-BX₂) intense π-holes were chosen to evaluate their capacities to attract one and two HCN molecules (Lewis bases). We discovered that the energetic characteristics of π-hole dimers severely conflict with the monomers MEP pattern since the weakest π-hole monomer forms a dimer characterized by interaction energy compared to those created by the monomers with noticeably greater power in the π-hole region. This outcome is due to the deformation of the weakest π-hole donor. Furthermore, the MEP analysis for monomers in the geometry of respective dimers revealed a "residual π-hole" site that was able to drive second ligand attachment, giving rise to the two "unusual trimers" examined further by the NCI and QTAIM analyses. Apart from them, the π-hole/π-hole and σ-hole/π-hole trimers have also been obtained throughout this study and described using energetic and geometric parameters. The SAPT approach revealed details of the bonding in one of the "unusual trimers". Finally, Born-Oppenheimer Molecular Dynamics (BOMD) simulations were carried out to investigate the time evolution of the interatomic distances of the studied complexes as well as their stability.

R•-hole interactions of group IV-VII radical-containing molecules: A comparative study

For the first time, the potentiality of the sp²-hybridized group IV-VII radical (R^•)-containing molecules to participate in R^•-hole interactions was comparatively assessed using ^•SiF_3,^•POF₂, ^•SO₂F, and ^•ClO₃ models in the trigonal pyramidal geometry. In that spirit, a plethora of quantum mechanical calculations was performed at the MP2/aug-cc-pVTZ level of theory. According to the results, all the investigated R^•-containing molecules exhibited potent versatility to engage in R^•-hole … Lewis base interactions with significant negative binding energies for the NCH-based complexes. The strength of R^•-hole interactions was perceived to obey the ^•ClO₃ … > ^•SO₂F … > ^•POF₂ … > ^•SiF₃ … Lewis base order, outlining an inverse correlation between the binding energy and the atomic size of the R^•-hole donor. Benchmarking of the binding energy at the CCSD/CBS(T) computational level was executed for all the explored interactions and addressed an obvious similarity between the MP2 and CCSD energetic findings. QTAIM analysis critically unveiled the closed-shell nature of the explored R^•-hole interactions. SAPT-EDA proclaimed the reciprocal contributions of electrostatic and dispersion forces to the total binding energy. These observations demonstrate in better detail the nature of R^•-hole interactions, leading to a convincing amelioration for versatile fields relevant to materials science and drug design.

Yet another perspective on hole interactions

Hole interactions are known by different names depending on the key atom of the bond (halogen bond, chalcogen bond, hydrogen bond, etc.), and the geometry of the interaction (σ if in line, π if perpendicular to the Lewis acid plane). However, its origin starts with the creation of a Lewis acid by an underlying covalent bond, which forms an electrostatic depletion and a virtual antibonding orbital, which can create non-covalent interactions with Lewis bases. In this (maybe subjective) perspective, we will claim that hole interactions must be defined via the molecular orbital origin of the molecule. Under this premise we can better explore the richness of such bonding patterns. For that, we will study old, recent and new systems, trying to pinpoint some misinterpretations that are often associated with them. We will use as exemplars the triel bonds, a couple of metal complexes, a discussion on convergent σ-holes, and many cases of anti-electrostatic hole interactions.

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.

Competition between Inter and Intramolecular Tetrel Bonds: Theoretical Studies Complemented by CSD Survey

Crystal structures document the ability of a TF₃ group (T=Si, Ge, Sn, Pb) situated on a naphthalene system to engage in an intramolecular tetrel bond (TB) with an amino group on the adjoining ring. Ab initio calculations evaluate the strength of this bond and evaluate whether it can influence the ability of the T atom to engage in a second, intermolecular TB with another nucleophile. A very strong CN^- anionic base can approach the T either along the extension of a T-C or T-F bond and form a strong TB with an interaction energy approaching 100 kcal/mol, although this bond is weakened a bit by the presence of the internal T⋅⋅⋅N bond. The much less potent NCH base engages in a correspondingly longer and weaker TB, less than 10 kcal/mol. Such an intermolecular TB is weakened by the presence of the internal TB, to the point that it only occurs for the two heavier tetrel atoms Sn and Pb.

Experimental and Theoretical Studies of Dimers Stabilized by Two Chalcogen Bonds in the Presence of a N···N Pnicogen Bond

The structure of the 5,6-dichloro-2,1,3-benzoselenadiazole homodimer, obtained by adding the ligand, 4,5-dichloro-o-phenylenediamine, to the methanolic solution of SeCl₄, was determined by X-ray crystallography, augmented by Fourier transform infrared, Raman, and NMR spectroscopy. The binding motif involves a pair of Se···N chalcogen bonds, with a supplementary N···N pnicogen bond. Quantum calculations provide assessments of the strengths of the individual interactions as well as their contributing factors. All together, these three bonds compose a total interaction energy between 5.4 and 16.8 kcal/mol, with the larger chalcogen atom associated with the strongest interactions. Replacement of the Se atoms by S and Te analogues allows analysis of the dependence of these forces on the nature of the chalcogen atom. Calculations also measure the importance to the binding of the presence of a second N atom on each diazole unit as well as the substituted phenyl ring to which it is fused.

Carbon (sp3) tetrel bonding mediated BODIPY supramolecular assembly via unprecedented synergy of Csp3⋯N and Csp3⋯F pair interactions

Here, we present the first example of sp³ hybridized carbon centered (C_sp3) tetrel bonding mediated 3D BODIPY assembly via the exceptional synergy of C_sp3⋯N and C_sp3⋯F pair interactions.

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.

Dichotomous Si-H Bond Activation by Alkoxide and Alcohol in Base-Catalyzed Dehydrocoupling of Silanes

The activation of silanes in dehydrogenative coupling with alcohols under general base catalysis was studied experimentally (using multinuclear NMR, IR, and UV-visible spectroscopies) and computationally (at DFT M06/6-311++G(d,p) theory level) on the example of Ph_4-nSiH_n (n = 1-3) interaction with (CF₃)₂CHOH in the presence of Et₃N. The effect of the phenyl groups' number and H^- substitution by the electron-withdrawing (CF₃)₂CHO^- group on Si-H bond hydricity (quantified as hydride-donating ability, HDA) and Lewis acidity of silicon atom (characterized by maxima of molecular electrostatic potential) was accessed. Our results show the coordination of Lewis base (Y = Me₃N, ROH, OR^-) leads to the increased hydricity of pentacoordinate hypervalent Ph_4-nSi(Y)H_n complexes and a decrease of the reaction barrier for H₂ release. The formation of tertiary complexes [Ph_4-nSi(Y)H_n]···HOR is a critical prerequisite for the dehydrocoupling with alkoxides being ideal activators. The latter can be external or internal, generated by in situ HOR deprotonation. The mutual effect of tetrel interaction and dihydrogen bonding in tertiary complexes (RO^-)Ph_4-nSiH_n···HOR leads to dichotomous activation of Si-H bond promoting the proton-hydride transfer and H₂ release.

Nature of the Interaction of Pyridines with OCS. A Theoretical Investigation

Ab initio calculations were carried out to investigate the interaction between para-substituted pyridines (X-C₅H₄N, X=NH₂, CH₃, H, CN, NO₂) and OCS. Three stable structures of pyridine.OCS complexes were detected at the MP2=full/aug-cc-pVDZ level. The A structure is characterized by N…S chalcogen bonds and has binding energies between −9.58 and −12.24 kJ/mol. The B structure is bonded by N…C tetrel bond and has binding energies between −10.78 and −11.81 kJ/mol. The C structure is characterized by π-interaction and has binding energies between −10.76 and −13.33 kJ/mol. The properties of the systems were analyzed by AIM, NBO, and SAPT calculations. The role of the electrostatic potential of the pyridines on the properties of the systems is outlined. The frequency shift of relevant vibrational modes is analyzed.

Dual Geometry Schemes in Tetrel Bonds: Complexes between TF₄ (T = Si, Ge, Sn) and Pyridine Derivatives

When an N-base approaches the tetrel atom of TF₄ (T = Si, Ge, Sn) the latter molecule deforms from a tetrahedral structure in the monomer to a trigonal bipyramid. The base can situate itself at either an axial or equatorial position, leading to two different equilibrium geometries. The interaction energies are considerably larger for the equatorial structures, up around 50 kcal/mol, which also have a shorter R(T··N) separation. On the other hand, the energy needed to deform the tetrahedral monomer into the equatorial structure is much higher than the equivalent deformation energy in the axial dimer. When these two opposite trends are combined, it is the axial geometry which is somewhat more stable than the equatorial, yielding binding energies in the 8–34 kcal/mol range. There is a clear trend of increasing interaction energy as the tetrel atom grows larger: Si < Ge < Sn, a pattern which is accentuated for the binding energies.

Tetrel Bond between 6-OTX₃-Fulvene and NH₃: Substituents and Aromaticity

Carbon bonding is a weak interaction, particularly when a neutral molecule acts as an electron donor. Thus, there is an interesting question of how to enhance carbon bonding. In this paper, we found that the –OCH₃ group at the exocyclic carbon of fulvene can form a moderate carbon bond with NH₃ with an interaction energy of about −10 kJ/mol. The –OSiH₃ group engages in a stronger tetrel bond than does the –OGeH₃ group, while a reverse result is found for both –OSiF₃ and –OGeF₃ groups. The abnormal order in the former is mainly due to the stronger orbital interaction in the –OSiH₃ complex, which has a larger deformation energy. The cyano groups adjoined to the fulvene ring not only cause a change in the interaction type, from vdW interactions in the unsubstituted system of –OCF₃ to carbon bonding, but also greatly strengthen tetrel bonding. The formation of tetrel bonding has an enhancing effect on the aromaticity of the fulvene ring.