24小时热门版块排行榜    

CyRhmU.jpeg
查看: 5603  |  回复: 30
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

whywhy01

金虫 (小有名气)

[求助] VASP加电场的问题已有1人参与

在金属薄膜体系加电场,有15A的真空层,电场沿Z轴方向
开电场用的下面三个参数:
  EFIELD = 0.05
  IDIPOL = 3
  LDIPOL = .TRUE.
自洽计算电荷密度时发现只有LDIPOL的开关对计算过程和结果有影响,电场的正负和大小都不能改变CHGCAR和WAVECAR,感觉这个结果很不符合常理
这是完整的INCAR:
    PREC = Accurate
   ENCUT = 500
  IBRION = -1
  NELMIN = 2
   EDIFF = 1.0e-05
VOSKOWN = 1
  NBLOCK = 1
    NELM = 60
    ALGO = Normal (blocked Davidson)
   ISPIN = 2
  INIWAV = 1
  ISTART = 0
  ICHARG = 2
ADDGRID = .FALSE.
  ISMEAR = -5
   SIGMA = 0.2
   LREAL = A
  LORBIT = 11
  MAGMOM = 4*0.749 4*0.813 4*1.974
LMAXMIX = 4
    NPAR = 6
  EFIELD = 0.3
  LDIPOL = .TRUE.
  IDIPOL = 3

其他参数应该都没有问题。希望大家能帮我看一看问题出在哪儿。。。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)

引用回帖:
19楼: Originally posted by valenhou001 at 2014-12-04 10:24:09
加减电子数,还不是模拟加静电场。 通过偶极矩、锯齿型势、平行板电容器加静电场会改变上下两个表面上的电荷分布,但是不改变体系的总电荷。

你试过把电场强度增大一些,比如0.2, 0.3或更大一点。...

侯老师,我把金属层加厚到了6层原子,发现电场能比较明显地起作用了。现在看可能是太薄的金属层会限制电场的效果。
29楼2014-12-05 17:35:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 31 个回答

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流。。。 2014-11-30 17:14:04
你怎么判断”电场的正负和大小都不能改变CHGCAR和WAVECAR,“。
你将静电势输出到LOCPOT,然后对LOCPOT处理一下,看加电场后真空区间的势的变化,是不是斜的。
2楼2014-11-29 19:06:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)

引用回帖:
2楼: Originally posted by valenhou001 at 2014-11-29 19:06:10
你怎么判断”电场的正负和大小都不能改变CHGCAR和WAVECAR,“。
你将静电势输出到LOCPOT,然后对LOCPOT处理一下,看加电场后真空区间的势的变化,是不是斜的。

我简单的看了一下CHGCAR和WAVECAR的大小,发现是一样的,可能这个方法不够准确。我试一下用LOCPOT处理的方法,谢谢!
3楼2014-11-30 14:33:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whywhy01

金虫 (小有名气)

引用回帖:
2楼: Originally posted by valenhou001 at 2014-11-29 19:06:10
你怎么判断”电场的正负和大小都不能改变CHGCAR和WAVECAR,“。
你将静电势输出到LOCPOT,然后对LOCPOT处理一下,看加电场后真空区间的势的变化,是不是斜的。

是用vtotav处理LOCPOT出来的VLINE文件看吗?加正负方向电场的曲线是完全重合的
4楼2014-11-30 15:21:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见