晶胞几何优化计算出错求助

lora2021

大家好:      今天遇到个奇怪的问题，一个小晶胞提交计算，第一遍几何优化，一切正常，收敛很快，精度是1E-4，第二遍，啥也没动，调整精度为1E-5，然后就报错：
VERY BAD NEWS! internal error in subroutine IBZKPT:
Reciprocal lattice and k-lattice belong to different class of lattices. Often results are still useful...     168


网上查了半天，修改了ALGO=48， NCORE =1，之类的，都没用
感到很莫名，啥也没动啊，就调了个精度而已啊？
请教大家，怎么办？附上我的INCAR

#   ALGO = Fast
   EDIFF  = 1E-5
   ISMEAR = -5
   SIGMA = 0.2
   NCORE = 4     ## Important to reduce calculate time, equal sqrt(cpu cores)

#   LDAU=.TRUE.  ## For LDA + U set
#   LDAUTYPE=2
#   LDAUL=2 -1
#   LDAUU=3  0
#   LDAUJ=0  0

  ISYM = 0
#  SYMPREC = 1E-8

#   ISPIN = 1    ## For magnetism set
#  MAGMOM = 3*1
Ionic Relaxation
   ICHARG = 2
   EDIFFG = -0.02        force
   NSW    = 300          max steps for IOM
   IBRION = 2
   POTIM = 0.5
   ISIF   = 3            suface:2 crystal:3
   PSTRESS= 0
   LWAVE  =  T


KPOINTS是： 5 10 5， 改成10 21 10，也一样的报错。

lora2021

Auto
0
G
5 10 5
0.  0.  0.

sobereva

lora2021 发表于 2022-8-4 11:10
Auto
0
G

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

lora2021

问题解决了，是我POSCAR的问题，重算了一遍，结果合理。

lora2021

sobereva 发表于 2022-8-5 10:09
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员 ...

好的，我是新手，一定注意。

huanglei

lora2021 发表于 2022-8-6 09:52
问题解决了，是我POSCAR的问题，重算了一遍，结果合理。

您好，请问POSCAR哪里出问题会导致这个原因呢？

lora2021

huanglei 发表于 2022-8-6 17:57
您好，请问POSCAR哪里出问题会导致这个原因呢？

我也不知道是哪里问题，没找出来。最初出错，我第二次几何优化的时候，是把第一次优化完的CONTCAR直接拷贝成POSCAR计算，估计就是这里出问题了。后来我直接把两次几何优化合成一次，就没问题了。所幸我的体系不大，12个原子，计算速度还是比较快的。