一元溶剂体系β-HMX球形化结晶形貌的分子动力学模拟

陈 芳, 周 涛

(中北大学化学工程与技术学院,太原 030051)

1 引 言

奥克托金(HMX)为典型的八元杂环硝胺类猛炸药, 因其密度高、熔点高、爆轰性能优良和热安定性好, 使其成为国内外现用炸药中综合性最好的单质炸药.奥克托金有四种晶型, 其中β-HMX在常温常压下是最稳定的. 有关β-HMX晶体的弹性性质[1]、内部缺陷[2]、热膨胀特性[3]、高压性能[4]、热分解[5]已有相关研究.炸药的晶体形貌在很大程度上影响其安定性能、流散性和能量输出, 对感度的影响尤为明显. 研究发现, 棱角分明的炸药晶体在机械刺激下尤其敏感, 立方体状或球形化炸药晶体机械感度较低[6]. 影响晶体晶习主要有晶体的内部结构和外部因素, 在外部因素中, 溶剂对晶习的影响起着很关键的作用[7]. 实验和计算模拟是研究结晶形貌最主要的方法, 而实验研究是比较昂贵和耗时的, 目前基于附着能模型(attachment energy, AE)以及修正附着能模型(modified attachment energy, MAE)[8-17], 应用分子动力学模拟方法成功的预测了含能材料在溶剂中的结晶行为.β-HMX几乎不溶于水, 但极易溶解在有机溶剂中. 应用修正的附着能模型(modified attachment energy, MAE)对β-HMX在乙腈(ACN)[10]和丙酮(AC)溶剂[16]中的结晶形貌进行了理论预测, 与实验结果相吻合. 但有关N,N-二甲基甲酰胺(DMF)、环己酮(CH)、γ-丁内酯(GBL)、二甲基亚砜(DMSO)、N-甲基吡咯烷酮(NMP)五种溶剂作用下β-HMX结晶形貌的变化规律没有被系统的研究.

为了能够更全面的从理论上剖析β-HMX的晶体生长过程, 本文运用分子动力学方法(MD)模拟 DMF、CH、GBL、DMSO、NMP五种溶剂对β-HMX各晶面的影响, 采用修正的附着能模型……