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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】VASP加电场的具体计算步骤是怎么样的？
51/1返回列表
查看: 12814  |  回复: 31
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-09-28 15:13:04
已阅   回复此楼   关注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
赞一下(5人) 回复此楼
10楼2010-09-30 08:24:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 32 个回答

215215

银虫 (小有名气)
★
chk0521(金币+1):谢谢参与
学习学习！静候高人...
赞一下(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中进行结构优化. 加不加电场的基态肯定是不一样的.
赞一下(2人) 回复此楼
3楼2010-09-29 05:18:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chk0521

金虫 (初入文坛)
谢谢您的回复，我知道基态应该不一样。但是我的问题是电子步来回震荡，没有趋于一个稳定值，这是为什么呢？我看到手册上说加电场后电子步收敛的比较慢，但是我这个都算了200多步了还是不收敛，每当快要收敛到一个稳定值，能量就上升到一个很高的值，然后来回震荡。
赞一下(1人) 回复此楼
4楼2010-09-29 09:31:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 32 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +10 Kittylucky 2026-08-27 11/550 2026-08-28 10:45 by Kittylucky
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
信息提示
关闭
请填处理意见
关闭