abacus计算静电势没有平缓部分

Wednesday

计算了外加电场的Ni(111)slab的静电势，切了4层，扩胞到4x4，一共64个原子，k点考虑的441，然后用abacus自带的后处理分析发现没有平缓期，其中锯齿电势加在分数坐标0.0-0.1，Ni最下方一层的原子分数坐标0.16，锯齿电势是正确加在真空区域的，切面时的真空层也有25埃，但是看这个图像还是受到了周期性影响，请教下各位老师这是什么问题

ultramanm87

偶极校正和消除周期性没啥关系，要拉平加上compensating charge。https://github.com/MCresearch/ab ... /abacus-surface6.md。

Wednesday

ultramanm87 发表于 2025-10-16 00:07
偶极校正和消除周期性没啥关系，要拉平加上compensating charge。https://github.com/MCresearch/abacus-us ...

感谢老师，还有几个问题想请教下，我这个体系肯定是上下不对称的，还需要开偶极矫正么？偶极矫正我看是通过外电场形式去抵消的，那是不是开启偶极矫正时我想要添加的外电场会失效？

ultramanm87

Wednesday 发表于 2025-10-16 08:18
感谢老师，还有几个问题想请教下，我这个体系肯定是上下不对称的，还需要开偶极矫正么？偶极矫正我看是通 ...

例子里说：

• efield_amp：类型是 Real，电场强度，默认值为 0，单位为 a.u.，其中 1 a.u. = 51.4220632*10^10 V/m。在偶极修正时需要设为 0，不为 0 则说明同时添加了外加电场和偶极修正
  

那么应该是可以同时的

Wednesday

ultramanm87 发表于 2025-10-16 09:36
例子里说：

efield_amp：类型是 Real，电场强度，默认值为 0，单位为 a.u.，其中 1 a.u. = 51.422063 ...

感谢感谢

Wednesday

ultramanm87 发表于 2025-10-16 09:36
例子里说：

efield_amp：类型是 Real，电场强度，默认值为 0，单位为 a.u.，其中 1 a.u. = 51.422063 ...

老师，增加了参数后还是有问题，能麻烦您再帮我看一下是哪里有问题么？

ultramanm87

你图里的曲线上的凸起应该是偶极和带电平板的位置，那我觉得应该是它们位置的原因。你可以该方法的参考原文献：Electrochemical doping of few-layer ZrNCl from first principles: Electronic and structural properties in field-effect configuration和QE的例子（与abacus中实现一致）https://gitlab.com/QEF/q-e/-/blo ... ngle_+0.10.bands.in （ zgate       = 0.011, emaxpos （对应 efield_pos_max）    = 0.000001,）。据我观察应该是它们可以放在一起而不是分开（参考qe例子的设置和文献的图二以及方法部分的最后一段，要不要加block也可以尝试）。

Wednesday

ultramanm87 发表于 2025-10-16 14:38
你图里的曲线上的凸起应该是偶极和带电平板的位置，那我觉得应该是它们位置的原因。你可以该方法的参考原文 ...

老师，我最后看到了相关的文献，去掉了外电场使用补偿电荷去模拟恒电位了，但是我发现最后静电势的结果有两个平台区，12-20的是补偿电荷板所在，21以后是偶极矫正位置，请问计算功函数时应该选取哪个值作为真空电势呢？