VASP计算前的各种测试
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
BatchDoc Word文档批量处理工具
（计算前的）验证
一、检验赝势的好坏：
（一）方法：对单个原子进行计算；
（二）要求：1、对称性和自旋极化均采用默认值；
2、ENCUT要足够大；
3、原胞的大小要足够大，一般设置为15 ?足矣，对某些元素还可以取得更小一些。
（三）以计算单个Fe原子为例：
1、INCAR文件：
SYSTEM = Fe atom
ENCUT = 450.00 eV
NELMDL = 5 ! make five delays till charge mixing，详细意义见注释一
ISMEAR = 0
SIGMA=0.1
2、POSCAR文件：
atom
15.00
1.00 0.00 0.00
0.00 1.00 0.00
0.00 0.00 1.00
1
Direct
0 0 0
3、KPOINTS文件：（详细解释见注释二。）
Automatic
Gamma
1 1 1
0 0 0
4、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 up
wavefunctions are optimized iteratively so that they get closer to the exact
a optimized wavefunctions wavefunctions of Hamiltonian. this From the
old with density charge is calculated, the which is then mixed new
Manual P105input-charge density. A brief flowchart is given below.（参自页）
是比较离谱的，在前一般情况下，the initial guessed wavefunctions
不变、保持初始的density次非自洽迭代过程中保持NELMDLcharge
哈密顿量不变，只对wavefunctions进行优化，在得到一个与the exact
文档批量处理工具BatchDoc Word
BatchDoc 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 described
more accurately.
（三）计算任务执行方法：
输入：vasp
（四）赝势好的判断标准：计算得到的OUTCAR文件中的“energy without entropy”能量值在?0.001～?0.01 eV之间。
文档批量处理工具BatchDoc Word
BatchDoc Word文档批量处理工具
二、筛选合适的ENCUT大小：
（一）输入文件：
1、用脚本程序optencut.sh代替INCAR文件：
rm WAVECAR
for i in 100 150 200 250 300 350 400 450 500 550 600
do
cat > INCAR <
SYSTEM = bcc Fe
ENCUT = $i
ISTART = 0 ; ICHARG = 2
ISMEAR = -5
PREC = Accurate
!
echo ENCUT = $i eV; time vasp
E=$(grep TOTEN OUTCAR | tail -1| awk '{printf _x0012_.6f \n, $5}')
echo $i $E >> comment
done2、POSCAR文件：
bcc Fe
2.866
1.00 0.00 0.00
0.00 1.00 0.00
0.00 0.00 1.00
2
Direct
0 0 0
0.5 0.5 0.5
3、KPOINTS文件：
Automatic generation
Monkhorst-Pack
9 9 9
0.0 0.0 0.0
4、POTCAR文件：（略）
（二）计算任务执行方法：
输入：dos2unix optencut.sh
bash optencut.sh
（三）判别标准：计算完成后得到comment文件，它列出了在每个ENCUT时计算得到的相应的总能，只要总能变化在0.001 eV左右就足够了。
文档批量处理工具BatchDoc Word
BatchDoc Word文档批量处理工具