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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP优化不成功，找不到原因
51/1返回列表
查看: 2351  |  回复: 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

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

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的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +8 医学老男孩 2026-08-13 19/950 2026-08-13 20:15 by yq199112088
[考博] 澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域 +25 AI8RaGaPaSCR 2026-08-07 26/1300 2026-08-13 20:15 by YSLfXhguwSL6
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
信息提示
关闭
请填处理意见
关闭