Gas-Phase Mechanism Study of Methane Nonoxidative Conversion by ReaxFF Method

Yuan Liu Zenghui Duan Jun Li Chunran Chang

Citation:  Liu Yuan, Duan Zenghui, Li Jun, Chang Chunran. Gas-Phase Mechanism Study of Methane Nonoxidative Conversion by ReaxFF Method[J]. Acta Physico-Chimica Sinica, 2021, 37(11): 201101. doi: 10.3866/PKU.WHXB202011012 shu

基于ReaxFF的甲烷无氧转化气相机理研究

    通讯作者: 常春然, changcr@mail.xjtu.edu.cn
  • 基金项目:

    中国博士后科学基金 2018M630139

    国家自然科学基金(91645203, 22078257), 中国博士后科学基金(2018T111034, 2018M630139), 中央高校基本科研业务费(xtr0218016, cxtd2017004), 陕西省科技创新团队支持计划(2019TD-039)及王宽诚教育基金会资助项目

    国家自然科学基金 22078257

    中央高校基本科研业务费 xtr0218016

    国家自然科学基金 91645203

    中央高校基本科研业务费 cxtd2017004

    陕西省科技创新团队支持计划 2019TD-039

    中国博士后科学基金 2018T111034

摘要: 晶格限域的Fe©SiO2催化剂在甲烷无氧直接转化生成乙烯的反应中表现出优异的性能。但由于反应条件苛刻,对该反应的分子机理研究一直存在较大的挑战。本文采用反应力场的方法模拟近反应条件下甲烷无氧直接转化气相机理,发现当气相只有甲基自由基存在时,很难产生高选择性乙烯产物。当在气相中加入氢自由基时,虽能在一定程度上增强甲烷的活化,但同样较难生成乙烯。高温下热裂解C10H12分子能同时产生氢自由基和乙烯分子,能合理地解释实验中加入C10H12分子可以在一定程度上提高乙烯选择性和甲烷转化率的现象。总之,甲烷无氧直接转化高选择性生成乙烯很难通过单纯的气相反应机理来实现,进而推断催化剂表面在甲烷活化和转化的整个过程中起着至关重要的作用。

English

  • 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.

    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

    Figure 1.  The model C of CH4/CH3∙/H2/H∙ mixture for ReaxFF simulation, i.e., 80 methane molecules, 80 methyl radicals, 20 hydrogen molecules and 40 hydrogen radicals in a 50 Å × 50 Å × 50 Å box.

    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

    Table 1.  Comparison of energetics calculated by ReaxFF and quantum chemistry methods of typical elementary reactions.
    下载: 导出CSV
    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

    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

    Figure 2.  (a) Products distribution of ReaxFF simulations beginning from gas-phase CH4(80)/CH3∙(80)/H2(40) mixture as a function of time, (b) Products distribution of ReaxFF simulations beginning from 200 gas-phase CH3∙ radicals as a function of time.

    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.

    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

    Figure 3.  (a) Products distribution of ReaxFF simulations beginning from gas-phase CH4(80)/CH3∙(80)/H2(20)/H∙(40) mixture as a function of time, (b) Products distribution of ReaxFF simulations beginning from CH4(100)/H∙(100) atoms as a function of time.

    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.

    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

    Figure 4.  (a) Products distribution of ReaxFF simulations beginning from gas-phase C10H12(30) as a function of time at 1500 K, (b) Products distribution of ReaxFF simulations beginning from C10H12(30) as a function of time at 2000 K.

    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.

    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. [1]

      Amos, R. D. Mol. Phys. 1979, 38, 33. doi: 10.1080/00268977900101511

    2. [2]

      Zhan, C. G.; Nichols, J. A.; Dixon, D. A. J. Phys. Chem. A 2003, 107, 4184. doi: 10.1021/jp0225774

    3. [3]

      Luo, Y. R. Comprehensive Handbook of Chemical Bond Energies; CRC Press: Boca Raton, 2007; pp. 19–145.

    4. [4]

      Lunsford, J. H. Catal. Today 2000, 63, 165. doi: 10.1016/S0920-5861(00)00456-9

    5. [5]

      Schwarz, H. Angew. Chem. Int. Ed. 2011, 50, 10096. doi: 10.1002/anie.201006424

    6. [6]

      Tang, P.; Zhu, Q. J.; Wu, Z. X.; Ma, D. Energy Environ. Sci. 2014, 7, 2580. doi: 10.1039/c4ee00604f

    7. [7]

      Weaver, J. F.; Hakanoglu, C.; Antony, A.; Asthagiri, A. Chem. Soc. Rev. 2014, 43, 7536. doi: 10.1039/c3cs60420a

    8. [8]

      Spivey, J. J.; Hutchings, G. Chem. Soc. Rev. 2014, 43, 792. doi: 10.1039/c3cs60259a

    9. [9]

      Horn, R.; Schlogl, R. Catal. Lett. 2015, 145, 23. doi: 10.1007/s10562-014-1417-z

    10. [10]

      Zhao, Z. J.; Chiu, C. C.; Gong, J. L. Chem. Sci. 2015, 6, 4403. doi: 10.1039/c5sc01227a

    11. [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. [12]

      Schwach, P.; Pan, X. L.; Bao, X. H. Chem. Rev. 2017, 117, 8497. doi: 10.1021/acs.chemrev.6b00715

    13. [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. [14]

      York, A. P. E.; Xiao, T. C.; Green, M. L. H. Top. Catal. 2003, 22, 345. doi: 10.1023/A:1023552709642

    15. [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. [16]

      Li, D. L.; Nakagawa, Y.; Tomishige, K. Appl. Catal. A 2011, 408, 1. doi: 10.1016/j.apcata.2011.09.018

    17. [17]

      Pakhare, D.; Spivey, J. Chem. Soc. Rev. 2014, 43, 7813. doi: 10.1039/c3cs60395d

    18. [18]

      Keller, G. E.; Bhasin, M. M. J. Catal. 1982, 73, 9. doi: 10.1016/0021-9517(82)90075-6

    19. [19]

      Ito, T.; Wang, J. X.; Lin, C. H.; Lunsford, J. H. J. Am. Chem. Soc. 1985, 107, 5062. doi: 10.1021/ja00304a008

    20. [20]

      Hutchings, G. J.; Scurrell, M. S.; Woodhouse, J. R. Chem. Soc. Rev. 1989, 18, 251. doi: 10.1039/cs9891800251

    21. [21]

      Lunsford, J. H. Angew. Chem. Int. Ed. 1995, 34, 970. doi: 10.1002/anie.199509701

    22. [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. [23]

      Palkovits, R.; Antonietti, M.; Kuhn, P.; Thomas, A.; Schuth, F. Angew. Chem. Int. Ed. 2009, 48, 6909. doi: 10.1002/anie.200902009

    24. [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. [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. [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. [27]

      Sushkevich, V. L.; Palagin, D.; Ranocchiari, M.; van Bokhoven, J. A. Science 2017, 356, 523. doi: 10.1126/science.aam9035

    28. [28]

      Wang, P. W.; Zhao, G. F.; Wang, Y.; Lu, Y. Sci. Adv. 2017, 3, e1603180. doi: 10.1126/sciadv.1603180

    29. [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. [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. [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. [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. [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. [34]

      Xu, Y. D.; Lin, L. W. Appl. Catal. A 1999, 188, 53. doi: 10.1016/S0926-860x(99)00210-0

    35. [35]

      Liu, S. T.; Wang, L.; Ohnishi, R.; Ichikawa, M. J. Catal. 1999, 181, 175. doi: 10.1006/jcat.1998.2310

    36. [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. [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. [38]

      Xu, Y. D.; Bao, X. H.; Lin, L. W. J. Catal. 2003, 216, 386. doi: 10.1016/S0021-9517(02)00124-0

    39. [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. [40]

      Ismagilov, Z. R.; Matus, E. V.; Tsikoza, L. T. Energy Environ. Sci. 2008, 1, 526. doi: 10.1039/b810981h

    41. [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. [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. [43]

      Tan, P. L. J. Catal. 2016, 338, 21. doi: 10.1016/j.jcat.2016.01.027

    44. [44]

      Lai, Y.; Veser, G. Catal. Sci. Technol. 2016, 6, 5440. doi: 10.1039/c5cy02258d

    45. [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. [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. [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. [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. [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. [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. [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. [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. [53]

      Martyna, G. J.; Klein, M. L.; Tuckerman, M. J. Chem. Phys. 1992, 97, 2635. doi: 10.1063/1.463940

    54. [54]

      Chenoweth, K.; van Duin, A. C. T.; Goddard, W. A. J. Phys. Chem. A 2008, 112, 1040. doi: 10.1021/jp709896w

    55. [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. [56]

      Zhao, Y.; Truhlar, D. G. Theor. Chem. Acc. 2008, 120, 215. doi: 10.1007/s00214-007-0401-8

    57. [57]

      Dunning, T. H. J. Chem. Phys. 1989, 90, 1007. doi: 10.1063/1.456153

    58. [58]

      Purvis, G. D.; Bartlett, R. J. J. Chem. Phys. 1982, 76, 1910. doi: 10.1063/1.443164

    59. [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.

  • 加载中
计量
  • PDF下载量:  51
  • 文章访问数:  1713
  • HTML全文浏览量:  451
文章相关
  • 发布日期:  2021-11-15
  • 收稿日期:  2020-11-03
  • 接受日期:  2020-11-26
  • 修回日期:  2020-11-25
  • 网络出版日期:  2020-12-03
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

/

返回文章