基于ReaxFF的甲烷无氧转化气相机理研究
English
Gas-Phase Mechanism Study of Methane Nonoxidative Conversion by ReaxFF Method
-
1. Introduction
With the constant depletion of oil resources, using natural gas/shale gas to produce chemical raw materials has been receiving more and more attention. As the main component of natural gas and shale gas, the activation and conversion of methane is a challenging task due to its inherent inertia, such as zero dipole moment 1, low electron affinity 2, and high C-H bond energy (439 kJ∙mol−1) 3.
A number of methods have been developed to activate and convert methane, which are mainly classified into indirect routes and direct routes 4-12. Indirect conversion of methane involves first converting methane into syngas and then converting syngas into other products 13-17, the efficiency of which is relatively low because of the complexity of the transformation process. Direct conversion of methane includes oxidative and nonoxidative pathways. Oxidative conversion of methane such as selective oxidation of methane (SOM) and oxidative coupling of methane (OCM) cannot avoid the overoxidation of methane to carbon dioxide and thus leads to low carbon-atom efficiency and CO2 emission 18-29. In contrast, the nonoxidative conversion of methane because of short process, high carbon-atom efficiency and less CO2 emission is considered to be a more promising way 30-47.
Since 1993, Wang et al. have firstly reported that Mo doped ZSM-5 catalysts can transform methane to aromatic products 30, but the catalyst was easily deactivated due to carbon deposition. Until 2014, a breakthrough work was reported by Guo and colleagues that methane can be transformed to olefins, aromatics, and hydrogen (MTOAH) under nonoxidative conditions 48. The conversion of methane reached at 48.1% and the selectivity of ethylene in the products peaked at 48.4% at 1363 K. Importantly, no coke deposition and no catalyst deactivation were observed during a 60-hour test. However, the mechanism of MTOAH is still under discussion. Through preliminary theoretical calculations and VUV-SPI-MBMS experimental analysis 48, it was proposed that methane first dissociated at the Fe1©SiC2 active center with the methyl and hydrogen atom adsorbed at Fe and C site, respectively. Then the methyl released to the gas phase and collided with each other to form the final product (ethylene and aromatics). In order to further demonstrate the possibility of gas-phase mechanism, in 2019, Hao et al. reported that methane conversion to olefins and aromatics over an iron-coated catalytic quartz reactor could be enhanced by hydrogen radicals which were produced by 1, 2, 3, 4-tetrahydronaphthalene (THN) and benzene 49. Kim's group studied the mechanism of nonoxidative methane coupling on a single-atom iron catalyst by density functional theory (DFT) and microkinetic analyses 50. They concluded that ethylene was formed by joint surface and gas-phase mechanisms, where ethane molecules generated in the gas phase were activated by surface active center to form ethylene molecules. Recently, our group also studied the mechanism of this system by utilizing ab initio molecular dynamics (AIMD) simulations and DFT methods 51, we found that ethylene can be easily formed on the surface active center through the quasi-MvK mechanism. From these mechanism studies, the main argument is that whether the ethylene is formed in the gas phase or on the catalyst surface during MTOAH reaction.
In this work, the reactive force field (ReaxFF) simulation method was used to study the gas-phase mechanism of methane nonoxidative conversion under near experimental conditions. It is found that ethylene can be hardly formed under sole gas-phase reaction conditions. The main simulated product via gas-phase mechanism is ethane. The C2H4 can be produced by thermal cracking of THN molecules at high temperature, which may influence the conclusion of experiment that THN could promote the conversion of methane to olefins and aromatics.
2. Models and computational methods
Five different models were built to study the gas-phase reaction mechanism based on the experimental speculation. (A) CH4/CH3∙/H2 mixture with the ratio of 2 : 2 : 1 (i.e., 80 methane molecules, 80 methyl radicals and 40 hydrogen molecules); (B) 200 methyl radicals; (C) CH4/CH3∙/H2/H∙ mixture with the ratio of 4 : 4 : 1 : 2 (i.e., 80 methane molecules, 80 methyl radicals, 20 hydrogen molecules and 40 hydrogen radicals); (D) CH4/H∙ mixture with the ratio of 1 : 1 (i.e., 100 methane molecules, 100 hydrogen radicals); (E) 30 C10H12 (THN) molecules. The periodic cubic box size for all the five models were set as 50 Å × 50 Å × 50 Å (1 Å = 0.1 nm), as shown in Fig. 1. Models A and B were used to study the mechanism of methyl radical collision reaction in the gas phase. Models C and D were used to study the effect of hydrogen radical in the gas-phase reaction. Model E was used to study the thermal cracking behavior of THN at high temperature. Three parallel simulations were performed for each model, one typical simulation result of products distribution is shown in the text, the others are shown in Figs. S1–S24 (see Supporting Information).
Figure 1
Canonical NVT ensemble molecular dynamics simulations were performed using ReaxFF package in Amsterdam Density Functional (ADF) software suite 52. The time step was set as 0.25 fs. The simulation started from a pre-equilibrium process at the first 2.5 ps and the total simulation time was set as 10 ns. The simulation and data collection temperature was 1363 K for model A/B/C/D, in agreement with the experimental temperature. The temperature 1500 and 2000 K were set for model E to accelerate THN cracking. The temperature was controlled using a Nosé-Hoover chains 53. The initial force field parameters for C and H were obtained from the literature 54. The identification and quantification of elementary reactions of ReaxFF trajectories were based on ChemTraYzer scripts 55.
In order to verify the reliability of reactive force field parameters in our model systems, the reaction energies (ΔE) and reaction barriers (Ea) of the representative elementary reactions were calculated by ReaxFF method in ADF software suite and high-level quantum chemistry method (CCSD(T)/cc-pVTZ//M06-2x/cc-pVTZ) 56-58 in Gaussian-09 software 59, which is shown in Table 1. The maximum difference in energy values is within 0.40 eV, therefore ReaxFF method could describe the reaction trends of our model systems qualitatively and quantitatively.
Table 1
Elementary step ΔECCSD(T)/ΔEReaxFF (eV) EaCCSD(T)/EaReaxFF (eV) 2CH3∙ → C2H6 −1.25/−0.85 N/A CH3∙ + H∙ → CH4 −2.45/−2.43 N/A CH3∙ + H2 → CH4 + H∙ 0.48/0.39 2.04/2.07 C2H6 + H∙ → C2H5∙ + H2 −0.82/−0.81 1.54/1.82 C2H6 + CH3∙ → C2H5∙ + CH4 −0.35/−0.42 2.09/2.40 C2H5∙ → C2H4 + H∙ 0.36/0.51 1.77/2.03 3. Results and discussion
3.1 The transformation of methyl radicals
To examine the possibility of the gas-phase reaction mechanism proposed by the original literature 48. ReaxFF simulations were performed from the CH4(80)/CH3∙(80)/H2(40) mixture in a 50 Å × 50 Å × 50 Å box, in agreement with the theoretical products distribution upon methane conversion at 48.1% in the experiment. The population of the products obtained from 10 ns simulation is presented in Fig. 2a. It shows that the number of methyl radicals drop rapidly from 80 to 2 mainly due to the collision between methyl radicals to produce ethane. Hydrogen molecules decrease from 40 to 25 in the first 3 ns and methane as one of the original reactants increase from 80 to 115, which owes to the reactions CH3∙ + H2 → CH4 + H∙ and CH3∙ + H∙ → CH4. After 3 ns, the number of species and the potential energy curve (Fig. S1) become steady. The main products obtained are ethane and methane with the number of 18 and 34, respectively. In this simulation, C2H5∙ radicals, C3H8 and C2H4 molecules are also detected, but the number of each species is all less than 5. The reaction mechanism of the simulated trajectories is revealed in detail in the inset of Fig. 2a. Other two parallel simulations with the same reactants composition were also performed, the results of which are similar to those of the simulation above (Figs. S2 and S3).
Figure 2
Another gas-phase reaction system with 200 pure methyl radicals in a 50 Å × 50 Å × 50 Å box was also simulated by ReaxFF. As shown in Fig. 2b, ethane and methane increase rapidly as methyl radical consumes at the beginning of the simulation. But at the end of the simulation, ethane is the only dominant product with the number of 78 due to the intensive collision between methyl radicals. When the reaction system reaches to the equilibrium at 10 ns, only four ethylene molecules are detected due to the dehydrogenation of ethyl radicals, similar to the reaction mechanism revealed by DFT calculations 48. Guo et al. 48 proposed that the carbon chain might grow through C2H3∙ radicals that produced via the reaction of C2H4 with CH3∙ or H∙. However, at the end of the simulation the number of C2H4, CH3∙ and H∙ is only 4, 2 and 0, respectively, indicating the formation of large amount of C2H3∙ radicals seems unlikely.
The ReaxFF simulations of the two sets of reaction systems demonstrate that the gas-phase reactions between CH3∙ radicals or CH4/CH3∙/H2 mixture cannot yield the main products (ethylene and aromatics) reported in the experiment 48. Instead, the main products only stop at ethane/methane and a tiny proportion of ethylene/propane, suggesting that the gas-phase reaction mechanism between desorbed CH3∙ radical and others is not the only mechanism for the methane conversion on single-atom iron catalyst.
3.2 The role of hydrogen radical
Hydrogen radical is considered to be a critical intermediate to activate C-H bond and finally form C2H4 products in the gas-phase mechanism. In this section, two models containing some hydrogen radicals as the initial reactants were used to explore the role of hydrogen radicals in the reaction. The first model is CH4(80)/CH3∙(80)/H2(20)/H∙(40) mixture in a 50 Å × 50 Å × 50 Å box. As shown in Fig. 3a, the number of hydrogen radicals falls steeply to zero in the first 16 ps. By exploring the elementary reaction, the H∙ radicals mainly combine with CH3∙ and H∙ to form CH4 and H2 molecules respectively, which is demonstrated by the blue and celeste lines in Fig. 3a that the methane and hydrogen molecules increased rapidly in the first 16 ps. Specifically, 20 methane and 9 hydrogen molecules are produced. Compared with the role of H∙ radicals in activating C-H bond, H∙ radicals are more likely to combine with other radicals to form stable molecules, which can also be concluded from Table 1. The combination of CH3∙ and H∙ is about −2.45 eV exothermally by quantum chemistry calculation at 1363 K. While the activation of C-H bond of ethane molecule by H∙ needs to overcome a barrier of 1.54 eV. At the same time, methyl radicals can combine with themselves to form ethane molecules, as shown by the green line in Fig. 3a. After 75 ps of the simulation, the number of hydrogen molecules has a downward trend because methyl radicals can react with hydrogen molecules to form methane. The whole reaction reach balance after 3 ns simulation, as shown in Fig. S4. In the whole simulation process, only one ethylene molecule is detected, demonstrating H radicals do not play a key role to produce ethylene molecules with high selectivity.
Figure 3
In order to study the role of hydrogen radicals in activating methane, a second model with CH4(100)/H∙(100) mixture in a 50 Å × 50 Å × 50 Å box was constructed, as shown in Fig. 3b. At the first 70 ps of the simulation, the number of hydrogen radicals decrease rapidly and consume completely (Fig. 3b). Fifty-five hydrogen molecules are produced at the same time, about forty of them come from the H∙ + H∙ → H2 reaction and the others are produced from the CH4 + H∙ → CH3∙ + H2 reaction. The hydrogen radical does play an important role to activate methane and produce methyl radical because of the increased collision probability of H∙ and CH4 in this model system. However, the methyl radicals can hardly combine with themselves to form ethane molecules because of their relatively small quantity, whereas they are easier to collide with H2 to form CH4 again due to the large number of H2 molecules around, which is shown by the blue line in Fig. 3b that methane molecules start to increase after 100 ps simulation. After 2 ns simulation, the reaction reach balance as shown in Fig. S5. In the whole simulation process, we only detect one ethylene molecule, which is formed by the dehydrogenation of ethane molecules. From this simulation, the large amount of H radicals could also hardly to promote the formation of ethylene molecules.
From these two simulations, the H∙ radicals indeed help to activate methane to some extent. But compared with the activation of methane, the H∙ radicals are more likely to combine with themselves and methyl radicals to form the H2 and methane. Therefore, the presence of hydrogen radical has very limited effect on the selective production of ethylene.
3.3 Thermal cracking of C10H12
According to the experimental report that hydrogen radicals produced by THN thermal decomposition can promote methane conversion to olefins and aromatics, the model with 30 THN molecules in a 50 Å × 50 Å × 50 Å box was built. In order to speed up the thermal cracking process, two high temperatures, 1500 and 2000 K, were adopted. The simulation results are shown in Fig. 4a and Fig. 4b, respectively. After 10 ns simulation, the number of C10H12 molecules decreased from 30 to 6 due to the thermal cracking in Fig. 4a. The maximum number of H radical is only 3, which is not as many as expected. Interestingly, ethylene molecules can be formed through thermal cracking (C10H12 → C8H8 + C2H4) from the simulation, as shown by the orange line in Fig. 4a. The maximum number of ethylene is 8, which may be one of the reasons for the promoted olefins production by adding THN in the experiment. The acetylene and hydrogen molecules can also be detected in the THN thermal cracking process from the simulation (Fig. 4a).
Figure 4
The other simulation at 2000 K shows similar tendency and phenomenon with those at 1500 K. Due to the higher simulating temperature, thirty C10H12 molecules completely cracked in the first 5 ns. The maximum number of hydrogen radicals is also below 5. The maximum number of C2H4, C2H2 and H2 reach 20, 48 and 33 respectively, which are the main cracking products. From these two simulations, the C10H12 thermal cracking could produce C2H4 at high temperature, which thus promotes the production of ethylene.
4. Conclusions
Five models were used to study the gas-phase mechanism of methane nonoxidative conversion by ReaxFF force field methods. It is found that only gas-phase methyl radicals released from the surface after methane dissociation can hardly produce ethylene. Instead, they are more inclined to combine with themselves to mainly produce ethane. The hydrogen radicals in the gas phase do play a role in activating methane to some extent, but they are also more likely to combine with themselves to form H2 and difficult to promote ethylene production. The C2H4 could be formed at high temperature by thermal cracking of C10H12 molecules. Overall, the highly selective production of ethylene products can hardly be achieved via sole gas-phase mechanism, thus the catalyst surface is not only important in the activation of methane but also critical in the following transformation steps.
Supporting Information: available free of charge via the internet at http://www.whxb.pku.edu.cn.
-
-
[1]
Amos, R. D. Mol. Phys. 1979, 38, 33. doi: 10.1080/00268977900101511
-
[2]
Zhan, C. G.; Nichols, J. A.; Dixon, D. A. J. Phys. Chem. A 2003, 107, 4184. doi: 10.1021/jp0225774
-
[3]
Luo, Y. R. Comprehensive Handbook of Chemical Bond Energies; CRC Press: Boca Raton, 2007; pp. 19–145.
-
[4]
Lunsford, J. H. Catal. Today 2000, 63, 165. doi: 10.1016/S0920-5861(00)00456-9
-
[5]
Schwarz, H. Angew. Chem. Int. Ed. 2011, 50, 10096. doi: 10.1002/anie.201006424
-
[6]
Tang, P.; Zhu, Q. J.; Wu, Z. X.; Ma, D. Energy Environ. Sci. 2014, 7, 2580. doi: 10.1039/c4ee00604f
-
[7]
Weaver, J. F.; Hakanoglu, C.; Antony, A.; Asthagiri, A. Chem. Soc. Rev. 2014, 43, 7536. doi: 10.1039/c3cs60420a
-
[8]
Spivey, J. J.; Hutchings, G. Chem. Soc. Rev. 2014, 43, 792. doi: 10.1039/c3cs60259a
-
[9]
Horn, R.; Schlogl, R. Catal. Lett. 2015, 145, 23. doi: 10.1007/s10562-014-1417-z
-
[10]
Zhao, Z. J.; Chiu, C. C.; Gong, J. L. Chem. Sci. 2015, 6, 4403. doi: 10.1039/c5sc01227a
-
[11]
Olivos-Suarez, A. I.; Szecsenyi, A.; Hensen, E. J. M.; Ruiz-Martinez, J.; Pidko, E. A.; Gascon, J. ACS Catal. 2016, 6, 2965. doi: 10.1021/acscatal.6b00428
-
[12]
Schwach, P.; Pan, X. L.; Bao, X. H. Chem. Rev. 2017, 117, 8497. doi: 10.1021/acs.chemrev.6b00715
-
[13]
Vernon, P. D. F.; Green, M. L. H.; Cheetham, A. K.; Ashcroft, A. T. Catal. Today 1992, 13, 417. doi: 10.1016/0920-5861(92)80167-L
-
[14]
York, A. P. E.; Xiao, T. C.; Green, M. L. H. Top. Catal. 2003, 22, 345. doi: 10.1023/A:1023552709642
-
[15]
Jones, G.; Jakobsen, J. G.; Shim, S. S.; Kleis, J.; Andersson, M. P.; Rossmeisl, J.; Abild-Pedersen, F.; Bligaard, T.; Helveg, S.; Hinnemann, B.; et al. J. Catal. 2008, 259, 147. doi: 10.1016/j.jcat.2008.08.003
-
[16]
Li, D. L.; Nakagawa, Y.; Tomishige, K. Appl. Catal. A 2011, 408, 1. doi: 10.1016/j.apcata.2011.09.018
-
[17]
Pakhare, D.; Spivey, J. Chem. Soc. Rev. 2014, 43, 7813. doi: 10.1039/c3cs60395d
-
[18]
Keller, G. E.; Bhasin, M. M. J. Catal. 1982, 73, 9. doi: 10.1016/0021-9517(82)90075-6
-
[19]
Ito, T.; Wang, J. X.; Lin, C. H.; Lunsford, J. H. J. Am. Chem. Soc. 1985, 107, 5062. doi: 10.1021/ja00304a008
-
[20]
Hutchings, G. J.; Scurrell, M. S.; Woodhouse, J. R. Chem. Soc. Rev. 1989, 18, 251. doi: 10.1039/cs9891800251
-
[21]
Lunsford, J. H. Angew. Chem. Int. Ed. 1995, 34, 970. doi: 10.1002/anie.199509701
-
[22]
Groothaert, M. H.; Smeets, P. J.; Sels, B. F.; Jacobs, P. A.; Schoonheydt, R. A. J. Am. Chem. Soc. 2005, 127, 1394. doi: 10.1021/ja047158u
-
[23]
Palkovits, R.; Antonietti, M.; Kuhn, P.; Thomas, A.; Schuth, F. Angew. Chem. Int. Ed. 2009, 48, 6909. doi: 10.1002/anie.200902009
-
[24]
Kwapien, K.; Paier, J.; Sauer, J.; Geske, M.; Zavyalova, U.; Horn, R.; Schwach, P.; Trunschke, A.; Schlogl, R. Angew. Chem. Int. Ed. 2014, 53, 8774. doi: 10.1002/anie.201310632
-
[25]
Grundner, S.; Markovits, M. A. C.; Li, G.; Tromp, M.; Pidko, E. A.; Hensen, E. J. M.; Jentys, A.; Sanchez-Sanchez, M.; Lercher, J. A. Nat. Commun. 2015, 6, 7546. doi: 10.1038/ncomms8546
-
[26]
Ikuno, T.; Zheng, J.; Vjunov, A.; Sanchez-Sanchez, M.; Ortuno, M. A.; Pahls, D. R.; Fulton, J. L.; Camaioni, D. M.; Li, Z. Y.; Ray, D.; et al. J. Am. Chem. Soc. 2017, 139, 10294. doi: 10.1021/jacs.7b02936
-
[27]
Sushkevich, V. L.; Palagin, D.; Ranocchiari, M.; van Bokhoven, J. A. Science 2017, 356, 523. doi: 10.1126/science.aam9035
-
[28]
Wang, P. W.; Zhao, G. F.; Wang, Y.; Lu, Y. Sci. Adv. 2017, 3, e1603180. doi: 10.1126/sciadv.1603180
-
[29]
Xie, J. J.; Jin, R. X.; Li, A.; Bi, Y. P.; Ruan, Q. S.; Deng, Y. C.; Zhang, Y. J.; Yao, S. Y.; Sankar, G.; Ma, D.; et al. Nat. Catal. 2018, 1, 889. doi: 10.1038/s41929-018-0170-x
-
[30]
Wang, L. S.; Tao, L. X.; Xie, M. S.; Xu, G. F.; Huang, J. S.; Xu, Y. D. Catal. Lett. 1993, 21, 35. doi: 10.1007/BF00767368
-
[31]
Weckhuysen, B. M.; Wang, D. J.; Rosynek, M. P.; Lunsford, J. H. Angew. Chem. Int. Ed. 1997, 36, 2374. doi: 10.1002/anie.199723741
-
[32]
Zhang, C. L.; Li, S. A.; Yuan, Y.; Zhang, W. X.; Wu, T. H.; Lin, L. W. Catal. Lett. 1998, 56, 207. doi: 10.1023/A:1019046104593
-
[33]
Weckhuysen, B. M.; Wang, D. J.; Rosynek, M. P.; Lunsford, J. H. J. Catal. 1998, 175, 338. doi: 10.1006/jcat.1998.2010
-
[34]
Xu, Y. D.; Lin, L. W. Appl. Catal. A 1999, 188, 53. doi: 10.1016/S0926-860x(99)00210-0
-
[35]
Liu, S. T.; Wang, L.; Ohnishi, R.; Ichikawa, M. J. Catal. 1999, 181, 175. doi: 10.1006/jcat.1998.2310
-
[36]
Ma, D.; Shu, Y. Y.; Han, X. W.; Liu, X. M.; Xu, Y. D.; Bao, X. H. J. Phys. Chem. B 2001, 105, 1786. doi: 10.1021/jp002011k
-
[37]
Su, L. L.; Ma, D.; Liu, X. M.; Xu, Y. D.; Bao, X. H. Chin. J. Catal. 2002, 23, 41. doi: 10.3321/j.issn:0253-9837.2002.01.010
-
[38]
Xu, Y. D.; Bao, X. H.; Lin, L. W. J. Catal. 2003, 216, 386. doi: 10.1016/S0021-9517(02)00124-0
-
[39]
Su, L. L.; Liu, L.; Zhuang, J. Q.; Wang, H. X.; Li, Y. G.; Shen, W. J.; Xu, Y. D.; Bao, X. H. Catal. Lett. 2003, 91, 155. doi: 10.1023/B:CATL.0000007149.48132.5a
-
[40]
Ismagilov, Z. R.; Matus, E. V.; Tsikoza, L. T. Energy Environ. Sci. 2008, 1, 526. doi: 10.1039/b810981h
-
[41]
Gao, J.; Zheng, Y. T.; Jehng, J. M.; Tang, Y. D.; Wachs, I. E.; Podkolzin, S. G. Science 2015, 348, 686. doi: 10.1126/science.aaa7048
-
[42]
Sun, C. Y.; Fang, G. Z.; Guo, X. G.; Hu, Y. L.; Ma, S. Q.; Yang, T. H.; Han, J.; Ma, H.; Tan, D. L.; Bao, X. H. J. Energy Chem. 2015, 24, 257. doi: 10.1016/S2095-4956(15)60309-6
-
[43]
Tan, P. L. J. Catal. 2016, 338, 21. doi: 10.1016/j.jcat.2016.01.027
-
[44]
Lai, Y.; Veser, G. Catal. Sci. Technol. 2016, 6, 5440. doi: 10.1039/c5cy02258d
-
[45]
Sun, K. D.; Ginosar, D. M.; He, T.; Zhang, Y. L.; Fan, M. H.; Chen, R. P. Ind. Eng. Chem. Res. 2018, 57, 1768. doi: 10.1021/acs.iecr.7b04707
-
[46]
陈强, 姜利学, 李海方, 陈娇娇, 赵艳霞, 何圣贵.物理化学学报, 2019, 35, 1014. doi: 10.3866/PKU.WHXB201811039Chen, Q.; Jiang, L. X.; Li, H. F.; Chen, J. J.; Zhao, Y. X.; He, S. G. Acta Phys. -Chim. Sin. 2019, 35, 1014. doi: 10.3866/PKU.WHXB201811039
-
[47]
王丹, 丁迅雷, 廖珩璐, 戴佳钰.物理化学学报, 2019, 35, 1005. doi: 10.3866/PKU.WHXB201809006Wang, D.; Ding, X. L.; Liao, H. L.; Dai, J. Y. Acta Phys. -Chim. Sin. 2019, 35, 1005. doi: 10.3866/PKU.WHXB201809006
-
[48]
Guo, X. G.; Fang, G. Z.; Li, G.; Ma, H.; Fan, H. J.; Yu, L.; Ma, C.; Wu, X.; Deng, D. H.; Wei, M. M.; et al. Science 2014, 344, 616. doi: 10.1126/science.1253150
-
[49]
Hao, J. Q.; Schwach, P.; Fang, G. Z.; Guo, X. G.; Zhang, H. L.; Shen, H.; Huang, X.; Eggart, D.; Pan, X. L.; Bao, X. H. ACS Catal. 2019, 9, 9045. doi: 10.1021/acscatal.9b01771
-
[50]
Kim, S. K.; Kim, H. W.; Han, S. J.; Lee, S. W.; Shin, J.; Kim, Y. T. Commun. Chem. 2020, 3, 58. doi: 10.1038/s42004-020-0306-1
-
[51]
Liu, Y.; Liu, J. C.; Li, T. H.; Duan, Z. H.; Zhang, T. Y.; Yan, M.; Li, W. L.; Xiao, H.; Wang, Y. G.; Chang, C. R.; et al. Angew. Chem. Int. Ed. 2020, 59, 18586. doi: 10.1002/anie.202003908
-
[52]
van Duin, A. C. T.; Goddard, W. A.; Islam, M. M.; van Schoot, H.; Trnka, T.; Yakovlev, A. L. ReaxFF, 2017, SCM, Theoretical Chemistry, Vrije Universiteit, Amsterdam, The Netherlands, http://www.scm.com
-
[53]
Martyna, G. J.; Klein, M. L.; Tuckerman, M. J. Chem. Phys. 1992, 97, 2635. doi: 10.1063/1.463940
-
[54]
Chenoweth, K.; van Duin, A. C. T.; Goddard, W. A. J. Phys. Chem. A 2008, 112, 1040. doi: 10.1021/jp709896w
-
[55]
Dontgen, M.; Przybylski-Freund, M. D.; Kroger, L. C.; Kopp, W. A.; Ismail, A. E.; Leonhard, K. J. Chem. Theory Comput. 2015, 11, 2517. doi: 10.1021/acs.jctc.5b00201
-
[56]
Zhao, Y.; Truhlar, D. G. Theor. Chem. Acc. 2008, 120, 215. doi: 10.1007/s00214-007-0401-8
-
[57]
Dunning, T. H. J. Chem. Phys. 1989, 90, 1007. doi: 10.1063/1.456153
-
[58]
Purvis, G. D.; Bartlett, R. J. J. Chem. Phys. 1982, 76, 1910. doi: 10.1063/1.443164
-
[59]
Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; et al. Gaussian 09, Revision D.01; Gaussian Inc.: Wallingford, CT, 2013.
-
[1]
计量
- PDF下载量: 51
- 文章访问数: 1713
- HTML全文浏览量: 451