在差分电荷计算中，自洽计算的nsw=0，那么该如何判断收敛

Drizzy

如题所示，我有如下几个问题，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

UW_0728.

感觉你这些问题属于是计算化学还没有完全入门，你没有完全弄明白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”这个是结构优化收敛的提示，静态自洽计算不会出现。

Drizzy

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

好的，谢谢老师