VASP计算AgGaS2能带及态密度及光学性质.docx
- 文档编号:4181181
- 上传时间:2022-11-28
- 格式:DOCX
- 页数:16
- 大小:189.26KB
VASP计算AgGaS2能带及态密度及光学性质.docx
《VASP计算AgGaS2能带及态密度及光学性质.docx》由会员分享,可在线阅读,更多相关《VASP计算AgGaS2能带及态密度及光学性质.docx(16页珍藏版)》请在冰豆网上搜索。
VASP计算AgGaS2能带及态密度及光学性质
VASP计算AgGaS2能带及态密度及光学性质
第一步:
构型优化
1.准备四个输入文件
POSCARINCARPOTCARKPOINTS
POSCAR:
从ms中导入AgGaS2结构,选择CASTEP,file,save,并保存成原包。
这样,得到一隐藏文件。
cell,将它用编辑器打开,从中的到vasp所需的POSCAR信息,修改得到POSCAR。
AgGaS2bulk
1.000000000000000
-2。
79349994659424102。
79349994659423705。
2045001983642580
2。
7934999465942370—2。
79349994659423905.2045001983642590
2。
79349994659423802。
7934999465942380-5.2045001983642580
422
Direct
0。
37500000000000000.40720000863075260.5322000086307526
0.87500000000000000。
84279999136924740。
4677999913692474
。
.。
。
。
.。
.。
。
。
。
2.POTCAR:
用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=。
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-pgivasp.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 就可以得到坐标 手动复制到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 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 1 5 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 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单胞体积 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 能带 密度 光学 性质