不同模型下HMX晶体结构和性能的MD研究.pdfVIP

不同模型下 HMx晶体结构和性能的MD研究 765 文章编号 ：1006—9941(2013)06—0765—06 不同模型下 HMX晶体结构和性能的MD研究 刘冬梅 ，肖继军，陈 军 ，姬广富，朱 伟。，赵 锋 ，吴 强，肖鹤鸣 (1．南京理工大学化工学院 分子与材料计算研究所，江苏 南京 210094；2．中国工程物理研究院流体物理研究所 ，四川 绵阳 621900；3．嘉兴学院生物与化学工程学院，浙江 嘉兴 314001) 摘 要 ：用分子动力学 (MD)方法 ，在 COMPASS力场和 295K—NPT系综下，对奥克托今 (HMX)(4×2×4)和(4X4X4)超 晶胞及 其沿各 自(100)晶面构成 的共 6种模型进行 了模拟 ，求得它们 的平衡结构、引发键键长分布、引发键连双原子作用能和弹性力学 性能。结果表明，引发键呈近似对称的高斯型分布，平均键长(L…)与实验值很接近；随着体系中原子数 目增加，引发键连双原子 作用能 ￡一 单调增加 ；除柯西压外，各模型的力学性能比较稳定一致 。 关键词 ：物理化学 奥克托今 (HMx)；分子动力学 (MD)模拟 ；模型；结构 ；相互作用能；力学性能 中图分类号 ：T／55；O64 文献标识码 ：A DOI：10．3969／j．issn．1006-9941．2013．06．015 原子作用能)和力学性能及其对感度的关联等方面 ， 1 引 言 提供模型对 MD模拟结果影响的基础信息，希望可对 后续高能材料结构和性能 (特别是安全性和力学性 奥克托今 (HMX，CH。O。N。)是当前综合性能最 能)的MD模拟工作有所助益。 好的单体炸药，已作为重要组分广泛应用于混合炸药 和固体推进剂等复合含能材料 中。国内外已有较多关 2 模型搭建和模拟细节 于 HMx晶体 和以HMX为基混合炸药如高聚物 粘结炸药 (PBX) 。‘等的分子动力学 (MD)模拟研 2．1 模型搭建 究。但其 中专 门讨论建立模型的文献很少 。肖继军 以口．HMX中子衍射获得的晶体数据 ¨为依据，搭 等。。’’在 PBXs的MD模拟中，曾经建议过 “吸附包 建 HMX晶胞模 型 。选取 (4×2X4)超 晶胞 ，含 64 个 HMX分子 ，共 1792个原子 ，由此得到的纯 HMX周 覆”、“渗透添加”和 “切割分面”三种模型，主要用于研 期箱记为模型I；对 (4×2×4)超晶胞沿 (100)晶面方 究 TATB和 HMx等基炸药与氟聚物之间的结合能和 力学性能 。而对于基炸药晶体模 型的讨论 ，特别是其 向进行 “切割”，并使 z轴平行晶轴C向量 ，同时垂直于 (100)晶面，晶面 (100)面积为 2．21×3．48nm ，c方 模型选取对涉及安全性 (如感度)判别 的MD模拟结 向上取深度 1即OC=2．38nm其上真空层高度设置 果 的影响问题 ，至今 尚未见有报道 。因为力场选定之 为0，这样得到的晶体模拟周期箱 ，记作模型 Ⅱ(原子 后 ，影响MD模拟结果的关键 因素就是模型建立 ，故 数亦为 1792)。用 同样 的方法在 c方向取深度 1．5即 本工作选择 HMx晶体为研究对象，通过 MD模拟 ，对 OC：3．46nm，亦得到另一种晶体模拟周期箱 ，记作模