24小时热门版块排行榜    

查看: 2334  |  回复: 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的回帖

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的回帖
查看全部 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 15/750 2026-08-10 19:29 by 要乐观耀哥
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 好奇怪的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
[基金申请] 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 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见