计算化学公社

标题: 晶胞几何优化计算出错求助 [打印本页]

作者
Author: lora2021 时间: 2022-8-4 11:03
标题: 晶胞几何优化计算出错求助
大家好:    今天遇到个奇怪的问题，一个小晶胞提交计算，第一遍几何优化，一切正常，收敛很快，精度是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，也一样的报错。

作者
Author: lora2021 时间: 2022-8-4 11:10
Auto
0
G
5 10 5
0. 0. 0.

作者
Author: sobereva 时间: 2022-8-5 10:09

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

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

作者
Author: lora2021 时间: 2022-8-6 09:52
问题解决了，是我POSCAR的问题，重算了一遍，结果合理。

作者
Author: lora2021 时间: 2022-8-6 09:54

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

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

作者
Author: huanglei 时间: 2022-8-6 17:57

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

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

作者
Author: lora2021 时间: 2022-8-7 09:45

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

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

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