24小时热门版块排行榜    

查看: 34685  |  回复: 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的回帖
★ ★
youzhizhe(金币+2): 谢谢交流。 2011-10-30 10:22:03
引用回帖:
9楼: Originally posted by 顺水 at 2011-10-28 22:40:38:
KPOINTS:(此处为HSE计算时的2-3倍)
为什么?请给出参考文献。
不能任意加K点的。
你用什么算的?
波函数啊!!!!
怎么能够半路加K点。
贴出你的输出信息让大家看看发生了什么。
前 ...

这个我也很怀疑,这是一个印度人告诉我的,不置可否
running on    1 nodes
distr:  one band on    1 nodes,    1 groups
vasp.5.2.2 15Apr09 complex
POSCAR found type information on POSCAR  Si
POSCAR found :  1 types and    2 ions
LDA part: xc-table for Pade appr. of Perdew
found WAVECAR, reading the header
  number of k-points has changed, file:    36 present:   112
  trying to continue reading WAVECAR, but it might fail
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...(           1 )
reading WAVECAR
the WAVECAR file was read sucessfully
charge-density read from file: Si-Diamond                              
entering main loop
       N       E                     dE             d eps       ncg     rms          ort
现在强烈要求贴出你的方法,谢谢
19楼2011-10-28 23:23:09
已阅   回复此楼   关注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的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿重庆大学085700资源与环境专硕,总分308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-21 00:39 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 286分人工智能专业请求调剂愿意跨考! +3 lemonzzn 2026-03-17 4/200 2026-03-20 11:04 by lemonzzn
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
信息提示
请填处理意见