《v做光学性质计算.doc》由会员分享,可在线阅读,更多相关《v做光学性质计算.doc(5页珍藏版)》请在金锄头文库上搜索。
1、关于vasp5.2计算光学性质的问题作者: ljx2006 收录日期: 2010-05-20 发布日期: 2010-05-10两个问题:1。做单点计算直接就可以得到电子的能带信息,态密度,和光学性质?我们往常不是做完静态之后再算能带,态密度的吗?单点和静态有什么联系和区别吗?2。 用vasp5.2如何计算材料的光学性质。我的计算INCAR文件如下:System=1-4-optical-1ALGO=Fast LOPTICS=.TRUE.ISTART=1 ISMEAR=-5 SIGMA=0.1 IBRION=-1 EDIFF=1.0E-05 EDIFFG=-0.02NSW=0 ENCUT=400
2、NPAR=1 CSHIFT=0.1LNABLA=.FALSE.LEPSILON=.FALSE.LRPA=.TRUE. VOSKOWN=1 LREAL=Auto LWAVE=.FALSE. LCHARG=.FALSE. 输出文件中有个OPTIC文件,但是怎么都打不开,请高手指点,应该要怎么做输入文件才能的到材料的光学性质:介电常数,吸收系数。折射和反射率? Last edited by wuli8 on 2010-5-10 at 17:07 作者:wuli8【重点讨论帖】vasp处理f电子体系http:/ http:/ MS1、【重点讨论帖】castep 如何设置FM和AFM http:/ Ba
3、nd的用处和设置http:/ 中得到,包含两部分,实部和虚部 都是一个矩阵的形式。frequency dependent IMAGINARY DIELECTRIC FUNCTION (RPA, no local field effects) E(ev) X Y Z XY YZ ZX- 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.092884 0.001782 0.001782 0.001782 0.000000 0.000000 0.000000 0.185768 0.003575 0.003575 0.00
4、3575 0.000000 0.000000 0.000000 0.278652 0.005388 0.005388 0.005388 0.000000 0.000000 0.000000 0.371535 0.007232 0.007232 0.007232 0.000000 0.000000 0.000000 0.464419 0.009119 0.009119 0.009119 0.000000 0.000000 0.000000 0.557303 0.011059 0.011059 0.011059 0.000000 0.000000 0.000000 0.650187 0.01306
5、7 0.013067 0.013067 0.000000 0.000000 0.000000 0.743071 0.015156 0.015156 0.015156 0.000000 0.000000 0.000000 0.835955 0.017341 0.017341 0.017341 0.000000 0.000000 0.000000 0.928839 0.019640 0.019640 0.019640 0.000000 0.000000 0.000000 1.021722 0.022073 0.022073 0.022073 0.000000 0.000000 0.000000 1
6、.114606 0.024661 0.024661 0.024661 0.000000 0.000000 0.000000 1.207490 0.027431 0.027431 0.027431 0.000000 0.000000 0.000000作者:ljx2006你好,我的INCAR文件是在得到静态电荷密度和波函数之后开始调用光学性质关键词的,我自己就是看手册算的,所以不是很肯定。但是,我在我的OUTCAR文件里也看到类似的信息了,就是你提到的介电常数,但是我不清楚的怎么进一步得到折射率,反射率,还有吸收系数与能量的谱图,所以还请多指点,希望和你再交流!作者:vasp001QUOTE:Or
7、iginally posted by ljx2006 at 2010-05-13 11:37:43:你好,我的INCAR文件是在得到静态电荷密度和波函数之后开始调用光学性质关键词的,我自己就是看手册算的,所以不是很肯定。但是,我在我的OUTCAR文件里也看到类似的信息了,就是你提到的介电常数,但是我不清楚 . 我看到他们有人用这种方法做过,求的是频率和吸收强度的关系。至于反射率和折射率和介电常数的关系,你需要看一下光学课本,以前学过,都忘记了。作者:vasp001补充一点资料,一篇SCIENCE 的资料:www.sciencemag.org SCIENCE VOL 293 13 JULY 20
8、01关于对这个矩阵的处理,看来是将每行的xyz 取平均了,其他的三个分量为零,因此不用考虑,下面是他画的图像,将能量通过E=hc/a ,c 为光速,a为光的波长,这样可以将能量转换成波长,可以做图了。这是这篇science里的图像。希望有用。照片传的不好看,重新上传一个文件: Last edited by vasp001 on 2010-5-13 at 18:01 作者:westmonster为什么要三个方向取平均呢? 如果你是个低维体系的话。作者:hlshi立方结构对称性时,介电张量只有一个独立的分量xx=yy=zz,其他的为零。在做介电谱虚部时只需用xx分量即可;至于吸收谱、能量耗散函数等可以参考文献中的公式,编个小程序即可可以参照http:/arxiv.org/abs/0911.3733