CN114464263B 一种分子晶体能量计算方法、装置及存储介质 （深圳晶泰科技有限公司）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN114464263B(45)授权公告日2025.07.04

(21)申请号202111600787.1

(22)申请日2021.12.24

(65)同一申请的已公布的文献号申请公布号CN114464263A

(43)申请公布日2022.05.10

(73)专利权人深圳晶泰科技有限公司

地址518000广东省深圳市福田区福保街

道福保社区红柳道2号顺丰工业厂房3

层

(72)发明人曾群李柏辰付文博

(74)专利代理机构深圳市科吉华烽知识产权事务所(普通合伙)44248

专利代理师胡吉科

(51)Int.CI.

G16C10/00(2019.01)

G16C20/30(2019.01)

(56)对比文件

CN108959842A,2018.12.07CN112086132A,2020.12.15

审查员李椰

权利要求书4页说明书15页附图3页

(54)发明名称

一种分子晶体能量计算方法、装置及存储介

质

获取目标晶体结构，并从目标晶体结构中确定中心分子、M

获取目标晶体结构，并从目标晶体结构中确定中心分子、M个团簇和M个壳结构

采用量子化学方法分别计算中心分子的分子能量、每一团簇的团簇能量和每一壳结构的壳能量

根据分子能量、团簇能量和壳能量，计算得到目标品体结构的品格能量

根据分子能量和晶格能量，计算得到目标晶体结构的晶体总能量

S110

S120

S130

S140

CN114464263B本申请提供一种分子晶体能量计算方法、装置及存储介质。其中，该方法包括：获取目标晶体结构，并从目标晶体结构中确定中心分子、M个团簇和M个壳结构；采用量子化学方法分别计算中心分子的分子能量、每一团簇的团簇能量和每一壳结构的壳能量；根据分子能量、团簇能量和壳能量，计算得到目标晶体结构的晶格能量；根据分子能量和晶格能量，计算得到目标晶体结构的晶体总能量。本申请的技术方案，通过建立“核-

CN114464263B

CN114464263B权利要求书1/4页

2

1.一种分子晶体能量计算方法，其特征在于，包括：

获取目标晶体结构，并从所述目标晶体结构中确定中心分子、M个团簇和M个壳结构，所述M为大于或等于1的整数；

采用量子化学方法分别计算所述中心分子的分子能量、每一所述团簇的团簇能量和每一所述壳结构的壳能量；

根据所述分子能量、所述团簇能量和所述壳能量，计算得到所述目标晶体结构的晶格能量；

根据所述分子能量和所述晶格能量，计算得到所述目标晶体结构的晶体总能量；

其中，所述从所述目标晶体结构中确定中心分子、M个团簇和M个壳结构，包括：

从所述目标晶体结构中选取一个分子作为中心分子；

以所述中心分子的几何中心为球心，且以M个不同的预设半径作为截断半径分别对所述目标晶体结构进行原子截取，选取原子位于所述截断半径内的分子构建团簇，得到M个所述团簇；

将所述中心分子从每一所述团簇中进行删除，得到对应的M个壳结构。

2.如权利要求1所述的方法，其特征在于，所述选取原子位于所述截断半径内的分子构建团簇，包括：

选取所有原子均位于所述截断半径内的分子构建团簇；或者，

若有分子仅部分原子位于所述截断半径内，对所述分子进行补全处理，并利用补全后的分子以及所有原子均位于所述截断半径内的分子构建团簇。

3.如权利要求1所述的方法，其特征在于，所述采用量子化学方法分别计算所述中心分子的分子能量、每一所述团簇的团簇能量和每一所述壳结构的壳能量，包括：

以所述中心分子的分子结构作为输入，采用量子化学方法计算得到所述中心分子的分子能量以及构成所述分子能量的各分项能量；

分别以每一所述团簇的结构作为输入，采用所述量子化学方法计算得到每一所述团簇的团簇能量以及构成所述团簇能量的各分项能量；

分别以每一所述壳结构作为输入，采用所述量子化学方法计算得到每一所述壳结构的壳能量以及构成所述壳能量的各分项能量。

4.如权利要求3所述的方法，其特征在于，所述量子化学方法包括基于密度泛函理论的紧束缚DFTB方法，能量计算公式为：

能=Eo近n+Esc+Erep+Eme+Edisp_cor;

其中，各分项能量E0近似Escc、Erep、Enuc、Edisp_corr分别为：零阶近似下的轨道能，二/三阶静电能，价键的近程排斥能，核排斥能，长程的色散校正能。

5.如权利要求3所述的方法，其特征在