3.1 ZSM-5(4T)@MIL-53(Al) core-shell catalyst
From the above model (Fig. 1(c)), it can be observed that Cat. had 5 active sites (marked by purple circles) locally. In order to avoid unnecessary repetition of similar reaction mechanism of each active site, this research mainly took active site 1 as an example to discuss the catalytic mechanism of disproportionation between M1 and M3 in order to prepare M2. The reaction equations of Cat. catalytic system was shown in Fig. 2(a), the corresponding reaction process and its key atomic numbers were shown in Fig. 2(b).
The disproportionation reaction was divided into two channels. In the first channel, Cat. and (CH3)3SiCl adsorbed each other and reacted assigning to intermolecular force, and the intermediate I1 and CH4 were generated through the transition state TS1. Then, I1 continued to adsorbed with CH3SiCl3, when the intermediate II and (CH3)2SiCl2 were formed through the transition state TS2. Finally, the intermediate II adsorbed with CH4 though the transition state TS3, the reduced Cat. and (CH3)2SiCl2 were formed. Thus, the production of main product (CH3)2SiCl2 indicated that the first channel was the main reaction route.
In the second channel, Cat. and CH3SiCl3 adsorbed each other and reacted assigning to intermolecular force, and the intermediate III and CH4 were generated through the transition state TS4. Then, III continued to adsorb with CH3SiCl3, and the intermediate IV and SiCl4 were formed through the transition state TS5. Finally, the intermediate IV adsorbed with CH4 though the transition state TS6, the reduced Cat. and Si(CH3)4 were formed. Thus, the by-product SiCl4 and Si(CH3)4 indicated that the second channel was the side reaction route.
Figure 2 (a) The reaction equations of ZSM-5(4T)@MIL-53(Al) catalytic system. (b) The reaction process catalyzed by ZSM-5(4T)@MIL-53(Al) core-shell catalyst (No. 1 active site) and its key atomic numbers. Rx, TSx, Ix, and Px (x = 1–6) represent the reactant, the transition state, the intermediate, and the product, respectively. Cat. represents ZSM-5(4T)@MIL-53(Al) core-shell catalyst
Based on the reaction mechanism described in Fig. 2, geometric optimized calculations were performed for all reactants, products and transition states involved in the disproportionation reaction on all active sites of Cat. The virtual vibrational modes of transition states were obtained by vibrational frequency analysis. Furthermore, the structures of the transition states TSx(x = 1–6) were verified. The key parameters and virtual vibrational modes of transition states in Cat. system were shown in Fig. 3 (active site 1) and Fig.S1-S4 (active sites 2–5) in the supplementary material. Obviously, the atomic vibrations of the transition states in each step mainly came from relative atomic stretching vibrations such as O1, H1, Si1, Si2, Cl1, Cl2, C1, and C2, and each transition state had a simultaneous reaction to its corresponding trends in the vibrations of reactants and products. In detail, the key atomic vibrational directions of the magenta arrows tend to move towards the product, while those of the purple arrow tend to move towards the reactant. Combined with Fig. 2 and Fig. 3, take the main reaction step TS1 of disproportionation reaction catalyzed by Cat. as an example: during the reaction between Cat. with (CH3)3SiCl, the relatively low energy Si1-C1 bond of (CH3)3SiCl and O1-H1 bond of Cat. and (CH3)3SiCl broke apart and move away from each other, while the Brønsted acidic H1 atom got closer and bonded with C1 atom of the methyl group to form CH4. The O1 atom of catalyst and Si1 atom groups were bonded to form intermediate I. Those were the product P1. Observing the reverse process of this step, starting from the transition state TS1, the Si1 atom and C1 atom in the methyl group gradually approached and bonded to form (CH3)3SiCl, and the H1 atom combined with O1 atom to bond and form Cat., those were the corresponding reactant R1. The analysis of other transition states was similar to TS1, and the results showed that the frequency vibrational directions and trends of all key atoms of the transition states were consistent with the aforementioned reaction mechanism. In summary, the reliability of the transition state can verify the credibility of the reaction mechanism.
Figure 3 The vibrational modes of the transition states of active site 1, its key atomic number and distance of each reactant, product and transition state
In order to further test the reliability of the transition state structure of each step, the intrinsic response coordinates (IRC) method, which could track both the forward and reverse directions of the transition state was used. Thus, IRC calculation was carried out for the transition state structure of each step. Figure 4 showed the variational tendency of the spacing of key atoms on active sites 1, 2 and 3 of Cat. along IRC. In addition, relative IRC analyses of active sites 4 and 5 were shown in Fig.S5 of supplementary materials. The capital letter S in English stood for the measure of IRC, and the positive and negative value ranges of S represented the regions of product and reactant, respectively. When S pointed to 0, it meant that the structure tended to transition state. It can be seen from the figure that the disproportionation reaction started from the reactants and passed through the transition state to generate the product, and the distance between the key atoms showed a trend of lengthening or shortening. It was consistent with the vibrational trends of the transition state in Fig. 3, which proved the correctness of the reaction mechanism.
Figure 4 Variation trend of key atom spacing along IRC in ZSM-5(4T)@MIL-53(Al) catalytic system; (a)-(f) represent the IRC of transition states TS1-TS6, respectively
The energy calculation of the disproportionated structure was carried out by using the M06-2X/Def2-TZVP hybrid functional method of the Minnesota series, and ZPE correction was used at the same time. Figure 5 and Table 1 are the reaction energy diagram of Cat (active sites 1–5) disproportionate synthesis on M2 and the numerical statistics table of the related activation energy. Taking active site 1 (indicated by the black line) as an example, the activation energy barrier data (1Ea1-1Ea6) of each step of the disproportionation reaction were: 122.76, 21.89, 152.14, 76.02, 20.32 and 112.55 kJ•mol− 1, respectively. From the above numerical comparison, it can be seen that the largest activation energy barrier in the main reaction was the third step and its activation energy of the rate-determining step was 152.14 kJ•mol− 1. Similarly, the activation energy of the rate-determining step of the side reaction was 112.55 kJ•mol− 1. The rate-determining step energy barrier of the main reaction was larger than that of the side reaction, indicating that the active site 1 of Cat. mainly produces by-products. Other active sites (2–5) were analyzed in the same way as active site 1. From the above data, it can be seen that the activation energy barrier data of the active sites 1, 2, 3, 4, and 5 in the main reaction rate-determining step were 152.14, 118.06, 152.74, 162.69 and 185.31 kJ•mol− 1 (marked in red in the table), the activation energy barrier data for the rate-determining step of the side reaction were 112.55, 122.85, 131.34, 162.47 and 157.62 kJ•mol− 1 (in the table marked in blue). The rate-determining step activation energy of the main reaction at sites 1, 3, 4 and 5 were all larger than those of the side reactions, indicating that these sites had little competitive advantage in synthesizing the target product (CH3)2SiCl2. However, the activation energy of the main reaction rate-determining step of the active site No. 2 was lower than that of the side reaction, and the active site No. 2 had a lower energy barrier than other sites. It was obvious that active site 2 was more conducive to the synthesis of M2 and showed a better disproportionation effect.
Table 1
Activation energy data of all active sites (1–5) of ZSM-5(4T)@MIL-53(Al) core-shell catalytic system and AlCl3/ZSM-5@MIL-53(Al) catalytic system (kJ•mol− 1)
Catalyst and Active Center | Cat. | Cat.& |
No.1 | No.2 | No.3 | No.4 | No.5 |
Activation energy (kJ•mol− 1) | Ea1 | 122.76 | 118.06 | 152.74 | 140.04 | 185.31 | 61.69 |
Ea2 | 21.89 | 20.37 | 30.66 | 33.53 | 100.74 | 86.92 |
Ea3 | 152.14 | 96.03 | 29.02 | 162.69 | 124.26 | 137.94 |
Ea4 | 76.02 | 82.28 | 12.65 | 162.47 | 64.79 | 41.66 |
Ea5 | 20.32 | 18.21 | 51.47 | 30.15 | 157.62 | - |
Ea6 | 112.55 | 122.85 | 131.34 | 136.51 | 101.64 | - |
Figure 5 Active sites (1–5) of ZSM-5(4T)@MIL-53(Al) core-shell catalyst and activation energy of AlCl3/ZSM-5@MIL-53(Al) catalytic system
Table 1 Activation energy data of all active sites (1–5) of ZSM-5(4T)@MIL-53(Al) core-shell catalytic system and AlCl3/ZSM-5@MIL-53(Al) catalytic system (kJ•mol− 1)
3.2 AlCl3/ZSM–5(4T)@MIL-53(Al) catalyst
The reaction equations of Cat.& catalytic system were shown in Fig. 7(a). The disproportionation reaction was divided into two channels: main reaction and side reaction. In the first channel, Cat.& and (CH3)3SiCl adsorbed each other and reacted assigning to intermolecular force, and the intermediate I& and (CH3)2SiCl2 were generated through the transition state TS1. Then, I& continued to adsorb with CH3SiCl3 through the transition state TS2. The reduced Cat.& and (CH3)2SiCl2 were formed. Thus, the main product (CH3)2SiCl2 indicated that the first channel was the main reaction route. In the second channel, Cat.& and CH3SiCl3 adsorbed each other and reacted assigning to intermolecular force, and the intermediate I& and CH4 were generated through the transition state TS3. Finally, the intermediate I& adsorbed with (CH3)3SiCl through the transition state TS4, the reduced Cat.& and Si(CH3)4 were formed. Thus, the by-product SiCl4 and Si(CH3)4 indicated that the second channel was the side reaction route.
Figure 7 (a) Reaction equations of AlCl3/ZSM-5(4T)@MIL-53(Al) catalytic system. (b) Reaction process catalyzed by AlCl3/ZSM-5(4T)@MIL-53(Al) catalyst and its key atomic numbers. Rx, TSx, Ix, and Px(x = 1–4) represent reactants, transition states, intermediates, and products, respectively. Cat.& represents AlCl3/ZSM-5(4T)@MIL-53(Al) core-shell catalyst
Based on the reaction flow diagram and its key atomic numbers of Cat.& for disproportionation synthesis on M2 as shown in Fig. 7(b), the virtual frequency vibrational analysis of the disproportionation reaction was carried out. Figure 8 showed the imaginary frequency vibration of the transition state of the reaction channel and its key atomic numbers. It can be seen from the figure that the atomic vibrations of the transition states of the disproportionation steps mainly came from the relative atomic stretching vibrations such as Si1, Si2, Cl1, Cl2, C1, C2, Al& and Cl&. Taking TS1 as an example: Cat.& and (CH3)3SiCl approached and interacted with each other, the Si1-C1 bond of (CH3)3SiCl and Al&-Cl& of the catalyst were broken respectively, and the C1 atom of methyl group and Al& atom were close to each other and formed a bond, resulting in the forming of the intermediate I&. The Si1 atom of (CH3)3SiCl and Cl& atom were close to each other and formed a bond to generate the product (CH3)2SiCl2. Those were the product P1. Observing the reverse process of this reaction and another vibrational direction of TS1, Si1 and C1 atoms were close to each other and form a bond, and Al& and Cl& atoms were close to each other and formed a bond. Those were the reactant R1. The vibrations of other transition states were also analyzed similarly. It was shown from the results that the distance between the key atoms increased or decreased with the reaction progresses, indicating that the reaction path was correct and each transition state was in the force constant matrix. There was positive frequency for all stable structure (reactants, products and intermediates) or only one imaginary frequency for each transition state, and the frequency vibrational direction and trend of key atoms were unified with the aforementioned reaction mechanism. Not only those above showed that the transition state structure was real and reliable, but also indicated that the reaction mechanism was credible.
Figure 8 The vibrational modes of the transition states of AlCl3/ZSM-5(4T)@MIL-53(Al) catalyst, its key atomic numbers and distances of each reactant, product and transition state
In order to further test the reliability of the transition state structure of each step, the IRC method, which can track both the forward and reverse directions of the transition state was used. Thus, IRC calculation was carried out for the transition state structure of each step. The change trend of the spacing of important atoms of Cat.& along IRC was shown in Fig. 9. It can be seen from the figure that disproportionation reaction started from the reactant and passed through the transition state to produce products, and its distance between key atoms tended to lengthen or shorten. That was consistent with the above virtual vibrational analysis and reaction mechanism, which proved the correctness of the reaction mechanism.
Figure 9 The variational trend of key atomic spacing along IRC in modified AlCl3/ZSM-5(4T)@MIL-53(Al) catalytic system
The energy calculation of the structure involved in the Cat.& catalytic system was carried out by using the M06-2X/Def2-TZVP hybrid functional method of the Minnesota series, and ZPE correction was used at the same time. Combined with Fig. 5 (the brown line) and Table 1, we can saw that the rate-determining step of the main reaction occurs in the second step of disproportionation, and its activation energy barrier datum was 86.92 kJ•mol− 1. That was to say, the activation energy barrier of the rate-determining step in the main reaction was 86.92 kJ•mol− 1. Similarly, the activation energy barrier of the rate-determining step in the side reaction was 137.94 kJ•mol− 1. It was obvious that the activation energy of the rate-determining step in the side reaction was higher than that of the main reaction, indicating that the Cat.& was beneficial to the disproportionation synthesis on M2. Compared with the activation energy data before loading (Fig. 6), the activation energy barrier value (86.92 kJ•mol− 1) of the rate-determining step in the main reaction of Cat.& catalytic system rate-determining step was the smallest. Moreover, the activation energy barrier difference value between the main reaction and side reaction was quite large, indicating that the main reaction of the disproportionation reaction involving the Cat.& was more competitive, which made it easier to synthesize the target product (CH3)2SiCl2 through a lower energy barrier and enhanced the difficulty in producing by-products at the same time.
In order to further study the activity of Cat.&, the catalyst was analyzed in combination with bond order, LOL and ELF. After calculation by using the Multiwfn software, the MBO bond order of the Al&-Cl& bond of Cat.& was 0.7114; the LBO bond order was 0.1984, and the Mulliken bond order was 0.3313. Figure 10 showed the ELF and LOL diagrams of Cat.&. As shown in Fig. 10(a) (ELF), the orange-yellow inner core of the Cl& atom indicated that its strong localization. The color gradually changed outward a very narrow blue, indicating that the orbital electrons of Cl& atom inner layer the orbital could be shared. However, the outer orbital of the entire Cl& atom was occupied by broad band orange and narrow band yellow, indicating that the valence orbital electrons of the outer layer were well localized. In addition, its electron-withdrawing characteristics were in line with Lewis acidic characters. Therefore, Cl& atom could attract the electrons provided by the methyl group on the silane and reacted each other. The orange inner area of Al& atom indicated its robust electron localization. Besides, the outer layer hue rapidly transformed to a broadband dark blue, indicating that its strong delocalization and electronic sharing with other atoms. From the region of Al& and Cl& atoms, it could be observed that the Al&-Cl& bond exhibited obvious polar and ionic bond characteristics.
LOL analysis in Fig. 10(b) was similar to the said ELF analysis. It could be seen from the figure that the lighting yellow center of the Cl& atomic nucleus showed well localization. The color gradually changed outward a very narrow blue, indicating that the orbital electrons of Cl& atom inner layer the orbital could be shared. However, the outer orbital of the entire Cl& atom was occupied by broad band orange and narrow band yellow, indicating that the valence orbital electrons of the outer layer were well localized. In addition, its electron-withdrawing characteristics were in line with Lewis acidic characters. Therefore, Cl& atom could attract the electrons provided by the methyl group on the silane and reacted each other. The light green inner area of Al& atom indicated its weak electron localization. Besides, the outer layer hue rapidly transformed to a broadband dark blue, indicating that its strong delocalization and electronic sharing with other atoms. In summary, the catalytic activity of Cat.& was further improved after loading AlCl3. This can be well matched with the previous activation energy data and mechanism analysis.
Figure 10 ELF(a) and LOL(b) analyses of AlCl3/ZSM-5(4T)@MIL-53(Al) catalyst