24小时热门版块排行榜    

查看: 4877  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

clifflee

金虫 (小有名气)

[求助] VASP5.3.2用Berry phase 计算铁电自发极化结果怎么处理?

INCAR文件:
ENCUT=500
GGA=PS
ISTART=0
ICHARG=2
ISMEAR=0
SIGMA=0.05
ALGO=Normal
LREAL=Auto

EDIFF=1E-6
PREC=Accurate

DIPOL=0 0 0
LCALCPOL=.TRUE.
计算结果:
0%
Ionic dipole moment:                 p[ion]=( 0.00000 0.00000 -0.00002 ) electrons Angst
Total electronic dipole moment: p[elc]=( 0.00000 -0.00000 -0.00004 ) electrons Angst
10%
Ionic dipole moment:                 p[ion]=( 119.28048 -0.00000 -0.00001 ) electrons Angst
Total electronic dipole moment: p[elc]=( 0.92280 0.00000 -0.00016 ) electrons Angst
20%
Ionic dipole moment:                 p[ion]=( 107.91264 -0.00000 -0.00001 ) electrons Angst
Total electronic dipole moment: p[elc]=( 1.85015 0.00000 -0.00016 ) electrons Angst
30%
Ionic dipole moment:                 p[ion]=( 96.54480 -0.00000 -0.00001 ) electrons Angst
Total electronic dipole moment: p[elc]=( 2.78106 0.00000 -0.00016 ) electrons Angst
40%
Ionic dipole moment:                 p[ion]=( 85.17696 -0.00000 -0.00001 ) electrons Angst
Total electronic dipole moment: p[elc]=( 3.71535 0.00000 -0.00016 ) electrons Angst
50%
Ionic dipole moment:                 p[ion]=( 73.80912 -0.00000 -0.00000 ) electrons Angst
Total electronic dipole moment: p[elc]=( 4.65097 0.00000 -0.00016 ) electrons Angst
60%
Ionic dipole moment:                p[ion]=( 62.44129 -0.00000 -0.00002 ) electrons Angst
Total electronic dipole moment: p[elc]=( -5.29986 0.00000 -0.00016 ) electrons Angst
70%
Ionic dipole moment:                p[ion]=( 51.07344 -0.00000 -0.00001 ) electrons Angst
Total electronic dipole moment: p[elc]=( -4.36483 0.00000 -0.00017 ) electrons Angst
80%
Ionic dipole moment:                 p[ion]=( 39.70561 -0.00000 0.00001 ) electrons Angst
Total electronic dipole moment: p[elc]=( -3.43232 0.00000 -0.00017 ) electrons Angst
90%
Ionic dipole moment:                 p[ion]=( 28.33777 -0.00000 0.00000 ) electrons Angst
Total electronic dipole moment: p[elc]=( -2.50352 0.00000 -0.00017 ) electrons Angst
100%
Ionic dipole moment:                 p[ion]=( 16.96993 -0.00000 0.00000 ) electrons Angst
Total electronic dipole moment: p[elc]=( -1.58067 0.00000 -0.00017 ) electrons Angst
100%表示铁电相,0%表示对应的中心对称相。
100方向p[ion]+ p[elc]:
0%               0
10%             120.20328
20%             109.76279
30%             99.32586
40%             88.89231
50%             78.46009
60%             57.14143
70%             46.70861
80%             36.27329
90%             25.83425
100%           15.38926
如果这样处理得出下图,
VASP5.3.2用Berry phase 计算铁电自发极化结果怎么处理?
比较诡异,不知道是否处理正确。望各位指教。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lightgjx

木虫 (正式写手)

引用回帖:
5楼: Originally posted by clifflee at 2016-03-28 21:11:48
对于没有磁性的结构,极化量子是否为2erR/V?...

请问这个2erR/V?是从哪里得到的?谢谢
12楼2016-10-03 22:10:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

iamikaruk

木虫 (著名写手)

★ ★ ★ ★ ★
fzx2008: 金币+5, 谢谢指导,代发金币 2016-03-31 19:37:23
50%-60%的地方跳跃好厉害啊。
把除了0%外的所有极化值减去极化量子即可。
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
3楼2016-03-28 21:01:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
3楼: Originally posted by iamikaruk at 2016-03-28 21:01:05
50%-60%的地方跳跃好厉害啊。
把除了0%外的所有极化值减去极化量子即可。

跳跃有可能是要减去不同的极化量子
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
4楼2016-03-28 21:02:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clifflee

金虫 (小有名气)

引用回帖:
4楼: Originally posted by iamikaruk at 2016-03-28 21:02:56
跳跃有可能是要减去不同的极化量子...

对于没有磁性的结构,极化量子是否为2erR/V?
5楼2016-03-28 21:11:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 能否退出参与的面上项目解除限项 +22 koalala 2026-08-24 25/1250 2026-08-26 11:14 by dxcharlary
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 8/400 2026-08-26 11:03 by zjlp
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 出来了 +8 trojank 2026-08-26 8/400 2026-08-26 09:00 by gxc
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +24 木水思豆 2026-08-25 27/1350 2026-08-26 08:26 by yanchen918
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见