计算的体系经过结构优化后,需要进行静态计算,才能进行的后续电子结构分析计算。静态计算时,原子位置不动,对电子进行优化,使体系能量达到最低。
计算准备
POSCAR
复制结构优化得到CONTCAR 作为静态计算的POSCAR
KPOINTS , POTCAR
与结构优化的KPOINTS和POTCAR保持一致
INCAR
修改结构优化INCAR的部分参数,保存好得到的电子密度和波函数文件,提供给后续电子结构计算,提高计算速度。
System =
ISTART = 0
ICHARG = 2
ENCUT = 450
PREC = Medium
ISIF = 2
EDIFF = 1E-4 (电子收敛标准 1E-4 或者1E-5)
ISPIN = 2
NSW = 0 (最大离子步优化设置为0)
POTIM = 0.2
EDIFFG = -0.03
IVDW = 11
GGA = RP
LCHGCAR=.TRUE (输出电荷密度文件, CHG和CHGCAR)
LWAVE=.TRUE (输出波函数文件,WAVECAR)
NELM=60 (电子优化最大步数)
计算结构
静态计算得到的CHGCAR, CHG, WAVECAR 可以直接复制到后续的计算中。但是,静态计算的能量没有意义。
实际操作
提交结构优化计算
mkdir OPT STA #建立结构优化文件夹OPT和静态计算文件夹STA
导入输入文件
cd OPT
ls #打开OPT文件夹,检查输入文件
文件夹里输入问价齐全后可以提交计算
2.提交静态计算文件
结构优化完成后,可以进行静态优化。
cd.. #返回上一级
cd STA #打开静态计算文件夹
导入静态计算INCAR-static
cp ../OPT/KPOINTS . #复制结构优化文件OPT中的KPOINTS
cp ../OPT/POTCAR . #复制结构优化文件OPT中的POTCAR
cp ../OPT/opt.lsf sta.lsf #复制结构优化文件OPT中的脚本文件,并重新命名
vi sta.lsf #查看脚本文件内容
输入i,
将cp INCAR-OPT INCAR 修改与输入文件名称一致 cp INCAR-static
:wq (保存修改)
确认修改无误
cp ../OPT/CONTCAR POSCAR #复制结构优化文件夹OPT中的CONTCAR为POSCAR
输入文件准备完毕后,提交计算即可。
本文涉及的文件可以点击下载
本文由一嘟向日葵投稿。