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

道福保社区红柳道2号顺丰工业厂房3层质簇和M个壳结构；采用量子化学方法分别计算中2采用量子化学方法分别计算所述中心分子的分子能量、每一所述团簇根据所述分子能量、所述团簇能量和所述壳能量，计算得到所述目从所述目标晶体结构中选取一个分子作为中心以所述中心分子的几何中心为球心，且以M个不同的预设半径作为截断半径分别对所若有分子仅部分原子位于所述截断半径内，对所述分子进行补分别以每一所述团簇的结构作为输入,采用所述量子化学方法计算得到每一所述团簇分别以每一所述壳结构作为输入，采用所述量子化学5.如权利要求3所述的方法，其特征在于，所述3根据所述分子能量的各分项能量、每一所述团簇能量的各分项能根据所述中心分子与每一所述壳结构之间的相互作用能的各分根据所述晶格能量的各分项能量以及各分项能量的贡献系数，–E所述根据所述中心分子与每一所述壳结构之间的相互作用能的各分项能对所述中心分子与每一所述壳结构之间的相互作用能的各分项能量利用所述第二预中心分子与第2个壳结构之间的相互作用能的第i个分项能量，所述R1为构建第1个团簇时；所述中心分子与壳结构之间的相互作用能的第i4获取多个参考晶体结构以及每一所述参考晶体结分别计算每一所述参考晶体结构的晶格能量的各分项能量；设定所述晶格能量的各根据每一所述参考晶体结构的晶格能量的各分项能量以及设定的贡献利用每一所述参考晶体结构的参考晶格能量的各分项能量和预测晶格能量的各分项利用如权利要求1-13任一项所述的方法对每一根据每一所述晶体结构的晶体总能量，确定所述至少两个晶体结构第一计算模块，用于采用量子化学方法分别计算所述中第三计算模块，用于根据所述分子能量和所述晶格能量，计算得从所述目标晶体结构中选取一个分子作为中心分子；以所述中心分子的几何中心为球心，且以M个不同的预设半径作为截断半径分别对所述目标晶体结构进行原子截取，选取原子计算模块，用于利用如权利要求1-13任一项所述的确定模块，用于根据每一所述晶体结构的晶体总能量，确56[0004]由于计算资源限制，现有计算方案通常会按照一定的规[0013]以所述中心分子的几何中心为球心，且以M个不同的预设半径作为截断半径分别7全后的分子以及所有原子均位于所述截断半径内的分[0020]分别以每一所述团簇的结构作为输入,采/预设公式计算得到所述目标晶体结构的晶格能量–E述中心分子与所述第n个壳结构之间的相互作用能的第i+A*(Rn)-B；8[0039]对所述中心分子与每一所述壳结构之间的相互作用能的各分项能量利用所述第分子与第2个壳结构之间的相互作用能的第i个分项能量，所述R1为构建第1个团簇时采用；所述中心分子与壳结构之间的相互作用能的第i[0051]获取多个参考晶体结构以及每一所述参考晶体结构的参考晶格能量的各分项能[0055]利用每一所述参考晶体结构的参考晶格能量的各分项能量和预测晶格能量的各[0061]利用如本申请第一方面提供的所述分子晶体能量计算方法对每一所述晶体结构9处理器执行如本申请第一方面提供的所述分子晶体能量计算方法或者如本申请第二方面分子晶体能量计算方法或者如本申请第二方面提供的所述分[0086]在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申含一个或多个相关联的列出项目的任何或所体结构中选取相互之间最靠近(分子几何中心之间的距离最短)的分子A和分子B共同作为分子的分子能量以及构成分子能量的各分项能量；分别以每一团簇的结构作为输入,采用[0101]在一实施方式中，量子化学方法可以包括基于密度泛函理论的紧束缚(Density[0104]在一实施方式中，量子化学方法可以包括密度泛函理论(DensityFunctional–E-3)n个–E2个壳结构(壳2[0119]考虑到不同截断半径(Rn)下的壳结构与中心分子之间的相互作用能可能存在衰EEE[0128]对于公式(1-4)，通常可以将B看作是一个物理确定的参数，其理论值可以设为–E-6)-7)以及这些参考晶体结构的参考晶格能量Eref,j，j为1～N。参考晶体结构是指利用理论方法理论方法(如DFT)预测得到的晶体晶格能量或者试验。[0157]5)利用每一参考晶体结构的参考晶格能量的各分项能量和预测晶格能量的各分E[0164]1)选择4个等价分子中的一个作为中心分子，计算E分子,0近似、E分子,scc、E分子,rep、Aspirin晶体中的中心分子的分子结构如图2所示。(1-5)进行外包含4个等价分子。本实例选择密度泛函理论(DFT)方法中正常精度的PBE-D3BJ方法和6-[0180]4)使用PBE-D3BJ方法和def-SV(P)对中心分子、团簇1和壳1结构进行各分项能量3所示：表3为PBE-D3BJ和def-SV(P)下得到的各个分量和[0184]5)选择预先针对PBE-D3BJ方法和def-SV(P)拟合的KPBE-D3BJ/6-31G*,壳1,使用公式(1--117.1kJ/mol,也非常接近文献报道的高精度计算的结果。度的方法(如MP2/aug-cc-pVQZ方法)对中心分子的分子能量进行重新计算，以替换掉前面[0190]本申请使用中心分子-外围团簇进行晶体能量评估的技术，兼容主流的量子化学[0206]可选的，获取模块610选取原子位于截断半径内的分子构建团簇的实施方式可以[0209]第二计算单元，用于分别以每一团簇/[0217]第五计算单元，用于根据中心分子与每一壳结构之间的相互作用能的各分项能En个团簇的+A*(Rn)-B；；量；利用每一参考晶体结构的参考晶格能量的各分项能量和预测晶格能量的各分项能量，通用处理器、数字信号处理器(DigitalSignalProcessor，DSP)、专用集成电路(ApplicationSpecificIntegratedCircuit，ASIC)、现场可编程门阵列(Field-分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器存储媒介不包含载波和通过无线或有线传输的瞬算机程序或计算机程序产品包括用于执行本申请的上述方法中部分或全部步骤的计算机[0247]或者，本申请还可以实施为一种计算机可读存储介质(或