请教区分polar surface和polar slab以及pymatgen模块处理slab的相关问题

Jiongci-21

社友们好，

首先我感觉 polar surface 应该和固体物理和表面物理相关，所以斗胆直接在物理学区开贴，如果有什么不对的地方还请批评指正。

我最近在用pymatgen.core.surface module 处理生成一些rock salt晶体的表面。以经典的NaCl为例，当我切了(1 1 1)面并生成了一定厚度的surface slab model（如图1，slab上下面非对称），调用is_polar()函数的时候，结果显示是极性的（is_polar()介绍说这个是检查slab是否是极性的）。但是当我添加了symmetrize=True 参数后，生成的slab被强制为对称的上下面（如图2），is_polar()的结果就是非极性了。我能理解NaCl的(1 1 1)面由于表面存在net dipole density肯定是极性面，但是我不太理解 polar surface 和 polar slab 之间的关系和区别。J. Neugebauer and M. Scheffler, Phys. Rev. B 46, 16067 (1992) 这篇文章里提到了asymmetric slab，是否只要是非对称的slab都是polar slab？


我接下来要对这些slab进行DFT计算优化，请问我该怎么判别哪些slab需要添加dipole correction？我最大的疑惑是，是否只要是polar surface 就必须添加dipole correction？还是说我做了symmetrize=True的处理后，得到non-polar slab，就不需要添加dipole correction了呢？


希望有懂的社友多多赐教，提供一些阅读的材料，万分感谢。


   
图 1 symmetrize=False 以及生成的NaCl (111) slab。


图 2 symmetrize=False 以及生成Na-terminal 和 Cl-terminal 的两个NaCl (111) 上下对称的slabs。