24小时热门版块排行榜    

查看: 3825  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

leonidasn

新虫 (初入文坛)

[求助] 求问,VASP计算born有效电荷后估算极化强度的方法

新手入门,刚接触VASP,很多地方都不懂
最近用vasp算BFO体系,照着文献算极化强度。算出来移动Bi原子后的有效电荷值是4.83,跟文献值的4.38差距不大,但是始终无法按照文献里的公式折算出极化强度,还望各位高人指点一二下面是我移动前和移动后的POSCAR文件
移动之前
  O   Fe   Bi
1.0
        5.6970000267         0.0000000000         0.0000000000
        2.9141183529         4.8952756337         0.0000000000
        2.9141183529         1.6566271578         4.6064422269
    6    2    2
Cartesian
     6.896384774         5.220795805         1.780389925
     6.485353118         4.165162336         4.310247989
     8.013711728         2.776994286         2.460761496
     5.165211715         1.636018466         4.083611039
     4.047884760         5.061454305         2.007026875
     6.419734778         2.189528260         0.157540383
     2.572432762         1.462384660         1.028157875
     8.335050785         4.738335860         3.331378851
     0.000000000         0.000000000         0.000000000
     5.762618366         3.275951396         2.303221113

移动之后
O   Fe   Bi
1.0
        5.6970000267         0.0000000000         0.0000000000
        2.9141183529         4.8952756337         0.0000000000
        2.9141183529         1.6566271578         4.6064422269
    6    2    2
Cartesian
     6.896384774         5.220795805         1.780389925
     6.485353118         4.165162336         4.310247989
     8.013711728         2.776994286         2.460761496
     5.165211715         1.636018466         4.083611039
     4.047884760         5.061454305         2.007026875
     6.419734778         2.189528260         0.157540383
     2.572432762         1.462384660         1.028157875
     8.335050785         4.738335860         3.331378851
     0.020000000         0.020000000         0.020000000
     5.762618366         3.275951396         2.303221113
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

vasp4.6支持DFPT计算BEC吗?
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
5楼2012-10-16 16:35:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
6楼: Originally posted by leonidasn at 2012-10-17 15:13:21
不支持。。。不过最近换5.2用LEPSILON算出了BEC,但是怎么从这个BEC得到极化数据呢?而且如何判断正负极化呢?...

你可以用BEC张量矩阵乘以原子的相对位移估算得到极化矢量,由矢量的方向判断正负极化
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
7楼2012-10-17 15:14:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
8楼: Originally posted by leonidasn at 2012-10-17 20:05:59
额,BEC乘以位移再除以体积可以估算极化这个我知道,但是用LEPSILON算出来的BEC好像没有给出对应的原子位移啊...

BEC本来就不会给出原子位移啊……它是极化对位移的偏微分
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
9楼2012-10-17 21:58:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
10楼: Originally posted by leonidasn at 2012-10-18 15:49:10
额,不好意思,那请问一下你说的用BEC张量矩阵乘以原子的相对位移估算得到极化矢量具体应该怎么操作呢?实在不是太懂啊。我看文献里面给出的公式一般都是我说那个意思啊...

[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
12楼2012-10-19 08:48:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)

引用回帖:
13楼: Originally posted by leonidasn at 2012-10-23 14:56:06
问题就是这个Uij如何得到啊?...

你的铁电相相对于顺电相的位移
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
14楼2012-10-23 20:44:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leonidasn 的主题更新
信息提示
请填处理意见