计算六方晶体结构的弹性常数一直出现K点报错

yangxunscience

请问在使用VASP, IBRION=6 计算六方晶体的结构的弹性常数一直出现如下报错如何解决，谢谢大家！

晶体结构使用原胞计算，结构正确。不能关闭ISYM否则计算速度慢而且计算出的弹性常数不满足六方晶体结构的对称性。


internal ERROR in RE_READ_KPOINTS_RD: the new k-point set for the       |
reduced symmetry case                                                   |
does not contain all original k-points. Try to switch off symmetry     |
                                                                           |
---->  I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <----



POSCAR :内容

CONTCAR
   1.0
           7.9028000832         0.0000000000         0.0000000000
           4.9780868173         6.1378254125         0.0000000000
           4.9780868173         2.3720903699         5.6609264322
      Te   Cr   Si
       6    2    2
   Direct
       0.721279979         0.448960006         0.080269992
       0.278719991         0.551039994         0.919730008
       0.080269992         0.721279979         0.448960006
       0.919730008         0.278719991         0.551039994
       0.448960006         0.080269992         0.721279979
       0.551039994         0.919730008         0.278719991
       0.665589988         0.665589988         0.665589988
       0.334410012         0.334410012         0.334410012
       0.944400012         0.944400012         0.944400012
       0.055599999         0.055599999         0.055599999


INCAR： 内容

SYSTEM=CrSiTe3
   ENCUT=600
   ISTART=0
   ICHARG=2
   ISMEAR=0
   SIGMA=0.1
   ALGO=Fast
   GGA=PE
  #NELM=1000
  EDIFF=1E-7
  EDIFFG=-0.001
  POTIM=0.1
  NSW=1000
  PREC=Accurate
  ISYM=2
  ISIF=3
  PREC=Accurate
  IBRION=6
  LCHARG=.F.
  LWAVE=.F.

yangxunscience

使用VASPKIT下基于应变和应变能201得到的六方弹性常数也不满足六方晶体的对称性。仔细检查了结构，没有问题。

vaspkit

贴出vaspkit寻找的输出信息。

yangxunscience

vaspkit 发表于 2021-11-15 10:25
贴出vaspkit寻找的输出信息。

+-------------------------- Summary ----------------------------+
Crystal Class: -3
Space Group: R-3
Crystal System: Trigonal/Rhombohedral II system
Including Point group classes: 3, -3.
There are 7 independent elastic constants
    C11  C12  C13  C14 -C15    0
    C12  C11  C13 -C14  C15    0
    C13  C13  C33    0    0    0
    C14 -C14    0  C44    0 -C15
    C15 -C15    0    0  C44  C14
      0    0    0 -C15  C14  C66

Stiffness Tensor C_ij (in GPa):
     58.621     16.459     20.296     -0.326    -12.553      0.000
     16.459     58.621     20.296      0.326     12.553      0.000
     20.296     20.296     86.873      0.000      0.000      0.000
     -0.326      0.326      0.000     17.428      0.000     12.553
    -12.553     12.553      0.000      0.000     17.428     -0.326
      0.000      0.000      0.000     12.553     -0.326     21.081


根据文献六方晶体的独立弹性常数有5个，C11 C33 C44 C12 C13。(PRB,76,054115,2007)
VASPKIT给出的独立弹性常数7个。
不知道问题出在那个地方? 非常感谢！

vaspkit

yangxunscience 发表于 2021-11-17 20:44
+-------------------------- Summary ----------------------------+
Crystal Class: -3
Space Group:  ...

Trigonal/Rhombohedral II也可以是六角晶格。当然也可以在~/.vaspkit文件中降低寻找对称性精度SYMPREC  1E-4试试独立弹性常数是否有变化。

yangxunscience

vaspkit 发表于 2021-11-17 22:35
Trigonal/Rhombohedral II也可以是六角晶格。当然也可以在~/.vaspkit文件中降低寻找对称性精度SYMPREC  1 ...

好的，谢谢！

yangxunscience

R3有7各独立的晶格常数。

fjy

您好，请问您的问题解决了吗，我也遇到同样的问题，
internal ERROR in RE_READ_KPOINTS_RD: the new k-point set for the reduced symme
try case
    does not contain all original k-points. Try to switch off symmetry
，没有---->  I REFUSE TO CONTINUE WITH THIS SICK JOB ... BYE!!! <----
，这是我的INCAR
SYSTEM = AB
ISTART = 0
ICHARG = 2
ENCUT = 450
EDIFF = 1E-5
EDIFFG = -1E-2
ISMEAR = 1
SIGMA = 0.2
NSW = 200
ISIF = 3
IBRION = 6
NFREE = 4
SYMPREC = 1E-3
这是优化后的POSCAR
AB                                   
   1.00000000000000     
     6.5138484608820528   -3.7607721623506798   -0.0000000000000000
    -0.0000000000000001    7.5215443247013614    0.0000000000000000
    -0.0000000000000000    0.0000000000000000   10.1728737622478409
   A   B
     6    18
Direct
  0.0000000000000000  0.0000000000000000  0.2725901154470381
  0.6666666666666643  0.3333333333333357  0.6059234487803737
  0.3333333333333357  0.6666666666666643  0.9392567821137023
  0.0000000000000000  0.0000000000000000  0.7274098845529620
  0.6666666666666643  0.3333333333333357  0.0607432178862976
  0.3333333333333357  0.6666666666666643  0.3940765512196262
  0.5184988578880706  0.4815011421119293  0.2798341422480040
  0.1851655245547349  0.8148344754452651  0.6131674755813399
  0.8518321912214063  0.1481678087785937  0.9465008089146685
  0.5184988578880706  0.0369977157761413  0.2798341422480040
  0.1851655245547349  0.3703310491094769  0.6131674755813399
  0.8518321912214063  0.7036643824428056  0.9465008089146685
  0.9630022842238587  0.4815011421119293  0.2798341422480040
  0.6296689508905230  0.8148344754452651  0.6131674755813399
  0.2963356175571945  0.1481678087785937  0.9465008089146685
  0.4815011421119293  0.5184988578880706  0.7201658577519958
  0.1481678087785937  0.8518321912214063  0.0534991910853317
  0.8148344754452651  0.1851655245547349  0.3868325244186603
  0.0369977157761413  0.5184988578880706  0.7201658577519958
  0.7036643824428056  0.8518321912214063  0.0534991910853317
  0.3703310491094769  0.1851655245547349  0.3868325244186603
  0.4815011421119293  0.9630022842238587  0.7201658577519958
  0.1481678087785937  0.2963356175571945  0.0534991910853317
  0.8148344754452651  0.6296689508905230  0.3868325244186603

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00

alpha

看看K点文件

alpha

fjy 发表于 2022-1-13 11:10
您好，请问您的问题解决了吗，我也遇到同样的问题，
internal ERROR in RE_READ_KPOINTS_RD: the new k-poi ...

提供你的INCAR、KPOINTS和VASP版本。