24小时热门版块排行榜    

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

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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +12 医学老男孩 2026-08-13 31/1550 2026-08-14 07:47 by 医学老男孩
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 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
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +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 布布和一二
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
信息提示
请填处理意见