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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » berry phase 参考态
363/4返回列表上一页1234下一页
查看: 2030  |  回复: 35
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 的主题更新
363/4返回列表上一页1234下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭