计算化学公社

标题: abacus计算静电势没有平缓部分 [打印本页]
作者
Author: Wednesday    时间: 2025-10-15 17:34
标题: abacus计算静电势没有平缓部分
本帖最后由 Wednesday 于 2025-10-15 17:35 编辑

计算了外加电场的Ni(111)slab的静电势，切了4层，扩胞到4x4，一共64个原子，k点考虑的441，然后用abacus自带的后处理分析发现没有平缓期，其中锯齿电势加在分数坐标0.0-0.1，Ni最下方一层的原子分数坐标0.16，锯齿电势是正确加在真空区域的，切面时的真空层也有25埃，但是看这个图像还是受到了周期性影响，请教下各位老师这是什么问题 (, 下载次数 Times of downloads: 1)
作者
Author: ultramanm87    时间: 2025-10-16 00:07
偶极校正和消除周期性没啥关系，要拉平加上compensating charge。https://github.com/MCresearch/ab ... /abacus-surface6.md。
作者
Author: Wednesday    时间: 2025-10-16 08:18
ultramanm87 发表于 2025-10-16 00:07
偶极校正和消除周期性没啥关系，要拉平加上compensating charge。https://github.com/MCresearch/abacus-us ...

感谢老师，还有几个问题想请教下，我这个体系肯定是上下不对称的，还需要开偶极矫正么？偶极矫正我看是通过外电场形式去抵消的，那是不是开启偶极矫正时我想要添加的外电场会失效？
作者
Author: ultramanm87    时间: 2025-10-16 09:36
本帖最后由 ultramanm87 于 2025-10-16 09:44 编辑
Wednesday 发表于 2025-10-16 08:18
感谢老师，还有几个问题想请教下，我这个体系肯定是上下不对称的，还需要开偶极矫正么？偶极矫正我看是通 ...
例子里说：

那么应该是可以同时的
作者
Author: Wednesday    时间: 2025-10-16 10:07
ultramanm87 发表于 2025-10-16 09:36
例子里说：
  • efield_amp：类型是 Real，电场强度，默认值为 0，单位为 a.u.，其中 1 a.u. = 51.422063 ...

    • 感谢感谢
    作者
    Author:     Wednesday    时间: 2025-10-16 10:21
    ultramanm87 发表于 2025-10-16 09:36
    例子里说：
  • efield_amp：类型是 Real，电场强度，默认值为 0，单位为 a.u.，其中 1 a.u. = 51.422063 ...

    • 老师，增加了参数后还是有问题，能麻烦您再帮我看一下是哪里有问题么？
    作者
    Author:     ultramanm87    时间: 2025-10-16 14:38
    你图里的曲线上的凸起应该是偶极和带电平板的位置，那我觉得应该是它们位置的原因。你可以该方法的参考原文献：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也可以尝试）。
    作者
    Author:     Wednesday    时间: 2025-10-20 08:18
    ultramanm87 发表于 2025-10-16 14:38
    你图里的曲线上的凸起应该是偶极和带电平板的位置，那我觉得应该是它们位置的原因。你可以该方法的参考原文 ...

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



    欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3