本文介绍了如何利用VASP进行金刚石硬度的计算，通过设置ENCUT为500和K mesh为8*8*8进行高精度计算。采用Murnaghan状态方程，通过不同细胞尺寸对应的能量值拟合得出体积弹性模量Bulk modulus。详细步骤包括脚本修改、数据获取、二次曲线拟合和使用murn.f程序确定关键参数V0和E0。
K Mesh的选取
这里我们正式计算一个体积弹性模量bulk modulus，用精度较高的ENCUT＝500和K mesh 8*8*8
那个run.sh脚本改成这个样子#!/bin/bash
for i in 1.71 1.72 1.73 1.74 1.75 1.76 1.77 1.78 1.79 1.8 1.81 1.82 1.83 1.84
do
cat > POSCAR > /dev/null
E=`grep TOTEN OUTCAR | tail -n 1 | awk '{print($5)}'`
echo $i $E
done复制代码运行 run.sh | tee run.out
可以得到每个cell尺寸对应的体系能量值，根据Murnaghan状态方程 http://en.wikipedia.org/wiki/Bir ... n_equation_of_state
E=E0(B/2V)(V-V0)^2 (*)
即可拟合出B的数值，即为体系的Bulk modulus，体积随压强的变化性能。其中V0，E0是平衡点的体积和压强值。
该系列点近似是一个二次曲线，像下面这个图一样。先求出这个曲线的对称轴和最低点吧，也就是V0和E0，图里的B点是bay自己求出的，方法见下-17.9 ---A-------------------------------------------------------
"run.out" A
| |
-17.95
| |
| A
最低0.47元/天 解锁文章
9031
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
