VASP的输入输出文件课件.ppt

VASPVASP的输入输出文件的输入输出文件侯侯 柱柱 锋锋复旦大学物理系复旦大学物理系2005, 12, 3北京，宏剑Hou Zhufeng2内容(1)输入文件输入文件INCAR :计算控制参数文件计算控制参数文件POSCAR :描述体系结构的文件描述体系结构的文件KPOINTS :k点取样设置文件点取样设置文件POTCAR :赝势文件赝势文件输出文件输出文件OUTCAR :最主要的输出文件最主要的输出文件DOSCAR : 电子态密度文件电子态密度文件EIGENVAL :本征值文件本征值文件OSZICAR :每次迭代或原子迟豫每次迭代或原子迟豫(或或MD)的信息的信息CHG和和CHGCAR :电荷密度文件电荷密度文件WAVECAR :波函数文件波函数文件3内容(2)输出文件输出文件(续续)CONTCAR :原子迟豫或原子迟豫或MD后的体系结构文件后的体系结构文件IBZKPT :布里渊区中的布里渊区中的k点点PCDAT :对关联函数对关联函数XDATCAR :在在MD时时, 原子位置变化的跟踪文件原子位置变化的跟踪文件PROCAR和和PROOUT : 波函数投影或分解的文件波函数投影或分解的文件LOCPOT :总的局域势总的局域势ELFCAR :电子局域函数电子局域函数4输入文件输入文件5INCAR(1)设置计算的精度设置是做什么计算设置交换关联函数设置优化的算法和收敛标准设置MD的步长、温度、时间设置每个波函数或电子的占有数(也就是 smearing方法及相关的参数)等等自由格式共100多个关键词，一般都有合适的默认值。

通常需设置的有10个左右即使INCAR文件无任何内容也可以进行计算6INCAR(2)：例子静态计算静态计算(计算总能和自洽的电荷密度计算总能和自洽的电荷密度)懒：删掉计算目录下的懒：删掉计算目录下的WAVECAR，INCAR中不作任何设置中不作任何设置推荐：自己手动设置推荐：自己手动设置ENCUT, TITLE, ISTART, ICHARG, ISMEAR, EDIFF, PREC等等SYSTEM = SiliconENCUT = 350ISTART = 0ICHARG = 2ISMEAR = -5EDIFF = 1E-5PREC = Accurate设置标题，以说明所计算的体系设置标题，以说明所计算的体系设置平面波切断动能设置平面波切断动能(不采用默认值不采用默认值)说明这次计算是一次全新的计算说明这次计算是一次全新的计算按体系中的原子构造初始的原子密度按体系中的原子构造初始的原子密度采用四面体方法采用四面体方法电子迭代的收敛标准是电子迭代的收敛标准是1E-5精度为精度为Accurate可以在一行设置多个关键词的值，但是每个关键值之间用分号可以在一行设置多个关键词的值，但是每个关键值之间用分号(;)隔开。

如隔开如ISMEAR= 0; SIGMA= 0.2当想不用当想不用INCAR中某个关键词的值时，在该行前面加上井号中某个关键词的值时，在该行前面加上井号(#)注释掉，注释掉，如如#ISMEAR=0; SIGMA = 0.27POSCAR(1)描述了体系的结构文件: 原(或晶)胞的基矢, 原子的位置, 原子是否移动，原子的初始速度等)分数坐标和卡笛尔坐标Cubic SiC3.570.0 0.5 0.50.5 0.0 0.50.5 0.5 0.0 1 1Direct 0.00 0.00 0.000.25 0.25 0.25例子例子:Cubic SiC3.570.0 0.5 0.50.5 0.0 0.50.5 0.5 0.0 1 1Cartesian 0.00 0.00 0.000.25 0.25 0.25设置体系的名称设置体系的名称晶格常数或缩放系数晶格常数或缩放系数原原(或晶或晶)胞的基矢胞的基矢每类原子的个数每类原子的个数确定按何种坐标来写原子位置确定按何种坐标来写原子位置第一类原子的第一个坐标第一类原子的第一个坐标第二类原子的第一坐标第二类原子的第一坐标第七行以字母第七行以字母D开头表示下面的是分数坐标，如果是以开头表示下面的是分数坐标，如果是以C或或K开头表示下面的坐标开头表示下面的坐标是卡笛尔坐标。

是卡笛尔坐标8POSCAR(2)对原子位置进行选择性的迟豫(Selective dynamic)Cubic SiC3.570.0 0.5 0.50.5 0.0 0.50.5 0.5 0.0 1 1Selective dynamicDirect 0.000 0.000 0.000 F F F0.248 0.248 0.248 T T T例子例子:表示对体系中的原子进行限制性的迟豫表示对体系中的原子进行限制性的迟豫或优化，这一行以字母或优化，这一行以字母S开头如何写出具有复杂结构的晶体的如何写出具有复杂结构的晶体的POSCAR：a) 查到该晶体的晶格常数、空间群查到该晶体的晶格常数、空间群和乌科夫和乌科夫(Wyckoff)坐标坐标; b)用用Material Studio中中Crystal builder或其他的晶体学软或其他的晶体学软件画出晶体，并得到各个原子的坐标件画出晶体，并得到各个原子的坐标;或或c)根据空间群国际表写出各等价位置根据空间群国际表写出各等价位置; d)晶晶体数据库体数据库固定第一类原子的位置，固定第一类原子的位置，对第二类原子的位置进行优化对第二类原子的位置进行优化9KPOINTS(1)设置布里渊区k点取样或k点的坐标,计算能带时k点的高对称线设置方式:手动输入所有的k点，按Monkhorst-Pack方法自动产生格式：普通，四面体，M-P网格，Line模式卡笛尔坐标系和倒格子坐标系10KPOINTS(2)按普通格式手动输入各个按普通格式手动输入各个k点点(推荐用在计算能带时推荐用在计算能带时)k-points entered manually5Reciprocal 0.500 0.000 0.500 1.00 0.475 0.000 0.475 1.00 0.450 0.000 0.450 1.00 0.425 0.000 0.425 1.00 0.400 0.000 0.400 1.00标题或注释行，无特别意义标题或注释行，无特别意义K点的数目点的数目以字母以字母R开头表示开头表示k点是按倒格子坐标系点是按倒格子坐标系前三个数是前三个数是k点的坐标，最后一个数是相应点的坐标，最后一个数是相应k点的权重点的权重(下面共下面共5个个k点点)如果是以卡笛尔坐标系来写k点坐标，则第三行以字母C开头。

按四面体格式手动输入各个按四面体格式手动输入各个k点点(复杂，不推荐用复杂，不推荐用)Example file 4 Cartesian 0.0 0.0 0.0 1. 0.0 0.0 0.5 1. 0.0 0.5 0.5 2. 0.5 0.5 0.5 4. Tetrahedra 1 0.183333333333333 6 1 2 3 4 标题或注释行，无特别意义K点的数目以字母C开头表示k点是按卡笛尔坐标系前三个数是k点的坐标，最后一个数是相应k点的坐标(下面共4个k点)当smearing方面采用的是四面体方法时，以字母T开头表示下面的是描述各个四面体的信息四面体的个数，每个四面体的体积权重每个四面体顶角的权重，每个四面体顶角k点坐标的序号当当Smearing时不采用四面体方法，则时不采用四面体方法，则Tetrahedra后面的内容不起作用后面的内容不起作用11KPOINTS(2)按按Line模式手动输入各个高对称模式手动输入各个高对称k点点(用在计算能带中，用在计算能带中，4.6以上版本支持以上版本支持)k along high symmetry lines 10 Line-mode Rec 0.00 0.00 0.00 ! gamma 0.50 0.50 0.00 ! X0.50 0.50 0.00 ! X 0.50 0.75 0.25 ! W0.50 0.75 0.25 ! W 0.00 0.00 0.00 ! gamma 标题或注释行，无特别意义标题或注释行，无特别意义每对高对称点之间产生每对高对称点之间产生10个个k点点以字母以字母L开头表示按开头表示按line模式产生模式产生k点点以字母以字母R开头表示开头表示k点按倒格子坐标系点按倒格子坐标系每对高对称点的坐标。

共三对，表明沿着每对高对称点的坐标共三对，表明沿着Gamma X W Gamma计算能带）计算能带）如果是以卡笛尔坐标系来写k点坐标，则第四行以字母C开头写出各种布拉菲格子对应的倒格子中的高对称点标记和坐标：写出各种布拉菲格子对应的倒格子中的高对称点标记和坐标：a)查空间群查空间群国际表国际表; b)用用xcrysden软件画出布里渊区，并选取高对称点软件画出布里渊区，并选取高对称点12KPOINTS(3)按按Monkhorst-Pack网格方法自动产生k点Automatic mesh0Auto10 标题或注释行，无特别意义标题或注释行，无特别意义0表示自动产生表示自动产生k点点以字母以字母A开头表示完全自动以开头表示完全自动以点为中心产生点为中心产生k点点长度，间接确定了沿倒格子各方向上网格点的数目长度，间接确定了沿倒格子各方向上网格点的数目 完全自动的方式 (1)Automatic mesh0Gamma1010 10 0.0 0.0 0.0 完全自动的方式 (2)标题或注释行，无特别意义标题或注释行，无特别意义0表示自动产生表示自动产生k点点以字母以字母G开头表示按开头表示按M-P方法以方法以点为中心点为中心自动产生自动产生k点点确定沿倒格子各方向上网格点的数目确定沿倒格子各方向上网格点的数目对所按网格分割产生的对所按网格分割产生的k点进行平移的量（这里不平移）点进行平移的量（这里不平移）第三行也可以字母M开头，则不以点为中心13KPOINTS(4)按按Monkhorst-Pack网格方法自动产生k点 完全自动的方式 (2) (续)Automatic mesh0Monkhorst-Pack1010 10 0.5 0.5 0.5 完全自动的方式 (3)标题或注释行，无特别意义标题或注释行，无特别意义0表示自动产生表示自动产生k点点以字母以字母M开头表示按开头表示按M-P方法自动产生方法自动产生k点点确定沿倒格子各方向上网格点的数目确定沿倒格子各方向上网格点的数目对所按网格分割产生的对所按网格分割产生的k点进行平移的量（这里不平移点进行平移的量（这里不平移)如果第三行以字母G开头，则以点为中心产生k点14KPOINTS(5) 对六角晶系的结构，强烈推荐采用以对六角晶系的结构，强烈推荐采用以为中心按为中心按M-P网格产生网格产生k点点 K点网格大小点网格大小N1, N2, N3按倒格子基矢的长度进行合适的选取按倒格子基矢的长度进行合适的选取一般选取一般选取N1, N2和和N3为奇数，以便产生的为奇数，以便产生的k点包含了点包含了点点 设置不同的设置不同的N1,N2和和N3，看它们对总能的影响，以确定能保证一定，看它们对总能的影响，以确定能保证一定计算精度的计算精度的k点数目点数目#/bin/shfor i in 3 5 7 9 11docat KPOINTS SUMMARYdone15POTCAR按按POSCAR文件中，每类原子的顺序把各类原子相应的赝文件中，每类原子的顺序把各类原子相应的赝势合并到同一个文件势合并到同一个文件(POTCAR)中。

比如在中比如在POSCAR中原中原子坐标按子坐标按Si, O的顺序写出给原子的坐标，则到赝势库目录的顺序写出给原子的坐标，则到赝势库目录中使用下面的命令中使用下面的命令: zcat Si/POTCAR.Z POTCAR zcat O/POTCAR.Z POTCAR每类原子的赝势类型每类原子的赝势类型(PAW或超软或超软)一致一致每类原子的赝势类型每类原子的赝势类型(交换关联交换关联)要与要与INCAR中交换关联设置中交换关联设置一致一致1。