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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问OUTCAR中的这段玻恩有效电荷的信息怎么分析呀
151/2返回列表12下一页
查看: 2686  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了，能看出什么问题？ +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +50 医学老男孩 2026-08-13 109/5450 2026-08-18 16:56 by 医学老男孩
[基金申请] 重要消息，中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[教师之家] 两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 7/350 2026-08-18 16:42 by LNP@mRNA
[基金申请] 明天放榜？ +4 Shxjjxjkx 2026-08-18 4/200 2026-08-18 15:59 by abcabc1133
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
关闭
请填处理意见
关闭