surface-construct 0.12__tar.gz → 0.12.1__tar.gz

{surface_construct-0.12/surface_construct.egg-info → surface_construct-0.12.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: surface_construct
-Version: 0.12
+Version: 0.12.1
 Summary: Surface construction and surface reaction sampling tools.
 Home-page: https://gitee.com/pjren/surface_construct/
 Author: ren
@@ -32,7 +32,9 @@ Dynamic: summary
 
 催化剂构效关系是建立表面位点结构和反应过程势能面的关联关系，这依赖对表面上各种可能位点的充分分析和采样，找到这些位点对关键反应步骤的影响情况。从催化剂的角度看，研究的目标是从表面上存在的位点中找到表面最适合反应的位点，即使反应势能面最低的位点；从反应的角度看，研究目标是对于一系列反应通道，找到反应势能面上的卡点，找到能够突破卡点的催化环境。
 
-本软件同时进行表面位点的构建、采样、计算、分析，和反应过程的分析采样，即过渡态和分子构象，由此可以实现分子在表面反应的全过程采样。其中的技术关键是采样的效率，我们针对性地使用综合使用了多样化的采样策略，保证表面各类型位点和关键位点的充分采样。
+本软件同时进行表面位点的构建、采样、计算、分析，和反应过程的分析采样，即过渡态和分子构象，由此可以实现分子在表面反应的全过程采样。不仅包括平板模型的表面，还包括分子筛孔道內、团簇表面等。
+
+表面反应过程采样的技术关键是采样的效率，我们针对性地使用综合使用了多样化的采样策略，保证表面各类型位点和关键位点的充分采样。
 
 ### 应用场景
 * 催化剂表面最可能反应位点采样
@@ -41,15 +43,17 @@ Dynamic: summary
   * 优化 C、R、Z （z方向，距离表面的距离）空间，固定COM
   * 返回所采结构中最低能量和对应的结构和能量
   
-* 催化剂表面全势能面采样
+* 催化剂表面全势能面采样（位点、分子构象、旋转）
   * 适用于催化剂位点全局考察
   * 适用于表面分子全局的构象和姿态
   * 适用于反应全过程采样
   * 不优化，或者仅优化Z空间
   * 返回全部能量和结构，以及全局的fitting结果
   * 使用能量截断（dEmax<2eV）
+* 原子力显微镜（AFM）谱图模拟
+* 分子筛孔道內分子构象采样
+
 
-<img src="docs/birds_view.png" alt="表面位点全局分析" style="zoom:50%;" />
 
 ## 程序流程图  Program Workflow
 
@@ -96,24 +100,6 @@ Dynamic: summary
 * MDMin
 * FIRE
 
-各种优化算法的对比，参考[链接](https://wiki.fysik.dtu.dk/gpaw/devel/ase_optimize/ase_optimize.html)
-
-**使用方法**
-
-修改 `surface_reaction_sample.py` 其中的一行 
-
-```
-from ase.optimize import BFGS
-```
-
-改为
-
-```
-from ase.optimize import XXX as BFGS
-```
-
-注意：目前这只是权宜之计，后面会把相应的设置加入到 `parameter.py`
-
 ### Gaussian Process Regression 方法
 
 高斯过程回归 GPR 的优点：
@@ -166,6 +152,8 @@ $$
   * v 0.9.2: debug 支持三斜晶系的表面格点化
   * v 0.9.3: debug GridGenerator
 * v 0.10: 表面终结生成完整并入
+* v0.11：分子构象和旋转采样
+* v0.12：AFM 模拟
 
 **TODO**
 
@@ -174,13 +162,16 @@ $$
 * Low-Energy Region Explorer 方法借鉴[^LoreX]: 划分能量区域， 进行 Delaunay 划分，合并较小的区域为较大的区域。然后针对低能的区域进行额外采样。也可以适用不同的方法，低精度进行低能高能区域识别，然后高精度方法加强。
 * 表面格点使用对称性进行降维。在最初期就引入对称性。参考 adsorbate rotation 思路
 * 表面位点数据库
-* 多原子体系（内坐标受限体系）
+* 覆盖度问题
+  * 表面上每增加一个吸附分子，就更新一次格点（可以仅仅更新 vip points，或者dock points）
+  * 表面上新增 dock points 属性，根据吸附能、位点类型或者用户进行定义。默认是所有的位点，根据能量分析或者用户更新
+  * sg 增加 del_grid(indices), 删除对应的 points， vectors，Dga 等等
 * 完善用户界面、例子、教程
 
 
 ## Reference
 
-[^Duvenaud]: [The Kernel Cookbook: Advice on Covariance functions](https://www.cs.toronto.edu/~duvenaud/)]
+[^Duvenaud]: [The Kernel Cookbook: Advice on Covariance functions](https://www.cs.toronto.edu/~duvenaud/)
 [^BASC]: Shane Carr, Roman Garnett, Cynthia LoBASC: Applying Bayesian Optimization to the Search for Global Minima on Potential Energy Surfaces.
 
 [^向量空间转化]: 计算实空间和向量空间的相邻格点距离的映射系数，根据此系数将实空间的距离转化为向量空间距离。

{surface_construct-0.12 → surface_construct-0.12.1}/README.md

@@ -2,7 +2,9 @@
 
 催化剂构效关系是建立表面位点结构和反应过程势能面的关联关系，这依赖对表面上各种可能位点的充分分析和采样，找到这些位点对关键反应步骤的影响情况。从催化剂的角度看，研究的目标是从表面上存在的位点中找到表面最适合反应的位点，即使反应势能面最低的位点；从反应的角度看，研究目标是对于一系列反应通道，找到反应势能面上的卡点，找到能够突破卡点的催化环境。
 
-本软件同时进行表面位点的构建、采样、计算、分析，和反应过程的分析采样，即过渡态和分子构象，由此可以实现分子在表面反应的全过程采样。其中的技术关键是采样的效率，我们针对性地使用综合使用了多样化的采样策略，保证表面各类型位点和关键位点的充分采样。
+本软件同时进行表面位点的构建、采样、计算、分析，和反应过程的分析采样，即过渡态和分子构象，由此可以实现分子在表面反应的全过程采样。不仅包括平板模型的表面，还包括分子筛孔道內、团簇表面等。
+
+表面反应过程采样的技术关键是采样的效率，我们针对性地使用综合使用了多样化的采样策略，保证表面各类型位点和关键位点的充分采样。
 
 ### 应用场景
 * 催化剂表面最可能反应位点采样
@@ -11,15 +13,17 @@
   * 优化 C、R、Z （z方向，距离表面的距离）空间，固定COM
   * 返回所采结构中最低能量和对应的结构和能量
   
-* 催化剂表面全势能面采样
+* 催化剂表面全势能面采样（位点、分子构象、旋转）
   * 适用于催化剂位点全局考察
   * 适用于表面分子全局的构象和姿态
   * 适用于反应全过程采样
   * 不优化，或者仅优化Z空间
   * 返回全部能量和结构，以及全局的fitting结果
   * 使用能量截断（dEmax<2eV）
+* 原子力显微镜（AFM）谱图模拟
+* 分子筛孔道內分子构象采样
+
 
-<img src="docs/birds_view.png" alt="表面位点全局分析" style="zoom:50%;" />
 
 ## 程序流程图  Program Workflow
 
@@ -66,24 +70,6 @@
 * MDMin
 * FIRE
 
-各种优化算法的对比，参考[链接](https://wiki.fysik.dtu.dk/gpaw/devel/ase_optimize/ase_optimize.html)
-
-**使用方法**
-
-修改 `surface_reaction_sample.py` 其中的一行 
-
-```
-from ase.optimize import BFGS
-```
-
-改为
-
-```
-from ase.optimize import XXX as BFGS
-```
-
-注意：目前这只是权宜之计，后面会把相应的设置加入到 `parameter.py`
-
 ### Gaussian Process Regression 方法
 
 高斯过程回归 GPR 的优点：
@@ -136,6 +122,8 @@ $$
   * v 0.9.2: debug 支持三斜晶系的表面格点化
   * v 0.9.3: debug GridGenerator
 * v 0.10: 表面终结生成完整并入
+* v0.11：分子构象和旋转采样
+* v0.12：AFM 模拟
 
 **TODO**
 
@@ -144,13 +132,16 @@ $$
 * Low-Energy Region Explorer 方法借鉴[^LoreX]: 划分能量区域， 进行 Delaunay 划分，合并较小的区域为较大的区域。然后针对低能的区域进行额外采样。也可以适用不同的方法，低精度进行低能高能区域识别，然后高精度方法加强。
 * 表面格点使用对称性进行降维。在最初期就引入对称性。参考 adsorbate rotation 思路
 * 表面位点数据库
-* 多原子体系（内坐标受限体系）
+* 覆盖度问题
+  * 表面上每增加一个吸附分子，就更新一次格点（可以仅仅更新 vip points，或者dock points）
+  * 表面上新增 dock points 属性，根据吸附能、位点类型或者用户进行定义。默认是所有的位点，根据能量分析或者用户更新
+  * sg 增加 del_grid(indices), 删除对应的 points， vectors，Dga 等等
 * 完善用户界面、例子、教程
 
 
 ## Reference
 
-[^Duvenaud]: [The Kernel Cookbook: Advice on Covariance functions](https://www.cs.toronto.edu/~duvenaud/)]
+[^Duvenaud]: [The Kernel Cookbook: Advice on Covariance functions](https://www.cs.toronto.edu/~duvenaud/)
 [^BASC]: Shane Carr, Roman Garnett, Cynthia LoBASC: Applying Bayesian Optimization to the Search for Global Minima on Potential Energy Surfaces.
 
 [^向量空间转化]: 计算实空间和向量空间的相邻格点距离的映射系数，根据此系数将实空间的距离转化为向量空间距离。

{surface_construct-0.12 → surface_construct-0.12.1}/setup.py

@@ -15,7 +15,7 @@ install_requires = [
 
 setup(
     name='surface_construct',
-    version='0.12',
+    version='0.12.1',
     packages=find_packages(),
     url='https://gitee.com/pjren/surface_construct/',
     license='GPL',

{surface_construct-0.12 → surface_construct-0.12.1}/surface_construct/structures/surface_grid.py

@@ -98,11 +98,12 @@ class MLATVectorGen:
             index_array = grid_dist.argsort(axis=-1)
             grid_dist = np.take_along_axis(grid_dist, index_array, axis=-1)
             # 设定距离向量长度
-            vlen = self.vlen
             r0 = self.kwargs.get('r0_dict')[atomnum]  # r0_dict 是一个字典，放着 r0 值
-            v = grid_dist[:, :vlen]
             kwargs.update({'r0': r0})
-            v_w = self.weight_func(v, **kwargs)
+            vlen = min(self.vlen,grid_dist.shape[-1])
+            rt = grid_dist[:, :vlen]
+            v_w = np.zeros((grid_dist.shape[0],self.vlen))
+            v_w[:,:vlen] = self.weight_func(rt, **kwargs)
             vector.append(v_w)
         vector = np.concatenate(vector, axis=1)
         return vector

{surface_construct-0.12 → surface_construct-0.12.1/surface_construct.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: surface_construct
-Version: 0.12
+Version: 0.12.1
 Summary: Surface construction and surface reaction sampling tools.
 Home-page: https://gitee.com/pjren/surface_construct/
 Author: ren
@@ -32,7 +32,9 @@ Dynamic: summary
 
 催化剂构效关系是建立表面位点结构和反应过程势能面的关联关系，这依赖对表面上各种可能位点的充分分析和采样，找到这些位点对关键反应步骤的影响情况。从催化剂的角度看，研究的目标是从表面上存在的位点中找到表面最适合反应的位点，即使反应势能面最低的位点；从反应的角度看，研究目标是对于一系列反应通道，找到反应势能面上的卡点，找到能够突破卡点的催化环境。
 
-本软件同时进行表面位点的构建、采样、计算、分析，和反应过程的分析采样，即过渡态和分子构象，由此可以实现分子在表面反应的全过程采样。其中的技术关键是采样的效率，我们针对性地使用综合使用了多样化的采样策略，保证表面各类型位点和关键位点的充分采样。
+本软件同时进行表面位点的构建、采样、计算、分析，和反应过程的分析采样，即过渡态和分子构象，由此可以实现分子在表面反应的全过程采样。不仅包括平板模型的表面，还包括分子筛孔道內、团簇表面等。
+
+表面反应过程采样的技术关键是采样的效率，我们针对性地使用综合使用了多样化的采样策略，保证表面各类型位点和关键位点的充分采样。
 
 ### 应用场景
 * 催化剂表面最可能反应位点采样
@@ -41,15 +43,17 @@ Dynamic: summary
   * 优化 C、R、Z （z方向，距离表面的距离）空间，固定COM
   * 返回所采结构中最低能量和对应的结构和能量
   
-* 催化剂表面全势能面采样
+* 催化剂表面全势能面采样（位点、分子构象、旋转）
   * 适用于催化剂位点全局考察
   * 适用于表面分子全局的构象和姿态
   * 适用于反应全过程采样
   * 不优化，或者仅优化Z空间
   * 返回全部能量和结构，以及全局的fitting结果
   * 使用能量截断（dEmax<2eV）
+* 原子力显微镜（AFM）谱图模拟
+* 分子筛孔道內分子构象采样
+
 
-<img src="docs/birds_view.png" alt="表面位点全局分析" style="zoom:50%;" />
 
 ## 程序流程图  Program Workflow
 
@@ -96,24 +100,6 @@ Dynamic: summary
 * MDMin
 * FIRE
 
-各种优化算法的对比，参考[链接](https://wiki.fysik.dtu.dk/gpaw/devel/ase_optimize/ase_optimize.html)
-
-**使用方法**
-
-修改 `surface_reaction_sample.py` 其中的一行 
-
-```
-from ase.optimize import BFGS
-```
-
-改为
-
-```
-from ase.optimize import XXX as BFGS
-```
-
-注意：目前这只是权宜之计，后面会把相应的设置加入到 `parameter.py`
-
 ### Gaussian Process Regression 方法
 
 高斯过程回归 GPR 的优点：
@@ -166,6 +152,8 @@ $$
   * v 0.9.2: debug 支持三斜晶系的表面格点化
   * v 0.9.3: debug GridGenerator
 * v 0.10: 表面终结生成完整并入
+* v0.11：分子构象和旋转采样
+* v0.12：AFM 模拟
 
 **TODO**
 
@@ -174,13 +162,16 @@ $$
 * Low-Energy Region Explorer 方法借鉴[^LoreX]: 划分能量区域， 进行 Delaunay 划分，合并较小的区域为较大的区域。然后针对低能的区域进行额外采样。也可以适用不同的方法，低精度进行低能高能区域识别，然后高精度方法加强。
 * 表面格点使用对称性进行降维。在最初期就引入对称性。参考 adsorbate rotation 思路
 * 表面位点数据库
-* 多原子体系（内坐标受限体系）
+* 覆盖度问题
+  * 表面上每增加一个吸附分子，就更新一次格点（可以仅仅更新 vip points，或者dock points）
+  * 表面上新增 dock points 属性，根据吸附能、位点类型或者用户进行定义。默认是所有的位点，根据能量分析或者用户更新
+  * sg 增加 del_grid(indices), 删除对应的 points， vectors，Dga 等等
 * 完善用户界面、例子、教程
 
 
 ## Reference
 
-[^Duvenaud]: [The Kernel Cookbook: Advice on Covariance functions](https://www.cs.toronto.edu/~duvenaud/)]
+[^Duvenaud]: [The Kernel Cookbook: Advice on Covariance functions](https://www.cs.toronto.edu/~duvenaud/)
 [^BASC]: Shane Carr, Roman Garnett, Cynthia LoBASC: Applying Bayesian Optimization to the Search for Global Minima on Potential Energy Surfaces.
 
 [^向量空间转化]: 计算实空间和向量空间的相邻格点距离的映射系数，根据此系数将实空间的距离转化为向量空间距离。

{surface_construct-0.12 → surface_construct-0.12.1}/tests/test_surface_grid.py

@@ -118,4 +118,12 @@ class TestSurfaceGrid:
         atoms = ase.io.read('atoms_files/Fe28C12_115_12_POSCAR')
         sg_obj = SurfaceGrid(atoms, interval=0.1, ads_num=6, lpca=False)
         sg_obj.gridize(subtype='slab')
-        print(sg_obj.points.shape)
+        print(sg_obj.points.shape)
+
+    def test_ads_C4H4O(self):
+        atoms = ase.io.read('atoms_files/C4H4O_Cu.vasp')
+        sg_obj = SurfaceGrid(atoms, interval=0.2, ads_num=6, lpca=False)
+        sg_obj.gridize(subtype='slab')
+        sg_obj.vectorize()
+        assert sg_obj.vector.shape==(6582, 40)
+        pass