24小时热门版块排行榜    

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

chk0521

金虫 (初入文坛)

[交流] 【求助】VASP加电场的具体计算步骤是怎么样的?已有14人参与

VASP加电场的具体计算步骤是否是:

1.先在未加电场的情况下进行几何优化
2.几何优化完后,进行静态自洽计算得到WAVECAR和CHGCAR
3.WAVECAR和CHGCAR的基础上,在INCAR中加入
LDIPOL=.TRUE.
IDIPOL=1/2/3
EFIELD= ×× (units eV/A)再进行几何优化。

请高人指点一下我上面的步骤是否正确?

我用上面的步骤计算时,在第3步加入电场后几何优化时电子步的能量来回震荡不收敛,这是什么原因呢?INCAR中还需要特别加入其他什么参数吗?我用的vasp版本是vasp.4.6.28。

请高人帮忙解答,万分感谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhsh

木虫 (小有名气)

这是相应的力的收敛情况:

grep 'FORCE t' ../OUTCAR
  FORCE total and by dimension    1.553955    0.387786
  FORCE total and by dimension    3.782150    0.840370
  FORCE total and by dimension    0.868240    0.220130
  FORCE total and by dimension    1.089236    0.224087
  FORCE total and by dimension    2.489097    0.557312
  FORCE total and by dimension    1.109065    0.230254
  FORCE total and by dimension    1.211541    0.255748
  FORCE total and by dimension    0.640010    0.146316
  FORCE total and by dimension    0.592315    0.201554
  FORCE total and by dimension    0.588112    0.198947
  FORCE total and by dimension    0.350382    0.097598
  FORCE total and by dimension    0.375239    0.104847
  FORCE total and by dimension    0.207250    0.051444
  FORCE total and by dimension    0.195840    0.049750
  FORCE total and by dimension    0.137772    0.031733
12楼2010-09-30 08:29:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 32 个回答

215215

银虫 (小有名气)


chk0521(金币+1):谢谢参与
学习学习!静候高人...
nevergiveup!
2楼2010-09-28 15:15:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)

★ ★
chk0521(金币+1):谢谢参与
chk0521(金币+5): 2010-09-29 09:25:34
zxzj05(金币+1):谢谢回帖交流 2010-10-01 20:58:50
我加电场的时候, 就直接把这些参数加入incar中进行结构优化. 加不加电场的基态肯定是不一样的.
3楼2010-09-29 05:18:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chk0521

金虫 (初入文坛)

谢谢您的回复,我知道基态应该不一样。但是我的问题是电子步来回震荡,没有趋于一个稳定值,这是为什么呢?我看到手册上说加电场后电子步收敛的比较慢,但是我这个都算了200多步了还是不收敛,每当快要收敛到一个稳定值,能量就上升到一个很高的值,然后来回震荡。
4楼2010-09-29 09:31:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见