[发明专利]一种NTO晶体晶面间结合能的计算方法

权利要求书

1.一种NTO晶体晶面间结合能的计算方法，其特征在于，包括：

对获得的α型NTO晶胞的NTO原子进行调整，得到待计算的NTO晶胞；

对所述待计算的NTO晶胞进行生长预测，得到待计算的NTO晶胞的5个主要晶面及所述主要晶面占总面积的百分比数值；

对各个所述主要晶面进行切面和扩胞，得到各个所述主要晶面的晶面模型；

对任意两个所述晶面模型进行晶格参数的近似匹配，得到作用模型；所述作用模型中的晶体取向以所述百分比数值最大的晶面为基准；

对所述作用模型的c参数进行调整，以添加用于分子动力学计算的真空层；

对调整后的作用模型依次进行几何优化和分子动力学模拟，得到分子动力学结果；

对所述分子动力学结果进行能量分析，并根据分析结果计算结合能；

所述分子动力学模拟的模拟参数具体包括：温度为298k，步数为100000，每50步输出一帧；

对所述分子动力学结果进行能量分析，并根据分析结果计算结合能，包括：

对所述分析结果进行排序，选择后50000步中能量最低步对应三帧分别进行体系总能量E_tot和两个晶面超胞能量E_l1、E_l2的计算；

根据公式E_int=E_tot-E_l1-E_l2进行计算，分别得到三帧的结合能，并将三帧对应的结合能的平均值作为最终结果；其中，E_int为所述结合能。

2.根据权利要求1所述的NTO晶体晶面间结合能的计算方法，其特征在于，对获得的α型NTO晶胞的NTO原子进行调整，得到待计算的NTO晶胞，包括：

从晶体数据中心数据库中获取所述α型NTO晶胞；

根据所述α型NTO晶胞的各原子的成键方式和化学环境调整各原子力场，得到所述待计算的NTO晶胞。

3.根据权利要求1所述的NTO晶体晶面间结合能的计算方法，其特征在于，所述对任意两个所述晶面模型进行晶格参数的近似匹配，得到作用模型，包括：

将两个所述晶面模型的晶胞的a、b、γ参数分别取平均值作为作用模型的新的晶胞参数，以使得所述作用模型容纳两个晶面超胞结构，以得到所述作用模型。

4.根据权利要求1所述的NTO晶体晶面间结合能的计算方法，其特征在于，所述真空层的厚度为20Å。

5.根据权利要求1所述的NTO晶体晶面间结合能的计算方法，其特征在于，在对所述分子动力学结果进行能量分析，并根据分析结果计算结合能之后，还包括：

根据各个作用模型的结合能计算结果进行对比，确定引发NTO聚集现象的主要晶面。

6.根据权利要求5所述的NTO晶体晶面间结合能的计算方法，其特征在于，所述引发NTO聚集现象的主要晶面为（0 0 1）。