24小时热门版块排行榜    

查看: 5541  |  回复: 12

elliot.cc

铁虫 (初入文坛)

[交流] VASP自发极化强度怎么计算 已有5人参与

想请教一下大家,自发极化强度是怎么计算的?比如BaTiO3,居里温度两侧的铁电相和顺电相的空间群分别是p4mm和pm3m,用berry phase方法计算极化强度。是用BaTiO3的p4mm结构的极化减去pm3m结构的极化?还是BaTiO3 p4mm实际畸变结构的计划减去完美 p4mm结构的计划? 一直不得其解,请大家指教。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elliot.cc

铁虫 (初入文坛)

另外,想问下,计算极化强度的时候不是不能把原子移出晶胞之外嘛,那BaTiO3的氧原子不就移出晶胞外了吗,这该怎么办?
2楼2014-03-27 15:37:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

werhaohao

金虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
似乎就是减去一个无畸变的参考相的值,构造一个中心对称的P4mm结构计算极化。
Come on!
3楼2014-04-10 15:36:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zique

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
相对于中心对称的P4mm结构;移出晶胞的氧原子应该没问题吧,由于有周期性边界条件,晶胞中实际上是有对应的氧原子的
4楼2014-04-10 15:46:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elliot.cc

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by werhaohao at 2014-04-10 15:36:16
似乎就是减去一个无畸变的参考相的值,构造一个中心对称的P4mm结构计算极化。

我发现这个无畸变的p4mm比较好构建,但是实际畸变的p4mm移出晶胞外的氧原子VASP好像没有计算在内,导致c轴方向极化强度很大,一直解决不了
5楼2014-04-16 12:00:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elliot.cc

铁虫 (初入文坛)

引用回帖:
4楼: Originally posted by zique at 2014-04-10 15:46:57
相对于中心对称的P4mm结构;移出晶胞的氧原子应该没问题吧,由于有周期性边界条件,晶胞中实际上是有对应的氧原子的

但是原来上下两个氧原子在面上,各占0.5,现在一个出去了,一个进来了,分别占0和1,从计算的结果上来看也能反映出来,比如c轴方向我算出来的值是2.6左右,比实验值大了10倍,实在找不出其他原因了,苦恼
6楼2014-04-16 12:02:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

werhaohao

金虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by elliot.cc at 2014-04-16 12:00:02
我发现这个无畸变的p4mm比较好构建,但是实际畸变的p4mm移出晶胞外的氧原子VASP好像没有计算在内,导致c轴方向极化强度很大,一直解决不了...

我PBE优化出来的铁电相算出的极化强度大致是用实验上的结构参数算出来的2倍。
Come on!
7楼2014-04-16 22:30:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

werhaohao

金虫 (初入文坛)

引用回帖:
6楼: Originally posted by elliot.cc at 2014-04-16 12:02:54
但是原来上下两个氧原子在面上,各占0.5,现在一个出去了,一个进来了,分别占0和1,从计算的结果上来看也能反映出来,比如c轴方向我算出来的值是2.6左右,比实验值大了10倍,实在找不出其他原因了,苦恼...

能给看一下你计算的参考相和铁电相里OUTCAR里极化强度么
Come on!
8楼2014-04-16 22:32:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by elliot.cc at 2014-04-16 12:00:02
我发现这个无畸变的p4mm比较好构建,但是实际畸变的p4mm移出晶胞外的氧原子VASP好像没有计算在内,导致c轴方向极化强度很大,一直解决不了...

还没有解决这个问题啊?
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
9楼2014-04-17 02:17:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elliot.cc

铁虫 (初入文坛)

引用回帖:
9楼: Originally posted by iamikaruk at 2014-04-17 02:17:26
还没有解决这个问题啊?...

解决了!多谢你的帮助哈!
10楼2014-05-07 20:28:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 elliot.cc 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天放榜吗? +12 布布和一二 2026-08-19 13/650 2026-08-19 14:33 by 陈晨晨陈啊
[基金申请] filecode,4个jtjc了 +8 ziyangfang 2026-08-19 9/450 2026-08-19 13:08 by blacksandy
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +52 医学老男孩 2026-08-13 115/5750 2026-08-19 12:09 by 医学老男孩
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +4 Tide man 2026-08-12 5/250 2026-08-19 07:41 by 201120518022
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
信息提示
请填处理意见