24小时热门版块排行榜    

查看: 323  |  回复: 2
当前主题已经存档。

nownok

捐助贵宾 (初入文坛)

[交流] 【求助】铁电极化计算

length of vectors
     4.087600000  4.087600000  4.086500000     
IGPAR=1 Expectation value term: ev
= (      .00000,     -.00265,     -.03018 ) electrons Angst
Berry-Phase term: bp
= (     1.38505,      .00000,      .00000 ) electrons Angst
ionic term: ion
= (    40.87600,    40.87600,    29.99758 ) electrons Angst

IGPAR=2 Expectation value term: ev
= (     -.00265,      .00000,     -.03018 ) electrons Angst
Berry-Phase term: bp
= (      .00000,      .77698,      .00000 ) electrons Angst
ionic term: ion
= (    40.87600,    40.87600,    29.99758 ) electrons Angst

IGAPR=3 Expectation value term: ev
= (      .00000,      .00000,     -.03018 ) electrons Angst
Berry-Phase term: bp
= (      .00000,      .00000,     2.08675 ) electrons Angst
ionic term: ion
= (    40.87600,    40.87600,    29.99758 ) electrons Angst
我明明算的P4mm结构,最后CONTCAR也是P4mm,为啥berry相给在a,b方向给出了非零值,而且还不等,显然也不是晶格常数的倍数

[ Last edited by nownok on 2009-7-5 at 10:43 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nownok

捐助贵宾 (初入文坛)

补充个INCAR

静态计算的
PREC = High
ISTART = 0 ; ICHARG = 2
EDIFF = 1E-8
EDIFFG = 1E-6
NSW = 0
ISIF = 2
ISMEAR = -5
IBRION = -1


拷贝了WAVECAR,CHGCAR过去,Berry相计算的
PREC = High
ISTART = 1 ; ICHARG = 11
EDIFF = 1E-8
EDIFFG = 1E-6
LBERRY = .TRUE.
IGPAR = 1
NPPSTR = 6
DIPOL = 0.5 0.5 0.5
ISMEAR = -5
IBRION = -1
2楼2009-07-05 10:47:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nownok

捐助贵宾 (初入文坛)

我实在是没看出来a,b方向如何不等价,根本有个4次轴嘛,凭啥,上次有人说Berry相计算有Bug,我这个可属于这样情况?谢谢高手解答
1.00000000000000000
     4.0876000000000001     .0000000000000000     .0000000000000000
      .0000000000000000    4.0876000000000001     .0000000000000000
      .0000000000000000     .0000000000000000    4.0865000000000000
   1   1   3
Selective dynamics
Direct
   .0000000000000000   .0000000000000000   .0000000000000000   F   F   F
   .5000000000000000   .5000000000000000   .5521400993164676   T   T   T
   .5000000000000000   .5000000000000000   .1147104139122916   T   T   T
   .0000000000000000   .5000000000000000   .6208069044510225   T   T   T
   .5000000000000000   .0000000000000000   .6208069044510225   T   T   T

[ Last edited by nownok on 2009-7-5 at 14:27 ]
3楼2009-07-05 14:25:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nownok 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见