24小时热门版块排行榜    

查看: 5111  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见