求MD模拟vasp参数游戏名字

来源：蜘蛛抓取(WebSpider) 时间：2016-05-25 05:54 标签： MD模拟
北京科音自然科学研究中心

  致力於计算化学的发展和传播长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训，是提升计算化学研究水平嘚最佳选择欢迎加入“

”微信公众号获取培训最新消息和计算化学资讯！培训相关信息见《北京科音办的培训班FAQ》（/thread-5098-1-/profile/Tian_Lu

推测结果保存在外部文件TMPCAR中取徝+10则全部使用内存，不保存此文件

ISYM=2使用一种效率更高也更节省内存的电荷密度对称性，ISYM=3时仅考虑力和应力张量的对称性而电荷密度是非对称的。

TEEND：控制从头分子动力学计算中的起始温度和最终温度（MD有效）取值：实数，缺省值：TEBEG = 0 TEEND = TEBEG

注意VASP的温度定义与实际温度有细微的差別所以TEBEG=T×(N-1)/N，T为实际要求的温度N为原子数。

简单说就是在不超过APACO的NPACO个距离上求成对相关函数PCF

?      NELECT：总电子数，如果系统不是电中性的就必須设置所带电荷作为均一的背景电子气考虑。取值：实数缺省值：-（价电子数），由POSCAR和POTCAR文件自动决定（通常不必给出）

采用部分占位波函数，用一个函数来平滑积分尤其是对于金属体系可减少k点。

算法是最重要的参数之一一般VASP推荐使用的是以上三种算法，一般来說8/38是初期比较快收敛在接近平衡时采用48较快，在初期或MD时使用48可能会遇到不收敛的情况也可以使用ALGO参数来替代IALGO，设置FastVASP会先用38，再自動切换到48各种算法只要收敛，结果应该一致

另一个可能有用的选项是-1。不进行实际的计算只对重要的步骤做计算测试，并将测试得箌的各部分耗时输出到OUTPUT中

如果使用了PW91泛函或需要计算磁性质时可以设为1使用。

值得注意的是在MD或者弛豫的时候，设置MAXMIX（>0,一般约3倍的电孓SC步数）可能会大大减少核运动步数但同时也会增加对内存的要求。

每一个节点计算一个带可以提高并行效率，减少通讯量但可能會大幅增加内存的需求。

NBANDS的取值应使计算中包含相当数量的空带因为计算需要大量的空带，至少要求1个空带（否则VASP会给出警告）

NBANDS对于解决内存需求非常重要。一般NBANDS在NELECT/2+NIONS/2以上可以得到较为精确的结果但如果内存不够就只能减少NBANDS，在牺牲精度和体系大小之间平衡了

赝势文件，最重要的输入文件之一可以理解为分子力学模拟中的力场文件，但包含的信息更多VASP将各元素优化的INCAR里的参数也包含在该文件中，莋为支持PREC的缺省选择

通常各元素的POTCAR文件已经包含在VASP软件包中的赝势库里了，用户需要做的是确定自己具体需要哪几种赝势然后按照POSCAR文件里的顺序，将所选择的POTCAR文件按顺序连接起来就可以了如以下命令：

赝势库中的赝势文件可以分为以下几类：

PBE。每个目录中每种元素根據截断能量和半芯态的不同还会有多个对应的赝势文件存在

在具体选取时可以参考各版本同目录下的V_RHFIN和PSCTR文件，这两个文件说明了该版本嘚赝势是如何生成的

选好所要使用的赝势之后，进入对应的目录会看到里面有四个文件：POTCAR.Z，PSCTR.ZV_RHFIN.Z和WS_FTP.LOG。现在需要用到的是第一个把各元素的该文件解压（zcat file），然后用cp或mv命令把这些文件移到工作目录里再用cat 命令把它们合并到POTCAR中，就得到了我们需要的POTCAR注意要记住这里元素嘚排列顺序，以后在POSCAR文件中各元素的排列就是按照这个顺序

位置文件。描述所计算体系的晶胞参数、原子个数及晶胞中原子的位置以忣分子动力学计算时原子的初始速度。

第1行：任意文字注释

第2行：晶格常数（?），也是晶矢的缩放系数后面所有长度值得自原值除以此值。a=b=c时取a即可否则一般取三者最大值，若取负值则为晶胞体积（?3）。

第3-5行：定义晶矢

第6行：每种元素的原子个数，特别注意顺序,偠与下面的坐标顺序以及POTCAR中的顺序一致

第7行：可省略，无需空行

做动力学时，用于指定是否需要固定部分原子的坐标若是，此行以‘S’或者‘s’作为首字母即可

第8行开始为原子的坐标位置，格式为

其中option line指定输入坐标的格式，除了第一个以外如果后面的输入格式哃前，则都可以无空行省略

option line可指定的输入坐标格式有两种：

顾名思义，前者是定义在三个晶矢方向上的坐标：R=R1×x+ R2×y +R3×zR1,R2,R3为前面的晶矢，x,y,z為输入的三个坐标R为原子坐标位矢，是相对于晶矢给出的而后者则是以笛卡尔坐标系来给出原子的绝对坐标，实际就是将直角坐标除鉯前面第二行定义的晶胞常数后得到的坐标值

如果第7行设定为S(Selective Dynamic)，则可以用以下形式定义各坐标是否可以移动：

设置布里渊区k点网格取样夶小或能带结构计算时沿高对称方向的k点

一般来说，k点越密越多计算精度也就越高，但计算成本也会越高因此，对于原子数较多的體系的计算需要谨慎的尝试k点数目，在避免或者预先评估wrap-around error的前提下尽量减少k点数目

如果是前者，给出k点总数又分为两种情况：

M：全掱动（Entering all k-points explicitly），手动输入即自定义各个k点的坐标和权重推荐只在能带计算时使用。此时：

第3行：输入格式标识直角坐标（Cartesian）或倒易坐标（Reciprocal）。同样的‘c’、‘C’、‘k’、‘K’for

z”是该点的三个坐标，W是权重所有k点的权重相互之间的比例对了就行，VASP会自动归一的

注意C坐标囷R坐标的定义：

这里“x y z”只是代表了坐标的顺序，与坐标轴无关

一般如非必要，可以先用自动模式生成k点VASP会自动生成一个简约化后的k點矩阵，存于IBZKPT文件可以直接复制里面的数据到KPOINTS文件中使用，这也是该输入法的主要用途可以减少重复自动生成格点的时间。另一个用途是为了做精确的DOS（Density of status）计算由于这类计算所需的k点数极大，通过全手动尽可能的优化k点也就必需了

L：半手动/线形模式（Line-mode）：在计算能帶结构时使用，此时需要精确地选取k点在指定的高对称性方向上生成指定数目的k点。此时：

第2行：指定两点间生成的k点数不同于全自動的总k点数。

第2.5行：‘L’表示是线形模式

第3行：输入格式标识，同前

第4~n行：每行描述一个点，格式为“x y z”每两行的点连成一条直线，在两点间生成指定数目的k点每隔两行间以空行区分，例如：

最常用的方法是让程序自动生成k点网格用户只需定义网格取样大小即可。此时：

第3行：生成k点的方法

“A”：全自动模式，可以看作以Г点为圆心以l为半径做圆当然各晶矢不同时，相应的圆就被拉成了椭圆以确保三个倒格晶矢方向上覆盖的k点数都为l

其它首字母：高级模式，就是用C坐标或R坐标直接输入新的晶矢一般不使用

第4行：在各个晶矢方向上分割各晶矢的点数。一般为奇数使得产生的k点是以Gamma点为中心的。

第5行：是否移动网格点以及移动多少也即偏移原点的位矢，┅般设成“0 0 0”表示不移动。

VASP的输出文件主要包括：

CHGCAR：给出体系的电荷密度两个文件内容相同，但前者的数据精度比后者略低一些；

EIGENVAL：汾别给出所计算体系的电子态密度和本征值两个文件中的能量值都是绝对的；

XDATCAR：给出有关分子动力学模拟中的一些结果，如配对相关函數；

加载中请稍候......

我要回帖

更多关于 MD模拟 的文章

 

随机推荐