24小时热门版块排行榜    

查看: 2339  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
回复此楼

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

第一性原理相关文档

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 9/450 2026-08-11 22:12 by 医学老男孩
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见