| 查看: 36496 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
» 本帖已获得的红花（最新10朵）
donggaomin» 猜你喜欢
时间戳变了，能看出什么问题？
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有80人回复
-
感觉是下周放榜了
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp5.2光学性质计算的疑问
已经有5人回复
- vasp5.2如何计算DOS和band structure？
已经有10人回复
- VASP 5.2.2INCAR中的LEPSILON问题
已经有3人回复
- VASP中如何创建 INCAR , POSCAR, KPOINTS和POTCAR文件
已经有7人回复
- GW计算的错误提示【问题已解决】
已经有4人回复
- vasp5.2计算光学性质
已经有22人回复
- 半导体中的基本性质（课件 北京大学 微电子学研究所）【已查重】
已经有505人回复
- vasp5.2光学计算中怎么加入scissor operation?
已经有13人回复
- vasp5.2计算误差太大
已经有11人回复
- 介电函数的虚部与能带结构
已经有5人回复
- 【求助】VASP计算光学性质时的GW方法（RPA, no local field effects)
已经有5人回复
- 【讨论】VASP 光学性质计算结果和实验图的比较，以及一个问题
已经有4人回复
- 【求助】vasp5 光学计算结果问题
已经有4人回复
- 【求助】vasp 5.2中用HSE06计算不能进入主循环
已经有7人回复
- 【求助】vasp5.2光学性质问题
已经有18人回复
- 【求助】GW计算错误
已经有6人回复
25