VASP计算实例Word格式文档下载.docx
《VASP计算实例Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《VASP计算实例Word格式文档下载.docx（8页珍藏版）》请在冰豆网上搜索。
12
晶格常数计算结果
13
三、
Si元素单原子能量计算
14
由内聚能倒推单原子能量
15
单原子能量计算
四、
Si的VASP力学常数计算
16
计算所需文件
计算与数据处理
17
计算所用到的公式：
18
五、
SI晶体的电子结构
19
采用VASP计算能带的步骤
电荷分布计算结果
20
能带计算和结果
21
态密度计算和结果
六、
Si晶体介电函数和光学性质的计算
22
计算步骤
用到的文件
23
计算结果
26
七、
VASP的声子谱计算
29
30
声子谱、声子态密度计算和结果
33
热学性质计算和结果
34
八、
化合物
键长计算
35
一、氢气分子H2键长的计算
1.基本文件
准备基本文件INCAR、POTCAR、POSCAR、KPOINT以及脚本文件encut、k、optimize
2.赝势类型的选择
输入文件如下其中参数要靠经验初选
INCAR:
System=F2
ISTART=0
ICHARG=2
NELMDL=5
ISMEAR=0
SIGMA=0.1
PREC=Accurate
KPOINTS：
Automaticmesh
0
M
1
1
1
0
0
0
POSCAR：
O
1
15.0
0.00
0.00
0.00
14.0
13.0
D
分别选用五个贋势文件进行计算。
计算方法：
把INCAR、POSCAR、POTCAR、KPOINTS四个文件分别放入同一个文件夹中，在终端中输入“VASP”得到五个能量。
GGA-91
E0=-.36620018E+01
GGA-CA
E0=-.35466314E+01
PAW
E0=-.34415634E+01
PAW-GGA
E0=-.34948621E+01
PBE
E0=-.35753416E+01
选择单原子能量最低的，晶体最稳定，故选择GGA-91赝势文件POTCAR。
3.截断能ENCUT参数的选择
准备文件KPOINTS、POSCAR、POTCAR以及脚本文件encut并放入同一个文件夹中
encut：
#!
/bin/sh
rm
WAVECAR
for
i
in
550
600
650
700
750800850900
do
cat
>
INCAR
<
!
System
=
H2
ISTART=0
ISMEAR=0;
SIGMA=0.1
ICHARG
2
ENCUT=$i
PREC=Accurate
echo
"
ENCUT
$i
eV"
;
vasp
E=`grep
TOTEN"
OUTCAR
|
tail
-1
awk
'
{printf
%12.6f\n"
$5}'
`
KP=`grep
irreducible"
%5i
\n"
$2
}'
$KP
$E
comment
done
Automatic
mesh
M
1
F2
1.209
4.000.00
0.004.00
0.000.00
4.00
Direct
0.25
选用上次计算是最低能量的POTCAR
运行脚本文件encut，得到comment文件内容如下：
5501-3.654087
6001-3.654701
6501-3.651275
7001-3.652794
7501-3.652555
8001-3.652600
8501-3.653120
9001-3.652554
9501-3.652669
总能变化在0.001eV左右就足够了，因此我们选择encut=700eV
4.KPOINTS参数选择
准备文件POSCA、INCAR、POTCAR以及脚本文件k，POSCAR、POTCAR与上面相同
脚本文件
k：
3
5
7
9
11
1315
KPOINTS
$i
$i
k
mesh
x
$i"
%12.6f
$5
INCAR：
O2
700
ISTART
ISMEAR
SIGMA
0.1
PREC
Accurate
运行脚本文件k，得到comment文件：
36-3.652382
518-3.652411
740-3.652388
975-3.652386
总能变化在
0.001eV
左右就足够了，因此我们选择
作为
点值。
5.对晶格常数进行优化
准备文件INCAR、KPOINTS、POTCAR以及脚本文件optimize，其余文件与上一样
Optimize：
BIN=/usr/bin/vasp
-rf
1.161.211.261.311.361.411.461.511.561.611.66;
POSCAR
$i
4.00
4.00
a=
$BIN
E=`tail
OSZICAR`;
RESULT.fcc
RESULT.fcccat
运行脚本文件optimize，得到RESULT文件如下：
1.161F=-.15006150E+01E0=-.15006150E+01dE=-.179658E-31
1.211F=-.24084278E+01E0=-.24084278E+01dE=-.145599E-16
1.261F=-.30094870E+01E0=-.30094870E+01dE=-.247133E-21
1.311F=-.33821346E+01E0=-.33821346E+01dE=-.403048E-16
1.361F=-.35847564E+01E0=-.35847564E+01dE=-.800387E-69
1.411F=-.36607958E+01E0=-.36607958E+01dE=-.270236E-21
1.461F=-.36424044E+01E0=-.36424044E+01dE=-.866906E-24
1.511F=-.35553122E+01E0=-.35553122E+01dE=-.289397E-27
1.561F=-.34180694E+01E0=-.34180694E+01dE=-.998984E-16
1.611F=-.32473112E+01E0=-.32473112E+01dE=-.696853E-23
1.661F=-.30547759E+01E0=-.30547759E+01dE=-.540502E-15
利用计算数据拟合二次曲线最低点就是计算出的晶格常数。
拟合后的曲线如下图：
计算出晶格常数是1.41实验室是同样是1.41故参数选择合适计算准确。
二、Si晶体晶格常数计算
1.赝势类型选择
分别带入五个赝势库文件带入计算去能量最小值的POTCAR。
输入文件：
System=Si
ISTART=0
=
240
NELM=100
EDIFF=1E-04
EDIFFG=1E-03
NSW=5
IBRION=2
升级会员 