CL-20/DNB共晶高溫?zé)峤獾腞eaxFF反應(yīng)分子動(dòng)力學(xué)模擬

苗瑞珍, 劉偉帥, 王 建, 康志鵬, 楊潞霞, 荊興斌, 付一政, 劉亞青

(1. 中北大學(xué)山西省高分子復(fù)合材料工程技術(shù)研究中心, 山西 太原 030051; 2. 中北大學(xué)材料科學(xué)與工程學(xué)院, 山西 太原 030051; 3. 山西大學(xué)商務(wù)學(xué)院 信息學(xué)院, 山西 太原 030051)

1 引 言

近年來(lái),由Nielsen A T[1]首次合成的六硝基六氮雜異伍茲烷(HNIW,CL-20)由于具有高的能量密度比以及優(yōu)異的爆速和爆壓性能,引起研究者廣泛的關(guān)注。然而,較高的感度制約了它在武器方面的應(yīng)用,如何平衡能量與安全性之間的矛盾已成為研究難點(diǎn)。有研究表明,通過(guò)制備共晶可降低CL-20感度,如CL-20/三硝基甲苯(TNT)[2-3],CL-20/環(huán)四亞甲基四硝胺(HMX)[4]和CL-20/苯并三氧化呋咱(BTF)[5],其中王玉平[6]合成的摩爾比為1∶1的CL-20/1,3-二硝基苯(DNB)共晶兼具了CL-20高能和DNB鈍感廉價(jià)的優(yōu)勢(shì),具有優(yōu)異的使用價(jià)值。

目前,關(guān)于CL-20熱解機(jī)理的研究比較多,Isayev[7]運(yùn)用從頭算分子動(dòng)力學(xué)方法分別研究了單分子和凝聚態(tài)CL-20的熱解機(jī)理,結(jié)果表明二者熱解路徑并不完全相同,其中NO2的生成勢(shì)壘決定了整個(gè)反應(yīng)的熱解動(dòng)力學(xué)。張力等[8]運(yùn)用ReaxFF力場(chǎng)對(duì)不同溫度和密度條件下CL-20單分子和超晶胞的初始熱解路徑進(jìn)行了研究,結(jié)果表明不同密度和溫度條件只會(huì)對(duì)熱解的速率有影響,熱解機(jī)理不受影響。Patil和Brill[9]通過(guò)恒溫?zé)崾е胤ê虵TIR光譜法研究了CL-20熱解動(dòng)力學(xué),結(jié)果表明低溫時(shí)N—O2鍵的斷裂反應(yīng)決定其熱解速率,其他氣體產(chǎn)物的濃度比是溫度的函數(shù)。由于DNB比較鈍感,一般不單獨(dú)使用,所以對(duì)其熱解機(jī)理的研究比較少,對(duì)CL-20/DNB共晶熱解機(jī)理的研究更少。含能材料主要通過(guò)在極端條件下發(fā)生復(fù)雜的化學(xué)反應(yīng)釋放能量,反應(yīng)速度非?？?實(shí)驗(yàn)研究物質(zhì)的熱解機(jī)理存在難度。反應(yīng)分子動(dòng)力學(xué)模擬可以描述原子水平和飛秒尺度上化學(xué)反應(yīng)的詳細(xì)信息,尤其是近期發(fā)展起來(lái)的ReaxFF反應(yīng)分子動(dòng)力學(xué)方法,計(jì)算精度與量子力學(xué)相當(dāng),速度快,可以處理上百萬(wàn)級(jí)原子的體系,時(shí)間尺度可達(dá)飛秒級(jí)別,有助于在原子水平上認(rèn)識(shí)含能材料在極端條件下的反應(yīng)機(jī)理,得到其能量釋放規(guī)律[10]。此力場(chǎng)在很多含能材料的沖擊起爆、沖擊點(diǎn)火和爆轟反應(yīng)等反應(yīng)機(jī)理的研究方面得到廣泛應(yīng)用[11-23]。本研究基于ReaxFF力場(chǎng),運(yùn)用LAMMPS (large-scale atomic/molecular massively parallel simulator)軟件[24]研究了2000,2500 K以及3000 K條件下CL-20/DNB的熱解過(guò)程,通過(guò)分析產(chǎn)物隨模擬時(shí)間的變化得到其熱解引發(fā)機(jī)理,以期為新材料的設(shè)計(jì)、開(kāi)發(fā)、合成、運(yùn)輸、起爆、安全性和風(fēng)險(xiǎn)評(píng)估等實(shí)際問(wèn)題提供理論指導(dǎo)。

2 ReaxFF力場(chǎng)簡(jiǎn)介與模擬細(xì)節(jié)

2.1 ReaxFF力場(chǎng)簡(jiǎn)介

ReaxFF[10]反應(yīng)力場(chǎng)基于第一性原理,由Adri van Duin和William A. Goddard提出,力場(chǎng)勢(shì)能函數(shù)以原子間的鍵級(jí)BOij(Bond Order)為基礎(chǔ)。ReaxFF力場(chǎng)的勢(shì)能表達(dá)式如下,其中主要分為鍵、角、二面角、共軛、庫(kù)侖、范德華及調(diào)整項(xiàng)等(單位: kJ·mol-1)。

Esystem=Ebond+Elp+Eover+Eunder+Eval+Epen+Ecoa+

EC2+Etriple+Etors+Econj+EH-bond+Evdwaals+ECoulomb

(1)

該系統(tǒng)中各能量項(xiàng)的具體意義為:Ebond表示鍵能:Elp是和孤對(duì)電子相關(guān)的能量補(bǔ)償項(xiàng);Eover和Eunder分別是過(guò)配位的能量矯正項(xiàng)和低配位的能量矯正項(xiàng);Eval、Ecoa、Epen是價(jià)角能量項(xiàng)；EC2是C2修正項(xiàng);Etriple是三鍵修正項(xiàng);Etors和Econj是四體作用項(xiàng);EH-bond是氫鍵作用項(xiàng);Evdwaals和ECoulomb分別是范德華和庫(kù)倫作用項(xiàng)。其中,Ebond、Eval、Etors是依賴(lài)于鍵級(jí)的價(jià)鍵相互作用。

2.2 模擬

CL-20/DNB共晶X-ray衍射數(shù)據(jù)來(lái)源于劍橋晶體學(xué)數(shù)據(jù)庫(kù)(CCDC NO.940129),由CL-20(圖1a)和DNB(圖1b)按摩爾比為1∶1組成,屬于正交晶系,Pbca空間群,晶格參數(shù)a=9.4703 ?,b=13.4589 ?,c=33.620 ?。首先運(yùn)用MS 6.0(Materials Studio 6.0)軟件[25]構(gòu)建CL-20/DNB的單晶胞模型,如圖2a所示,為了避免溫度和升溫速率波動(dòng)對(duì)反應(yīng)的影響,構(gòu)建了2×2×1的超晶胞模型,共1664個(gè)原子,如圖2b所示。

a. CL-20b. DNB

圖1CL-20和DNB的分子結(jié)構(gòu)

Fig.1Molecular structures for CL-20 and DNB

運(yùn)用Lammps軟件[24]對(duì)CL-20/DNB超晶胞進(jìn)行優(yōu)化以獲得合理的初始構(gòu)型,在300 K下進(jìn)行500 ps等溫等容分子動(dòng)力學(xué)模擬(Canonical ensemble,NVT),接著采用(iso-thermal-isobaric,NPT)系綜對(duì)其壓力進(jìn)行10 ps弛豫,壓力為0 GPa。為了驗(yàn)證ReaxFF力場(chǎng)的實(shí)用性,將CL-20/DNB共晶在常溫下的晶胞參數(shù)與相應(yīng)實(shí)驗(yàn)值[6]進(jìn)行了對(duì)比,結(jié)果如表1所示,由表1可得,ReaxFF力場(chǎng)作用下的模擬值與實(shí)驗(yàn)值具有較好的一致性,說(shuō)明ReaxFF力場(chǎng)可用于CL-20/DNB共晶的研究。最后分別在2000,2500 K和3000 K下進(jìn)行50 ps的NVT分子動(dòng)力學(xué)計(jì)算,選擇上述高溫條件主要是為了模擬常規(guī)高能量密度材料爆轟中柴普曼-柔格(CJ)點(diǎn)的溫度波動(dòng),得到不同溫度對(duì)共晶熱解反應(yīng)的影響,并且能夠兼顧計(jì)算機(jī)的運(yùn)算能力與所用時(shí)間(在不影響反應(yīng)路徑的條件下提高溫度可以加快反應(yīng),節(jié)省時(shí)間)。采用Berendsen方法[26]對(duì)溫度和壓力進(jìn)行控制,耦合參數(shù)為100 fs,使其在設(shè)定的范圍內(nèi)波動(dòng),采用周期性邊界條件,步長(zhǎng)為0.1 fs,每隔50 fs記錄一次原子軌跡以及鍵級(jí)信息,鍵級(jí)為0.3。

a. unit cell

b. supercell

圖2CL-20/DNB共晶單胞結(jié)構(gòu)和超晶胞結(jié)構(gòu)

Fig.2The unit cell and supercell structures of CL-20/DNB cocrystal

表1常溫下CL-20/DNB共晶晶胞參數(shù)模擬值與實(shí)驗(yàn)值

Table1The lattice parameters of the cocrystal of CL-20/DNB calculated by ReaxFF force field compared with the experiment value

latticeparameterexperimentvalue[6]ReaxFFerror/%a/? 9．4703 9．5606 0．95b/? 13．4589 13．5872 0．95c/? 33．620 33．9406 0．95density/g·cm-3 1．880 1．82679 -2．8

3 結(jié)果與討論

3.1 溫度對(duì)勢(shì)能和總物種數(shù)的影響

圖3給出了不同溫度下體系勢(shì)能和總物種數(shù)隨時(shí)間的變化情況。從圖3a看不出明顯的平衡和誘導(dǎo)期,體系迅速吸收熱量發(fā)生初級(jí)熱解反應(yīng),隨著反應(yīng)的進(jìn)行,熱解生成一系列中間及最終產(chǎn)物,并放出大量的熱,因而勢(shì)能發(fā)生明顯的下降; 溫度越高,系統(tǒng)達(dá)到反應(yīng)平衡所需時(shí)間越短,熱解反應(yīng)越完全,生成的產(chǎn)物越多,勢(shì)能下降也越多。

由圖3b可知,初始晶胞內(nèi)只有CL-20和DNB兩種物質(zhì),幾乎在0.01 ps時(shí)就有產(chǎn)物生成,熱解溫度為2500 K和3000 K時(shí)物種總量分別在25 ps和6.9 ps內(nèi)達(dá)到了最大值,隨著反應(yīng)的進(jìn)一步發(fā)生呈現(xiàn)一定程度的衰減。而在較低的溫度(2000 K)時(shí),前5 ps內(nèi)物種迅速增加到50種,而后緩慢增加。由此可知,溫度越高,共晶分解速率越快,初期產(chǎn)生的新物種數(shù)量越多。

3.2 共晶熱解的引發(fā)反應(yīng)機(jī)理

圖4為2000,2500 K和3000 K下CL-20及DNB的消耗情況和主要產(chǎn)物的分布。由圖4a可知,三個(gè)溫度下CL-20在0.05 ps時(shí)均分解完畢,而DNB消耗較慢,隨著溫度的升高其消耗速率顯著加快。由圖4b,圖4c,圖4d可知,熱解產(chǎn)物主要有NO2、N2、NO、H2O、HNO3、HONO、HON以及CO2等。每個(gè)溫度下最早生成的產(chǎn)物均為NO2,其含量在0.05 ps時(shí)就達(dá)到了最大值,它是由CL-20中N—O2鍵和DNB中C—NO2鍵斷裂共同作用的結(jié)果,—NO2鍵的斷裂是CL-20/DNB熱解的主要引發(fā)反應(yīng)。對(duì)于CL-20這種典型的環(huán)狀硝胺類(lèi)化合物來(lái)說(shuō)N—NO2鍵的斷裂是其熱解的特征引發(fā)反應(yīng)[27],這與文獻(xiàn)[28-30]中的模擬結(jié)果一致。此外,也與Naik和Gore[31]得出N—NO2鍵的斷裂是CL-20熱解過(guò)程中生成其他小分子產(chǎn)物的前提這一實(shí)驗(yàn)結(jié)論相互吻合。DNB中的C—NO2鍵斷裂生成NO2是其主要引發(fā)反應(yīng),因?yàn)镃—NO2鍵是硝基芳香化合物中最弱的鍵[11]。

a. potential energy

b. total species

圖3不同溫度下體系勢(shì)能以及反應(yīng)過(guò)程中總物種數(shù)量隨模擬時(shí)間的變化

Fig.3Time evolution of potential energy and total species at various temperatures

a. CL-20 and DNBb. 2000 K

c. 2500 Kd. 3000 K

圖4不同溫度下CL-20及DNB的消耗情況和主要產(chǎn)物隨模擬時(shí)間的變化

Fig.4Time evolution of consumption of CL-20 and DNB and main products at various temperatures

3.3 一些主要產(chǎn)物的生成機(jī)理以及溫度的影響

3.2節(jié)的分析表明,NO2是共晶熱解引發(fā)反應(yīng)所產(chǎn)生,不同溫度下,各主要產(chǎn)物隨時(shí)間的變化如圖5所示,由圖5a可知,NO2數(shù)量在達(dá)到最大值后會(huì)迅速減小,這是由于次級(jí)反應(yīng)消耗掉大量的NO2,生成NO、HONO和HNO3等物質(zhì)。然而在5 ps(2000 K),6 ps(2500 K)和10 ps(3000 K)左右時(shí)NO2數(shù)量略有增加,這主要是由于DNB在此時(shí)開(kāi)始快速分解引起。

NO主要的生成途徑有以下幾種: (1)兩個(gè)NO2自由基結(jié)合生成N2O4,接著繼續(xù)熱解生成NO3和NO; (2)對(duì)于DNB來(lái)說(shuō),由硝基重排形成亞硝酸鹽結(jié)構(gòu)(C—NO2→C—ONO),隨后O—NO鍵斷裂也會(huì)生成NO,這與具有相似結(jié)構(gòu)的TNT生成NO的機(jī)理相同[11]; (3)由CL-20中N—N鍵和DNB中C—N鍵斷裂生成的硝基自由基重排布形成亞硝酸鹽結(jié)構(gòu),通過(guò)奪取一個(gè)質(zhì)子生成HONO,接著分子中N—O鍵斷裂生成NO,這與Isayev[7]等對(duì)凝聚態(tài)CL-20熱解機(jī)理的研究所得結(jié)果相吻合; (4)NO2與含碳分子片段生成NO。NO和NO2的分布非常相似,也是先快速增大然后減少,大多參與反應(yīng)生成N2[11],且溫度越高,NO2和NO的消耗速率越快,生成N2的量越多,此過(guò)程中N2的量一直在不斷增多,如圖5c所示,也進(jìn)一步驗(yàn)證了消耗掉的大部分NO2和NO生成了N2。

a. NO2b. NO

c. N2d. HONO

e. H2Of. CO2

圖5主要產(chǎn)物在不同溫度下的對(duì)比分布曲線

Fig.5Comparative distribution curves of main products at various temperatures

HONO的主要來(lái)源是CL-20和DNB中—NO2鍵斷裂生成的硝基自由基重排布形成亞硝酸鹽結(jié)構(gòu),通過(guò)奪取一個(gè)質(zhì)子所生成,如圖5d所示隨著反應(yīng)的進(jìn)行HONO的數(shù)量先增多后減少,主要是由于HONO熱解消耗引起。

由圖5e可以看出H2O在整個(gè)熱解過(guò)程中數(shù)量先增加然后基本保持不變,溫度越高生成速率越快。水的生成路徑有: (1)兩個(gè)分子間發(fā)生質(zhì)子的轉(zhuǎn)移使硝基生成—NOOH,然后繼續(xù)奪取一個(gè)質(zhì)子生成—NOOH—H,接著脫去一分子的H2O; (2)NO2自由基重排生成的ONO自由基取鄰近分子中的質(zhì)子生成HONO,隨著反應(yīng)的進(jìn)行,HONO會(huì)進(jìn)一步熱解產(chǎn)生HO自由基,HO自由基繼續(xù)奪取質(zhì)子生成H2O,這與反應(yīng)過(guò)程中HONO先快速積累而后減少是一致的; (3)HNO3的熱解也可以產(chǎn)生H2O。

由圖5f可以看出CO2數(shù)量隨著反應(yīng)的進(jìn)行不斷增多,且溫度越高,CO2生成越早,在2000 K時(shí)一直到16.85 ps才有CO2生成,而在3000 K時(shí)在2.85 ps就有CO2產(chǎn)生,這是由于C—C鍵斷裂所需要的能量較多,一般是在環(huán)上基團(tuán)脫去或者相互反應(yīng)后,主環(huán)上的C—C鍵才開(kāi)始斷裂,進(jìn)而反應(yīng)生成CO2。溫度為3000 K時(shí),在40 ps后CO2數(shù)量基本保持不變,這與2000 K和2500 K的現(xiàn)象不同,其可能的原因在于,一方面高溫下隨著反應(yīng)的進(jìn)行,生成的一些小分子片段運(yùn)動(dòng)速率加快,發(fā)生碰撞繼續(xù)聚合生成大分子含碳團(tuán)簇的幾率變大,使得CO2數(shù)量減少; 另一方面隨著能量的積累,大的團(tuán)簇又發(fā)生熱解,最終使得總的數(shù)量基本保持不變。通過(guò)對(duì)產(chǎn)物的分析可知,產(chǎn)物中確實(shí)有少部分大的含碳團(tuán)簇,這是含能材料爆轟過(guò)程中的常見(jiàn)現(xiàn)象,一些富碳炸藥尤為明顯[32-35]。

4 結(jié) 論

(1)隨著溫度的升高,系統(tǒng)達(dá)到平衡時(shí)間越短,熱解反應(yīng)越完全,勢(shì)能下降越多,最終總產(chǎn)物生成越多。三個(gè)溫度下CL-20均先于DNB分解完畢,升高溫度DNB的消耗速率以及產(chǎn)物的生成顯著加快。CL-20/DNB共晶熱解的主要終產(chǎn)物為NO2、NO、N2、H2O、HNO3、HON、HONO、CO2等,還有一些重要的中間產(chǎn)物如N2O5、N2O4、OH自由基、N2O等。

(2)熱解反應(yīng)過(guò)程中,CL-20分子中N—NO2鍵和DNB分子中C—NO2鍵的斷裂生成NO2是最主要的引發(fā)反應(yīng),接著兩個(gè)NO2分子結(jié)合或者重排生成ONO自由基,進(jìn)一步生成HONO、HON、NO、NO3、H2O等,其中大部分的NO2和NO繼續(xù)反應(yīng)生成N2。3000 K時(shí)對(duì)產(chǎn)物的分析表明在模擬后期會(huì)發(fā)生小分子片段的聚合生成大分子碳團(tuán)簇的現(xiàn)象。

(3)本研究的熱解產(chǎn)物和引發(fā)機(jī)理與前人的一些結(jié)論相互吻合,說(shuō)明了ReaxFF力場(chǎng)的實(shí)用性與正確性。

參考文獻(xiàn):

[1] Nielsen A T, Chan M L,Kraeutle C K, et al. Polynitropolyazacaged explosives, part 7, NWC TP7200[M]. China Lake: Navel Weapons Center, 1989.

[2] 楊宗偉, 張艷麗, 李洪珍, 等. CL-20/TNT 共晶炸藥的制備, 結(jié)構(gòu)與性能[J]. 含能材料, 2013, 20(6): 674-679.

YANG Zong-wei, ZHANG Yan-li, LI Hong-zhen, et al. Preparation, structure and properties of CL-20/TNT cocrystal[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2013, 20(6): 674-679.

[3] Bolton O,Matzger A J. Improved stability and smart-material functionality realized in an energetic cocrystal[J].AngewandteChemieInternationalEdition, 2011, 50(38): 8960-8963.

[4] Bolton O,Simke L R, Pagoria P F, et al. High power explosive with good sensitivity: a 2:1 cocrystal of CL-20: HMX[J].CrystalGrowth&Design, 2012, 12(9): 4311-4314.

[5] YANG Zong-wei, LI Hong-zhen, ZHOU Xiao-qing, et al. Characterization and properties of a novel energetic-energetic cocrystal explosive composed of HNIW and BTF[J].CrystalGrowth&Design, 2012, 12(11): 5155-5158.

[6] WANG Yu-ping, YANG Zong-wei, LI Hong-zhen, et al. A novel cocrystal explosive of HNIW with good comprehensive properties[J].Propellants,Explosives,Pyrotechnics, 2013, 35: 1-7.

[7] Olexandr I, Gorb L, Qasim M, et al. Ab initio molecular dynamics study on the initial chemical events in nitramines: thermal decomposition of CL-20[J].TheJournalofPhysicalChemistryB, 2008, 112(35): 11005-11013.

[8] 張力, 陳朗, 王晨等. CL-20初始熱分解反應(yīng)機(jī)理的分子動(dòng)力學(xué)計(jì)算[J]. 火炸藥學(xué)報(bào), 2012 (4): 5-9.

ZHANG Li, CHEN Lang, WANG Chen, et al. Mechanism of the initial thermal decomposition of CL-20 via molecular dynamic simulation[J].ChineseJournalofExplosives&Propellants, 2012 (4): 5-9.

[9] Patil D G, Brill T B. Thermal decomposition of energetic materials 53 kinetics and mechanism of thermolysis of hexanitrohexazaisowurtzitane[J].CombustionandFlame, 1991, 87(2): 145-151.

[10] VanDuin A C T, Dasgupta S, Lorant F, et al. ReaxFF: a reactive force field for hydrocarbons[J].TheJournalofPhysicalChemistryA, 2001, 105(41): 9396-9409.

[11] LIU H, DONG X, HE Y H. Reactive molecular dynamics simulations of carbon-containing clusters formation during pyrolysis of TNT[J].ActaPhysico-ChimicaSinica, 2014, 30(2): 232-240.

[12] Alejandro, Strachan, Edward M, Kober, Adri C T, van Duin, et al. Thermal decomposition of RDX from reactive molecular dynamics[J].JournalofChemicalPhysics, 2005, 122(5):054502-054502-10.

[13] Strachan A, vanDuin A C T, Chakraborty D, et al. Shock waves in high-energy materials: the initial chemical events in nitramine RDX[J].PhysicalReviewLetters, 2003, 91(9): 098301.

[14] Zhou T, Liu L, Goddard III W A, et al. ReaxFF reactive molecular dynamics on silicon pentaerythritol tetranitrate crystal validates the mechanism for the colossal sensitivity[J].PhysicalChemistryChemicalPhysics, 2014, 16(43): 23779-23791.

[15] 周婷婷, 石一丁, 黃風(fēng)雷. 高壓下β-HMX 熱分解機(jī)理的 ReaxFF 反應(yīng)分子動(dòng)力學(xué)模擬[J]. 物理化學(xué)學(xué)報(bào), 2012, 28(11): 2605-2615.

ZHOU Ting-ting, SHI Yi-ding, HUANG Feng-lei. Thermaldecomposition mechanism ofβ-HMX under high pressuresviaReaxFF reactive molecular dynamics simulations[J].ActaPhysico-ChimicaSinica, 2012, 28(11): 2605-2615.

[16] Zhou TT, Huang F L. Effects of defects on thermal decomposition of HMX via ReaxFF molecular dynamics simulations[J].TheJournalofPhysicalChemistryB, 2010, 115(2): 278-287.

[17] Zhang L,Zybin S V, van Duin A C T, et al. Carbon cluster formation during thermal decomposition of octahydro-1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazocine and 1, 3, 5-triamino-2, 4, 6-trinitrobenzene high explosives from ReaxFF reactive molecular dynamics simulations[J].TheJournalofPhysicalChemistryA, 2009, 113(40): 10619-10640.

[18] vanDuin A C T, Zeiri Y, Dubnikova F, et al. Atomistic-scale simulations of the initial chemical events in the thermal initiation of triacetonetriperoxide[J].JournaloftheAmericanChemicalSociety, 2005, 127(31): 11053-11062.

[19] Han S, vanDuin A C T, Goddard III W A, et al. Thermal decomposition of condensed-phase nitromethane from molecular dynamics from ReaxFF reactive dynamics[J].TheJournalofPhysicalChemistryB, 2011, 115(20): 6534-6540.

[20] Rom N,Zybin S V, van Duin A C T, et al. Density-dependent liquid nitromethane decomposition: molecular dynamics simulations based on ReaxFF[J].TheJournalofPhysicalChemistryA, 2011, 115(36): 10181-10202.

[21] 張力, 陳朗. 高壓下固相硝基甲烷分解的分子動(dòng)力學(xué)計(jì)算[J]. 物理學(xué)報(bào), 2013, 62(13): 138201-138201.

ZHANG Li, CHEN Lang. The effect of pressure on thermal decomposition of solid nitromethane via MD simulation[J].ActaPhysicalSinica, 2013, 62(13): 138201-138201.

[22] Chenoweth K, vanDuin A C T, Dasgupta S, et al. Initiation mechanisms and kinetics of pyrolysis and combustion of JP-10 hydrocarbon jet fuel[J].TheJournalofPhysicalChemistryA, 2009, 113(9): 1740-1746.

[23] Guo F, F G, Guo F , et al. ReaxFF molecular dynamics study of initial mechanism of JP-10 combustion[J].CombustionScienceandTechnology, 2012, 184:1233-1243.

[24] Plimpton S J, Large-scale atomic/molecular massively parallel simulator[CP]. Sandia National Laboratories: 2003.

[25] Materials Studio 6.0[CP]. Accelrys Software Inc, USA: 2013.

[26] Berendsen H J C, Postma J P M, van Gunsteren W F, et al. Molecular dynamics with coupling to an external bath[J].TheJournalofChemicalPhysics, 1984, 81(8): 3684-3690.

[27] Oxley J C, Kooh A B, Szekeres R, et al. Mechanisms of nitramine thermolysis[J].TheJournalofPhysicalChemistry, 1994, 98(28): 7004-7008.

[28] 張力, 陳朗, 王晨, 等. 不同晶型 CL-20 熱分解反應(yīng)機(jī)理計(jì)算分析[J]. 爆炸與沖擊, 2014, 34(2): 188-194.

ZHANG Li, CHEN Lang, WANG Chen, et al. Molecular dynamic simulation on thermal decomposition mechanism of CL-20 with different polymorphs[J].ExplosionandShockWaves, 2014, 34(2): 188-194.

[29] 張力, 陳朗, 王晨, 等. 水分子對(duì)α相CL-20熱分解機(jī)理影響的分子動(dòng)力學(xué)研究[J]. 物理化學(xué)學(xué)報(bào), 2013, 29(6): 1145-1153.

ZHANG Li, CHEN Lang, WANG Chen, et al. Molecular dynamics study of the effect of H2O on the thermal decomposition ofαphase CL-20[J].ActaPhysico-ChimicaSinica, 2013, 29(6): 1145-1153..

[30] 劉海, 李啟楷, 何遠(yuǎn)航. CL20-TNT共晶高溫?zé)峤獾腞eaxFF/lg反應(yīng)力場(chǎng)分子動(dòng)力學(xué)模擬[J]. 物理學(xué)報(bào), 2013, 62(20): 208202.

LIU Hai, LI Qi-kai, HE Yuan-hang.Pyrolysis of CL20-TNT cocrystal from ReaxFF/lg reactive molecular dynamics simulations[J].ActaPhysicalSinica, 2013, 62(20): 208202.

[31] Naik N H, Gore G M, Gandhe B R, et al. Studies on thermal decomposition mechanism of CL-20 by pyrolysis gas chromatography-mass spectrometry (Py-GC/MS)[J].JournalofHazardousMaterials, 2008, 159(2): 630-635.

[32] Mironov E V, Petrov E A, Korets A Y. From analysis of the structure of ultrafine diamond to the problem of its formation kinetics[J].Combustion,ExplosionandShockWaves, 2004, 40(4): 473-476.

[33] Chevrot G, Sollier A, Pineau N. Molecular dynamics and kinetic study of carbon coagulation in the release wave of detonation products[J].TheJournalofChemicalPhysics, 2012, 136(8): 084506.

[34] Ten K A, Aulchenko V M, Lukjanchikov L A, et al. Application of introduced nano-diamonds for the study of carbon condensation during detonation of condensed explosives[J].NuclearInstrumentsandMethodsinPhysicsResearchSectionA:Accelerators,Spectrometers,DetectorsandAssociatedEquipment, 2009, 603(1): 102-104.

[35] Krüger A, Kataoka F, Ozawa M, et al. Unusually tight aggregation in detonation nanodiamond: identification and disintegration[J].Carbon, 2005, 43(8): 1722-1730.