24小时热门版块排行榜    

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

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

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

铁虫 (初入文坛)

引用回帖:
8楼: Originally posted by werhaohao at 2014-04-16 22:32:11
能给看一下你计算的参考相和铁电相里OUTCAR里极化强度么...

这是参考相的
Total electronic dipole moment: p[elc]=(     0.00000     0.00000     0.00000 ) e*Angst

            ionic dipole moment: p[ion]=(    31.95600    31.95600    32.26800 ) e*Angst
这是铁电相的
Total electronic dipole moment: p[elc]=(     0.00000     0.00000     2.27151 ) e*Angst

            ionic dipole moment: p[ion]=(    31.95600    31.95600    31.17895 ) e*Angst
11楼2014-05-07 20:29:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 elliot.cc 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见