(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111673770.9 (22)申请日 2021.12.31 (71)申请人 深圳晶泰科技有限公司 地址 518000 广东省深圳市福田区福保街 道福保社区红柳道2号顺丰工业厂房3 层 (72)发明人 方利文　孙广旭　 (74)专利代理 机构 深圳市科吉华烽知识产权事 务所(普通 合伙) 44248 专利代理师 胡吉科 (51)Int.Cl. G16C 10/00(2019.01) G06F 30/20(2020.01) G06F 119/14(2020.01) (54)发明名称 基于分子动力学的晶体熔点的计算方法、 装 置及存储介质 (57)摘要 本申请提供了一种基于分子动力学的晶体 熔点的计算方法、 装置及存储介质。 其中， 该方法 包括： 获取目标晶体结构， 并根据目标晶体结构 构建周期性结构模型； 利用周 期性结构模型， 生 成含有空洞的多个目标结构模型， 其中， 每个目 标结构模型包含的空洞的体积不同； 对每个目标 结构模型进行分子动力学升温模拟， 分别获得每 个目标结构 模型在预设温度范围下的结构信息； 根据每个目标结构模型的结构信息在预设温度 范围内的变化状态， 得到目标晶体结构的熔点 值。 本申请的技术方案， 能够提高熔点的计算精 度。 权利要求书3页 说明书13页 附图8页 CN 114464264 A 2022.05.10 CN 114464264 A 1.一种基于分子动力学的 晶体熔点的计算方法， 其特 征在于， 包括： 获取目标晶体结构， 并根据所述目标晶体结构 构建周期性结构模型； 利用所述周期性结构模型， 生成含有空洞的多个目标结构模型， 其中， 每个所述目标结 构模型包 含的所述空洞的体积不同； 对每个所述目标结构模型进行分子动力学升温模拟， 分别获得每个所述目标结构模型 在预设温度范围下的结构信息； 根据每个所述目标结构模型的结构信 息在所述预设温度 范围内的变化状态， 得到所述 目标晶体结构的熔点 值。 2.根据权利要求1所述的方法， 其特 征在于， 所述方法还 包括： 对所述周期性结构模型进行扩胞， 得到超胞周期性结构模型； 所述利用所述周期性结构模型， 生成含有空洞的多个目标 结构模型， 包括： 利用所述超胞周期性结构模型， 生成含有空洞的多个目标 结构模型。 3.根据权利要求2所述的方法， 其特征在于， 所述利用所述超胞周期性结构模型， 生成 含有空洞的多个目标 结构模型， 包括： 分别将所述超胞周期性结构模型中的不同预设阈值的原子或分子进行删除， 生成含有 空洞的多个目标 结构模型； 其中， 不同的所述预设阈值按照大小进行排序后， 任意相邻的两个所述预设阈值之间 符合预设变化 规律。 4.根据权利要求1所述的方法， 其特 征在于， 所述方法还 包括： 对每个所述目标 结构模型进行 结构优化， 分别得到优化后的所述目标 结构模型； 所述对每 个所述目标 结构模型进行分子动力学升温模拟， 包括： 对每个优化后的所述目标 结构模型进行分子动力学升温模拟。 5.根据权利要求 4所述的方法， 其特 征在于， 所述方法还 包括： 根据所述周期性结构模型， 确定目标力场； 所述对每个所述目标结构模型进行结构优化， 分别得到优化后的所述目标结构模型， 包括： 利用与所述目标力场相对应的预设算法对每个所述目标结构模型进行结构优化， 使得 所述目标 结构模型的能量 最小化， 分别得到优化后的所述目标 结构模型。 6.根据权利要求1所述的方法， 其特征在于， 所述结构信 息包括所述目标结构模型的势 能、 密度、 原子坐标和原子坐标的均方位移中的至少一种。 7.根据权利要求1 ‑6任一所述的方法， 其特征在于， 所述根据每个所述目标结构模型的 结构信息在所述预设温度范围内的变化状态， 得到所述目标晶体结构的熔点 值， 包括： 根据每个所述目标结构模型的结构信 息在所述预设温度 范围内的变化状态， 确定每个 所述目标 结构模型的结构信息在所述预设温度范围内的突变点； 根据每个所述目标结构模型的突变点对应的温度值， 确定每个所述目标结构模型的熔 点值； 根据每个所述目标 结构模型的熔点 值， 确定所述目标晶体结构的熔点 值。 8.根据权利要求7所述的方法， 其特征在于， 所述根据每个所述目标结构模型的结构信 息在所述预设温度范围内的变化状态， 确定每个所述目标结构模型的结构信息在所述预设权　利　要　求　书 1/3 页 2 CN 114464264 A 2温度范围内的突变点， 包括： 根据每个所述目标结构模型的结构信 息在所述预设温度 范围内的变化状态， 确定在所 述预设温度范围内每个所述目标结构模型的结构信息的变化量超过预设变化量的突变温 度区间； 根据每个所述目标 结构模型的突变温度区间， 确定对应的突变点。 9.根据权利要求7所述的方法， 其特征在于， 所述根据每个所述目标结构模型的熔点 值， 确定所述目标晶体结构的熔点 值， 包括： 对多个所述目标 结构模型按照空洞的体积大小 进行排序； 根据排序 结果， 确定是否存在连续排列的预设数量或预设比例的所述目标结构模型的 熔点值满足预设条件； 若存在， 根据所述连续排列的预设数量或预设比例的所述目标结构模型的熔点值， 确 定所述目标晶体结构的熔点 值。 10.根据权利要求9所述的方法， 其特征在于， 所述根据排序 结果， 确定是否存在连续排 列的预设数量或预设比例的所述目标 结构模型的熔点 值满足预设条件， 包括： 根据排序 结果， 确定是否存在连续排列的预设数量或预设比例的所述目标结构模型中 任意两个所述目标 结构模型的熔点 值之间的差异值小于第一预设值； 或者， 根据排序 结果， 确定是否存在连续排列的预设数量或预设比例的所述目标结构模型的 熔点值之间的均方误差或均方根 误差小于第二预设值。 11.根据权利要求9所述的方法， 其特征在于， 所述根据所述连续排列的预设数量或预 设比例的所述目标 结构模型的熔点 值， 确定所述目标晶体结构的熔点 值， 包括： 计算所述连续排列的预设数量或预设比例的所述目标结构模型的熔点值的平均熔点 值； 将所述平均熔点 值确定为所述目标晶体结构的熔点 值。 12.根据权利要求9所述的方法， 其特 征在于， 所述方法还 包括： 若不存在连续排列的预设数量或预设比例的所述目标结构模型的熔点值满足所述预 设条件， 利用所述周期性结构模型， 生成更多数量的含有空洞的目标结构模型； 或者利用所 述周期性结构模型， 生成空洞的体积差异更小的多个目标结构模型； 或者利用所述周期性 结构模型， 生成空洞的体积差异更小的且数量更多的目标结构模型； 或者对所述周期性结 构模型进行扩胞， 生成更 大尺寸的超胞周期性结构模型。 13.一种晶体筛 选方法， 其特 征在于， 包括： 获取至少两个晶体结构； 利用如权利要求1 ‑12任一项所述的方法对所述至少两个晶体结构进行晶体熔点计算， 得到每一所述晶体结构的熔点 值； 根据每一所述晶体结构的熔点 值， 从所述至少两个晶体结构中确定出候选晶体结构。 14.一种基于分子动力学的 晶体熔点的计算装置， 其特 征在于， 包括： 获取模块， 用于获取目标晶体结构； 构建模块， 用于根据所述目标晶体结构 构建周期性结构模型； 生成模块， 用于利用所述周期性结构模型， 生成含有空洞的多个目标结构模型， 其中， 每个所述目标 结构模型包 含的所述空洞的体积不同；权　利　要　求　书 2/3 页 3 CN 114464264 A 3