vasp计算压电系数_求助DFTP算出来的压电系数 - 第一原理 - 小木虫 - 学术 科研 互动社区...

我用的VASP5.2 下的LEPSILON-tag算的压电系数和BEC等东西，之前发过类似的贴，不过这次想一起贴上来在好好求助一下：

我的INCAR

SYSTEM=PLPK

ENCUT=400

PREC=Normal

ISMEAR=0

SIGMA=0.05

EDIFF=1E-3

ISTART=0

ICHARG=2

IBRION=8

POTIM=0.1

LEPSILON=.TRUE.

LDAU=.TRUE.

LDAUTYPE=2

LDAUL=-1 -1 3 -1 -1 -1 -1 -1 -1

LDAUU=0 0 6 0 0 0 0 0 0

LDAUJ=0 0 0 0 0 0 0 0 0

LDAUPRINT=2

算的是一个20个原子的四方相结构，含极性不连续的。主要算它的极化强度和压电特性在应变条件下的变化，现在只是OUTCAR中的结果还差那么点儿分析，想多请教下算过相关内容的前辈和同行，希望不吝赐教。

比如0应变的结构，我从OUTCAR中找到的有用内容如下：(目前提取出来了弹性模量，BEC和压电张量矩阵，还有个介电张量)

MACROSCOPIC STATIC DIELECTRIC TENSOR (including local field effects in DFT)

---------------------------------------------

7.562     0.000     0.00