24小时热门版块排行榜    

查看: 11199  |  回复: 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

木虫 (小有名气)

8 F= -.42852002E+03 E0= -.42852002E+03  d E =-.180923E+00  mag=     6.8367
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.428493520979E+03   -0.42849E+03   -0.14089E+00  3448   0.754E+00    0.549E-01
DAV:   2    -0.428537349688E+03   -0.43829E-01   -0.43742E-01  3448   0.138E+00    0.423E-01
DAV:   3    -0.428583118184E+03   -0.45768E-01   -0.28996E-02  3992   0.607E-01    0.716E-01
DAV:   4    -0.428431343877E+03    0.15177E+00   -0.14322E-02  3928   0.364E-01    0.474E-01
DAV:   5    -0.428534043777E+03   -0.10270E+00   -0.22263E-03  3640   0.162E-01    0.203E-01
DAV:   6    -0.428544061511E+03   -0.10018E-01   -0.19082E-03  3544   0.184E-01    0.217E-01
DAV:   7    -0.428544207702E+03   -0.14619E-03   -0.25122E-03  3736   0.209E-01    0.137E-01
DAV:   8    -0.428553761990E+03   -0.95543E-02   -0.11437E-03  3480   0.149E-01    0.101E-01
DAV:   9    -0.428551606202E+03    0.21558E-02   -0.26775E-04  2816   0.772E-02    0.377E-02
DAV:  10    -0.428555932216E+03   -0.43260E-02   -0.10605E-04  2328   0.352E-02
   9 F= -.42855593E+03 E0= -.42855593E+03  d E =-.359095E-01  mag=     6.9330
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.428559886058E+03   -0.42856E+03   -0.96774E-05  2264   0.521E-02    0.361E-02
DAV:   2    -0.428563769600E+03   -0.38835E-02   -0.39471E-04  2720   0.108E-01    0.925E-02
DAV:   3    -0.428564431416E+03   -0.66182E-03   -0.25486E-04  2512   0.779E-02    0.270E-02
DAV:   4    -0.428566650894E+03   -0.22195E-02   -0.44587E-05  2184   0.207E-02    0.314E-02
DAV:   5    -0.428570210090E+03   -0.35592E-02   -0.80388E-05  2272   0.214E-02    0.341E-02
DAV:   6    -0.428575877792E+03   -0.56677E-02   -0.92622E-05  2320   0.275E-02    0.227E-02
DAV:   7    -0.428582738357E+03   -0.68606E-02   -0.88411E-05  2408   0.336E-02    0.253E-02
DAV:   8    -0.428588961256E+03   -0.62229E-02   -0.14081E-04  2424   0.338E-02    0.128E-02
DAV:   9    -0.428593463499E+03   -0.45022E-02   -0.98978E-05  2376   0.253E-02    0.581E-03
DAV:  10    -0.428595180851E+03   -0.17174E-02   -0.32054E-05  2280   0.183E-02
  10 F= -.42859518E+03 E0= -.42859518E+03  d E =-.751582E-01  mag=     7.0652
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.428565754789E+03   -0.42857E+03   -0.11942E+00  3328   0.688E+00    0.502E-01
DAV:   2    -0.428594201604E+03   -0.28447E-01   -0.36405E-01  3384   0.123E+00    0.413E-01
DAV:   3    -0.428607626468E+03   -0.13425E-01   -0.34906E-02  3968   0.515E-01    0.963E-01
DAV:   4    -0.428597426387E+03    0.10200E-01   -0.17902E-02  3872   0.335E-01    0.188E-01
DAV:   5    -0.428598518418E+03   -0.10920E-02   -0.12241E-03  3056   0.133E-01    0.167E-01
DAV:   6    -0.428599256484E+03   -0.73807E-03   -0.15638E-03  3376   0.178E-01    0.142E-01
DAV:   7    -0.428599967441E+03   -0.71096E-03   -0.15082E-03  3432   0.165E-01    0.107E-01
DAV:   8    -0.428600169659E+03   -0.20222E-03   -0.41109E-04  2824   0.921E-02    0.309E-02
DAV:   9    -0.428600251055E+03   -0.81395E-04   -0.58568E-05  2200   0.435E-02
  11 F= -.42860025E+03 E0= -.42860025E+03  d E =-.507020E-02  mag=     7.0681
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.428599598521E+03   -0.42860E+03   -0.28477E-02  3344   0.106E+00    0.787E-02
DAV:   2    -0.428600150280E+03   -0.55176E-03   -0.88255E-03  3312   0.197E-01    0.102E-01
DAV:   3    -0.428600412113E+03   -0.26183E-03   -0.15557E-03  2848   0.107E-01    0.184E-01
DAV:   4    -0.428600207899E+03    0.20421E-03   -0.63007E-04  2448   0.571E-02    0.337E-02
DAV:   5    -0.428600352055E+03   -0.14416E-03   -0.47989E-05  2048   0.324E-02    0.293E-02
DAV:   6    -0.428600437960E+03   -0.85905E-04   -0.33050E-05  2168   0.284E-02
  12 F= -.42860044E+03 E0= -.42860044E+03  d E =-.525711E-02  mag=     7.0692
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.428595653637E+03   -0.42860E+03   -0.29825E-01  3296   0.320E+00    0.222E-01
DAV:   2    -0.428600295553E+03   -0.46419E-02   -0.69096E-02  3344   0.558E-01    0.153E-01
DAV:   3    -0.428601654412E+03   -0.13589E-02   -0.27325E-03  3464   0.217E-01    0.175E-01
DAV:   4    -0.428601419702E+03    0.23471E-03   -0.11103E-03  2976   0.142E-01    0.952E-02
DAV:   5    -0.428601654821E+03   -0.23512E-03   -0.64752E-04  2520   0.666E-02    0.130E-01
DAV:   6    -0.428601853013E+03   -0.19819E-03   -0.60945E-04  2648   0.904E-02    0.695E-02
DAV:   7    -0.428602078456E+03   -0.22544E-03   -0.33366E-04  2600   0.846E-02    0.554E-02
DAV:   8    -0.428602118303E+03   -0.39848E-04   -0.15261E-04  2456   0.702E-02
  13 F= -.42860212E+03 E0= -.42860212E+03  d E =-.168034E-02  mag=     7.0692
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.428602122040E+03   -0.42860E+03   -0.98232E-04  3280   0.180E-01    0.145E-02
DAV:   2    -0.428602161965E+03   -0.39925E-04   -0.22469E-04  2536   0.350E-02
  14 F= -.42860216E+03 E0= -.42860216E+03  d E =-.172401E-02  mag=     7.0692
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.428600638057E+03   -0.42860E+03   -0.80203E-02  3272   0.172E+00    0.126E-01
DAV:   2    -0.428602040838E+03   -0.14028E-02   -0.21843E-02  3320   0.316E-01    0.138E-01
DAV:   3    -0.428602758851E+03   -0.71801E-03   -0.32073E-03  3184   0.160E-01    0.278E-01
DAV:   4    -0.428602237849E+03    0.52100E-03   -0.16957E-03  2824   0.923E-02    0.566E-02
DAV:   5    -0.428602484087E+03   -0.24624E-03   -0.12020E-04  2120   0.498E-02    0.408E-02
DAV:   6    -0.428602596061E+03   -0.11197E-03   -0.99174E-05  2184   0.494E-02    0.373E-02
DAV:   7    -0.428602666672E+03   -0.70611E-04   -0.84209E-05  2176   0.479E-02
  15 F= -.42860267E+03 E0= -.42860267E+03  d E =-.504707E-03  mag=     7.0693
11楼2010-09-30 08:24:13
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见