24小时热门版块排行榜    

Znn3bq.jpeg
查看: 34935  |  回复: 106
【奖励】 本帖被评价33次,作者liliangfang增加金币 22.6
本帖产生 1 个 1ST强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] vasp5.2计算半导体Si光学性质的总结

前一段时间计算了Si的光学性质,现在整理下帖出来,希望对光学方面的虫友有帮助,当然有不当的地方请大家指教。由于DFT方法计算出的半导体带隙不太准确,所以这里在用HSE杂化计算能带后直接计算光学性质
就顺水提出的问题,做下面说明,标志红色
1.建模结构优化这里就不说了
2..优化后DFT计算(HSE杂化需要DFT的支持)这里需要说明的是后面计算要用到这步的波函数
INCAR:
System = diamond Si
ENCUT = 400
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.1
PREC = Accurate
EDIFF = 1.0E-05
EDIFFG = -0.001
NSW = 0
IBRION = -1
ISIF = 2
GGA = PS(此处如果用GGA = PS请在POTCAR中修改LEXCH  = PS)
KPOINTS:
Automatic mesh
0
M
7   7   7
0   0   0

POSCAR:
Si-Diamond:                             
   5.43800000000000     
     0.0000000000000000    0.5000000000000000    0.5000000000000000
     0.5000000000000000    0.0000000000000000    0.5000000000000000
     0.5000000000000000    0.5000000000000000    0.0000000000000000
   Si
   2
Selective dynamics
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000   T   T   T
  0.2500000000000000  0.2500000000000000  0.2500000000000000   T   T   T

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
运行vasp
3.HSE06计算(此步在DFT基础上)这步用到上步的波函数

INCAR:
SYSTEM = Si-Diamond
ISTART = 1
ICHARG = 2
EDIFF = 0.00001
EDIFFG = -0.001
GGA = PS
ISMEAR = 0
SIGMA = 0.1
ENCUT = 400
ENAUG = 800
LREAL = .FALSE.
LWAVE = .TRUE.
LCHARG = .TRUE.
NELM = 200
NSW = 0
IBRION = -1
LMAXMIX = 4
NSIM = 4
ISYM = 0
LHFCALC = .TRUE.
HFSCREEN = 0.2
ALGO = Damped
TIME = 0.4
ENCUTFOCK = 0
AEXX = 0.25
NBANDS = 24(此处为DFT计算是的2—3倍)根据说明要有足够的空带
其他文件从DFT文件夹拷贝过来,运行vasp计算

4.optic计算(这步用到上面的波函数)
INCAR:
SYSTEM = Si-Diamond
ISTART = 1
ICHARG = 2
EDIFF = 0.00001
EDIFFG = -0.01
GGA = PS
ENCUT = 400
ENAUG = 800
NELM = 200
NSW = 0
IBRION = 1
LMAXMIX = 4
ISMEAR = -5
NSIM = 4
ISYM = 0
LHFCALC = .TRUE.
HFSCREEN = 0.2
ALGO = N
TIME = 0.4
ENCUTFOCK = 0
AEXX = 0.25
NELMIN = 5
LOPTICS = TRUE
CSHIFT = 0.1
NBANDS = 24
NPAR = 1
NEDOS = 2000
LRPA = .FALSE.              
LREAL = .FALSE.
LWAVE = .FALSE.   
LCHARG = .FALSE.

KPOINTS:(此处为HSE计算时的2-3倍,关于这一点加密K点是为保证介电谱收敛)
Automatic mesh
0
M
9   9   9
0   0   0

其他全部文件从HSE06计算文件夹中拷贝
主要此处NPAR = 1,vasp运行时不可并行运行这里要说的是NPAR = 1按能带依次处理。如果此处还用并行那么vasp会要求你设置NPAR = number ,与NPAR = 1不合
运行vasp

5.结果(与实验值做比较)
介电函数实部:

介电函数虚部:

介电函数:

折射率:

还有差别,欢迎大家指正
这里感谢顺水提出的问题,如大家还有问题,请指出并指出问题的关键及解决方法,不要让大家都感到很迷惑谢谢

[ Last edited by liliangfang on 2012-4-23 at 11:25 ]
回复此楼

» 收录本帖的淘帖专辑推荐

资源收集 仿真建模与计算 VASP and MS 材料计算模拟实用技巧
第一性原理相关文档 第一性原理计算-tg 光学工程新思想 VASP软件从入门到高深
第一原理资料贴 第一性原理计算经验 表界面计算 VASP第一性原理的模拟计算
VASP资料 第一性原理 VASP算例 vasp光学
first principle HSE & 杂化函数 计算化学-光学性质 计算资料
@个人收集@ 量化软件 VASP计算资料收集byZF 博士发论文的思考体会
杂化泛函方面的研究 第一原理 VASP相关 纳米及第一性原理方向
网贴收集 模拟,理论 科研 vasp
vasp good resources 光学性质计算 化工专业课资料
vasp vasp-光学 vasp讨论帖子 SIESTA学习
VASP soft study 光学 半导体

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)


★★★★★ 五星级,优秀推荐

引用回帖:
6楼: Originally posted by liliangfang at 2011-10-28 21:57:39:
首先感谢指出错误,我是使用5.2的说明书的,我们的计算机一共是128个核的。还有啊请指出错误的地方,谢谢

受内伤了用四核并行算任务的孩子伤不起啊。
楼主很细心,总结的也很好,向你学习!加油
28楼2011-11-02 21:19:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 107 个回答

顺水

金虫 (小有名气)


★ ★
youzhizhe(金币+2): 谢谢交流。 2011-10-30 10:18:52
你这有很多是错的。
建议要重新核实一下下。
看看vasp的说明书。
5.2就是5.2不是4.6。
3楼2011-10-28 20:42:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

顺水

金虫 (小有名气)


★ ★
youzhizhe(金币+2): 谢谢交流。 2011-10-30 10:18:58
不该加的K点不能加啊!
不然会万劫不复,特别是你的计算机计算能力弱的时候!
5楼2011-10-28 20:46:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youzhizhe(金币+1): 谢谢交流。 2011-10-30 10:22:17
引用回帖:
5楼: Originally posted by 顺水 at 2011-10-28 20:46:04:
不该加的K点不能加啊!
不然会万劫不复,特别是你的计算机计算能力弱的时候!

首先感谢指出错误,我是使用5.2的说明书的,我们的计算机一共是128个核的。还有啊请指出错误的地方,谢谢
6楼2011-10-28 21:57:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +11 汪!?! 2026-04-10 13/650 2026-04-10 23:16 by Ftglcn90
[考研] 080500求调剂 +16 黄宇博 2026-04-06 16/800 2026-04-10 22:33 by Ftglcn90
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 298求调剂 +8 钉叮咚冬瓜 2026-04-07 10/500 2026-04-10 11:25 by 高维春
[考研] 求调剂 +16 熊二想上岸 2026-04-04 16/800 2026-04-10 11:24 by may_新宇
[考研] 求调剂 +11 翩翩一书生 2026-04-09 13/650 2026-04-10 10:27 by liuhuiying09
[考研] 22408 352分求调剂0854类 +3 努力的夏末 2026-04-09 3/150 2026-04-10 10:12 by 314126402
[考研] 0854调剂 +7 950824he@ 2026-04-09 7/350 2026-04-10 09:10 by Delta2012
[考研] 085500求调剂材料 +8 易11122 2026-04-09 8/400 2026-04-09 23:15 by parmtree
[考研] 274求调剂 +5 山阿蔓 2026-04-07 5/250 2026-04-09 15:28 by 18828373951
[考研] 291求调剂 +7 关忆北. 2026-04-09 8/400 2026-04-09 15:17 by 探123
[考研] 086000调剂 +4 十七sa 2026-04-07 4/200 2026-04-09 09:27 by 壹往無前
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 求调剂 +9 月@163.com 2026-04-07 11/550 2026-04-08 14:48 by qlm5820
[考研] 315求调剂 +17 欣喜777 2026-04-04 18/900 2026-04-08 13:54 by hangsimei
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 一志愿211生物学280分 求调剂 +5 李rien 2026-04-05 5/250 2026-04-06 10:30 by zhyzzh
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 能动调剂326专硕 +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
信息提示
请填处理意见