VASP计算AgGaS2能带及态密度及光学性质Word格式.docx

VASP计算AgGaS2能带及态密度及光学性质Word格式.docx

《VASP计算AgGaS2能带及态密度及光学性质Word格式.docx》由会员分享，可在线阅读，更多相关《VASP计算AgGaS2能带及态密度及光学性质Word格式.docx（16页珍藏版）》请在冰豆网上搜索。

用PBE-GGA的赝势，提取，AgGaO的赝势合并成一个赝势。

（一般VASP有自带）

3。

下面是INCAR

SYSTEM=optimizationofAgGaS2

LPLANE=。

TRUE.

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=。

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文件

mpiexec—np8~/bin/vasp。

4.5-mk-mp—pgi〈/dev/null〉vasp。

out&

（运行命令）

第三部：

能带的计算。

建立新的计算目录

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

mpiexec—np8～/bin/vasp。

4.5—mk—mp-pgi〈/dev/null>

vasp。

out＆（运行命令）

在运行

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=。

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

LCHARG=.TRUE。

LVTOT=。

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即非线性光学性质）

编辑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

创建一个新目录，将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<

out＆

正常运行将产生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单胞体积

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。

d（1，6）SHG（realpart）=-2。

46438and（imagpart）0.00000（pm/V）

d（2,1）SHG（realpart）=-2。

46438and（imagpart）0。

d（2，2）SHG（realpart）=0.34388and（imagpart）0。

d（2，3）SHG（realpart）=—1。

40097and（imagpart）0。

d（2,4）SHG（realpart）=-0。

42303and（imagpart）0。

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。

d（3，3）SHG（realpart）=1.21340and（imagpart）0.00000（pm/V）

d（3，4）SHG（realpart）=—1。

d（3，5）SHG（realpart）=-2.72264and（imagpart）0.00000（pm/V）

d（3,6）SHG（realpart）=10.98424and（imagpart）0。

由上可得最大倍频系数为d（1。

4）

动态倍频系数计算运行vasp

需要输入以下数据

momentum_matrix存放跃迁矩阵元的文件名

42。

34单胞体积

2.26实验带隙

1起始能带

26终止能带

T为True时表示要计算动态倍频系数

10考察动态倍频系数时的能量范围

500格点数

14倍频系数分量d14