vasp做光学性质计算

关于vasp5.2计算光学性质的问题

作者: ljx2006 收录日期: 2010-05-20 发布日期: 2010-05-10

两个问题：

1。做单点计算直接就可以得到电子的能带信息，态密度，和光学性质？ 我们往常不是做完静态之后再算能带，态密度的吗？单点和静态有什么联系和区别吗？

2。 用vasp5.2如何计算材料的光学性质。我的计算INCAR文件如下： System=1-4-optical-1 ALGO=Fast

LOPTICS=.TRUE. ISTART=1 ISMEAR=-5 SIGMA=0.1

IBRION=-1 EDIFF=1.0E-05 EDIFFG=-0.02 NSW=0 ENCUT=400 NPAR=1 CSHIFT=0.1

LNABLA=.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://emuch.net/bbs/viewthread.php?tid=1953056&fpage=1

【重点讨论】关于vasp5.2计算光学性质的问题

http://emuch.net/bbs/viewthread.php?tid=2037690&fpage=1 About MS

1、【重点讨论帖】castep 如何设置FM和AFM

http://emuch.net/bbs/viewthread.php?tid=1703336&fpage=1 2、【重点讨论】Castep中的Empty Band的用处和设置

http://emuch.net/bbs/viewthread.php?tid=1744052&fpage=1 3、【重点讨论】一个氧空位表面charge是不是要设为2

http://emuch.net/bbs/viewthread.php?tid=1749642&fpage=1 4、【重点讨论帖】如何将两种结构整合在一起

http://emuch.net/bbs/viewthread.php?tid=1849550&fpage=1 5、【重点讨论】掺杂原子化合价怎么确定

http://emuch.net/bbs/viewthread.php?tid=1942922&fpage=1

作者：cbFeng200881 OUTCAR

作者：ljx2006

看到你的帖子也学习了，希望用VASP计算光学性质的虫友也一起参与讨论，互相学习。

作者：ljx2006

强烈期待用vasp计算过光学性质的或者对此有见解的虫友发表看法，共同讨论，互相学习！

作者：vasp001

介电常数，可以在VASP的OUTCAR 中得到，包含两部分，实部和虚部 都是一个矩阵的形式。

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.003575 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.013067 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.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文件里也看到类似的信息了，就是你提到的介电常数，但是我不清楚的怎么进一步得到折射率，反射率，还有吸收系数与能量的谱图，所以还请多指点，希望和你再交流！

作者：vasp001 QUOTE:

Originally posted by ljx2006 at 2010-05-13 11:37:43:

你好，我的INCAR文件是在得到静态电荷密度和波函数之后开始调用光学性质关键词的，我自己就是看手册算的，所以不是很肯定。

但是，我在我的OUTCAR文件里也看到类似的信息了，就是你提到的介电常数，但是我不清楚 ...

我看到他们有人用这种方法做过，求的是频率和吸收强度的关系。

至于反射率和折射率和介电常数的关系，你需要看一下光学课本，以前学过，都忘记了。

作者：vasp001

补充一点资料，一篇SCIENCE 的资料：

www.sciencemag.org SCIENCE VOL 293 13 JULY 2001

关于对这个矩阵的处理，看来是将每行的x y z 取平均了，其他的三个分量为零，因此不用考虑，下面是他画的图像，将能量通过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

本文来源：https://www.bwwdw.com/article/gbbh.html