CN108268750B - 基于枚举Wyckoff位置组合的假想无机晶体结构预测方法 - Google Patents

Abstract

本发明提供了一种基于枚举Wyckoff位置组合的假想无机晶体结构预测方法，第一步，选择进行结构预测的目标空间群，确定晶胞参数和骨架密度的范围；第二步，枚举出Wyckoff位置组合；第三步，进行结构模型的搭建；本发明在特定的空间群下，通过枚举出的Wyckoff位置组合，将原子固定在Wyckoff位置上进行假想无机晶体结构模型搭建，最终完成假想无机晶体结构预测，此方法能够均衡地搜索整个构象空间，提高结构模型的搭建效率，预测出更多合理的无机晶体结构。

Description

基于枚举Wyckoff位置组合的假想无机晶体结构预测方法

技术领域

本发明属于无机晶体结构预测领域，通过枚举Wyckoff位置组合的方法进行假想无机晶体的结构预测。

背景技术

近年来，由于在无机晶体中发现了许多独特的现象，并在技术上得到了广泛的应用，无机晶体领域得到了迅猛的发展。如MOFs和沸石分子筛等结构，拥有着极大的比表面积和大而均一的孔径，在分离科学和催化领域发挥着关键的作用。将无机晶体的结构、生长和性质作为一个统一的综合性问题，不仅研究晶体的理想结构，而且研究带有各种缺陷的实际结构，从而指导发现具有良好性能的无机晶体结构。然而，无机晶体结构的开发仍面临着极大地挑战，预测具有新颖拓扑类型的无机晶体结构一直是该领域迫切的需求。国际上主流的结构预测方法在搭建模型的过程中，骨架原子随机进行移动，以原子融合的方式进入Wyckoff位置，以实现结构模型的搭建。这种方法的劣势在于：原子无法均衡地搜索整个构象空间，忽略了大量可能存在的结构模型，极大地降低了结构的预测效率。本发明枚举出所有符合总骨架原子数目范围的Wyckoff位置组合，将每个原子安排在空间群中固定的Wyckoff位置上，进行假想无机晶体结构模型的搭建。这种方法的优势在于：所有骨架原子被安排在预先设定的Wyckoff位置上，实现对整个构象空间均衡搜索，从而得到数量较多的假想无机晶体结构，提高结构预测效率。

发明内容

本发明目的在于提供一种均衡高效的无机晶体结构预测方法，采用枚举Wyckoff位置组合的方式进行假想无机晶体结构模型的搭建，从而对无机晶体结构进行预测，本发明是采用如下具体技术方案实现的：

基于枚举Wyckoff位置组合的假想无机晶体结构预测方法，其特征在于，根据总骨架原子数目范围枚举出Wyckoff位置组合，将每个原子安排在空间群中固定的Wyckoff位置上进行假想无机晶体结构的预测，具体步骤如下：

第一步，选择进行假想无机晶体结构预测的目标空间群，确定晶胞参数和骨架密度FD的范围：

在选定目标空间群后，晶胞参数中α、β、γ、A、B和C的晶胞参数特征随之被确定，晶胞参数中A、B和C的初始值A₀、B₀和C₀可任意设定，随后选取增长步幅m_A、m_B和m_C沿晶轴方向增加A、B和C的值，从而逐步扩大晶胞，如式(1)所示：

A＝A₀+m_A×n_A；B＝B₀+m_B×n_B；C＝C₀+m_C×n_C； (1)

其中，n_A、n_B和n_C为非负整数，m_A、m_B和m_C的值可任意设定，可得到[(n_A+1)×(n_B+1)×(n_C+1)]种晶胞参数组合方式；删去不符合目标空间群下晶胞参数特征的晶胞参数组合方式；每种组合方式下的晶胞形状、大小均已被α、β、γ、A、B和C所确定，可据此计算出晶胞的体积V；骨架密度FD的范围可任意设定，骨架密度FD的最小值用FD_min表示，骨架密度FD的最大值用FD_max表示，根据式(2)，可求出每种晶胞参数组合方式下向目标晶胞中投放的总骨架原子数目T的最小值T_min和最大值T_max，即总骨架原子数目T的范围：

T_min＝V×FD_min；T_max＝V×FD_max (2)

第二步，枚举出Wyckoff位置组合：

在确定每种晶胞参数组合方式下需投放的总骨架原子数目T的范围后，定义P为选定的目标空间群下存在的Wyckoff位置种类数，如式(3)所示，可枚举出每种晶胞参数组合方式下的Wyckoff位置组合方式：

N_x为位于第x个Wyckoff位置的独立骨架原子数目，g_x为相应Wyckoff位置的多重度，T_min、T_max和g_x均为已知，可得到若干组N_x的组合，即用来进行假想无机晶体结构预测的Wyckoff位置组合。

第三步，进行结构模型的搭建：

确定晶胞的形状、大小和向其中投放的总骨架原子数目范围后，将所有骨架原子按照枚举出的Wyckoff位置组合方式固定在Wyckoff位置上，定义势函数E为搭建出的假想无机晶体结构与理想无机晶体结构在拓扑和几何上的差异，如式(4)所示：

E＝w₁E₁+w₂E₂+w₃E₃ (4)

w₁、w₂和w₃分别为E₁、E₂和E₃的权重，其数值由操作者根据E₁、E₂和E₃各自的重要程度进行指定，E₁为假想无机晶体结构与理想无机晶体结构在配位数上的偏差，E₂为假想无机晶体结构与理想无机晶体结构在键长上的偏差，E₃为假想无机晶体结构与理想无机晶体结构在键角上的偏差，由E的定义可知，E的数值恒大于零，且E值越小，代表搭建出的假想无机晶体结构在原子排布上越接近理想无机晶体结构，随机修改骨架原子的坐标，使所有骨架原子在整个晶胞内移动，骨架原子的移动导致E₁、E₂和E₃的改变，实时计算每次骨架原子移动后的假想无机晶体结构的E值，控制骨架原子坐标向降低E值的方向移动，并以操作者指定的概率允许骨架原子向E值升高的方向移动，骨架原子经过操作者指定的移动次数进行移动，当E值不再变化后，存储所有骨架原子的原子坐标和对应的E值，完成该Wyckoff位置组合的结构预测，重复此步骤，直至完成所有枚举出的Wyckoff位置组合计算，实现目标空间群下对假想无机晶体结构的预测。

附图说明

下面结合附图和实例对本发明进一步详细说明：

图1为本发明所述的基于枚举Wyckoff位置组合的假想无机晶体结构预测方法流程图；

图2为本方法预测出的已知分子筛结构JRY的结构模型图；

图3为本方法预测出的假想分子筛结构H024ddd的结构模型图：

具体实施方式

下面结合附图和实例对本发明进一步详细说明：

本实例是在I2₁2₁2₁(NO.24)空间群下，通过枚举Wyckoff位置组合的方法进行假想无机分子筛结构的预测：

第一步，确定晶胞参数和骨架密度FD的范围。I2₁2₁2₁(NO.24)空间群属于正交晶系，晶胞参数特征为α＝β＝γ＝90°，设定晶胞参数中A、B和C的初始值A₀、B₀、C₀为

将增长步幅定为

沿晶轴方向将A、B和C的值逐渐增长到

如式(1)所示，可得到125种不同的晶胞参数组合方式。从中选取5组晶胞参数组合方式作为示例进行下一步计算，每种晶胞参数组合方式下的晶胞形状、大小已经确定，可据此计算出晶胞的体积V＝A×B×C。将骨架密度范围定为

根据式(2)，可求出每种晶胞参数组合方式下的总骨架原子数目T的范围，如表1所示。

第二步，枚举出Wyckoff位置组合。由I2₁2₁2₁(NO.24)空间群下存在的4个Wyckoff位置(a、b、c、d)所对应的多重度(g₁＝4、g₂＝4、g₃＝4、g₄＝8)，依据表1中每组晶胞参数组合下的总骨架原子数目范围，枚举出Wyckoff位置组合方式，如式(3)所示，T_min、T_max、和多重度数g_x的值均为已知，从而可得到若干组符合式(3)的N₁、N₂、N₃和N₄的组合，据此可枚举出Wyckoff位置组合，如表2所示。

第三步，进行结构模型的搭建：

确定晶胞的大小和向其中投放的总骨架原子数目范围后，将所有骨架原子按照枚举出的Wyckoff位置组合方式固定在Wyckoff位置上，定义势函数E为搭建出的假想无机分子筛结构与理想无机分子筛结构在拓扑和几何上的差异，如式(4)所示，设E₁的权重w₁＝5、E₂的权重w₂＝20、E₃的权重w₃＝5，理想分子筛结构的配位数、相邻骨架原子间距、次相邻骨架原子间距分别为4、

据此计算E1、E2、E3，由E的定义可知，E的数值恒大于零，且E值越小，代表搭建出的假想分子筛结构越合理，随机修改骨架原子的坐标，使所有骨架原子在整个晶胞内移动，实时计算每次骨架原子移动后的假想无机分子筛结构的E值，以一组Wyckoff位置为例，表3记录了晶胞参数为

Wyckoff位置组合为bd时的E值变化，在该实例下，操作者设定原子移动步数为10000，以50％的概率接受骨架原子的移动，每2000步记录一次E值，当E值不再变化后，所有骨架原子不再移动且坐标被存储，完成该Wyckoff位置组合下的假想分子筛结构预测，重复此步骤，完成所有枚举出的Wyckoff位置组合计算，用此方法，共得到了114种Wyckoff位置组合，按枚举出的Wyckoff位置组合固定原子进行结构搭建，得到了225种假想无机晶体结构，其中当晶胞参数设定为

时，成功搭建已知分子筛结构JRY(附图1)，同时得到了对称性较高的假想分子筛结构H024ddd(附图2)

A＝10+3×n_A；B＝10+3×n_B；C＝10+3×n_C (1)

其中，n_A、n_B、n_C为非负整数

T_min＝V×FD_min；T_max＝V×FD_max (2)

T_min≤N₁×4+N₂×4+N₃×4+N₄×8≤Tmax (3)

其中，N₁、N₂、N₃和N₄为位于Wyckoff位置a、b、c和d的独立骨架原子数目；

E＝5 E₁ +20E₂+ 5E₃ (4)

表1

表2

表3

Wyckoff位置组合	骨架原子移动次数	E值
			bd	2000	7.200
bd	4000	7.193
			bd	6000	7.193
bd	8000	7.193
			bd	10000	7.193

Claims (1)

1.基于枚举Wyckoff位置组合的假想无机晶体结构预测方法，其特征在于，根据总骨架原子数目范围枚举出Wyckoff位置组合，将每个原子安排在空间群中固定的Wyckoff位置上进行假想无机晶体结构的预测，具体步骤如下：

第一步，选择进行假想无机晶体结构预测的目标空间群，确定晶胞参数和骨架密度FD的范围：

在选定目标空间群后，晶胞参数中α、β、γ、A、B和C的晶胞参数特征随之被确定，晶胞参数中A、B和C的初始值A₀、B₀和C₀可任意设定，随后选取增长步幅m_A、m_B和m_C沿晶轴方向增加A、B和C的值，从而逐步扩大晶胞，如式(1)所示：

A＝A₀+m_A×n_A；B＝B₀+m_B×n_B；C＝C₀+m_C×n_C； (1)

其中，n_A、n_B和n_C为非负整数，m_A、m_B和m_C的值可任意设定，可得到[(n_A+1)×(n_B+1)×(n_C+1)]种晶胞参数组合方式；删去不符合目标空间群下晶胞参数特征的晶胞参数组合方式；每种组合方式下的晶胞形状、大小均已被α、β、γ、A、B和C所确定，可据此计算出晶胞的体积V；骨架密度FD的范围可任意设定，骨架密度FD的最小值用FD_min表示，骨架密度FD的最大值用FD_max表示，根据式(2)，可求出每种晶胞参数组合方式下向目标晶胞中投放的总骨架原子数目T的最小值T_min和最大值T_max，即总骨架原子数目T的范围：

T_min＝V×FD_min；T_max＝V×FD_max (2)

第二步，枚举出Wyckoff位置组合：

在确定每种晶胞参数组合方式下需投放的总骨架原子数目T的范围后，定义P为选定的目标空间群下存在的Wyckoff位置种类数，如式(3)所示，可枚举出每种晶胞参数组合方式下的Wyckoff位置组合方式：

N_x为位于第x个Wyckoff位置的独立骨架原子数目，g_x为相应Wyckoff位置的多重度，T_min、T_max和g_x均为已知，可得到若干组N_x的组合，即用来进行假想无机晶体结构预测的Wyckoff位置组合；

第三步，进行结构模型的搭建：

确定晶胞的形状、大小和向其中投放的总骨架原子数目范围后，将所有骨架原子按照枚举出的Wyckoff位置组合方式固定在Wyckoff位置上，定义势函数E为搭建出的假想无机晶体结构与理想无机晶体结构在拓扑和几何上的差异，如式(4)所示：

E＝w₁E₁+w₂E₂+w₃E₃ (4)

w₁、w₂和w₃分别为E₁、E₂和E₃的权重，其数值由操作者根据E₁、E₂和E₃各自的重要程度进行指定，E₁为假想无机晶体结构与理想无机晶体结构在配位数上的偏差，E₂为假想无机晶体结构与理想无机晶体结构在键长上的偏差，E₃为假想无机晶体结构与理想无机晶体结构在键角上的偏差，由E的定义可知，E的数值恒大于零，且E值越小，代表搭建出的假想无机晶体结构在原子排布上越接近理想无机晶体结构，随机修改骨架原子的坐标，使所有骨架原子在整个晶胞内移动，骨架原子的移动导致E₁、E₂和E₃的改变，实时计算每次骨架原子移动后的假想无机晶体结构的E值，控制骨架原子坐标向降低E值的方向移动，并以操作者指定的概率允许骨架原子向E值升高的方向移动，骨架原子经过操作者指定的移动次数进行移动，当E值不再变化后，存储所有骨架原子的原子坐标和对应的E值，完成该Wyckoff位置组合的结构预测，重复此步骤，直至完成所有枚举出的Wyckoff位置组合计算，实现目标空间群下对假想无机晶体结构的预测。

Images