武宗凱，舒遠杰，劉　寧，丁小勇，吳敏杰，王　可，王　波，盧瑩瑩

(西安近代化學研究所， 陜西 西安 710065)

?

CL-20/FOX-7共晶的分子動力學模擬

武宗凱，舒遠杰，劉寧，丁小勇，吳敏杰，王可，王波，盧瑩瑩

(西安近代化學研究所， 陜西 西安 710065)

摘要：依據(jù)共晶形成的氫鍵規(guī)則，搭建了7種CL-20/1,1-二氨基-2,2-二硝基乙烯(FOX-7)共晶模型; 采用分子動力學(MD)模擬研究了CL-20/FOX-7共晶形成的可能性;通過徑向分布函數(shù)(RDF)考察了共晶模型內(nèi)分子的相互作用力;采用X-射線粉末衍射(XRD)模擬分析了共晶模型與純組分間衍射峰的區(qū)別。結(jié)果表明，F(xiàn)OX-7和CL-20分子間存在較強的氫鍵和范德華力，7 種共晶模型的結(jié)合能大小順序為:Eb(1 0 1)>Eb(1 1 -1)>Eb(隨機晶面)>Eb(0 1 1)>Eb(0 0 2)>Eb(1 1 0)>Eb(1 0 -1);7種共晶模型的分子間作用力以FOX-7取代CL-20(1 0 1)、(1 1 -1)晶面較強；7種共晶模型的XRD衍射峰相較純組分CL-20或FOX-7區(qū)別較大。據(jù)此推測在制備CL-20/FOX-7共晶過程中，F(xiàn)OX-7取代CL-20(1 0 1)晶面的共晶模型易于形成。

關(guān)鍵詞：共晶；分子動力學；六硝基六氮雜異伍茲烷；CL-20；1,1-二氨基-2,2-二硝基乙烯；FOX-7

引 言

隨著軍事科學技術(shù)的不斷發(fā)展，現(xiàn)代戰(zhàn)爭和新型武器對含能材料提出了高能量、高安全性的要求。新一代高能鈍感含能材料成為國內(nèi)外含能材料領(lǐng)域研究的熱點[1]。隨著共晶技術(shù)在含能材料領(lǐng)域的應(yīng)用，共晶成為研究新型高能鈍感含能材料的重要途徑[2-7]。六硝基六氮雜異伍茲烷(CL-20)是目前能量較高的單質(zhì)炸藥之一，但由于其機械感度高，使其應(yīng)用受到限制[8-9]。1,1-二氨基-2,2-二硝基乙烯(FOX-7)是一種新型高能鈍感炸藥[10]，其撞擊感度約為CL-20的5倍。根據(jù)共晶的形成原理[11]，若能將CL-20與FOX-7實現(xiàn)共晶，形成同時具有高能和低感特性的獨特結(jié)構(gòu)，將大大拓展CL-20的應(yīng)用范圍。

近年來，分子模擬計算在共晶含能材料的形成預(yù)測上起到了重要作用。衛(wèi)春雪等[12]模擬研究了HMX/TATB共晶體系的結(jié)構(gòu)與性能，結(jié)果表明對HMX 的5個主要生長面，TATB 在比表面能越小、生長速度越慢的晶面上取代形成的共晶能量越低, 結(jié)構(gòu)越穩(wěn)定；林鶴等[13]模擬研究了HMX/FOX-7共晶體系的結(jié)構(gòu)與性能，預(yù)測了共晶形成的可能性，并得出最易形成共晶的模型為FOX-7取代HMX(1 1 -1)晶面的模型；文國等[14]模擬研究了CL-20的晶習，并通過分子動力學模擬的方法推測了CL-20/TATB共晶體系的可行性；楊文升等[15]采用密度泛函方法對HMX/NQ共晶分子間的相互作用進行了模擬分析，結(jié)果表明HMX/NQ共晶的分子間作用本質(zhì)是一系列弱氫鍵和范德華力的共同作用。

本研究搭建了7種CL-20和FOX-7的共晶模型，并運用Materials studio軟件[16]中的Forcite模塊和Discover模塊對其進行幾何優(yōu)化和能量最小化計算，然后進行分子動力學模擬，通過Reflex模塊[17]進行X射線粉末衍射圖譜計算。通過考察共晶模型的結(jié)合能、徑向分布函數(shù)、X射線粉末衍射圖譜來判斷共晶是否生成，為CL-20/FOX-7共晶的形成提供理論依據(jù)。

1計算方法

1.1共晶模型的建立

根據(jù)共晶的形成原理，共晶體系中氫鍵是共晶形成的最重要的作用力[18]，CL-20含有極性的NO2基團，F(xiàn)OX-7含有NH2基團，二者可以形成N-O…H型氫鍵，考慮CL-20與FOX-7分子之間可能形成的分子間氫鍵，CL-20/FOX-7共晶中FOX-7作為氫鍵的給體，CL-20作為氫鍵的受體，構(gòu)建初始超分子結(jié)構(gòu)，建立摩爾比為1∶1 的CL-20/FOX-7共晶超分子結(jié)構(gòu)?；?CL-20的中子衍射單晶數(shù)據(jù)，建立ε-CL-20的晶胞和超晶胞(3×2×2)模型，共計48個CL-20分子。FOX-7分子對CL-20分子的取代分為隨機取代和在不同生長晶面上的取代，主要生長面模型采取切割分面法切割ε-CL-20超晶胞，用24個FOX-7分子隨機取代24個CL-20分子，以真空層厚度為200nm來構(gòu)建，共建立7種共晶模型。CL-20主要生長面的確定采用Morphology模塊中的Growth Morphology方法。

1.2共晶模型的優(yōu)化及分子動力學模擬

初始模型通過Forcite模塊下的Anneal過程進行200～500℃升、降溫退火循環(huán)處理，以消除超晶胞中不合理的能量和構(gòu)象，此過程步數(shù)為5000 步；之后通過Discover模塊下Minimizer過程將初始模型進行能量最小化；最后在Discover模塊下對前面優(yōu)化后的7 種共晶結(jié)構(gòu)模型進行NVT系綜下的MD模擬，采用Compass力場，NVT系綜，Andersen控溫方法，溫度設(shè)為298K，總時間為200ps，步長為1fs，范德華和靜電作用分別用Atom based和Ewald方法，此過程可充分弛豫超晶胞，為MD模擬提供合理、平衡的幾何構(gòu)象。FOX-7隨機取代CL-20(1 1 0)晶面的共晶模型如圖1所示。

圖1　FOX-7 隨機取代CL-20(1 1 0) 晶面的共晶模型Fig. 1　Co-crystal model of CL-20 (1 1 0) crystalsubstituted by FOX-7 randomly

1.3X射線粉末衍射

運用MS軟件包中的Reflex模塊, 基于CL-20/FOX-7共晶的平衡構(gòu)型進行X射線粉末衍射圖模擬，選擇Cu靶，掃描范圍5 °～75 °，步長為0.02 °。

2結(jié)果與討論

2.1共晶體系平衡的判別

選擇CL-20/FOX-7共晶(1 1 0)面為例，以同時滿足能量平衡和溫度平衡來判別分子動力學體系是否平衡。要達到能量平衡，要求體系平衡后能量沿恒定值變化波動較??；要達到溫度平衡，要求溫度上下波動差小于20K。體系平衡后溫度和能量隨時間波動的曲線如圖2所示。從圖2(a)可知，體系平衡后的溫度上下波動不超過10%；從圖2(b)可知,體系平衡后的能量波動較小，表明CL-20/FOX-7共晶模擬體系已達到平衡，其他幾種共晶模型具有相同的模擬結(jié)果。由圖2(b)還可知，CL-20/FOX-7共晶模型存在較強的非鍵能，由此可推測CL-20/FOX-7共晶存在較強的分子間作用力，有利于CL-20/FOX-7共晶的形成。

圖2　體系平衡后溫度和能量隨時間波動的曲線Fig.2　Curves of temperature and energy vs.simulation time

2.2CL-20 晶體主要生長面的預(yù)測

采用Mrophology模塊中Growth Morphology方法，使用Compass力場和力場自帶電荷，對CL-20晶體主要生長晶面進行預(yù)測,得出CL-20晶體有6個主要生長面，具體參數(shù)見表1。

從表1中可看出，晶面附著能主要由范德華力和靜電力構(gòu)成，且總附著能由大到小順序為Ea(0 1 1)>Ea(1 0 -1)>Ea(1 1 0)>Ea(0 0 2)>Ea(1 1 -1)>Ea(1 0 1)；晶面間距dhkl相對大小順序為：d(0 1 1)>d(1 0 -1)>d(1 1 0)>d(1 1 -1)>d(0 0 2)>d(1 0 1)；最大顯露面積的晶面為(0 1 1)面，其占總面積的比例為38.24%，相應(yīng)的附著能為-151.71 kJ/mol。

2.3CL-20/FOX-7 共晶的結(jié)合能

通過CL-20 的硝基和FOX-7的氨基在分子間形成氫鍵組裝，可以實現(xiàn)CL-20/FOX-7共晶。二者之間的結(jié)合能可定義為：

ΔE=EC-(ECL-20+EFOX-7)

(1)

Eb=-ΔE

(2)

式中:EC為CL-20/FOX-7共晶平衡后的單點能；ECL-20為CL-20/FOX-7共晶平衡后去除FOX-7之后的單點能；EFOX-7為CL-20/FOX-7共晶平衡后去除CL-20之后的單點能；結(jié)合能Eb定義為二者之間作用能ΔE的負值。

表1　真空中CL-20晶體的重要晶面及參數(shù)

注：M為多重度；dhkl為晶面間距；S為表面積；Etotal為總附著能；Eedw為范德華附著能；Eelc為靜電附著能；D為面心距離；Stotal為總顯露面積。

由Curie-Wulff原理可知，當晶體處于平衡狀態(tài)時，其總表面能最小，且各晶面的比表面能與生長速率成正比。

由于(1 0 1)晶面的生長速率最慢，因而具有最小的比表面能，所以(1 0 1)最有利于共晶的形成，具體結(jié)果見表2。

表2　CL-20各晶面與FOX-7的結(jié)合能參數(shù)

由表2可知，7 種共晶模型均有較強的結(jié)合能，其結(jié)合能大小順序為:Eb(1 0 1)>Eb(1 1 -1)>Eb(隨機晶面)>Eb(0 1 1)>Eb(0 0 2)>Eb(1 1 0)>Eb(1 0 -1),其中(1 0 1)晶面結(jié)合能最大。綜合以上兩點推測，CL-20/FOX-7在(1 0 1)晶面上總表面能最小，結(jié)合能最大，體系能量最低，更加穩(wěn)定，有利于CL-20/FOX-7共晶的形成。

2.4CL-20/FOX-7 共晶體系的徑向分布函數(shù)

徑向分布函數(shù)( RDF) 為系統(tǒng)區(qū)域密度與平均密度的比值。一般是指定某個特定粒子的坐標，在空間中其他粒子的分布幾率(離特定粒子的距離),用g(r)表示。因而徑向分布函數(shù)不僅可以研究物質(zhì)的有序性，同時也可以進行電子相關(guān)性的描述。通常情況下，分子間作用力包括氫鍵和范德華力。氫鍵鍵長在0.21～0.31nm之間，強范德華力相互作用鍵長范圍為0.31～0.50nm，弱范德華力相互作用鍵長大于0.50nm。

CL-20/FOX-7共晶結(jié)構(gòu)的徑向分布函數(shù)如圖3所示。

圖3　CL-20/FOX-7的7種共晶結(jié)構(gòu)的徑向分布函數(shù)Fig. 3　The radial distribution function g(r) of seven kinds of CL-20/FOX-7 co-crystal structures

由圖3可以看出，CL-20/FOX-7的7種共晶模型在r=0.18nm附近均有出峰，這表明CL-20的氧原子和 FOX-7的氫原子形成了分子間氫鍵，且(1 0 1)晶面和(1 1 -1)晶面模型的峰值較強，隨機取代模型的峰值最弱，(1 1 0)、(0 0 2)、(0 1 1)和(1 0 -1) 晶面模型的峰值處在兩者之間，說明(1 0 1)、(1 1 -1)晶面成氫鍵的作用力最強；在r=0.32nm處，7種共晶模型都有不同的出峰，這說明7種共晶模型同時存在著強范德華力作用，但峰值比形成氫鍵的峰值要小，說明共晶模型分子間作用力主要以氫鍵為主；由于(1 0 1)、(1 1 -1)晶面模型在氫鍵及范德華力的峰面積比其他模型大，由此可以預(yù)測(1 0 1)、(1 1 -1)晶面共晶模型存在著較強分子間作用力，較其他模型更容易形成共晶。

2.5XRD圖譜分析

經(jīng)分子動力學平衡后，7種共晶模型用Reflex模塊對其XRD圖譜進行了模擬，如圖4所示。

圖4　CL-20/FOX-7共晶模型的XRD粉末衍射圖譜Fig. 4　X-ray power diffraction patterns of CL-20/FOX-7co-crystal models

由圖4可以看出，純組分CL-20和FOX-7的特征峰強度在CL-20/FOX-7共晶模型XRD衍射圖譜中有所下降或消失，且7種共晶模型在2θ為5°附近均有新的衍射峰出現(xiàn)。這些特征衍射峰消失或新的衍射峰出現(xiàn)不是CL-20和FOX-7衍射峰的簡單疊加。根據(jù)這些特征峰消失、衍射峰強度下降和新的衍射峰出現(xiàn)，證明了CL-20/FOX-7共晶模型是一種有別于純組分CL-20和FOX-7的晶體，具有新的晶胞參數(shù)。主要是由于CL-20和FOX-7形成共晶時，分子間氫鍵及范德華力對FOX-7和CL-20的結(jié)構(gòu)進行破壞，形成了一種新的結(jié)構(gòu)。

3結(jié)論

(1)FOX-7和CL-20共晶模型存在較強的結(jié)合能，其大小順序為：Eb(1 0 1)>Eb(1 1 -1)>Eb(隨機晶面)>Eb(0 1 1)>Eb(0 0 2)>Eb(1 1 0) >Eb(1 0 -1)。

(2)通過徑向分布函數(shù)研究得到，F(xiàn)OX-7和CL-20分子間存在較強的氫鍵和范德華力作用，且氫鍵是形成共晶的主要動力。

(3)隨機取代模型的XRD衍射峰與純組分CL-20的衍射峰相近，6種主要生長面模型的XRD衍射峰位置和強度均有別于純組分CL-20或FOX-7，且均有新的衍射峰出現(xiàn)，表明晶格參數(shù)發(fā)生變化，形成了一種新的結(jié)構(gòu)。

參考文獻：

[1]Agrawal J P. Recent trends in high-energy materials[J]. Progress Energy and Combustion Science, 1998, 24: 1-7.

[2]舒遠杰, 武宗凱, 劉寧, 等. 晶形控制及形成共晶：含能材料改性研究的重要途徑[J]. 火炸藥學報, 2015, 38(5): 1-9.SHU Yuan-jie, WU Zong-kai, LIU Ning, et al. Crystal control and cocrystal formation：important route of modification research of energetic materials[J]. Chinese Journal of Explosives and Propellants, 2015, 38(5): 1-9.

[3]Landenberger K B, Matzger A J. Cocrystal engineering of a prototype energetic material: supramolecular chemistry of 2,4,6-trinitrotoluene[J]. Crystal Growth & Design, 2010, 10: 5341-5345.

[4]宋小蘭, 王毅, 宋朝陽, 等. CL-20/DNT 共晶炸藥的制備及其性能研究[J]. 火炸藥學報, 2016, 39(1): 23-27.

SONG Xiao-lan, WANG Yi, SONG Zhao-yang, et al. Preparation of CL-20/DNT cocrystal explosive and study on its performance[J]. Chinese Journal of Explosives and Propellants, 2016, 39(1): 23-27.

[5]Bolton O, Matzger A J. Improved stability and smart-material functionality realized in an energetic coccrystal[J]. Angewandte Chemie International Edition, 2011, 50:8960-8963.

[6]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]. Crystal Growth & Design,2012, 12, 5155-5158.

[7]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]. Crystal Growth & Design, 2012, 12, 4311-4314.

[8]徐金江, 蒲柳, 劉渝,等. HTPB 基粘結(jié)體系中ε-CL-20 的晶型轉(zhuǎn)變規(guī)律[J]. 含能材料, 2015, 23(2): 113-119.

XU Jin-jiang, PU Liu, LIU Yu, et al. Polymorphic transformation of ε-CL-20 in different HTPB-based composite systems[J]. Chinese Journal of Energetic Materials, 2015, 23(2): 113-119.

[9]Simpson R L, Urtiew P A, Ornellas D L, et al. CL-20 performanceexceeds that of HMX and its sensitivity is moderate[J]. Propellants, Explosives, Pyrotechnics, 1997, 22( 5) : 249-255．

[10] Badgujar D M, Talawar M B, Asthana S N, et al. Advances in science and technology of modern energetic materials: an overview[J]. Journal of Hazardous Materials, 2008, 151(2/3):289-305.

[11] Thomas S J M. Crystal engineering: origins early adventures and some current trend [J].Cryst Eng Comm, 2011, 13: 4304-4306.

[12] 衛(wèi)春雪, 段曉惠, 劉成建, 等. 環(huán)四甲撐四硝胺/1,3,5-三氨基-2,4,6-三硝基苯共晶炸藥的分子模擬研究[J]. 化學學報, 2009, 67: 2822-2826.WEI Chun-xue, DUAN Xiao-hui, LIU Cheng-jian, et al. Molecular simulation on co-crystal structure of HMX/TATB[J]. Acta Chimica Sinica, 2009, 67: 2822-2826.

[13] 林鶴, 張琳, 朱順官, 等. HMX/FOX-7共晶炸藥分子動力學模擬[J]. 兵工學報, 2012, 33(9): 1026-1030.

LIN He, ZHANG Lin, ZHU Shun-guan, et al. Molecular dynamic simulation of cyclotetramethylene tetranitramine /1,1-diamino-2, 2-dinitroethylene co-crystal explosive [J]. Acat Armamentarii, 2012, 33(9): 1026-1030.

[14] 文國. CL-20晶習研究及CL-20/TATB共晶的分子動力學模擬[D]. 太原: 中北大學, 2014.

WEN Guo. Morpho1ogy of CL-20 and molecular dynamic simulation on co-crystal of CL-20/TATB [D]. Taiyuan: North University of China, 2014.

[15] 楊文升, 茍瑞君, 張樹海, 等. HMX/NQ 共晶分子間相互作用的密度泛函理論研究[J]. 火炸藥學報, 2015, 38(6): 72-86.

YANG Wen-sheng, GOU Rui-jun, ZHANG Shu-hai, et al. Study on the intermolecular interaction of HMX/NQ cocrystal explosive by density functional theory[J]. Chinese Journal of Explosives and Propellants, 2015, 38(6): 72-86.

[16] Accelrys Inc. Materials Studio5.5[CP/CD]. San Diego, CA: Accelrys Inc., 2010.

[17] Engel G E, Wilke S, Kontig O, et al. Powder solve a complete package for crystal structure solution from powder diffraction patterns[J]. Applied Crystallography, 1999, 32: 1169-1179.

[18] Jayasankar A, Somwangthana A,Shao Z J, et al. Co-crystal formation during cogrinding and strorage is mediated by amorphous phase[J]. Pharm aceutical Research, 2006, 23: 2381-2392.

Molecular Dynamics Simulation of CL-20/FOX-7 Co-crystal

WU Zong-kai, SHU Yuan-jie, LIU Ning, DING Xiao-yong, WU Min-jie, WANG Ke, WANG Bo, LU Ying-ying

(Xi′an Modern Chemistry Research Institute, Xi′an 710065, China)

Abstract:Seven kinds of models of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) /1,1-diamino-2,2-dinitroethylene (FOX-7) co-crystals were constructed based on the hydrogen bond rules of co-crystal formation. The probability of CL-20/FOX-7 co-crystal formation was investigated using the molecular dynamics (MD) simulation. The radical distribution function (RDF) was used to study the acting force between the molecules in co-crystal model. X-ray powder diffraction (XRD) was used to simulate and analyze the difference between XRD peaks of co-crystal models and pure components. The results show that there are strong hydrogen bonds and Van Der Waals forces between CL-20 and FOX-7.The binding energies of seven kinds of co-crystal models decrease in the order ofEb(1 0 1)>Eb(1 1 -1)>Eb(randomness surface)>Eb(0 1 1)>Eb(0 0 2)>Eb(1 1 0)>Eb(1 0 -1).The intermolecular forces of seven kinds of co-crystal models are stronger when some CL-20 molecules are substituted by FOX-7 on the surface (1 0 1) and (1 1 -1) . The XRD peaks of seven kinds of co-crystal models are quite different from that of pure component CL-20 or FOX-7. Therefore, it can be predicted that CL-20 molecule is more easily substituted by FOX-7 on the surface (1 0 1) to obtain the stable co-crystal structure in the preparation process of the CL-20/FOX-7 co-crystal.

Keywords:co-crystal; molecular dynamics; 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane;CL-20; 1,1-diamino-2, 2-dinitroethylene;FOX-7

DOI：10.14077/j.issn.1007-7812.2016.03.007

收稿日期:2015-11-29；修回日期:2016-01-04

基金項目:國家自然科學基金資助項目(51373159)；國際(地區(qū))合作與交流項目(51511130036)

作者簡介:武宗凱(1989-)，男，助理工程師，從事含能材料合成、性能及理論計算研究。E-mail：zkwu204@sina.com

中圖分類號：TJ55;O641.3

文獻標志碼：A

文章編號：1007-7812(2016)03-0037-06