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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP优化不成功，找不到原因
91/1返回列表
查看: 2284  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

csjfding

新虫 (小有名气)

[求助] VASP优化不成功，找不到原因

我用VASP来做驰豫，incar中的参数是下面这样设置的，运行了好几天，计算中没报错，可从驰豫后的OUTCAR中发现驰豫好像没最后完成。不知道原因出在哪？希望大家帮帮忙。谢谢！
sys=Mo4O12

Startparameter for this Run:
   NWRITE =      2;   LPETIM=.TRUE.    write-flag & timer
   ISTART =      0    job   : 0-new  1-cont  2-samecut
   LCHARG =      T
   LWAVE  =      F
   ISPIN  =      1
   ISYM   =      2
   PREC   =      A
   ENCUT  =      520

Electronic Relaxation
   NELM   =      60;
   EDIFF  =      1E-04   stopping-criterion for ELM
   LREAL  =      A   #Auto: N_ion>20 Test RPOT; F(alse) : N_ion < 20

Ionic Relaxation
   IBRION =      2    ionic relax: 0-MD 1-quasi-New 2-CG
   ISIF   =      2    relax z axiz
   NSW    =      500
   EDIFFG =      -1E-03
   POTIM  =      0.50

DOS related values:
   ISMEAR =     0
   SIGMA  =     0.20

Algorithm:
   ALGO  =     FAST  #Normal: Default;  use Fast for N_ion > 20
回复此楼

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

第一性原理相关文档

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-02-23 09:30:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csjfding

新虫 (小有名气)
另外，OUTCAR的最后是这样的：
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      1.98482      6.14018      0.92410         0.005673      0.000510      0.000000
      2.06219      1.25279      0.92410        -0.001134     -0.002683      0.000000
      0.12291      3.14644      0.92410        -0.001420      0.002226      0.000000
      3.95321      0.91244      2.77230        -0.005605      0.000354      0.000000
      0.06950      5.72390      2.77230         0.006083     -0.000140      0.000000
      2.08853      3.80669      2.77230         0.003194     -0.003781      0.000000
      1.88994     12.60393      2.77230        -0.010795      0.000467      0.000000
      3.82052     10.70233      2.77230        -0.002636      0.000654      0.000000
     -0.00099     12.95302      0.92410        -0.002106     -0.000253      0.000000
      3.78707      8.09622      0.92410         0.005029     -0.000106      0.000000
      1.85596     10.00875      0.92410         0.000718     -0.000804      0.000000
      0.30614      1.45067      0.92410        -0.003371     -0.004212      0.000000
      2.25699      5.49980      2.77230         0.003039      0.003812      0.000000
      3.64718     12.39233      2.77230         0.004993     -0.001382      0.000000
      1.61945      8.32384      0.92410        -0.001661      0.005337      0.000000
-----------------------------------------------------------------------------------
    total drift:                               -0.005357     -0.009776      0.000000
赞一下 回复此楼
2楼2012-02-23 09:31:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxf_dxf_dxf

银虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
uuv2010(金币+1): 多谢提示 2012-02-23 12:51:41
uuv2010(金币+2): 详细的解答，赞！ 2012-02-23 12:52:21
你设置了NSW=500，如果500步内还没有打到收敛标准，计算会自动停止。
建议：
先把ECUT、KPOINTS设的小一点，还有EDIFF和EDIFFG设的大一点，总之对精度要求先小一点，收敛之后得到的CONTCAR拷贝称POSCAR，把计算精度提高，继续优化。
分几个台阶来进行优化，这样收敛的希望更大一些。

附：VASP官方网站上Workshop内的提示
What can one do when convergence is bad?
1 Start from charge density of non-spin-polarized calculation, using
ISTART = 0 (or remove WAVECAR)
ICHARG = 1
2 Linear mixing
BMIX = 0.0001 ; BMIX MAG = 0.0001
3 Mix slowly, i.e., reduce AMIX and AMIX MAG
4 Reduce MAXMIX, the number of steps stored in the Broyden mixer
(default = 45)
5 Restart from partly converged results
(stop a calculation after say 20 steps and restart from the WAVECAR)
6 Use constraints to stabilize the magetic configuration
7 Pray
赞一下 回复此楼
3楼2012-02-23 11:17:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)
uuv2010(金币+1): 多谢提示 2012-02-23 12:51:50
查看下例子优化到底进行了几步，如果是500，说明是没有完全优化。
不过这样也说明初始结构很差。。
赞一下 回复此楼
timeflies..
4楼2012-02-23 12:20:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guodonglin

至尊木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
uuv2010(金币+1): 多谢提示 2012-02-23 12:51:56
500步还没有收敛，可能结构有问题，请仔细检查POSCAR
赞一下 回复此楼
5楼2012-02-23 12:23:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csjfding

新虫 (小有名气)
谢谢各位的回复。我的结构和POSCAR都是直接通过MS得到的。
刚看了一下OUTCAR，应该是循环了217步结束的。
里面还多次出现这句话：aborting loop because EDIFF is reached
不知道怎么粘附件，不然粘上来让大家帮我看看
赞一下 回复此楼
6楼2012-02-23 14:31:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

llh2010

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
liliangfang(金币+1): 多谢指教 2012-02-24 13:50:56
要把运行输出文件，就是开头是VSAP版本及读取输入文件内容然后是电子步和离子步循环的那个输出文件（不是OUTCAR）贴出来看看。即使是没有完成结构优化，程序正常结束的话，OUTCAR的最后几行是给出程序总运行时间和使用的内存情况。从你的贴出的OUTCAR看应该是程序没有正常结束。
赞一下 回复此楼
知识引导人生，学习成就未来
7楼2012-02-24 11:47:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rxd337

木虫 (正式写手)

养殖场经理

【答案】应助回帖

感谢参与，应助指数 +1
liliangfang(金币+1): 多谢指教 2012-02-24 18:38:48
看OUTCAR，看弛豫到第几部了。你的情况可能是程序崩溃了，还是任务被杀掉了？反正没有自然停止。如果是程序崩溃，可能是计算量太大，你要重新配置计算节点
赞一下 回复此楼
都是假设！
8楼2012-02-24 18:22:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
结构优化一般60步，超过了一般就很难收敛，你可以调小一点收敛精度
赞一下 回复此楼
9楼2012-02-24 18:41:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 csjfding 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 5个%2F，非固定段22个字母 +9 jaasxulu 2026-07-15 10/500 2026-07-19 18:58 by wowsunflower
[基金申请] 小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子 +15 botar 2026-07-14 16/800 2026-07-19 18:23 by wowsunflower
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 bjm2e4yior 2026-07-18 3/150 2026-07-19 16:43 by 0u3sh5hhfi
[论文投稿] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 6/300 2026-07-19 16:23 by 0u3sh5hhfi
[考博] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 4/200 2026-07-19 16:10 by 0u3sh5hhfi
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 4/200 2026-07-19 15:54 by 0u3sh5hhfi
[基金申请] 国社科通讯评审 +5 潜水冠军123 2026-07-18 6/300 2026-07-19 12:28 by lfy8008
[基金申请] 生命口C13面上会评时间 +7 luckinging 2026-07-14 11/550 2026-07-17 20:43 by 虫卜1988
[基金申请] 生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？ +6 碧水00 2026-07-16 9/450 2026-07-17 15:21 by jaasxulu
[论文投稿] 投英文期刊如何查重 10+5 不知道取什么名� 2026-07-15 6/300 2026-07-17 09:07 by bobvan
[有机交流] 伯胺上甲基求助 100+4 002tmumu 2026-07-14 7/350 2026-07-16 19:47 by 002tmumu
[论文投稿] 跨出版社商投稿 20+5 倪好520 2026-07-13 5/250 2026-07-16 17:24 by 无视军团
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[基金申请] H口，面上，时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗？ 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[基金申请] 没收到消息，再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
信息提示
关闭
请填处理意见
关闭