RDX及其衍生物高温热解的反应分子动力学模拟

引用本文: 彭莉娟, 姚倩, 王静波, 李泽荣, 朱权, 李象远. RDX及其衍生物高温热解的反应分子动力学模拟[J]. 物理化学学报, 2017, 33(4): 745-754. doi: 10.3866/PKU.WHXB201701161 shu

Citation: PENG Li-Juan, YAO Qian, WANG Jing-Bo, LI Ze-Rong, ZHU Quan, LI Xiang-Yuan. Pyrolysis of RDX and Its Derivatives via Reactive Molecular Dynamics Simulations[J]. Acta Physico-Chimica Sinica, 2017, 33(4): 745-754. doi: 10.3866/PKU.WHXB201701161 shu

RDX及其衍生物高温热解的反应分子动力学模拟

彭莉娟 ^1, ,
姚倩 ^2, ,
王静波 ^3, ,
李泽荣 ^2, ,
朱权 ^3, ,
李象远 ^3,

1.
四川大学空天科学与工程学院, 成都 610065;
2.
四川大学化学学院, 成都 610064;
3.
四川大学化学工程学院, 成都 610065
基金项目:
国家自然科学基金（91441132）和教育部新世纪优秀人才支持计划（NCET-13-0398）资助项目

摘要: 本文采用ReaxFF 反应分子动力学方法，研究了RDX及其衍生物晶体在高温条件下（2000、2500 和3000 K）的热分解机理以及主要产物随时间的变化情况。结果表明：RDX及其衍生物晶体热解的第一步反应均为N―NO₂键断裂生成NO₂分子，随后反应由于六元环上和侧链基团的不同导致侧链脱除或开环反应的顺序不同。在这四种单胞体系热解中，NO₂和NO分子为共同的中间产物，形成之后迅速发生次级反应并最终生成N₂。各体系热解终产物一致，均为N₂、H₂O和CO₂，其中N₂分子数最多，大于20 个。由于原始分子结构和组成中C/N比、H/O比的不同，各体系热解后H₂O和CO₂分子数目相差较大。不同温度下，各单胞体系热分解生成的最大含碳团簇中C原子数均较小。在进一步超胞体系的模拟中，RDX和RDX-D2体系生成的含碳团簇中C原子数分别达到约30 和16 个，远高于单胞模拟，且受温度影响较大；而RDX-D1 和RDX-D3 单胞或超胞模拟结果相近，均未生成含碳团簇，仅存在小分子含碳碎片。因此，初始分子的结构和元素比对含碳团簇的生成有明显影响。

关键词:

English

Pyrolysis of RDX and Its Derivatives via Reactive Molecular Dynamics Simulations

1.
School of Aeronautics & Astronautics, Sichuan University, Chengdu 610065, P. R. China;
2.
College of Chemistry, Sichuan University, Chengdu 610064, P. R. China;
3.
School of Chemical Engineering, Sichuan University, Chengdu 610065, P. R. China
Received Date: 22 October 2016
Revised Date: 16 January 2017
Available Online: 16 January 2017
Fund Project:
The project was supported by the National Natural Science Foundation of China (91441132) and Program for New Century Excellent Talents in University of Ministry of Education of China (NCET-13-0398).

Abstract: The thermal decomposition mechanisms of cyclotrimethylene trinitramine (RDX) and its derivatives are investigated using ReaxFF reactive molecular dynamics simulation at high temperatures (2000, 2500, and 3000 K). It is shown that the first pyrolysis step of RDX and its derivatives is the N―NO₂ homolytic cleavage to produce NO₂; however, the subsequent reaction mechanism is completely different due to the different sixmembered rings and side chain groups. In these four model systems, NO₂ and NO molecules are common intermediates in the thermal decomposition process, which eventually transform to N₂ in the following reactions. The most stable products are N₂, H₂O, and CO₂ in the thermal process, among which N₂ has the maximum molecular number (more than 20). The numbers of cracked H₂O and CO₂ molecules are very different in the four model systems due to the different C/N and H/O ratios. At different temperatures, the maximum numbers of carbon atoms in the carbon clusters are all small in the unit cell simulations of the four systems. In the further super cell simulation for RDX and RDX-D2 the numbers of carbon atoms reach about 30 and 16, respectively, which are higher than those from the unit cell simulation. In both the unit and super cell simulations for RDX-D1 and RDX-D3, the carbon cluster cannot be formed and there exist only small carbon molecular fragments. Therefore, the initial molecular structure and elemental ratio have a large effect on the production of carbon clusters.

Key words:

1. [1]
  Shaw, R.W.; Brill, T. B.; Thompson, D. L. Overviews of Recent Research on Energetic Materials ;World Scientific: NJ, 2005, Vol. 16.
2. [2]
  Teipel, U. Energetic Materials: Particle Processing and Characterization; Wiley-VCH: Weinheim, Germany, 2005.
3. [3]
  Chakraborty, D.; Muller, R. P.; Dasgupta, S.; Goddard, W. A. J. Phys. Chem. A 2000, 104, 2261. doi: 10.1021/jp9936953
4. [4]
  Xu, J.; Zhao, J.; Sun, L. Mol. Simul. 2008, 34, 961. doi: 10.1080/08927020802162892
5. [5]
  Schweigert, I. V. J. Phys. Chem. A 2015, 119, 2747. doi: 10.1021/jp510034p
6. [6]
  Irikura, K. K. J. Phys. Chem. A 2013, 117, 2233. doi: 10.1021/jp310247z
7. [7]
  Behrens, R., Jr.; Bulusu, S. J. Phys. Chem. 1992, 96, 8877. doi: 10.1063/1.462245
8. [8]
  Van Duin, A. C.; Dasgupta, S.; Lorant, F.; Goddard, W. A. J. Phys. Chem. A 2001, 105, 9396. doi: 10.1021/jp004368u
9. [9]
  Cohen, R.; Zeiri, Y.; Wurzberg, E.; Kosloff, R. J. Phys. Chem. A 2007, 111, 11074. doi: 10.1021/jp072121s
10. [10]
  Strachan, A.; van Duin, A. C.; Chakraborty, D.; Dasgupta, S.; Goddard, W. A., III Phys. Rev. Lett. 2003, 91, 098301. doi: 10.1103/PhysRevLett.91.098301
11. [11]
  Zhou, T. T.; Shi, Y. D.; Huang, F. L. Acta Phys. -Chim. Sin. 2012, 28 (11), 2605. [周婷婷, 石一丁, 黄风雷. 物理化学学报, 2012, 28 (11), 2605.]. doi: 10.3866/PKU.WHXB201208031
12. [12]
  Zhang, L.; Chen, L.; Wang, C.; Wu, J. Y. Acta Phys. -Chim. Sin. 2013, 29 (6), 1145. [张力, 陈朗, 王晨, 伍俊英. 物理化学学报, 2013, 29 (6), 1145.]. doi: 10.3866/PKU.WHXB201303221
13. [13]
  Rom, N.; Zybin, S. V.; van Duin, A. C.; Goddard, W. A., III.; Zeiri, Y.; Katz, G.; Kosloff, R. J. Phys. Chem. A 2011, 115, 10181. doi: 10.1021/jp202059v
14. [14]
  Liu, H.; Dong, X.; He, Y. H. Acta Phys. -Chim. Sin. 2014, 30, 232. [刘海, 董晓, 何远航. 物理化学学报, 2014, 30, 232.] doi: 10.3866/PKU.WHXB201312101
15. [15]
  Qian, H. J.; van Duin, A. C.; Morokuma, K.; Irle, S. J. Chem. Theory Comput. 2011, 7, 2040. doi: 10.1021/ct200197v
16. [16]
  Liu, L.; Bai, C.; Sun, H.; Goddard, W. A., III. J. Phys. Chem. A 2011, 115, 4941. doi: 10.1021/jp110435p
17. [17]
  Ge, N. N.; Wei, Y. K.; Ji, G. F.; Chen, X. R.; Zhao, F.; Wei, D.Q. J. Phys. Chem. B 2012, 116, 13696. doi: 10.1021/jp309120t
18. [18]
  Zhang, L. Z.; Zybin, S. V.; Van Duin, A. C. T.; Goddard, W. A., III. J. Energ. Mater. 2010, 28, 92. doi: 10.1080/07370652.2010.504682
19. [19]
  Strachan, A.; Kober, E. M.; van Duin, A. C.; Oxgaard, J.; Goddard, W. A., III. J. Chem. Phys. 2005, 122, 054502. doi: 10.1063/1.1831277
20. [20]
  Zhang, L.; Zybin, S. V.; van Duin, A. C.; Dasgupta, S.; Goddard, W. A., III; Kober, E. M. J. Phys. Chem. A 2009, 113, 10619. doi: 10.1021/jp901353a
21. [21]
  Hakey, P.; Ouellette, W.; Zubieta, J.; Korter, T. Acta Crystallographica Section E: Structure Reports Online 2008, 64
22. [22]
  Qiu, L.; Gong, X.; Ju, X.; Xiao, H. Sci. China Ser. B: Chem. 2008, 51, 1231. doi: 10.1007/s11426-008-0141-1
23. [23]
  Berendsen, H. J.; Postma, J. V.; van Gunsteren, W. F.; DiNola, A. R. H. J.; Haak, J. R. J. Chem. Phys. 1984, 81, 3684. doi: 10.1063/1.448118
24. [24]
  Botcher, T. R.; Wight, C. A. J. Phys. Chem. 1994, 98, 5441. doi: 10.1021/j100072a009
25. [25]
  Oyumi, Y.; Brill, T. B. Combust. Flame 1985, 62, 213. doi: 10.1016/0010-2180(85)90147-6
26. [26]
  Kohno, Y.; Ueda, K.; Imamura, A. J. Phys. Chem. 1996, 100, 4701. doi: 10.1021/jp9503223
27. [27]
  Prabhakaran, K. V.; Bhide, N. M.; Kurian, E. M. Thermochim. Acta 1995, 249, 249. doi: 10.1016/0040-6031(95)90704-1
28. [28]
  Oxley, J. C.; Kooh, A. B.; Szekeres, R.; Zheng, W. J. Phys. Chem. 1994, 98, 7004. doi: 10.1021/j100079a019
29. [29]
  Brill, T. B.; Oyumi, Y. J. Phys. Chem. 1986, 90, 6848. doi: 10.1021/j100284a028
30. [30]
  Mironov, E. V.; Petrov, E. A.; Korets, A. Y. Combust. Explos. Shock Waves. 2004, 40, 473. doi: 10.1023/B:CESW.0000033571.82326.6a
31. [31]
  Anisichkin, V. F.; Dolgushin, D. S.; Petrov, E. A. Combust. Explos. Shock Waves. 1995, 31, 106. doi: 10.1007/BF00755966
32. [32]
  Titov, V. M.; Anisichkin, V. F.; Mal'kov, I. Y. Combust. Explos. Shock Waves. 1989, 25, 372. doi: 10.1007/BF00788819
33. [33]
  Baidakova, M. J. Phys. D: Appl. Phys. 2007, 40, 6300. doi: 10.1088/0022-3727/40/20/S14
34. [34]
  Mal'Kov, I. Y.; Filatov, L. I.; Titov, V. M.; Litvinov, B. V.; Chuvilin, A. L.; Teslenko, T. S. Combust. Explos. Shock Waves. 1993, 29, 542. doi: 10.1007/BF00782983
35. [35]
  Chevrot, G.; Sollier, A.; Pineau, N. J. Chem. Phys. 2012, 136, 084506. doi: 10.1063/1.3686750
36. [36]
  Hobbs, M. L.; Kaneshige, M. J.; Gilbert, D.W.; Marley, S. K.; Todd, S. N. J. Phys. Chem. A 2009, 113, 10474. doi: 10.1021/jp906134f

扫一扫看文章

计量

PDF下载量: 4
文章访问数: 857
HTML全文浏览量: 116

通讯作者: 陈斌, bchen63@163.com

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

RDX及其衍生物高温热解的反应分子动力学模拟

English

Pyrolysis of RDX and Its Derivatives via Reactive Molecular Dynamics Simulations

CCS Chemistry：嵌段共聚物自组装成 “三明治”二维有机材料，实现纳米级压力传感功能

CCS Chemistry：颜徐州、鲍哲南：你的皮肤可能是一片SPMs

CCS Chemistry：工作高效不疲劳，只须让催化剂动起来

CCS Chemistry：静电组装导向聚合- 聚电解质纳米凝胶量化制备新策略

CCS Chemistry：金黄色葡萄球菌醛基脱氢酶是如何识别特异性底物的？

计量

通讯作者: 陈斌, bchen63@163.com

中国化学会秘书处

RDX及其衍生物高温热解的反应分子动力学模拟

English

Pyrolysis of RDX and Its Derivatives via Reactive Molecular Dynamics Simulations

CCS Chemistry：嵌段共聚物自组装成 “三明治”二维有机材料，实现纳米级压力传感功能

CCS Chemistry：颜徐州、鲍哲南： 你的皮肤可能是一片SPMs

CCS Chemistry：工作高效不疲劳，只须让催化剂动起来

CCS Chemistry：静电组装导向聚合- 聚电解质纳米凝胶量化制备新策略

CCS Chemistry：金黄色葡萄球菌醛基脱氢酶是如何识别特异性底物的？

计量

文章相关

通讯作者: 陈斌, bchen63@163.com

中国化学会秘书处

CCS Chemistry：颜徐州、鲍哲南：你的皮肤可能是一片SPMs