《高能量密度材料3,3'-偶氮-1,2,4,5-四嗪衍生物的分子设计.pdf》由会员分享,可在线阅读,更多相关《高能量密度材料3,3'-偶氮-1,2,4,5-四嗪衍生物的分子设计.pdf(10页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、高能量密度材料3,3-偶氮-1,2,4,5-四嗪衍生物的分子设计张婧婧1高洪伟2卫涛2,*王朝杰2(1温州医学院仁济学院,浙江 温州325035;2温州医学院药学院,浙江 温州325035)摘要:运用密度泛函理论(DFT)方法,计算系列3,3-偶氮-1,2,4,5-四嗪衍生物的生成热.结果显示:N3取代基在增加3,3-偶氮-1,2,4,5-四嗪衍生物的生成热方面起了非常重要的作用.通过分析标题化合物的最弱键离解能发现:NH2或N3取代基非常有利于增加衍生物的热稳定性.计算的爆速(D)和爆压(p)数值表明:NO2或NF2取代基有利于提高3,3-偶氮-1,2,4,5-四嗪衍生物的爆轰性能.综合爆轰
2、性能和热稳定性的计算结果,3种3,3-偶氮-1,2,4,5-四嗪衍生物可以作为潜在的品优高能量密度材料(HEDM)候选物.关键词:密度泛函理论;3,3-偶氮-1,2,4,5-四嗪;生成热;键离解能;爆轰性能中图分类号:O641Molecular Design of 3,3-Azobis-1,2,4,5-tetrazine-Based High-EnergyDensity MaterialsZHANG Jing-Jing1GAO Hong-Wei2WEI Tao2,*WANG Chao-Jie2(1Renji College,Wenzhou Medical College,Wenzhou3250
3、35,Zhejiang Province,P.R.China;2Department of Pharmacy,Wenzhou Medical College,Wenzhou325035,Zhejiang Province,P.R.China)Abstract:We systematically studied the heats of formation(HOFs)for a series of 3,3-azobis-1,2,4,5-tetrazine derivatives by density functional theory(DFT).The results show that the
4、 N3group plays avery important role in increasing the HOFs for these derivatives.An analysis of the bond dissociationenergies for the weakest bonds indicates that the attachment of NH2or N3group to 3,3-azobis-1,2,4,5-tetrazine is favorable in enhancing its thermal stability.The calculated detonation
5、 velocities(D)andpressures(p)indicates that NO2or NF2largely enhances the detonation performance of thederivatives.Considering the detonation performance and the thermal stability,the three derivatives may beregarded to be promising candidates for high-energy density materials(HEDMs).Key Words:Densi
6、ty functional theory;3,3-Azobis-1,2,4,5-tetrazine;Heat of formation;Bond dissociation energy;Detonation propertyA物理化学学报(Wuli Huaxue Xuebao)Acta Phys.Chim.Sin.,2010,26(12):3337-3344DecemberReceived:July 19,2010;Revised:September 19,2010;Published on Web:November 1,2010.Corresponding author.Email:;Tel
7、:+86-577-86689984.The project was supported by the Natural Science Foundation of Zhejiang Province,China(Y5080043).浙江省自然科学基金(Y5080043)资助项目 Editorial office ofActa PhysicoChimica Sinica六元杂氮环成为传统含能材料的基本单元,特别是S-四嗪环(1,2,4,5-四嗪环).近些年,S-四嗪衍生物作为高能量密度材料(HEDM)无论在理论还是在实验合成方面都得到了广泛关注1-7,归因于它们具有高密度、高正生成热以及较好的热稳
8、定性等特点.S-四嗪衍生物中3,3-偶氮-(6-氨基-1,2,4,5-四嗪)(DAAT)以及类似衍生物的合成和性能研究已有相关报道8-11,具有良好的应用前景,一些可作为潜在的高爆炸物、火箭推进剂以及烟火型气体发生剂原料等.尽管不含氧,但DAAT展现出了较好的爆轰性能、高正生成热以及感度低、稳定性好等特点.DAAT预测生成热高达 1000 kJmol-1 8,实验测量值在3337万方数据R 1NC 1R 1=N H C 2R 1=N H N H C 3:R 1 _ N。,C 4:R 1=N H N O C 5:R 1=N F,NR 2=N O R 2=N O R 2=N O R 2=N O R
9、 2=N O 一弋NVol.26Acta Phys.Chim.Sin.2010862-1034 kJmol-1 9-10之间,实验测量爆速为 7.40km s-1(=1.64 g cm-3)、爆压为24.00 GPa(=1.61 g cm-3)12.Wilcox等8对DAAT以及含有NO2取代3,3-偶氮-1,2,4,5-四嗪衍生物的生成热和分子结构进行了理论研究.虽然DAAT已经表现出较好的综合特性,但同时也存在密度较低和氧平衡偏低等不足.尽管DAAT已经成功合成并已有一些相关报道,但一些爆轰性能以及热力学性能仍然缺乏,另外其它类似衍生物也尚未合成.因此,迫切需要对DAAT及其类似衍生物的设
10、计和合成进一步研究,寻找具有高爆轰性能、较高稳定性和高密度的品优3,3-偶氮-1,2,4,5-四嗪类HEDMs.由于一些化合物的合成较为困难且具有一定危险性,因此在设计HEDM方面,理论预测就成为行之有效的方法.理论研究不仅为筛选候选化合物提供了可能性,而且可以揭示其结构与性能之间的关系,因此可以为实验设计提供更多、更有效的基础性数据13-14.本文报道了系列3,3-偶氮-1,2,4,5-四嗪衍生物的气相生成热、热稳定性以及含能性质.首先运用密度泛函理论(DFT)15-18,通过设计的等键反应预测标题化合物的生成热,然后利用键离解能评估了化合物的热稳定性,最后在计算的生成热和密度基础上预测了它
11、们的爆速和爆压.希望这些结果能为设计理想的HEDM提供必要信息,并且能够对实验合成以及性能研究起到实际的指导意义.1计算方法运用DFT方法结合设计的等键反应,在B3LYP/6-311G*和B3P86/6-311G*水平上预测了标题化合物的生成热.先前的报道4表明6-311G*基组的计算结果能够较好地符合实验值.基于总能量计算的等键反应方法已成功应用于预测化合物的生成热4-5,19-25.在等键反应中,因反应物和生成物电子环境相近,由电子相关能造成的误差可以相互抵消,使求得的生成热误差大为降低23.在设计的等键反应中1,2,4,5-四嗪的基本环骨架未破坏,选取1,2,4,5-四嗪和CH3NNCH
12、3为参考物质,降低了计算误差,该方法已被证明非常可靠4,5,19-21.目标化合物被分成3组(A、B和C三个系列),见图1.先前的工作已初步判定带有NO2的衍生物其性能较为理想4-5,因此C系列衍生物一端取代基定为NO2,用以探讨当另一侧引入其它基团后对其结构以及性能的影响.在298 K下,用于计算标题化合物生成热的等键反应设计如下:XnR+(n+2)CH42Tz+nCH3R+CH3NNCH3(n2)(1)XR1R2+4CH42Tz+CH3R1+CH3R2+CH3NNCH3(2)其中R(或R1,R2)=NH2、NHNH2、NO2、N3、NHNO2和NF2.X 表示 3,3-偶氮-1,2,4,5
13、-四嗪,Tz表示1,2,4,5-四嗪.对于等键反应(1)和(2),在298K时的反应热(H298 K)可表示如下:H298 K=Hf,p-Hf,r(3)式中Hf,r和Hf,p分别表示反应物和生成物在298 K时的生成热.参考物质 Tz、CH4、CH3NNCH3、CH3NH2、CH3NHNH2和CH3NO2的气相生成热是已知的26.由于CH3N3、CH3NHNO2和CH3NF2缺乏实验值,因此需要借助于G2方法计算精确生成热27-28.通过设计原子化反应(4)和(5),可以精确计算CH3N3和CH3NF2的生成热,而CH3NHNO2生成热的计算基于等键反应(6).CH3N3(g)C(g)+3H(
14、g)+3N(g)(4)图13,3-偶氮-1,2,4,5-四嗪衍生物的分子结构Fig.1Molecular frameworks of 3,3-azobis-1,2,4,5-tetrazine derivatives3338万方数据No.12张婧婧等:高能量密度材料3,3-偶氮-1,2,4,5-四嗪衍生物的分子设计CH3NF2(g)C(g)+3H(g)+N(g)+2F(g)(5)CH3NHNO2(g)+CH4(g)CH3NH2(g)+CH3NO2(g)(6)当已知反应热H298 K时,则可按照式(3)求得标题化合物的生成热.H298 K可以按照如下公式进行计算:H298 K=E298 K+(pV
15、)=E0+ZPE+HT+nRT(7)式中E0表示 0 K 时产物与反应物的总能量之差,ZPE为0 K时产物与反应物的零点能(ZPE)之差;HT是 0-298 K 的温度校正项.在理想气体反应中,(pV)=nRT.对于等键反应(1)和(2),n=0,因此(pV)=0.化学键的强度可用键离解能(BDE)来评价,它对于了解化学反应过程十分重要29.在 UB3LYP/6-311G*水平上计算标题化合物的键离解能,用以对比化合物的键强度以及热稳定性.在 298 K 和101.325 kPa 下化学键均裂的能量等于反应 AB(g)A (g)+B (g)的反应焓,该反应焓可以定义为AB的键离解焓30.文献报
16、道许多有机化合物的键离解能(BDE)和键离解焓通常是近似相等的31.因此,在0 K下物质的键离解能可以依据式(8)进行计算31:BDE0(AB)=E0(A )+E0(B )-E0(AB)(8)运用式(9)可以计算零点能(ZPE)校正的键离解能:BDE(AB)ZPE=BDE0(AB)+ZPE(9)其中ZPE表示生成物与反应物的零点能之差.运用Kamlet和Jacobs于1968年提出的半经验K-J方程32估算其爆速(D)、爆压(p)值:D=(1.011+1.312)1/2(10)p=1.5582(11)=NM1/2Q1/2(12)式中p为爆压(GPa),D为爆速(km s-1),为密度(g cm
17、-3),N为每克炸药爆轰生成气体的摩尔数(molg-1),M为气体产物的平均摩尔质量(gmol-1),Q为每克炸药的爆轰化学能(J g-1),结合等键反应所得生成热后,参照文献32的方法求得Q.计算化合物的密度必须已知摩尔体积,基于0.001 e bohr-3的等电子密度面所包围的体积空间,应用Monte Carlo方法求算摩尔体积,这种方法已经成功应用于高氮化合物33.对标题化合物在B3LYP/6-31G*理论水平上进行几何优化获得稳定构型,对每一个优化稳定构型执行100次单点计算取平均值.计算方法为B3LYP/6-311G*和B3P86/6-311G*,所有计算使用Gaussian 09量
18、子化学软件包34完成.收敛精度取程序默认值,振动分析表明所有得到的优化构型在此计算水平上没有虚频率,都是势能面上的局域极小点.2结果与讨论2.1生成热生成热是化合物的基本热力学性质,是其含能高低的标志,也是衡量 HEDM 爆轰性能的重要参数,因此精确计算生成热对设计、筛选以及合成新型品优HEDM具有重要意义.一些高氮化合物标准生成热数据仍然缺少文献值或不易查找,其标准生成热的估算无论在理论上还是在实验中均具有重要意义.本文通过设计等键反应对标题化合物的生成热进行了计算.图1给出了3,3-偶氮-1,2,4,5-四嗪衍生物的分子结构.表1列出等键反应(1)和(2)中9种参考物质的在B3LYP/6-
19、311G*和B3P86/6-311G*水平Compd.TzCH4CH3NNCH3CH3NH2CH3NHNH2CH3NO2CH3N3CH3NHNO2CH3NF2B3LYP/6-311G*E0(a.u.)-296.390906-40.533744-189.328011-95.888439-151.222209-245.081673-204.148401-300.434462-294.298331ZPE(kJ mol-1)134.40117.09220.62167.56213.44130.50131.67176.48122.75HT/(kJ mol-1)13.8410.0316.2911.5613.
20、7414.1514.4116.3413.85B3P86/6-311G*E0(a.u.)-297.134547-40.713980-189.900331-96.213717-151.692042-245.649182-204.679786-301.147794-294.895029ZPE(kJ mol-1)135.41117.42221.53168.31214.55131.51132.53177.78123.62HT/(kJ mol-1)13.8710.0416.2111.5513.7013.9914.3916.3013.75HOF(kJ mol-1)464.00b-74.60b141.00b-
21、22.50b94.70b-80.80b296.54c-15.20c-115.23c表1计算的参考物质总能量(E0)、零点能(ZPE)、温度校正项(HT)和生成热(HOFs)aTable 1Calculated total energies(E0),zero-point energies(ZPE),thermal corrections(HT),and heats of formation(HOFs)of the reference compoundsaaThe scaling factor for ZPE is 0.98 and the scaling for HTis 0.9635;bDat
22、a are the experimental values taken from Ref.26;cThe values are calculated at the G2 level.3339万方数据Vol.26Acta Phys.Chim.Sin.2010上计算的总能量、零点能和温度校正项.对振动频率以校正因子0.96进行校正(参考Scott和Radom的研究结果35).在表1中参考物质的气相实验生成热参考文献26(包括 Tz、CH4、CH3NNCH3、CH3NH2、CH3NHNH2和CH3NO2).CH3NF2和CH3N3的精确生成热分别参考文献计算值22,25.另外,在G2水平上运用式(7
23、)计算了反应(6)的等键反应反应热(H298 K)27-28.CH4(g)、CH3NH2(g)和CH3NO2(g)的生成热实验值是已知的26,故可以计算 CH3NHNO2的精确生成热(Hf).为了验证理论计算的可靠性,我们在G2水平上计算了NH2NO2的生成热.结果显示,通过设计等键反应:NH2NO2(g)+CH4(g)NH3(g)+CH3NO2(g),G2方法计算的生成热为-3.20 kJmol-1,与实验值(-3.20 kJmol-1 20)完全一致.由于 CH3NHNO2和NH2NO2具有相似的结构特点,说明G2方法计算的CH3NHNO2生成热应该是可靠的.表2列出3,3-偶氮-1,2,
24、4,5-四嗪衍生物计算的总能量、零点能、温度校正项以及生成热.结果显示,B3P86/6-311G*和B3LYP/6-311G*这两种方法计算的生成热结果十分接近,具有很好的线性关系:HOFB3P86=1.0003HOFB3LYP-6.7479,R2=0.9997.这说明两种方法对标题化合物生成热的计算结果都较为理想.两种方法对于B1(DAAT)的HOF计算结果与实验值11都比较接近,前者略大于后者约50 kJ mol-1且在误差范围内(48 kJ mol-1),见表2.同时与先前报道的B1理论计算结果1148.51 kJ mol-1 8相比,我们计算的结果更加接近实验值.这也表明了两种方法用以
25、预测标题化合物的生成热数据是可靠的.从表2可见,所有的3,3-偶氮-1,2,4,5-四嗪衍生物都具有高正生成热,超过1000kJ mol-1,其中B4具有最高的生成热大约1830 kJ mol-1.符合先前的研究报道,含能高氮杂环化合物具有高正生成热1-2.标题化合物的生成热数值远高于我们之前研究的1,2,4,5-四嗪衍生物4,原因可能为前者分子结构中含有偶氮基以及两个四嗪环.NHNH2、NO2、N3或NHNO2取代的衍生物生成热比未取代的3,3-偶氮-1,2,4,5-四嗪生成热高,而NH2或NF2取代的衍生物生成热相对较低.在带有相同取代基数量的情况下,带有N3取代基的衍生物生成热最大,这与
26、先前的研究结果一致:N3基团可以显著增加分子的生成热,是已知含能材料官能团中增加生成热显著的一类取代基,大约增加 293-320 kJmol-1 4-5,9,36.通过对比单取代和双取代衍生物的生成热,我们发现3,3-偶氮-1,2,4,5-四嗪衍生物的生成热并不符合简单的基团加和规律.2.2热稳定性键级是基于分子轨道系数进行Mulliken布居数分析所得的重要电子结构参数.一般而言,原子间键级越大,表明其间电子云密度较大或重叠较多,预示表2计算的3,3-偶氮-1,2,4,5-四嗪衍生物计算的总能量(E0)、零点能(ZPE)、温度校正项(HT)和生成热(HOFs)aTable 2Calculat
27、ed total energies(E0,),zero-point energies(ZPE),thermal corrections(HT),and heats of formation(HOFs)of the 3,3-azobis-1,2,4,5-tetrazine derivatives at the B3LYPand B3P86 levelsaCompd.AA1A2A3A4A5A6B1B2B3B4B5B6C1C2C3C4C5B3LYP/6-311G*E0(a.u.)-701.045715-756.442385-811.773000-905.578657-864.678942-960.9
28、50646-954.797844-811.836975-922.498331-1110.112202-1028.312054-1220.854751-1208.549265-960.977999-1016.302184-1069.213191-1165.484393-1159.331006ZPE(kJ mol-1)239.42283.74329.05243.68247.11289.26235.82327.74418.74248.38254.82338.68232.05288.44332.77251.63293.51240.36HT(kJ mol-1)30.6634.5438.2637.9537
29、.5640.6938.3838.5445.8145.0644.4551.0146.2741.6045.3444.7447.9945.59HOF(kJ mol-1)1200.571138.781265.261228.001518.721243.281187.471082.271335.091254.101837.181288.021176.241159.401301.821542.861268.591214.52B3P86/6-311G*E0(a.u.)-702.745830-758.289307-813.764950-907.666232-866.731179-963.184447-956.9
30、15433-813.830611-924.782021-1112.587199-1030.716396-1223.622082-1211.084383-963.2124065-1018.681585-1071.652870-1168.105706-1161.836066ZPE(kJ mol-1)241.86286.68332.38246.86250.23292.77238.96331.17422.95252.27258.60343.21235.87292.10336.95255.50297.77244.17HT(kJ mol-1)30.6334.3938.1337.8437.4540.5838
31、.2338.2945.5644.9144.2750.8646.0641.4145.0844.5647.7845.44HOF(kJ mol-1)1198.161131.731257.021225.171513.851239.381182.551070.811321.251250.951829.891283.041168.701151.871293.241537.641264.131209.16HOF(kJ mol-1)103548b1148.51c1314.86c1220.85caThe scaling factor for ZPE is 0.98 and the scaling factor
32、for HTis 0.9635;bData are the experimental values taken from Ref.11;cThe theoretical HOFs were taken from Ref.8.3340万方数据No.12张婧婧等:高能量密度材料3,3-偶氮-1,2,4,5-四嗪衍生物的分子设计该键较强、不易断裂.表3给出标题化合物分子的CN(四嗪环)、NN(四嗪环)、N7N8、CR(R=NH2、NHNH2、NO2、N3、NHNO2、NF2)以及NR(R=NH2、N2、NO2)的重叠布居数.结果显示,CR 或 NR 键的键级大都比 CN、NN和N7N8键的键级小(衍
33、生物带有NR 键时不比较 CR 键),只有 A1、B1 和 B2 的 CR 或NR 键的键级比N7N8键级大.A3、A6、B3、B6以及C1C5的CR键级在所有化学键中最小,表明在这些衍生物中CR键较弱,在热解或撞击下可能最易断裂,是热解引发键.A2、A4、A5、B4和B5的NR 键具有最小的布居数,因此热解时可能生成R 自由基.然而,A1、B1和B2的N7N8键级最小,表明热解时连接两个四嗪环的偶氮键可能作为引发键发生断裂.我们发现,大多数衍生物NN键的键级比 3,3-偶氮-1,2,4,5-四嗪大,除了 A3、A6、B3、B4、B6、C3和C5,说明取代基的引入可以增加3,3-偶氮-1,2,
34、4,5-四嗪NN键的强度.在所有标题化合物中,A1最弱键的键级(0.2615)最大.根据判别撞击感度相对大小的“最小键级原理(PSBO)”19,预示A1的感度最小、热力学稳定性最好.选择 Mulliken 布居数最小的 CR、NN 或NR 键 作 为 可 能 的 热 解 引 发 键,用 UB3LYP/Compd.AA1A2A3A4A5A6B1B2B3B4B5B6C1C2C3C4C5PCN0.46400.46270.46060.45480.45940.45750.45580.46400.45960.44680.45510.45330.45040.45260.45230.45020.44850.4
35、475PNN0.28390.29610.29460.28130.28400.29050.28190.30740.30460.28130.28370.29660.28050.29440.28960.28250.28960.2813PN7N80.27380.26150.26020.27940.26660.26740.27140.25220.24940.27360.25990.26200.26980.25960.26040.26610.26700.2725PCR0.41040.38560.36040.14610.40850.27320.21890.38180.35760.14670.40700.27
36、230.22100.14490.14630.14710.14750.1512PCR0.38890.35550.41170.27410.2200PNR0.2601(NHNH2)0.2068(NN2)0.1617(NHNO2)0.2609(NHNH2)0.2075(NN2)0.1576(NHNO2)0.2605(NHNH2)0.2035(NN2)0.1558(NHNO2)表33,3-偶氮-1,2,4,5-四嗪衍生物在B3LYP/6-311G*水平上计算的部分键平均重叠布居数(P)Table 3Calculated average bond overlap populations(P)of part
37、 bonds of the 3,3-azobis-1,2,4,5-tetrazinederivatives at the B3LYP/6-311G*levelCompd.AA1A2A3A4A5A6B1B2B3B4B5B6C1C2C3C4C5RDXHMXPNN0.28390.29610.28400.30740.2837PN7N80.27380.26150.25220.2494PCR0.14610.21890.14670.22100.14490.14630.14710.14750.1512PNR0.2601(NHNH2)0.2068(NN2)0.1617(NHNO2)0.2609(NHNH2)0.
38、2075(NN2)0.1576(NHNO2)0.1558(NHNO2)0.1376(NNO2)0.1382(NNO2)BDE0(kJ mol-1)302.24(409.46)b335.80(401.98)b262.59241.24325.91(735.51)b105.11270.22318.59(389.05)b262.29(387.93)239.34322.22(733.40)b103.27267.14244.46245.06242.0299.61(241.20)b240.11166.19178.77BDEZPE(kJ mol-1)283.25(384.69)b318.90(378.12)b
39、230.90224.10309.02(709.19)b86.82253.39304.23(366.39)b230.36(364.66)222.26305.89(706.97)b85.43250.60227.10227.83224.8581.85(224.17)b222.95145.62160.41aBDE0denotes the bond dissociation energies without zero-point energy corrections,while BDEZPEdenotes the bond dissociation energies includingzero-poin
40、t energy corrections,bThe values in parentheses are the bond dissociation energies of the bonds with the smallest bond overlap populations.表4在B3LYP/6-311G*水平上计算的3,3-偶氮-1,2,4,5-四嗪衍生物以及RDX和HMX的最弱键重叠布居数以及键离解能(BDE)aTable 4Calculated bond dissociation energies(BDE)and bond overlap populations of the weak
41、est bonds for the3,3-azobis-1,2,4,5-tetrazine derivatives together with RDX and HMX at the B3LYP/6-311G*levela3341万方数据Vol.26Acta Phys.Chim.Sin.20106-311G*方法计算这些键断裂所需的离解能(BDE),标题化合物最弱键的重叠布居数以及键离解能(BDE)列于表4.结果表明没有零点能校正的键离解能(BDE0)通常比校正的键离解能(BDEZPE)大,但零点能不影响化合物的稳定性次序.当3,3-偶氮-1,2,4,5-四嗪被NH2和N3基团取代后BDEZPE
42、有所增加,而其它衍生物的BDEZPE比未取代的3,3-偶氮-1,2,4,5-四嗪小.计算的BDE值可以用来比较含能材料的相对热稳定性顺序37-38.因此,可以推断出NH2和N3基团有利于增加3,3-偶氮-1,2,4,5-四嗪的热稳定性.通过与常用炸药RDX(黑索今,环三亚甲基三硝胺)以及HMX(奥克托今,1,3,5,7-四硝基-1,3,5,7-四氮杂环辛烷)比较,我们发现大多数衍生物具有较高的 BDEZPE,这表明它们具有较为理想的热稳定性.分析结果发现A4的NR 具有相对较低的键级(0.2068),但在所有取代3,3-偶氮-1,2,4,5-四嗪衍生物的最弱键中,A4具有最大的BDEZPE值(
43、709.19kJmol-1).然而,虽然A4四嗪环上NN键具有较高的键级(0.2840),但却具有相对较低的BDEZPE值(309.02 kJmol-1).B4的计算结果与A4相似,表明A4和B4在受热或撞击下应该发生四嗪环断裂.由此说明,在判定标题化合物热稳定性时不能只是简单依据化学键的布居数,有必要借助键离解能来判定.从表4可以看出除了A1和B1,B系列衍生物的BDE与A系列衍生物的BDE非常接近,说明增加取代基数量对衍生物的稳定性影响较小.同时与之前研究的单环1,2,4,5-四嗪衍生物结果4相比较,我们发现相同取代基对于不同体系稳定性的影响有所不同,在单环1,2,4,5-四嗪衍生物中其稳
44、定性随着取代基数目的增加而增加.另外,B1四嗪环上NN的键级(0.3074)比N7N8的键级(0.2522)高,但前者的BDE比后者低,说明B1热解时首先发生四嗪环断裂,与先前的报道相符11.同时说明在键离解能数据没有相关实验值的情况下,其数值与传统炸药在相同计算水平UB3LYP/6-311G*的所得数据相比,可用以判定新设计的偶氮四嗪衍生物的热稳定性是否理想,先前的工作已证明该计算方法是可行的4-5.2.3爆轰性能含能材料的爆速和爆压是两个重要的性能参数,预测的爆速(D)和爆压(p)列于表5.为了进行比较,已知爆轰性能实验值的RDX和HMX也列于表539.先前的研究表明4,5,40-42,将
45、半经验的K-J方程应用于预测高氮化合物的爆轰性能是可靠的.Compd.AA1A2A3A4A5A6B1B2B3B4B5B6C1C2C3C4C5RDXbHMXbV/(cm3 mol-1)114.12126.38136.56132.81140.50142.39133.66134.31151.80153.32161.98173.48149.84144.66153.86152.70160.12150.70124.91157.50N/(mol g-1)0.030.040.040.030.030.030.030.030.040.030.030.030.030.030.030.030.030.030.030.
46、03M/(g mol-1)23.6721.5020.2032.0028.0029.4728.4422.4321.5029.0426.1427.7327.0327.7327.1330.0030.6330.1227.2027.20Q/(J g-1)1506.791162.931262.661739.461607.591630.511771.461319.061365.161668.931565.951603.871688.601520.231602.261672.071679.891755.241597.391633.88/(g cm-3)1.671.621.611.771.641.761.801
47、.64(1.6412)1.651.831.681.791.951.731.721.811.841.901.78(1.8239)1.88(1.9139)D/(km s-1)7.817.617.808.387.858.348.567.58(7.4012)8.028.858.038.709.398.148.218.628.899.158.87(8.7539)9.28(9.1039)p/GPa25.7624.0825.2030.8525.8330.3932.5624.04(24.1012)27.0135.0327.4233.4240.9428.6729.0933.0535.5938.3534.67(3
48、4.0039)39.19(39.0039)aData in the parentheses are the experimental values taken from Refs.12,39;bThe calculated values of RDX and HMX are taken from Ref.4;V:average volume from 100 single-point calculations at the B3LYP/6-31G*level,N:the moles of detonation gases per gram explosive,M:theaverage mole
49、cular weight of these gases,Q:the heat of detonation(J g-1),:the loaded density of explosives(g cm-3),D:the detonation velocity(km s-1),p:the detonation pressure(GPa)表53,3-偶氮-1,2,4,5-四嗪衍生物以及RDX和HMX的爆轰性能aTable 5Predicted detonation properties of the 3,3-azobis-1,2,4,5-tetrazine derivatives together w
50、ith RDX and HMXa3342万方数据AA1A 2 A 3 A 4 A 5 A 6B1B 2B3B 4B5B 6C1C 2C3C 4C5HM XMolecularnumber ing5O5O9988(L-E邑Q5O5O574332lB卫oCyOOOOOO2O5O5O32211(|-|oE王一山QNo.12张婧婧等:高能量密度材料3,3-偶氮-1,2,4,5-四嗪衍生物的分子设计如表5所示,B1、RDX以及HMX的理论预测值(、D和p)与实验值12,39比较接近,说明该方法用于预测标题化合物的爆轰性能是相对可靠的.结果显示,带有不同取代基团的3,3-偶氮-1,2,4,5-四嗪衍生物其密