应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问OUTCAR中的这段玻恩有效电荷的信息怎么分析呀
151/2返回列表12下一页
查看: 2280  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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，求大神赐教
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
主做电弧增材制造方向
1楼 2014-07-08 09:40:41
已阅   回复此楼   关注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朵）

jk射手座
[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射手座 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 5/250 2026-03-19 11:08 by kakakapanpan
[考研] 332求调剂 +3 ydfyh 2026-03-17 3/150 2026-03-19 10:14 by 功夫疯狂
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 334求调剂 +3 志存高远意在机� 2026-03-16 3/150 2026-03-18 08:34 by lm4875102
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 0856专硕279求调剂 +5 加油加油！? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 复试调剂 +3 呼呼？~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 328化工专硕求调剂 +4 。，。，。，。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
信息提示
关闭
请填处理意见
关闭