2. 北京农业职业学院,北京 102442
3. 河南理工大学 物理化学学院,河南 焦作 454000
2. Beijing Vocational College of Agriculture, Beijing 102442, China
3. School of Physics and Chemistry, Henan Polytechnic University, Jiaozuo 454000, China
高能量密度材料(HEDM)作为武器系统、民用爆破、喷涂加工等领域的能源材料,其性能的改进将对一个国家的军事、生产和生活产生重大影响,所以兼具高能、低耗、安全和环保特性的高能量密度化合物的设计与合成已成为各国争先研究的焦点.全氮化合物是一种仅靠自身分解即可产生高能、无毒无污染的超高能量密度物质,其能量来源于分子键能.N—N的键能是159.8 kJ/mol,NN的键能是418.0 kJ/mol,N≡N的键能是954.0 kJ/mol,全氮化合物完全分解成N2时,由于键接方式改变释放出的能量是现有炸药的几倍,甚至几十倍,且最终产物是环境友好的N2.
1891年,Curtius和Radenhausen发现了单质N2之外的全氮物质N3-离子[1].之后,国内外学者对各种全氮分子设计、研究、预测的步伐从未停歇[2-10].Christe等成功合成了蕴含巨大爆炸能量的链状结构N5+,并通过理论计算预测其稳定性[2-3].Cacace等通过电子轰击氮气得到了亚稳态的N4+[4].根据已合成的N3-和N5+,科研工作者预测今后最可能合成的物质是N8,并进行了相关研究[5-10],例如Hirshberg等报道了理论预测的不同压力下N8晶体的存在形式[7].本文拟从分子设计的角度出发,采用高精度的量子化学计算方法预测不同N8分子的结构、热稳定性和爆轰性能,筛选出6种优秀的高能材料候选分子.
1 计算方法本文采用Gaussian 09程序包[11]中的三种量子化学计算方法B3LYP,MP2和CCSD结合6-31G*和6-311++G**基组对20种设计新颖的N8分子进行优化和振动频率分析,对筛选出的6种稳定构型进行了性质预测,包括Wiberg键级、分子轨道、分子静电势、能量和密度、生成热以及爆速和爆压等.
生成热、爆速和爆压是衡量高能材料爆轰性能的重要参数.生成热可通过原子化能方法[12]分两步求得,首先计算分子在0 K时的生成热:
其次计算分子在298.15 K时的生成热:
其中:M表示分子;X表示组成分子M的每一种元素;x表示分子M中元素X的原子个数;ε0(M)表示分子的总能量;xε0(X)表示元素X的总能量;εZPE(M)表示分子的零点能;Δf
其中:D为炸药的爆速(m/s);P为炸药的爆压(GPa);ρ0为装药密度(g/cm3);∑NW为炸药的修正氮当量;Pi为每克炸药爆轰时所形成第i种产物的摩尔数;NPi为第i种产物的修正氮当量系数;Bk为第k种化学键在分子中出现的次数;NBk为第k种化学键的修正氮当量系数;Gj为第j种基团在分子中出现的次数;NGj为第j种基团的修正氮当量系数.
2 结果与讨论 2.1 几何结构利用三种方法B3LYP,MP2和CCSD优化均可得到图 1所示的6种稳定的N8同分异构体,并通过频率分析确认它们是N8势能面上的极小值点(以B3LYP数据为例进行结果讨论,其他方法可得到相同结论).其中,结构1为平面双五元环,D2h点群,Wiberg键级分析表明N2—N5最弱(1.063 8),为单键,N1—N3,N3—N4,N6—N7和N6—N8四键相等且最强(1.479 8),介于单键和双键之间,表明分子中存在共轭π键.结构2是C2h点群的梯状结构,N4—N7和N5—N6键级相等(1.873 3)且接近于双键,其余皆为单键.N2—N1—N5和N3—N8—N7的角度分别是84.59°和84.65°,远远小于N原子sp3等性杂化角度107.3°,表明分子内部存在较大张力.结构3和4是C2v点群的笼状结构,各键的键级值均接近1,为单键氮化物.结构5近似于链状,C2v点群,N1—N4,N2—N4,N6—N7和N7—N8四键相等且最弱(0.840 1),N1—N2和N6—N8为双键.结构6为C3点群的单键立方结构,分子张力较大.通过上述对分子化学键强度的分析可以预测爆炸反应发生的最可能途径.
表 1和图 2给出了6种N8分子在B3LYP/6-311++G**计算水平下得到的能量、前线轨道及静电势.前线轨道的能级差ε为分子的最低空轨道LUMO和最高占据轨道HOMO的能量之差,ε越大,电子越难转移和跃迁,分子越稳定.因此,经过对比发现结构1的能量最低且ε最大,表明其发生电子转移和跃迁比较困难,相对热稳定性最高.
分子中与核距离不等的空间各点处的静电势是不同的.分子静电势图(图 2)可以清晰直观地描述出分子中电子的空间排布情况,蓝色和红色分别代表正电性和负电性.通过对分子中电子富集和缺失区域的观察,能够预测分子的活性位点.位点附近的电子密度越大,负电性越高,越容易发生化学反应.
表 2给出了6种N8分子在B3LYP/6-31G*计算水平下得到的生成热,数值大小顺序为:结构1<5<2<4<3<6.结合图 1中各分子结构,发现含有共轭键的分子的生成热较小,不含共轭键且张力较大的分子的生成热较大.例如,结构1的双五元环中存在较强的离域π键,分子结构稳定、张力较小,生成热最小;立方烷结构6完全由单键构成且存在很大张力,从而导致分子不稳定,生成热最大.
表 3列出了采用修正氮当量公式在B3LYP/6-311++G**计算水平下[10]得到的6种N8分子的密度和爆速爆压,并与黑索金(RDX)和奥克托今(HMX)的计算值[14-15]和实验值[16]进行了比较.结果表明,6种N8分子的性质(ρ,D,P)排序为:笼状结构(3,4,6)>环状结构(1,2)>链状结构(5).N8分子的爆速和爆压大于RDX和HMX的爆速和爆压,因此,这6种N8分子是很有潜力的优秀高能量密度材料分子.
本文设计并筛选出6种新颖的N8同分异构体,采用高精度的量子化学计算方法预测其热稳定性和爆轰性能,发现6种N8分子的爆速和爆压大于爆轰性能优越的RDX和HMX,是非常理想的环保型高能量密度材料的候选目标.本文为此类新型全氮高能材料的研发提供了可靠的理论依据,针对该类分子的热裂解途径和机理等后续研究正在进行中.
[1] | Curtius T, Radenhausen R. Zur kenntniss der stickstoffwasserstoffsäure[J]. Journal für Praktische Chemie, 1891, 43 (1) : 207 –208. (0) |
[2] | Christe K O, Wilson W W, Sheehy J A, et al. N5+:a novel homoleptic polynitrogen ion as a high energy density material[J]. Angewandte Chemie International Edition, 1999, 38 (13/14) : 2004 –2009. (0) |
[3] | Vij A, Pavlovich J G, Wilson W W, et al. Experimental detection of the pentaazacyclopentadienide (pentazolate) anion,cyclo-N5-[J]. Angewandte Chemie International Edition, 2002, 41 (16) : 3051 –3054. (0) |
[4] | Cacace F, De Petris G, Troiani A. Experimental detection of tetranitrogen[J]. Science, 2002, 295 (5554) : 480 –481. (0) |
[5] | Lauderdale W J, Stanton J F, Bartlett R J. Stability and energetics of metastable molecules:tetraazatetrahedrane (N4),hexaazabenzene (N6),and octaazacubane (N8)[J]. The Journal of Physical Chemistry, 1992, 96 (3) : 1173 –1178. (0) |
[6] | Chung G, Schmidt M W, Gordon M S. An ab initio study of potential energy surfaces for N8 isomers[J]. The Journal of Physical Chemistry A, 2000, 104 (23) : 5647 –5650. (0) |
[7] | Hirshberg B, Gerber R B, Krylov A I. Calculations predict a stable molecular crystal of N8[J]. Nature Chemistry, 2014, 6 (1) : 52 –56. (0) |
[8] | Wang L J, Li S, Li Q S. Theoretical studies on a possible synthesis reaction pathway on N8 (CS) clusters[J]. Journal of Computational Chemistry, 2001, 22 (13) : 1334 –1339. (0) |
[9] | Li Y C, Qi C, Li S H, et al. 1,1’-Azobis-1,2,3-triazole:a high-nitrogen compound with stable N8 structure and photochromism[J]. Journal of the American Chemical Society, 2010, 132 (35) : 12172 –12173. (0) |
[10] |
牛晓庆, 张建国, 王颖, 等. 叠氮唑类高氮含能化合物的理论研究[J].
化学学报, 2011, 69 (6) : 610 –616.
( Niu Xiao-qing, Zhang Jian-guo, Wang Ying, et al. Theoretical studies of azide-azole nitrogen-rich energetic compounds[J]. Acta Chimica Sinica, 2011, 69 (6) : 610 –616. ) (0) |
[11] | Frisch M J, Trucks G W, Schlegel H B, et al. Gaussian 09,revision A.02[M]. Wallingford: Gaussian Inc, 2009 . (0) |
[12] | Curtiss L A, Raghavachari K, Redfern P C, et al. Assessment of Gaussian-2 and density functional theories for the computation of enthalpies of formation[J]. Journal of Chemical Physics, 1997, 106 (3) : 1063 –1079. (0) |
[13] | Curtiss L A, Redfern P C, Smith B J, et al. Gaussian-2 (G2) theory:reduced basis set requirements[J]. Journal of Chemical Physics, 1996, 104 (13) : 5148 –5152. (0) |
[14] |
国迂贤, 张厚生. 炸药爆轰性质计算的氮当量公式及修正氮当量公式:炸药爆速的计算[J].
爆炸与冲击, 1983, 3 (3) : 56 –66.
( Guo Yu-xian, Zhang Hou-sheng. Nitrogen equivalent (NE) and modified nitrogen equivalent (MNE) equations for predicting detonation parameters of explosives:prediction of detonation velocity of explosives[J]. Explosion and Shock Waves, 1983, 3 (3) : 56 –66. ) (0) |
[15] |
张厚生. 用氮当量公式及修正氮当量公式计算炸药爆轰压力[J].
爆炸与冲击, 1984, 4 (2) : 79 –82.
( Zhang Hou-sheng. Prediction of detonation pressures of explosives by nitrogen equivalent (NE) and modified nitrogen equivalent (MNE) equations[J]. Explosion and Shock Waves, 1984, 4 (2) : 79 –82. ) (0) |
[16] | Wang R H, Xu H Y, Guo Y, et al. Bis [3-(5-nitroimino-1,2,4-triazolate)]-based energetic salts:synthesis and promising properties of a new family of high-density insensitive materials[J]. Journal of the American Chemical Society, 2010, 132 (34) : 11904 –11905. (0) |