BatchDoc Word文档批量处理工具（计算前的）验证一、检验赝势的好坏：（一）方法：对单个原子进行计算；（二）要求：1、对称性和自旋极化均采用默认值；2、ENCUT要足够大；3、原胞的大小要足够大，一般设置为15 ?足矣，对某些元素还可以取得更小一些。
（三）以计算单个Fe原子为例：1、INCAR文件：SYSTEM = Fe atomENCUT = 450.00 eVNELMDL = 5 ! make five delays till charge mixing，详细意义见注释一ISMEAR = 0SIGMA=0.12、POSCAR文件：atom15.001.00 0.00 0.000.00 1.00 0.000.00 0.00 1.001Direct0 0 03、KPOINTS文件：（详细解释见注释二。
）AutomaticGamma1 1 10 0 04、POTCAR文件：（略）注释一：关键词“NELMDL”:A）此关键词的用途：指定计算开始时电子非自洽迭代的步数（即NELMDL gives the number of non-selfconsistent steps at the beginning），文档批量处理工具BatchDoc Word文档批量处理工具BatchDoc Word densitycharge fastermake calculations 。
目的是“非自洽”指的是保持“非自Charge density is used to set up the Hamiltonian, 所以不变，由于洽”也指保持初始的哈密顿量不变。
:B）默认值（default value）(时) 当ISTART=0, INIWANELMDL = -5 V=1, and IALGO=8 ) ISTART=0, INIWA V=1, and IALGO=48( NELMDL = -12 时当)其他情况下NELMDL = 0 (NELMDL might be positive or negative.ionic each applied means A positive number that after a delay is（movement -- in general not a convenient option. ）在每次核运动之后（只在A negative value results in a delay only for the start-configuration.第一步核运动之前）NELMDL”为什么可以减少计算所需的时间？C）关键词“the the is Charge density used Hamiltonian, to set then upwavefunctions are optimized iteratively so that they get closer to the exacta optimized wavefunctions wavefunctions of Hamiltonian. this From theold with density charge is calculated, the which is then mixed newManual P105input-charge density. A brief flowchart is given below.（参自页）是比较离谱的，在前一般情况下，the initial guessed wavefunctions不变、保持初始的density次非自洽迭代过程中保持NELMDLcharge哈密顿量不变，只对wavefunctions进行优化，在得到一个与the exact文档批量处理工具BatchDoc WordBatchDoc Word文档批量处理工具wavefunctions of initial Hamiltonian较为接近的wavefunctions后，再开始同时优化charge density。
这样一来，计算时间要比一开始就同时优化charge density 和wavefunctions短得多。
注释二：为什么这里只需要一个k点？For atoms and molecules, the Bloch theorem does not apply, hence there is no need to use more than one single k-point. When more k-points are used,only the interaction between the atoms (which should be zero) is describedmore accurately.（三）计算任务执行方法：输入：vasp（四）赝势好的判断标准：计算得到的OUTCAR文件中的“energy without entropy”能量值在?0.001～?0.01 eV之间。
文档批量处理工具BatchDoc WordBatchDoc Word文档批量处理工具二、筛选合适的ENCUT大小：（一）输入文件：1、用脚本程序optencut.sh代替INCAR文件：rm WAVECARfor i in 100 150 200 250 300 350 400 450 500 550 600docat > INCAR <<!SYSTEM = bcc FeENCUT = $iISTART = 0 ; ICHARG = 2ISMEAR = -5PREC = Accurate!echo ENCUT = $i eV; time vaspE=$(grep TOTEN OUTCAR | tail -1| awk '{printf _x0012_.6f \n, $5}')echo $i $E >> commentdone2、POSCAR文件：bcc Fe2.8661.00 0.00 0.000.00 1.00 0.000.00 0.00 1.002Direct0 0 00.5 0.5 0.53、KPOINTS文件：Automatic generationMonkhorst-Pack9 9 90.0 0.0 0.04、POTCAR文件：（略）（二）计算任务执行方法：输入：dos2unix optencut.shbash optencut.sh（三）判别标准：计算完成后得到comment文件，它列出了在每个ENCUT时计算得到的相应的总能，只要总能变化在0.001 eV左右就足够了。
文档批量处理工具BatchDoc WordBatchDoc Word文档批量处理工具三、选择合适的k点数目：（一）输入文件：1、INCAR文件：SYSTEM = bcc FeENCUT = 450.00 eVISTART = 0 ; ICHARG = 2ISMEAR = -5PREC = Accurate2、POSCAR文件：bcc Fe2.8661.00 0.00 0.000.00 1.00 0.000.00 0.00 1.002Direct0 0 00.5 0.5 0.53、用脚本程序optkpoints.sh代替KPOINTS文件：rm WAVECARfor i in 5 7 9 11 13 15docat > KPOINTS <<!Automatic generationMonkhorst-pack$i $i $i0.0 0.0 0.0!echo k mesh = $i x $i x $i; time vaspE=$(grep TOTEN OUTCAR | tail -1| awk '{printf _x0012_.6f \n, $5}')KP=$(grep irreducible OUTCAR | tail -1| awk '{printf %5i \n, $2}')echo $i $KP $E >> commentdone4、POTCAR文件：（略）（二）计算任务执行方法：输入：dos2unix optkpoints.shbash optkpoints.sh（三）判别标准：计算完成后得到comment文件，它列出了在k点数目与总能的对应值，只要总能变化在0.001 eV左右就非常足够了。
文档批量处理工具BatchDoc WordBatchDoc Word文档批量处理工具四、优化选择合适的SIGMA值（展宽σ值）：（一）为什么要优化SIGMA值？若展宽σ太小，则计算难以收敛；若展宽σ太大，则会产生多余的熵(entropy)，因此必须选择合适的σ值。
(Too large smearing-parameters might result in a wrong total energy, small smearing parameters require a large k-point mesh.)（二）ISMEAR和SIGMA：1、ISMEAR和SIGMA这两个关键词要联合起来使用，前者用来指定smearing 的方法，后者用来指定smearing的展宽——σ值。
2、ISMEAR和SIGMA的默认值分别为1和0.2。
3、ISMEAR可能的取值为-5，-4，-3，-2，-1，0，N (N表示正整数)：ISMEAR＝-5，表示采用Blochl修正的四面体方法；ISMEAR＝-4，表示采用四面体方法，但是没有Blochl修正；ISMEAR＝-1，表示采用Fermi-Dirac smearing方法；ISMEAR＝0，表示采用Gaussian smearing方法；ISMEAR＝N，表示采用Methfessel-Paxton smearing方法，其中N是表示此方法中的阶数，一般情况下N取1或2, 但是In most casesandleads to very similar。
results4、σ值一般在0.1～0.3 eV范围内。
5、ISMEAR取值的一些经验：（1）一般说来，无论是对何种体系，进行何种性质的计算，采用ISMEAR＝0并选择一个合适的SIGMA值，都能得到合理的结果。
（2）在进行静态计算（能量单点计算, ）或态密度计算且k点no relaxation in metals 数目大于4时，取ISMEAR＝-5。
（3）当原胞较大而k点数目较小（小于4个）时，取ISMEAR＝0，并选择一个合适的SIGMA值。
(if the cell is too large (or if you use only a single or two k-points) use ISMEAR=0 in combination with a small SIGMA=0.05)（4）对半导体或绝缘体，不论是静态还是结构优化计算，都取ISMEAR＝-5。