vasp计算介电张量的结果文件

体系微观介电常数epsilon(r,r')是一个二变量的函数,对其作傅立叶变换后,傅立叶系数可以写成一个矩阵的形式epsilon( G, G' ),G和G'是倒格矢。宏观介电常数由微观介电常数在原胞内积分取均值获得,可以证明最终结果等于其傅立叶变换系数的(0,0)项,也就是epsilon(0,0)。

但是一般程序在计算时算的不是epsilon(G,G'),而是epsilon^-1(G,G')。这样必须先取整个epsilon^-1(G,G')矩阵的逆,再取其(0,0)元,这样做时就包含了局域场效应。

如果为了计算简便,仅对epsilon^-1(G,G')的对角元取逆,就忽略了局域场效应。
一般来说,局域场效应对inhomogeneous系统影响较大。

4楼: Originally posted by liqizuiyang at 2013-11-05 14:21:04体系微观介电常数epsilon(r,r')是一个二变量的函数,对其作傅立叶变换后,傅立叶系数可以写成一个矩阵的形式epsilon( G, G' ),G和G'是倒格矢。
宏观介电常数由微观介电常数在原胞内积分取均值获得,可以证明最 ...

这位专家,你点应助啊,我好给你金币!

回答的很专业,学到很多,但是我还是对局域场效应不懂,你能通俗的解释一下这个效应吗?是不是大多体系在计算介电性质的时候都应该考虑这种效应才对?多谢指导!

6楼: Originally posted by future_wl at 2013-11-06 09:11:12回答的很专业,学到很多,但是我还是对局域场效应不懂,你能通俗的解释一下这个效应吗?是不是大多体系在计算介电性质的时候都应该考虑这种效应才对?
多谢指导!...

这个效应就是由于电子分布不均一而对介电函数产生了影响,导致通常采用的一种简便计算方法不再适用。一般来说,电子分布不均一的系统需要考虑。如果计算能力许可,最好还是选考虑局域场效应的计算结果。

7楼: Originally posted by liqizuiyang at 2013-11-06 09:35:56这个效应就是由于电子分布不均一而对介电函数产生了影响,导致通常采用的一种简便计算方法不再适用。
一般来说,电子分布不均一的系统需要考虑。如果计算能力许可,最好还是选考虑局域场效应的计算结果。...

哦,那局域场效应同随机相近似(Random Phase Approximation)又有什么联系吗? 在vasp计算中如果不考虑RPA那只是在DFT层面上计算局域场效应,考虑了RPA就会在Hartree层面上计算局域场效应,我这种理解对不对呢?

8楼: Originally posted by future_wl at 2013-11-06 20:31:06哦,那局域场效应同随机相近似(Random Phase Approximation)又有什么联系吗? 在vasp计算中如果不考虑RPA那只是在DFT层面上计算局域场效应,考虑了RPA就会在Hartree层面上计算局域场效应,我这种理解对不对呢?...

RPA是一种计算节点函数的理论,局域场效应是计算介电函数过程中所涉及的一种偏数学的近似手段。这里VASP手册里写得我也看得不是太明白,按我理解应该是这样的:
算介电函数时先算的是极化函数X,然后由X计算eps^-1。

由前面DFT计算得到的波函数和本征值可以根据RPA算出X0,更精确的应该是由X0算出的X。考不考虑RPA的区别主要在X上。
你的理解应该是对的。
推荐去看下DP和exciting这两个TDDFT程序的主页,上面有讲解TDDFT理论的PPT,里面对极化率节点函数讲得很清楚。

:sweat:抱歉,上面打错了。应该是DP和yambo。

10楼: Originally posted by liqizuiyang at 2013-11-07 09:46:59:sweat:抱歉,上面打错了。应该是DP和yambo。

你好,还想请教一下,如果设置LEPSILON=.TRUE和LRPA=.F,结果文件还会给出电子和离子两部分对介电常数的贡献,MACROSCOPIC STATIC DIELECTRIC TENSOR (including local field effects in DFT)
MACROSCOPIC STATIC DIELECTRIC TENSOR IONIC CONTRIBUTION
如果要和实验上的静态介电常数去比较的话,是不是应该将这两项加起来?
多谢:hand: