24小时热门版块排行榜    

查看: 2337  |  回复: 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 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +16 Lanmanbaby 2026-08-09 26/1300 2026-08-11 17:16 by 天神眷顾
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 yPv6ndDtiRxQ 2026-08-11 3/150 2026-08-11 17:15 by CLKTSgjM0Uqj
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +6 布布和一二 2026-08-10 7/350 2026-08-11 16:57 by Ai19911123
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +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 +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +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
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
信息提示
请填处理意见