计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 301 | 评论数: 7 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-10-15 17:34

正文摘要:

本帖最后由 Wednesday 于 2025-10-15 17:35 编辑 计算了外加电场的Ni(111)slab的静电势,切了4层,扩胞到4x4,一共64个原子,k点考虑的441,然后用abacus自带的后处理分析发现没有平缓期,其中锯齿电势加在分数坐 ...

回复 Reply

ultramanm87 发表于 Post on 2025-10-16 14:38:46
你图里的曲线上的凸起应该是偶极和带电平板的位置,那我觉得应该是它们位置的原因。你可以该方法的参考原文献: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 发表于 Post on 2025-10-16 10:07:42
ultramanm87 发表于 2025-10-16 09:36
例子里说:
  • efield_amp:类型是 Real,电场强度,默认值为 0,单位为 a.u.,其中 1 a.u. = 51.422063 ...

  • 感谢感谢
    ultramanm87 发表于 Post on 2025-10-16 09:36:01
    本帖最后由 ultramanm87 于 2025-10-16 09:44 编辑
    Wednesday 发表于 2025-10-16 08:18
    感谢老师,还有几个问题想请教下,我这个体系肯定是上下不对称的,还需要开偶极矫正么?偶极矫正我看是通 ...
    例子里说:
    • efield_amp:类型是 Real,电场强度,默认值为 0,单位为 a.u.,其中 1 a.u. = 51.4220632*10^10 V/m。在偶极修正时需要设为 0,不为 0 则说明同时添加了外加电场和偶极修正

    那么应该是可以同时的
    Wednesday 发表于 Post on 2025-10-16 08:18:28
    ultramanm87 发表于 2025-10-16 00:07
    偶极校正和消除周期性没啥关系,要拉平加上compensating charge。https://github.com/MCresearch/abacus-us ...

    感谢老师,还有几个问题想请教下,我这个体系肯定是上下不对称的,还需要开偶极矫正么?偶极矫正我看是通过外电场形式去抵消的,那是不是开启偶极矫正时我想要添加的外电场会失效?
    ultramanm87 发表于 Post on 2025-10-16 00:07:48
    偶极校正和消除周期性没啥关系,要拉平加上compensating charge。https://github.com/MCresearch/ab ... /abacus-surface6.md

    手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

    GMT+8, 2026-1-24 16:41 , Processed in 0.176743 second(s), 25 queries , Gzip On.

    快速回复 返回顶部 返回列表 Return to list