VASP原包的计算过程及原理

1、by Tan SuVASP计算流程概括：首先在Materiale Studio(MS)里导入原包模型，即纯金属模型、氧化物模型、原子模型等把原包坐标转换成VASP坐标，用FTP转入所要计算的服务器（a）然后在SecureCRT中进入服务器a，建立文件夹（mkdir 文件名，如sut）, 也可以直接用ftp键入好文件夹，在所建好的原包文件件，如Au中键入四个指标（vi 名称，名称必需大写！）四个指标分别是：INCAR，KPOINTS，POSCAR，POTCAR复制VASP脚本，即：cp VASP脚本 交作业： qsub 目录（如所在即是要算的目录，可不加） VASP脚本 作业运行：qstat一般

2、会进入队列排队（R），也可能出错，请检查OUTCAR文件多余作业取消：qdel 作业代号23原包的计算点击左键4保存为*.cif和*.car文件根据所需金属的不同，建立好文件夹5通过Gauss程序输入*.cif格式，保存转换成*.gjf格式另一转换方法见表面部分讲解6用FTP进入所要计算的服务器：ljy01 ： 06ljy03a：02ljy07a：02ljy01b-ljy04b：53pdr06b: 53ljy01c：04ljy01e-ljy02e: 192.9.2

3、07.240ljy01g：8ljy01h-ljy02h：9st01d: 06建立属于自己的文件夹：建立属于自己的文件夹：sut(sut(姓的全拼和名拼音首字母姓的全拼和名拼音首字母) )在上一文件夹中建立新文件夹，以所要计算的原包命名：在上一文件夹中建立新文件夹，以所要计算的原包命名：AuAu在上一文件夹中建立新文件夹，以所计算的点或者面，原包在上一文件夹中建立新文件夹，以所计算的点或者面，原包选择：选择：0 0或或000, 000, 同理选择所计算的面命名，如：同理选择所计算的面命名，如：110110在上一文件夹中建立四个参数

4、：在上一文件夹中建立四个参数： INCAR INCAR，KPOINTSKPOINTS，POSCARPOSCAR，POTCARPOTCAR（也可以用第（也可以用第2 2页的方法添加，稍后说明）页的方法添加，稍后说明）并添加当前服务器所用并添加当前服务器所用VASPVASP程序，即左图程序，即左图vasp.pbsvasp.pbs。注意：每个服务器的注意：每个服务器的VASPVASP版本不同，一定要确认好！版本不同，一定要确认好！78首次进入按照上图设置，下次进入会出现左下图快捷单，选择即可通过”cd 文件夹”命令读取所需文件夹，进入错了输入”cd .”命令退回上一级，”ls”命令检查文件夹所包含的

5、子目录9原包的计算除了可以在FTP中把已有四个参数放入相应文件夹外，还可以在SecureCRT面板中进行：如下图所示，进入你相应的文件夹在0文件夹中输入：“vi INCAR”（参数） 注意：大小写！字母”I”键进入编辑输入左边黑色字符，作出相应修改；Esc键退出修改模式输入”：wq”保存修改，没有修改直接退出键入“：q”;修改错了取消修改，强制退出编辑状态键入”：q!”此外复制键入“cp”, ctrl+C是不管用的！“shift+G”键直接跳到段尾SYSTEM =AuENCUT=400PREC =MediumEDIFF = 1E-05EDIFFG = -0.02ISPIN =2NELM = 1

6、00NELMIN =6ISYM =2VOSKOWN = 1LREAL = AISMEAR =1SIGMA = 0.2ISIF =3IALGO = 8IBRION =2POTIM = 0.1NSW =300LWAVE = .FALSE.LCHARG = .FALSE.计算体系名称.常规为400，有时用300.根据精度选择，一般第一遍用低精度: low.为收敛标准.一般用-0.02-0.03.1不自旋 2自旋（原包选2）.为步数.对称性.选择实/倒空间计算.电子占据数.选择维数：2表面，3原包. 8 原包的计算38 特定快分解法（Kosugi算法）48 预处理后残量最小化法1过渡态 2稳态（原包）

7、 5频率步长：低精度0.5，高精度0.1圈数：稳态300，过渡态1000是否输出波函数是否输出电荷密度10原包的计算在0文件夹中输入：“vi KPOINTS”（K点） 注意大写！ 输入下列字符行：Automatic mesh0M11 11 110. 0. 0.M角度为90时，即方形；G菱形一般用高精度的即可如果过没有KPOINTS文件，在INCAR中加入两个参数，则自动生成K点11原包的计算在0文件夹中输入：“vi POSCAR”（结构坐标） 输入下列字符行：Au 1.000000000000000 4.07830000 0.00000000 0.00000000 0.00000000 4.0

8、7830000 0.00000000 0.00000000 0.00000000 4.07830000 4 原子数 C 笛卡尔坐标（D相对坐标） 0.00000000 0.00000000 0.00000000 0.00000000 2.03915000 2.03915000 2.03915000 0.00000000 2.03915000 2.03915000 2.03915000 0.00000000 基矢缩放系数基矢缩放系数12原包的计算“POTCAR”（贋势）不需要输入，从文件夹中选取即可，一般常用的是paw_gga方法中的贋势 贋势就是把离子实的内部势能用假想的势能取代真实的势能，但

9、在求解波动方程时，不改变能量本征值和离子实之间区域的波函数。由赝势求出的波函数叫赝波函数，在离子实之间的区域真实的势和赝势给出同样的波函数。13交作业如果VASP版本没有在该目录下，需要写路径在交完作业后，一定要键入”qstat”, 检查是否成功运行运行完查看OUTCAR文档查看是否正确键入“qdel 作业号”删除多余的作业，以免浪费机时作业号14表面的计算表面的计算和原包的计算步骤大致相同，只是在原包计算的基础上多出了开始的切出不同的表面，以Au（1 1 1）面为例：确认确认MS中晶胞长度中晶胞长度a，b，c及角度与原包的及角度与原包的*.OUTCAR中相同中相同坐标调为坐标调为Origin

10、al模式（见原包部分）模式（见原包部分）15表面的计算所选表面所选表面原子数为原子数为五个为准，五个为准，即五层即五层先切出所要研究面，厚先切出所要研究面，厚度选择为五层；再建立度选择为五层；再建立超包，超包，2*2或者或者3*316表面的计算建立真空层，一般固气相要选择15以上即够了。液固相需要根据两者晶胞长度，要比晶胞长些。存成*.cif和*.car文件17表面的计算在VESTA软件中打开所要转化的*.cif文件，如右图所示，转化为*.VASP文件另外一种方式见过渡态计算部分18表面的计算固定最上两层，即最大数值的两个，所以区域为7-9.5之间即可自动生成四个参数，需自动生成四个参数，需要检查参数是否正确要检查参数是否正确19表面的计算元素元素原子数原子数选择性动力学选择性动力学笛卡尔坐标笛卡尔坐标20表面的计算M方形；G菱形算算表面表面时先用低精度时先用低