应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP加电场的具体计算步骤是怎么样的？
322/4返回列表上一页1234下一页
查看: 11040  |  回复: 31
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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的回帖

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的回帖

valenhou001

至尊木虫 (职业作家)
★ ★ ★
chk0521(金币+1):谢谢参与
chk0521(金币+3): 2010-10-01 03:32:29
zxzj05(金币+2):谢谢回帖交流 2010-10-01 21:00:00
一开始电场不要加的太大，电场的强度逐渐加上去。比如一开始做一个电场强度很小的计算，得到了收敛的电荷密度和波函数之后，再增加电场强度，读入上一步的电荷和波函数进行计算。。。。。这样的方式逐渐增加电场强度。

另外也要注意，这里所加的电场所能处理的体系，即电场的方向上是要有真空层的。
赞一下(9人) 回复此楼
13楼2010-09-30 08:34:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freer007

木虫 (正式写手)
★ ★
chk0521(金币+1):谢谢参与
chk0521(金币+2): 2010-10-01 03:32:38
zxzj05(金币+1):谢谢回帖交流 2010-10-01 21:00:27
楼主是给什么材料加的电场？能得到什么结果呢？想得到什么结果呢？我以前尝试给块体晶体材料加电场，也做成了加有真空层的slab模型，但是没有得到有价值的结果。（我原本是想看看电场作用下结构有什么变化。）希望楼主能有得到想要的结果。
赞一下(2人) 回复此楼
14楼2010-09-30 10:08:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★
zxzj05(金币+2):谢谢回帖交流 2010-10-01 21:00:39
可以先找找相关的文献，理解vasp里面所实现加电场的方法。它的这种方法处理晶体材料是不行的。
赞一下(2人) 回复此楼
15楼2010-09-30 12:12:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chk0521

金虫 (初入文坛)
多谢各位专家的指点，随着粒子步的进行，电子步确实也逐渐收敛了。真的非常感谢各位！
赞一下 回复此楼
16楼2010-10-01 03:34:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljx2006

木虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
vasp加电场的时候这种情况还是比较常见的，一直也很疑惑，即使一个很小的体系，只要加上电场，就半天收敛不了，有时候勉强计算完成了，但是最后的几个能量就又开始变得比较高。有没有人处理过类似的情况呀，交流一下呀！谢谢了
赞一下(2人) 回复此楼
17楼2011-06-28 21:53:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ffzx

金虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
10楼: Originally posted by xhsh at 2010-09-30 08:24:02:
一个建议，做结构优化时，电子步可以取小一点，比如20步，或30步。不需要等到电子完全自洽，因为刚开始时几何结构不大合理，收敛会非常慢。但是你会发现，随着离子步的进行，电子步的自洽会越来越快，这是因为几何 ...

这样真的可以吗？一般是觉得应该电子步在60步以内自洽的，这样才好。
赞一下(1人) 回复此楼
18楼2011-12-22 13:10:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sg18408926
19楼2011-12-22 15:27:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superest

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
13楼: Originally posted by valenhou001 at 2010-09-30 08:34:32:
一开始电场不要加的太大，电场的强度逐渐加上去。比如一开始做一个电场强度很小的计算，得到了收敛的电荷密度和波函数之后，再增加电场强度，读入上一步的电荷和波函数进行计算。。。。。这样的方式逐渐增加电场强 ...

加电场后,是否还需要弛豫计算?
赞一下 回复此楼
20楼2012-02-16 16:23:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chk0521 的主题更新
322/4返回列表上一页1234下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭