VASP计算AgGaS2能带及态密度及光学性质Word文档下载推荐.docx
- 文档编号:19736447
- 上传时间:2023-01-09
- 格式:DOCX
- 页数:16
- 大小:188.34KB
VASP计算AgGaS2能带及态密度及光学性质Word文档下载推荐.docx
《VASP计算AgGaS2能带及态密度及光学性质Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《VASP计算AgGaS2能带及态密度及光学性质Word文档下载推荐.docx(16页珍藏版)》请在冰豆网上搜索。
NPAR=8
Elecronicminimisation
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
M
444
000至此,四文件已准备好,进行计算,
mpiexec-np8~/bin/vasp.4.5-mk-mp-pgi<
/dev/null>
vasp.out&
(运行命令)
(问题:
计算前对某些参数的测试,比如截断能,晶格参数等,标准时什么?
vasp运行后,不知道结果是否满足要求呢?
)
第二部:
将构型优化后产生的CONTCAR文件拷贝为POSCAR文件
第三部:
能带的计算.建立新的计算目录
mkdirband创建目录
cdband进入目录
cp../INCAR.复制
cp../PO*.
cp../CHG*.
对INCAR文件作如下修改:
ICHARG=11
NBANDS=120与上面构型优化时一致
建立新的KPOINTS手动定义K点
AgGaS2
51
Rec(关键字要对齐第一行)
0.50000.50000.50001.0000
0.45000.45000.45001.0000
......
K点坐标获得途径
布里渊区各特殊点的坐标可利用CASTEP获得
然后编辑
viinp
lin1特殊点
lin2特殊点之间的取点数
lin3第一个K点的坐标
lin4第二个K点的坐标
....(同上..)
....
保存
运行vasp_kpoints<
inp
就可以得到坐标
手动复制到KPOINTS
运行VASP
在运行
Vaspband
DrawthebandstructureofVASP
Twofiles,EIGENVALandOUTCARareneeded
----------------------------------------------------------------
Spin-polarizedcalculation(T/F)?
(是否自旋极化计算)
f
NumberoflinesalongtheBZ:
(布里渊区特殊线数目)
5
Kpointsequenceusedtoseparatelines:
(各线段始终点序号)
11121314151
SettheFermileveltozero(T/F)?
(是否费米设置为零点)
t
可以得到相应的band.dat用origin绘图
然后标注各个特殊点
Fermienergy=1.2182vexp=2.76ev(可能是参数设计太粗糙…)不过放大还算清楚..
(用MS作图K网格相同)
第四部:
DOS的计算及其绘图
其POSCARPOTCARKPOINTS可以相同于第一步优化
设置INCAR
LPLANE=.TRUE.
ISTART=3(从WAVECAR读入波函数)
LREAL=.FALSE.
PREC=Medium
EDIFFG=-0.02
NELMIN=3
ISYM=1
ISIF=3
ISMEAR=0
SIGMA=0.1
NBANDS=120
LORBIT=11(投影到s.p.d)
NEDOS=501(DOS的点数)
ICHARG=11(从CHG读入电荷密度)
EMIN=-20
EMAX=20
POTIM=0.05
运行vasp命令同上
然后编写DOS
F(是否自旋极化计算)
9(轨道数)
6(要计算的DOS)
-1(为负数时,表示所考察的为原子的DOS)
1(原子编号)
-1(为正正数时,表示为原子轨道的DOS)
2(下面都同上)
-1
3
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–np8~/bin/vasp.4.5-mk-mp-pgi<
正常运行将产生2个文件
EPS----存放介电函数数据
JDOS----存放JDOS数据
在运行vasp_lo
用origin软件根据EPS数据绘图
在此基础上才可以得到其他线性光学性质:
(1)折射率(refractiveindex)
(2)吸收系数(adsorptionconstant)
(3)能量损失系数(energy-losscoefficient):
(4)消光系数(extinctioncoefficient)
(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.37106and(imagpart)0.00000(pm/V)
d(1,3)SHG(realpart)=-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.46438and(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(imagpart)0.00000(pm/V)
d(2,5)SHG(realpart)=10.98424and(imagpart)0.00000(pm/V)
d(2,6)SHG(realpart)=-2.37106and(imagpart)0.00000(pm/V)
d(3,1)SHG(realpart)=-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(realpart)=-2.72264and(imagpart)0.00000(pm/V)
d(3,6)SHG(realpart)=10.98424and(imagpart)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 能带 密度 光学 性质