VASP(1)_参数测试_截断能(ENCUT)测试

必要输入文件：

encut.sh
POTCAR
KPOINTS和POSCAR可以在encut.sh里设置，也可以独立给出。

encut.sh脚本参考

#!/bin/bash#生成KPOINTS文件，KPOINTS不能与INCAR中的KSPACING参数共用，二选一
#注：KPOINTS文件每行一定要顶格写，不能前面不能空格，不然可能会报错。
cat > KPOINTS <<!
A
0
M
9  9  9
0  0  0
!#产生计算所需POSCAR
cat > POSCAR<<!
Si8
1.00000000005.4687280000        0.0000000000        0.00000000000.0000000000        5.4687280000        0.00000000000.0000000000        0.0000000000        5.4687280000
Si
8
Direct0.0000000000        0.0000000000        0.00000000000.2500000000        0.7500000000        0.75000000000.5000000000        0.0000000000        0.50000000000.0000000000        0.5000000000        0.50000000000.5000000000        0.5000000000        0.00000000000.7500000000        0.2500000000        0.75000000000.7500000000        0.7500000000        0.25000000000.2500000000        0.2500000000        0.2500000000
!for i in $(seq 400 50 700) #截断能从450-800，步数为50
do#生成vasp静态计算的INCAR
cat > INCAR_static <<!
Global ParametersISTART =  0            (Read existing wavefunction; if there)ICHARG =  2            (Non-self-consistent: GGA/LDA band structures)LREAL  = .FALSE.       (Projection operators: automatic)ENCUT  =  $i           (Cut-off energy for plane wave basis set, in eV)PREC   =  Accurate     (Precision level)  LWAVE  = .FALSE.       (Write WAVECAR or not)LCHARG = .FALSE.       (Write CHGCAR or not) ADDGRID= .TRUE.        (Increase grid; helps GGA convergence)Electronic RelaxationISMEAR =  -5           (Gaussian smearing; metals:1)#SIGMA  =  0.05        (Smearing value in eV; metals:0.2)NELM   =  60           (Max electronic SCF steps)  NELMIN =  4            (Min electronic SCF steps)EDIFF  =  1E-06        (SCF energy convergence; in eV)GGA  =  PE             (PBEsol exchange-correlation)Ionic RelaxationISIF   =  2            (Stress/relaxation: 2-Ions, 3-Shape/Ions/V，4-Shape/Ions)EDIFFG = -0.001        (Ionic convergence; eV/AA)KSPACING = 0.10
!cp INCAR_static INCAR
echo "ENCUT = $i eV" ; time mpirun -np 16 vasp_std #vasp并行运行命令，根据系统自行修改rm INCAR_static
rm INCAR
#提取计算得到的能量
E=$(grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.9f \n", $5 }')
echo $i $E >>encut_energy.out
done

提交脚本任务运行，成功结束后可以得到输出文件encut_enengy.out 如下

400 -43.398923980
450 -43.399878560
500 -43.399694690
550 -43.399861720
600 -43.400281930
650 -43.400531580
700 -43.400560950

可以使用OriginPro等画图工具画出截断能-能量关系.

当能量的变化在大约0.001 eV/atom 时差不多认为是收敛的。这里选择硅的截断能 500 eV 作为后续的计算。

VASP(1)_参数测试_截断能(ENCUT)测试相关推荐

VASP(1)_参数测试_晶格常数测试
晶格常数测试 (Equation of state method) 必要输入文件 run_a0.sh POTCAR INCAR 和 KPOINTS可在run_a0.sh中直接设置,也可以单独给出. E ...
python table类中的可以作为初步观察工具_Python程序设计基础_课程2020最新章节测试网课课后答案...
Python程序设计基础_课程2020最新章节测试网课课后答案更多相关问题科学探究的过程一般要经过______.作出假设.制定方案.实施方案,得出结论等五个步骤．在生物学研究中,当实验的结果与原 ...
21天Jmeter打卡Day17 后置处理器_JSON_正则表达式_边界提取器_完成删除场景模拟测试
1天Jmeter打卡Day17 后置处理器_JSON_正则表达式_边界提取器_完成删除场景模拟测试 https://www.jianshu.com/p/5c1d64e5d724 Json提取器见Day ...
软件测试过程与方法_第1周小测试
软件测试过程与方法_第1周小测试 <软件测试过程与方法>作业第一题:单项选择题(每小题5分,共60分) 三.判断题(每小题5分,共40分) <软件测试过程与方法>作业第一题 ...
复习Java异常处理_异常分类_自定义异常_线程初步了解
复习Java异常处理_异常分类_自定义异常_线程主要内容异常.线程教学目标第一章异常 1.1 异常概念异常,就是不正常的意思.在生活中:医生说,你的身体某个部位有异常,该部位和正常相比有点 ...
李宏毅_机器学习_作业1（详解）_COVID-19 Cases Prediction (Regression)
李宏毅_机器学习_作业1:COVID-19 Cases Prediction (Regression) 本文旨在,读懂hw1代码.李宏毅老师作业的完成有一个默认前提,掌握python基础语法,pyto ...
ava_212_反射机制_动态操作_构造器_方法_属性_练习
package java_212_反射机制_动态操作_构造器_方法_属性_练习; /** 属性 */ public class User { private int id; private int a ...
Python_机器学习_算法_第4章_4.决策树算法
Python_机器学习_算法_第4章_4.决策树算法文章目录 Python_机器学习_算法_第4章_4.决策树算法决策树算法学习目标 4.1 决策树算法简介学习目标小结 4.2 决策树分类原 ...
停车场管理_充电桩_停车收费_物业管理_源码
[声明]:本项目里面的代码没有任何私jar包,本来就是开源项目,代码没有丢包或者故意丢代码导致项目报错的情况,但是不保证能商用,个人拿来作为入门研究和快速二次开发是个不错的选择!这项目代码肯定能跑起来 ...

VASP(1)_参数测试_截断能(ENCUT)测试

VASP(1)_参数测试_截断能(ENCUT)测试相关推荐

最新文章

热门文章