Next Article in Journal
Conformational and Supramolecular Aspects in Chirality of Flexible Camphor-Containing Schiff Base as an Inducer of Helical Liquid Crystals
Next Article in Special Issue
On the Nature of the Partial Covalent Bond between Noble Gas Elements and Noble Metal Atoms
Previous Article in Journal
Bioactive Molecules Derived from Plants in Managing Dengue Vector Aedes aegypti (Linn.)
Previous Article in Special Issue
Unicorns, Rhinoceroses and Chemical Bonds
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 5
10.3390/molecules28052385
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Substituent Effects in Tetrel Bonds Involving Aromatic Silane Derivatives: An ab initio Study

by
Sergi Burguera
,
Antonio Frontera
and
Antonio Bauzá
*
Department of Chemistry, Universitat de les Illes Balears, Ctra. de Valldemossa km 7.5, 07122 Palma de Mallorca, Baleares, Spain
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(5), 2385; https://doi.org/10.3390/molecules28052385
Submission received: 17 February 2023 / Revised: 2 March 2023 / Accepted: 3 March 2023 / Published: 5 March 2023
(This article belongs to the Special Issue Fundamental Aspects of Chemical Bonding)
Browse Figures
Versions Notes

Abstract

:
In this manuscript substituent effects in several silicon tetrel bonding (TtB) complexes were investigated at the RI-MP2/def2-TZVP level of theory. Particularly, we have analysed how the interaction energy is influenced by the electronic nature of the substituent in both donor and acceptor moieties. To achieve that, several tetrafluorophenyl silane derivatives have been substituted at the meta and para positions by several electron donating and electron withdrawing groups (EDG and EWG, respectively), such as –NH2, –OCH3, –CH3, –H, –CF3 and –CN substituents. As electron donor molecules, we have used a series of hydrogen cyanide derivatives using the same EDGs and EWGs. We have obtained the Hammett’s plots for different combinations of donors and acceptors and in all cases we have obtained good regression plots (interaction energies vs. Hammet’s σ parameter). In addition, we have used the electrostatic potential (ESP) surface analysis as well as the Bader’s theory of atoms in molecules (AIM) and noncovalent interaction plot (NCI plot) techniques to further characterize the TtBs studied herein. Finally, a Cambridge Structural Database (CSD) inspection was carried out, retrieving several structures where halogenated aromatic silanes participate in tetrel bonding interactions, being an additional stabilization force of their supramolecular architectures.

Graphical Abstract

1. Introduction

Noncovalent interactions (NCIs) are of crucial importance in molecular recognition and set the foundation of modern chemistry [1,2,3,4]. From the beginning of the century, NCIs have earned a general recognition that has moved the field far beyond classical and ubiquitous forces such as hydrogen bonding and π–π stacking interactions. Concretely, the σ-/π-hole crown of interactions [5], composed by halogen [6], chalcogen [7], pnictogen [8], tetrel [9] and aerogen bonds [10], has become a reliable and robust resource for driving molecular recognition [11,12] and self-assembly processes [13,14,15] with the assistance of computational chemistry [16,17], which is crucial for achieving a full exploitation of the NCIs in several fields (e.g., materials science, catalysis, biochemistry, etc.) [18,19,20,21,22,23].
The number of Investigations in NCIs [24,25] encompassing elements from groups 14–17 have exponentially grown during the past decade, leading to the creation and exploitation of novel fields of research [26,27,28]. Several computational studies [5,6,7,8,9,10,24,25] have demonstrated the anisotropic nature of the electron density distribution in Halogen (Ha), Chalcogen (Ch) and Pnictogen (Pn) elements owing to the coexistence of lone pairs (LPs) and σ–holes on the same atom. More precisely, in monovalent Ha, divalent Ch, trivalent Pn and tetravalent Tt atoms, the LP number goes from three to none, while the number of σ-holes goes from one to four (moving from group 17 to group 14). This is important, since the number and position of the positive and negative molecular electrostatic surface potentials on these atoms is closely related to the number and location of the covalent bonds the atoms are involved in, as well as the location of the LPs.
Among this number of NCIs, tetrel bonds (TtBs), which involve a σ-hole located on an element from group IV and a Lewis base, were theoretically described by the groups of Frontera [29] and Arunan [30] and have nowadays achieved great recognition in the fields of supramolecular chemistry and crystal engineering, as well as in biological systems [31,32,33,34,35,36,37,38,39]. In spite of this, Johnson and co-workers first showed the stability of carbon-Lewis base adducts [40]. They demonstrated that the geometrical disposition of the [H2O–CO2] dimer was ruled by a C⋯O interaction that was stronger than the dimer based on a H⋯O hydrogen bond (HB). This finding was confirmed by the group of Klemperer et al. several years later via microwave spectral analysis [41]. In more detail, they demonstrated the prevalence of a supramolecular complex featuring a O2C⋯OH2 tetrel bond instead of a hydrogen bond (HO–H⋯O=C=O). Simple calculations at B97D/6-31 + G* yielded a more stable tetrel bonded complex by 1.8 kcal/mol compared to the HB geometry. In addition, it was also demonstrated during the 1980s that other HB donor molecules, such as HBr [42] and HCN [43], formed stronger tetrel bonding complexes with carbon dioxide.
Later in 2001, Alkorta and collaborators calculated a series of supramolecular complexes between Si derivatives (SiX4, X = halogen) and lone pair donors (Z = NH3, H2O, etc.). The Si⋯Z interaction distances obtained from the calculations range between 2.1 and 4.1 Å, while the strength of the interaction reached up to −10 kcal/mol for the stronger complexes [44]. In a parallel work, Alkorta evaluated the conformational trend of aminopropylsilanes [45] by using NMR data from experiments (15N and 29Si chemical shifts and JN–Si coupling constants) in combination with theoretical calculations. The results obtained provided evidence of the existence of an equilibrium between an open chain structure (dominated by the entropic component) and the supramolecular cycle (governed by the enthalpic term). Finally, the energetic and geometric features of TtBs in neutral and protonated RTtF3 (R = pyridinyl and furanyl) systems have been investigated by Scheiner and collaborators [46] for C, Si and Ge, using NH3 as a Lewis base.
The physical nature of the interaction is based on two main factors [47,48]. First, the polarizability of the tetrel atom (Tt), which increases upon descending in the group. Consequently, the electropositive region of the σ–hole increases if the EWG–Tt bond (EWG: electron withdrawing group) is more polarized, resulting in a strengthening of the NCI. Secondly, another way to polarize the EWG–Tt bond is by increasing the electron-withdrawing ability of the EWG. Therefore, the combination of heavy elements and a strong EWG increases the positive potential and size of the σ–hole, thus reinforcing the NCI (by increasing the contribution of electrostatics). In terms of steric demands, TtBs are a world apart compared to other σ-hole interactions. For instance, while the approaching LP donor molecule in a typical halogen bond (HalB) complex is usually separated by 180° from the EWG covalently attached to the halogen atom (with the minimal amount of steric repulsion), the tetrahedral spacing of the four substituents bound to the Tt atom dramatically increases the steric crowding and, at the same time, augments the steric repulsion with the Lewis base.
In this context, silicon is probably more prone to participate in noncovalent TtBs than the heavier tetrels due to their tendency to expand their valence, thus engaging in covalent/coordination chemistry. Moreover, tetrahedral lead(IV) compounds are quite uncommon, and lead(II) has a rich coordination chemistry [49,50,51,52] and it is considered a metal. Nonetheless, hypervalent species of Si are also known [53,54,55,56] although less abundant. The aim of this study is to investigate the geometric and energetic features of tetrel bonded complexes involving tetrafluorophenyl silane moieties and to study substituent effects in the meta and para positions of the ring. A similar study encompasses the work from Franconetti and collaborators [57], where substituent effects on a series of “like–like” Sn···Sn Tetrel bond complexes were evaluated at the RI-MP2/def2-TZVP level of theory. The authors used para-substituted phenylstannane derivatives as TtB donors and several tin-containing benzene derivatives as TtB acceptor species. In addition, the study from An and co-workers [58] shed light on the basis set influence (using double and triple zeta Dunning’s basis sets), substitution and competition in Tetrel bonding interactions involving phenyltrifluorosilane and dimethyl sulfoxide molecules. Building upon these previous works, we were interested in (i) evaluating meta-substitution effects on the aromatic silane moiety, (ii) including substitution on the electron donor molecule and (iii) inspecting the Cambridge Structural Database (CSD) [59] for specific examples of Si TtBs involving aromatic silanes.
To achieve that, we have used several electron-donating and electron-withdrawing groups (EDG and EWG, respectively), such as –NH2, –OCH3, –CH3, –H, –CF3 and –CN substituents, incorporated at the meta and para positions of a fluorinated aromatic silane derivative (see Figure 1). As electron donor molecules, we have used a series of hydrogen cyanide derivatives using the same EDGs and EWGs (see ESI for cartesian coordinates of complexes 1 to 72 in supplementary material). In addition, the topological analysis of the electron density was carried out using the quantum theory of atoms in molecules (QTAIM) and noncovalent interactions plot (NCI plot) methodologies. Finally, a CSD inspection was performed, revealing some examples where halogenated aromatic silanes undergo TtBs in their solid-state structure, thus giving reliability to the results obtained from the calculations. We expect that the evidence reported herein will be useful for those scientists devoted to the fields of supramolecular chemistry and crystal engineering.

2. Results and Discussion

2.1. ESP Analysis

We started our study by carrying out an electrostatic potential surface (ESP) analysis on the tetrel bond donor and acceptor molecules used herein (see Figure 2 for some representative examples). Compounds involving –CN, –H and –NH2 substituents are discussed in more detail while the complete list of ESP values is gathered in Table 1. Firstly, in the case of the three hydrogen cyanide derivatives, the basicity of the sp N atom increased continually from –CN to –NH2 substituents, following their electron-donating or electron-withdrawing nature. Secondly, in the case of the TtB donor moieties four σ-holes were found over the Si atom, due to the four established covalent bonds. Our interests focused on the σ-hole belonging to the C–Si bond, which is depicted in Figure 2, thus, the presence of a EWG or EDG group resulted in an increase or a decrease of the ESP value over the Si σ-hole, respectively. This was observed in both the meta- and para-substituted aromatic silanes, achieving the latter slightly more positive (meta-CN (+34.5 kcal/mol) and para-CN (+35.8 kcal/mol)) or less positive (meta-NH2 (+23.8 kcal/mol) and para-NH2 (+21.9 kcal/mol)) ESP values. These results are in line with previous studies on the substituent effects in halogen bonding interactions [36], pointing out that resonance effects are not an important factor affecting the Si σ-hole’s potential and that TtB complexes involving either meta- or para-substituted silanes are expected to show a similar strength. On the other hand, using –H as substituent (either on the TtB donor or acceptor molecule) middle point results were reported, as expected. Finally, the rest of the ESP values included in Table 1 follow the expected trends for both the TtB donor and acceptor molecules, depending on their electron-donating or electron-withdrawing nature.

2.2. Energetic Study

The values of the binding energies for all complexes of the study are included in Table 2 and Table 3, following the same behaviour independently of the position of the substituent in the aromatic ring. Firstly, in all the cases negative and attractive interaction energy values were obtained, ranging between −1.0 and −3.6 kcal/mol. If one focuses on the same electron donor molecule, those complexes involving EWGs (–CN and –CF3) resulted in more favourable binding energy values than those involving EDGs (–NH2, –OCH3 and –CH3), as expected from the ESP analysis discussed above. Besides, the position of the substituent on the ring (either meta or para) does not have a noticeable influence on the strength of the interaction, which is also in agreement with the ESP analysis shown above, indicating that the resonance effects do not play an important role in the stabilization of the Si tetrel bond complexes studied herein. Figure 3 shows a graphical representation of the ESP values vs. the binding energy values involving meta- and para-substituted aromatic silanes.
As noted, in both cases we found excellent correlations, thus reinforcing the ESP as a powerful predicting tool of the strength of these Si Tetrel bond complexes. Finally, in Figure 4 some optimized geometries at the RI-MP2/def2-TZVP level of theory are shown for six representative complexes exhibiting both meta and para substitution. In all the cases, a bond critical point (BCP) connects the sp N atom from the electron donor molecule with the Si atom of the aromatic silane moiety. In addition, a greenish isosurface can also be found along the bond path that connects both molecules, indicating the presence of the tetrel bonding interaction. The value of the density at the bond critical point that characterizes the TtBs is also given in Figure 4, while the rest of the values are also gathered in Table 2 and Table 3.

2.3. Hammett’s Representations

Figure 5 and Figure 6 gather several Hammett regression plots, where the σm and σp constants were plotted against the binding energies of those meta and para complexes involving –CN, –H and –NH2 as electron donor substituents. As noted, very good correlations were found in those representations that separately involve the meta- and para-substituted complexes (Figure 5, top). In addition, if these two plots were combined (Figure 5, bottom) the correlation was still very good, with R values of 0.918 (–CN), 0.856 (–H) and 0.874 (–H), thus showing that resonance effects do not play a significant role in the substituent effects involving these systems, similarly to previously published results [60].

2.4. CSD Search

As the last stage of our study, the CSD was inspected (December 2022) and five structures were found exhibiting tetrel bonding interactions between halogenated aromatic silane moieties. Three of them are shown in Figure 6, which correspond to AHESIR [61] and WIZHUI01 [62] X-ray structures. Firstly, the AHESIR structure corresponds to (E) 4,4′-Dimethyl-2,2′-disilylazobenzene, which exhibits discrete monomeric units interacting through several NCIs (e.g., π-π stacking, CH-π) on its solid state architecture. In addition to these, Si TtBs are undertaken between a F atom from the –SiF3 moiety and the Si–F σ-hole from a neighbouring unit. The WIZHUI01 structure involves a polymorph of the compound 2,6-Mes2C6H3SiF3, which shows discrete units self-assembling through F···Si TtBs bonds in the solid state architecture, forming a one-dimensional columnar pattern. In Figure 6 the interaction strength as well as the TtB distances are indicated, which lie within the same range as the results obtained using fully optimized models, thus giving reliability to the binding energies obtained.

3. Conclusions

The above computations indicate the ability of tetrafluorophenyl silanes to engage in TtB interactions. The substituent effects were evaluated by using a series of EDGs and EWGs on the meta and para positions of the ring and on the electron donor molecule and by correlating the Hammett substituent constants with the strength of the interaction, obtaining excellent correlations in all cases. These indicate that resonance effects do not play a crucial role and that substituent effects can be explained on the basis of inductive effects. In addition, a survey of the CSD revealed several X-ray structures where TtBs between halogenated aromatic silane moieties were undertaken, contributing to their stabilization of the supramolecular architectures. We expect that the results reported herein will be useful for those scientists devoted to the fields of supramolecular chemistry and crystal engineering of silicon-based compounds.

4. Materials and Methods

4.1. General Considerations

The interaction energies of all the complexes used were calculated at the RI-MP2 [63]/def2-TZVP [64] level of theory. This level of theory has achieved success to accurately represent interaction energies involving both neutral and charged electron donors [65]. The calculations were performed using the TURBOMOLE program, version 7.0. [66]. The interaction energies (ΔEBSSE, in kcal/mol) were calculated as the energy difference between the complex and the isolated monomers following the supermolecule approximation (ΔEcomplex = Ecomplex − EmonomerA − EmonomerB). These values were corrected using the Boys and Bernardi counterpoise correction method [67]. Frequency calculations demonstrated the true minima nature of the TtB complexes studied herein. We have also calculated the Gibbs free energies for two representative complexes (18 and 54) involving –CH3 and –CN groups and resulting in positive values, thus pointing to the importance of entropic effects when dealing with this family of compounds. In the case of the X-ray crystal models of AHESIR and WIZHUI01 structures, single point calculations were carried out to obtain the interaction energy values at the RI-MP2/def2-TZVP level of theory. The MEP surfaces were computed using Gaussian 16 software [68]. The Bader’s AIM theory [69] was used to analyse and describe the interactions discussed in this work using the AIMall calculation package [70]. The RI-MP2/def2-TZVP level of theory was also used for the wavefunction analysis (also using Gaussian 16 software). The NCI plot [71] isosurfaces acknowledge the presence of both attractive and repulsive interactions, as denoted by the sign of the second density Hessian eigenvalue and characterized by the isosurface colour. The colour scheme is composed of a red–yellow–green–blue scale using red for repulsive (λ+cut) and blue for attractive (λ−cut) NCI interaction density. Weak repulsive and weak attractive interactions are identified by yellow and green surfaces, respectively.

4.2. CSD Survey

The CSD was inspected in December 22 (February 2021 update) using the Conquest program [72], retrieving five structures where tetrel bonds between halogenated aromatic silane moieties were undertaken. The following CSD codes from the structures were retrieved from the search: ASEHIR, EHIFAF, GOFSUR, PIDMOH and WIZHUI01. The geometrical criteria used to classify a noncovalent contact as TtB are as follows (see also Figure 7):
  • dSi···A ≤ sum of vdW radii + 0.5 Å.
  • C/X–Si···A between 160 and 180 degrees (X = C, F, Cl, Br and I and A = any atom).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28052385/s1, Cartesian coordinates of complexes 1 to 72.

Author Contributions

Conceptualization, A.B. and A.F.; methodology, S.B., A.B. and A.F.; formal analysis, S.B. and A.B.; investigation, S.B., A.F. and A.B; writing—original draft preparation, A.B. and A.F.; writing—review and editing, A.F. and A.B.; project administration, A.F.; funding acquisition, A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by MICIU/AEI, grant number PID2020-115637GB-I00 FEDER funds.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data to reproduce the results is gathered in the ESI in supplementary material.

Acknowledgments

We thank the CTI (UIB) for the computational facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lehn, J.M. Supramolecular Chemistry: Concepts and Perspectives, 1st ed.; Wiley VCH: Weinheim, Germany, 1995. [Google Scholar]
  2. Schneider, H.J. Supramolecular Systems in Biomedical Fields, 1st ed.; RSC Publishing: Cambridge, UK, 2013. [Google Scholar]
  3. Steed, A.W.; Atwood, J.L. Supramolecular Chemistry, 1st ed.; John Wily & Sons, Ltd.: Chichester, UK, 2009. [Google Scholar]
  4. Cragg, P.J. Supramolecular Chemistry: From Biological Inspiration to Biomedical Applications, 1st ed.; Springer: Dordrecht, The Netherlands, 2010. [Google Scholar]
  5. Politzer, P.; Murray, J. σ-Hole Interactions: Perspectives and Misconceptions. Crystals 2017, 7, 212. [Google Scholar] [CrossRef] [Green Version]
  6. Desiraju, G.R.; Ho, P.S.; Kloo, L.; Legon, A.C.; Marquardt, R.; Metrangolo, P.; Politzer, P.; Resnati, G.; Rissanen, K. Definition of the halogen bond (IUPAC Recommendations 2013). Pure Appl. Chem. 2013, 85, 1711–1713. [Google Scholar] [CrossRef]
  7. Aakeroy, C.B.; Bryce, D.L.; Desiraju, G.R.; Frontera, A.; Legon, A.C.; Nicotra, F.; Rissanen, K.; Scheiner, S.; Terraneo, G.; Metrangolo, P.; et al. Definition of the chalcogen bond (IUPAC Recommendations 2019). Pure Appl. Chem. 2019, 91, 1889–1892. [Google Scholar] [CrossRef]
  8. Zahn, S.; Frank, R.; Hey-Hawkins, E.; Kirchner, B. Pnicogen bonds: A new molecular linker? Chem. Eur. J. 2011, 17, 6034–6038. [Google Scholar] [CrossRef]
  9. Bauzá, A.; Mooibroek, T.J.; Frontera, A. Tetrel Bonding Interaction: Rediscovered Supramolecular Force? Angew. Chem. Int. Ed. 2013, 52, 12317–12321. [Google Scholar] [CrossRef] [PubMed]
  10. Bauzá, A.; Frontera, A. σ/π-Hole noble gas bonding interactions: Insights from theory and experiment. Coord. Chem. Rev. 2020, 404, 213112. [Google Scholar] [CrossRef]
  11. Bauzá, A.; Frontera, A.; Mooibroek, T.J. NO3− anions can act as Lewis acid in the solid state. Nat. Commun. 2017, 8, 14522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. García-Llinás, X.; Bauzá, A.; Seth, S.K.; Frontera, A. Importance of R–CF3···O Tetrel Bonding Interactions in Biological Systems. J. Phys. Chem. A 2017, 28, 5371–5376. [Google Scholar] [CrossRef] [PubMed]
  13. Bauzá, A.; Sharko, A.V.; Senchyk, G.A.; Rusanov, E.B.; Frontera, A.; Domasevitch, K.V. π–hole interactions at work: Crystal engineering with nitro-derivatives. Cryst. Eng. Comm. 2017, 19, 1933–1937. [Google Scholar] [CrossRef]
  14. Bauzá, A.; Frontera, A.; Mooibroek, T.J. π-hole interactions involving nitro compounds: Directionality of nitrate esters. Cryst. Growth Des. 2016, 16, 5520–5524. [Google Scholar] [CrossRef]
  15. Hazari, A.; Das, L.K.; Bauzá, A.; Frontera, A.; Ghosh, A. The influence of H-bonding on the ‘ambidentate’coordination behaviour of the thiocyanate ion to Cd (II): A combined experimental and theoretical study. Dalton Trans. 2014, 43, 8007–8015. [Google Scholar] [CrossRef]
  16. Politzer, P.; Murray, J.S.; Clark, T. Halogen bonding and other σ-hole interactions: A perspective. Phys. Chem. Chem. Phys. 2013, 15, 11178–11189. [Google Scholar] [CrossRef] [PubMed]
  17. Bauza, A.; Mooibroek, T.J.; Frontera, A. The Bright Future of Unconventional σ/π-Hole Interactions. Chem. Phys. Chem. 2015, 16, 2496–2517. [Google Scholar] [CrossRef]
  18. Riwar, L.J.; Trapp, N.; Root, K.; Zenobi, R.; Diederich, F. Supramolecular Capsules: Strong versus Weak Chalcogen Bonding. Angew. Chem. Int. Ed. 2018, 57, 17259–17264. [Google Scholar] [CrossRef]
  19. Wonner, P.; Dreger, A.; Vogel, L.; Engelage, E.; Huber, S.M. Chalcogen Bonding Catalysis of a Nitro-Michael Reaction. Angew. Chem. Int. Ed. 2019, 58, 16923–16927. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Borissov, A.; Marques, I.; Lim, J.Y.C.; Félix, V.; Smith, M.D.; Beer, P.D. Anion Recognition in Water by Charge-Neutral Halogen and Chalcogen Bonding Foldamer Receptors. J. Am. Chem. Soc. 2019, 141, 4119–4129. [Google Scholar] [CrossRef] [PubMed]
  21. Mallada, B.; Gallardo, A.; Lamanec, M.; de la Torre, B.; Špirko, V.; Hobza, P.; Jelinek, P. Real-space imaging of anisotropic charge of σ-hole by means of Kelvin probe force microscopy. Science 2021, 374, 863–867. [Google Scholar] [CrossRef] [PubMed]
  22. Pascoe, J.; Ling, K.B.; Cockroft, S.L. The Origin of Chalcogen-Bonding Interactions. J. Am. Chem. Soc. 2017, 139, 15160–15167. [Google Scholar] [CrossRef] [Green Version]
  23. Macchione, M.; Goujon, A.; Strakova, K.; Humeniuk, H.V.; Licari, G.; Tajkhorshid, E.; Sakai, N.; Matile, S. A Chalcogen-Bonding Cascade Switch for Planarizable Push-Pull Probes. Angew. Chem. Int. Ed. 2019, 58, 15752–15756. [Google Scholar] [CrossRef]
  24. Politzer, P.; Murray, J.S. An overview of strengths and directionalities of noncovalent interactions: σ-holes and π-holes. Crystals 2019, 9, 165. [Google Scholar] [CrossRef] [Green Version]
  25. Politzer, P.; Murray, J.S.; Clark, T.; Resnati, G. The σ-hole revisited. Phys. Chem. Chem. Phys. 2017, 19, 32166–32178. [Google Scholar] [CrossRef]
  26. Cavallo, G.; Metrangolo, P.; Pilati, T.; Resnati, G.; Terraneo, G. Naming interactions from the electrophilic site. Cryst. Growth Des. 2014, 14, 2697–2702. [Google Scholar] [CrossRef]
  27. Terraneo, G.; Resnati, G. Bonding Matters. Cryst. Growth Des. 2017, 17, 1439–1440. [Google Scholar] [CrossRef] [Green Version]
  28. Scheiner, S. Steric crowding in tetrel bonds. J. Phys. Chem. A 2018, 122, 2550–2562. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Bauzá, A.; Mooibroek, T.J.; Frontera, A. Tetrel Bonding Interactions. Chem. Rec. 2016, 16, 473–487. [Google Scholar] [CrossRef] [PubMed]
  30. Mani, D.; Arunan, E. The X–C⋯Y (X = O/F, Y = O/S/F/Cl/Br/N/P) ‘carbon bond’ and hydrophobic interactions. Phys. Chem. Chem. Phys. 2013, 15, 14377–14383. [Google Scholar] [CrossRef]
  31. Sethio, D.; Oliveira, V.; Kraka, E. Quantitative Assessment of Tetrel Bonding Utilizing Vibrational Spectroscopy. Molecules 2018, 23, 2763. [Google Scholar] [CrossRef] [Green Version]
  32. Heywood, V.L.; Alford, T.P.J.; Roeleveld, J.J.; Deprez, S.J.L.; Verhoofstad, A.; van der Vlugt, J.I.; Domingos, S.R.; Schnell, M.; Davis, A.P.; Mooibroek, T.J. Observations of tetrel bonding between sp3-carbon and THF. Chem. Sci. 2020, 11, 5289–5293. [Google Scholar] [CrossRef]
  33. Zhang, Y.; Wang, W.; Wang, Y.-B. Tetrel bonding on Graphene. Comp. Theor. Chem. 2019, 1147, 8–12. [Google Scholar] [CrossRef]
  34. Taylor, M.S. Anion recognition based on halogen, chalcogen, pnictogen and tetrel bonding. Coord. Chem. Rev. 2020, 413, 213720. [Google Scholar] [CrossRef]
  35. Bauzá, A.; Seth, S.K.; Frontera, A. Tetrel bonding interactions at work: Impact on tin and lead coordination compounds. Coord. Chem. Rev. 2019, 384, 107–125. [Google Scholar] [CrossRef]
  36. Mundlapati, V.R.; Sahoo, D.K.; Bhaumik, S.; Jena, S.; Chandrakar, A.; Biswal, H.S. Noncovalent Carbon-Bonding Interactions in Proteins. Angew. Chem. Int. Ed. 2018, 57, 16496. [Google Scholar] [CrossRef]
  37. Dutta, J.; Sahoo, D.K.; Jena, S.; Tulsiyan, K.D.; Biswal, H.S. Non-covalent interactions with inverted carbon: A carbo-hydrogen bond or a new type of hydrogen bond? Phys. Chem. Chem. Phys. 2020, 22, 8988–8997. [Google Scholar] [CrossRef]
  38. Murray, J.S.; Lane, P.; Clark, T.; Riley, K.E.; Politzer, P. σ-Holes, π-holes and electrostatically-driven interactions. J. Mol. Mod. 2012, 18, 541–548. [Google Scholar] [CrossRef]
  39. Burgi, H.B.; Dunitz, J.D.; Shefter, E. Geometrical reaction coordinates. II. Nucleophilic addition to a carbonyl group. J. Am. Chem. Soc. 1973, 95, 5065–5067. [Google Scholar]
  40. Jönsson, B.; Karlström, G.; Wennerström, H. Ab initio molecular orbital calculations on the water-carbon dioxide system: Molecular complexes. Chem. Phys. Lett. 1975, 30, 58–59. [Google Scholar] [CrossRef]
  41. Peterson, K.I.; Klemperer, W. Structure and internal rotation of H2O−CO2, HDO−CO2 and D2O−CO2 van der Waals complexes. J. Chem. Phys. 1984, 80, 2439–2445. [Google Scholar] [CrossRef]
  42. Peng, Y.P.; Sharpe, S.W.; Shin, S.K.; Wittig, C.; Beaudet, R.A. Infrared spectroscopy of CO2–D(H)Br: Molecular structure and its reliability. J. Chem. Phys. 1992, 97, 5392. [Google Scholar]
  43. Leopold, K.R.; Fraser, G.T.; Klemperer, W. Rotational spectrum and structure of the complex HCNCO2. J. Chem. Phys. 1984, 80, 1039. [Google Scholar] [CrossRef]
  44. Alkorta, I.; Rozas, I.; Elguero, J. Molecular Complexes between Silicon Derivatives and Electron-Rich Groups. J. Phys. Chem. A 2001, 105, 743–749. [Google Scholar] [CrossRef]
  45. Alkorta, I.; Elguero, J.; Fruchier, A.; Macquarrie, D.J.; Virgili, A. Aminopropylsilanes versus silatranes: An experimental and theoretical study. J. Organomet. Chem. 2001, 625, 148–153. [Google Scholar] [CrossRef]
  46. Liu, M.; Li, Q.; Scheiner, S. Comparison of tetrel bonds in neutral and protonated complexes of pyridineTF3 and furanTF3 (T= C, Si, and Ge) with NH3. Phys. Chem. Chem. Phys. 2017, 19, 5550–5559. [Google Scholar] [CrossRef] [Green Version]
  47. Murray, J.S.; Lane, P.; Politzer, P. Expansion of the sigma-hole concept. J. Mol. Model. 2009, 15, 723–729. [Google Scholar] [CrossRef] [PubMed]
  48. Murray, J.S.; Concha, M.C.; Politzer, P. Molecular surface electrostatic potentials as guides to Si-O-N angle contraction: Tunable σ-holes. J. Mol. Model. 2011, 7, 2151–2157. [Google Scholar] [CrossRef] [PubMed]
  49. Patia, S.; Rappoport, Z. (Eds.) The Chemistry of Functional Groups. In The Chemistry of Organic Germanium, Tin and Lead Compounds; Wiley: Hoboken, NJ, USA, 1995; Volume 19. [Google Scholar]
  50. Parr, J. Comprehensive Coordination Chem II; McCleverty, J.A., Meye, T., Jr., Eds.; Elsevier Pergamon: Oxford, UK, 2004; Volume 3, p. 545. [Google Scholar]
  51. Sato, T. Comprehensive Organometallic Chem II—8—Tin; Abel, E.W., Stone, F.G.A., Wilkinson, G., Eds.; Pergamon Press: Oxford, UK, 1995; Volume 11, p. 355. [Google Scholar]
  52. Pinhey, J.T. 11—Lead. In Comprehensive Organometallic Chem II; Abel, E.W., Stone, F.G.A., Wilkinson, G., Eds.; Pergamon Press: Oxford, UK, 1995; Volume 11, p. 461. [Google Scholar]
  53. Greenberg, A.; Wu, G. Structural relationships in silatrane molecules. Struct. Chem. 1990, 1, 79–85. [Google Scholar] [CrossRef]
  54. Hencsei, P. Evaluation of silatrane structures by correlation relationships. Struct. Chem. 1991, 2, 21–26. [Google Scholar] [CrossRef]
  55. Voronkov, M.G.; Barishok, V.P.; Petukhov, L.P.; Rahklin, R.G.; Pestunovich, V.A. 1-Halosilatranes. J. Organomet. Chem. 1988, 358, 39–55. [Google Scholar] [CrossRef]
  56. Lukevics, E.; Dimens, V.; Pokrovska, N.; Zicmane, I.; Popelis, J.; Kemme, A. Addition of nitrile oxides to 2,3-dihydrofurylsilanes. Crystal and molecular structure of tetrahydrofuro-[2,3-d]-isoxazolylsilanes. J. Organomet. Chem. 1999, 586, 200–207. [Google Scholar] [CrossRef]
  57. Franconetti, A.; Frontera, A. “Like–like” tetrel bonding interactions between Sn centres: A combined ab initio and CSD study. Dalton Trans. 2019, 48, 11208–11216. [Google Scholar] [CrossRef] [PubMed]
  58. An, X.; Yang, X.; Li, Q. Tetrel Bonds between Phenyltrifluorosilane and Dimethyl Sulfoxide: Influence of Basis Sets, Substitution and Competition. Molecules 2021, 26, 7231. [Google Scholar] [CrossRef] [PubMed]
  59. Groom, C.R.; Bruno, I.J.; Lightfoot, M.P.; Ward, S.C. The Cambridge Structural Database. Acta Cryst. 2016, 72, 171–179. [Google Scholar] [CrossRef]
  60. Bauzá, A.; Quiñonero, D.; Frontera, A.; Deyà, P.M. Substituent effects in halogen bonding complexes between aromatic donors and acceptors: A comprehensive ab initio study. Phys. Chem. Chem. Phys. 2011, 13, 20371–20379. [Google Scholar] [CrossRef]
  61. Kano, N.; Yamamura, M.; Kawashima, T. 2,2′-Disilylazobenzenes featuring double intramolecular nitrogen⋯silicon coordination: A photoisomerizable fluorophore. Dalton Trans. 2015, 44, 16256–16265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Schröder, A.; Lork, E.; Beckmann, J. A monoclinic polymorph of 2,6-Mes2 C6 H3 SiF3. Main Group Met. Chem. 2014, 37, 153. [Google Scholar] [CrossRef]
  63. Weigend, F.; Häser, M. RI-MP2: First derivatives and global consistency. Theor. Chem. Acc. 1997, 97, 331–340. [Google Scholar] [CrossRef]
  64. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef]
  65. Bauzá, A.; Alkorta, I.; Frontera, A.; Elguero, J. On the Reliability of Pure and Hybrid DFT Methods for the Evaluation of Halogen, Chalcogen, and Pnicogen Bonds Involving Anionic and Neutral Electron Donors. J. Chem. Theory Comput. 2013, 9, 5201–5210. [Google Scholar] [CrossRef] [Green Version]
  66. Ahlrichs, R.; Bar, M.; Haser, M.; Horn, H.; Kolmel, C. Electronic Structure Calculations on Workstation Computers—The Program System turbomole. Chem. Phys. Lett. 1989, 162, 165–169. [Google Scholar] [CrossRef]
  67. Boys, S.B.; Bernardy, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar]
  68. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision B.01; Gaussian, Inc.: Wallingford, UK, 2016. [Google Scholar]
  69. Bader, R.F.W. A quantum theory of molecular structure and its applications. Chem. Rev. 1991, 91, 893–928. [Google Scholar] [CrossRef]
  70. Todd, A.; Keith, T.K. AIMAll, version 13.05.06; Gristmill Software: Overland Park, KS, USA, 2013. [Google Scholar]
  71. Contreras-García, J.; Johnson, E.R.; Keinan, S.; Chaudret, R.; Piquemal, J.-P.; Beratan, D.N.; Yang, W. NCIPLOT: A Program for Plotting Noncovalent Interaction Regions J. Chem. Theory Comput. 2011, 7, 625–632. [Google Scholar] [CrossRef] [PubMed]
  72. Bruno, I.J.; Cole, J.C.; Edgington, P.R.; Kessler, M.; Macrae, C.F.; McCabe, P.; Pearson, J.; Taylor, R. New software for searching the Cambridge Structural Database and visualising crystal structures. Acta Cryst. 2002, 58, 389–397. [Google Scholar] [CrossRef] [PubMed]
Molecules 28 02385 g001 550
Figure 1. Compounds and complexes 1 to 72 studied in this work.
Figure 1. Compounds and complexes 1 to 72 studied in this work.
Molecules 28 02385 g001
Molecules 28 02385 g002 550
Figure 2. Electrostatic potential surfaces (ESP) of some selected compounds used in this study. The energy values given at selected points in the surface are in kcal/mol using an isovalue of 0.001 a.u. The ESP values corresponding to the rest of the compounds are gathered in Table 1.
Figure 2. Electrostatic potential surfaces (ESP) of some selected compounds used in this study. The energy values given at selected points in the surface are in kcal/mol using an isovalue of 0.001 a.u. The ESP values corresponding to the rest of the compounds are gathered in Table 1.
Molecules 28 02385 g002
Molecules 28 02385 g003 550
Figure 3. Electrostatic surface potential values (ESP) vs. the BSSE-corrected binding energies (ΔEBSSE) of several representative meta- (left) and para-substituted (right) complexes.
Figure 3. Electrostatic surface potential values (ESP) vs. the BSSE-corrected binding energies (ΔEBSSE) of several representative meta- (left) and para-substituted (right) complexes.
Molecules 28 02385 g003
Molecules 28 02385 g004 550
Figure 4. Distribution of critical points and bond paths in several tetrel bonding complexes. The intermolecular bond critical point (BCP) is denoted by a red dot. The bond path connecting the bond critical point is also represented. The value of the density at the BCP characterizing the TtB (ρ × 102, a.u.) is also indicated. NCI plot colour range –0.002 a.u. ≤ (signλ2) ρ ≤ –0.002 a.u.
Figure 4. Distribution of critical points and bond paths in several tetrel bonding complexes. The intermolecular bond critical point (BCP) is denoted by a red dot. The bond path connecting the bond critical point is also represented. The value of the density at the BCP characterizing the TtB (ρ × 102, a.u.) is also indicated. NCI plot colour range –0.002 a.u. ≤ (signλ2) ρ ≤ –0.002 a.u.
Molecules 28 02385 g004
Molecules 28 02385 g005 550
Figure 5. Hammett’s plots (ΔEBSSE vs. σ) of meta- and para-substituted complexes involving –CN, –H, –NH2 groups (top left and right, respectively). Hammett’s plots (ΔEBSSE vs. σ) including meta- and para-substituted complexes in the same representation (bottom).
Figure 5. Hammett’s plots (ΔEBSSE vs. σ) of meta- and para-substituted complexes involving –CN, –H, –NH2 groups (top left and right, respectively). Hammett’s plots (ΔEBSSE vs. σ) including meta- and para-substituted complexes in the same representation (bottom).
Molecules 28 02385 g005
Molecules 28 02385 g006 550
Figure 6. Partial views of the X-ray structures corresponding to AHESIR (left) and WIZHUI01 (middle). The tetrel bond distances are given in Å. The CSD codes are also indicated.
Figure 6. Partial views of the X-ray structures corresponding to AHESIR (left) and WIZHUI01 (middle). The tetrel bond distances are given in Å. The CSD codes are also indicated.
Molecules 28 02385 g006
Molecules 28 02385 g007 550
Figure 7. Conceptual models of the tetrel bonding interactions found in the CSD.
Figure 7. Conceptual models of the tetrel bonding interactions found in the CSD.
Molecules 28 02385 g007
Table
Table 1. Values of the electrostatic potential surface (ESP, in kcal/mol) for both tetrel bond donors (TtB donor, measured at the Si atom) and acceptor (TtB acceptor, measured at the N atom) molecules used in this study.
Table 1. Values of the electrostatic potential surface (ESP, in kcal/mol) for both tetrel bond donors (TtB donor, measured at the Si atom) and acceptor (TtB acceptor, measured at the N atom) molecules used in this study.
TtB Donor (para substitution)ESP (kcal/mol)
–CF3+33.3
–CN+35.8
–H+28.2
–CH3+25.9
–NH2+21.9
–OCH3+26.4
TtB donor (meta substitution)ESP (kcal/mol)
–CF3+32.0
–CN+34.5
–H+27.0
–CH3+25.2
–NH2+23.8
–OCH3+26.4
TtB acceptorESP (kcal/mol)
–CF3−23.3
–CN−19.5
–H−32.6
–CH3−38.3
–NH2−40.2
–OCH3−39.2
Table
Table 2. Values of the uncorrected and BSSE-corrected interaction energies (ΔE and ΔEBSSE, in kcal/mol), equilibrium distances (d, in Å) and value of the density (ρ × 100, in a.u.) at the bond critical point (BCP) that characterizes the tetrel bonding interaction in complexes 1 to 36.
Table 2. Values of the uncorrected and BSSE-corrected interaction energies (ΔE and ΔEBSSE, in kcal/mol), equilibrium distances (d, in Å) and value of the density (ρ × 100, in a.u.) at the bond critical point (BCP) that characterizes the tetrel bonding interaction in complexes 1 to 36.
ComplexΔEΔEBSSEdρ × 100
1 (X = –NH2, Y = –NH2)−1.9−1.23.1590.72
2 (X = –NH2, Y = –OCH3)−2.4−1.73.1140.77
3 (X = –NH2, Y = –CH3)−2.4−1.73.1190.77
4 (X =–NH2, Y = –H)−2.7−1.93.0950.79
5 (X = –NH2, Y = –CF3)−3.4−2.63.0300.89
6 (X = –NH2, Y = –CN)−3.7−2.93.0240.93
7 (X = –OCH3, Y = –NH2)−1.9−1.23.1730.69
8 (X = –OCH3, Y = –OCH3)−2.4−1.73.1420.72
9 (X = –OCH3, Y = –CH3)−2.3−1.63.1320.74
10 (X = –OCH3, Y = –H)−2.6−1.93.1080.77
11 (X = –OCH3, Y = –CF3)−3.3−2.53.0460.86
12 (X = –OCH3, Y = –CN)−3.6−2.83.0240.89
13 (X = –CH3, Y = –NH2)−1.9−1.23.1820.69
14 (X = –CH3, Y = –OCH3)−2.3−1.63.1390.74
15 (X = –CH3, Y = –CH3)−2.3−1.63.1430.74
16 (X = –CH3, Y = –H)−2.5−1.83.1210.77
17 (X = –CH3, Y = –CF3)−3.2−2.43.0600.85
18 (X = –CH3, Y = –CN)−3.5−2.73.0390.88
19 (X = –H, Y = –NH2)−1.5−1.03.2590.42
20 (X = –H, Y = –OCH3)−1.9−1.33.2210.64
21 (X = –H, Y = –CH3)−1.9−1.33.2180.64
22 (X = –H, Y = –H)−2.1−1.53.2060.66
23 (X = –H, Y = –CF3)−2.6−2.03.1550.71
24 (X = –H, Y = –CN)−2.8−2.13.1380.74
25 (X = –CF3, Y = –NH2)−1.5−1.03.3100.37
26 (X = –CF3, Y = –OCH3)−1.8−1.23.2680.58
27 (X = –CF3, Y = –CH3)−1.7−1.13.2830.57
28 (X = –CF3, Y = –H)−1.8−1.33.2670.58
29 (X = –CF3, Y = –CF3)−2.2−1.63.2280.62
30 (X = –CF3, Y = –CN)−2.3−1.73.2140.47
31 (X = –CN, Y = –NH2)−1.6−1.03.3240.69
32 (X = –CN, Y = –OCH3)−1.8−1.23.2970.74
33 (X = –CN, Y = –CH3)−1.8−1.23.2990.74
34 (X = –CN, Y = –H)−1.8−1.33.2860.77
35 (X = –CN, Y = –CF3)−2.1−1.53.2500.85
36 (X = –CN, Y = –CN)−2.2−1.63.2380.89
Table
Table 3. Values of the uncorrected and BSSE-corrected interaction energies (ΔE and ΔEBSSE, in kcal/mol), equilibrium distances (d, in Å) and value of the density (ρ × 100, in a.u.) at the bond critical point (BCP) that characterizes the tetrel bonding interaction in complexes 37 to 72.
Table 3. Values of the uncorrected and BSSE-corrected interaction energies (ΔE and ΔEBSSE, in kcal/mol), equilibrium distances (d, in Å) and value of the density (ρ × 100, in a.u.) at the bond critical point (BCP) that characterizes the tetrel bonding interaction in complexes 37 to 72.
ComplexΔEΔEBSSEdρ × 100
37 (X = –NH2, Y = –NH2)−2.4−1.73.1350.74
38 (X = –NH2, Y = –OCH3)−2.4−1.73.1120.77
39 (X = –NH2, Y = –CH3)−2.6−1.93.1310.75
40 (X = –NH2, Y = –H)−2.5−1.83.1620.71
41 (X = –NH2, Y = –CF3)−3.2−2.43.0480.87
42 (X = –NH2, Y = –CN)−3.6−3.63.0140.92
43 (X = –OCH3, Y = –NH2)−2.3−1.63.1600.71
44 (X = –OCH3, Y = –OCH3)−2.7−2.03.1280.75
45 (X = –OCH3, Y = –CH3)−2.6−1.83.1410.73
46 (X = –OCH3, Y = –H)−2.4−1.73.1420.73
47 (X = –OCH3, Y = –CF3)−3.1−2.43.0730.82
48 (X = –OCH3, Y = –CN)−3.5−2.73.0560.84
49 (X = –CH3, Y = –NH2)−2.0−1.33.1760.70
50 (X = –CH3, Y = –OCH3)−2.3−1.63.1640.71
51 (X = –CH3, Y = –CH3)−2.1−1.53.1570.72
52 (X = –CH3, Y = –H)−2.4−1.73.1210.77
53 (X = –CH3, Y = –CF3)−3.0−2.33.0700.84
54 (X = –CH3, Y = –CN)−3.4−2.63.0450.88
55 (X = –H, Y = –NH2)−1.7−1.13.2580.60
56 (X = –H, Y = –OCH3)−1.9−1.33.2320.63
57 (X = –H, Y = –CH3)−2.1−1.53.2380.62
58 (X = –H, Y = –H)−1.9−1.43.2230.64
59 (X = –H, Y = –CF3)−2.4−1.83.1680.70
60 (X = –H, Y = –CN)−2.7−2.73.1380.74
61 (X = –CF3, Y = –NH2)−1.9−1.33.3070.54
62 (X = –CF3, Y = –OCH3)−1.7−1.23.2870.56
63 (X = –CF3, Y = –CH3)−2.0−1.43.2900.56
64 (X = –CF3, Y = –H)−1.8−1.23.2800.57
65 (X = –CF3, Y = –CF3)−2.1−1.53.2290.62
66 (X = –CF3, Y = –CN)−2.3−1.73.2140.64
67 (X = –CN, Y = –NH2)−1.6−1.13.3410.51
68 (X = –CN, Y = –OCH3)−1.7−1.23.3020.55
69 (X = –CN, Y = –CH3)−1.7−1.13.3180.53
70 (X = –CN, Y = –H)−1.8−1.23.3080.54
71 (X = –CN, Y = –CF3)−2.0−1.43.2770.57
72 (X = –CN, Y = –CN)−2.2−1.63.2450.60
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Burguera, S.; Frontera, A.; Bauzá, A. Substituent Effects in Tetrel Bonds Involving Aromatic Silane Derivatives: An ab initio Study. Molecules 2023, 28, 2385. https://doi.org/10.3390/molecules28052385

AMA Style

Burguera S, Frontera A, Bauzá A. Substituent Effects in Tetrel Bonds Involving Aromatic Silane Derivatives: An ab initio Study. Molecules. 2023; 28(5):2385. https://doi.org/10.3390/molecules28052385

Chicago/Turabian Style

Burguera, Sergi, Antonio Frontera, and Antonio Bauzá. 2023. "Substituent Effects in Tetrel Bonds Involving Aromatic Silane Derivatives: An ab initio Study" Molecules 28, no. 5: 2385. https://doi.org/10.3390/molecules28052385

APA Style

Burguera, S., Frontera, A., & Bauzá, A. (2023). Substituent Effects in Tetrel Bonds Involving Aromatic Silane Derivatives: An ab initio Study. Molecules, 28(5), 2385. https://doi.org/10.3390/molecules28052385

Article Metrics

Back to TopTop