vasp攻略

打包压缩命令：tar zcvf 文件名.tar.gz 源文件名
采用link方式避免重复的文件浪费内存：ln –s 源文件命名
Eg：ln –s ../optic/MME ./EuO.mme
交互式绘图工具gnuplot: 命令行打:gnuplot进入格式:plot “文件名”退出:quit
求磁矩:getmag
如果由于节点掉线在提交任务后秒退,不输出outcar可以指定节点提交任务
LJRS -l nodes=c0104:ppn=4
chmod +x 名称——使脚本可以执行
构造potcar，以A和B元素为例：
如果是以Z为拓展名的文件：zcat A/POTCAR.Z B/POTCAR.Z > POTCAR
如果是解压后的potcar文件：cat A/POTCAR B/POTCAR > POTCAR
INCAR中的RWIGS通过POTCAR文件获得
单位晶胞体积：grep “vol”OUTCAR
自动计算加应力情况下的最优化情况：
C/a ：vaspcaopt softmode-e -=*(画曲线）
如果要看某种材料是FM还是AFM，需要以相同结构计算一次，看能量哪个低。

将POSCAR/CONTCAR/CHGCAR装换成xcrysden(进入xcrysden文件执行./xcrysden)可读取的形式Eg：（BFO)v2xsf CHGCAR -1 83 -2 26 -3 8 -d
v2xsf POSCAR -1 83 -2 26 -3 8 -d
v2xsf CONTCAR -1 83 -2 26 -3 8 -d
晶胞放大时，k点需要等比例的缩小
画出曲线图：
ISMEAR如果是半导体/绝缘体取-5，如果是金属取1，SIGMA=0.2(一般不改变)
如果求DOS，则ISMEAR=-5.如果求band，则ISMEAR=1.
1、Relax ISTART=0，ICHARG=2, ISIF=3，NSW=200，EDIFFG=-1*10-3，IBRION=2
2、Scf NSW=0(关闭结构优化) 用Relax后的CONTCAR替换POSCAR
3、DOS 添加scf后的CHGCAR，ISTART=1，ICHARG=11，k放大一倍
DOSCAR 第六行：Emax Emin Emin与Emax之间点的数目Ef
第七行：能量总的态密度（spin up）总的态密度（spin down）态密度积分（up）态密度积分（down）
后面按原子分：能量s轨道态密度（spin up）s轨道态密度（spin down）p轨道态密度（spin up）p轨道态密度（spin down）d轨道态密度（spin up）d轨道态密度（spin
down）f轨道态密度（spin up）f轨道态密度（spin down）
Split_dos 对dos按原子区分
4、band 与DOS所需文件一致，KPOINTS需要使用line模式手动输入。

(getbandk进行修正）高
对称点的坐标通过文献确定。

一定要正确判断结构。

EIGENVAL 第六行：总的价电子数（nelect）k点的个数能带的条数（一般来说
nbands=nelect/2+nions/2，要考虑附加更多空带nbands=nelect/2+nions/2+nions*2）
Split_band 对能带按自旋方向区分
查看总能grep 'TOTEN' OUTCAR
5、加应力情况（省略relax步骤）
加应力会导致dos图向低能方向移动，费米面处能隙逐渐消失，使得半导体向金属转变
ISIF=2。

如果是四方情况，修改POSCAR，给a=b加不大于6%的应变（，
是稳定结构下的a），调整c直到获得使得总能最低的最优化c。

6、Berry Phase情况
计算par（顺电）情况：首先做自洽（LBERRY=.F.）NSW=0
LBERRY=.T. 如果是非正交体系，需要计算三次（IGPAR=1,2,3）
NPPSTR可以取6-8，一般取8，如果是一个较大的体系，则取6
DIPOL与所取原胞的原点不重合即可，一般取0.15 0.15 0.15
ISTART=0，ICHARG=11，ISMEAR=0
计算ferr（铁电）情况：与顺电情况类似，但需要对POSCAR进行修改（参见vasp手册88页）
POSCAR中的原子位移设定为0.05 eg:EuO的poscar中Eu原子位置
0.05 0.05 -0.05 T T T ;O原子位置0.0 0.0 0.0 F F F
作弛豫其中NSW=200,ISIF=2,进行结构优化.
如果弛豫后的CONTCAR中原子的位移几乎是0,就表示在该方向是不存在极化的.(EuO+5%应变
情况认为极化沿z轴方向设置Eu原子位置0.05 0.05 -0.05,relax后几乎是0,表示z方向不存在极化)
之后与顺电情况一致,进行自洽后求BF
计算结果分析（Total electronic dipole moment&ionic dipole moment）：以NaF为例
铁电情况下：P1[elc]+ P12[elc]+ P3[elc]+P[ion]——（-10.64900，-10.64899，-10.60335）顺电情况下：P1[elc]+ P12[elc]+ P3[elc]+P[ion]= P[ion]（P1[elc]=P12[elc]=P3[elc]=0）——（-10.649，-10.649，-10.649）
将铁电算出的结果减去顺电算出的结果就是总的偶极矩。

（此时的单位是e*A）——（0,0,0.04565）
进行单位换算——*16.02/单位晶胞的体积（grep “vol”OUTCAR，)得到最终结果。

（此时的单位是C/m2）——（0,0,0.03）
如果要计算模：e*R/V 其中R是极化轴的长度(POSCAR)，计算过程中可以直接采用16.02*R/单位晶胞的体积来进行计算，计算结果的单位仍然是C/m2 ——R=C=4.63*0.5,模是1.5
7、计算声子谱
INCAR: 关闭SIGMA,关闭NSW,关闭ISIF,关闭NPAR,关闭IALGO
ISYM=1,IBRION=6,NFREE=2,EDIFF=1.D-6
OUTCAR中倒数几行出现f/i.如果f/i那行的第四个数字接近0表示是零频,否则表示存在铁电软模即负频。

如果dx有值，则极化方向是x方向;如果dy有值，则极化方向是y方向; 如果dz有值，
则极化方向是z方向.
8、计算lamda
1 将铁电相{做自洽那个文件夹中的}POSCAR重命名为POSCAR-fer，顺电相的重名为POSCAR.ini
2 用脚本getxyzdiff得到铁电相和顺电相的坐标差值g etxyzdiff POSCAR-fer POSCAR.ini
3 得到xyzdiff.out并重命名为XYZDIFF.OUT
4 在job.pbs中加入vasplamdaBP -0.2 1.2 0.1 $NPROCS $LJRS_NODEFILE 4
在“mpirun -v -machinefile $LJRS_NODEFILE -np $NPROCS ~/vasp”前加#注释掉
5 将前面自洽和BP的INCAR复制进来，分别命名为INCAR-scf和INCAR-Berry
6 备齐“INCAR-Berry INCAR-scf POSCAR-fer POSCAR.ini POTCAR KPOINTS ZDIFF.OUT job.pbs”共8个文件就可以提交任务了。

对应于IGPAR=1，2，3各需要计算一次。

设置铁电相POSCAR时选定的原子位移不要出现0.9x这样的要设置成-0.0x以确保XYZDIFF.OUT 中不出现过大的值,一般XYZDIFF.OUT中的值都是正负0.0x
使用"softmode-e? -=*",输出文件softmode-e.dat得到能量E-序参量关系，看到在顺电相时能量最高，铁电相时能量最低。

用之前的方法获得极化强度P-序参量关系(getpolar-new-total OUTCAR，方便数据收集)
如果得到的计划强度P-序参量关系有跳变,通过软模进行修正.软模:feR/V,其中R是极化轴的长度.对于自旋极化体系,f=1;对于非自旋极化体系,f=2.
9、计算光学性质
第一步：
算自洽scf得到CHGCAR，INCAR（ISTART=1，ICHARGE=11，LOPTICS=.T.，NBANDS的设置看自洽OUTCAR中NBANDS的值，这里取它的两倍，要有足够的导带以供跃迁），POSCAR，POTCAR，KPOINTS(k点放大一倍，即与dos的KPOINTS一致)，job.pbs中采用vaspNLO而非vasp5。

得到MME 和STRUC
第二步：
需要四个文件——MME重命名---元素.mme, STRUC重命名---元素.strucin，文件opticpack.def 文件内只有一行是元素名，文件---元素.opticin（eg：EuO.mme EuO.opticin EuO.strucin opticpack.def）
对于文件.opticin：如果第一行选择1线性dielectricze第二行选择0或者1，如果第一行选择2SHG 则第二行选择0，如果第一行选择3SHG II则第二行选择1；第四行的第一个参数是费米能级，通过第一步中grep ‘fermi’OUTCAR 得到，第三个参数是自旋，如果是磁性体系为2,如果是非磁性体系为0；第五行开始设置的时候是NELECT/2（对于磁性体系,通过第一步中grep
‘NELECT’ OUTCAR得到NELECT,如果NELECT是23，就设置成12 11;对于非磁性体系,设置成NELECT/2,如果NELECT是32,就设置成16即可）；第六行是能量范围，其中第五个参数要与第一步中INCAR的NBANDS一致；第七行开启“剪刀“修正，如果无需修正第八行设置为0.
四个文件备齐后执行“opticpack-FM $”(非磁性系统则是opticpack $) .opticout文件，spin1中能级少于费米能级的个数与spin2中能级少于费米能级的个数就是第五行中需要的设置。

(对于非磁性体系与磁性体系,OPTICIN文件有两行是对换的,需要注意)。