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 选取参考点，任意选的，但是不要和离子重合
BP计算极化过程中，Dipole设置问题
设置在计算离子的dipole时的参考点即设置 DIPOL（注意的是，它的 设置需要使得原子移动前后的原子都在这个参考点的一侧。比如这个例 子中 Al处于(0,0,0)，As处于(0.25, 0.25, 0.25)位置，而将DIPOL设置为( 0.5, 0.5, 0.5)和(0.125, 0.125, 0.125)都是可以的，但是在考虑移动Al原子时，不要将
NELM = 40 # maximum of 40 electronic steps
杂化泛函的计算 --HF Functional
GGA=PS(选用的赝势文件为PBBEsol, =PE为PBE的方法) LHFCALC = .True. PRECFOCK = Normal # NKRED = 2 (设置此参数容易报错，不知为何) TIME = 0.4 HFSCREEN = 0.2 AEXX = 0.25 #the exact exchange is used
2、ALGO, IALGO, LDIAG If the self-consistency loop does not converge within 40 steps, it will probably not converge at all. In this case you should reconsider the tags IALGO, LDIAG, and the mixing-parameters. 一般情况下，或使用IALGO=48时遇到收敛问题的话，可以考虑设IALGO为38， 或设置ALGO=Normal or Fast (in VAS P.4.5 and later versions)。 Default ALGO = Normal 3、NELMDL NELMDL gives the number of non-selfconsistent steps at the beginning In some cases (for instance MD’s, or ionic relaxation) you might set NELMIN to a larger value (4 to 8)
用VASP计算DOS
ISTA-5 LORBIT=11
将上一步自洽计算得到的CHG、CHGCAR拷贝至同一目录下 准备好KPOINTS文件，增加k点网格
用VASP计算能带
ISTART= 1 ICHARG = 11
简立方的高对 称点 NAME 20 Line-mode Rec 0.0 0.0 0.0 0.0 0.5 0.0 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
首先，判断一个体系是否有可能存在极化 1.看晶胞所属的群是否是极化群，磁性有影响时要考虑进去。 MS查看对称性，找到所对应的群。 http://homepage.univie.ac.at/nikos.pinotsis/spacegroup.html#14 找到对应的点群 铁电体物理，钟维烈著，P651。附录一，30个晶体点群。 10个极性点群是非常有可能产生极化的，11个非极性中兴对称群是可能有 极化的，11个中心对称点群是没有极化的。
此段从别处复制的，侵删
加真空-用MS
1.Build-surface-Cleave surface (0 0 1)指沿c方向加真空层，要 是在其他方向上加对应修改一 下就可 点击Cleave就行
2.Build-Crystal-Build Vacuum Slab Vacuum thickness就是真空层的厚度 点击 Build 即可
点电荷估计： 晶格某一方向所有原子坐标*离子价态相加，与优化后相加得到的值对比。
Spin orbital coupling (SOC)的计算 （LMCO为例）
ISTART=1 ICHARG=1 要读取WAVCAR 最好选取1而不用11，对结果影响是比较大的
LSORBIT=.TRUE. 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 磁矩的方向(见vasp手册，有两种设置方法) NBANDS=408 能带数是线性计算的二倍 ISYM=0 计算soc最好去掉对称性 GGA_COMPAT=.FALSE.
原子移动原胞之外即偏移量为负数；另外也不要将DIPOL设置在所要移动的原子上，
如果是这样 的话，则会导致移动该原子后，该原子不在DIPOL的同一侧，使得原子 移动之后的极化的Berry-phase项要比平衡态时的大很多。） 要注意的就是DIPOL 的设置，以及原子移动的选取：不要将DIPOL设置在某个原子 位置上，原子移动的选取要保证原子在移动前后都是DIPOL的同一侧
最后可以将.cif的文件放到vesta里面。 导出.vasp的文件，就能放到vasp中运算 了
Partial charge density（PCD）计算-VASP
/blog/static/2094090822012892478779/ 定义： Partial charge density计算或称为Band decomposed charge density计算，即计算 特定的某个(或某些)k点和本征值(这些k点和本征值是相互对应的)所对应的本征波函 数的平方(也就是电荷密度)。 第一种Partial Charge分析的INCAR ISTART = 1 ICHARG = 1 LPARD=.TRUE. IBAND= 20 21 22 23 KPUSE= 1 2 3 4 LSEPB=.TRUE. LSEPK=.TRUE.
vasp计算总结
加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= 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.
bnd000.dat文件的第一列数据是k点距离的绝对值，第二列数据是以 Ferim level为参考的本征值。
简立方的高对称点
能量本征值文件EIGENVAL
第一行，前三个整数无意义，第四个整数，如果是2, 表 示是自旋极化的计算，如果是1, 表示非自旋极化的计算。
60 60 1 1 0.1158924E+02 0.5522500E-09 0.2280000E-08 0.5522500E-09 0.5000000E-15 1.000000000000000E-004 第2至5行的数据含义不大明确，可以不管它。 CAR YFTO 492 32 300 第6行的数据表示：第一个数表示体系总的价电子数目，第二个数 表示的计算能带时总的k点数目，第三个数表示的是计算能带时计 算了多少条能带。 0.3742002E-15 0.1262379E-14 0.3742002E-15 0.1388889E-01 1 -46.511682 第8行的前三个数是k点的坐标，第四个数是相 2 -46.511280 应k点的权重。 3 -46.451176 4 -46.449273 第9行给出的是该k点对应的本征值的序号(即第几 5 -46.199455 条能带)，及相应的本征值。 6 -46.197355
NPAR=16(使用核数为4,16,64此种) ALGO=ALL
Berry Phase计算极化，老方法
一般来说，现在都用新方法来计算极化，因为用起来方便。 老方法的话还要把三个方向分来来计算，等于要跑三个任务，相对复杂，但是
准确性可能好一些。
LBERRY=.TRUE. IGPAR=3 #1,2,3对应x, y, z方向 NPPSTR=7 DIPOL=0.5 0.5 0.5 EDIFF=1E-5
找最低能量的晶格参数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)
! 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