不同模型下HMX晶体结构和性能的MD研究-含能材料.PDFVIP

不同模型下 ＨＭＸ晶体结构和性能的ＭＤ研究 ７６５ 文章编号：１００６９９４１（２０１３）０６０７６５０６ 不同模型下 ＨＭＸ晶体结构和性能的ＭＤ研究 １ １ ２ ２ ３ ２ ２ １ 刘冬梅 ，肖继军 ，陈　军 ，姬广富 ，朱　伟 ，赵　锋 ，吴　强 ，肖鹤鸣 （１．南京理工大学化工学院分子与材料计算研究所，江苏 南京 ２１００９４；２．中国工程物理研究院流体物理研究所，四川 绵阳 ６２１９００；３．嘉兴学院生物与化学工程学院，浙江嘉兴 ３１４００１） 摘　要：用分子动力学（ＭＤ）方法，在ＣＯＭＰＡＳＳ力场和２９５ＫＮＰＴ系综下，对奥克托今（ＨＭＸ）（４×２×４）和（４×４×４）超晶胞及 其沿各自（１００）晶面构成的共６种模型进行了模拟，求得它们的平衡结构、引发键键长分布、引发键连双原子作用能和弹性力学 性能。结果表明，引发键呈近似对称的高斯型分布，平均键长（Ｌ ）与实验值很接近；随着体系中原子数目增加，引发键连双原子 ａｖｅ 作用能 Ｅ 单调增加；除柯西压外，各模型的力学性能比较稳定一致。 Ｎ—Ｎ 关键词：物理化学奥克托今（ＨＭＸ）；分子动力学（ＭＤ）模拟；模型；结构；相互作用能；力学性能 中图分类号：ＴＪ５５；Ｏ６４ 文献标识码：Ａ ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１００６９９４１．２０１３．０６．０１５ 原子作用能）和力学性能及其对感度的关联等方面， １　引　言 提供模型对 ＭＤ模拟结果影响的基础信息，希望可对 后续高能材料结构和性能（特别是安全性和力学性 　　奥克托今（ＨＭＸ，ＣＨＯ Ｎ）是当前综合性能最 ４ ８ ８ ８ 能）的 ＭＤ模拟工作有所助益。 好的单体炸药，已作为重要组分广泛应用于混合炸药 和固体推进剂等复合含能材料中。国内外已有较多关 ２　模型搭建和模拟细节 于 ＨＭＸ晶体［１－６］和以 ＨＭＸ为基混合炸药如高聚物 ［７－１４］ ２．１　模型搭建 粘结炸药（ＰＢＸ） 等的分子动力学（ＭＤ）模拟研 　　以 ＨＭＸ中子衍射获得的晶体数据［１７］为依据，搭 究。但其中专门讨论建立模型的文献很少。肖继军 β ［１８］ 等［７，１５，１６］在 ＰＢＸｓ的ＭＤ模拟中，曾经建议过“吸附包 建 ＨＭＸ晶胞模型 。选取（４×２×４）超晶胞，含６４ 覆”、“渗透添加”和“切割分面”三种模型，主要用于研 个 ＨＭＸ分子，共 １７９２个原子，由此得到的纯 ＨＭＸ周 期箱记为模型 ；对（４×２×４）超晶胞沿（１００）晶面方 究 ＴＡＴＢ和 ＨＭＸ等基炸药与氟聚物之间的结合能和 Ⅰ 力学性能。而对于基炸药晶体模型的讨论，特别是其 向进行“切割”，并使Ｚ轴平行晶轴 ｃ向量，同时垂直于 模型选取对涉及安全性（如感度）判别的 ＭＤ模拟结 （１００）晶面，晶面（１００）面积为２．２１×３．４８ｎｍ２，ｃ方 向上取深度１即ＯＣ＝２．３８ｎｍ其上真空层高度设置 果的影响问题，至今尚未见有报道。因为力场选定之 为０，这样得到的晶体模拟周期箱，记作模型 （原子