[求助] 请问LCALCPOL和LBERRY计算的极化有什么区别

36楼: Originally posted by jk射手座 at 2014-08-26 22:20:07
谢谢师兄！我明白我错在哪儿了，最后一个问题呀，参考相所有的原子移动相同的距离吗？...

38楼: Originally posted by jk射手座 at 2014-08-27 08:59:43
不好意思师兄，说的太不专业了，就是原子的坐标移动相同的量，比如Ti(0.5 0.5 0.5)到Ti(0.5 0.5 0.51)的移动；还有说到要保证绝缘，最近确实在算多层结构的时候出现了错误提示说系统不绝缘，我还没想明白是什么原因 ...

40楼: Originally posted by jk射手座 at 2014-08-29 08:43:37
嗯师兄，我昨天照你说的做了，而且铁电相用的弛豫结果，不过有一部分提示铁电相不绝缘，然后参考相提示：系统不绝缘，我换了ISMEAR=-5也不管用，这个是不是原子坐标的问题呢？...

42楼: Originally posted by jk射手座 at 2014-08-30 08:52:19
只有1:1的两层超晶格结构是-5,2：2以上的到5:5我都用的0，不然总会出现BRMIX: very serious problems
the old and the new charge density differ的错误提示，目前连3:3以上的PbTiO3/KNbO3 结构弛豫都会出这个错 ...

45楼: Originally posted by jk射手座 at 2014-10-04 15:40:26
刚马上又想到一个可能我出错的地方。。。
师兄，是不是算一个体系的极化强度，直接一个POSCAR就可以了？？比如用自洽后的POSCAR，然后INCAR里设置LCALCPOL让它直接跑出来？（还有是否需要设置IDIPOL，因 ...

47楼: Originally posted by jk射手座 at 2014-10-05 08:19:42
哦这样啊，可是如果不设置IDIPOL，LCALCPOL算出来的好像是xyz三个方向极化的矢量和，而我只需要得到c方向的就可以呀，这样也不用设置IDIPOL=3吗？！



另外一个问题，就是POSCAR，还是不太确定，参考的顺电相 ...

49楼: Originally posted by jk射手座 at 2014-10-05 10:31:23
For LCALCPOL=.TRUE., however, VASP calculates the electronic contribution to the polarization, along the three reciprocal lattice vectors  (i.e.  ) in a single run (unlike LBERRY=.TRUE.).
  这段话 ...

51楼: Originally posted by jk射手座 at 2014-10-05 23:23:13

意思是……我提取‘moment’ 得到
Ionic dipole moment: p=(   -44.39000   -44.39000    -8.35428 ) electrons Angst
Total electronic dipole moment: p=(     0.00000     0.00000     4.78781  ...

53楼: Originally posted by jk射手座 at 2014-10-06 09:20:42

太谢谢你了师兄，我也知道自己问的问题挺低级的。。。谢谢你这么不厌其烦地指教，嘿嘿！

最后确认一件事情，参考相的那部分算出来的极化还有什么用吗？
OUTCAR2里的是铁电相， OUTCAR1里是参考相
#  ...