• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

      2016-05-09 02:42:04苗瑞珍劉偉帥康志鵬楊潞霞荊興斌付一政劉亞青
      含能材料 2016年2期
      關(guān)鍵詞:力場(chǎng)共晶勢(shì)能

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

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

      猜你喜歡
      力場(chǎng)共晶勢(shì)能
      基于物理模型的BaZrO3鈣鈦礦機(jī)器學(xué)習(xí)力場(chǎng)
      “動(dòng)能和勢(shì)能”知識(shí)鞏固
      力場(chǎng)防護(hù)罩:并非只存在于科幻故事中
      作 品:景觀設(shè)計(jì)
      ——《勢(shì)能》
      文化縱橫(2022年3期)2022-09-07 11:43:18
      “動(dòng)能和勢(shì)能”知識(shí)鞏固
      調(diào)性的結(jié)構(gòu)力場(chǎng)、意義表征與聽(tīng)覺(jué)感性先驗(yàn)問(wèn)題——以貝多芬《合唱幻想曲》為例
      “動(dòng)能和勢(shì)能”隨堂練
      Cr12Mo1V1鍛制扁鋼的共晶碳化物研究
      模具制造(2019年3期)2019-06-06 02:11:04
      脫氧核糖核酸柔性的分子動(dòng)力學(xué)模擬:Amber bsc1和bsc0力場(chǎng)的對(duì)比研究?
      《含能材料》“含能共晶”征稿
      含能材料(2017年1期)2017-03-04 15:46:20
      岚皋县| 轮台县| 娱乐| 景洪市| 巴彦县| 敖汉旗| 贵港市| 彭州市| 高平市| 龙门县| 连山| 新巴尔虎右旗| 奉节县| 龙口市| 璧山县| 皋兰县| 海兴县| 长宁区| 岳阳县| 弥渡县| 兴海县| 中卫市| 游戏| 黄骅市| 浮山县| 家居| 正镶白旗| 永嘉县| 慈利县| 阳新县| 永城市| 上饶县| 泰顺县| 繁峙县| 台北县| 云南省| 云安县| 亳州市| 武清区| 青神县| 墨竹工卡县|