外电场影响HMX/MDNI复合物感度的理论研究

韩 刚， 苟瑞君， 张树海， 武春磊， 朱双飞

(1. 中北大学化工与环境学院， 山西 太原 030051; 2. 应用物理化学国家级重点实验室， 陕西 西安 710061)

1 引 言

现代武器装备的发展对含能材料的性能提出了越来越高的要求，高能量密度化合物感度与能量之间的矛盾严重。如何改善炸药的性能尤其是感度特性一直是国内外学者研究的重点。

在含能材料领域，已有研究证明外电场可改变分子的结构从而改变某些性能(如感度)。Politzer等[1-3]研究发现，不同大小和方向的外加电场对几种典型含能分子的结构和性能都有较大的影响。Ren[4]的研究表明，外加电场对CH3NO2的N—O键和NH2NO2的N—N键影响显著，并且有可能改变硝基甲烷的分解反应过程、炸药感度及爆速和爆压[5]。基于此，可知外电场对单体含能分子的性能影响较大。实际应用中，两种或多种含能组分形成复合体系的情况很多，然而，外电场对两种(或多种)含能分子所构成的复合物的影响还未见报道。若对复合体系施加适当的外电场，则也有可能改善炸药性能，使其感度降低。因此研究外电场对炸药复合体系性能的影响具有实际意义。

奥克托今(HMX)是一种高能量密度材料，但感度较高因而其应用受到了一定的限制，目前已有报道表明含HMX共晶炸药的感度相对于纯HMX有所降低[6-8]。1-甲基-4,5-二硝基咪唑(MDNI)是一种新型钝感咪唑类炸药，易于合成，有可能与高能敏感炸药形成共晶从而降低高能组分的感度。本课题组前期对HMX/MDNI共晶炸药进行了初步的理论研究[9]，发现共晶结构中—CH3(MDNI)和—NO2(HMX)形成的分子间氢键可使HMX的感度降低。由于实验操作的复杂性及危险性，本研究借助计算机模拟技术对不同外电场下HMX/MDNI复合体系进行了优化，然后从结构、分子间相互作用、引发键强度、硝基基团电荷、电子密度及分子表面静电势等方面研究了不同外电场对复合体系的影响，以期为含能材料的感度调节和爆炸特性实验研究提供一定理论基础。

2 计算方法

本研究首先借助Gaussian 09软件包[10]，采用DFT-B3LYP方法，在6-311++G(d,p)水平下优化得到了3种稳定的HMX/MDNI复合物(Ⅰ, Ⅱ, Ⅲ)。由于与MDNI形成分子间氢键的HMX的引发键(N—NO2)是我们关注的重点，所以外电场正方向定义为由N指向NO2，负方向则由NO2指向N。对所得的3种稳定构型分别施加不同大小的外电场 (±0.005和±0.010 a.u.)，在相同水平下获得了外电场作用时的稳定几何构型。然后采用M06-2X和B3LYP两种方法计算了几种稳定构型的相互作用能 (Eint)和引发键解离能(EBDE); 对施加外电场和未施加外电场的稳定构型进行电子密度拓扑和电子密度转移分析，再借助Multiwfn软件[11]分析了不同外加电场下硝基基团电荷 (QNO2)及分子结构的表面静电势的变化情况，以研究外电场对炸药感度的影响。

3 结果与讨论

3.1 结构分析

在B3LYP/6-311++G(d,p)水平下分别对3种初始构型施加不同大小的外电场进行优化(-0.010～0.010 a.u.)，以获得最稳定的构型(共15个)，未施加外电场时优化所得的稳定构型如图1所示。构型Ⅰ为—NO2(HMX)与—CH/—CH3(MDNI)作用; 构型Ⅱ为—NO2(HMX)与—CH3(MDNI)，—CH2(HMX)与—NO2(MDNI)作用; 构型Ⅲ为—NO2(HMX)与—CH(MDNI)，—CH2(HMX)与—NO2(MDNI)作用。由图1可知，3种复合构型中两分子间的H…O距离均在氢键作用范围内，分别为2.363～2.549 Å，2.344～2.501 Å和2.334～2.633 Å。在外电场作用下，H…O间距离发生了不同程度的变化，分别变为2.170～2.670 Å，2.199～2.606 Å和2.217～2.662 Å，表明外电场对HMX/MDNI复合物中的氢键作用影响程度不同。

图1 未施加外电场时3种复合物优化后的稳定构型

Fig.1 Optimized configurations of 3 complexes without external electric field

将所得的15个复合物构型中的HMX分子引发键键长数据列于表1中。从表1可以看出，对于每种构型，随着正向外电场的增加，Ⅰ、Ⅱ和Ⅲ的引发键键长都有所减小，分别由1.374, 1.392, 1.392 Å减少至1.351, 1.377, 1.378 Å，表明此时引发键键长有变短的趋势; 负向时则由未施加外电场时的键长分别增加至1.445, 1.430, 1.451 Å。引发键长度越短则强度越强，键长分析表明正向外电场可使引发键强度增加，增强炸药稳定性，降低感度; 负向外电场时则可能削弱其稳定性，这与Ren[4]对单质炸药所得的研究结果一致。

表1 不同大小外电场下复合物中HMX引发键键长

Table 1 Trigger bond lengths of HMX in complexes with different external electric fields

externalelectricfield/a．u．L/ÅⅠⅡⅢ-0．0101．4451．4301．451-0．0051．3901．4171．41601．3741．3921．3920．0051．3661．3821．3920．0101．3511．3771．378

3.2 相互作用能和解离能

已有研究表明，相对于B3LYP和M06-2X，MP2方法计算所得键解离能误差较大[12]，因此采用B3LYP和M06-2X两种方法计算了3种构型不同外电场作用下的相互作用能(经BSSE校正)和各自的键解离能 (Eint)。其中，分子间相互作用能可表示为:

ΔE=-Eint=-(EHMX/MDNI-EHMX-EMDNI+EBSSE)

式中，Eint为相互作用能，kJ·mol-1;EHMX/MDNI为复合物的总能量，kJ·mol-1;EHMX(EMDNI) 为去除MDNI(HMX)后的总能量，kJ·mol-1;EBSSE为矫正能，kJ·mol-1。

由Tan[13]的研究可知，键解离能可表示为:

EBDE=E·R+E·NO2-ER-NO2

式中，EBDE为键解离能，kJ·mol-1;E·R，E·NO2，ER-NO2分别表示·R，·NO2，R—NO2的能量，kJ·mol-1。

由上述公式所得Eint和EBDE的数据列于表2中。

表2 不同大小外电场下的分子间相互作用能(经BSSE校正)和HMX的键解离能

Table 2 Intermolecular interaction energy (with BSSE corrected) and bond dissociation energy of HMX with different external electric fields

externalelectricfield/a．u．methodΔE/kJ·mol-1ⅠⅡⅢEBDE/kJ·mol-1ⅠⅡⅢ-0．010B3LYP29．7346．5846．97161．08174．50185．52M06⁃2X45．2373．0061．99216．31235．89245．83-0．005B3LYP37．5933．7638．24179．41191．00201．08M06⁃2X54．0259．0158．90238．07257．36261．880B3LYP16．9129．9139．78200．83214．26206．02M06⁃2X21．4747．4964．26259．21282．01268．740．005B3LYP30．4626．3531．87219．74228．74217．40M06⁃2X38．8431．1440．44281．04294．47282．430．010B3LYP47．2248．9648．53246．81237．44228．03M06⁃2X57．1958．8857．98310．54303．30293．55

图2 HMX的键解离能与外电场关系图

Fig.2 The relationship between bond dissociation energy (EBDE) of HMX and external electric field

3.3 硝基基团电荷

已有研究表明，可通过硝基基团所带负电荷多少来揭示炸药感度的大小[15]。根据以下公式，计算了B3LYP/6-311++G(d,p)水平下不同外电场作用时形成分子间氢键的硝基基团的电荷并列于表3中。

QNO2=QN+QO1+QO2

式中，QNO2代表硝基基团的电荷，e;QN、QO1、QO2分别代表硝基基团中N和O原子上的电荷，e。

硝基基团所带负电荷越多，炸药越钝感，负电荷越少，炸药越敏感。从表3可以看出，相对于未施加外电场的情况，负方向外电场可使硝基基团所带负电荷减少: 对于Ⅰ，从-0.1276降低至0.0068; 对于Ⅱ，从-0.1351减少至-0.0458; 对于Ⅲ，从-0.1502减少至-0.0059，这表明负方向外电场可使炸药的感度增加。外电场为正值时，Ⅰ的硝基基团负电荷增加了，表明此时更多的负电荷聚集在硝基周围，炸药感度降低; 但Ⅱ和Ⅲ的硝基所带负电荷变化复杂: 随正向外电场增加，有升有降，无法通过硝基基团电荷多少判断炸药感度的变化，这也说明了外电场对炸药感度影响的复杂性，需要综合其它因素分析考虑。

表3 复合物中的硝基基团电荷

Table 3 Charges of nitro groups in complexes

externalelectricfield/a．u．ⅠⅡⅢ-0．0100．0068-0．0458-0．0059-0．0050．0296-0．0800-0．07600-0．1276-0．1351-0．15020．005-0．1528-0．1594-0．03730．010-0．2079-0．0503-0．0954

3.4 电子密度拓扑

3.5 电子密度转移

电子密度分布的改变可以被用于研究键强度的变化[18]，近年已有研究者借助电子密度转移来分析复合物与单体的引发键强度的变化，从而分析炸药的感度变化[19]，但还未见外电场作用下炸药分子电子密度转移的报道。本研究通过分析不同外加电场下各复合物的电子密度转移情况，得出了外电场对感度的影响。为了直观的观察炸药感度的变化，以构型Ⅰ为例，将Ⅰ的电子密度在不同外电场下的变化情况画于图3中进行分析，紫色代表电子聚集，黄色代表失去电子。

从图3可以看出，不施加外电场时，氢原子(C—H)附近呈黄色而硝基氧(N—NO2)被紫色围绕，这是因为HMX与MDNI形成了分子间氢键，发生了电子转移(氢原子失去电子)所致。HMX的引发键(N—NO2)周围出现了紫色区域，表明电子密度也转移到了引发键上。值得一提的是，在不同外电场下，构型Ⅰ中引发键周围紫色区域不同，表明不同外电场对电子密度转移的影响不同。正向外电场下，HMX的引发键周围紫色区域比负向外电场时明显，表明此处聚集了较多的电子(构型Ⅱ和构型Ⅲ也有相同现象)。两原子间电子越集中，其重叠的几率越大，引发键的强度就越高。因此，正向外电场下复合物的引发键强度比负向外电场的高，炸药感度比负向外电场时的低。这与解离能分析所得出的结论一致。

3.6 分子表面静电势

分子表面静电势(ESP)源自分子的总净电荷，是分子的全域性特征，因此分析炸药分子的表面静电势变化对研究炸药的感度变化更具有指导意义。Politzer[20]指出，不仅是引发键，分子表面静电势也应该被用于分析炸药的感度。已有研究利用ESP分析共晶炸药的感度变化[9]，但外电场作用下炸药的ESP相关研究还未见报道。图4以不同外电场下构型Ⅰ的分子表面静电势变化为例进行分析，相对于无外电场的静电势，外电场存在时，正/负表面静电势最大/小值数值都变大，如表4所示，且外电场强度越大变化越大，与图4中静电势图颜色变深一致，表明强外电场对电荷迁移影响较大。图4中分子连接处的蓝色变浅，正向外电场时甚至变为白色，表明正向外电场对电荷转移的影响较大，构型Ⅱ和Ⅲ也有同样的现象。这主要是因为硝基氧和C—H的正负静电势重合，改变了两分子间的静电势，说明外电场作用下炸药的感度有可能发生变化。

a. -0.010 a.u.b. -0.005 a.u.c. 0 a.u.

d. 0.005 a.u.e. 0.010 a.u.

图3 不同外电场下构型Ⅰ的电子密度转移图

Fig.3 Shifts of electron density for configuration Ⅰ with different external electric fields

a. -0.010 a.u.b. -0.005 a.u.c. 0 a.u.

d. 0.005 a.u.e. 0.010 a.u.

图4 不同外电场下构型Ⅰ的静电势图

Fig.4 Electrostatic surface potentials for configuration Ⅰ with different external electric fields

表4 复合物的最大/小分子表面静电势

Table 4 The maximum and minimum electrostatic surface potentials of complexes

externalelectricfield/a．u．Vmax/VminⅠⅡⅢ-0．010297．591/-219．124324．101/-212．426315．394/-205．146-0．005235．371/-170．837275．131/-167．385255．174/-157．4980261．454/-141．750215．074/-118．646206．125/-118．8130．005302．838/-186．364283．161/-175．854267．345/-189．3010．010349．184/-224．844297．336/-228．007293．901/-227．413

Note:VmaxandVminare the maximum and minimum electrostatic surface potentials, respectively.

表5 N—NO2键的局域正静电势极值

Table 5 The maximum local positive electrostatic surface potential of selected N—NO2bond

externalelectricfield/a．u．Vs，max/kJ·mol-1ⅠⅡⅢ-0．010154．674(-4．059)182．489(16．619)177．900(20．184)-0．005119．591(-39．142)172．180(6．31)196．870(39．154)0158．733165．870157．7160．005140．859(-17．874)162．054(-3．816)147．478(-10．238)0．010125．457(-33．276)143．227(-22．643)73．166(-84．55)

Note:Vs,maxis the maximum local positive electrostatic surface potential, values in parentheses are increments compared withVs,maxof no external electric field.

此外，Politzer等[21]指出，C/N—NO2键周围的局域正静电势极值点Vs,max与感度存在正相关关系，即Vs,max值越大，感度越高，将涉及到氢键的N—NO2键的Vs,max列于表5中。由表中数据可以看出，相对于未施加外电场时的极值，在不同大小外电场作用下，构型Ⅰ的Vs,max都减小了，表明感度降低; 对于构型Ⅱ和Ⅲ，虽然变化幅度不同，但负向外电场都使Vs,max增大，且正向外电场都使相应的Vs,max减小，说明负向外电场使复合物Ⅱ和Ⅲ中的HMX感度增加了，而正向外电场则使其感度降低，这与前面的分析一致。然而，必须指出: 炸药的感度变化复杂，受到多种因素的影响，必须全面考虑。

4 结 论

在M06-2X∥B3LYP/6-311++G(d,p)水平下对3种HMX/MDNI复合物施加了不同的外加电场(-0.010～0.010 a.u.)，分析了不同外电场下复合物中的引发键键长、解离能、电子密度拓扑、电子密度转移和分子表面静电势的变化，得到如下结论:

(1) 随着正向外电场的增加，HMX/MDNI复合物中HMX的引发键变短，负向外电场增加时，引发键则变长; 引发键解离能与外加电场存在正相关关系，正向外电场增加则引发键解离能增加，炸药的感度相应降低。

(2) 电子密度拓扑分析表明，当施加负向外电场时，3种HMX/MDNI复合物中均出现了色散作用O…O/N，但ρ值仅为约0.0040 a.u.，说明此时色散作用极其微弱; 而O…H对应的ρ值介于0.002～0.04 a.u.，表明仍然存在微弱的氢键作用。

(3) 电子密度转移分析一方面说明了两分子间电子转移受到了外加电场的影响; 另一方面，复合物在正向外电场作用下，其引发键周围紫色区域(电子聚集)更加明显，电子重叠几率增加，引发键强度增强。这与引发键解离能分析一致。

(4) 不同外电场下HMX/MDNI复合体系的Vmax和Vmin变化表明强外电场对电子迁移的影响较大; HMX引发键(N—NO2)的局域正静电势极值点Vs,max的变化说明，负向外电场趋于增加炸药感度，而正向外电场则可降低炸药的感度。

参考文献:

[1]Politzer P, Murray J S, Concha M C, et al. Effects of electric fields upon energetic molecules: nitromethane and dimethylnitramine [J].CentralEuropeanJournalofEnergeticMaterials, 2007, 4(4): 3-21.

[2]Politzer P, Murray J S. Computed effects of electric fields upon the C—NO2and N—NO2bonds of nitromethane and dimethylnitramine [J].InternationalJournalofQuantumChemistry, 2009, 109(1): 3-7.

[3]Politzer P, Murray J S, Lane P. Computational determination of effects of electric fields upon “trigger linkages” of prototypical energetic molecules [J].InternationalJournalofQuantumChemistry, 2009, 109(3): 534-539.

[4]REN Fu-de, CAO Duan-lin, SHI Wen-jing, et al. A theoretical prediction of the possible trigger linkage of CH3NO2and NH2NO2in an external electric field [J].JournalofMolecularModeling, 2015, 21: 145.

[5]REN Fu-de, CAO Duan-lin, SHI Wen-jing. A dynamics prediction of nitromethane → methyl nitrite isomerization in external electric field [J].JournalofMolecularModeling, 2016, 22: 96.

[6]LIN He, ZHU Shun-guan, ZHANG Lin, et al. Intermolecular interactions, thermodynamic properties, crystal structure, and detonation performance of HMX/NTO cocrystal explosive [J].InternationalJournalofQuantumChemistry, 2013, 113(10): 1591-1599.

[7]LIN He, ZHU Shun-guan, LI Hong-zhen, et al. Structure and detonation performance of a novel HMX/LLM-105 cocrystal explosive [J].JournalofPhysicalOrganicChemistry, 2013, 26(11): 898-907.

[8]LIN He, ZHU Shun-guan, ZHANG Lin, et al. Synthesis and first principles investigation of HMX/NMP cocrystal explosive [J].JournalofEnergeticMaterials, 2014, 31(4): 261-272.

[9]HAN Gang, GOU Rui-jun, REN Fu-de, et al. Theoretical investigation into the influence of molar ratio on binding energy, mechanical property and detonation performance of 1,3,5,7-tetranitro-1,3,5,7-tetrazacyclo octane (HMX)/1-methyl-4,5-dinitroimidazole (MDNI) cocrystal explosive [J].ComputationalandTheoreticalChemistry, 2017, 1109: 27-35.

[10]Frisch M J, Trucks G W, Schlegel H B, et al. Gaussian 09 [CP]∥Revision A.01, Gaussian Inc, Wallingford, CT, 2009.

[11]LU Tian, CHEN Fei-wu. Multiwfn: A multifunctional wavefunction analyzer [J].JournalofComputationalChemistry, 2012, 33(5): 580-592.

[12]Budyka M F, Zyubina T S, Zarkadis A K. Correlating ground and excited state properties: a quantum chemical study of the photodissociation of the C—N bond in N-substituted anilines [J].JournalofMolecularStructure:Theochem, 2002, 594(1-2): 113-125.

[13]TAN Bi-sheng, LONG Xin-ping, PENG Ru-fang, et al. Two important factors influencing shock sensitivity of nitro compounds: bond dissociation energy of X—NO2(X=C, N, O) and mulliken charges of nitro group [J].JournalofHazardousMaterials, 2010, 183(1-3): 908-912.

[14]SONG Xiao-shu, CHENG Xin-lu, YANG Xiang-dong, et al. Relationship between the bond dissociation energies and impact sensitivities of some nitro-explosives [J].Propellants,Explosives,Pyrotechnics, 2006, 31(4): 306-310.

[15]ZHANG Chao-yang, SHU Yuan-jie, HUANG Yi-gang, et al. Investigation of correlation between impact sensitivities and nitro group charges in nitro compounds[J].TheJournalofPhysicalChemistryB, 2005, 109(18): 8978-8982.

[16]Bader R F W. Atoms in Molecules: A quantum Theory [M]. New York: Oxford University Press, 1990.

[17]Nishiyama Y, Langan P, Wada M, et al. Looking at hydrogen bonds in cellulose [J].ActaCrystallographicaSectionD, 2010, 66: 1172-1177.

[18]Scheiner S, Kar T. Red-versus blue-shifting hydrogen bonds: are there fundamental distinctions? [J].TheJournalofPhysicalChemistryA, 2002, 106(9): 1784-1789.

[19]QIU Wei, REN Fu-de, SHI Wen-jing. A theoretical study on the strength of the C—NO2bond and ring strain upon the formation of the intermolecular H-bonding interaction between HF and nitro group in nitrocyclopropane, nitrocyclobutane, nitrocyclopentane or nitrocyclohexane [J].JournalofMolecularModeling, 2015, 21: 114.

[20]Politzer P, Murray J S. Impact sensitivity and crystal lattice compressibility/free space [J].JournalofMolecularModeling, 2014, 20: 222-2230.

[21]Politzer P, Murray J S. Relationships between dissociation energies and electrostatic potentials of C—NO2bonds: applications to impact sensitivities [J].JournalofMolecularStructure, 1996, 376: 419-424.