24小时热门版块排行榜    

查看: 4459  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 317求调剂 +8 申子申申 2026-03-19 13/650 2026-03-21 00:09 by 刘国森
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +4 一瓶苯甲酸 2026-03-14 4/200 2026-03-20 20:36 by fen_rao
[考研] 0703化学调剂 ,六级已过,有科研经历 +13 曦熙兮 2026-03-15 13/650 2026-03-20 19:35 by Dream007008
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
信息提示
请填处理意见