24小时热门版块排行榜    

查看: 2139  |  回复: 35

tuhaidan

银虫 (小有名气)

引用回帖:
19楼: Originally posted by iamikaruk at 2012-09-13 20:39:36
1.0000和0.0000是一样的....

优化前:
    POSITION         DIRECTION 3           BORN EFFECTIVE CHARGE    (rigid.aug.   ionic)
      1.99600      1.99600      0.00000     0.00000  0.00000 -6.10018 ( -0.62733  6.00000)
      0.00000      1.99600      2.01750     0.00000  0.00000 -2.14138 ( -0.62823  6.00000)
      1.99600      0.00000      2.01750     0.00000  0.00000 -2.14138 ( -0.62823  6.00000)
      1.99600      1.99600      2.01750     0.00000  0.00000  7.66991 ( -0.93586  4.00000)
      0.00000      0.00000      0.00000     0.00000  0.00000  2.71799 (  1.70517 10.00000)
ion    1
    1    -2.17589     0.00000     0.00000
    2     0.00000    -2.16014     0.00000
    3     0.00000     0.00000    -6.10018
ion    2
    1    -6.27563     0.00000     0.00000
    2     0.00000    -2.14016     0.00000
    3     0.00000     0.00000    -2.14138
ion    3
    1    -2.13176     0.00000     0.00000
    2     0.00000    -6.27949     0.00000
    3     0.00000     0.00000    -2.14138
ion    4
    1     7.86272     0.00000     0.00000
    2     0.00000     7.85854     0.00000
    3     0.00000     0.00000     7.66991
ion    5
    1     2.72494     0.00000     0.00000
    2     0.00000     2.72557     0.00000
    3     0.00000     0.00000     2.71799
优化后:
POSITION         DIRECTION 3           BORN EFFECTIVE CHARGE    (rigid.aug.   ionic)
-----------------------------------------------------------------------------------------
      1.99600      1.99600      3.93009     0.00000  0.00000 -4.86448 ( -0.62850  6.00000)
      1.99600      0.00000      1.96908     0.00000  0.00000 -1.94027 ( -0.62888  6.00000)
      0.00000      1.99600      1.96908     0.00000  0.00000 -1.94027 ( -0.62888  6.00000)
      1.99600      1.99600      2.09820     0.00000  0.00000  5.92345 ( -0.95221  4.00000)
      0.00000      0.00000      0.00000     0.00000  0.00000  2.82578 (  1.70531 10.00000)
PIEZOELECTRIC TENSOR  for field in x, y, z        (e  Angst)
          XX          YY          ZZ          XY          YZ          ZX
  --------------------------------------------------------------------------------
  x    -0.00117     0.00075    -0.00012     0.00000     0.00000     0.53845
  y     0.00117    -0.00179     0.00024     0.00000     0.53878     0.00000
  z     1.66535     1.66535     0.25263     0.00000     0.00000     0.00000

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.13417
  y     0.00000     0.00000     0.00000     0.00000     0.13421     0.00000
  z     0.41495     0.41495     0.06295     0.00000     0.00000     0.00000
  ion    1
    1    -1.98675     0.00000     0.00111
    2     0.00000    -1.98679     0.00110
    3     0.00000     0.00000    -4.86448
ion    2
    1    -2.15716     0.00000    -0.00020
    2     0.00000    -5.91181     0.00096
    3     0.00000     0.00000    -1.94027
ion    3
    1    -5.91189     0.00000     0.00114
    2     0.00000    -2.15722    -0.00002
    3     0.00000     0.00000    -1.94027
ion    4
    1     7.35215     0.00000    -0.00150
    2     0.00000     7.35211    -0.00149
    3     0.00000     0.00000     5.92345
ion    5
    1     2.70938     0.00000    -0.00004
    2     0.00000     2.70921    -0.00006
    3     0.00000     0.00000     2.82578
应该如何处理这些数据啊?原来前面我真的算错了
21楼2012-09-14 11:10:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
21楼: Originally posted by tuhaidan at 2012-09-14 11:10:06
优化前:
    POSITION         DIRECTION 3           BORN EFFECTIVE CHARGE    (rigid.aug.   ionic)
      1.99600      1.99600      0.00000     0.00000  0.00000 -6.10018 ( -0.62733  6.00000)
       ...

什么晶系?把坐标做平移之后相减就可以了。
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
22楼2012-09-14 17:11:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuhaidan

银虫 (小有名气)

引用回帖:
22楼: Originally posted by iamikaruk at 2012-09-14 17:11:22
什么晶系?把坐标做平移之后相减就可以了。...

POSCAR的是direct的  上面的是cartesian的,这个我会,然后的所谓的born effctive charge张量指的是这个吗?
ion    1
    1    -1.98675     0.00000     0.00111
    2     0.00000    -1.98679     0.00110
    3     0.00000     0.00000    -4.86448
它乘以原子位移,如:
ion    1
    1    -1.98675     0.00000     0.00111
    2     0.00000    -1.98679     0.00110
    3     0.00000     0.00000    -4.86448    *【(1.99600      1.99600      3.93009)- (1.99600      1.99600      0.00000)】=P是这么理解的吗?
23楼2012-09-14 17:57:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuhaidan

银虫 (小有名气)

引用回帖:
22楼: Originally posted by iamikaruk at 2012-09-14 17:11:22
什么晶系?把坐标做平移之后相减就可以了。...

其实  这个可以用berry phase直接算电极化的,然后我也想用born effective charge的方法估算一下
24楼2012-09-14 18:02:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
23楼: Originally posted by tuhaidan at 2012-09-14 17:57:02
POSCAR的是direct的  上面的是cartesian的,这个我会,然后的所谓的born effctive charge张量指的是这个吗?
ion    1
    1    -1.98675     0.00000     0.00111
    2     0.00000    -1.98679     0.00110 ...

3.93009那个不对,你平移了么?
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
25楼2012-09-14 21:18:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuhaidan

银虫 (小有名气)

引用回帖:
25楼: Originally posted by iamikaruk at 2012-09-14 21:18:13
3.93009那个不对,你平移了么?...

这个就是优化前为0.00    优化后为0.9829307947555620的,我也觉这个地方有问题,但是不明白你指的平移是什么,能不能一次说明白啊
26楼2012-09-14 22:14:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
26楼: Originally posted by tuhaidan at 2012-09-14 22:14:07
这个就是优化前为0.00    优化后为0.9829307947555620的,我也觉这个地方有问题,但是不明白你指的平移是什么,能不能一次说明白啊...

……晶体平移对称性啊
1-0.9829307xxxxx就是位移了
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
27楼2012-09-15 00:02:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuhaidan

银虫 (小有名气)

引用回帖:
27楼: Originally posted by iamikaruk at 2012-09-15 00:02:56
……晶体平移对称性啊
1-0.9829307xxxxx就是位移了...

计算Born effective charge张量时,是不是采用铁电相的POSCAR?  POSCAR里是Direct坐标,求得的Born effective charge 张量是什么坐标?
位移量应该用什么坐标?
我记得应该用cartesian吧,我怎么觉得我现在越来越糊涂了呢
你能说一下完整的过程吗?
我说一下我的做法吧:
用铁电相的POSCAR求的Born effective charge张量,同时在OUTCAR中得到它的cartesian原子坐标。
然后求得perfect POSCAR的cartesian原子坐标同样是OUTCAR中得到。
前后两个cartesian原子坐标相减,得到位移。
最后就是Born effective charge张量*位移=P
请指出错误在哪里,好吗
28楼2012-09-15 09:45:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)


xueht987: 金币+1, 鼓励交流! 2012-09-16 10:33:20
引用回帖:
28楼: Originally posted by tuhaidan at 2012-09-15 09:45:40
计算Born effective charge张量时,是不是采用铁电相的POSCAR?  POSCAR里是Direct坐标,求得的Born effective charge 张量是什么坐标?
位移量应该用什么坐标?
我记得应该用cartesian吧,我怎么觉得我现在越来 ...

用顺电相的POSCAR求Born effective charge张量
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
29楼2012-09-15 11:01:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tuhaidan

银虫 (小有名气)

引用回帖:
29楼: Originally posted by iamikaruk at 2012-09-15 11:01:03
用顺电相的POSCAR求Born effective charge张量...

啊?那其他的有没有问题?  那个位移的算法呢,错没有错啊?
顺电相应该是1.00还是0.00?
30楼2012-09-15 12:51:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tuhaidan 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 寻找调剂 +3 倔强芒? 2026-03-21 3/150 2026-03-22 14:27 by 30660438
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 306求调剂 +4 来好运来来来 2026-03-22 4/200 2026-03-22 12:43 by 素颜倾城1988
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 299求调剂 +4 某某某某位 2026-03-21 4/200 2026-03-21 16:30 by barlinike
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
请填处理意见