Next Article in Journal
2-Deoxyglucose-Modified Folate Derivative: Self-Assembling Nanoparticle Able to Load Cisplatin
Next Article in Special Issue
Tetrel Interactions from an Interacting Quantum Atoms Perspective
Previous Article in Journal
Characterization and Antibacterial Activity of 7S and 11S Globulins Isolated from Cowpea Seed Protein
Previous Article in Special Issue
Calculation of VS,max and Its Use as a Descriptor for the Theoretical Calculation of pKa Values for Carboxylic Acids
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 24
Issue 6
10.3390/molecules24061083
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

Identification of the Tetrel Bonds between Halide Anions and Carbon Atom of Methyl Groups Using Electronic Criterion

by
Ekaterina Bartashevich
1,*,
Yury Matveychuk
1 and
Vladimir Tsirelson
1,2
1
Research Laboratory of Multiscale Modelling of Multicomponent Functional Materials, REC Nanotechnology, South Ural State University, 454080 Chelyabinsk, Russia
2
Quantum Chemistry Department, D.I. Mendeleev University of Chemical Technology, 125047 Moscow, Russia
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(6), 1083; https://doi.org/10.3390/molecules24061083
Submission received: 26 February 2019 / Revised: 15 March 2019 / Accepted: 16 March 2019 / Published: 19 March 2019
(This article belongs to the Special Issue Tetrel Bonds)
Browse Figures
Versions Notes

Abstract

:
The consideration of the disposition of minima of electron density and electrostatic potential along the line between non-covalently bound atoms in systems with Hal−···CH3–Y (Hal− = Cl, Br; Y = N, O) fragments allowed to prove that the carbon atom in methyl group serves as an electrophilic site provider. These interactions between halide anion and carbon in methyl group can be categorized as the typical tetrel bonds. Statistics of geometrical parameters for such tetrel bonds in CSD is analyzed. It is established that the binding energy in molecular complexes with tetrel bonds correlate with the potential acting on an electron in molecule (PAEM). The PAEM barriers for tetrel bonds show a similar behavior for both sets of complexes with Br− and Cl− electron donors.

Graphical Abstract

1. Introduction

The problem of categorizing non-covalent interactions in molecular crystals and complexes is now a focus of attention [1,2]. Nowadays, the systematization of the halogen, chalcogen, pnictogen, and tetrel bonds already exists [3]; however, in most cases, only the simplistic geometrical approach underlies the analysis of such types of interactions. In this context, the types of non-covalent interactions are traditionally discussed in terms of interatomic distances and angles, which specify the mutual orientation of pivotal chemical bonds [4,5]. However, due to the pronounced and specific electrostatic component of such non-covalent bonds [6,7], more careful analysis of the electronic features of the halogen, chalcogen, pnictogen, and tetrel bonds is required. Such analysis needs to focus on the features of valence electron shells and related anisotropy of the electrostatic potential of interacting atoms.
The estimation of the binding energy for molecules in the Y4T···Hal− complexes, where a tetrel atom T = C, Si, Ge, Sn, as well as the description of electron density characteristics for tetrel bonds, were presented in References [8,9,10,11,12,13]. The carbon atom in the СH3-group is fairly often noted as the owner of σ-hole, and the fact that the oxygen atom can act as an electron-rich center in the CH3···O tetrel bonding has been confirmed in studies [14,15]. Note that in these early works such non-covalent interaction has been referred as a “carbon bond”. Pal et al. [16] described the CH3···N tetrel bonding in a Co(II) coordination polymeric system using the analysis of calculated electron density. The typical tetrel bonds formed by СH3-group in crystals have been also been observed by high-precision X-ray diffraction method using the analysis of the experimental electron density [17,18].
The rows of binding energy were calculated for the series of complexes in which the compounds NH3 [19,20,21,22], PH3 and AsH3 [23], benzene and unsaturated hydrocarbons [24], HCN and pyrazine [25] acted as the tetrel bond acceptors. According to these results, the tetrel bond strength depends on the tetrel atom, the Lewis base that acts as a tetrel bond acceptor, and the fragments covalently bound with this atom. Therefore, in the complexes with ammonia at transition from CF4 to SnF4, the binding energy increases from −0.82 to −25.53 kcal/mol (MP2/aug-cc-pVDZ) [21], that is, the tetrels of higher periods form stronger interactions in complexes. Decrease of the Lewis base strength in series from X = NH3 to AsH3 reduces the binding energy from −3.23 to −1.28 kcal/mol (MP2/aug-cc-pVDZ) in complexes X···CH3Cl [23]. Therefore, the values of binding energy estimated in neutral complexes with tetrel bonds formed by the carbon in a methyl group are commonly very small. The binding energy grows more sharply if the hydrogen atoms are replaced by halogens, which leads to enhancement of electron acceptor properties of the tetrel atom providing the σ-hole for bonding. This fact is illustrated by complexes H4Sn···NH3 (−2.44 kcal/mol) and F4Sn···NH3 (−25.53 kcal/mol) with the tetrel bond lengths of 3.17 and 2.28 Å (MP2/aug-cc-pVDZ), respectively [21]. If the anion acts as the tetrel bond acceptor, the estimated binding energy can be greater by an order of magnitude as compared to the neutral complexes [9,10,11]. The highest value of binding energy, −93.58 kcal/mol, has been recorded for the SnF4···F− complex (MP2/aug-ccpVTZ) [11]; replacing the F− ion with Cl− and Br− reduces this energy by 2–3 times.
Therefore, a variety of functional groups can act as the tetrel bond donors in molecular complexes and a methyl group as an electrophilic site deliverer is quite often in the focus of attention for the study of tetrel bond properties. The tetrel bonds formed by the methyl carbon atom occur in crystalline systems as well. There are at least two studies [17,26] based on high-precision X-ray diffraction data, which established the participation of the carbon atom of a methyl group in non-covalent interactions. These facts motivate us to focus on the search for evidential electronic criterion for recognizing the type of non-covalent bonding and for systematization of the electrophilic site features for the carbon atom of a methyl group in molecular complexes and crystals.
Quantum Theory of Atoms in Molecules (QTAIM) [27] suggests zero-flux conditions for electron-density gradient [28] and electrostatic-potential gradient [29] to determine the boundaries of chemically bonded atoms and electrically neutral atomic fragments, respectively. For non-covalent interactions with a significant electrostatic component, the boundaries of the electron density basins (ρ-basins) do not coincide with the boundaries of the electrostatic potential (ESP) basins (φ-basins). It means that electrons formally belonging to the electron donor atom can be electrostatically attracted to the nucleus of an electron acceptor atom along a specific direction. The superposition of the ρ- and φ-basins has been discussed in the literature [30,31,32], and the features of zero-flux surface in electrostatic potential in solids have been studied using the experimental charge density [29,33].
Previously we have proposed the following electronic criterion for recognizing the atom that prescribes the name of non-covalent bonding [34,35]. The minimum of electrostatic potential along the interaction line is located at the side of the atom that donates electrons; the minimum of electron density is closer to the atom that delivers its electrophilic site for bonding. More explicitly, the latter atom prescribes the name of bonding. This observation opens up the broad possibilities for identifying the role of atoms involved in non-covalent interactions. For example, atoms of the Group 14 of the Periodic Table are able to deliver their electrophilic sites for non-covalent interactions with electron donors for the tetrel bond formation, and the suggested electronic criterion can specify the electron acceptor role of these atoms. Along the Hal−···C interatomic line, the 1D minimum of ESP should be located closer to the electron donor atom, while the minimum of the electron density will be found closer to the atom that is an electron acceptor. Such disposition of minima indicates that the fraction of electron density from the atomic ρ-basin of an electron donor is electrostatically attracted to the nucleus of an electron acceptor atom. In this case, the disposition of the 1D minima of both functions can unambiguously indicate that only one of the pair of atoms provides its electrophilic site. Namely, the minimum of electron density on the interatomic line is always closer to the atom that has delivered its own electrophilic site for bonding.
We have recently proposed the potential acting on an electron in molecule (PAEM [36]) [37] as a function that not only characterizes the properties of non-covalent bonding with a significant electrostatic component, but also allows us to observe the quantitative relationship with the interaction energy in complexes. Unlike ESP, the PAEM contains both Coulomb and exchange components. The first of them has a classic nature and the second one is the two-electron contribution of the quantum exchange-correlation potential. PAEM was examined [38] for the halogen and chalcogen bonding characterization and its usefulness was confirmed.
The aim of the present study is to demonstrate the efficiency and productivity of the above-mentioned electronic criterion and the PAEM for analysis of tetrel bonds between the carbon atom of methyl groups and halide anions, which sometimes occur in molecular crystals. We also try to understand to what extent the characteristics of the gas-phase complexes are suitable for describing Hal−···CH3Y (Hal− = Cl, Br; Y = N, O) tetrel bonds in crystals.

2. Results and Discussion

2.1. Population of Hal−···CH3Y Tetrel Bonds in Crystals

The search of short contacts between a halide anion and the methyl carbon, Hal−···CH3–Y, in Cambridge Structural Database (CSD) v.5.39 (The Cambridge Crystallographic Data Centre, Cambridge, UK) [39] was performed with the following restrictions: organic derivatives without disordered, polymeric, powder, organometallic and repeating structures have been considered. The main condition was to choose interactions in which the halide anion, Hal− = Cl, Br, I, was placed on the extension of the C–Y covalent bond, where Y = C, N, O. We set that condition using the angle θ (Hal−–C–Y), the value of which was in the range from 160° to 180°. At the same time, we selected the structures with interatomic distances, d(Hal, C), falling into the range (rvdw(C) + rvdw(Hal) ± 0.2 Å), where rvdw is the Bondi atomic radius [40]. The total number of selected structures that satisfied those conditions was 164. The analysis of the obtained sample has shown that the Y atom covalently bound with the CH3 group is nitrogen in most cases. We have found 43 cases of Cl−···CH3–N interactions, 36 cases of Br−···CH3–N interactions and 53 cases of I−···CH3–N interactions. Oxygen and carbon are involved in such covalent bonds much less frequently and approximately equally. In these cases, the methyl group forms more interactions with the Cl− anion than with the Br− and I− anions taken together. It can be concluded that the polarity of the covalent bond CH3–Y affects the probability and strength of the tetrel bond formation. It should be noted that the distances d(Hal, C) for the cases of Hal−···CH3–C interactions everywhere exceed the sum of van der Waals radii (Figure 1). If the CH3-group is bound with oxygen (Hal−···CH3–O), the distances d(Hal, C) can be less than the sum of van der Waals radii, though all of these distances are more than this sum for I− cases.
Therefore, we conclude that the studied type of interactions, Hal−···CH3–Y, are not widely spread within crystals listed in CSD, but they are not exceptional.

2.2. Evidence of Electrophilic Sites for the CH3-Groups Bound in Tetrel Bonds

Let us now look at the same examples of halide crystal structures containing tetrel bonds (Figure 2). All the results considered in this chapter were obtained for crystal structures by the calculations with the periodic boundary conditions. In the crystalline N,N,N′,N′-tetramethylchloroformamide chloride, LONGEB [41], the Cl− anion forms non-covalent interactions with the Cl, H, C atoms, which are characterized by interatomic distances smaller than the sums of van der Waals radii. The Cl−···Cl–C non-covalent interaction with a distance of 3.122 Å refers to a typical charge-assisted halogen bond; the next five hydrogen bonds, Сl−···H–C, are characterized by interatomic distances ranging from 2.924 to 2.664 Å. Finally, the Cl−···CH3 interaction of 3.425 Å can be called a tetrel bond. In the crystalline dimethylmethyleneammonium chloride VAPREJ [42], the chloride anion forms eight Сl−···H–C interactions, which are shorter than the sum of van der Waals radii and two Cl−···CH3 tetrel bonds. In the dimethylmethylenimine bromide crystal, LILLOH [43], in addition to multiple Br−···H–C interactions, there are two tetrel bonds: Br−···CH3 (3.533 Å) and Br−···CH2 (3.503 Å). Quantum-topological analysis of the electron density in all considered crystals have confirmed the presence of the Hal−···C bond path and bcp of electron density (Table 1). Our series of tetrel bonds in the considered crystals does not vary much, and we have observed the small changes in electron density at the bond critical points, ρ(rbcp), which are in the range 0.0042–0.0070 a.u. for Cl−···CH3–Y and 0.0060–0.0068 a.u. for Br−···CH3–Y.
It is possible to demonstrate the electrophilic site on a carbon atom using the electrostatic potential (ESP) mapped on the isosurface of electron density or the distribution of Electron Localization Function (ELF) [44] for CH3-group, which participate in a tetrel bond. For example, relatively higher positive values of ESP on the isosurface of electron density (0.003 a.u.) we can see in the region of the σ-hole, which belongs to the C atomic basin in trimethylammonium cation (Figure 3a,b).
Moving strictly along the line linking the Сl(1)− and C(4) atoms of CH3-group (Figure 4а), we reach the extension of the covalent bond formed by CH3-group. Along this line the ELF is less than 0.5 near the carbon atom, showing the region of the reduced probability of electron pairing. It can be considered as the manifestation of the carbon atom σ-hole. The values of ELF (rbcp) at the bond critical points do not exceed 0.05. This excludes the hypothesis about significant covalent character of the Hal−···CH3Y tetrel bonds. At about 0.8 Å from the chloride anion nucleus the maximum values of ELF are distributed around the circumference. Figure 4b, depicting the tetrel bonds formed by the bromide anion in LILLOH crystal, shows a similar ELF distribution. Near the carbon atom and along the Br(1)− ···C(2) line the ELF does not attain high values, but it increases sharply, affecting the basins of hydrogen atoms, if we slightly deviate from this line.
Now let us consider and evaluate how the electronic criterion works for the cases of non-covalent interactions formed by the carbon atoms of methyl groups in halide crystals. In Figure 5a it can be seen that in the LONGEB crystal, the Cl− anion forms two non-covalent interactions at least, as follows from the presence of corresponding bcp. In both cases, Cl(2)···Cl(1)−···C(4), the one-dimensional ESP minimum is closer to electron donating anion, Cl(1)− in a crystal. The electron density minima along the Cl(2)···Cl(1)− and Cl(1)−···C(4) lines are located on the side of the Cl(2) and C(4) atoms. They indicate the electrophilic site providers and dictate the name of the non-covalent bonding.
According to the proposed electronic criterion, the first interaction, Cl(2)···Cl(1)−, can be categorized as a charge-assisted halogen bond, and the second one, Cl(1)−···C (4), is a tetrel bond enhanced by charges. In Figure 5b, the minimum of electron density along the Br(1)−···C(2) line is located on the side of the C(4) atom, while the minimum of ESP is closer to Br(1)−. Such disposition of minima shows that the carbon atom accepts electrons along the Br(1)−···C(2) line and that interaction can be called a tetrel bond.

2.3. Binding Energy in Molecular Complexes with the Hal−···CH3 Tetrel Bonds

The determination of the equilibrium geometry for ion pairs “halide anion–cation” extracted from the crystal environment is not a straightforward procedure. In general, the retention of the halide anion position strictly on the extension of a covalent bond of CH3-group is rather difficult in the gas phase state. This task requires us to dwell on the level of gas phase calculations, different from those used for crystal structures. Nevertheless, this step allows us to obtain the stationary state for the maximal number of complexes, for which the tetrel bonds in crystalline state have attracted our attention. Some relative estimations and the features of electronic properties can be quite useful for understanding the nature of charge-assisted tetrel bonds.
The binding energy, Eb, between the halide anion and cation in the considered complexes varies from −52.28 to −82.67 kcal/mol (Table 2). These values do not fall out of the range that is determined in similar studies [8,9,10,11]. The BSSE correction, ΔEBSSE, is negligible and influences the energy values of the third decimal place of kcal/mol units.
The considered tetrel bonds exhibit significantly shorter lengths in the models of complexes extracted from the crystalline environment. On average, the observed Hal−···C bond lengths in such complexes differ by ~17% from those in crystal structures. As a result, the different approaches for complexes and crystals calculations lead to the values of ρ(rbcp) that are almost twice higher in crystals. Obviously, the direct transfer of tetrel bond properties in isolated complexes to the crystals, neglecting the rest of interactions between a halide anion and crystalline environment, is not entirely correct. Comparing the properties of the CH3–Y (Y=N, O) covalent bonds in complexes and isolated cations, we see that the participation of CH3-group in tetrel bond with Hal− weakens the CH3–Y covalent bond. Therefore, the tetrel bonding elongates the covalent bond of a methyl group by 0.01–0.04 Å, and the values of ρ(rbcp) for the CH3–Y bonds decrease by ~8%.
It is useful to understand how the electronic properties of Hal−···CH3Y tetrel bonds in complexes are related to the strength of complexes. In our opinion, the tetrel bonds belong to electrostatically driven interactions. In addition, the binding energy between two oppositely charged ions is much higher in comparison with neutral molecules. For this reason, the electrostatic properties of tetrel bonds have been analyzed first of all.
We found that the properties of both ESP and PAEM for the Hal−···CH3Y tetrel bonds are linearly correlated with the binding energy, Eb, in complexes, as shown in Figure 6. The correlation coefficient for the minima of the electrostatic potential, ESPmin, on the line between Hal− and C atoms is 0.917. For the maximum of PAEM along this line, PAEMmax, or PAEM barrier, the correlation coefficient is 0.985. It is important to note that in the relationship “PAEMmax vs Eb“, PAEMmax for the tetrel bonds formed by Br− and Cl− fits strictly on the common line. This is a rare case among the established relationships between local properties of non-covalent bonds and the binding energy for bound fragments. For example, the electronic potential and kinetic energy densities at bcp do not allow constructing a good common relationship for Br− and Cl− rows (Figure S1). This finding has been discussed by us earlier for the halogen bonds formed by different atoms or fragments that play the role of halogen acceptors [45]. This fact has recently been illustrated in detail in Reference [46], where the large series of non-covalent interactions with different halide anions have been studied.
Note that the extreme values of ESPmin and PAEMmax slightly differ from their local values at the bond critical points of tetrel bonds (Table S4). Nevertheless, the PAEM(rbcp) values correlate with the binding energy better than ESP(rbcp) (Figure S2). This is probably due to the fact that in our series the maximum of PAEM is closer to the tetrel bond critical point than the minimum of ESP (Figure 7). The relative location of ESPmin and PAEMmax in the common projection is demonstrated by examples of the weakest Cl−···CH3 (ZENJAD) and the strongest Br−···CH3 (FADXIR) tetrel bonds in our set. Though the gap between the PAEMmax and the minimum of electron density is larger for the Cl−···CH3 tetrel bond, and ESP has a lower negative minimum, it can be seen that the PAEM barrier is higher in absolute value. It means that the Cl−···CH3 tetrel bond is weaker, and this is confirmed by the linear correlation between PAEMmax and Eb. Moreover, the relative positions of ESPmin and PAEMmax along the tetrel bond line allow us to distinguish, which atom is the acceptor of electrons. As it has been noted earlier [45], and as can be seen from the above, PAEMmax position is located closer to the electrophilic site, while the position of ESPmin is closer to the electron donor.

3. Materials and Methods

The structure optimization of molecular complexes consisting of organic cations and halide anions was carried out at M06-2X/aug-cc-pVDZ level [47,48,49] in GAMESS (v. 2017 R2, Mark Gordon’s Quantum Theory Group, Ames Laboratory, Iowa State University, Ames, IA, USA [50]) with gradient convergence that equaled 0.00001. The optimized structures were tested for the absence of imaginary frequencies. The binding energy between cations and halide anions in electrically neutral complexes was estimated as Eb = Ecom − (EHal + Ecat) − ΔEBSSE, where Ecom, Ecat, EHal were the total energies of the optimized complex, relaxed isolated organic cation and halide anion. BSSE correction, ΔEBSSE, was carried out taking into account the phantom orbitals in complexes calculated for compounds without energy relaxation, see Table S1 in Supporting Information.
All calculations with periodic boundary conditions were performed using CRYSTAL14 (v. 1.0.4, CRYSTAL Theoretical Chemistry Group, Chemistry Department, University of Turin, Turin, Italy [51]) at the B3LYP/6-31G** level for C, N, O, H atoms and DZVP basis set for halogen atoms [52,53] with Grimme dispersion correction D2 [54]. The structure relaxation was carried out with the atomic coordinate optimization only, with the fixed unit cell parameters for the purpose of maximum conformity to experimental data. The following convergence parameters have been used for all calculation: TOLDEG (root-mean-square on gradient) is less than 0.0001 a.u., TOLDEX (root-mean-square on estimated displacements) is less than 0.0003 a.u., TOLDEE (energy change between optimization steps threshold) is less than 10−10 a.u., TOLINTEG (truncation criteria for bielectronic integrals: overlap threshold for Coulomb integrals; penetration threshold for Coulomb integrals; overlap threshold for HF exchange integrals; pseudo-overlap for g and n HF exchange series) are 10, 10, 10, 10 and 16, respectively. The number of k-points in the Pack–Monkhorst net (in the irreducible part of Brillouin zone) was 125 or 170 depending on crystals; the number of k-points in the Gilat net was 729 or 1170, that corresponded to the set SHRINK 8 16 values. All calculations for isolated cations were performed using CRYSTAL17 (v. 1.0.2, CRYSTAL Theoretical Chemistry Group, Chemistry Department, University of Turin, Turin, Italy [55]) at the B3LYP/6-31G** level with the Grimme dispersion correction D2 and DOPING option to account for the cation positive charge.
The QTAIM analyses of electron density and electrostatic potential were carried out in TOPOND [56] in crystals and in AIMAll software package [57] for the complexes. PAEM and ESP distributions were computed using Multiwfn [58] program (Beijing Kein Research Center for Natural Sciences, Beijing, China).
The reported calculations were performed using the supercomputer resources of the South Ural State University [59].

4. Conclusions

In this computational study, the charge-assisted tetrel bonds in the crystals formed between halide anions and the methyl groups of organic cations such as Hal−···CH3Y (Hal− = Cl, Br; Y = N, O) have been considered. The bond paths between the Hal− and C atoms confirm the existence of these uncommon bonds in both the crystal structures and gas phase complexes. To define the type of Hal−···CH3Y bonding in crystals more precisely, we have suggested using the order of one-dimensional minima of electron density and electrostatic potential along the interatomic lines between the carbon atom of СH3-group and the halide anion. This allowed us to apply a simple criterion which reveals that the carbon atom provides its electrophilic site for a typical tetrel bond formation.
The strong correlation between the binding energy in complexes and the extreme values of potential acting on an electron in a molecule calculated along the lines between the Hal− and C atoms has been obtained. Therefore, PAEM extends and enforces the electronic criterion for revealing electrophilic sites and sheds some light on the nature of tetrel bonds. We may speculate that its application will be useful for the other electrostatically driven non-covalent interactions as well.

Supplementary Materials

The following are available online, Figure S1: Binding energy (kcal/mol) in complexes vs the potential (a) and kinetic (b) energy density (a.u.) at the bond critical point of tetrel bonds, Figure S2: Binding energy (kcal/mol) in complexes vs the electrostatic potential (a.u.) (a) and (b), potential acting on an electron in molecule (a.u.) at the bond critical point of tetrel bonds, Figure S3: ESP in the trimethylammonium chloride on the isosurface of electron density of 0.02 a.u, Table S1: The energy characteristics of Hal−···CH3–YR (Hal− = Cl, Br) complexes taken from crystal structures with listed refcodes, Table S2. Experimental and calculated tetrel and C–Y bond lengths D (Å), angles Hal−···C–Y and electron density (a.u.) at bond critical points for considered crystal and cation structures calculated in CRYSTAL code. Table S3, Bond lengths D(Å), the characteristics of electron density, potential and kinetic energy densities (a.u.), electrostatic potential (a.u.), potential acting on an electron in molecule PAEM at bcp (a.u.) for Hal−···CH3 and Y–C bonds in complexes and cations calculated in GAMESS code.

Author Contributions

Conceptualization, V.T. and E.B.; investigation, Y.M.; formal analysis and project administration E.B.

Funding

This research was funded by the Russian Foundation for Basic Research, grant No. 17-03-00406, Ministry of Science and Higher Education of the Russian Federation (grant No. 4.1157.2017/4.6) and by the Government of the Russian Federation, Act 211, contract No. 02.A03.21.0011.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Legon, A.C. Tetrel, Pnictogen and Chalcogen Bonds Identified in the Gas Phase before they had Names: A Systematic Look at Non-covalent Interactions. Phys. Chem. Chem. Phys. 2017, 19, 14884–14896. [Google Scholar] [CrossRef] [PubMed]
  2. 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]
  3. Terraneo, G.; Resnati, G. Bonding Matters. Cryst. Growth Des. 2017, 17, 1439–1440. [Google Scholar] [CrossRef]
  4. Desiraju, G.R.; Parthasarathy, R. The Nature of Halogen···Halogen Interactions: Are Short Halogen Contacts Due to Specific Attractive Forces or Due to Close Packing of Nonspherical Atoms? J. Am. Chem. Soc. 1989, 111, 8725–8726. [Google Scholar] [CrossRef]
  5. 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]
  6. Politzer, P.; Murray, J.S.; Clark, T. Halogen Bonding: An Electrostatically-Driven Highly Directional Noncovalent Interaction. Phys. Chem. Chem. Phys. 2010, 12, 7748–7757. [Google Scholar] [CrossRef]
  7. Politzer, P.; Murray, J.S.; Clark, T.; Resnati, G. The σ-hole Revisited. Phys. Chem. Chem. Phys. 2017, 19, 32166–32178. [Google Scholar] [CrossRef]
  8. Del Bene, J.E.; Alkorta, I.; Elguero, J. Anionic Complexes of F− and Cl− with Substituted Methanes: Hydrogen, Halogen, and Tetrel Bonds. Chem. Phys. Lett. 2016, 655–656, 115–119. [Google Scholar] [CrossRef]
  9. Scheiner, S. Tetrel Bonding as a Vehicle for Strong and Selective Anion Binding. Molecules 2018, 23, 1147. [Google Scholar] [CrossRef]
  10. Esrafili, M.D.; Asadollahi, S.; Mousavian, P. Anionic Tetrel Bonds: An ab Initio Study. Chem. Phys. Lett. 2018, 691, 394–400. [Google Scholar] [CrossRef]
  11. Esrafili, M.D.; Mousavian, P. Strong Tetrel Bonds: Theoretical Aspects and Experimental Evidence. Molecules 2018, 23, 2642. [Google Scholar] [CrossRef]
  12. Grabowski, S.J. Tetrel bond–σ-hole bond as a preliminary stage of the SN2 reaction. Phys. Chem. Chem. Phys. 2014, 16, 1824–1834. [Google Scholar] [CrossRef]
  13. 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]
  14. Bauza, A.; Frontera, A. RCH3···O Interactions in Biological Systems: Are They Trifurcated H-Bonds or Noncovalent Carbon Bonds? Crystals 2016, 6, 26. [Google Scholar] [CrossRef]
  15. Garcia-LLinas, X.; Bauza, A.; Seth, S.K.; Frontera, A. Importance of R−CF3···O Tetrel Bonding Interactions in Biological Systems. J. Phys. Chem. A 2017, 121, 5371–5376. [Google Scholar] [CrossRef] [PubMed]
  16. Pal, P.; Konar, S.; Lama, P.; Das, K.; Bauza, A.; Frontera, A.; Mukhopadhyay, S. On the Importance of Noncovalent Carbon Bonding Interactions in the Stabilization of a 1D Co(II) Polymeric Chain as Precursor of a Novel 2D Coordination Polymer. J. Phys. Chem. B 2016, 120, 6803–6811. [Google Scholar] [CrossRef]
  17. Tsirelson, V.G.; Stash, A.I.; Potemkin, V.A.; Rykounov, A.A.; Shutalev, A.D.; Zhurova, E.A.; Zhurov, V.V.; Pinkerton, A.A.; Gurskaya, G.V.; Zavodnik, V.E. Molecular and Crystal Properties of Ethyl-4,6-dimethyl-2-thioxo-1,2,3,4-tetrahydropyrimidine-5-carboxylate from Experimental and Theoretical Electron Densities. Acta Cryst. B 2006, B62, 676–688. [Google Scholar] [CrossRef]
  18. Thomas, S.P.; Pavan, M.S.; Row, T.N.G. Experimental Evidence for ‘Carbon Bonding’ in the Solid State from Charge Density Analysis. Chem. Commun. 2014, 50, 49–51. [Google Scholar] [CrossRef]
  19. Scheiner, S. Steric Crowding in Tetrel Bonds. J. Phys. Chem. A 2018, 122, 2550–2562. [Google Scholar] [CrossRef]
  20. Zierkiewicz, W.; Michalczyk1, M.; Scheiner, S. Comparison between Tetrel Bonded Complexes Stabilized by σ and π Hole Interactions. Molecules 2018, 23, 1416. [Google Scholar] [CrossRef]
  21. Scheiner, S. Systematic Elucidation of Factors That Influence the Strength of Tetrel Bonds. J. Phys. Chem. A 2017, 121, 5561–5568. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. 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] [PubMed]
  23. Laconsay, C.J.; Galbraith, J.M. A Valence Bond Theory Treatment of Tetrel Bonding Interactions. Comp. Theor. Chem. 2017, 1116, 202–206. [Google Scholar] [CrossRef]
  24. Grabowski, S.J. Tetrel Bonds with π-Electrons Acting as Lewis Bases—Theoretical Results and Experimental Evidences. Molecules 2018, 23, 1183. [Google Scholar] [CrossRef] [PubMed]
  25. Zierkiewicz, W.; Michalczyk1, M.; Scheiner, S. Implications of Monomer Deformation for Tetrel and Pnicogen Bonds. Phys. Chem. Chem. Phys. 2018, 20, 8832–8841. [Google Scholar] [CrossRef] [PubMed]
  26. Stash, A.I.; Chen, Y.S.; Kovalchukova, O.V.; Tsirelson, V.G. Electron Density, Electrostatic Potential, and Spatial Organization of Ammonium Hydrooxalate Oxalic Acid Dihydrate Heteromolecular Crystal from Data of Diffraction Experiment at 15 K Using Synchrotron Radiation and Theoretical Calculations. Russ. Chem. Bull. 2013, 62, 1752–1763. [Google Scholar] [CrossRef]
  27. Bader, R.F.W. Atoms in Molecules: A Quantum Theory; Clarendon Press: Oxford, UK, 1990; pp. 1–438. [Google Scholar]
  28. Bader, R.F.W.; Carroll, M.T.; Cheeseman, J.R.; Chang, C. Properties of Atoms in Molecules: Atomic Volumes. J. Am. Chem. Soc. 1987, 109, 7968–7979. [Google Scholar] [CrossRef]
  29. Tsirelson, V.G.; Avilov, A.S.; Lepeshov, G.G.; Kulygin, A.K.; Stahn, J.; Pietsch, U.; Spence, J.C.H. Quantitative Analysis of the Electrostatic Potential in Rock-Salt Crystals Using Accurate Electron Diffraction Data. J. Phys. Chem. 2001, B105, 5068–5074. [Google Scholar] [CrossRef]
  30. Tsirelson, V.G.; Shishkina, A.V.; Stash, A.I.; Parsons, S. The Experimental and Theoretical QTAIMC Study of the Atomic and Molecular Interactions in Dinitrogen Tetroxide. Acta Crystallogr. Sect. B Struct. Sci. Cryst. Eng. Mater. 2009, B65, 647–658. [Google Scholar] [CrossRef] [PubMed]
  31. Mata, I.; Molins, E.; Alkorta, I.; Espinosa, E. Topological Properties of the Electrostatic Potential in Weak and Moderate N···H Hydrogen Bonds. J. Phys. Chem. 2007, A111, 6425–6433. [Google Scholar] [CrossRef]
  32. Bartashevich, E.V.; Yushina, I.D.; Kropotina, K.K.; Muhitdinova, S.E.; Tsirelson, V.G. Testing the Tools for Revealing and Characterizing the Iodine-Iodine Halogen Bond in Crystals. Acta Crystallogr. Sect. B Struct. Sci. Cryst. Eng. Mater. 2017, B73, 217–226. [Google Scholar] [CrossRef] [PubMed]
  33. Pathak, R.K.; Gadre, S.R. Maximal and Minimal Characteristics of Molecular Electrostatic Potentials. J. Chem. Phys. 1990, 93, 1770–1773. [Google Scholar] [CrossRef]
  34. Bartashevich, E.V.; Yushina, I.D.; Stash, A.I.; Tsirelson, V.G. Halogen Bonding and Other Iodine Interactions in Crystals of Dihydrothiazolo(oxazino)quinolinium Oligoiodides from the Electron-Density Viewpoint. Cryst. Growth Des. 2014, 14, 5674–5684. [Google Scholar] [CrossRef]
  35. Bartashevich, E.V.; Yushina, I.D.; Muhitdinova, S.E.; Tsirelson, V.G. Electronic Criterion for Categorizing the Chalcogen and Halogen Bonds: Sulfur—Iodine Interactions in Crystals. Acta Crystallogr. Sect. B Struct. Sci. Cryst. Eng. Mater. 2019, in press. [Google Scholar] [CrossRef]
  36. Zhao, D.-X.; Gong, L.-D.; Yang, Z.-Z. The Relations of Bond Length and Force Constant with the Potential Acting on an Electron in a Molecule. J. Phys. Chem. A 2005, 109, 10121–10128. [Google Scholar] [CrossRef] [PubMed]
  37. Bartashevich, E.V.; Tsirelson, V.G. A Comparative View on the Potential Acting on an Electron in a Molecule and the Electrostatic Potential through the Typical Halogen Bonds. J. Comput. Chem. 2018, 39, 573–580. [Google Scholar] [CrossRef] [PubMed]
  38. Bartashevich, E.V.; Mukhitdinova, S.E.; Tsirelson, V.G. Characterizing the Halogen and Chalcogen Bonds in Crystals: PAEM vs ESP. In Book of Abstracts of International Union of Crystallography (IUCr)’s Sagamore XIX Conference on Quantum Crystallography, Halifax, NS, Canada; Mount Saint Vincent University’s Printshop: Halifax, Canada, 2018; pp. 93–95. [Google Scholar]
  39. Groom, C.R.; Bruno, I.J.; Lightfoot, M.P.; Ward, S.C. The Cambridge Structural Database. Acta Crystallogr. Sect. B Struct. Sci. Cryst. Eng. Mater. 2016, 72, 171–179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Bondi, A. Van der Waals Volumes and Radii. J. Phys. Chem. 1964, 68, 441–451. [Google Scholar] [CrossRef]
  41. Tiritiris, I.; Kantlehner, W. Crystal Structure of N,N,N′,N′-tetramethylchloroformamidinium chloride, [C5H12N2Cl]Cl. Z. Kristallogr. New Cryst. Struct. 2008, 223, 345–346. [Google Scholar] [CrossRef]
  42. Burg, A.B. Restudy of the Action of Sulfur Dioxide on Dry Trimethylamine Oxide: Iodine Oxidation and Lewis Acid Chemistry of the Most Reactive Product, (CH3)2(H)NCH2SO3. Inorg. Chem. 1989, 28, 1295–1300. [Google Scholar] [CrossRef]
  43. Clark, G.R.; Shaw, G.L.; Surman, P.W.J.; Taylor, M.J.; Steele, D. Preparation, Structure and Vibrational Spectrum of the Dimethylmethyleniminium Ion, Including the Role of Cationic Polymers in its Formation. J. Chem. Soc. Faraday Trans. 1994, 90, 3139–3144. [Google Scholar] [CrossRef]
  44. Silvi, B.; Savin, A. Classification of Chemical Bonds Based on Topological Analysis of Electron Localization Functions. Nature 1994, 371, 683–686. [Google Scholar] [CrossRef]
  45. Bartashevich, E.V.; Tsirelson, V.G. Interplay between Non-covalent Interactions in Complexes and Crystals with Halogen Bonds. Russ. Chem. Rev. 2014, 83, 1181–1203. [Google Scholar] [CrossRef]
  46. Kuznetsov, M.L. Can Halogen Bond Energy be Reliably Estimated from Electron Density Properties at Bond Critical Point? The Case of the (A)nZ–Y···X− (X, Y = F, Cl, Br) Interactions. Int. J. Quantum Chem. 2018. [Google Scholar] [CrossRef]
  47. Dunning, T.H., Jr. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  48. Wilson, A.K.; Woon, D.E.; Peterson, K.A.; Dunning, T.H., Jr. Gaussian Basis Sets for Use in Correlated Molecular Calculations. IX. The Atoms Gallium through Krypton. J. Chem. Phys. 1999, 110, 7667–7676. [Google Scholar] [CrossRef]
  49. Zhao, Y.; Truhlar, D.G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-class Functionals and 12 Other Functionals. Theoret. Chem. Acc. 2008, 120, 215–241. [Google Scholar] [CrossRef]
  50. Schmidt, M.W.; Baldridge, K.K.; Boatz, J.A.; Elbert, S.T.; Gordon, M.S.; Jensen, J.H.; Koseki, S.; Matsunaga, N.; Nguyen, K.A.; Su, S.; et al. General Atomic and Molecular Electronic Structure System. J. Comput. Chem. 1993, 14, 1347–1363. [Google Scholar] [CrossRef]
  51. Dovesi, R.; Saunders, V.R.; Roetti, C.; Orlando, R.; Zicovich-Wilson, C.M.; Pascale, F.; Civalleri, B.; Doll, K.; Harrison, N.M.; Bush, I.J.; et al. CRYSTAL14 User’s Manual; University of Torino: Torino, Italy, 2014. [Google Scholar]
  52. Becke, A.D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  53. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B Condens. Matter Mater. Phys. 1988, 37, 785–789. [Google Scholar] [CrossRef]
  54. Grimme, S. Semi-empirical GGA-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef] [PubMed]
  55. Dovesi, R.; Saunders, V.R.; Roetti, C.; Orlando, R.; Zicovich-Wilson, C.M.; Pascale, F.; Civalleri, B.; Doll, K.; Harrison, N.M.; Bush, I.J.; et al. CRYSTAL17 User’s Manual; University of Torino: Torino, Italy, 2017. [Google Scholar]
  56. Gatti, C.; Casassa, S. Topond14 User’s Manual; University of Torino: Torino, Italy, 2014. [Google Scholar]
  57. Keith, T.A. AIMALL, Version 12.06.03, 2012 Professional. Available online: http://aim.tkgristmill.com (accessed on 20 February 2019).
  58. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  59. Kostenetskiy, P.; Semenikhina, P. SUSU Supercomputer Resources for Industry and Fundamental Science. In Proceedings of the Global Smart Industry Conference (GloSIC), Chelyabinsk, Russia, 13–15 November 2018; pp. 1–7. [Google Scholar] [CrossRef]
Molecules 24 01083 g001 550
Figure 1. The distribution of the Hal−···CH3–Y interactions in crystals for Cl− (green), Br− (brown), I− (violet); the blue lines mark the sum of van der Waals radii.
Figure 1. The distribution of the Hal−···CH3–Y interactions in crystals for Cl− (green), Br− (brown), I− (violet); the blue lines mark the sum of van der Waals radii.
Molecules 24 01083 g001
Molecules 24 01083 g002 550
Figure 2. Fragments of structures with tetrel bonds and other non-covalent interactions in crystals LONGEB (a), VAPREJ (b), POSTUM02 (c), LILLOH (d). The interatomic distances are given in Angstroms.
Figure 2. Fragments of structures with tetrel bonds and other non-covalent interactions in crystals LONGEB (a), VAPREJ (b), POSTUM02 (c), LILLOH (d). The interatomic distances are given in Angstroms.
Molecules 24 01083 g002
Molecules 24 01083 g003 550
Figure 3. (a) ESP in the trimethylammonium cation on the isosurface of electron density of 0.003 a.u.; (b) contour map of ESP in the plane N-C-H, red point indicates the maximum of ESP on the van der Waals surface (red line) and belongs to C atomic basin.
Figure 3. (a) ESP in the trimethylammonium cation on the isosurface of electron density of 0.003 a.u.; (b) contour map of ESP in the plane N-C-H, red point indicates the maximum of ESP on the van der Waals surface (red line) and belongs to C atomic basin.
Molecules 24 01083 g003
Molecules 24 01083 g004 550
Figure 4. The ELF for (a) the halogen and tetrel bonds Cl(2)···Cl(1)−···C(4) in LONGEB crystal; (b) the tetrel bond in LILLOH crystal.
Figure 4. The ELF for (a) the halogen and tetrel bonds Cl(2)···Cl(1)−···C(4) in LONGEB crystal; (b) the tetrel bond in LILLOH crystal.
Molecules 24 01083 g004
Molecules 24 01083 g005 550
Figure 5. The disposition of electron density and electrostatic potential minima (a.u.) (a) along two interatomic lines Cl(2)···Cl(1)−···C(4), (Å), in LONGEB; (b) along the tetrel bond in LILLOH. The arrows point to the electrophilic site provider.
Figure 5. The disposition of electron density and electrostatic potential minima (a.u.) (a) along two interatomic lines Cl(2)···Cl(1)−···C(4), (Å), in LONGEB; (b) along the tetrel bond in LILLOH. The arrows point to the electrophilic site provider.
Molecules 24 01083 g005
Molecules 24 01083 g006 550
Figure 6. Binding energy in complexes vs minimum of electrostatic potential (a) and maximum of potential acting on an electron in a molecule, (b) along the line of tetrel bonds.
Figure 6. Binding energy in complexes vs minimum of electrostatic potential (a) and maximum of potential acting on an electron in a molecule, (b) along the line of tetrel bonds.
Molecules 24 01083 g006
Molecules 24 01083 g007 550
Figure 7. Potential acting on an electron in a molecule, а.u., (blue), electrostatic potential, а.u., (red) and electron density, а.u., (black) along the tetrel bond in (a) ZENJAD and (b) FADXIR complexes.
Figure 7. Potential acting on an electron in a molecule, а.u., (blue), electrostatic potential, а.u., (red) and electron density, а.u., (black) along the tetrel bond in (a) ZENJAD and (b) FADXIR complexes.
Molecules 24 01083 g007
Table
Table 1. Experimental and calculated tetrel and C–Y bond lengths, D (Å), angles Hal−···C–Y, θ and calculated electron density, ρ(rbcp) (a.u.) at bond critical points for considered crystals.
Table 1. Experimental and calculated tetrel and C–Y bond lengths, D (Å), angles Hal−···C–Y, θ and calculated electron density, ρ(rbcp) (a.u.) at bond critical points for considered crystals.
CrystalBondCrystal Dexp, θ (Hal−···C–N)expCrystal Dcalc, θ (Hal−···C–N)calcρ(rbcp), Crystal
GETQIFСl(3)− ···C(2)3.4584
169.08		3.4260
163.91		0.0056
C(2)–N(1)1.48151.49580.2441
LONGEBСl(1)− ···C(4)3.4251
175.28		3.4087
166.64		0.0068
C(4)–N(2)1.47221.47400.2458
VAPREJСl(1)− ···C(1)3.417
164.88		3.4385
164.49		0.0064
C(1)–N(1)1.4661.47470.2480
TMHYZCСl(1)− ···C(2)3.4374
174.96		3.4280
176.34		0.0062
C(2) ···N(1)1.49761.50800.2406
ZENJADСl(1)− ···C(7)3.5111
170.58		3.4644
171.72		0.0056
C(7)–O(2)1.44711.44680.2306
LILLOHBr(1)− ···C(2)3.533
167.25		3.5664
166.86		0.0061
C(2)–N(1)1.4741.47350.2490
FADXIRBr(1)− ···C(6)3.6014
170.87		3.5722
173.15		0.0058
C(6)–N(1)1.50251.50820.2371
POSTUM02Br(1)− ···C(1)3.7012
175.21		3.6667
174.15		0.0048
C(1)–N(1)1.48521.49280.2432
ZZZGVM01Br(1)− ···C(2)3.742
168.65		3.7283
169.04		0.0042
C(2)–N(1)1.4741.49680.2437
ZZZUQO03Br(1)− ···C(1)3.685
171.12		3.6819
171.11		0.0049
C(1)–N1.4871.50390.2411
Table
Table 2. Calculated binding energy Eb (kcal/mol), bond lengths Dcalc (Å), electron density ρ(rbcp) at bcp (a.u.) for Hal−···CH3 minimum of electrostatic potential ESPmin (а.u.) and maximum of potential acting on an electron in molecule PAEMmax (а.u.) along the line of tetrel bonds in complexes.
Table 2. Calculated binding energy Eb (kcal/mol), bond lengths Dcalc (Å), electron density ρ(rbcp) at bcp (a.u.) for Hal−···CH3 minimum of electrostatic potential ESPmin (а.u.) and maximum of potential acting on an electron in molecule PAEMmax (а.u.) along the line of tetrel bonds in complexes.
RefcodeTetrel BondEbDcalcρ(rbcp)ESPminPAEMmax
GETQIFСl−···CH3Y−81.072.82620.019−0.038−0.4819
LONGEBСl−···CH3Y−71.052.87820.017−0.063−0.4454
TMHYZCСl−···CH3Y−82.672.82480.019−0.034−0.4874
VAPREJСl−···CH3Y−81.802.82260.019−0.034−0.4870
ZENJADСl−···CH3Y−52.282.92680.015−0.107−0.3883
FADXIRBr−···CH3Y−81.092.98550.017−0.015−0.4741
LILLOHBr−···CH3Y−78.672.98960.016−0.020−0.4664
POSTUM02Br−···CH3Y−79.822.98020.017−0.019−0.4699
ZZZGVM01Br−···CH3Y−77.972.99490.0165−0.024−0.4612
ZZZUQO03Br−···CH3Y−76.633.00150.0162−0.028−0.4560

Share and Cite

MDPI and ACS Style

Bartashevich, E.; Matveychuk, Y.; Tsirelson, V. Identification of the Tetrel Bonds between Halide Anions and Carbon Atom of Methyl Groups Using Electronic Criterion. Molecules 2019, 24, 1083. https://doi.org/10.3390/molecules24061083

AMA Style

Bartashevich E, Matveychuk Y, Tsirelson V. Identification of the Tetrel Bonds between Halide Anions and Carbon Atom of Methyl Groups Using Electronic Criterion. Molecules. 2019; 24(6):1083. https://doi.org/10.3390/molecules24061083

Chicago/Turabian Style

Bartashevich, Ekaterina, Yury Matveychuk, and Vladimir Tsirelson. 2019. "Identification of the Tetrel Bonds between Halide Anions and Carbon Atom of Methyl Groups Using Electronic Criterion" Molecules 24, no. 6: 1083. https://doi.org/10.3390/molecules24061083

APA Style

Bartashevich, E., Matveychuk, Y., & Tsirelson, V. (2019). Identification of the Tetrel Bonds between Halide Anions and Carbon Atom of Methyl Groups Using Electronic Criterion. Molecules, 24(6), 1083. https://doi.org/10.3390/molecules24061083

Article Metrics

Back to TopTop