VASP计算AgGaS2能带及态密度及光学性质文档格式.docx
- 文档编号:22587912
- 上传时间:2023-02-04
- 格式:DOCX
- 页数:19
- 大小:192.95KB
VASP计算AgGaS2能带及态密度及光学性质文档格式.docx
《VASP计算AgGaS2能带及态密度及光学性质文档格式.docx》由会员分享,可在线阅读,更多相关《VASP计算AgGaS2能带及态密度及光学性质文档格式.docx(19页珍藏版)》请在冰豆网上搜索。
ISTART= 0
LREAL=.FALSE.
PREC =Low
EDIFF =1e-4
EDIFFG=-0.03
IALGO =48
NELMIN=4
ISYM =0
GGA= PBE
ISPIN=1
NBANDS=120
OUTPUTCONTROL
LCHARG =.TRUE.
LWAVE =.TRUE.
LVTOT =.FALSE.
IONICRELAXATION
NBLOCK =1
NSW=1
IBRION=-1
DOSRELATED(disabled)
ISMEAR =0 (tetrahedron/gaussian/m-p)
SIGMA=0.05
4.再下面KPOINTS
A
0
M
444
000至此,四文件已准备好,进行计算,
mpiexec-np8~/bin/vasp.4.5-mk-mp-pgi <
/dev/null>
vasp.out&
(运行命令)
(问题:
计算前对某些参数的测试,比如截断能,晶格参数等,标准时什么?
vasp运行后,不知道结果是否满足要求呢?
)
第二部:
将构型优化后产生的CONTCAR文件拷贝为POSCAR文件
mpiexec-np8~/bin/vasp.4.5-mk-mp-pgi<
/dev/null>
vasp.out&(运行命令)
第三部:
能带的计算.建立新的计算目录
mkdir band 创建目录
cd band 进入目录
cp../INCAR.复制
cp ../PO*.
cp../CHG*.
对INCAR文件作如下修改:
ICHARG=11
NBANDS=120 与上面构型优化时一致
建立新的KPOINTS 手动定义K点
AgGaS2
51
Rec(关键字要对齐第一行)
0.5000 0.50000.50001.0000
0.4500 0.4500 0.4500 1.0000
... ...
K点坐标获得途径
布里渊区各特殊点的坐标可利用CASTEP获得
然后编辑
viinp
lin1 特殊点
lin2特殊点之间的取点数
lin3第一个K点的坐标
lin4 第二个K点的坐标
.. ..(同上..)
. ...
保存
运行 vasp_kpoints<
inp
就可以得到坐标
手动复制到KPOINTS
运行VASP
mpiexec -np8 ~/bin/vasp.4.5-mk-mp-pgi<
/dev/null >
vasp.out&
在运行
Vaspband
DrawthebandstructureofVASP
Twofiles, EIGENVALand OUTCARare needed
----------------------------------------------------------------
Spin-polarized calculation(T/F)?
(是否自旋极化计算)
f
Numberof linesalongtheBZ:
(布里渊区特殊线数目)
5
K pointsequenceusedtoseparatelines:
(各线段始终点序号)
11121314151
Setthe Fermileveltozero(T/F)?
(是否费米设置为零点)
t
可以得到相应的band.dat用origin绘图
然后标注各个特殊点
Fermienergy=1.2182v exp=2.76ev(可能是参数设计太粗糙…) 不过放大还算清楚..
(用MS作图 K网格相同)
第四部:
DOS的计算及其绘图
其POSCARPOTCARKPOINTS可以相同于第一步优化
设置INCAR
SYSTEM= optimizationofAgGaS2
NPAR=8
LPLANE= .TRUE.
Elecronicminimisation
ISTART=3 (从WAVECAR读入波函数)
LREAL=.FALSE.
PREC =Medium
EDIFF =1e-4
EDIFFG=-0.02
IALGO =48
NELMIN= 3
ISYM =1
ISIF=3
ISPIN=1
ISMEAR=0
SIGMA=0.1
NBANDS=120
LORBIT=11(投影到s.p.d)
NEDOS=501 (DOS的点数)
ICHARG=11 (从CHG读入电荷密度)
EMIN=-20
EMAX=20
OUTPUTCONTROL
LCHARG=.TRUE.
LWAVE =.TRUE.
LVTOT= .FALSE.
IONICRELAXATION
NBLOCK= 1
NSW = 1
IBRION=-1
POTIM =0.05
运行vasp 命令同上
然后编写DOS
F (是否自旋极化计算 )
9(轨道数)
6(要计算的DOS)
-1 (为负数时,表示所考察的为原子的DOS)
1 (原子编号)
-1 (为正正数时,表示为原子轨道的DOS)
2 (下面都同上)
-1
3
1
1
8
保存运行vaspdos(如果没有设定路径的话。
要自己指定路径)
可以得到fort.25文件
运行cdos98不过需要加一条DOS数
就可以得到DOS1DOS2DOS3. .. . ..
用origin分别绘图…..
(可能是DOS编号出现错误了的原因 导致Ga(p)和S(p)没有画出..)
第五部:
光学性质的计算
(一)线性的计算
1.也是进行构型优化(同能带计算时结构优化)..只是计算中INCAR文件需包含关键词:
ISIF =3
所有计算均需采用PAW型赝势
2.构型优化后将CONTCAR拷贝成POSCAR
在INCAR里面
NPAR=1(计算光学性质时,必须按能带依次处理)
LOPTICS=.TRUE.(计算光学性质)
NSW=1
IBRION=-1(无需进行构型优化)
计算结束后(正常),会得到OPTIC(用于计算线性光学性质)和momentum_matrix
(用于计算NLO即非线性光学性质)
3.编辑OPTCRT文件常用关键词及其含义
ISYMM=2
OMMIN=0
OMMAX =20
NEDOS =4000
NBCON =200
LJDOS= .TRUE.
LDOS =.TRUE.
LKRAMERS =.TURE.
GAMMA=0.002
LSEARCH=.TRUE.
EMINSEARCH = 0
EMAXSEARCH=20
SCISSOR=0.5
4.创建一个新目录,将OPTICOPTCTR,KPOINTS,POSCAR复制到其下
命令cp../ OPTIC. OPTCTR,KPOINTS,POSCAR.
修改KPOINTS
创建个临时目录复制vasp计算所需要的四个输入文件.运行单机版vasp。
(由于56机子没有装单机版,所以这一步COPY到其他的机子上做的)
单机运行命令
~/bin/vasp.4.5-mk-sp-pgi<
/dev/null>vasp.out&
将其中的IBZKPT复制为KPOINTS
运行vasp
命令:
mpiexec –np 8 ~/bin/vasp.4.5-mk-mp-pgi<
/dev/null>
vasp.out&
正常运行将产生2个文件
EPS----存放介电函数数据
JDOS----存放JDOS数据
在运行vasp_lo
用origin软件根据EPS数据绘图
在此基础上才可以得到其他线性光学性质:
(1)折射率(refractiveindex)
(2)吸收系数(adsorption constant)
(3)能量损失系数(energy-loss coefficient) :
(4)消光系数(extinction coefficient)
(5)双折射率曲线
(6)光导电变化曲线
(二)非线性光学性质
1.非线的运行需要以下三个文件
KPOINTS、EIGENVAL和momentum_matrix三个文件,
其中KPOINTS文件内容与计算线性光学性质相同.
EIGENVAL文件为vasp计算输出文件
运行vasp_nlo, 需输入以下数据:
momentum_matrix 存放跃迁矩阵元的文件名
177.32 单胞体积
2.76 实验带隙(若为0则无需对带隙进行校正)
1 起始能带(一般均为1)
38 终止能带
F 为False时表示只计算静态倍频系数
2.运行后会得到静态倍频系数的计算结果.(此数据好像没有自己保存)
d(1, 1)SHG (realpart)= 1.60120and(imagpart) 0.00000(pm/V)
d(1, 2)SHG (realpart)= -2.37106 and(imag part) 0.00000(pm/V)
d(1, 3)SHG(real part)= -2.72264and (imagpart) 0.00000(pm/V)
d(1,4) SHG (realpart)= 10.98424and (imagpart) 0.00000(pm/V)
d(1,5)SHG (realpart)= -0.00510and (imagpart) 0.00000(pm/V)
d(1,6)SHG(realpart)= -2.46438 and(imagpart) 0.00000(pm/V)
d(2,1) SHG(realpart)= -2.46438and (imagpart) 0.00000(pm/V)
d(2,2) SHG(realpart)= 0.34388and(imagpart) 0.00000(pm/V)
d( 2,3)SHG (realpart)= -1.40097and(imagpart) 0.00000(pm/V)
d(2,4) SHG(realpart)= -0.42303and (imag part) 0.00000(pm/V)
d(2, 5)SHG(real part)= 10.98424 and(imagpart) 0.00000(pm/V)
d(2,6)SHG(realpart)= -2.37106and(imagpart)0.00000(pm/V)
d(3,1)SHG(real part)= -0.00510and(imagpart) 0.00000(pm/V)
d( 3,2) SHG(realpart)= -0.42303and(imagpart) 0.00000(pm/V)
d(3,3)SHG (realpart)= 1.21340and(imagpart) 0.00000(pm/V)
d(3,4)SHG(realpart)= -1.40097and(imagpart) 0.00000(pm/V)
d( 3,5)SHG (real part)= -2.72264and(imagpart) 0.00000(pm/V)
d(3, 6) SHG(realpart)= 10.98424and(imag part) 0.00000(pm/V)
由上可得最大倍频系数为d(1.4)
动态倍频系数计算 运行vasp
需要输入以下数据
momentum_matrix存放跃迁矩阵元的文件名
42.34 单胞体积
2.26 实验带隙
1 起始能带
26 终止能带
T 为True时表示要计算动态倍频系数
10 考察动态倍频系数时的能量范围
500 格点数
14 倍频系数分量d14
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- VASP 计算 AgGaS2 能带 密度 光学 性质