vasp做大体系结构优化

1、Vasp对大体系的结构优化ENCUTvasp运行时间主要消耗在对角化上。运行时间正比于 nbands*nplNwttr者是能带数目，后者是平面波数目；由于nplnwencuT2,故运行时间正比于encuT运行时间也正比于 nelecT，电子数目三次方。IALGO选择对角化算法：对于小体系，用 IALGO=38 D Davidson algorithm ）;对于大体系，用 IALGO=48 （RMM-DIIS）。可以设置 ALGO=veryfast or fast。RMM-DIIS 并行效率比 Davidson algorithm 高一些。NPAR如果 IALGO=38 贝U NPA做 1。对I

2、ALGO=48,影响不是特另I大，可选 2或4,可选节点数，取值越大，内存占用越多。并行效率总是低于 线性叠加效率的，核越多并行效率越低。所以对于一定核（如 20）, 一定作业（如2个），同时算（每个作 业10个核）比先后算（每个作业 20个核）要更节约时间。LREAL控制震势能量的非局域部分如何计算（k空间或实空间）。大体系，如果用集群算，可以考虑让NPAR节点的CPU数。LREAL对k-空间，计算量正比于平面波数目（8ENCUT3/2*aJa2*a3）。小于25个原子，可用K空间。 实空间计算依赖体系大小。对大体系用 LREAL = Auto or LREAL = .Ture.。KPAR设

3、置KPARJ计算节点数或K点数。KPARffi来设置K点的计算并行度。每个 K点用N/KPAR个核来计算，N 为总核数。核数很多时（100）,这个参数的影响比较大。NELEcTNELECT= real DefaultNELECT =- （number of valence electrons） 价电子数NELECT = number of electrons 电子数Usually you should not set this line - the number of electrons is determined automatically from POTCAR （ZVAL of the

4、element） and POSCAR （number of the atoms of the respective atom type）.通常这个参数是不需要手动设置的一电子数将自动的从POTCAR元素的类型）和 POSCA即确定（每种类型原子数）If the number of electrons is not compatible 不匹配 with the number derived from 来源于 the valence 化合价 and the number of atoms a homogeneous 均匀的、 同类的 background-charge is assumed.

5、如果电子数与来源于化合价的电子数不一样，就假定这些数目的原子具有均匀的背景电荷。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.如果在POSCA珅指定的离子数是0, NELECT=n然后计算一个同质 LDA电子气 在计算中遇到算不动的情况，需要适当减小精度。可以减小收敛精度的一些参数设置：k点密度 减少K点大体系，一般 ALGO=Very_Fast,IALGO=48SIGM

6、AN认 0.2改变迭代算法（ALGO 提高高斯展宽（SIGMA曾加）设置自洽延迟（NELMDL截断能ENCUT在一开始计算时初始化的过程中电子非自洽迭代的步数 确定平面波的切断动能。PREC确定计算的精度，它决定了 ENCU丽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、越大越好。这是一个两相竞争的问题。一般用POTCA呻的ENMAXiR可，若做静态等计算，可适当调大些。控制参与计算的平面波的个数，截断能越大，所参与的平面波就越多，计算的精度当然会随之提高，然而， 在超过一定的数值后，计算的结果没有实质的变化，反而计算的时间却增加了。这个一定要测试的。还是要把计算理论好好看看的。后面需要用到的次数可能还会很多，科研不是做一天 两天的事情。完全有必要找一个合理的截断能。否则写文章的时候，你直接写个很大的截断能（准确性没 问题），但是会让人怀疑你对计算的理解，其他的结果都可能被怀疑了。一般，认为截断能越大，k点设置越多，收敛精度越高，收敛会越慢，但是准确度不一定越高。因此有必要做一些此类参数的测试。一般会这些参数大概会有个能量收敛的趋势。就是设置到某个临界参数，能量 几乎