计算化学公社

标题: 在差分电荷计算中，自洽计算的nsw=0，那么该如何判断收敛 [打印本页]

作者
Author: Drizzy 时间: yesterday 14:24
标题: 在差分电荷计算中，自洽计算的nsw=0，那么该如何判断收敛
如题所示，我有如下几个问题，1.为什么nsw=0，在OSZICAR中还跑了20个电子步和一个离子步，这是我非常不解的，其次如果IBRION=-1固定了原子，那么此时是不是无法使用能量和力来判断收敛，该如何判断收敛，亦或者说只需要得到CHGCAR文件就能满足需求。
2.在②的结构优化中，关于能量的收敛是看dE和d eps还是d E,我的EDIFF统一设置为1E-6,如果是看d E,那为什么即使我的d E没有达到这个值，在OUTCAR中依然显示了reached required accuracy，但是在vaspview中，我的力在最后达到了收敛，这是不是说只需要力收敛就行了。
3.关于vaspview和qvasp，接着第二个问题，如果在vaspview中我的能量差达到了一个平衡趋势，但是没有达到我的设定值，是否也会判断为收敛；使用qvasp -e作为判断收敛的方法可不可信，在问题1中，我的scf自洽计算在OUTCAR中没有显示reached required accuracy字样，但是使用qvasp -e其会显示get converged :YES!
恳请各位老师解答，谢谢
① N    E                   dE          d eps    ncg    rms       rms(c)DAV: 1    0.265335398225E+04 0.26534E+04 -0.11251E+05 12160 0.166E+03
DAV: 2    0.138468769993E+03 -0.25149E+04 -0.23711E+04 13584 0.430E+02
DAV: 3 -0.326252082037E+03 -0.46472E+03 -0.45756E+03 13432 0.198E+02
DAV: 4 -0.353808320272E+03 -0.27556E+02 -0.27259E+02 14232 0.536E+01
DAV: 5 -0.354596908309E+03 -0.78859E+00 -0.78597E+00 15600 0.894E+00 0.363E+01
RMM: 6 -0.323428990432E+03 0.31168E+02 -0.80218E+01 14310 0.224E+01 0.208E+01
RMM: 7 -0.321531452066E+03 0.18975E+01 -0.26096E+01 14665 0.144E+01 0.883E+00
RMM: 8 -0.321209962810E+03 0.32149E+00 -0.53326E+00 15083 0.703E+00 0.643E+00
RMM: 9 -0.321072155280E+03 0.13781E+00 -0.18081E+00 14463 0.469E+00 0.146E+00
RMM:  10 -0.321053949960E+03 0.18205E-01 -0.55313E-01 15006 0.196E+00 0.792E-01
RMM:  11 -0.321063186894E+03 -0.92369E-02 -0.13058E-01 15050 0.114E+00 0.413E-01
RMM:  12 -0.321065248860E+03 -0.20620E-02 -0.39982E-02 15302 0.530E-01 0.167E-01
RMM:  13 -0.321068682860E+03 -0.34340E-02 -0.14140E-02 15228 0.295E-01 0.172E-01
RMM:  14 -0.321068942137E+03 -0.25928E-03 -0.47307E-03 14515 0.195E-01 0.469E-02
RMM:  15 -0.321069025842E+03 -0.83705E-04 -0.14924E-03 17048 0.962E-02 0.465E-02
RMM:  16 -0.321068861300E+03 0.16454E-03 -0.40173E-04 14854 0.489E-02 0.307E-02
RMM:  17 -0.321068825395E+03 0.35905E-04 -0.14834E-04 14669 0.322E-02 0.136E-02
RMM:  18 -0.321068815745E+03 0.96494E-05 -0.53332E-05 15317 0.179E-02 0.871E-03
RMM:  19 -0.321068814257E+03 0.14887E-05 -0.14132E-05 13674 0.905E-03 0.443E-03
RMM:  20 -0.321068813356E+03 0.90053E-06 -0.43199E-06 11051 0.516E-03
1 F= -.32518044E+03 E0= -.32518044E+03  d E =-.381524E-12
②.  N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.321065474727E+03 0.37187E-04 -0.50349E-02  2480 0.614E-01 0.489E-02
RMM: 2 -0.321065517617E+03 -0.42891E-04 -0.62989E-04  2451 0.745E-02 0.304E-02
RMM: 3 -0.321065516790E+03 0.82745E-06 -0.10400E-04  2452 0.305E-02 0.122E-02
RMM: 4 -0.321065517763E+03 -0.97314E-06 -0.25295E-05  2464 0.143E-02 0.458E-03
RMM: 5 -0.321065517922E+03 -0.15859E-06 -0.42737E-06  1912 0.695E-03
372 F= -.32517637E+03 E0= -.32517778E+03  d E =-.203398E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.321065664127E+03 -0.14636E-03 -0.23144E-02  2544 0.424E-01 0.344E-02
RMM: 2 -0.321065690288E+03 -0.26161E-04 -0.36054E-04  2432 0.568E-02 0.215E-02
RMM: 3 -0.321065689724E+03 0.56436E-06 -0.56213E-05  2439 0.236E-02 0.896E-03
RMM: 4 -0.321065690169E+03 -0.44551E-06 -0.14933E-05  2362 0.113E-02 0.336E-03
RMM: 5 -0.321065690248E+03 -0.79173E-07 -0.25799E-06  1567 0.558E-03
373 F= -.32517651E+03 E0= -.32517792E+03  d E =-.141241E-03

作者
Author: UW_0728. 时间: yesterday 16:23
本帖最后由 UW_0728. 于 2026-1-22 16:59 编辑

感觉你这些问题属于是计算化学还没有完全入门，你没有完全弄明白INCAR里必要的计算参数以及计算输出或OSZICAR中的东西各自代表什么含义，有空一定要多到VASP Wiki上勤搜索。
1. NSW=0等于是直接对当前结构做静态自洽（算单点能），显然需要进行SCF迭代以收敛到稳定的波函数；另外VASP wiki写得明明白白：应当避免IBRION=-1搭配NSW>0以免完全相同的结构跑n多次无意义的SCF迭代，且NSW=0时IBRION默认为-1（见https://www.vasp.at/wiki/IBRION）。此时没有离子步收敛标准一说，而电子步收敛标准（即SCF收敛标准）就是EDIFF所指定的能量变化值，对应需要检测的是电子步迭代过程中的“dE”和“d eps”两个值（假如EDIFF设为1.0E-6，那就是要求这两个值都达到1.0E-6以下），而不是电子步完成后最下面输出里面的“d E”。你的输出1中已经SCF正常收敛（达到了EDIFF的要求），没别的问题的话可以直接拿得到的CHGCAR用就行了。
2. 结构优化的收敛标准看EDIFFG，不看EDIFF（SCF收敛标准），一般设为EDIFFG = -0.02（加负号表示看受力而非能量变化；0.02表示最大原子受力小于0.02eV/angstrom，这也是结构优化任务一般推荐的收敛限）；如果编译VASP时集成了VTST模块，那么通过在任务所在目录运行“grep FORCES: OUTCAR”可以很方便地监控受力变化（包括最大受力和均方根受力）。所以你说得没错，只需要力收敛就行了。
3. 其一，同2，电子步收敛标准和离子步收敛标准不是看的一个东西；其二，“reached required accuracy”这个是结构优化收敛的提示，静态自洽计算不会出现。

作者
Author: Drizzy 时间: yesterday 17:41

UW_0728. 发表于 2026-1-22 16:23
感觉你这些问题属于是计算化学还没有完全入门，你没有完全弄明白INCAR里必要的计算参数以及计算输出或OSZIC ...

好的，谢谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)