[发明专利]基于分子间相互作用的粗粒化方法

权利要求书

1.基于分子间相互作用的粗粒化方法，其特征在于：包括以下步骤：

步骤1、对晶体材料进行量子力学计算：根据晶体材料的晶向对晶体材料的单体分子进行建模，并用QM方法进行计算，计算完成后，改变两个单体分子之间的距离，重新计算；经过一系列的计算后，得到距离-能量矩阵；

步骤2、对三斜晶系进行校正：对于三斜晶系材料，对其晶体结构进行正交化处理，得到斜方晶系；

步骤3、建立伪化学键模型：根据步骤2的斜方晶系建立伪化学键模型，具体为，模型x、y、z三个方向的距离依据校正后的斜方晶系晶格常数确定，然后将每一个晶体材料单体粗粒化为一个珠子，珠子与珠子之间，使用伪化学键连接，伪化学键不是真正意义上的化学键，他的长度远远超出了化学键的长度，表示成一种分子间的相互作用；

步骤4、对Morse函数参数进行拟合：对步骤1计算得到的距离-能量矩阵进行Morse函数参数拟合，得到D、α、r₀的参数矩阵；

步骤5、将步骤4获得的参数矩阵和步骤3建立的模型带入MD框架进行模拟；

步骤6、输出模拟结果：输出需要的计算结果或者可视化轨迹。

2.根据权利要求1所述基于分子间相互作用的粗粒化方法，其特征在于：步骤2所述对于三斜晶系材料的正交化处理为，a′＝a,b′＝b×cos(γ-90),c′＝c×cos(α-90),a、b、c、α、β、γ是指原三斜晶系中的6个晶体参数，其中a、b、c是指晶胞的三组棱长，α、β、γ是指三组棱相互间的夹角，带′的是正交化为斜方晶系后的6个晶体参数。

3.根据权利要求1所述基于分子间相互作用的粗粒化方法，其特征在于：所述伪化学键的长度远远超出化学键的长度，且不表示成对电子，而表示为一种分子间的相互作用，且不同方向的伪化学键的键长和类型不同。

4.根据权利要求1-3任一项所述基于分子间相互作用的粗粒化方法，其特征在于：步骤5MD框架模拟中将键角设置为90°或180°。

5.根据权利要求4所述基于分子间相互作用的粗粒化方法，其特征在于：在MD模拟中，在一个晶体颗粒内，使用伪化学键描述分子之间的相互作用，而在不同的晶体颗粒之间使用Morse对势来描述相互作用。