计算化学公社

标题: vasp结构优化一直跑，POSCAR原文件改变 [打印本页]

作者
Author: Drizzy 时间: 2026-1-6 13:10
标题: vasp结构优化一直跑，POSCAR原文件改变
这是INCAR文件，由VASPKIT的Lattice Relaxation生成，因为并没有收敛，所以我参照论坛里的方法，修改参数之后就收敛了。我有三个问题，第一个问题是①是哪个参数有问题导致它一直在超算跑，跑了两天还没有停止，直到我手动用scancel命令才停止，会不会是我没加电子步收敛参数EDDIF的原因，还是LREAL没有开Auto导致的，其次我的ENCUT截断能参数是按照POTCAR里ENMAX的1.3倍设置的，这个参数合理吗？第二个问题是我的②有没有哪些参数设置的不合理。第三个问题是我计算的是ZnCO3晶胞，且我固定了底部的原子，为什么在①中我的POSCAR原文件被改变了，而收敛的②中是正常的，请老师们解答。
POTCAR和KPOINT均由VASPKIT生成，①用的是Gamma点，②用的是M点。相关文件已上传

①
Global Parameters
ISTART =  0          (Read existing wavefunction, if there)
ISPIN  =  1          (Non-Spin polarised DFT)
# ICHARG =  11       (Non-self-consistent: GGA/LDA band structures)
LREAL  = .FALSE.    (Projection operators: automatic)
ENCUT  =  520       (Cut-off energy for plane wave basis set, in eV)
# PREC =  Accurate (Precision level: Normal or Accurate, set Accurate when perform structure lattice relaxation calculation)
LWAVE  = .FALSE.       (Write WAVECAR or not)
LCHARG = .FALSE.       (Write CHGCAR or not)
ADDGRID= .TRUE.       (Increase grid, helps GGA convergence)
LASPH  = .TRUE.       (Give more accurate total energies and band structure calculations)
PREC = Accurate    (Accurate strictly avoids any aliasing or wrap around errors)
# LVTOT  = .TRUE.    (Write total electrostatic potential into LOCPOT or not)
# LVHAR  = .TRUE.    (Write ionic + Hartree electrostatic potential into LOCPOT or not)
# NELECT =          (No. of electrons: charged cells, be careful)
# LPLANE = .TRUE.    (Real space distribution, supercells)
# NWRITE = 2          (Medium-level output)
# KPAR = 2          (Divides k-grid into separate groups)
# NGXF = 300       (FFT grid mesh density for nice charge/potential plots)
# NGYF = 300       (FFT grid mesh density for nice charge/potential plots)
# NGZF = 300       (FFT grid mesh density for nice charge/potential plots)

Lattice Relaxation
NSW =  300       (number of ionic steps)
ISMEAR =  0          (gaussian smearing method )
SIGMA  =  0.05       (please check the width of the smearing)
IBRION =  2          (Algorithm: 0-MD, 1-Quasi-New, 2-CG)
ISIF =  3          (optimize atomic coordinates and lattice parameters)
EDIFFG = -2E-02    (Ionic convergence, eV/A)

K-POINTS
Auto
0
G
2 2 1
0.  0.  0.

②
Global Parameters
ISTART =  0          (Read existing wavefunction, if there)
ISPIN  =  1          (Non-Spin polarised DFT)
# ICHARG =  11       (Non-self-consistent: GGA/LDA band structures)
LREAL  = Auto    (Projection operators: automatic)
ENCUT  =  520       (Cut-off energy for plane wave basis set, in eV)
PREC =  Accurate (Precision level: Normal or Accurate, set Accurate when perform structure lattice relaxation calculation)
LWAVE  = .FALSE.       (Write WAVECAR or not)
LCHARG = .FALSE.       (Write CHGCAR or not)
ADDGRID= .TRUE.       (Increase grid, helps GGA convergence)
LASPH  = .TRUE.       (Give more accurate total energies and band structure calculations)
PREC = Accurate    (Accurate strictly avoids any aliasing or wrap around errors)
# LVTOT  = .TRUE.    (Write total electrostatic potential into LOCPOT or not)
# LVHAR  = .TRUE.    (Write ionic + Hartree electrostatic potential into LOCPOT or not)
# NELECT =          (No. of electrons: charged cells, be careful)
# LPLANE = .TRUE.    (Real space distribution, supercells)
# NWRITE = 2          (Medium-level output)
# KPAR = 2          (Divides k-grid into separate groups)
# NGXF = 300       (FFT grid mesh density for nice charge/potential plots)
# NGYF = 300       (FFT grid mesh density for nice charge/potential plots)
# NGZF = 300       (FFT grid mesh density for nice charge/potential plots)

Lattice Relaxation
NSW =  1000       (number of ionic steps)
ISMEAR =  0          (gaussian smearing method )
SIGMA  =  0.05       (please check the width of the smearing)
IBRION =  2          (Algorithm: 0-MD, 1-Quasi-New, 2-CG)
ISIF =  3          (optimize atomic coordinates and lattice parameters)
EDIFF  =  1E-06
POTIM  = 0.05
EDIFFG =  -0.02       (Ionic convergence, eV/A)

K-POINTS
K-Spacing Value to Generate K-Mesh: 0.020
0
Monkhorst-Pack
7  11 1
0.0  0.0  0.0

作者
Author: shinnashiori 时间: 2026-1-6 22:11
本人也在学习中，如有不对也请各位指出，谢谢
1. 应该就是你没设置电子步收敛的问题，截断能设置应该是没问担心改大点儿再跑跑
2. 应该是没有啥不合理的地方吧
3. 你说的 POSCAR 改变是指啥意思？是原子坐标改变还是说你固定的层数没固定住？

作者
Author: UW_0728. 时间: 2026-1-6 23:11
1. 看OSZICAR文件大致判断。EDIFF不设置的话默认为10E-4；LREAL用哪个更好与晶胞大小有关，但绝对不至于造成这个问题；ENCUT设置我认为合理，不放心的话可以做做收敛性测试
2. NSW太大，设小一点比如300免得振荡得太厉害了还在一直跑；k点数目有什么理由变动那么明显？
3. 不理解你说的“POSCAR改变”具体是什么意思。一般来说POSCAR作为晶体结构输入文件应该自始至终没有变化，后续的结构是在CONTCAR（最后一步结构）和XDATCAR（优化或动力学的轨迹）里的，除非你自己中途手动修改了。

作者
Author: Drizzy 时间: 2026-1-7 21:18

UW_0728. 发表于 2026-1-6 23:11
1. 看OSZICAR文件大致判断。EDIFF不设置的话默认为10E-4；LREAL用哪个更好与晶胞大小有关，但绝对不至于造 ...

老师您好，谢谢您的解答，第三个问题是我没看清楚文件。我又看了下我的OSZICAR文件，在即将到达第十五个离子步的时候，它就没有往后算了，我觉得应该和这个有关系。其次我想问一下ENCUT能不能直接根据ENMAX设置，还是说必须要进行收敛性测试；再者，如果是对于菱锌矿，我用不同的面去切，它的ENMAX会不会不同，比如我用101面和104面去切。
  N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.897682525044E+03 0.17203E+01 -0.22555E+03  5376 0.636E+01 0.169E+01
DAV: 2 -0.934429529914E+03 -0.36747E+02 -0.10097E+02  7344 0.176E+01 0.514E+01
DAV: 3 -0.899899381734E+03 0.34530E+02 -0.90171E+01  8304 0.147E+01 0.108E+01
DAV: 4 -0.900037397480E+03 -0.13802E+00 -0.58981E+00  6384 0.313E+00 0.674E+00
DAV: 5 -0.899926198766E+03 0.11120E+00 -0.13345E+00  8208 0.173E+00 0.413E+00
DAV: 6 -0.899986300521E+03 -0.60102E-01 -0.21937E+00  7632 0.190E+00 0.413E+00
DAV: 7 -0.899807148736E+03 0.17915E+00 -0.94680E-01  7776 0.181E+00 0.143E+00
DAV: 8 -0.899818058377E+03 -0.10910E-01 -0.15099E-01  6336 0.517E-01 0.101E+00
DAV: 9 -0.899813392553E+03 0.46658E-02 -0.59008E-02  7488 0.397E-01 0.504E-01
DAV:  10 -0.899814770507E+03 -0.13780E-02 -0.26989E-02  6912 0.196E-01 0.512E-01
DAV:  11 -0.899813247873E+03 0.15226E-02 -0.11145E-02  7152 0.149E-01 0.122E-01

作者
Author: Drizzy 时间: 2026-1-7 21:20

shinnashiori 发表于 2026-1-6 22:11
本人也在学习中，如有不对也请各位指出，谢谢
1. 应该就是你没设置电子步收敛的问题，截断能设置应该是没 ...

老师您好，感谢您的解答，对于①楼下的老师说EDIFF默认设置为10E-04，③是我自己的原始文件有问题，再次感谢您的解答

作者
Author: UW_0728. 时间: 2026-1-7 22:27

Drizzy 发表于 2026-1-7 21:18
老师您好，谢谢您的解答，第三个问题是我没看清楚文件。我又看了下我的OSZICAR文件，在即将到达第十五个 ...

ENCUT的设定原则上不得低于相应POTCAR中包含的所有元素的ENMAX的最大值，经验上讲设置为1.3倍较为稳妥；收敛性测试是更严谨的做法，但说不上必需。
ENMAX只取决于元素和使用的赝势/PAW方法，与具体体系毫无关系

作者
Author: Drizzy 时间: 2026-1-8 20:30

UW_0728. 发表于 2026-1-7 22:27
ENCUT的设定原则上不得低于相应POTCAR中包含的所有元素的ENMAX的最大值，经验上讲设置为1.3倍较为稳妥； ...

好的，感谢老师的解答

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