[实用参考]Vasp经验总结.ppt
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
vasp计算总结
wk.baidu.com
加U（LMCO为例）
LDAU=.TRUE. LDAUTYPE=2 #The simplified (rotationally invariant) approach to the LSDA+U
#specifies the l-quantum number for which the on-site LDAUL= -1 2 2 -1 interaction is added (-1=no on-site terms added, 1= p, 2= d, 3=
POSCAR.sh chmod +x POSCAR.sh 改为可执行文件 ./POSCAR.sh 运行
VASP中，用Berry Phase的方法计算极化值
INCAR中要添加的参数：
LCALCPOL=.TRUE. 打开计算Berry的开关
EDIFF=1E-6
精度
DIPOL=0.4 0.4 0.4 选取参考点，任意选的，但是不要和离子重合
磁矩的方向(见vasp手册，有两种设置方法)
NBANDS=408
能带数是线性计算的二倍
ISYM=0
计算soc最好去掉对称性
GGA_COMPAT=.FALSE.
把线性计算的 IBZKPT文件产生的K点copy成KPONITS文件，作为soc计 算的输入文件
YFTO 492 32 300 第6行的数据表示：第一个数表示体系总的价电子数目，第二个数
表示的计算能带时总的k点数目，第三个数表示的是计算能带时计
算了多少条能带。
0.3742002E-15 0.1262379E-14 0.3742002E-15 0.1388889E-01
1 -46.511682 2 -46.511280 3 -46.451176
用VASP计算DOS
ISTART=1 ICHARG=11 ISMEAR= -5 LORBIT=11
将上一步自洽计算得到的CHG、CHGCAR拷贝至同一目录下
准备好KPOINTS文件，增加k点网格
用VASP计算能带
简立方的高对 称点
ISTART= 1 ICHARG = 11
NBANDS一般可以用默认，有时候可以适当增加
f, Default: LDAUL=2)
LDAUU=0 5.5 4 0 #specifies the effective on-site Coulomb interaction parameters LDAUJ=0 1.5 1 0 #specifies the effective on-site Exchange interaction parameters.
ISMARE=-5 不可用。 半导体或绝缘体用0，金属用1
能带处理： 程序 band.cpp，g++编译，得到可执行文件a.out vasp计算得到EIGENVAL,删除前七行 grep -v E EIGENVAL > band.in ./a.out 运行，得到bnd000.dat
NAME 20 Line-mode Rec 0.0 0.0 0.0 0.0 0.5 0.0
LMAXMIX=6
计算soc的时候一定要有此参数，d电子4，f电子6
MAGMOM=6*0 0 0 4 0 0 4 0 0 4 0 0 -4 0 0 -4 0 0 -4 0 0 3 0 0 3 0 0 3 0 0 3 0 0 -3
0 0 -3 0 0 -3 0 0 -3 72*0
SAXIS=1 1 1
! echo "$i" ; mpirun -np 10 vasp >display.log E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'` mv display.log display.log-$i mv OUTCAR OUTCAR-$i mv DOSCAR DOSCAR-$i echo $i $E >>comment done
0.0 0.5 0.0 0.5 0.5 0.0
0.5 0.5 0.0 0.5 0.5 0.5
0.5 0.5 0.5 0.0 0.0 0.0
bnd000.dat文件的第一列数据是k点距离的绝对值，第二列数据是以 Ferim level为参考的本征值。
简立方的高对称点
能量本征值文件EIGENVAL
第一行，前三个整数无意义，第四个整数，如果是2, 表
第8行的前三个数是k点的坐标，第四个数是相 应k点的权重。
4 -46.449273 5 -46.199455 6 -46.197355
第9行给出的是该k点对应的本征值的序号(即第几 条能带)，及相应的本征值。
找最低能量的晶格参数C
#!/bin/sh rm WAVECAR for i in 7.6 7.7 7.8 7.9 8.0 8.1 8.2 8.3 8.4 8.5 8.6 （改为要计算的值） do cat > POSCAR <<! (POSCAR文件，要计算的量改为 $i)
极化的，11个中心对称点群是没有极化的。
点电荷估计： 晶格某一方向所有原子坐标*离子价态相加，与优化后相加得到的值对比。
Spin orbital coupling (SOC)的计算 （LMCO为例）
ISTART=1 ICHARG=1
要读取WAVCAR 最好选取1而不用11，对结果影响是比较大的
LSORBIT=.TRUE.
首先，判断一个体系是否有可能存在极化 1.看晶胞所属的群是否是极化群，磁性有影响时要考虑进去。 MS查看对称性，找到所对应的群。
http://homepage.univie.ac.at/nikos.pinotsis/spacegroup.html#14 找到对应的点群
铁电体物理，钟维烈著，P651。附录一，30个晶体点群。 10个极性点群是非常有可能产生极化的，11个非极性中兴对称群是可能有
60 60 1 1
示是自旋极化的计算，如果是1, 表示非自旋极化的计算。
0.1158924E+02 0.5522500E-09 0.2280000E-08 0.5522500E-09 0.5000000E-15
1.000000000000000E-004 CAR
第2至5行的数据含义不大明确，可以不管它。
wk.baidu.com
加U（LMCO为例）
LDAU=.TRUE. LDAUTYPE=2 #The simplified (rotationally invariant) approach to the LSDA+U
#specifies the l-quantum number for which the on-site LDAUL= -1 2 2 -1 interaction is added (-1=no on-site terms added, 1= p, 2= d, 3=
POSCAR.sh chmod +x POSCAR.sh 改为可执行文件 ./POSCAR.sh 运行
VASP中，用Berry Phase的方法计算极化值
INCAR中要添加的参数：
LCALCPOL=.TRUE. 打开计算Berry的开关
EDIFF=1E-6
精度
DIPOL=0.4 0.4 0.4 选取参考点，任意选的，但是不要和离子重合
磁矩的方向(见vasp手册，有两种设置方法)
NBANDS=408
能带数是线性计算的二倍
ISYM=0
计算soc最好去掉对称性
GGA_COMPAT=.FALSE.
把线性计算的 IBZKPT文件产生的K点copy成KPONITS文件，作为soc计 算的输入文件
YFTO 492 32 300 第6行的数据表示：第一个数表示体系总的价电子数目，第二个数
表示的计算能带时总的k点数目，第三个数表示的是计算能带时计
算了多少条能带。
0.3742002E-15 0.1262379E-14 0.3742002E-15 0.1388889E-01
1 -46.511682 2 -46.511280 3 -46.451176
用VASP计算DOS
ISTART=1 ICHARG=11 ISMEAR= -5 LORBIT=11
将上一步自洽计算得到的CHG、CHGCAR拷贝至同一目录下
准备好KPOINTS文件，增加k点网格
用VASP计算能带
简立方的高对 称点
ISTART= 1 ICHARG = 11
NBANDS一般可以用默认，有时候可以适当增加
f, Default: LDAUL=2)
LDAUU=0 5.5 4 0 #specifies the effective on-site Coulomb interaction parameters LDAUJ=0 1.5 1 0 #specifies the effective on-site Exchange interaction parameters.
ISMARE=-5 不可用。 半导体或绝缘体用0，金属用1
能带处理： 程序 band.cpp，g++编译，得到可执行文件a.out vasp计算得到EIGENVAL,删除前七行 grep -v E EIGENVAL > band.in ./a.out 运行，得到bnd000.dat
NAME 20 Line-mode Rec 0.0 0.0 0.0 0.0 0.5 0.0
LMAXMIX=6
计算soc的时候一定要有此参数，d电子4，f电子6
MAGMOM=6*0 0 0 4 0 0 4 0 0 4 0 0 -4 0 0 -4 0 0 -4 0 0 3 0 0 3 0 0 3 0 0 3 0 0 -3
0 0 -3 0 0 -3 0 0 -3 72*0
SAXIS=1 1 1
! echo "$i" ; mpirun -np 10 vasp >display.log E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'` mv display.log display.log-$i mv OUTCAR OUTCAR-$i mv DOSCAR DOSCAR-$i echo $i $E >>comment done
0.0 0.5 0.0 0.5 0.5 0.0
0.5 0.5 0.0 0.5 0.5 0.5
0.5 0.5 0.5 0.0 0.0 0.0
bnd000.dat文件的第一列数据是k点距离的绝对值，第二列数据是以 Ferim level为参考的本征值。
简立方的高对称点
能量本征值文件EIGENVAL
第一行，前三个整数无意义，第四个整数，如果是2, 表
第8行的前三个数是k点的坐标，第四个数是相 应k点的权重。
4 -46.449273 5 -46.199455 6 -46.197355
第9行给出的是该k点对应的本征值的序号(即第几 条能带)，及相应的本征值。
找最低能量的晶格参数C
#!/bin/sh rm WAVECAR for i in 7.6 7.7 7.8 7.9 8.0 8.1 8.2 8.3 8.4 8.5 8.6 （改为要计算的值） do cat > POSCAR <<! (POSCAR文件，要计算的量改为 $i)
极化的，11个中心对称点群是没有极化的。
点电荷估计： 晶格某一方向所有原子坐标*离子价态相加，与优化后相加得到的值对比。
Spin orbital coupling (SOC)的计算 （LMCO为例）
ISTART=1 ICHARG=1
要读取WAVCAR 最好选取1而不用11，对结果影响是比较大的
LSORBIT=.TRUE.
首先，判断一个体系是否有可能存在极化 1.看晶胞所属的群是否是极化群，磁性有影响时要考虑进去。 MS查看对称性，找到所对应的群。
http://homepage.univie.ac.at/nikos.pinotsis/spacegroup.html#14 找到对应的点群
铁电体物理，钟维烈著，P651。附录一，30个晶体点群。 10个极性点群是非常有可能产生极化的，11个非极性中兴对称群是可能有
60 60 1 1
示是自旋极化的计算，如果是1, 表示非自旋极化的计算。
0.1158924E+02 0.5522500E-09 0.2280000E-08 0.5522500E-09 0.5000000E-15
1.000000000000000E-004 CAR
第2至5行的数据含义不大明确，可以不管它。