Next Article in Journal
Two New Cytotoxic Steroidal Alkaloids from Sarcococca Hookeriana
Next Article in Special Issue
Differential Binding of Tetrel-Bonding Bipodal Receptors to Monatomic and Polyatomic Anions
Previous Article in Journal
Molecular Repolarisation of Tumour-Associated Macrophages
 
 
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 1
10.3390/molecules24010010
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

Tetrel Bond between 6-OTX3-Fulvene and NH3: Substituents and Aromaticity

by
Ming-Chang Hou
1,
Shu-Bin Yang
1,
Qing-Zhong Li
1,*,
Jian-Bo Cheng
1,
Hai-Bei Li
2,* and
Shu-Feng Liu
3
1
The Laboratory of Theoretical and Computational Chemistry, School of Chemistry and Chemical Engineering, Yantai University, Yantai 264005, China
2
School of Ocean, Shandong University, Weihai 264209, China
3
Shandong Key Laboratory of Biochemical Analysis, College of Chemistry and Molecular Engineering, Qingdao University of Science and Technology, Qingdao 266042, China
*
Authors to whom correspondence should be addressed.
Molecules 2019, 24(1), 10; https://doi.org/10.3390/molecules24010010
Submission received: 3 December 2018 / Revised: 16 December 2018 / Accepted: 17 December 2018 / Published: 20 December 2018
(This article belongs to the Special Issue σ and π Holes: A New Class of Non-Covalent Interactions)
Browse Figures
Versions Notes

Abstract

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

1. Introduction

Fulvenes serve not only as synthetic precursors to naturally occurring compounds [1,2], but also as starting materials for the synthesis of novel substituted titanocenes, which are bio-organometallic anti-cancer drugs [3]. As an isomer of benzene, the structures and properties of fulvene are different from those of benzene. Fulvene is a non-alternant hydrocarbon, although it also displays a planar structure. Unlike benzene, this planar molecule is unstable both thermodynamically and kinetically, with very low resonance energy [4] and aromatic stabilization energy [5]. Despite fulvene being a non-aromatic molecule, the realization of its aromaticity has attracted much interest [6,7,8,9,10]. At its lowest excited states, fulvene is moderately aromatic, with a substantial contribution from the dipolar aromatic structure [6]. The electronic structure of the conjugated system in fulvene suffers substantial changes when exocyclic substitution occurs; its aromaticity increases by electron-donating substituents, and decreases by electron-withdrawing groups [7,8,9]. Interestingly, when a lithium atom approaches a face of the fulvene molecule, the fulvene moiety becomes aromatic, with an ‘aromatic’ NICS (nuclear independent chemical shift) value of −11 ppm [10]. When a fulvene ring is fused into the phenyl ring of an O-hydroxy Schiff base, the strength of an intramolecular hydrogen bond is tuned by the substituent at the exocyclic carbon in fulvene [11]. An electron-donating group at the exocyclic carbon in fulvene weakens the intramolecular hydrogen bond in some of the fulvene’s Schiff bases [12]. Cyano groups adjoined to the fulvene ring greatly increase the acidity of 6-OH-fulvene, resulting in organic superacids [13].
Recently, tetrel bonding, an attractive interaction between a Group IV atom and a Lewis base, has been attracting much attention [14,15,16,17,18,19,20]. Similarly to hydrogen and halogen bonds, tetrel bonding has broad applications in supermolecular materials [21,22,23] and chemical reactions [24,25,26]. Owing to the greater electronegativity and lower polarization, carbon atoms seldom engages in tetrel bonding, wherein it is also called a carbon bond if the carbon atom acts as a Lewis acid. Even so, carbon bonding has been observed in different systems [27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43]. A weak Ar···C interaction was firstly detected in an Ar···propargyl alcohol complex by microwave spectroscopy and ab initio calculations [27]. Mani and Arunan then performed a theoretical study of carbon bonding in complexes of methanol and methyl fluoride, where a methyl group is adjoined with an electron-withdrawing group, and found that in some neutral bases the acidity of the carbon atom increases, and the interaction energy amounts to ~8 kJ/mol [28]. The ubiquity of such carbon bond was further witnessed in the solid state by means of Cambridge Structural Database (CSD) and charge density analyses [29]. Similarly, the –CF3 group in para-substituted ArCF3 is able to participate in carbon bonding [30]. The –CF3 group in pyridine–CF3 and furan–CF3 is favorable for a weak hydrogen bond with NH3, but the protonation on the N and O atoms in the pyridine and furan rings results in a moderate carbon bond with NH3 [31]. A similar protonation enhancing effect was also reported for the –CH3 group [32]. In addition to sp3-hybridized carbon being involved in carbon bonding, the sp2-hybridized carbon in C=C and C=O bonds also acts as a Lewis acid [33,34,35,36,37,38,39,40,41]. Moreover, the sp3-hybridized carbon forms a weaker carbon bond than does the sp2-hybridized carbon in most cases. How can a carbon bond be strengthened when a sp3-hybridized carbon is acting as a Lewis acid?
It has been demonstrated that the hydroxyl proton of 6-OH-fulvene is acidic, and that its acidity increases greatly when cyano groups adjoin to the fulvene ring [13]. Inspired by these findings, in this study, the hydroxyl H atom of 6-OH-fulvene was changed to –CX3 (X = H and F), which was able to form a stronger carbon bond with NH3. In order to compare the strength of the tetrel bonds formed by different tetrel atoms, the –CX3 group in 6-OCX3-fulvene was also changed to –SiX3 and –GeX3. The four H atoms in the 2–5 positions of fulvene were then replaced by four cyano groups to enhance the acidity of CF3 and SiF3. For comparison, the corresponding complexes of PhOTX3 and HOTX3 (T = C, Si, and Ge) were also studied. We focused on the effect of both TX3 and the formation of tetrel bonds on the aromaticity of the fulvene ring.

2. Theoretical Methods

The second-order Moller-Plesset perturbation theory (MP2) method was utilized to optimize the structures of the complexes and monomers, using the aug-cc-pVTZ basis set. Frequency calculations were performed at the same level, to ensure they were at minima on the potential energy surfaces. This level of the theory has been often used to study tetrel-bonded complexes [15,16,25,31]. The stability of the complex was estimated, using the interaction and binding energies, by the supramolecular method (subtracting the energies of the monomers from the energy of the complex). Both interaction and binding energies were corrected for the basis set superposition error (BSSE) with the counterpoise method proposed by Boys and Bernardi [44]. For the interaction energy, the geometry of the monomer was the one in its complex, while for the binding energy, in contrast, the monomers were fully optimized. The magnitude of their difference is the deformation energy (DE), which was used to measure the deformation of both monomers as they came close to each other and formed the complexes. All calculations were carried out using the Gaussian 09 package of codes [45].
The molecular electrostatic potentials (MEPs) were calculated on the 0.001 electrons bohr−3 isosurface at the MP2/aug-cc-pVDZ level, using the WFA-SAS program [46]. The topological parameters of the interaction, including electron density, its Laplacian, and the energy density at the bond critical point (BCP), were analyzed by the AIM2000 [47] program. The natural bond orbital (NBO) analysis was carried out at the HF/aug-cc-pVTZ level, while the charge transfer and second order perturbation energy were obtained by NBO3.1 version [48] implemented in the Gaussian 09 program. The orbital interaction contribution was obtained by means of an analysis of natural orbital for chemical valence (NOCV), using the ADF program [49]. Interaction energy was decomposed by the LMOEDA method [50] at the MP2/aug-cc-pVTZ level, using the GAMESS program [51].
The NICS(1)zz index [52] was calculated at the MP2/aug-cc-pVDZ level as the value of ZZ component of magnetic shielding tensor (taken with negative sign) at the point located 1 Å above a ghost atom located in the geometric center of the fulvene ring.

3. Results and Discussion

3.1. MEPs of 6-OTX3-Fulvene

The MEP maps of 6-OTX3-fulvene with two views are presented in Figure 1. It is evident from Figure 1 that four σ-holes were found on the surface of the –TX3 (T = C, Si, and Ge) group. The MEP values of the four σ-holes varied between the different –TX3 groups. For the –TH3 (T = C, Si and Ge) group, the σ-hole along the O–T bond had the largest MEP, followed by the one along the H–T bond, which was coplanar with the C(sp2)-H bond. The MEP values of the left two σ-holes were approximately equal, being attributed to the Cs symmetrical geometries of 6-OTH3-fulvene. The F substitution in 6-OTF3-fulvene caused a complicated change in the order of the σ-hole magnitude, which was primarily due to the conformation change of the –TF3 group. The MEP of the σ-hole at the O–T end was not largest for the –TF3 (T = Si and Ge) group, and the location of the largest σ-hole was different for the –SiF3 and –GeF3 groups. The σ-holes of the –CX3 groups had a similar distribution to that of the –TH3 (T = Si and Ge) groups. It was noted that the values of 0.037 a.u. and 0.039 a.u. of –CH3 group corresponded to the three H atoms of this group, while the value of 0.057 au of –CF3 group corresponded to the H atom of C(sp2)–H bond.
The MEP of the σ-hole at the O–T end is listed in the last column of Table 1. It is evident that the MEP of this σ-hole showed an increasing trend from C to Si to Ge with the increase of the T atomic mass, due to the lower electronegativity and greater polarizability of the heavier T atom. Furthermore, with the F substituent, the σ-hole was enlarged, owing to the electron-withdrawing ability of F atoms, with the exception of 6-OCF3-fulvene, where it was slightly smaller than that of 6-OCH3-fulvene.

3.2. Geometries and Interaction Energy of Complexes

Figure 2 shows the optimized structures of the 6-OTX3-fulvene···NH3 complexes (T = C, Si, and Ge; X = H and F). In most complexes except 6-OCF3-fulvene···NH3, the angle ∠O–T···N was close to 180°, and three T–X bonds displayed a stagger conformation with three N–H bonds. In 6-OCF3-fulvene···NH3, two H atoms of NH3 pointed to the lone pair of two F atoms of 6-OCF3-fulvene. Even so, both molecules were combined through vdW interactions in 6-OCF3-fulvene···NH3. According to studies on the complexes of α/β-furanCF3 and p-PyCF3 with NH3 [31], they still had the capacity to form a tetrel bond interaction, as in 6-OCF3-fulvene···NH3. The Si···N and Ge···N distances were much shorter than the sum of their van der Waals radii (3.6 Å and 3.9 Å, respectively), and the F substituents further shortened their distance, especially for the Ge···N interaction.
The second column of Table 1 is the interaction energy (Eint) of the complex calculated by the formula: Eint = E(complex) − E(monomer1) − E(monomer2), in which E(monomer) is the energy of the monomer within the complex structure. Eint was −9.93 kJ/mol in 6-OCH3-fulvene···NH3, which is twice as high as that of HOCH3···NH3. Thus, introducing a fulvene group to the carbon atom was an efficient method for enhancing the strength of tetrel bonding. In 6-OGeF3-fulvene···NH3, Eint was more negative than that in 6-OSiF3-fulvene···NH3. Without the F substituent, a reverse result was found, that is, the interaction in 6-OSiH3-fulvene···NH3 was stronger than that in 6-OGeH3-fulvene···NH3. This was inconsistent with the positive MEP of the Si and Ge atoms, indicating that electrostatic interaction is not the only dominant factor in the formation of a tetrel bond. Due to their strong electron-withdrawing ability, the F substituents in the tetrel donor molecule resulted in an increase of interaction energy by 87 and 112 kJ/mol for the Si and Ge complexes, respectively. It was found that the degree of interaction energy increased in that the effect of F substitution was more prominent for the Ge complex. The binding energy had a change similar to the interaction energy, although the former was smaller.
The geometrical deformation of tetrel donor molecules is often observed in tetrel-bonded complexes [53], where the tetrel atom is inclined to form a pentahedral structure, and thus to accommodate the approaching base. This geometrical deformation is measured by deformation energy (DE), which is the energy required to distort the interacting monomers from their equilibrium structures to the geometries in the complex. DE is negligible in 6-OCH3-fulvene···NH3 and 6-OCF3-fulvene···NH3, but is very large in the Si and Ge complexes, up to 99 kJ/mol in 6-OSiF3-fulvene···NH3. Dependent on the nature of the tetrel donor, the value of DE contributes to 24–72% of the interaction energy. It is evident from Table 1 that the silicon donor had larger DE than the Ge analogue. Furthermore, the F substituents in both Si and Ge tetrel donors led to a significant larger DE compared to that in 6-OTH3-fulvene···NH3.
The O–T bond was elongated during the formation of the complexes. The O–T bond elongation increased in the order of C < Ge < Si for a given –OTH3/–OTF3 group. This order was consistent with the interaction energy in 6-OTH3-fulvene···NH3, but an inconsistency was found in 6-OTF3-fulvene···NH3. This indicated that the O–T bond elongation is not a good indicator for the strength of a tetrel bond, different to that in H-bonding.
In order to get a deeper insight into the tetrel bond in 6-OTX3-fulvene···NH3, we introduced HOTX3 and PhOTX3 for comparison. We found that the binding modes with NH3 were the same for 6-OTX3-fulvene, HOTX3, and PhOTX3. The interaction strength was, in order, HOTX3 < PhOTX3 < 6-OTF3-fulvene, which was evidenced by the magnitudes of the intermolecular contact, the negative interaction energy, and the O–T bond elongation. This demonstrates that connecting with a benzene or fulvene ring is favorable to engage in a tetrel bond for a tetrel atom, especially with a fulvene group. It is apparent from Table 1 that the benzene or fulvene ring had a greater effect on the tetrel bond of an –OTH3 group than on that of an –OTF3 group. For example, when changing H to fulvene, the interaction energy of –OGeH3 complex increased by 175%, while that of –OGeF3 complex increased by 19%.

3.3. Substituents and AIM Analysis

The interaction mode between 6-OTX3-fulvene and NH3 was further determined by AIM analysis. Figure 3 shows the AIM diagrams of –CF3(4) and –SiF3(5). The other systems, except –CF3, were similar to 5, with only one BCP between the T atom and the N atom, confirming tetrel bond formation. Three BCP paths were found in 6-OCF3-fulvene···NH3, with one linear F···N path and two curved F···N paths in the vicinity of the N atom. The coexistence of three paths with a small interaction energy (~4 kJ/mol) implied that the interaction in this complex was a van der Waals force.
Table 2 presents the electron density, Laplacian, and total energy density at the BCP of the complexes, in which the mean values of the three F···N BCPs are given for 6-OCF3-fulvene···NH3. It confirms that the interaction type can be classified according to the sign of Laplacian and the total energy density [54]. Both were positive for 6-OCH3-fulvene···NH3 and 6-OCF3-fulvene···NH3, indicating, together with the small electron density, that both were weak close-shell interactions. For other complexes involving heavy tetrel donors, both the positive Laplacian and the negative total energy density demonstrated that the tetrel bond was a partially covalent interaction. The F substituents resulted in a further larger electron density, and thus gave rise to a stronger tetrel bond compared to the no substituent complexes. Moreover, the F substitution made the increase of electron density at the Ge···N BCP larger relative to the increase of electron density at the Si···N BCP. The enhancement of the electron density for the same T···N tetrel bond was accompanied by a greater Laplacian and more negative total energy density.
Through analysis of the optimized structure (Figure 2) and AIM maps, we can know that the –CF3 group of 6-OCF3-fulvene forms vdW interactions, not a tetrel bond, when it binds with NH3, like many molecules involving a –CF3 group [31]. Thus, a question is raised: is it possible for the –CF3 group of 6-OCF3-fulvene to form a tetrel bond with NH3? To answer it, the four H atoms on the five-membered ring of 6-OCF3-fulvene were replaced by the electron-withdrawing group, CN. The AIM analysis of the CN–substituent systems was calculated at the MP2/aug-cc-pVDZ level, together with the unsubstituted analogues—the corresponding AIM diagrams are also presented in Figure 3. Three intermolecular BCPs were found between the three F atoms of the –CF3 group and the N atom of NH3, which was also observed in the CF4···NCH tetrel-bonded complex [24]. This is clear evidence that the C···N tetrel bond exists in the complex of CN–substituent 6-OCF3-fulvene. Due to the strong electron-withdrawing effect of the CN group, the positive MEP on the σ-hole of –CF3 group was doubled from 0.039 a.u. to 0.072 a.u. Thus, the electron-withdrawing CN groups in the tetrel donor molecule can modulate the variation from vdW interactions to the tetrel bond. A similar effect was also realized by protonation of the tetrel donor molecule [31] or the methyl substituents of NH3, where the former similarly resulted in enlargement of the σ-hole of –CF3 group, and the latter not only led to an increase in basic electron donation, but also avoided the formation of the H-bond. The C···N separation was 3.291 Å in the CN–substituent complex, which was shortened by 0.161 Å relative to that in 6-OCF3-fulvene···NH3. As a consequence, the interaction energy increased from −2.98 to −9.63 kJ/mol, three times higher than that of vdW interactions.
To estimate the enhancing effect of CN substituents on the strength of tetrel bond, the four H atoms of the five-membered ring in 6-OSiF3-fulvene were also substituted by CN groups. This substitution caused an increase of the positive MEP on the σ-hole of the Si atom from 0.095 a.u. to 0.140 a.u., resulting in an increase in the acidity of Si atom. The tetrel bond was characterized by a Si···N BCP, with an electron density of 0.0679 a.u. and a negative total energy density. The Si···N distance was shortened from 2.035 Å to 1.992 Å, and the interaction energy increased from −143.62 to −190.13 kJ/mol. These data indicate that strong electron-withdrawing groups in tetrel donor molecules have a prominent enhancement effect on the strength of tetrel bonding.

3.4. NBO Analysis

Orbital interactions and charge transfer analysis were performed for the 6-OTX3-fulvene···NH3 complexes (Table 3). In 6-OCH3-fulvene···NH3, there was an orbital interaction of Lp(N)→σ*C-O, where Lp(N) denotes the lone pair orbital of the N atom and σ*C-O is the anti-bonding orbital of the C–O bond. This orbital interaction was very weak, with a perturbation energy of less than 6 kJ/mol. Additionally, the charge transfer was very small, consistent with the small interaction energy of this complex. Similarly, the orbital interaction Lp(F)→σ*N-H in 6-OCF3-fulvene···NH3 was negligible (Table 3). This confirms its weak vdW interactions. For 6-OXF3-fulvene···NH3 (X = Si and Ge), there were two types of orbital interactions (Lp(N)→σ*X-O and Lp(N)→σ*X-F), and the perturbation energy of Lp(N)→σ*X-F listed in Table 3 is the sum of the three orbital interactions. The contribution of the Lp(N)→σ*X-F orbital interaction was much larger than that of Lp(N)→σ*X-O, and both types of orbital interactions were stronger in 6-OGeF3-fulvene···NH3 than those in 6-OSiF3-fulvene···NH3, consistent with the interaction energy. Strong orbital interactions were related to high charge transfer (>0.16 e) in 6-OXF3-fulvene···NH3(X = Si and Ge). For 6-OXH3-fulvene···NH3 (X = Si and Ge), the Lp(N)→σ*X-O orbital interaction was changed into Lp(N)→p*X (the empty lone pair orbital on X atom), and Lp(N)→σ*X-H was not changed. The Lp(N)→p*X orbital interaction made a larger contribution than Lp(N)→σ*X-H, and both orbital interactions were stronger in 6-OSiH3-fulvene···NH3. Obviously, the F substituents significantly affected the type of orbital interaction as well as its strength.
In order to visually understand the contribution of orbital interactions to bonding energy, we performed an energy decomposition calculation, in conjunction with an analysis of natural orbital for chemical valence (NOCV), with the ADF 2008.01 program [49]. Figure 4 plots the deformation densities due to the pair-wise orbital interactions. The blue and red regions represent the increase and decrease of densities, respectively. It is obvious that NH3 was surrounded by red areas, while the group OTX3 was surrounded by blue regions. Thus, the charge flow shifted from NH3 (a base) to OTX3 (an acid). The largest contribution was from the orbital interaction between Lp(N) and σ*X-O in 6-OXH3-fulvene···NH3 (X = Si and Ge), while it was from that between Lp(N) and σ*X-F in 6-OXF3-fulvene···NH3 (X = Si and Ge). The former was larger for the SiH3 system, while the latter was larger for the GeF3 system. The stronger orbital interaction resulted in larger deformation densities. In TH3 complexes, the orbital interactions were relatively important; this can explain the stronger binding for Si than for Ge, although the positive MEP of the σ-hole for Ge was more positive. In TF3 complexes, the fluorines produced a stronger σ-hole, making the electrostatic interactions more important, and so the Ge was the stronger binder. This interplay between the electrostatic and orbital interactions has been previously suggested to be able to reverse the binding order [55].

3.5. Energy Decomposition Analysis

To gain a deeper insight into the characteristics of the tetrel bond in the complexes 6-OTX3-fulvene···NH3, its interaction energy was divided into five components: electrostatic energy (Eele), repulsion energy (Erep), exchange energy (Eex), polarization energy (Epol), and dispersion energy (Edisp). Because both Erep and Eex cancel each other out, only three attractive terms (Eele, Epol, and Edisp) are plotted in Figure 5, for ease of comparison. Each energy contribution was very small in the CX3 system, consistent with the weak characteristic of the interaction. In the SiX3 and GeX3 systems, the electrostatic energy was largest, followed by the polarization energy, and the dispersion energy was smallest. Both electrostatic and polarization contributions decreased from SiH3 to GeH3 systems, but increased from SiF3 to GeF3 systems (Figure 5), which illustrates the large effect of the F substituents on the nature of the tetrel bond. Owing to its smaller electronegativity, the Si atom in the SiH3 group was more easily polarized by a base than was the Ge atom in the GeH3 group. The greater polarization consequently resulted in a larger electrostatic interaction between both monomers, which was inconsistent with the magnitude of the positive MEP on the σ-hole along the O–T axis. The F substituents enlarged the MEP value of the σ-hole along the O–T axis. Thus, the electrostatic energy increased in the TF3 system, up to −90 kJ/mol in 6-OGeF3-fulvene···NH3, and was accompanied by a greater polarization. The large polarization energy means that the shape of the molecular orbital underwent a large change during the formation of the complexes, which is consistent with the deformation density of the pair-wise orbital interaction. This is consistent with the partially covalent nature of the tetrel bonding interactions of Si and Ge that was ascertained by analyzing the topological properties at the intermolecular BCPs.

3.6. Aromaticity of the Fulvene Ring

It has been confirmed that an electron-donating group at the exocyclic carbon atom of fulvene can increase its π-electron density [7,8,9]. Therefore, it is interesting to study the influence of a tetrel bond formation on the aromaticity of fulvene ring. Aromaticity is judged with NICS(1)zz at the point located 1 Å above the fulvene ring center. A more negative NICS(1)zz indicates a larger aromaticity. This value is collected in Table 4 for 6-OTX3-fulvene···NH3 and 6-OTX3-fulvene. It is evident that the NICS(1)zz value was more negative in the complex than that in the isolated molecule, indicating that the formation of a tetrel bond resulted in an increase of aromaticity for the fulvene ring. It was confirmed that the charge transfer in the formation of the tetrel bond is from the lone pair on the N atom of NH3 to 6-OTX3-fulvene. This leads to an increase in the aromaticity and electron density of the fulvene ring. That is, the formation of a tetrel bond increased the aromaticity of the fulvene ring, similar with the electron-donating group at the exocyclic carbon atom. Furthermore, the increase of the NICS(1)zz value was consistent with the charge transfer in SiX3 and GeX3 systems. For example, 6-OSiH3-fulvene···NH3 had a larger increase of NICS(1)zz value and charge transfer than 6-OGeH3-fulvene···NH3.

4. Conclusions

The complexes of 6-OTX3-fulvene···NH3 (T = C, Si, Ge; X = H, F) were studied using theoretical calculations, with regard to geometrics, energetics, charge transfer, orbital interactions, and AIM parameters. The main conclusions are summarized as:
(1) The interaction energy of 6-OCH3-fulvene···NH3 was up to −10 kJ/mol, which was twice as high as that in HOCH3···NH3, thus, carbon bonding is enhanced when –OCH3 combines with fulvene ring.
(2) The σ-hole at the T–O bond end was larger in the –OGeX3 group than that in the –OSiX3 group, thus, –OGeF3 forms a stronger tetrel bond than –OSiF3. However, –SiH3 engaged in a stronger tetrel bond than –OGeH3, which is inconsistent with the magnitude of the σ-hole on both atoms. This inconsistency can be partly attributed to the stronger orbital interaction in 6-OSiH3-fulvene···NH3, accompanied by prominent distortion of –SiH3 group.
(3) The cyano groups adjoined to the fulvene ring increased the positive MEP on the σ-hole of –CF3 and –SiF3 groups, thus, the vdW interactions in 6-OCF3-fulvene···NH3 (4) were changed to carbon bonding in CN-4, and the interaction energy of tetrel bond increased by 32%, up to −190 kJ/mol in CN-5. Such CN substitution has a prominent effect on the type and strength of interactions.
(4) The formation of a tetrel bond was accompanied by a charge transfer from the lone pair on the N atom of NH3 to 6-OTX3-fulvene, thus, it resulted in an increase in the electron density on the fulvene ring and its aromaticity.

Author Contributions

Conceptualization, Q.-Z.L. and H.-B.L.; methodology, S.-B.Y.; software, J.-B.C.; validation, M.-C.H., Q.-Z.L. and H.-B.L.; formal analysis, M.-C.H.; investigation, M.-C.H.; resources, J.-B.C.; data curation, S.-F.L.; writing—original draft preparation, M.-C.H.; writing—review and editing, Q.-Z.L. and H.-B.L.; visualization, S.-F.L.; supervision, J.-B.C.; project administration, S.-B.Y.; funding acquisition, Q.-Z.L., H.-B.L., and S.-F.L.

Funding

This work was supported by the National Natural Science Foundation of China (21573188 and 21403127) and the Open Subject of Faculty of Chemistry of QingDao University of Science and Technology (QUSTHX201807).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Stone, K.J.; Little, R.D. An exceptionally simple and efficient method for the preparation of a wide variety of fulvenes. Cheminform 1984, 15, 1849–1853. [Google Scholar] [CrossRef]
  2. Peloquin, A.J.; Stone, R.L.; Avila, S.E. Synthesis of 1,3-diphenyl-6-alkyl/aryl-substituted fulvene chromophores: Observation of π-π interactions in a 6-pyrene-substituted 1,3-diphenylfulvene. J. Org. Chem. 2012, 77, 6371–6376. [Google Scholar] [CrossRef] [PubMed]
  3. Strohfeldt, K.; Tacke, M. Bioorganometallic fulvene-derived titanocene anti-cancer drugs. Chem. Soc. Rev. 2008, 37, 1174–1187. [Google Scholar] [CrossRef] [PubMed]
  4. Schaad, L.J.; Hess, B.A., Jr. Hueckel molecular orbital π resonance energies. Question of the σ structure. J. Am. Chem. Soc. 1972, 94, 3068–3074. [Google Scholar] [CrossRef]
  5. Cyrański, M.K. Energetic aspects of cyclic pi-electron delocalization: Evaluation of the methods of estimating aromatic stabilization energies. Chem. Rev. 2005, 105, 3773–3811. [Google Scholar] [CrossRef]
  6. Malar, E.J.P.; Neumann, F.; Jug, K. Investigation of aromaticity in the excited states of fulvene. J. Mol. Struct. Theochem. 1995, 336, 81–84. [Google Scholar] [CrossRef]
  7. Stepień, B.T.; Krygowski, T.M.; Cyrański, M.K. Aromaticity strongly affected by substituents in fulvene and heptafulvene as a new method of estimating the resonance effect. Chem. Phys. Lett. 2001, 350, 537–542. [Google Scholar] [CrossRef]
  8. Krygowski, T.M.; Palusiak, M.; Oziminski, W.P. Aromaticity of substituted fulvene derivatives: Substituent-dependent ring currents. Phys. Chem. Chem. Phys. 2010, 12, 10740–10745. [Google Scholar] [CrossRef]
  9. Shakerzadeh, E.; Noorizadeh, S. Aromaticity study on tria-, penta-and hepta-fulvene derivatives. Comput. Theor. Chem. 2011, 964, 141–147. [Google Scholar]
  10. Krygowski, T.M.; Soncini, A.; Oziminski, W.P.; Fowler, P.W. Aromatization of fulvene by complexation with lithium. Org. Lett. 2010, 12, 4880–4883. [Google Scholar]
  11. Krygowski, T.M.; Zachara-horeglad, J.E.; Pelloni, S.; Lazzeretti, P.; Palusiak, M. Relation between π-electron localization/delocalization and H-bond strength in derivatives of omicron-hydroxy-schiff bases. J. Org. Chem. 2008, 73, 2138–2145. [Google Scholar] [CrossRef] [PubMed]
  12. Vatanparast, M.; Nekoei, A.R. An intramolecular hydrogen bond study in some Schiff bases of fulvene: A challenge between the RAHB concept and the σ-skeleton influence. New J. Chem. 2014, 38, 5886–5891. [Google Scholar]
  13. Vianello, R.; Maksić, Z.B. Tailoring of strong neutral organic superacids: DFT-B3LYP calculations on some fulvene derivatives. New J. Chem. 2004, 28, 843–846. [Google Scholar]
  14. 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]
  15. Guo, X.; Yang, X.; Li, W.Z.; Cheng, J.; Li, Q.Z. A σ-hole interaction with radical species as electron donors: Does single-electron tetrel bonding exist? Phys. Chem. Chem. Phys. 2014, 16, 11617–11625. [Google Scholar]
  16. Guo, X.; Liu, Y.W.; Li, W.Z.; Cheng, J.B.; Li, Q.Z. Competition and cooperativity between tetrel bond and chalcogen bond in complexes involving F2CX (X=Se and Te). Chem. Phys. Lett. 2015, 620, 7–12. [Google Scholar] [CrossRef]
  17. Bauzá, A.; Mooibroek, T.J.; Frontera, A. Tetrel bonding interactions. Chem. Rec. 2016, 16, 473–487. [Google Scholar] [CrossRef]
  18. Bauzá, A.; Mooibroek, T.J.; Frontera, A. Influence of ring size on the strength of carbon bonding complexes between anions and perfluorocycloalkanes. Phys. Chem. Chem. Phys. 2014, 16, 19192–19197. [Google Scholar] [CrossRef]
  19. Scheiner, S. Assembly of effective halide receptors from components. Comparing hydrogen, halogen, and tetrel Bonds. J. Phys. Chem. A 2017, 121, 3606–3615. [Google Scholar] [CrossRef]
  20. 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]
  21. Servati, G.M.; Stilinović, V.; Bauzá, A.; Frontera, A.; McArdle, P.; Derveer, D.V.; Mahmoudi, G. Design of lead(II) metal-organic frameworks based on covalent and tetrel bonding. Chem. Eur. J. 2015, 21, 17951–17958. [Google Scholar] [CrossRef] [PubMed]
  22. Bauzá, A.; Frontera, A.; Garczarek, P.; Stilinović, V.; Kirillov, A.M.; Kennedy, A.; Ruiz-Pérez, C.; Mahmoudi, G. Metal–organic and supramolecular lead (ii) networks assembled from isomeric nicotinoylhydrazone blocks: The effects of ligand geometry and counter-ion on topology and supramolecular assembly. CrystEngComm 2016, 18, 5375–5385. [Google Scholar]
  23. Bauzá, A.; Frontera, A.; Mahmoudi, G. Concurrent agostic and tetrel bonding interactions in lead(ii) complexes with an isonicotinohydrazide based ligand and several anions. Dalton Trans. 2016, 45, 4965–4969. [Google Scholar]
  24. 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] [PubMed]
  25. Liu, M.; Cheng, J.; Li, H.B.; Li, Q.Z. Tetrel bond of pseudohalide anions with XH3F (X=C, Si, Ge, and Sn) and its role in SN2 reaction. J. Chem. Phys. 2016, 145, 224310. [Google Scholar] [CrossRef] [PubMed]
  26. Elguero, J.; Del Bene, J.E.; Alkorta, I. Carbon–carbon bonding between nitrogen heterocyclic carbenes and CO2. J. Phys. Chem. A 2017, 121, 8136–8146. [Google Scholar]
  27. Mani, D.; Arunan, E. Microwave spectroscopic and atoms in molecules theoretical investigations on the Ar···propargyl alcohol complex: Ar···H-O, Ar···π, and Ar···C interactions. ChemPhysChem 2013, 14, 754–763. [Google Scholar] [CrossRef]
  28. 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]
  29. 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]
  30. Bauzá, A.; Mooibroek, T.J.; Frontera, A. Non-covalent sp3 carbon bonding with ArCF3 is analogous to CH-π interactions. Chem. Commun. 2014, 50, 12626–12629. [Google Scholar] [CrossRef]
  31. Liu, M.; Li, Q.; Scheiner, S. Comparison of tetrel bonds in neutral and protonated complexes of pyridine TF3 and furan TF3 (T=C, Si, and Ge) with NH3. Phys. Chem. Chem. Phys. 2017, 19, 5550–5559. [Google Scholar] [CrossRef] [PubMed]
  32. Bauzá, 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]
  33. Bučar, D.K.; Halasz, I.; Stilinović, V. V=O···C interactions in crystal structures of oxovanadium-coordination compounds. New J. Chem. 2013, 37, 619–623. [Google Scholar]
  34. Azofra, L.M.; Scheiner, S. Tetrel, chalcogen, and CH···O hydrogen bonds in complexes pairing carbonyl-containing molecules with 1, 2, and 3 molecules of CO2. J. Chem. Phys. 2015, 142, 034307. [Google Scholar] [CrossRef] [PubMed]
  35. Alkorta, I.; Elguero, J.; Del Bene, J.E. Azines as electron-pair donors to CO2 for N···C tetrel bonds. J. Phys. Chem. A 2017, 121, 8017–8025. [Google Scholar] [CrossRef] [PubMed]
  36. Cormanich, R.A.; Rittner, R.; O’Hagan, D.; Buhl, M. Inter- and intramolecular CF···C=O interactions on aliphatic and cyclohexane carbonyl derivatives. J. Comput. Chem. 2016, 37, 25–33. [Google Scholar] [CrossRef] [PubMed]
  37. Bauzá, A.; Frontera, A. Theoretical study on σ-and π-hole carbon···carbon bonding interactions: Implications in CFC chemistry. Phys. Chem. Chem. Phys. 2016, 18, 32155–32159. [Google Scholar] [CrossRef] [PubMed]
  38. Dong, W.; Yang, X.; Cheng, J.; Li, Q. Comparison for σ-hole and π-hole tetrel-bonded complexes involving F2C=CFTF3 (T=C, Si, and Ge): Substitution, hybridization, and solvation effects. J. Fluorine Chem. 2018, 207, 38–44. [Google Scholar] [CrossRef]
  39. Wei, Y.; Scheiner, S.; Li, Q. The π-tetrel bond and its influence on hydrogen bonding and proton transfer. ChemPhysChem 2018, 19, 736–743. [Google Scholar] [CrossRef] [PubMed]
  40. Xu, H.; Cheng, J.; Yang, X.; Li, Q. Comparison of σ-hole and π-hole tetrel bonds formed by pyrazine and 1,4-dicyanobenzene: The interplay between anion-π and tetrel bonds. ChemPhysChem 2017, 18, 2442–2450. [Google Scholar] [CrossRef] [PubMed]
  41. Quiñonero, D. Sigma-hole carbon-bonding interactions in carbon-carbon double bonds: An unnoticed contact. Phys. Chem. Chem. Phys. 2017, 19, 15530–15540. [Google Scholar] [CrossRef] [PubMed]
  42. Karim, A.; Schulz, N.; Andersson, H.; Nekoueishahraki, B.; Carlsson, A.C.C.; Sarabi, D.; Valkonen, A.; Rissanen, K.; Gräfenstein, J.; Keller, S.; et al. Carbon’s three-center-four-electron tetrel bond, treated experimentally. J. Am. Chem. Soc. 2018. [Google Scholar] [CrossRef] [PubMed]
  43. Pal, P.; Konar, S.; Lama, P.; Kinsuk, D.; Bauzá, A.; Frontera, A.; Mukhopadhyay, S. On the importance of noncovalent carbon-bonding interactions in the stabilization of a 1D Co(II) polymeric chain as a precursor of a novel 2D coordination polymer. J. Phys. Chem. B 2016, 120, 6803–6811. [Google Scholar] [CrossRef] [PubMed]
  44. Boys, S.F.; Bernardi, 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] [CrossRef]
  45. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A. Gaussian 09, Revision B.01; Gaussian Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  46. Bulat, F.A.; Toro-Labbé, A.; Brinck, T.; Murray, J.S.; Politzer, P. Quantitative analysis of molecular surfaces: Areas, volumes, electrostatic potentials and average local ionization energies. J. Mol. Model. 2010, 16, 1679–1691. [Google Scholar] [CrossRef]
  47. Bader, R.F.W. AIM2000 Program, V. 2.0; McMaster University: Hamilton, ON, Canada, 2000. [Google Scholar]
  48. Reed, A.E.; Curtiss, L.A.; Weinhold, F. Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint. Chem. Rev. 1988, 88, 899–926. [Google Scholar] [CrossRef]
  49. ADF2008.01, SCM, Theoretical Chemistry, VrijeUniversiteit, Amsterdam, The Netherlands. Available online: http://www.scm.com (accessed on 20 November 2018).
  50. Su, P.; Li, H. Energy decomposition analysis of covalent bonds and intermolecular interactions. J. Chem. Phys. 2009, 131, 014102. [Google Scholar] [CrossRef] [Green Version]
  51. 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]
  52. Schleyer, P.R.; Maerker, C.; Dransfeld, A.; Jiao, H.J.; van Eikema Hommes, N.J.R. Nucleus-independent chemical shifts: A simple and efficient aromaticity probe. J. Am. Chem. Soc. 1996, 118, 6317–6318. [Google Scholar] [CrossRef]
  53. Zierkiewicz, W.; Michalczyk, M.; Scheiner, S. Implications of monomer deformation for tetrel and pnicogen bonds. Phys. Chem. Chem. Phys. 2018, 20, 8832–8841. [Google Scholar] [CrossRef]
  54. Aronld, W.D.; Oldfield, E. The chemical nature of hydrogen bonding in proteins via NMR: J-couplings, chemical shifts, and AIM theory. J. Am. Chem. Soc. 2000, 122, 12835–12841. [Google Scholar] [CrossRef]
  55. Solel, E.; Kozuch, S. On the power of geometry over tetrel bonds. Molecules 2018, 23, 2742. [Google Scholar] [CrossRef] [PubMed]
Molecules 24 00010 g001 550
Figure 1. MEP maps of 6-OTX3-fulvene on the 0.001 electrons bohr−3 isosurface. Color ranges, in au, are: red, >0.03; yellow, 0.03–0.02; green, 0.02–0.00; blue, <0.
Figure 1. MEP maps of 6-OTX3-fulvene on the 0.001 electrons bohr−3 isosurface. Color ranges, in au, are: red, >0.03; yellow, 0.03–0.02; green, 0.02–0.00; blue, <0.
Molecules 24 00010 g001
Molecules 24 00010 g002 550
Figure 2. Optimized structures of the 6-OTX3-fulvene···NH3 complexes.
Figure 2. Optimized structures of the 6-OTX3-fulvene···NH3 complexes.
Molecules 24 00010 g002
Molecules 24 00010 g003 550
Figure 3. AIM diagrams before and after tetracyano substitution of 6-OTF3-fulvene···NH3 (T = C and Si). The units of distance and Eint are Å and kJ/mol, respectively.
Figure 3. AIM diagrams before and after tetracyano substitution of 6-OTF3-fulvene···NH3 (T = C and Si). The units of distance and Eint are Å and kJ/mol, respectively.
Molecules 24 00010 g003
Molecules 24 00010 g004 550
Figure 4. Plots of deformation densities of the pair-wise orbital interactions (Δρ) in 6-OTX3-fulvene···NH3 (T = Si, Ge; X = H, F) at the PBED3/TZ2P//MP2/aug-cc-pVTZ level. The associated orbital interaction energies are given in kJ/mol. The color code of the charge flow is red→blue, and the isovalue for Δρ(r) is 0.002 a.u.
Figure 4. Plots of deformation densities of the pair-wise orbital interactions (Δρ) in 6-OTX3-fulvene···NH3 (T = Si, Ge; X = H, F) at the PBED3/TZ2P//MP2/aug-cc-pVTZ level. The associated orbital interaction energies are given in kJ/mol. The color code of the charge flow is red→blue, and the isovalue for Δρ(r) is 0.002 a.u.
Molecules 24 00010 g004
Molecules 24 00010 g005 550
Figure 5. Dependence of three attractive energies on the TX3 group (T = C, Si, and Ge; X = H and F).
Figure 5. Dependence of three attractive energies on the TX3 group (T = C, Si, and Ge; X = H and F).
Molecules 24 00010 g005
Table
Table 1. Interaction energy (Eint, kJ/mol), binding energy (Eb, kJ/mol), deformation energy (DE, kJ/mol), intermolecular distance (R, Å), ∠O–T···N (θ, deg), change of O–T bond length (∆r, Å), and the MEP maximum (Vmax, a.u.) of the σ-hole at the O–T bond end.
Table 1. Interaction energy (Eint, kJ/mol), binding energy (Eb, kJ/mol), deformation energy (DE, kJ/mol), intermolecular distance (R, Å), ∠O–T···N (θ, deg), change of O–T bond length (∆r, Å), and the MEP maximum (Vmax, a.u.) of the σ-hole at the O–T bond end.
EintEbDERθ∆rVmax
6-OTX3-fulvene···NH3
–CH3−9.93−9.900.023.228178.30.0060.040
–SiH3−50.13−30.0720.392.353178.00.0510.074
–GeH3−39.30−29.739.562.565180.00.0440.084
–CF3−4.12−3.970.153.498169.30.0030.039
–SiF3−137.07−38.0499.032.052176.90.0630.095
–GeF3−151.31−70.3680.942.074176.70.0460.113
HOTX3···NH3
–CH3−4.76−4.700.063.384174.90.0030.019
–SiH3−18.24−14.893.352.792175.90.0170.051
–GeH3−18.13−15.892.242.872175.50.0180.059
–CF3−3.55−3.420.123.620166.60.0000.029
–SiF3−101.03−14.9686.072.120176.20.0410.078
–GeF3−126.75−49.4977.272.102176.00.0320.097
PhOTX3···NH3
–CH3−7.55−7.430.123.302176.00.0040.028
–SiH3−30.78−20.5610.222.539177.60.0320.057
–GeH3−27.44−21.965.482.697178.10.0300.067
–CF3−3.66−3.570.093.634166.60.0000.024
–SiF3−115.29−22.8392.462.088176.90.0500.080
–GeF3−134.93−−56.9279.172.094175.80.0380.096
Table
Table 2. Electron density (ρ, a.u.), Laplacian (∇2ρ, a.u.), and total energy density (H, a.u.) at the intermolecular BCP.
Table 2. Electron density (ρ, a.u.), Laplacian (∇2ρ, a.u.), and total energy density (H, a.u.) at the intermolecular BCP.
ρ∇2ρH
10.0060.0260.002
20.0330.079−0.008
30.0270.082−0.008
40.0040.0190.001
50.0630.251−0.020
60.0830.230−0.033
Table
Table 3. Charge transfer (CT, e) and second-order perturbation energies (E(2), kJ/mol) in the complexes.
Table 3. Charge transfer (CT, e) and second-order perturbation energies (E(2), kJ/mol) in the complexes.
CTOrbitalsE(2)OrbitalsE(2)
10.003Lp(N)→σ*C-O5.23
20.093Lp(N)→p*Si186.72Lp(N)→σ*Si-H60.69
30.062Lp(N)→p*Ge129.91Lp(N)→σ*Ge-H34.44
4−0.002Lp(F)→σ*N-H0.67
50.163Lp(N)→σ*Si-O102.66Lp(N)→σ*Si-F296.53
60.186Lp(N)→σ*Ge-O134.76Lp(N)→σ*Ge-F401.61
Note: CT is the sum of the charge on all atoms of NH3.
Table
Table 4. NICS(1)zz in 6-OTX3-fulvene···NH3 and 6-OTX3-fulvene (in parentheses), as well as their difference (Δ); all are in ppm.
Table 4. NICS(1)zz in 6-OTX3-fulvene···NH3 and 6-OTX3-fulvene (in parentheses), as well as their difference (Δ); all are in ppm.
NICS(1)zzΔ
1−6.0781 (−5.8122)−0.2659
2−6.7149 (−5.7074)−1.0125
3−6.9462 (−6.0684)−0.8778
4−5.2846 (−4.7814)−0.5032
5−6.0377 (−5.7074)−0.3303
6−6.2645 (−4.9984)−1.2661

Share and Cite

MDPI and ACS Style

Hou, M.-C.; Yang, S.-B.; Li, Q.-Z.; Cheng, J.-B.; Li, H.-B.; Liu, S.-F. Tetrel Bond between 6-OTX3-Fulvene and NH3: Substituents and Aromaticity. Molecules 2019, 24, 10. https://doi.org/10.3390/molecules24010010

AMA Style

Hou M-C, Yang S-B, Li Q-Z, Cheng J-B, Li H-B, Liu S-F. Tetrel Bond between 6-OTX3-Fulvene and NH3: Substituents and Aromaticity. Molecules. 2019; 24(1):10. https://doi.org/10.3390/molecules24010010

Chicago/Turabian Style

Hou, Ming-Chang, Shu-Bin Yang, Qing-Zhong Li, Jian-Bo Cheng, Hai-Bei Li, and Shu-Feng Liu. 2019. "Tetrel Bond between 6-OTX3-Fulvene and NH3: Substituents and Aromaticity" Molecules 24, no. 1: 10. https://doi.org/10.3390/molecules24010010

APA Style

Hou, M. -C., Yang, S. -B., Li, Q. -Z., Cheng, J. -B., Li, H. -B., & Liu, S. -F. (2019). Tetrel Bond between 6-OTX3-Fulvene and NH3: Substituents and Aromaticity. Molecules, 24(1), 10. https://doi.org/10.3390/molecules24010010

Article Metrics

Back to TopTop