24小时热门版块排行榜    

查看: 539  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳又变了8-15 +3 archvillain 2026-08-15 3/150 2026-08-15 11:22 by 医学老男孩
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 09:28 by sunzitan
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 不应该看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
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见