24小时热门版块排行榜    

查看: 461  |  回复: 1

yayanoba

新虫 (小有名气)

[求助] 为什么我得不出手册上计算NaF极化的值呢?Berry phase

我按着手册跑了一遍NaF的Berry phase计算,得到的结果和手册上基本相同了,但是我为什么就是得不出最后的结果E=0.04489.
我经过计算以后的各项值分别 如下:
r=(-0.0000033, 0, 0.00003)
r=(0.00001, 0.00001, 0.36054)
p[ion]=(0, 0, -0.94713)
以上的值是已经减过参考态的数值。我怎么加也加不出0.04489的值啊?请大侠赐教。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yayanoba

新虫 (小有名气)

补一下原始数据:
参考态:
                              e_ev=(     0.00000     0.00000     0.00000 ) e*Angst
                                e_bp=(     0.00000     0.00000     0.00000 ) e*Angst

Total electronic dipole moment: p[elc]=(     0.00000     0.00000     0.00000 ) e*Angst

            ionic dipole moment: p[ion]=(     2.25510     2.25510     2.25510 ) e*Angst
三个方向G1 G2 G3 都是以上结果。
极化:
G1:
                                e_ev=(     0.00000     0.00000     0.00003 ) e*Angst
                                e_bp=(     0.00000     0.18027     0.18027 ) e*Angst

Total electronic dipole moment: p[elc]=(     0.00000     0.18026     0.18030 ) e*Angst

            ionic dipole moment: p[ion]=(     2.25510     2.25510     1.93939 ) e*Angst
G2:
                                e_ev=(     0.00000     0.00000     0.00003 ) e*Angst
                                e_bp=(     0.18027     0.00000     0.18027 ) e*Angst

Total electronic dipole moment: p[elc]=(     0.18026     0.00000     0.18030 ) e*Angst

            ionic dipole moment: p[ion]=(     2.25510     2.25510     1.93939 ) e*Angst
G3:                                e_ev=(    -0.00001     0.00000     0.00003 ) e*Angst
                                e_bp=(    -0.18026    -0.18026     0.00000 ) e*Angst

Total electronic dipole moment: p[elc]=(    -0.18027    -0.18025     0.00003 ) e*Angst

            ionic dipole moment: p[ion]=(     2.25510     2.25510     1.93939 ) e*Angst
2楼2013-01-16 16:44:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yayanoba 的主题更新
信息提示
请填处理意见