VASP计算AgGaS2能带及态密度及光学性质Word文档下载推荐.docx

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

欢迎下载，资料仅供参考!

!

资料仅供参考!