| 查看: 36606 | 回复: 106 | ||
| 本帖产生 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» 猜你喜欢
可靠消息，8月21号(明天)放榜
已经有3人回复
-
只有每年这种时候来逛逛小木虫
已经有20人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
99楼2016-09-29 13:24:46
3楼2011-10-28 20:42:33
5楼2011-10-28 20:46:04
6楼2011-10-28 21:57:39
