VASP计算 各类输入文件说明

POSCAR 记录计算体系的结构信息;
Converted by xsd2pos.py        //备注文字
1.0                              //体系放大的比例
    12.29753000     0.00000000     0.00000000        //晶格基本三矢量的信息
     0.00000000    12.77997000     0.00000000
     0.00000000     0.00000000    15.00000000
  C  N  Fe                //原子的种类和数量
  54  4  1
Direct                          //各原子坐标的描述方式,下面是具体的原子坐标
     0.00046600     0.11056900     0.03455800
     0.09908200     0.27765400     0.03500900
     0.99991700     0.22169000     0.03378400
     0.10117600     0.05591200     0.03603000
     0.20128000     0.11235400     0.03757400
     0.29849900     0.28504400     0.03932900
     0.20036900     0.22459600     0.03736700
     0.30144900     0.05662300     0.03906400
     0.40000000     0.11546200     0.03967300
     0.50150100     0.28504400     0.03932900
     0.40000000     0.22974100     0.03996100
     0.49855100     0.05662300     0.03906400
     0.59872000     0.11235400     0.03757400
     0.70091800     0.27765400     0.03500900
     0.59963100     0.22459600     0.03736700
     0.69882400     0.05591200     0.03603000
     0.79953400     0.11056900     0.03455800
     0.90000000     0.27711600     0.03266100
     0.80008300     0.22169000     0.03378400
     0.90000000     0.05542900     0.03428600
     0.99834900     0.44461000     0.03201400
     0.09727700     0.61123000     0.03432800
     0.99834900     0.55539000     0.03201400
     0.09727700     0.38877000     0.03432800
     0.19497800     0.44417700     0.03705100
     0.19497800     0.55582300     0.03705100
     0.60502200     0.44417700     0.03705100
     0.70272300     0.61123000     0.03432800
     0.60502200     0.55582300     0.03705100
     0.70272300     0.38877000     0.03432800
     0.80165100     0.44461000     0.03201400
     0.90000000     0.61177900     0.03161300
     0.80165100     0.55539000     0.03201400
     0.90000000     0.38822100     0.03161300
     0.99991700     0.77831000     0.03378400
     0.10117600     0.94408800     0.03603000
     0.00046600     0.88943100     0.03455800
     0.09908200     0.72234600     0.03500900
     0.20036900     0.77540400     0.03736700
     0.30144900     0.94337700     0.03906400
     0.20128000     0.88764600     0.03757400
     0.29849900     0.71495500     0.03932900
     0.40000000     0.77025900     0.03996100
     0.49855100     0.94337700     0.03906400
     0.40000000     0.88453800     0.03967300
     0.50150100     0.71495500     0.03932900
     0.59963100     0.77540400     0.03736700
     0.69882400     0.94408800     0.03603000
     0.59872000     0.88764600     0.03757400
     0.70091800     0.72234600     0.03500900
     0.80008300     0.77831000     0.03378400
     0.90000000     0.94457100     0.03428600
     0.79953400     0.88943100     0.03455800
     0.90000000     0.72288400     0.03266100
     0.29355900     0.60644400     0.03990100
     0.29355900     0.39355600     0.03990100
     0.50644100     0.60644400     0.03990100
     0.50644100     0.39355600     0.03990100
     0.40000000     0.50000000     0.04213900
POTCAR 记录计算体系所用到的泛函;(文件过大,不宜展示,通常使用VASPKIT直接自动生成,注意赝势库的选用)
KPOINTS 记录用到的布里渊域的密度;(专业术语可能不太准确)
KPT-Resolved Value to Generate K-Mesh: 0.060        //备注
0                                                        //自动生成KPOINTS
Gamma                                                    //使用gamma点为中心
   1   1   1                                    //各晶格矢量下的K点的密度,越高计算结果越精确,但是计算量越大
0.0  0.0  0.0                                        //偏移量
INCAR 包含计算所用到的参数,决定VASP程序计算什么;
对于电子结构优化计算,INCAR的模版文件如下:
self-consistent                               //备注页,不读取

######Start Parameter for This Run
ISTART = 0                                  //决定是否读取WAVECAR来初始化轨道。https://www.vasp.at/wiki/index.php/ISTART
ICHARG = 2  //决定如何构建电子密度https://www.vasp.at/wiki/index.php/ICHARG
NWRITE = 0                                  //决定OUTCAR记录信息的详细程度
NPAR = 2                                    //并行参数

######Electronic Relaxation
ENCUT = 480                                //决定平面波基组的截断能量 https://www.vasp.at/wiki/index.php/ENCUT
NELM = 100                                 //电子迭代的最高步骤
EDIFF = 1E-5                               //电子迭代结束的收敛标准
GGA = PE                                   //期望使用的GGA类型    https://www.vasp.at/wiki/index.php/GGA
ISMEAR = 0                                 //设定每个轨道的密度的方式    https://www.vasp.at/wiki/index.php/ISMEAR
SIGMA = 0.05                               //决定电子结构迭代时费米能级展宽的宽度
PREC = Accurate                            //计算精确度的设定
ISPIN = 2                                  //是否计算自旋极化,1不计算,2计算
IBRION  = -1                                //不优化离子的几何结构

#for out                                    //决定各类输出文件是否被记录
LCHARG = .TRUE.                    
LWAVE = .FALSE.
LVTOT = .FALSE.
LELF = .FALSE.
对于几何结构的优化,INCAR的模版文件如下:
SYSTEM Relaxation    //同上

######Start Parameter for This Run
ISTART = 0          //同上 不再赘述
ICHARG = 2
NWRITE = 0
NPAR = 2

######Electronic Relaxation
ENCUT = 480        //同上 不再赘述
NELM = 120
EDIFF = 1E-5
GGA = PE
ISMEAR = 0
SIGMA = 0.05
PREC = Accurate
ISPIN = 2

######Ion Relaxation
IBRION = 2        //决定离子如何更新移动的
ISIF = 2          //体系几何结构优化的对象,2只优化几何结构,3也优化晶胞
NSW = 100         //几何优化离子的迭代步骤
EDIFFG = -0.01    //几何优化迭代的收敛判定
POTIM = 0.5       //几何优化每步的步长
ALGO = Fast       //一个方便的参数,确定电子迭代的算法
LREAL = A         //决定投影算符在时空间上计算或是在倒易空间中计算

#for out           //输出文件
LCHARG = .FALSE.
LWAVE = .FALSE.
LVTOT = .FALSE.
LELF = .FALSE.


全部评论

相关推荐

头像
11-18 16:08
福州大学 Java
影流之主:干10年不被裁,我就能拿别人一年的钱了,日子有盼头了
点赞 评论 收藏
分享
10-09 22:05
666 C++
找到工作就狠狠玩CSGO:报联合国演讲,报电子烟设计与制造
点赞 评论 收藏
分享
点赞 收藏 评论
分享
牛客网
牛客企业服务