0S.0I0G0M00A00=000.010000000 0.5000000000000000 0.4589585899816485
0#.5I0S0I0F00=02000000000 0.0000000000000000 0.9589585899816486
0I.5A0L0G00O00=04080000000 0.0000000000000000 0.5410414100183514
POTCAR输入文件: 赝势文件
Si 的一种势函数的部分内容
US Si 4.00000000000000000 parameters from PSCTR are: VRHFIN =Si: s2p2 LEXCH = CA EATOM = 115.7612 eV, 8.5082 Ry GGA = -1.4125 -1.4408 .0293 -.9884 eV
First-principles electronic structure calculations
• 简介 • 基本任务 • 输入文件 • 输出文件 • 程序举例 • 操作界面
简介—VASP是什么
• 全称Vienna Ab-inito Simulation Package • 是一个采用平面波赝势(或缀加投影波)方法进行
0.5000000000000000 0.0000000000000000 0.7500000000000000
(2). 优化晶格参数
• 运行VASP，从输出文件CONTCAR中查看优化之后的晶格 常数。
（3）固定晶格常数静态计算，得到后面计算要用到的电荷 密度。
TiO2
1S.0Y0S0T0E00M00=0T00iO0020
# ISIF = 3 IALGO = 48 NELMDL = -5
0.0000000000000000 0.5000000000000000 0.0410414100183514
0.5000000000000000 0.5000000000000000 0.2910414100183515 将CONTCAR
0.0000000000000000 0.0000000000000000 0.7910414100183514 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.5000000000000000 0.2500000000000000
basic characters of surface system 表面体系的基本性质的计算
......
输入文件
POTCAR KPOINTS POSCAR INCAR
pseudopotentail file（赝势文件，软件本身具有， 用时选择合理的即可）
Brillouin zone sampling（ k点取样设置文件）
TITEL = US Si LULTRA = T use ultrasoft PP ? IUNSCR = 1 unscreen: 0-lin 1-nonlin 2-no RPACOR = 1.580 partial core radius POMASS = 28.085; ZVAL = 4.000 mass and valenz RCORE = 2.480 outmost cutoff radius RWIGS = 2.480; RWIGS = 1.312 wigner-seitz radius (au A) ENMAX = 150.544; ENMIN = 112.908 eV EAUG = 241.945 …………
从头模拟的软件包。 • 是目前材料模拟和计算材料科学研究中非常流行
的商用软件之一。
基本任务
crystal structure 晶体结构基本性质计算
• Electronic Structure • Magnetic properties • Optimize parameters • Structural relaxation
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.5000000000000000 0.5000000000000000 0.5000000000000000
0.0000000000000000 0.5000000000000000 0.2500000000000000
中得到的数据 拷贝到
POSCAR中
0.5000000000000000 0.0000000000000000 0.7500000000000000
(4). 静态计算之后，进行能带计算
SYSTEM = TiO2 PREC = Accurate
NSW = 0 ENCUT = 400 ISTART = 1 ICHARG = 11 # IBRION = -1 # ISMEAR = -5
输出文件
OUTCAR :最主要的输出文件 DOSCAR : 电子态密度文件 EIGENVAL :本征值文件 OSZICAR :每次迭代或原子迟豫(或MD)的信息 CHG和CHGCAR :电荷密度文件 WAVECAR :波函数文件 CONTCAR :原子迟豫或MD后的体系结构文件 IBZKPT :布里渊区中的k点 PCDAT :对关联函数 XDATCAR :在MD时, 原子位置变化的跟踪文件 PROCAR和PROOUT : 波函数投影或分解的文件 LOCPOT :总的局域势 ELFCAR :电子局域函数
每类原子的个数 确定按何种坐标来写原子位置 第一个坐标 第二个坐标
TiO2_anatase.cell 隐藏文件
缩放倍数为1
命名为TiO2 基矢
体系中有8TS输入文件: 控制K点的选取方式
• 按普通格式手动输入各个k点(推荐用在计算能带时)
(1). 生成4个输入文件: POSCAR POTCAR KPOINTS INCAR
TiO2
1.00000000000000
3.776000000000000 0.000000000000000 0.000000000000000 0.00S0Y0S0T0E00M00=0T00iO002 3.776000000000000 0.000000000000000 D0008...i050r4000e0000ct.00000000PNEIII000CBS0RNS000T0HRE000WCA0IA000C0UOR000R=0TNT=0000G00020A===00000=c00004001c000020u00000ra0t0000e...0550000000000.000000000000000000000000KM0000040-000o0P000n0o000k0i40000nh0000to0sr0s000400t ...0P274a005c888k0009000.000400080006000000000000V各0000A0000种000S0000P0P000提0O0T供0C0的1AR 0.50000I0S0IF00=00300000 0.0000000000000000 0.9580000000000000 0.50000E0D00IF0F00G0=00-000.010.0000000000000000 0.5420000000000000 0.00000A0D00D0G00R0I0D0=00.Tr0u.e5.000000000000000 0.0420000000000000 0.50000I0S0M00E0A0R00=0000 0.5000000000000000 0.2920000000000000 0.00000S0I0G0M00A00=000.010 0.0000000000000000 -0.2080000000000000
structural data（描述体系结构的文件）
steering parameters（计算控制参数文件）
INCAR输入文件：设置程序控制参数的值
System 、ISTART、ENCUT、NELM、EDIFF 、EDIFFG、 GGA 、 NPAR、NSW、 IBRION 、 ISIF、 ISYM、 LWAVE 、LCHARG等
POSCAR输入文件：描述体系结构
例:SiC体系的POSCAR文件
Cubic SiC 3.57 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 11 Direct 0.00 0.00 0.00 0.25 0.25 0.25
设置体系的名称 晶格常数或缩放系数
原胞(或晶胞)的基矢
0.00000000000000000 9.518838639V6A3S40P2提63供的
Monkhorst Pack
各种POTCAR
444
0I.0B0R0I0O0N00=00-01000000 0.0000000000000000000 0.2089585899816485
0I.5S0M00E0A0R00=000000000 0.5000000000000000 0.7089585899816486
P3R.7E8C45=10A1c3c4u2ra2t5e9402 0.0000000000000000 0.0000000000000000
N0S.0W00=00000000000000 3.784510K1-3P4o2in2t5s9402 0.0000000000000000
E0N.0C00U0T00=0400000000000 8IST4ART = 0 DiIreCcHt ARG = 2
例：静态计算(计算总能和自洽的电荷密度)时
SYSTEM = Silicon ENCUT = 350 ISTART = 0 ICHARG = 2 ISMEAR = -5 EDIFF = 1E-5 PREC = Accurate