应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问LCALCPOL和LBERRY计算的极化有什么区别
202/2返回列表上一页12
查看: 4778  |  回复: 55
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jk射手座

新虫 (著名写手)

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

英文手册上的LCALCPOL和LBERRY 部分都看过，比较混乱，在豆丁上找到了相应部分的翻译版，还是有几个问题：
1、只是多知道一点LCALCPOL是算宏观电子极化，那这个是什么意思？LBERRY 算的是微观电子极化吗？？宏观、微观又是什么区别呢？
2、LCAL 和LBERRY 算出来的极化理论上值应该是一样或是相近的吗？
3、二者分别应该是在什么情况下优先选用呢？
4、LCAL的计算说明里提到手动移动原子（这里我一般两个POSCAR1，POSCAR2，分数坐标我把自洽之后的坐标当做是手动移动过原子了，这样做可以吗），这点上LBERRY 也是一样的吗？还是直接用模型的规矩的分数坐标POSCAR直接算？

另外有个公式想求证下各位前辈，极化P={（P2-P1）*16.022}/V， 这个公式算极化是正确的吗？单位由c/m2变为 uc/cm2 是直接乘以100就可以了吧？

问得比较多，最近的问题导致开始怀疑计算这一系列的方法的正确性，所以希望看到的路过的版主、大神们，当然还有偶尔出没的侯老师帮解下疑！小弟万分感谢啊！
回复此楼

» 猜你喜欢
主做电弧增材制造方向
1楼 2014-06-24 15:19:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

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

什么叫做移动同样的距离？你只要让它保持有中心反演对称就可以了。
不过事实上还要考虑算出来能不能保证是绝缘的……
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
37楼2014-08-26 22:25:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
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)的移动；还有说到要保证绝缘，最近确实在算多层结构的时候出现了错误提示说系统不绝缘，我还没想明白是什么原因 ...

移动不用一致，铁电相可以是弛豫的结果
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
39楼2014-08-28 17:08:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

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

弛豫时候用的是-5还是0？
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
41楼2014-08-29 14:52:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
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 结构弛豫都会出这个错 ...

这个就不知道了。。。没有做过
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
43楼2014-08-30 09:40:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

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

算macroscopic polarization需要用两个POSCAR，一个是参考的顺电相，一个是铁电相。
如果你的模拟体系不是slab的话，不用设置IDIPOL
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
46楼2014-10-04 22:11:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

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



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

LCALCPOL算出来的应该是个矢量，不是矢量和。
你最好先去看看现代极化理论的review文章，免得盲目计算。
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
48楼2014-10-05 09:25:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
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.).
  这段话 ...

当然不是矢量和了，只是沿三个方向分别计算而已。出来的是一个矢量不就知道c轴的极化了
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
50楼2014-10-05 22:43:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
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  ...

是的
回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
52楼2014-10-06 08:04:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
53楼: Originally posted by jk射手座 at 2014-10-06 09:20:42

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

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

macroscopic polarization=铁电相的formal polarization-顺电相的formal polarization
赞一下(1人) 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
54楼2014-10-06 10:32:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jk射手座 的主题更新
202/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭