24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3230  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

slow_winds

新虫 (初入文坛)

[交流] VASP计算压电张量不对称

请教各位一下,在下刚刚转来用VASP做钙钛矿结构压电材料的第一性原理计算,但是现在计算中出现了一些问题,我计算P4mm的四方相结构时,压电离子贡献24与15方向不相等,导致计算出来的e24也不等于e15,请教下是什么原因。
POSCAR
KNbO3                                   
   1.00000000000000     
     3.9780718153912340    0.0000000000000000    0.0000000000000000
     0.0000000000000000    3.9780718153912340    0.0000000000000000
     0.0000000000000000    0.0000000000000000    4.0697524641330913
   K    Nb   O
   1   1   3
Direct
  0.0000000000000000  0.0000000000000000  0.0208122828936224
  0.5000000000000000  0.5000000000000000  0.5370745926342191
  0.5000000000000000  0.5000000000000000  0.9956543127557665
  0.0000000000000000  0.5000000000000000  0.4982294058581956
  0.5000000000000000  0.0000000000000000  0.4982294058581956

INCAR
SYSTEM = KNbO3

  ISTART =  0
  ICHARG =  2
  PREC   =  Accurate
  EDIFF  =  1E-06
  EDIFFG =  -0.005
  NSW    =  1
  POTIM  =  0.2
  NFREE  =  2
  IBRION =  8
  LEPSILON = Ture
  ISIF   =  3
  ISPIN  =  1
  NELM   =  60
   ISMEAR = 0
   SIGMA  = 0.2
   LREAL = Auto
   LCHARGE = False
   LWAVE = False

计算结果
PIEZOELECTRIC TENSOR  for field in x, y, z        (C/m^2)
          XX          YY          ZZ          XY          YZ          ZX   
  --------------------------------------------------------------------------------
  x     0.00000     0.00000     0.00000     0.00000     0.00000     0.19952
  y     0.00000     0.00000     0.00000     0.00000     0.19957     0.00000
  z     0.16714     0.16714     0.06456     0.00000     0.00000     0.00000

PIEZOELECTRIC TENSOR IONIC CONTR  for field in x, y, z        (C/m^2)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x    -0.00020     0.00055    -0.00176     0.00000     0.00000    -5.47889
  y    -0.00076     0.00097    -0.00008     0.00000     1.73507     0.00000
  z    -0.46524    -0.46524     4.83952     0.00000     0.00000     0.00000

d33值利用d=e*s计算出来只有20pC/N多点。烦请各位告知哪里有问题
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

slow_winds

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by tsglss at 2018-02-20 21:03:38
在你的结果中在哪里能看出e24和e15呢?另外,d33和d等于es这个公式是什么意思?谢谢

xZX是15,yYZ是24,离子贡献上一个-5.4,一个1.7.
d=e*s,d和e是3x6矩阵,s是弹性柔顺系数,是弹性劲度系数c的6x6逆矩阵。
5楼2018-02-21 01:41:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

tsglss

铁杆木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
请问,为什么要e24和e15相等呢?根据哪个文献?

[ 发自手机版 http://muchong.com/3g ]
做好自己
2楼2018-02-20 15:40:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsglss

铁杆木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
在你的结果中在哪里能看出e24和e15呢?另外,d33和d等于es这个公式是什么意思?谢谢

[ 发自手机版 http://muchong.com/3g ]
做好自己
3楼2018-02-20 21:03:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

slow_winds

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by tsglss at 2018-02-20 15:40:41
请问,为什么要e24和e15相等呢?根据哪个文献?

4mm点群的x,y是相同的,15和24分量可以通过对称变换得出相等。
4楼2018-02-21 01:36:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见