vasp计算改变结构的对称性
case: graphite
POSCAR
graphite
2.45600000000000
1.0027447811380699 -0.0000000000000000 0.0000000000000000
-0.5013723905690349 0.8684024501626882 -0.0000000000000000
0.0000000000000000 -0.0000000000000000 2.6108453362728841
C
4
Direct
0.0000000000000000 -0.0000000000000000 0.0000000000000000
0.3333333347959524 0.6666666695919048 -0.0000000000000000
-0.0000000000000000 0.0000000000000000 0.5000000000000000
0.6666666673979762 0.3333333347959524 0.5000000000000000
做完静态计算后,我用vasp 转 Boltztrap的接口文件生成 Boltztrap的输入文件,如case.struct。
1. 从POSCAR里提取结构的对称信息,得到的结果如下
case
4.65391 -0.00000 0.00000
-2.32695 4.03040 -0.00000
0.00000 -0.00000 12.11737
24
(这个部分是24个对称操作的矩阵形式,太多我就不列举了,确实有24个之多)
2. 从OUTCAR里提取结构的对称信息,得到的结果如下
case
4.65391 -0.00000 0.00000
-2.32695 4.03040 -0.00000
0.00000 -0.00000 12.11737
1
1.0 0.0 0.0
0.0 1.0 0.0
0.0 0.0 1.0
(为什么只有一个对称操作呢,有两种可能。其一,OUTCAR里把结构转换成了P1,所以输出只有一个对称操作。其二,我在转换脚本里写了这么一句,如果找不到OUTCAR里的结构对称信息,就输出P1)
我更倾向从OUTCAR里提取结构的对称信息,感觉能带结构本征值的写法跟这个有关。一个很明显的区别就是,用两种方式提取出来的结构输入文件对boltztrap的计算结果有很大的影响。
问题:1. OUTCAR里这种改变结构对称性的结果,还能用于boltztrap的计算吗?
2. 有没有什么参数可以让OUTCAR里的对称性和POSCAR里的一样?
语言组织的太直白,大家担待一下 返回小木虫查看更多
谢谢~
如何“从OUTCAR里提取结构的对称信息”?
对于graphite,空间群为P 63 m c,对称操作应该是12个
从POSCAR的数据来看,VASP判断对称性的精度是够的(SYMPREC默认10^-5)
所以我猜测你计算的INCAR中设置了ISYM=0?
一般计算软件会认为graphite是No.194,即D6h对称性,有24个对称操作
上面的POSCAR文件,如果用VASP计算,OUTCAR给出的是24个对称操作
Subroutine GETGRP returns: Found 24 space group operations
你可以检查OUTCAR(同时还体现在IBZKPT上)
至于你的脚本从OUTCAR中提取对称操作,我想不出是什么操作的。估计正如你所说的,如果没有合理判断就置1吧,