24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2113  |  回复: 14

jk射手座

新虫 (著名写手)

[求助] 请问OUTCAR中的这段玻恩有效电荷的信息怎么分析呀已有1人参与

尝试了下用LCALCEPS算BEC, 一下是我的INCAR:
SYSTEM=LAO
ENCUT=400
PREC=Accurate
ISMEAR=-5
SIGMA=0.02
EDIFF=1E-3
ISTART=0
ICHARG=2
LCALCEPS=.TRUE.

LDAU=.TRUE.
LDAUTYPE=2
LDAUL=3 -1 -1
LDAUU=6 0 0
LDAUJ=0 0 0
LDAUPRINT=2

然后这个是OUTCAR中关于玻恩的一段信息,还有压电张量的,
BORN EFFECTIVE CHARGES (including local field effects)
------------------------------------------------------
ion    1
    1     4.52078     0.00000     0.00000
    2     0.00000     4.52078     0.00000
    3     0.00000     0.00000     4.52078
ion    2
    1     2.92932     0.00000     0.00000
    2     0.00000     2.92932     0.00000
    3     0.00000     0.00000     2.92932
ion    3
    1    -2.43128     0.00000     0.00000
    2     0.00000    -2.44923     0.00000
    3     0.00000     0.00000    -2.43128
ion    4
    1    -2.44923     0.00000     0.00000
    2     0.00000    -2.43128     0.00000
    3     0.00000     0.00000    -2.43128
ion    5
    1    -2.43128     0.00000     0.00000
    2     0.00000    -2.43128     0.00000
    3     0.00000     0.00000    -2.44923

PIEZOELECTRIC TENSOR (including local field effects) (e Angst)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x     0.00000     0.00000     0.00000     0.00000     0.00000     0.00000
  y     0.00000     0.00000     0.00000     0.00000     0.00000     0.00000
  z     0.00000     0.00000     0.00000     0.00000     0.00000     0.00000

PIEZOELECTRIC TENSOR (including local field effects) (C/m^2)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x     0.00000     0.00000     0.00000     0.00000     0.00000     0.00000
  y     0.00000     0.00000     0.00000     0.00000     0.00000     0.00000
  z     0.00000     0.00000     0.00000     0.00000     0.00000     0.00000

请教下,上面这段中的玻恩有效电荷要怎么分析?我想算的是铝酸镧的各原子的BEC,求大神赐教
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

主做电弧增材制造方向
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

La:4.52078
Al:2.92932
O垂直:-2.43128
O平行:-2.44923
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
2楼2014-07-08 11:52:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)

引用回帖:
2楼: Originally posted by iamikaruk at 2014-07-08 11:52:25
La:4.52078
Al:2.92932
O垂直:-2.43128
O平行:-2.44923

有两个小疑问师兄,:第一个是,我这里用的POSCAR我不知道合不合理,用的是优化之后的CONTCAR,不过Direct坐标部分基本是整数;还是应该用没优化之前的POSCAR呢?
第二个是,这个BEC的结果为啥我 grep 不出来呀,上面的结果是我打开OUTCAR里查找出来的。。。。
主做电弧增材制造方向
3楼2014-07-08 14:37:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)

引用回帖:
2楼: Originally posted by iamikaruk at 2014-07-08 11:52:25
La:4.52078
Al:2.92932
O垂直:-2.43128
O平行:-2.44923

我在算钛酸铅(PTO)的BEC 时候,跟师兄师姐的对比了下都偏小,所以怀疑是不是POSCAR用的不合理。。。
主做电弧增材制造方向
4楼2014-07-08 14:38:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)

引用回帖:
2楼: Originally posted by iamikaruk at 2014-07-08 11:52:25
La:4.52078
Al:2.92932
O垂直:-2.43128
O平行:-2.44923

还有个疑问,师兄,这是我刚算的PTO(钛酸铅)的BEC ,刚用的LEPSILON算出来的,怎么每个原子出来的数值不一样了呢?取平均值吗?
BORN EFFECTIVE CHARGES (in e, cummulative output)
-------------------------------------------------
ion    1
    1     3.93395     0.00000     0.00000
    2     0.00000     3.93425     0.00000
    3     0.00000     0.00000     3.93407
ion    2
    1     7.21256     0.00000     0.00000
    2     0.00000     7.21235     0.00000
    3     0.00000     0.00000     7.21237
ion    3
    1    -2.59456     0.00000     0.00000
    2     0.00000    -5.94618     0.00000
    3     0.00000     0.00000    -2.59451
ion    4
    1    -5.94537     0.00000     0.00000
    2     0.00000    -2.59464     0.00000
    3     0.00000     0.00000    -2.59451
ion    5
    1    -2.59456     0.00000     0.00000
    2     0.00000    -2.59464     0.00000
    3     0.00000     0.00000    -5.94593
主做电弧增材制造方向
5楼2014-07-08 15:29:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
3楼: Originally posted by jk射手座 at 2014-07-08 14:37:09
有两个小疑问师兄,:第一个是,我这里用的POSCAR我不知道合不合理,用的是优化之后的CONTCAR,不过Direct坐标部分基本是整数;还是应该用没优化之前的POSCAR呢?
第二个是,这个BEC的结果为啥我 grep 不出来呀 ...

看你关心的是什么结构的bec,跟优化不优化关系倒不是太大。
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
6楼2014-07-08 23:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
5楼: Originally posted by jk射手座 at 2014-07-08 15:29:32
还有个疑问,师兄,这是我刚算的PTO(钛酸铅)的BEC ,刚用的LEPSILON算出来的,怎么每个原子出来的数值不一样了呢?取平均值吗?
BORN EFFECTIVE CHARGES (in e, cummulative output)
------------- ...

因为vasp在对称性的处理上不是很严格吧。
如果是abinit的话这些数值是严格符合space group的要求的。
在你的这种情况中,Pb和Ti取平均,O垂直和O平行分别取平均

» 本帖已获得的红花(最新10朵)

[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
7楼2014-07-08 23:15:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)

送红花一朵
引用回帖:
7楼: Originally posted by iamikaruk at 2014-07-08 23:15:18
因为vasp在对称性的处理上不是很严格吧。
如果是abinit的话这些数值是严格符合space group的要求的。
在你的这种情况中,Pb和Ti取平均,O垂直和O平行分别取平均...

呀,还真猜对了,多谢师兄!
主做电弧增材制造方向
8楼2014-07-09 08:24:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)

引用回帖:
7楼: Originally posted by iamikaruk at 2014-07-08 23:15:18
因为vasp在对称性的处理上不是很严格吧。
如果是abinit的话这些数值是严格符合space group的要求的。
在你的这种情况中,Pb和Ti取平均,O垂直和O平行分别取平均...

我又对比了一下,师兄,LEPSILON这个微扰理论算的结果是不怎么对称,然后用LCALCEPS算的结果就比较对称,而且二者结果很相近,算是误差可接受范围内吧:
LEPSILON的:(铌酸钾KNO的BEC)
BORN EFFECTIVE CHARGES (in e, cummulative output)
-------------------------------------------------
ion    1
    1     1.14082     0.00000     0.00000
    2     0.00000     1.14027     0.00000
    3     0.00000     0.00000     1.14096
ion    2
    1     9.73272     0.00000     0.00000
    2     0.00000     9.76588     0.00000
    3     0.00000     0.00000     9.72262
ion    3
    1    -1.74789     0.00000     0.00000
    2     0.00000    -7.34884     0.00000
    3     0.00000     0.00000    -1.74115
ion    4
    1    -7.36245     0.00000     0.00000
    2     0.00000    -1.77111     0.00000
    3     0.00000     0.00000    -1.74115
ion    5
    1    -1.74789     0.00000     0.00000
    2     0.00000    -1.77111     0.00000
    3     0.00000     0.00000    -7.36584

LCALCEPS的:
BORN EFFECTIVE CHARGES (including local field effects)
------------------------------------------------------
ion    1
    1     1.14300     0.00000     0.00000
    2     0.00000     1.14300     0.00000
    3     0.00000     0.00000     1.14300
ion    2
    1     9.27829     0.00000     0.00000
    2     0.00000     9.27829     0.00000
    3     0.00000     0.00000     9.27829
ion    3
    1    -1.64676     0.00000     0.00000
    2     0.00000    -6.98256     0.00000
    3     0.00000     0.00000    -1.64676
ion    4
    1    -6.98256     0.00000     0.00000
    2     0.00000    -1.64676     0.00000
    3     0.00000     0.00000    -1.64676
ion    5
    1    -1.64676     0.00000     0.00000
    2     0.00000    -1.64676     0.00000
    3     0.00000     0.00000    -6.98256

这个差别是不是因为BORN EFFECTIVE CHARGES 后边括号里的那个注释呢,比如内电场作用。
主做电弧增材制造方向
9楼2014-07-09 08:50:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyqh

铁杆木虫 (著名写手)

新人顶顶楼主的帖子,话说我平时算BORN EFFECTIVE CHARGES时候,都是设置参数LEPSILON来实现的?
对于使用LCALCEPS还没有使用过,不知道两者算出来的结果差别大不?
雨后方能见彩虹
10楼2014-07-09 14:44:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jk射手座 的主题更新
信息提示
请填处理意见