24小时热门版块排行榜    

查看: 4838  |  回复: 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的回帖

clifflee

金虫 (小有名气)

2楼2016-03-28 15:48:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

clifflee

金虫 (小有名气)

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

多了个r
6楼2016-03-28 21:12:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

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

少了个整数倍
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
7楼2016-03-29 10:31:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clifflee

金虫 (小有名气)

引用回帖:
7楼: Originally posted by iamikaruk at 2016-03-29 10:31:23
少了个整数倍...

谢谢!
8楼2016-03-29 12:47:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qinshujiu

木虫 (小有名气)

9楼2016-03-30 13:11:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gooselin

木虫 (小有名气)

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

楼主你好,我也在学习铁电计算,这个极化量子2eR/V是什么意思啊?e和R代表什么,V是表示体积吗?非常感谢!
10楼2016-09-18 14:31:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 clifflee 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 28/1400 2026-08-13 17:51 by 医学老男孩
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看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
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
信息提示
请填处理意见