RDX及其衍生物高溫熱解的反應分子動力學模擬

彭莉娟 姚 倩 王靜波 李澤榮 朱 權,* 李象遠,*

RDX及其衍生物高溫熱解的反應分子動力學模擬

彭莉娟1姚 倩2王靜波3李澤榮2朱 權3,*李象遠3,*

(1四川大學空天科學與工程學院，成都610065；2四川大學化學學院，成都610064；3四川大學化學工程學院，成都610065)

本文采用ReaxFF反應分子動力學方法，研究了RDX及其衍生物晶體在高溫條件下(2000、2500和3000 K)的熱分解機理以及主要產物隨時間的變化情況。結果表明：RDX及其衍生物晶體熱解的第一步反應均為N―NO2鍵斷裂生成NO2分子，隨后反應由于六元環上和側鏈基團的不同導致側鏈脫除或開環反應的順序不同。在這四種單胞體系熱解中，NO2和NO分子為共同的中間產物，形成之后迅速發生次級反應并最終生成N2。各體系熱解終產物一致，均為N2、H2O和CO2，其中N2分子數最多，大于20個。由于原始分子結構和組成中C/N比、H/O比的不同，各體系熱解后H2O和CO2分子數目相差較大。不同溫度下，各單胞體系熱分解生成的最大含碳團簇中C原子數均較小。在進一步超胞體系的模擬中，RDX和RDX-D2體系生成的含碳團簇中C原子數分別達到約30和16個，遠高于單胞模擬，且受溫度影響較大；而RDX-D1和RDX-D3單胞或超胞模擬結果相近，均未生成含碳團簇，僅存在小分子含碳碎片。因此，初始分子的結構和元素比對含碳團簇的生成有明顯影響。

RDX及其衍生物；高溫熱分解；ReaxFF；含碳團簇；分子動力學

1 引言

含能材料是一類儲存高能量的特殊材料，在一定外界條件刺激下，自身能在短時間內發生氧化還原反應，釋放大量能量，在軍用和民用領域廣泛應用1。高溫熱分解過程和反應機理的研究，不僅可提高對含能材料反應特性的認識，而且還可在原子分子水平上進行改進和設計，從而研制出更好性能的炸藥和推進劑配方。

黑索金(RDX)是一種綜合性能優良的單質炸藥和推進劑。通常RDX存在α和β兩種晶型，常溫常壓下α-RDX性質穩定而β-RDX極不穩定，因此α-RDX具有更高的實用價值2，本文研究均采用α晶型。文獻中，針對RDX的熱分解過程開展了大量的理論和實驗研究。Goddard等3采用從頭算方法研究了單分子RDX熱分解過程，提出了三種主要初始分解途徑：協同分解反應，N―NO2的斷裂反應和連續HONO消去反應，其中連續HONO消去反應被認為是能量最低的反應途徑。趙繼軍等4研究了溫度為3000 K下RDX的熱分解行為，表明熱解后的主要產物為N2、H2O和CO2。Schweigert5采用從頭算分子動力學研究氣相單分子RDX在高溫2000 K下的詳細裂解反應機理。Irikura6采用密度泛函理論計算比較了RDX單分子和雙分子熱解反應途徑，為凝聚相反應機理的研究提供了理論依據。Behrens和Bulusu7從實驗上指出RDX的主要分解路徑為RDX分子去除HNO和兩個HONO分子生成1,3,5-三嗪(TAZ)。雖然，氣相單分子RDX熱解機理研究已經比較成熟，但凝聚相RDX由于高溫熱解中分子間復雜的相互作用，反應機理研究仍不夠深入。

ReaxFF反應力場8－10通過鍵長、鍵級和鍵能之間關系，以及原子間電荷的動態傳遞來模擬原子間化學鍵的斷裂和生成，為研究較大體系的物理化學性質提供了有效途徑。ReaxFF反應力場已經廣泛應用于熱分解11－14、燃燒15,16和爆轟17,18等快速反應過程。采用這一方法，Strachan等19研究了不同溫度和密度下RDX的熱分解過程，表明產物中CO和CO2分子數目很大程度上與系統密度有關，而N2和H2O分子的數目幾乎不受密度的影響。Zhang等20研究了不同溫度和密度下奧克托金(HMX)和三氨基三硝基苯(TATB)晶體的熱分解過程，表明HMX的初始裂解為N―NO2鍵的斷裂反應和HONO消去反應；TATB在低溫下幾乎沒有發現C―NO2鍵的斷裂，高溫下主要為C―NO2鍵斷裂反應和HONO消去反應。但在RDX的熱解模擬中，關于反應物分子結構和元素組成對裂解產物和反應途徑的影響，一直缺乏系統的對比分析，難以對含能材料的結構設計提供指導。

本工作根據RDX分子，改變六元環上和側鏈的取代基團設計衍生物分子，從而改變含能材料中碳、氮含量比，研究其對反應過程中能量釋放和含碳團簇形成過程的影響。通過RDX取代基團的逐步替代設計了三種化合物，它們的結構分別如圖1所示，這三種衍生物的分子式分別為：C2H5O6N7(RDX-D1)、C3H7O6N7(RDX-D2)、C2H4O8N8(RDXD3)。相對于RDX，RDX-D1采用NH取代RDX主環上的CH2，改變了主環結構，提高了N元素含量；相對于RDX-D1，RDX-D2用―CH3取代了主環上NH基團中的氫原子，改變了C元素含量和側鏈結構；相對于RDX-D1，RDX-D3用―NO2取代了主環上NH基團中的氫原子，進一步提高了N元素含量和改變了側鏈結構。本工作采用ReaxFFMD方法，在單胞水平上模擬研究了高溫熱解的初始反應過程和反應產物隨時間變化，在超胞水平上深入分析了反應過程中含碳團簇的演化情況，從而探討分子元素組成和分子結構對RDX類化合物高溫熱解性能的影響。

圖1 RDX及其衍生物的分子結構Fig.1 Structures of the RDX molecule and its derivatives

2 計算模型與方法

2.1 初始構型的構建

RDX模擬采用的初始體系為一個單胞，依據RDX的X射線衍射單晶數據21，運用Materials Studio(MS)程序包Visualizer模塊，構建單胞體系，如圖2(a)所示，該單胞結構包含8個RDX分子，共計168個原子，初始密度為1.86 g·cm－3。由于RDX-D1與RDX分子具有相似的分子結構，故采用和RDX相同的單晶參數，采用MS程序包Visualizer模塊構建單胞體系，密度為1.87 g· cm－3。依據肖鶴鳴等22采用密度泛函理論計算結果，RDX-D3密度選取為1.92 g·cm－3。RDX-D2與RDX-D3具有相似的分子結構，故選用相同的參數。各體系單胞結構如圖2所示，其中灰色、白色、藍色和紅色原子分別表示碳原子、氫原子、氮原子和氧原子，模擬均采用周期性邊界條件。

針對含碳團簇生成過程模擬，對RDX及其衍生物分別搭建了(2×2×1)超胞結構。RDX超胞結構共含有32個分子，共計672個原子；衍生物RDX-D1，RDX-D2和RDX-D3超胞結構所含原子個數分別為640個、736個和704個。

2.2 計算方法

分子動力學模擬采用ReaxFF力場20，選取Velocity Verlet算法，采用ReaxFF-MD方法分別模擬各晶胞體系的高溫熱裂解過程。在ReaxFF反應力場中，各初始晶胞體系分別在低溫10 K下進行幾何優化，獲得最小能量結構，然后在相同溫度采用等溫等容(NVT)系綜進行弛豫，使得模擬體系達到平衡狀態，弛豫時間為5 ps，時間步長選取為0.1 fs。其后，各體系得到的平衡結構分別在2000、2500、3000 K三個溫度下采取NVT系綜進行高溫分解模擬，模擬時間步長選取0.1 fs，總模擬時間為100 ps。由于升溫過程可能會導致化學鍵的斷裂，因此，本文不考慮升溫過程的模擬。在所有的NVT-MD模擬中，各系統溫度均用Berendsen方法23進行調節，溫度調整參數選取100 fs，使溫度在設定值附近波動。當任意原子對的鍵級大于0.3時，判定新化學鍵形成。

3 結果與討論

3.1 初始反應途徑

在3000 K條件下對四個單胞體系的熱分解進行模擬，通過分析原子成鍵情況及各體系在不同時刻下的裂解產物，最終得到高溫熱解的初始反應途徑，如圖3所示。RDX及其三種衍生物熱解過程，第一步反應均為N―NO2鍵發生斷裂生成NO2分子，表明主環和側鏈基團的變化并不影響側鏈硝基斷鍵的初始反應步驟。RDX分子第一步首先發生N―NO2鍵的斷裂，從而生成大量的NO2，與Botcher等在實驗中觀察到RDX快速熱解初期生成大量NO2一致24,25。在隨后反應中，RDX主要發生剩余N―NO2鍵斷裂反應和母體環斷裂反應，這一結果與Goddard等3采用密度泛函理論計算得到的RDX中N―N鍵斷裂(解離能約83.68 kJ·mol－1)優先于協同開環反應的結論一致。

圖2 RDX(a),RDX-D1(b),RDX-D2(c)和RDX-D3(d)單胞結構示意圖Fig.2 Unit cell structure of RDX(a),RDX-D1(b),RDX-D2(c)and RDX-D3(d) color online

圖3 RDX及其衍生物的初始反應途徑Fig.3 Initial reaction path for RDX and its derivatives

RDX由于分子本身的對稱性，側鏈第一和第二個NO2基的消去無選擇性，但在RDX-D1中，由于六元環上的CH2基團被NH基團取代，在熱解過程中優先脫去NH基團鄰位氮原子上的兩個側鏈NO2基，其后，不同于RDX的開環反應，RDX-D1發生HONO脫去反應，之后再開環。在RDX-D1基礎上將環上NH中的氫原子用甲基取代，RDXD2中兩個NO2基脫除也優先發生于N―CH3鄰位氮原子上，其后發生的是甲基脫除反應，然后發生開環反應。因此，與RDX-D1相比，CH3基的引入對鄰位N―NO2的斷鍵影響不大，而對對位NO2的脫除影響較大，使得后續的反應途徑完全不同。將RDX-D2中的甲基用硝基取代，RDX-D3中相鄰的兩個NO2基優先脫除，六元環上形成氮氮雙鍵，其后雙鍵間位上的NO2基發生脫除，然后發生開環反應。

在RDX及其衍生物晶胞熱解過程中主要觀察到三種初始分解機理：N―NO2鍵的斷裂反應，HONO消去反應和分子主環的開環反應，與Goddard等采用密度泛函理論計算的結果一致3,22。對比這四種單胞，第一步反應均為N―NO2鍵斷裂，表明這一通道是能量最低的反應通道，與Imamura等研究的其它硝胺化合物的熱分解機理一致26－29。然后才出現HONO消去、甲基取代基脫除和開環反應，進一步生成小分子碎片。在四種體系的高溫熱解中，側鏈取代基團對N―NO2第一步斷鍵無影響，主要改變了隨后的側鏈脫除或開環反應的順序。

3.2 溫度對勢能和總分子碎片數的影響

在不同溫度條件下，我們模擬了各單胞體系(經過5 ps初始弛豫)勢能隨時間的變化，如圖4所示。在初始熱解時，系統迅速吸收熱量發生初級熱解反應，使得系統內勢能升高，隨后的次級反應生成一系列中間及最終產物，并放出大量的熱，因而，系統勢能發生明顯下降。溫度越高，系統達到化學反應動態平衡所需時間越短，熱解反應越完全，生成的產物越多，勢能下降也越多。因而，各體系勢能的改變與其熱分解反應過程中的化學鍵斷裂和生成緊密相關30。

從圖4(a)可見，系統溫度越低，初始反應越慢，吸熱階段越長；隨著系統溫度的升高，各體系分子初始反應加劇，從而導致次級反應越早開始，系統達到化學反應動態平衡的時間縮短。當溫度為3000 K時，各晶胞體系勢能在20 ps后相對穩定；而當系統處于2000 K的較低溫度時，各體系勢能在100 ps內依然緩慢減少，需要更長時間才能達到穩定，故在低溫條件達到化學反應平衡的時間更長。

圖5為各單胞體系在不同溫度下熱解過程中總分子碎片數隨時間的變化情況。從圖5(a)可見，在高溫初始裂解階段，溫度越高RDX熱解釋放能量的速率越快，物種數增加越快、產生的總分子碎片數越多，在10 ps左右達到峰值后，這些快速熱解產生的分子碎片會進一步轉化，但整體數量維持在50左右。其他衍生物具有類似的變化規律，各晶胞體系的總分子碎片數在前10 ps內快速增加，隨后變化相對比較緩慢，且溫度越高，總分子碎片數增加速度越大。在50 ps后，各晶胞體系的總分子碎片數在3000 K溫度條件下基本維持恒定，而在較低的溫度2000 K條件下仍緩慢增加。綜上，隨熱解溫度升高，各晶胞體系產生的總分子碎片數變多，總分子碎片數增加到最大值的時間縮短，體系釋放的能量也更多。

3.3 主產物及中間產物分析

圖4 四個體系的勢能隨時間變化的關系Fig.4 Time evolution of potential energy for the four materials

圖5 四個體系的熱解總分子碎片數隨時間變化關系Fig.5 Time evolution of total fragments for the four materials

RDX及其衍生物在3000 K時的主要分解產物N2、H2O和CO2及中間產物NO和NO2分子數目隨時間的變化如圖6所示。在各體系主產物中，N2是生成數目最多的產物，主要來自次級反應中NO2和NO的轉化，其主要機理為NO2或NO與含碳分子碎片反應生成N2。在前20 ps內，NO2和NO均迅速增大然后減少，此過程中N2數目一直迅速增加到約20，說明消耗掉的大部分NO2和NO生成了N2。在這四種不同體系中，主產物分子數目相對較多的還有H2O和CO2，這三種主要產物數目的大小順序為：N2＞H2O＞CO2。

在四種不同的晶胞體系中，主要分解產物H2O和CO2分子數目相差較大。表1列出了各體系分子的元素組成比以及主產物CO2和H2O分子的最終數目。從表1可以清晰看出RDX和RDX-D2分子中C/N比大于等于0.5，高于RDX-D1和RDX-D3分子中的C/N比；RDX和RDX-D2體系中碳原子含量相多較高，熱解穩定后生成的CO2分子數均小于10，這說明體系中的碳原子的最終轉化途徑并不以生成穩定產物CO2為主，而是轉變成了其它的含碳分子碎片。與此對應，RDX和RDX-D2由于其分子構成的H/O元素比(≥1)高于另外兩種衍生物，熱解后生成H2O分子數大于20，H2O的生成消耗了大量的氧，同時也抑制了CO2的生成。說明在初始的分子設計中C/N比和H/O比，對熱解反應路徑和最終產物分布起決定作用。

仔細觀察圖6各體系中N2和CO2分子數目，發現它們幾乎在同一時間達到最大，這與N2和CO2的生成途徑有關；反應早期的主要步驟是N―NO2的斷裂生成NO2，所以剛開始體系中NO2數量最多，隨后NO2進一步反應生成N2，而NO2中的O則與中間物CO反應生成CO2，因而N2和CO2達到最大值的時間一致。同時還表明，高溫下N―N兩原子主要以N2的形式存在，與Strachan等10,19的研究結果一致。此外，在反應過程中還觀察到一些數量很少基團出現，如：H、HO、CO、HNO3、N2O和HONO等基團，這些小分子基團相互聚合或反應生成其它產物；因其數量較少，故在本文沒有對其進行詳細討論。

由于圖6很難清晰展示中間體NO2和NO的變化，而中間體的變化影響著整個熱分解反應過程，故研究了前20 ps中間體NO2和NO隨時間的變化，如圖7所示。RDX-D2中的NO2峰值最高，RDX-D3中的NO2峰值最低，RDX和RDX-D2中的NO2峰值(＞13)略高于RDX-D1和RDX-D3的NO2峰值(〈10)。從圖7(a)中可以看出，NO2的數量在前2.5 ps內迅速達到峰值，表明四種物質首先發生斷裂的鍵是N―NO2鍵，這與圖3中研究表明高溫下N―NO2鍵的斷裂是各體系最主要的初始反應途徑相一致。隨后NO2的數量開始減少，這是由于中間產物NO2參與次級反應不斷被消耗，生成NO，NO3和HONO等中間產物。

圖6 四個體系的主要產物分子(N2,H2O,CO2,NO2,NO)數目隨時間變化關系Fig.6 Time evolution of key produced molecules(N2,H2O,CO2,NO2,NO)for the four materials

同時，NO也是很活躍的中間產物，從圖7(b)可以看出，NO和NO2的分布曲線具有相似的演化特征，均是先增大后減小的過程，但各體系NO達到峰值前的增長速率以及隨后的衰減速率相比NO2較小。在整個反應過程中，NO數量整體上明顯低于NO2，且NO分子數達到最大值的時間相對NO2滯后一些，這是因為后續反應消耗掉一部分NO2生成NO。NO的產生主要有三種途徑，(1)N―NO2→N―ONO，然后NO―NO單鍵斷裂生成NO；(2)氫原子轉移反應生成的HONO，隨后分解為HO和NO分子；(3)由NO2和NO3等相互反應生成NO。在20 ps后，NO2和NO的數量很少，幾乎被完全消耗。

表1 各體系的C/N和H/O比值及主產物CO2和H2O的生成量Table 1 C/N and H/O ratios of each system and the amount of the main products of CO2and H2O

3.4 大分子含碳團簇的形成

實驗研究表明，富碳高能材料在爆轟過程中能形成大分子含碳團簇30－34。大分子含碳團簇的出現是深入了解含能材料性質的重要方面，目前采用分子動力學模擬，已經開展含能材料熱解過程含碳團簇形成動力學機制的研究35。為了深入了解含碳團簇形成過程，首先模擬了各單胞體系高溫熱解過程最大含碳團簇C原子數隨時間的演變，如圖8(a,b,c)所示。本文中我們將C原子數大于8定義為含碳團簇，隨著溫度升高，RDX和RDX-D2單胞體系熱解過程中最大含碳團簇C原子數有所增加，在3000 K時含碳團簇最大碳原子數可達19。而在RDX-D1和RDX-D3單胞體系熱解過程中，均未生成含碳團簇，只存在小分子含碳碎片(C原子數小于等于3)，且在不同溫度下幾乎保持不變。由于大分子含碳團簇的生成受初始模型的影響較大，因此，為了正確考慮其形成過程，我們進一步采用超胞體系進行模擬，如圖8(d)所示。

圖7 四個體系熱解過程中NO2(a)和NO(b)的分子數量變化關系Fig.7 Quantities of NO2(a)and NO(b)for the thermal pyrolysis of four systems

圖8 不同溫度下各體系單胞(a－c)和超胞(d)熱解過程中最大含碳團簇所含C原子數隨時間變化關系Fig.8 Time evolution of C atomic number in the maximum molecule formed in the thermal pyrolysis of unit cell(a－c)and super cell(d)under different temperatures

對比圖8(c,d)單胞和多胞體系可以看出，RDX和RDX-D2超胞體系中生成的最大含碳團簇C原子數增加迅速，最終C原子數目分別達到約30和16個；RDX-D1和RDX-D3單胞和超胞體系的模擬結果相近，均未生成大分子含碳團簇，生成的小分子含碳碎片中C原子數均小于等于3?？紤]初始分子結構中，RDX和RDX-D2分子中原子比C/N≥0.5，而RDX-D1和RDX-D3分子中C/N〈0.5，說明：當分子中C/N≥0.5時，生成的最大含碳團簇C原子數與模型大小有關；而對C/N〈0.5的體系，含碳團簇的形成過程基本不受模型大小的影響。

因而，初始設計的分子中碳、氮含量比和反應物分子結構對RDX及其衍生物中最大含碳團簇的形成過程影響較大，對最大含碳團簇中C原子數目多少也影響顯著，且C/N比越大越容易形成大分子含碳團簇。根據上述模擬結果，在熱分解的初級反應中，環上脫落NO2、NO和CH3等自由基后，剩余部分發生聚合形成大分子含碳團簇，這與實驗和理論上研究硝胺類化合物的高溫熱解過程是一致的20,36。

4 結論

采用ReaxFF反應力場分子動力學方法，模擬研究了RDX及其衍生物單胞和超胞在不同溫度下的熱分解行為。研究表明，這四種晶體高溫熱解的第一步均為N―NO2鍵的斷裂生成NO2，然后發生取代基的脫去和環的斷裂反應。RDX在斷裂兩個硝基后就發生開環反應，而RDX的衍生物在斷裂兩個硝基后并沒有發生開環反應，而是先發生取代基的脫去反應再發生開環反應，這是由于各體系分子六元環上基團和側鏈基團的改變使得初始反應途徑不同。

通過對各體系熱分解過程主要產物和中間產物進行詳細分析，發現NO2和NO是裂解過程中十分重要的中間產物。NO2主要來自N―NO2中的N―N鍵的斷裂產生；NO主要來自N―NO2→N―ONO，然后NO―NO單鍵斷裂生成NO，中間產物HONO分解為HO和NO以及次級反應生成NO分子。熱分解后最穩定的氣相產物為N2、H2O和CO2，其中分子數最多的是N2(＞20)，主要來源于NO2和NO參與的次級反應中。各分子體系中的C/ N比和H/O比不同，使得主要分解產物H2O和CO2分子數目相差較大。C/N和H/O比越大的體系，熱解后生成的H2O分子數越多，H2O的生成消耗了大量的氧，故生成穩定產物CO2分子數越少，生成含碳團簇分子量越大。模型分子設計中C/N比和H/O比，對于最終熱解過程和產物分布起決定作用。

單胞體系模擬分析表明：溫度升高有利于RDX和RDX-D2單胞體系熱解中最大含碳團簇的形成，在3000 K時初步形成大分子含碳團簇(C原子數約為10)。進一步超胞體系模擬結果表明：相對于單胞模擬，RDX和RDX-D2超胞熱解過程中最大含碳團簇C原子數顯著增加(約為30和16)，初始模型大小和反應溫度對含碳團簇的生成影響明顯。而RDX和RDX-D2體系中的C/N比明顯高于RDX-D1和RDX-D3體系，對于后兩者的模擬，單胞和超胞模擬結果相近，均未生成含碳團簇，僅存在小分子含碳碎片(碳原子數小于等于3)。含碳團簇的大量生成，不利于含能材料能量的釋放，故分子設計中減小C/N比能很好地改善含能材料的熱解性能。

(1) Shaw,R.W.;Brill,T.B.;Thompson,D.L.Overviews of Recent Research on Energetic Materials;World Scientific:NJ,2005, Vol.16.

(2) Teipel,U.Energetic Materials:Particle Processing and Characterization;Wiley-VCH:Weinheim,Germany,2005.

(3) Chakraborty,D.;Muller,R.P.;Dasgupta,S.;Goddard,W.A.J. Phys.Chem.A 2000,104,2261.doi:10.1021/jp9936953

(5) Schweigert,I.V.J.Phys.Chem.A 2015,119,2747. doi:10.1021/jp510034p

(6) Irikura,K.K.J.Phys.Chem.A 2013,117,2233.doi:10.1021/ jp310247z

(7) Behrens,R.,Jr.;Bulusu,S.J.Phys.Chem.1992,96,8877. doi:10.1063/1.462245

(8) Van Duin,A.C.;Dasgupta,S.;Lorant,F.;Goddard,W.A.J. Phys.Chem.A 2001,105,9396.doi:10.1021/jp004368u

(9) Cohen,R.;Zeiri,Y.;Wurzberg,E.;Kosloff,R.J.Phys.Chem.A 2007,111,11074.doi:10.1021/jp072121s

(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) 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) 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) 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) Liu,H.;Dong,X.;He,Y.H.Acta Phys.-Chim.Sin.2014,30, 232.[劉 海,董 曉,何遠航.物理化學學報,2014,30,232.] doi:10.3866/PKU.WHXB201312101

(15) Qian,H.J.;van Duin,A.C.;Morokuma,K.;Irle,S.J.Chem. Theory Comput.2011,7,2040.doi:10.1021/ct200197v

(16) Liu,L.;Bai,C.;Sun,H.;Goddard,W.A.,III.J.Phys.Chem.A 2011,115,4941.doi:10.1021/jp110435p

(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) 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) 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) 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) Hakey,P.;Ouellette,W.;Zubieta,J.;Korter,T.Acta Crystallographica Section E:Structure Reports Online 2008,64 (8),o1428.doi:10.1107/S1600536808019727

(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) 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) Botcher,T.R.;Wight,C.A.J.Phys.Chem.1994,98,5441. doi:10.1021/j100072a009

(25) Oyumi,Y.;Brill,T.B.Combust.Flame 1985,62,213. doi:10.1016/0010-2180(85)90147-6

(26)Kohno,Y.;Ueda,K.;Imamura,A.J.Phys.Chem.1996,100, 4701.doi:10.1021/jp9503223

(27)Prabhakaran,K.V.;Bhide,N.M.;Kurian,E.M.Thermochim. Acta 1995,249,249.doi:10.1016/0040-6031(95)90704-1

(28) Oxley,J.C.;Kooh,A.B.;Szekeres,R.;Zheng,W.J.Phys. Chem.1994,98,7004.doi:10.1021/j100079a019

(29) Brill,T.B.;Oyumi,Y.J.Phys.Chem.1986,90,6848. doi:10.1021/j100284a028

(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)Anisichkin,V.F.;Dolgushin,D.S.;Petrov,E.A.Combust. Explos.Shock Waves.1995,31,106.doi:10.1007/BF00755966

(32) Titov,V.M.;Anisichkin,V.F.;Mal′kov,I.Y.Combust.Explos. Shock Waves.1989,25,372.doi:10.1007/BF00788819

(33) Baidakova,M.J.Phys.D:Appl.Phys.2007,40,6300. doi:10.1088/0022-3727/40/20/S14

(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) Chevrot,G.;Sollier,A.;Pineau,N.J.Chem.Phys.2012,136, 084506.doi:10.1063/1.3686750

(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

Pyrolysis of RDX and Its Derivatives via Reactive Molecular Dynamics Simulations

PENG Li-Juan1YAO Qian2WANG Jing-Bo3LI Ze-Rong2ZHU Quan3,*LI Xiang-Yuan3,*
(1School of Aeronautics&Astronautics,Sichuan University,Chengdu 610065,P.R.China; 2College of Chemistry,Sichuan University,Chengdu 610064,P.R.China; 3School of Chemical Engineering,Sichuan University,Chengdu 610065,P.R.China)

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―NO2homolytic cleavage to produce NO2;however,the subsequent reaction mechanism is completely different due to the different sixmembered rings and side chain groups.In these four model systems,NO2and NO molecules are common intermediates in the thermal decomposition process,which eventually transform to N2in the following reactions. The most stable products are N2,H2O,and CO2in the thermal process,among which N2has the maximum molecular number(more than 20).The numbers of cracked H2O and CO2molecules 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.

RDX and its derivatives;Thermal decomposition;ReaxFF;Carbon clusters;Molecular dynamics

O643

Xu,J.;Zhao,J.;Sun,L.Mol.Simul.2008,34,961.

10.1080/ 08927020802162892

doi:10.3866/PKU.WHXB201701161

Received:October 22,2016;Revised:January 16,2017;Published online:January 16,2017.

*Corresponding authors.ZHU Quan,Email:qzhu@scu.edu.cn;Tel:+86-28-85402951.LI Xiang-Yuan,Email:xyli@scu.edu.cn.

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).國家自然科學基金(91441132)和教育部新世紀優秀人才支持計劃(NCET-13-0398)資助項目