计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1569|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] vasp计算能带报错,出现plane wave coefficients changed怎么办?

[复制链接 Copy URL]

240

帖子

0

威望

1005

eV
积分
1245

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
WAVECAR: different cutoff or change in lattice found
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...
reading WAVECAR
ERROR: while reading WAVECAR, plane wave coefficients changed       26059
       25966

先结构优化,在自洽,最后非自洽,非自洽的 INCAR 内容如下:Global Parameters
ISTART =  1            (Read existing wavefunction; if there)
# ISPIN =  2           (Spin polarised DFT)
ICHARG =  11         (Non-self-consistent: GGA/LDA band structures)
LREAL  = .FALSE.       (Projection operators: automatic)
# ENCUT  =  400        (Cut-off energy for plane wave basis set, in eV)
PREC   =  Normal       (Precision level)
LWAVE  = .TRUE.        (Write WAVECAR or not)
LCHARG = .TRUE.        (Write CHGCAR or not)
ADDGRID= .TRUE.        (Increase grid; helps GGA convergence)
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)
# NPAR   = 4           (Max is no. nodes; don't set for hybrids)
# NWRITE = 2           (Medium-level output)
# KPAR   = 2           (Divides k-grid into separate groups)
# NGX    = 500         (FFT grid mesh density for nice charge/potential plots)
# NGY    = 500         (FFT grid mesh density for nice charge/potential plots)
# NGZ    = 500         (FFT grid mesh density for nice charge/potential plots)
ICORELEVEL=1
NBANDS=368

Static Calculation
ISMEAR =  0            (gaussian smearing method)
SIGMA  =  0.01         (please check the width of the smearing)
LORBIT =  11           (PAW radii for projected DOS)
NEDOS  =  2001         (DOSCAR points)
NELM   =  60           (Max electronic SCF steps)
EDIFF  =  1E-08        (SCF energy convergence; in eV)



IBRION=-1
ISIF=2
NSW=0
EDIFFG=-0.05
NGX =    48
NGY =   48
NGZ =   90
# ISMEAR=-5
ALGO=Fast
LORBIT=10



1060

帖子

0

威望

3256

eV
积分
4316

Level 6 (一方通行)

2#
发表于 Post on 2024-5-15 17:35:18 | 只看该作者 Only view this author
你前面的计算里指定了 NGX, NGY, NGZ 三个参数没有?  如果没有设置的话,现在也不需要设置。

  1. IBRION=-1
  2. ISIF=2
  3. NSW=0
  4. EDIFFG=-0.05
  5. NGX =    48
  6. NGY =   48
  7. NGZ =   90
  8. # ISMEAR=-5
  9. ALGO=Fast
  10. LORBIT=10
复制代码

240

帖子

0

威望

1005

eV
积分
1245

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-5-16 09:03:36 | 只看该作者 Only view this author
谢谢大神

240

帖子

0

威望

1005

eV
积分
1245

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2024-5-16 09:03:50 | 只看该作者 Only view this author
乐平 发表于 2024-5-15 17:35
你前面的计算里指定了 NGX, NGY, NGZ 三个参数没有?  如果没有设置的话,现在也不需要设置。

2

帖子

0

威望

89

eV
积分
91

Level 2 能力者

5#
发表于 Post on 2024-9-23 22:12:17 | 只看该作者 Only view this author

麻烦问一下,plane wave coefficients changed这个报错应该怎么考虑去解决

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-23 20:55 , Processed in 0.206718 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list