24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2170  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 电气专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 02:43 by 秋豆菜芽
[考研] 求调剂,262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen~
[考研] 339求调剂 +7 hanwudada 2026-04-11 7/350 2026-04-11 23:06 by 852137818
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 22408 352分求调剂 +4 努力的夏末 2026-04-09 4/200 2026-04-11 10:42 by maddjdld
[考研] 还有化工二轮调剂的学校吗 5+14 化工人999 2026-04-09 48/2400 2026-04-11 10:27 by 89436494
[考研] 本科211 工科085400 280分求调剂 可跨专业 +11 LZH(等待调剂中 2026-04-10 11/550 2026-04-11 08:39 by zhq0425
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 一志愿2110,化学学硕310分,本科重点双非求调剂 +18 努力奋斗112 2026-04-08 18/900 2026-04-09 23:28 by wolf97
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 一志愿中国科学院上海有机所,有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 322求调剂,08工科 +3 今天是个小号 2026-04-08 3/150 2026-04-09 15:53 by wp06
[考研] 材料专硕322 +14 哈哈哈吼吼吼哈 2026-04-05 14/700 2026-04-09 13:25 by 5268321
[考研] 化工学硕 285求调剂 +26 Wisjxn 2026-04-07 26/1300 2026-04-08 14:42 by screening
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
信息提示
请填处理意见