TNAZ的第一原理分子力学力场.pdfVIP

维普资讯 第 l5卷 第 6期 含 能 材 料 Vo1．15，No．6 2007年 12月 CHINESEJOURNALOFENERGETICMATERIALS December，2007 文章编号 ：1006—9941(2007)06—0616-06 TNAZ的第一原理分子力学力场 王喜军，赵立峰，孙 淮 (上海交通大学化学化工学院，上海 200240) 摘要：分子力学力场的适用性和准确性是制约分子模拟应用 的瓶颈。本文依照 TEAM力场 的研发思想，推导了 1，3，3-三硝基氮杂环丁烷 (TNAZ)的全原子第一原理分子力学力场。键参数和电荷参数拟和 自第一性原理计算结果， 范德华参数迁移 自文献报道并通过模拟 TNAZ液相密度和蒸发焓进行了优化。通过分子力学、分子动力学模拟，计算 了TNAZ分子的气相、液相和晶相性质对参数进行了验证。验证结果显示气相分子的结构、能量、振动性质，液相的 密度、蒸发焓，晶体结构、晶格能等结果与实验测定一致。同时预测了TNAZ晶体的热胀性质、等温压缩等性质。 模拟得到的热胀 曲线在整个凝聚相 内显示出三个具有不同热胀率的区段。其中300K下TNAZ晶体体积较 400K 下减小约 13％，与实验观测吻合。300K下等温压缩到 10GPa，体积减小约28％，与文献报道的理论值一致。 关键词：物理化学；分子力学力场 ；分子动力学；1，3，3-三硝基氮杂环丁烷(TNAZ) 中图分类号：TJ55；TQ015．9；064 文献标识码：A 直接从量子化学计算结果推导其力场参数，避免在不 1 引 言 同的分子间迁移参数，大大提高力场参数的准确性。 含能材料的研究开发对国防、航天领域的发展有着 其次，采用数据库的方式存储大量中小分子的参数，对 不言而喻的重要性。可是，含能材料的实验研究成本昂 于大分子或高分子化合物，则将 已有的小分子或分子 贵、环境不友好且周期长。而采用理论计算方法辅助含 片断力场参数组合，得到大分子或高分子的力场。第 能材料的研究不仅能够揭示材料性能的微观本质，而且 三，范德华参数一般不能直接从量子化学计算得到，须 有可能有效地降低成本，缩短周期。将分子模拟用于预 通过 拟合大量分子 的凝聚相实验数据得到，因而 测含能材料的多晶形 ¨ 、熔点 、导热系数 、剪切 TEAM力场采用一套统一优化的范德华参数。 粘度 等性质的研究已有文献报道，这些研究对揭示含 表 1 已经报道 的含能材料 的分子力学力场研究 能材料的感度、初始爆轰反应的机理有着重要意义。 Tablel Reported forcefield researches011energeticmaterials 然而，从 目前对含能材料的模拟研究报道来看 ，存 在一个突出的问题，即缺乏准确、可靠和普适的分子力 学力场。表 1列出了文献 中部分含能材料的分子力学 力场 …。从表 中可 以看 出：已经报道的含能材料 的力场函数形式不统一 ，没有可迁移性 ，缺乏系统的对 含能材料的分子力学研究。 近年来，本课题组着力研究系统解决力场参数的 准确性、迁移性和可扩展性问题 ，提出了建立第一原理 传统的力场基于从一组分子的训练集获得参数， TEAM(Transferable，Extensible，Accurate，Manageable， 因此其力场