一种化合物晶体密度的预测方法技术

本发明专利技术公开了一种化合物晶体密度的预测方法，属于化合物晶体密度的计算方法技术。预测方法包括以下步骤：（１）应用量子化学计算软件，采用密度泛函Ｂ３ＬＹＰ／６－３１Ｇ＊方法优化化合物的可能构型的分子结构，确定能量最低的分子构型，得到该分子结构能量最高的电子占有轨道的能量ＥＨＯＭＯ和能量最低的空轨道的能量ＥＬＵＭＯ，然后采用Ｍｏｎｔｅ－Ｃａｒｌｏ方法计算该分子构型的摩尔体积，每一个分子结构进行５０－３００次，计算平均值得到摩尔体积Ｖ；（２）计算化合物理论密度ρｃａｌ为分子摩尔质量Ｍ除以平均摩尔体积Ｖ；（３）预测化合物晶体最大密度ρｅｓｔ＝０．２３２＋０．９７０ρｃａｌ＋０．５０２ＥＨＯＭＯ＋０．２２４ＥＬＵＭＯ，预测密度与实测密度的最大相对误差约为４％。

【技术实现步骤摘要】

本专利技术涉及化合物晶体密度的计算方法
，具体地说是一种化合物晶体密 度的预测方法。
技术介绍
1900年Planck提出量子论。1926年Schrodinger方程问世，标志着量子力学的建 立。1927年Heitler和London用量子力学基本原理研究氢分子结构，标志着量子化学新学 科的诞生。随着量子化学理论方法的发展和计算机技术的进步，80多年来，量子化学已渗透 应用到各学科各领域，在研究原子、分子和晶体结构以及解释物质结构与性质的关系方面， 取得了举世瞩目的成就。近年来发展迅速的密度泛函理论(DFT)方法因较方便、快速地处理电子相关问题 而得到极为广泛的关注和应用。DFT建立在两个著名定理的基础上。它们通常称为HK定 理,是在1964年Hohenberg和Kohn在研究均勻电子气Thomas-Fermi (TF)模型的理论基础 时提出的。根据HK定理，原子、分子和固体的基态性质可用电子密度加以描述。据此，Kohn和 Sham又导出了 DFT中的单电子自洽场方程，通常称为KS方程。HK定理和KS方程奠定了 DFT方法的坚实基础。对较大的体系，DFT耗时比传统的 超HF从头计算要少1 2个数量级。它可以处理有机、无机、金属、非金属体系，几乎可以 囊括周期表中所有元素的化合物。常用的“杂化” DFT方法有B3LYP、B3P86和B3PW91等，其中B3LYP杂化方法使用 密度泛函理论Beck三参数杂化方法，结合Lee-Yang-Parr非定域相关泛函。此方法充分考 虑了电子相关，保持了从头计算MO法的很多有点，又较为省时。相对于一组精确的实验数 据G2测试组，B3LYP的误差绝对值仅略大于8kJ. mo Γ1。因此，近年来B3LYP方法在分子的 电子结构计算中获得了广泛的应用和空前的成功，是目前最为流行的DFT方法之一。现今最受关注的高能量密度材料(HEDM，High Energy Density Material)是20 世纪80年代-90年代出现的一类新型含能材料，具有能量密度高，安全性能好，低易损性， 低特征型号等做为炸药的优良性能。就HEDC能量水平而言，在宏观上目前一致得出的结论 是取决于三个因素晶体密度、标准生成焓和氧平衡。在设计新的HEDC时，应采用有利于提 高上述三个因素的分子组成和分子结构(构型和构象)。基于此点，有望成功地合成且具应 用前景的HEDC可归纳为下述三类(1)分子组成与HMX或黑索金(RDX)相近的多环笼型多 硝铵；(2)笼形多硝基化合物；(3)高氮杂环化合物(包括四氮杂并环戊二烯衍生物、四嗪 衍生物和呋咱及氧化呋咱系化合物等)。炸药的密度与很多爆炸性能(爆速、爆压、爆热、猛度及爆容等)有密切关系。对 碳氢氧氮系炸药，其爆速随密度线性增长，而冲击波前沿爆压正比于密度的平方。提高炸药 的密度是提高炸药能量水平的重要途径。化学工作者一直都在不懈努力寻找性能更好的化合物作为炸药材料，在开发一种具有良好炸药性能的化合物之前，都是先预测其各种性质，如DOS值，带隙，感度，键离解 能，密度值等，并将其与已有的炸药材料的相关性质进行比较，最后确定有无合成价值。密 度是衡量炸药材料性质的很重要的一个标准，因此预测炸药晶体密度就具有了重大意义。早在1996 年，俄国化学家 Α. V. Dzyabchenko, T. S. Pivina 和 E. A. Arnautova 就在 Journal ofMolecular Structure 第 378 卷第 2 期 67-82 页发表了题为《Prediction of structure and density fororganic nitramines))白勺文章。文章介绍了通过软件包优化， 利用从头计算法预测晶体结构的方法在硝胺炸药上得到了巨大的发展。硝胺炸药是一类很 重要的能量材料。这种研究方法的主要特征是以各类有机晶体结构的统计数据为基础，选 择典型的空间官能团和对称点进行进一步的研究；选择超曲面能量对称点，以此决定单一 的研究范围；搜索程序广泛的的收敛性使得优化可以从任意一点开始，这样初步筛选模型 就没有必要。数值计算首次是被HMX的三种构型的计算所证明。在这个步骤中，硝胺片段 的力场参数被证明是可以通过其得到预测和实际分子构型之间最好的相关性。预测所得的 与X衍射所得的结构参数可以很好的符合。预测的晶格能不能很准确的预测观察到的升华 热和极性稳定的趋势。其次，这种方法还用来预测的伍兹系列分子结构，其结构在之前的一 个量子研究上已经被提出。最终，金刚烷和伍兹系列被确定其多种晶体构型里面，很可能含 有最小储能结构密度会高达2. 08g/cm_3和2. 04g/cm_3的构型。由于分子间的构型和晶格能 的相互作用，总的能量有一定的增加，但是分子的集合对称性受到了一些损失，正是由于这 样的原因，计算密度应该比实际密度要低一些。西安近代化学研究所和太原理工大学应用化学系的来蔚鹏等人于2007年12月 在《含能材料》第15卷第6期626-629页发表了 “用定量结构性质关系(QSPR)预测芳香 系炸药的密度”。文中以物质的电子、空间等结构性质为基础，运用GauSSian98和Cerius2 程序包对偶极距(Dipole)、最高占据轨道能量(Emb)、最低空轨道能量(Euitro)、分子总能 量(E)、旋转键(Rotlbonds)、最弱的R-NO2键长(R-N02bondlength，R为C或N)、氢键供体 (Hbonddonor)和中点势(Vmid)8种描述符进行了计算，采用Cerius2程序包中的QSI3R方法 建立了芳香系炸药密度与8种描述符之间的构效关系式，30个化合物所构成的训练集和15 个化合物所构成的预测密度集与实测密度之间的平均误差分别为3. 33%和2. 94%。寻找通用的、简便的、预测密度更准确的方法具有广泛的需求和重要意义。
技术实现思路
本专利技术的目的是提供，包括以下步骤(1)应用 量子化学计算软件，采用密度泛函B3LYP/6-31G*方法优化化合物的可能构型的分子结构， 确定能量最低的分子构型，得到该分子结构能量最高的电子占有轨道的能量ΕΗ(Μ^Π能量最 低的空轨道的能量Euik，然后采用Monte-Carlo方法计算该分子构型的摩尔体积，每一个分 子结构进行50-300次，计算平均值得到摩尔体积V ； (2)计算化合物理论密度P。al为分子 摩尔质量M除以平均摩尔体积V ； (3)预测化合物晶体最大密度P est = 0. 232+0. 970 P。al+ 0· 502Εηομο+0· 224Elumoo应用的量子化学计算软件为Gaussian量子化学计算软件。Gaussian量子化学计算软件包是目前通常使用的商业产品，也可采用其他商业或 自编软件程序进行量子化学计算。对选取的20个常见炸药化合物晶体进行预测，其预测密度与实测密度的最大相 对误差仅为4. 14%，平均相对误差仅为2. 15%。本专利技术的目的，技术方案及效果将结合实施例进行详细说明。附图说明附图1是采用B3LYP/6-31G*方法优的HNIW可能构型的分子结构HNIW1。附图2是采用B3LYP/6-31G*方法优的HNIW可能构型的分子结构本文档来自技高网...

【技术保护点】
一种化合物晶体密度的预测方法，包括以下步骤：（１）应用量子化学计算软件，采用密度泛函Ｂ３ＬＹＰ／６－３１Ｇ↑［＊］方法优化化合物的可能构型的分子结构，确定能量最低的分子构型，得到该分子结构能量最高的电子占有轨道的能量Ｅ↓［ＨＯＭＯ］和能量最低的空轨道的能量Ｅ↓［ＬＵＭＯ］，然后采用Ｍｏｎｔｅ－Ｃａｒｌｏ方法计算该分子构型的摩尔体积，每一个分子结构进行５０－３００次，计算平均值得到摩尔体积Ｖ；（２）计算化合物理论密度ρ↓［ｃａｌ］为分子摩尔质量Ｍ除以平均摩尔体积Ｖ；（３）预测化合物晶体最大密度ρ↓［ｅｓｔ］＝０．２３２＋０．９７０ρ↓［ｃａｌ］＋０．５０２Ｅ↓［ＨＯＭＯ］＋０．２２４Ｅ↓［ＬＵＭＯ］。

【技术特征摘要】

【专利技术属性】
技术研发人员：吴玉凯，陈冠玉，周智明，董海山，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：11[中国|北京]