《vasp做大体系结构优化(精品文档)》由会员分享，可在线阅读，更多相关《vasp做大体系结构优化(精品文档)（2页珍藏版）》请在金锄头文库上搜索。
1、Vasp对大体系的结构优化ENCUTvasp 运行时间主要消耗在对角化上。运行时间正比于 NBANDS*NPLNW 2,前者是能带数目，后者是平面波数目；由于 NPLNWENCUT 3/2,故运行时间正比于ENCUT 3。运行时间也正比于NELECT 3，电子数目三次方。IALGO选择对角化算法：对于小体系，用 IALGO=38 （ Davidson algorithm ） ；对于大体系，用 IALGO=48 （ RMM-DIIS ）。可以设置 ALGO=very fast or fast RMM-DIIS 并行效率比 Davidson algorithm 高一些。NPAR如果 IALGO=3
2、8 ，则 NPAR 取 1 。对 IALGO=48 ，影响不是特别大，可选2 或 4 ，可选节点数，取值越大，内存占用越多。并行效率总是低于线性叠加效率的，核越多并行效率越低。所以对于一定核（如 20） ，一定作业（如 2个） ，同时算（每个作业 10 个核）比先后算（每个作业20 个核）要更节约时间。 LREAL 控制赝势能量的非局域部分如何计算（k空间或实空间）。大体系，如果用集群算，可以考虑让NPAR=节点的CPU数。LREAL对k-空间，计算量正比于平面波数目（8ENCUT3/2*a1*a2*a3）。小于25个原子，可用K空间。实空间计算依赖体系大小。 对大体系用 LREAL = Au
3、to or LREAL = .Ture. 。KPAR设置 KPAR 为计算节点数或K 点数。 KPAR 用来设置 K 点的计算并行度。 每个 K 点用 N/KPAR 个核来计算，N 为总核数。核数很多时（ 100） ，这个参数的影响比较大。NELECTNELECT= realDefaultNELECT =- （number of valence electrons） 价电子数NELECT = number of electrons 电子数Usually you should not set this line - the number of electrons is determined au
4、tomatically from POTCAR （ZVALof the element） and POSCAR （number of the atoms of the respective atom type）.通常这个参数是不需要手动设置的 电子数将自动的从 POTCAR （元素的类型） 和 POSCAR 中确定 （每种类型原子数）If the number of electrons is not compatible 不匹配 with the number derived from 来源于 the valence 化合价 and the number of atoms a homogene
5、ous均匀的、 同类的 background-charge is assumed.如果电子数与来源于化合价的电子数不一样，就假定这些数目的原子具有均匀的背景电荷。If the number of ions specified 指定 in the POSCAR file is 0 and NELECT=n, then the energy of a homogeneous LDA-electron gas is calculated.如果在 POSCAR 中指定的离子数是0， NELECT=n ，然后计算一个同质LDA 电子气。在计算中遇到算不动的情况，需要适当减小精度。可以减小收敛精度的一些参
6、数设置：k 点密度 减少 K 点改变迭代算法（ALGO ）提高高斯展宽（SIGMA 增加）设置自洽延迟（NELMDL ）截断能 ENCUT大体系，一般 ALGO=Very_Fast ， IALGO=48SIGMA 默认 0.2在一开始计算时初始化的过程中电子非自洽迭代的步数确定平面波的切断动能。PREC确定计算的精度， 它决定了 ENCUT 和 ROPT 的默认值。 默认设为 Medium 中等的， VASP4.5以后的版本可设置的值为 Normal 普通和 Accurate 精确有两个提高并行效率的参数NPAR 和 KPAR大体系一般不进行收敛测试了，主要根据小体系的测试值和别人文章的使值来
7、选择比较合适的参数分步优化： 先采用低精度进行优化，比如增大离子步长，降低收敛精度，等收敛之后再提高精度进行精优化，这样相对比较可能会快一些对于大的体系，还是分步来优化好:1） 先把 EDIFF EDIFFG 弄大一点，设置 POTIM=0.1 ， NSW 的步数也设小一点， IALGO=48 ， ISIF=2 ，采 用 2*2*1 的 k 点先做一个优化。2） 然后把 EDIFF EDIFFG 弄小一点，设置 POTIM=0.2 ， NSW=1000 ，把 IALGO 去掉，采用大一点的 k 点， ISIF=4 （二维） ，进行第二次优化。截断能对计算速度的影响有多大？截断能取得大，电子波函
8、数的平面波展开中就要取更多的项。平面波越多，计算量就越大。截断能是表示平面波展开后取到多大能量的平面波，对于高能部分，展开后所占的比例非常小，而且影响计算速度，所以并不是截断能越大越好。这是一个两相竞争 的问题。一般用 POTCAR 中的 ENMAX 即可，若做静态等计算，可适当调大些。控制参与计算的平面波的个数，截断能越大，所参与的平面波就越多，计算的精度当然会随之提高，然而，在超过一定的数值后，计算的结果没有实质的变化，反而计算的时间却增加了。这个一定要测试的。还是要把计算理论好好看看的。后面需要用到的次数可能还会很多，科研不是做一天两天的事情。完全有必要找一个合理的截断能。否则写文章的时候，你直接写个很大的截断能（准确性没问题） ，但是会让人怀疑你对计算的理解，其他的结果都可能被怀疑了。一般，认为截断能越大， k 点设置越多，收敛精度越高，收敛会越慢，但是准确度不一定越高 。因此有必要做一些此类参数的测试。一般会这些参数大概会有个能量收敛的趋势。就是设置到某个临界参数，能量几乎变化很小。而参数再增加，只是增加计算量，而准确度没有变化。优化时可以用 ISMEAR=0 ，最后静态的时候再用-5 ，这样有些错误就不会发生。
