24小时热门版块排行榜    

CyRhmU.jpeg
查看: 11046  |  回复: 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。

请高人帮忙解答,万分感谢!
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 DFT加电场 第一性 计算—磁性
vasp VASP vasp电场 wly的好资源专辑
VASP计算资料收集byZF

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhsh

木虫 (小有名气)

★ ★ ★
chk0521(金币+1):谢谢参与
chk0521(金币+5): 2010-10-01 03:29:56
chk0521(金币+4): 2010-10-01 03:36:55
zxzj05(金币+2):谢谢回帖交流 2010-10-01 20:59:45
一个建议,做结构优化时,电子步可以取小一点,比如20步,或30步。不需要等到电子完全自洽,因为刚开始时几何结构不大合理,收敛会非常慢。但是你会发现,随着离子步的进行,电子步的自洽会越来越快,这是因为几何结构越来越合理了。当然,也许有人会反驳,只有完全自洽了才能得到这种构型下的真正的力,我不否认。但是在没有完全自洽的情况下让其走一步,尽管所得的力的方向不是“真正”的力的方向,让其动一动,调整一下“姿势”,结构相对前一步肯定会合理一些。这是我的经验,大家要是觉得有些道理,可以试一试。以下是我的一个例子: OSZICAR文件。我的电子步只取了10步。

       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.120638908206E+04    0.12064E+04   -0.80248E+04  4696   0.572E+02
DAV:   2    -0.229803875256E+03   -0.14362E+04   -0.13873E+04  4496   0.220E+02
DAV:   3    -0.453834997226E+03   -0.22403E+03   -0.21777E+03  3992   0.844E+01
DAV:   4    -0.464593354058E+03   -0.10758E+02   -0.99729E+01  4048   0.165E+01
DAV:   5    -0.464856331379E+03   -0.26298E+00   -0.46611E+00  4232   0.289E+00    0.533E+01
DAV:   6    -0.558904913888E+03   -0.94049E+02   -0.53982E+02  3960   0.845E+01    0.870E+01
DAV:   7    -0.432466329466E+03    0.12644E+03   -0.24277E+02  3968   0.983E+01    0.254E+01
DAV:   8    -0.439973166657E+03   -0.75068E+01   -0.61342E+01  4200   0.265E+01    0.467E+01
DAV:   9    -0.432947586448E+03    0.70256E+01   -0.43006E+01  4128   0.178E+01    0.233E+01
DAV:  10    -0.428479765682E+03    0.44678E+01   -0.63480E+00  4168   0.939E+00
   1 F= -.42847977E+03 E0= -.42847977E+03  d E =-.428480E+03  mag=     3.3222
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.429644510970E+03   -0.42964E+03   -0.22612E+01  3784   0.179E+01    0.212E+01
DAV:   2    -0.513735058730E+03   -0.84091E+02   -0.33368E+02  4264   0.508E+01    0.641E+01
DAV:   3    -0.437538753741E+03    0.76196E+02   -0.78457E+01  4152   0.404E+01    0.366E+01
DAV:   4    -0.430430830315E+03    0.71079E+01   -0.39526E+01  4336   0.212E+01    0.247E+01
DAV:   5    -0.428826961761E+03    0.16039E+01   -0.10748E+01  4136   0.932E+00    0.196E+01
DAV:   6    -0.428018879377E+03    0.80808E+00   -0.60342E+00  4144   0.549E+00    0.124E+01
DAV:   7    -0.428091220364E+03   -0.72341E-01   -0.58288E+00  3992   0.532E+00    0.109E+01
DAV:   8    -0.427890173081E+03    0.20105E+00   -0.15282E+00  3840   0.377E+00    0.471E+00
DAV:   9    -0.427928524555E+03   -0.38351E-01   -0.16457E-01  3864   0.136E+00    0.535E+00
DAV:  10    -0.427987697903E+03   -0.59173E-01   -0.31814E-02  3560   0.536E-01
   2 F= -.42798770E+03 E0= -.42798770E+03  d E =0.492068E+00  mag=     3.5314
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.427702311785E+03   -0.42770E+03   -0.65016E+00  3368   0.105E+01    0.554E+00
DAV:   2    -0.434220199992E+03   -0.65179E+01   -0.39630E+01  3904   0.114E+01    0.392E+01
DAV:   3    -0.427634428803E+03    0.65858E+01   -0.37155E+01  3944   0.925E+00    0.821E+00
DAV:   4    -0.428464850653E+03   -0.83042E+00   -0.35202E+00  3952   0.536E+00    0.122E+01
DAV:   5    -0.427849707030E+03    0.61514E+00   -0.16658E+00  3776   0.308E+00    0.367E+00
DAV:   6    -0.427808857524E+03    0.40850E-01   -0.24713E-01  3832   0.161E+00    0.273E+00
DAV:   7    -0.427747627481E+03    0.61230E-01   -0.30396E-02  3760   0.616E-01    0.348E+00
DAV:   8    -0.427655076314E+03    0.92551E-01   -0.21309E-02  3608   0.492E-01    0.411E+00
DAV:   9    -0.427643333254E+03    0.11743E-01   -0.32473E-03  3672   0.157E-01    0.420E+00
DAV:  10    -0.427640926731E+03    0.24065E-02   -0.10076E-03  3504   0.131E-01
   3 F= -.42764093E+03 E0= -.42764093E+03  d E =0.838839E+00  mag=     3.9323
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.427593204938E+03   -0.42759E+03   -0.59530E-02  3480   0.100E+00    0.373E+00
DAV:   2    -0.428282613001E+03   -0.68941E+00   -0.47857E+00  3848   0.471E+00    0.138E+01
DAV:   3    -0.427836820829E+03    0.44579E+00   -0.92126E+00  3904   0.542E+00    0.801E+00
DAV:   4    -0.427829538706E+03    0.72821E-02   -0.94823E-01  3896   0.356E+00    0.540E+00
DAV:   5    -0.427664222645E+03    0.16532E+00   -0.47511E-01  3760   0.226E+00    0.419E+00
DAV:   6    -0.427597829841E+03    0.66393E-01   -0.14340E-01  3952   0.118E+00    0.289E+00
DAV:   7    -0.427323642590E+03    0.27419E+00   -0.12967E-01  3688   0.120E+00    0.177E+00
DAV:   8    -0.427493460239E+03   -0.16982E+00   -0.64766E-02  3768   0.913E-01    0.162E+00
DAV:   9    -0.427428708618E+03    0.64752E-01   -0.30047E-02  3832   0.751E-01    0.123E+00
DAV:  10    -0.427551745208E+03   -0.12304E+00   -0.47309E-02  3728   0.704E-01
   4 F= -.42755175E+03 E0= -.42755175E+03  d E =0.928020E+00  mag=     5.1253
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.427465976634E+03   -0.42747E+03   -0.86630E+00  3528   0.146E+01    0.172E+00
DAV:   2    -0.427728718046E+03   -0.26274E+00   -0.24874E+00  3640   0.403E+00    0.642E+00
DAV:   3    -0.427586000266E+03    0.14272E+00   -0.96676E-01  4072   0.300E+00    0.208E+00
DAV:   4    -0.427489807897E+03    0.96192E-01   -0.10590E-01  3776   0.108E+00    0.112E+00
DAV:   5    -0.427463111443E+03    0.26696E-01   -0.37656E-02  3880   0.794E-01    0.823E-01
DAV:   6    -0.427504354524E+03   -0.41243E-01   -0.25528E-02  4104   0.622E-01    0.754E-01
DAV:   7    -0.427666933588E+03   -0.16258E+00   -0.37745E-02  3832   0.655E-01    0.655E-01
DAV:   8    -0.427788306719E+03   -0.12137E+00   -0.38314E-02  3880   0.691E-01    0.854E-01
DAV:   9    -0.427758917530E+03    0.29389E-01   -0.38824E-03  3728   0.271E-01    0.614E-01
DAV:  10    -0.427821171479E+03   -0.62254E-01   -0.20321E-02  4096   0.276E-01
   5 F= -.42782117E+03 E0= -.42782117E+03  d E =-.269426E+00  mag=     5.9283
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.427991523662E+03   -0.42799E+03   -0.37230E+00  3336   0.957E+00    0.960E-01
DAV:   2    -0.428117088376E+03   -0.12556E+00   -0.81263E-01  3608   0.250E+00    0.301E+00
DAV:   3    -0.428132054084E+03   -0.14966E-01   -0.32426E-01  3960   0.146E+00    0.136E+00
DAV:   4    -0.428133401584E+03   -0.13475E-02   -0.79143E-02  3928   0.824E-01    0.826E-01
DAV:   5    -0.428131342803E+03    0.20588E-02   -0.28406E-02  3904   0.559E-01    0.615E-01
DAV:   6    -0.428132126955E+03   -0.78415E-03   -0.13076E-02  3944   0.337E-01    0.497E-01
DAV:   7    -0.428270138552E+03   -0.13801E+00   -0.80236E-03  3768   0.384E-01    0.482E-01
DAV:   8    -0.428335816619E+03   -0.65678E-01   -0.53490E-03  3984   0.280E-01    0.359E-01
DAV:   9    -0.428219700913E+03    0.11612E+00   -0.12957E-03  3840   0.157E-01    0.440E-01
DAV:  10    -0.428339099893E+03   -0.11940E+00   -0.63678E-03  3648   0.235E-01
   6 F= -.42833910E+03 E0= -.42833910E+03  d E =-.787355E+00  mag=     6.3860
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.428188692737E+03   -0.42819E+03   -0.58311E+00  3424   0.142E+01    0.115E+00
DAV:   2    -0.428287075179E+03   -0.98382E-01   -0.15430E+00  3528   0.294E+00    0.258E+00
DAV:   3    -0.428376614244E+03   -0.89539E-01   -0.39444E-01  4096   0.167E+00    0.195E+00
DAV:   4    -0.428363706597E+03    0.12908E-01   -0.97121E-02  3904   0.851E-01    0.789E-01
DAV:   5    -0.428359816524E+03    0.38901E-02   -0.21249E-02  3952   0.552E-01    0.422E-01
DAV:   6    -0.428359153743E+03    0.66278E-03   -0.65648E-03  3912   0.361E-01    0.391E-01
DAV:   7    -0.428366027548E+03   -0.68738E-02   -0.89284E-03  3816   0.371E-01    0.361E-01
DAV:   8    -0.428389369844E+03   -0.23342E-01   -0.52421E-03  3888   0.335E-01    0.269E-01
DAV:   9    -0.428397677880E+03   -0.83080E-02   -0.17481E-03  3544   0.208E-01    0.113E-01
DAV:  10    -0.428419912503E+03   -0.22235E-01   -0.40205E-04  2992   0.669E-02
   7 F= -.42841991E+03 E0= -.42841991E+03  d E =-.808126E-01  mag=     6.7019
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.428452849579E+03   -0.42845E+03   -0.90597E-01  3400   0.561E+00    0.446E-01
DAV:   2    -0.428483565767E+03   -0.30716E-01   -0.24001E-01  3512   0.118E+00    0.105E+00
DAV:   3    -0.428515424415E+03   -0.31859E-01   -0.62955E-02  4016   0.670E-01    0.755E-01
DAV:   4    -0.428384970936E+03    0.13045E+00   -0.14913E-02  3952   0.371E-01    0.465E-01
DAV:   5    -0.428494186179E+03   -0.10922E+00   -0.41756E-03  3784   0.217E-01    0.236E-01
DAV:   6    -0.428504653158E+03   -0.10467E-01   -0.18715E-03  3568   0.188E-01    0.172E-01
DAV:   7    -0.428506388879E+03   -0.17357E-02   -0.15634E-03  3656   0.200E-01    0.174E-01
DAV:   8    -0.428397279822E+03    0.10911E+00   -0.13273E-03  3488   0.162E-01    0.300E-01
DAV:   9    -0.428510355965E+03   -0.11308E+00   -0.10393E-03  3272   0.123E-01    0.119E-01
DAV:  10    -0.428520022686E+03   -0.96667E-02   -0.79978E-04  3176   0.931E-02
10楼2010-09-30 08:24:02
已阅   回复此楼   关注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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见