KPOINTS
对于气体分子来说1 1 1足够
1.自己写KPOINTS
Automatic mesh #随便写不识别,但是得写点
0
Gamma #Gamma只写要给G就可以
N1 N2 N3
0. 0. 0.
晶格常数*N = 20~30
N的数值越大越精确
2.用qvasp书写
qvasp -k 0.03
0.03这个数字越校点位越多
POTCAR
1.自己书写POTCAR
以书写H2的POTCAR为例子
cat ~/pot/Al/POTCAR ~/pot/C/POTCAR ~/pot/H/POTCAR >POTCAR
其中~/pot/Al/POTCAR 这个是赝势存放的地址, 注意赝势调取的顺序应该和POSCAR中的原子顺序一样
2.通过qvasp书写
qvasp -pbe pseudo_potential_name
#查看生成情况
grep "PBE" POTCAR
这里再来认识一下各种赝势
这里的赝势种类很多,如何选择呢?手册的195页中的黑体字给出了每个元素的推荐选择
- Fe_GW(用的人少暂且不谈
- Fe_pv表示把p电子也作为价电子
- Fe_sv表示把s电子也作为价电子
- Ge_d表示把d电子也作为价电子
- _h, _s 应该是 hard和soft的缩写。带h的POTCAR中截断能比普通的要高出很多。带s的截断能要小很多。
认识一下赝势中几个常见的参数
- VRHFIN用来看元素的价电子排布
- TITEL 指的是哪个元素,以及POTCAR发布的时间;
- ENMAX 代表默认的截断能。与INCAR中的ENCUT这个参数相关。