Skip to content
BY-NC-ND 3.0 license Open Access Published by De Gruyter November 17, 2018

Theoretical study of the formation of a spiro-Sn-heterocyclic compound by cycloaddition reaction of Me2C=Sn: and ethene

  • Xiaojun Tan EMAIL logo and Xiuhui Lu EMAIL logo

Abstract

X2C=Sn: compounds (X=H, Me, F, Cl, Br, Ph, Ar) are new species. The cycloaddition reactions of X2C=Sn: are also a new study field of unsaturated stannylene chemistry. The mechanism of cycloaddition reaction between singlet Me2C=Sn: and ethene was investigated for the first time using the MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz) method. From the potential energy profile, it was predicted that the reaction has one dominant channel in which the 5p unoccupied orbital of Sn: in Me2C=Sn: and the π orbital of ethene form a π→p donor-acceptor bond in an intermediate product. Instability of the intermediate product results in its isomerization to a four-membered ring of stannylene. The four-membered stannylene further combines with ethene to form another intermediate product that further isomerizes to a spiro-Sn-heterocyclic ring compound.

Introduction

In 1997, Clouthier and co-workers [1] analyzed for the first time unsaturated germylidene (H2C=Ge:) derivatives. Subsequently, ab initio calculations predicted the molecular structure [2], [3], electronic spectrum [2], oscillatory fluorescence decay [2], the ground state structure [3] and stimulated emission pumping (SEP) spectra [4] of the first excited singlet state of germylidene or derivatives. Stogner and Grev have published many papers about ab initio calculations [5] on both germylidene and the trans-bent germyne HC≡GeH isomer. We have also been involved in this subject [6], [7], [8], [9]. As the studies progressed, research on the analogous species X2C=Sn: (X=H, Me, F, Cl, Br, Ph, Ar) has emerged. The stannylene X2C=Sn: is a new chemical species [10]. In this report, we have explored theoretically the cycloaddition reactions of X2C=Sn: with ethene.

Results and discussion

We used the method of second-order perturbation theory (MP2) and Gaussian 09 package to optimize the structure of Me2C=Sn:, and the parameters of cycloaddition reaction with ethene including the transition states and the intermediate products at the MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz) [11] theory level. In order to further confirm the correctness of the relevant species, a vibration analysis was included. Finally, the intrinsic reaction coordinate (IRC) [12], [13] was also calculated for all transition states to determine the reaction pathways.

The energies of the singlet and triplet states of Me2C=Sn: (R1) calculated by the MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz) method are −120.81220 and −120.77967 a.u, respectively. Accordingly, the ground state of Me2C=Sn: is a singlet state. Its cycloaddition reaction with ethene (R2) may have two possible outcomes. The geometrical parameters of the intermediate

products INT1 and INT2, transition states TS1, TS1.1, TS1.2 and TS2 and products P1, P1.1, P1.2 and P2 which appear in the schematic reactions are given in Figure 1, the energies are listed in Table 1 and the entropy, enthalpy and Gibbs free energy are listed in Table 2. The potential energy profile of the reactions is shown in Figure 2.

Figure 1 Optimized MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz) geometrical parameters and the atomic numbering for the species in cycloaddition reaction between Me2C=Sn: and ethene. Bond lengths are in (Å) and angles in (°).
Figure 1

Optimized MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz) geometrical parameters and the atomic numbering for the species in cycloaddition reaction between Me2C=Sn: and ethene. Bond lengths are in (Å) and angles in (°).

Table 1

The electronic structure energy (Eese, a.u.) and relative energies (ER, kJ/mol) for the species from MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz) method at 298 K and 101 325 Pa.

ReactionSpeciesMP2/GENECP
EeseER
aReaction (1)R1+R2−199.158720.0
INT1−199.17551−44.1
TS1 (INT1-P1)−199.16128−6.7
P1−199.22277−168.2
TS1.1 (P1-P1.1)−199.1326068.6
P1.1−199.16449−15.1
TS1.2 (P1-P1.2)−199.11668110.4
P1.2−199.1544611.2
bReaction (2)P1+R2−277.569300.0
INT2−277.58145−31.9
TS2 (INT2-P2)−277.5618719.5
P2−277.565948.8
  1. aER=EeseEese, (R1+R2), bER=Eese—Eese(P1+R2).

Table 2

Entropy (S, a.u.), enthalpy (H, a.u.) and Gibbs free energy (G, a.u.) for the species from MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz) methods at 298 K and 101 325 Pa.

ReactionSpeciesH, a.u.S, a.u.G, a.u.
Reaction (1)R1+R2−199.017032.09768×10−4−199.07675
INT1−199.031221.52654×10−4−199.07959
TS1 (INT1-P1)−199.018261.38714×10−4−199.05964
P1−199.077151.39870×10−4−199.11887
TS1.1 (P1-P1.1)−198.991801.42594×10−4−199.03433
P1.1−199.021121.50390×10−4−199.06597
TS1.2 (P1-P1.2)−198.977731.39018×10−4−199.01920
P1.2−199.013751.41544×10−4−199.05597
Reaction (2)P1+R2−277.368832.23496×10−4−277.42993
INT2−277.378291.73142×10−4−277.43549
TS2 (INT2-P2)−277.360031.54895×10−4−277.40622
P2−277.363101.61220×10−4−277.41119
Figure 2 The potential energy surface for the cycloaddition reaction between Me2C=Sn: and ethene with MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz).
Figure 2

The potential energy surface for the cycloaddition reaction between Me2C=Sn: and ethene with MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz).

The unique imaginary frequencies of the transition states TS1, TS1.1, TS1.2 and TS2 by vibrational analysis are 214.5i, 169.1i, 828.4i and 102.1i, respectively. Therefore, these transition states can be affirmed as the real ones. IRC (with the step size of 0.1 amu−1/2 bohr) analysis confirms that TS1 connects INT1 and P1, TS1.1 connects P1 and P1.1, TS1.2 connects P1 and P1.2 and TS2 connects INT2 and P2. According to Figure 2, the reaction (1) consists of four steps: in the first step, the two reactants R1 and R2 form an intermediate product INT1. According to Figure 2 and Table 2, the reaction is a barrier-free exothermic process; the molar constant volume heat of reaction (ΔrUm) and molar heat of reaction (ΔrHm) at normal temperature and pressure are −44.1 and −37.3 kJ/mol, and the molar Gibbs free energy of reaction (ΔrGm) is −7.5 kJ/mol. In the second step, the intermediate product INT1 undergoes isomerization to a four-membered ring of stannylene P1 via a transition state TS1 with an energy barrier of 37.4 kJ/mol. According to Figure 2 and Table 2, the reaction is exothermic, and the ΔrUm and ΔrHm values at normal temperature and pressure are 124.1 and 120.6 kJ/mol, and the ΔrGm value is −110.6 kJ/mol. In the third and fourth steps, the product P1 undergoes Me-transfer and H-transfer via transition states TS1.1 and TS1.2 with energy barriers of 236.8 and 278.6 kJ/mol, respectively. The products of Me-transfer and H-transfer are P1.1 and P1.2. According to Figure 2 and Table 2, both reactions are endothermic, and the ΔrUm and ΔrHm values at normal temperature and pressure are 153.1, 179.4 and 147.1, 166.5 kJ/mol, respectively. The ΔrGm values are 138.9 and 165.1 kJ/mol, respectively. Accordingly, the transitions P1→P1.1 and P1→P1.2 are thermodynamically prohibited at normal temperature and pressure, and the product of reaction (1) will be P1.

In reaction (2), the four-membered ring stannylene P1 undergoes a further reaction with ethene R2 to form a spiro-Sn-heterocyclic compound P2. According to Figure 2, in the reaction (2), the product P1 formed in reaction (1) undergoes a subsequent reaction with ethene R2 to form an intermediate INT2. According to Figure 2 and Table 2, the reaction is a barrier-free exothermic process with the ΔrUm and ΔrHm values at normal temperature and pressure of 31.9 and 24.8 kJ/mol, and the ΔrGm value is −14.6 kJ/mol. Then, the intermediate INT2 undergoes isomerization to a spiro-Sn-heterocyclic compound P2 via a transition state TS2 with an energy barrier of 51.4 kJ/mol. According to Figure 2 and Table 2, this reaction is an endothermic process with the ΔrUm and ΔrHm values at normal temperature and pressure of 40.7 and 39.9 kJ/mol, and the ΔrGm value is 63.8 kJ/mol. According to Figure 2, reaction (2) and the transformations P1→P1.1 and P1→P1.2 in reaction (1) are competitive processes. As the transitions P1→P1.1 and P1→P1.2 are thermodynamically forbidden, reaction (2) is the dominant reaction pathway. Considering that the ΔrGm value of INT2→P2 is 63.8 kJ/mol and the ΔrGm value of P1+R2→P2 is 49.2 kJ/mol, the transition P1+R2→INT2→P2 is a continuous reaction. Analysis of the following thermodynamic function shows that if the reaction is carried out under liquid phase conditions,

ΔG(p2)ΔG(p1)=p1p2ΔVdp

ΔV≈−1 and P2≈ΔG(P1)+P1=165 125 Pa. Accordingly, at the temperature of 298 K in liquid phase, for the reaction P1+R2→P2, the pressure of the reaction system must be greater than 165 125 Pa (1.6 atm).

According to all analyses, reaction (2) should be the dominant pathway of the cycloaddition between singlet Me2C=Sn: and ethene. It can be formulated as follows:

R1+R2INT1TS1P1+R2INT2TS2P2

The frontier molecular orbitals of R1, R2 and P1 are shown in Figure 3. As can be seen, the mechanism of the transformation can be explained by analysis of the frontier molecular orbital diagrams (Figures 4 and 5). Briefly (Figures 1 and 4), when Me2C=Sn: (R1) initially interacts with ethene (R2), the 5p unoccupied orbital of the Sn atom in Me2C=Sn: (R1) undergoes insertion of the π orbital of ethene to form a π→p donor-acceptor bond, leading to the formation of the intermediate product INT1. As the reaction progresses, INT1 isomerizes to a four-membered stannylene P1 via the transition state TS1. Because P1 is still an active molecule, it may further react with ethene to form a spiro-Sn-heterocyclic compound P2 through the intermediary of INT2 and TS2. The mechanism of this reaction can be explained by analysis of Figures 1 and 5. Before the transition state TS2, the covalent bonds are formed between Sn and C(6), and Sn and C(7). After the transition state TS2, the Sn atom hybridization changes to a sp3 hybrid orbital. Overall, the intermediate product INT2 undergoes isomerization to a spiro-Sn-heterocyclic compound P2 via the transition state TS2.

Figure 3 The frontier molecular orbitals of R1, R2 and P1.
Figure 3

The frontier molecular orbitals of R1, R2 and P1.

Figure 4 A schematic interaction diagram for the frontier orbitals of Me2C=Sn: (R1) and C2H4 (R2).
Figure 4

A schematic interaction diagram for the frontier orbitals of Me2C=Sn: (R1) and C2H4 (R2).

Figure 5 A schematic diagram for the frontier orbitals of P1 and C2H4 (R2).
Figure 5

A schematic diagram for the frontier orbitals of P1 and C2H4 (R2).

Conclusion

On the basis of the potential energy profile, the cycloaddition reaction between singlet species Me2C=Sn: and ethene obtained with the MP2/GENECP (C, H in 6-311++G**; Sn in LanL2dz) method was predicted. This reaction has one dominant channel that consists of four steps: (1) the two reactants first form an intermediate INT1 through a barrier-free exothermic reaction of 44.1 kJ/mol; (2) the intermediate INT1 isomerizes to a four-membered stannylene P1 via the transition state TS1 with an energy barrier of 37.4 kJ/mol; (3) the four-membered stannylene P1 further reacts with ethene R2 to form another intermediate product INT2 through a barrier-free exothermic reaction of 31.9 kJ/mol; and (4) the intermediate product INT2 undergoes isomerization to a spiro-Sn-heterocyclic compound P2 via the transition state TS2 with an energy barrier of 51.4 kJ/mol. At the temperature of 298 K in liquid phase, the pressure of the reaction system needs to be greater than 165 125 Pa (1.6 atm).

References

[1] Harper, E. W. H.; Ferrall, A.; Hilliard, R. K.; Stogner, S. M.; Grev, R. S.; Clouthier, D. J. Laser spectroscopic detection of the simplest unsaturated silylene and germylene. J. Am. Chem. Soc.1997, 119, 8361–8362.10.1021/ja9716012Search in Google Scholar

[2] Hostutler, D. A.; Smith, T. C.; Li, H. Y.; Clouthier, D. J. The electronic spectrum, molecular structure, and oscillatory fluorescence decay of jet-cooled germylidene (H2C=Ge:), the simplest unsaturated germylene. J. Chem. Phys.1999, 111, 950–958.10.1063/1.479187Search in Google Scholar

[3] Hostutler, D. A.; Clouthier, D. J.; Pauls, S. W. The ground state of germylidene (H2C=Ge:). J. Chem. Phys.2002, 116, 1417–1423.10.1063/1.1431274Search in Google Scholar

[4] He, S. G.; Tackett, B. S.; Clouthier, D. J. A stimulated emission pumping study of the first excited singlet state of germylidene (H2C=Ge:). J. Chem. Phys.2004, 121, 257–264.10.1063/1.1758699Search in Google Scholar PubMed

[5] Stogner, S. M.; Grev, R. S. Germyne, HC≡GeH, and the excited states of 1-germavinylidene H2C=Ge:. J. Chem. Phys.1998, 108, 5458–5464.10.1063/1.475934Search in Google Scholar

[6] Tian, C. L., Xu, Y. H., Lu, X. H. Theoretical studies on the mechanism of formation of a heteropolycyclic germanic compound between dimethyl-germylidene and ethylene. Int. J. Quant. Chem.2010, 110, 1675–1681.10.1002/qua.22337Search in Google Scholar

[7] Lu, X. H., Xu, Y. H., Yu, H. B., Wu, W. R. Ab initio study on the mechanism of forming a germanic hetero-polycyclic compound between alkylidenegermylene and ethylene. J. Phys. Chem. A2005, 109, 6970–6973.10.1021/jp0515075Search in Google Scholar PubMed

[8] Lu, X. H., Lian, Z. X., Li, Y. Q., Wang, Z. N. Ab initio study of mechanism of forming bis-heterocyclic compound with Si and Ge between dimethylsilylene germylidene (Me2Si=Ge:) and ethane. Struct. Chem.2012, 23, 809–813.10.1007/s11224-011-9929-2Search in Google Scholar

[9] Lu, X. H., Xu, Y. H., Shi, L. Y., Han, J. F., Lian, Z. X. Theoretical study on the mechanism of cycloaddition reaction between dimethyl germylidene and formaldehyde. J. Organomet. Chem.2009, 694, 4062–4068.10.1016/j.jorganchem.2009.08.023Search in Google Scholar

[10] Bundhun, A., Momeni, M. R., Shakib, F. A., Ramasami, P., Gaspar, P. P., Schaefer III, H. F. Toward unsaturated stannylenes Y2Z=Sn: and related compounds with triplet electronic ground states. RSC. Adv.2016, 6, 53749–53759.10.1039/C6RA00492JSearch in Google Scholar

[11] Curtis, L. A.; Raghavachari, K.; Pople, J. A. Gaussian-2 theory using reduced Møller-Plesset orders. J. Chem. Phys.1993, 98, 1293–1298.10.1063/1.464297Search in Google Scholar

[12] Fukui, K. A. Formulation of the reaction coordinate. J. Phys. Chem.1970, 74, 4161–4163.10.1142/9789812795847_0038Search in Google Scholar

[13] Ishida, K.; Morokuma, K.; Komornicki, A. The intrinsic reaction coordinates. An ab initio calculation for HNC to HCN and H+CH4 to CH4+H. J. Chem. Phys.1981, 66, 2153–2156.10.1063/1.434152Search in Google Scholar

Received: 2018-07-21
Accepted: 2018-09-10
Published Online: 2018-11-17
Published in Print: 2018-12-19

©2018 Walter de Gruyter GmbH, Berlin/Boston

This article is distributed under the terms of the Creative Commons Attribution Non-Commercial License, which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Downloaded on 24.4.2024 from https://www.degruyter.com/document/doi/10.1515/hc-2018-0129/html
Scroll to top button