Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 34993  |  回复: 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的回帖
引用回帖:
37楼: Originally posted by 小豆sumy at 2011-11-18 15:00:17:
你好,得到以下数值后,作图是用X下面的值还是Y呢?你能不能解释一下X、Y、Z、XY、YZ、ZX分别代表什么呢?谢谢你!

  frequency dependent IMAGINARY DIELECTRIC FUNCTION (RPA, no local field effects)
    ...

http://muchong.com/bbs/viewthread.php?tid=3561053
38楼2011-11-18 15:26:02
已阅   回复此楼   关注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的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408 312求调剂 +16 门路摸摸 2026-04-14 17/850 2026-04-14 23:59 by Xurambo2014
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 273求调剂 +4 白居不易. 2026-04-09 6/300 2026-04-14 19:54 by 逆水乘风
[考研] 药学305求调剂 +6 玛卡巴卡boom 2026-04-11 6/300 2026-04-14 19:48 by zhouxiaoyu
[考研] 材料专业344求调剂 +17 hualkop 2026-04-10 22/1100 2026-04-14 16:21 by sxdj2
[考研] 本科211,报考085601-310分 +16 ararak 2026-04-13 16/800 2026-04-14 14:55 by Delta2012
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 085404 298分求调剂 +11 呼啦呼啦呼呼呼 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[考研] 材料复试求调剂 +24 xhhdjdjsjks 2026-04-09 24/1200 2026-04-13 15:49 by 幸免 ..
[考研] 一志愿211 0703化学 346分求调剂 +26 土豆er? 2026-04-09 29/1450 2026-04-13 15:15 by 独醉梦孤城
[考研] 生物学调剂 +11 小冉要努力 2026-04-10 13/650 2026-04-13 11:46 by 电化学及催化
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 一志愿京区985,085401,与本科专业一致,电子信息工程, +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
信息提示
请填处理意见