视频1 视频21 视频41 视频61 视频文章1 视频文章21 视频文章41 视频文章61 推荐1 推荐3 推荐5 推荐7 推荐9 推荐11 推荐13 推荐15 推荐17 推荐19 推荐21 推荐23 推荐25 推荐27 推荐29 推荐31 推荐33 推荐35 推荐37 推荐39 推荐41 推荐43 推荐45 推荐47 推荐49 关键词1 关键词101 关键词201 关键词301 关键词401 关键词501 关键词601 关键词701 关键词801 关键词901 关键词1001 关键词1101 关键词1201 关键词1301 关键词1401 关键词1501 关键词1601 关键词1701 关键词1801 关键词1901 视频扩展1 视频扩展6 视频扩展11 视频扩展16 文章1 文章201 文章401 文章601 文章801 文章1001 资讯1 资讯501 资讯1001 资讯1501 标签1 标签501 标签1001 关键词1 关键词501 关键词1001 关键词1501 专题2001
vasp做光学性质计算
2025-09-24 13:55:20 责编:小OO
文档
关于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=17492&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.4419    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.0190    0.0190    0.0190    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

下载本文
显示全文
专题