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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问LCALCPOL和LBERRY计算的极化有什么区别
565/6返回列表上一页123456下一页
查看: 4692  |  回复: 55
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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的回帖

jk射手座

新虫 (著名写手)
引用回帖:
41楼: Originally posted by iamikaruk at 2014-08-29 14:52:49
弛豫时候用的是-5还是0？...

只有1:1的两层超晶格结构是-5,2：2以上的到5:5我都用的0，不然总会出现BRMIX: very serious problems
the old and the new charge density differ的错误提示，目前连3:3以上的PbTiO3/KNbO3 结构弛豫都会出这个错，原因还在找；而PbTiO3/LaAlO3的就不会。
赞一下 回复此楼
主做电弧增材制造方向
42楼2014-08-30 08:52:19
 已阅   回复此楼   关注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的回帖

jk射手座

新虫 (著名写手)
引用回帖:
11楼: Originally posted by iamikaruk at 2014-07-05 22:40:09
看了下tag，似乎lcalcpol是自动把lberry的所有工作给做了。

师兄，不好意思，这个就问题还有个地方需要请教您。
因为算的极化强度结果和文献一直有出入，我反复看了几遍LCALCPOL的tag之后想到一个问题，那个tag里的example是只移动的F原子，那么我在算整个体系的极化强度时要怎么移动原子，移动哪些原子呢？我目前用过的移动方式是，两个POSCAR一个用未自洽的整坐标，一个是自洽之后的坐标；另一种是POSCAR1 用自洽过后的坐标，POSCAR2是所有原子的z坐标移动0.01个单位（纯属模仿tag里的例子的方式），但是Pele 和Pion的结果处理完后得到的极化强度P太小了，所以我不知道是不是我移动原子的方式太盲目了，改怎么选定原子移动前后的POSCAR呢？
赞一下 回复此楼
主做电弧增材制造方向
44楼2014-10-04 15:34:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)
引用回帖:
43楼: Originally posted by iamikaruk at 2014-08-30 09:40:12
这个就不知道了。。。没有做过...

刚马上又想到一个可能我出错的地方。。。
师兄，是不是算一个体系的极化强度，直接一个POSCAR就可以了？？比如用自洽后的POSCAR，然后INCAR里设置LCALCPOL让它直接跑出来？（还有是否需要设置IDIPOL，因为我的内容里只关注c方向的极化）
赞一下 回复此楼
主做电弧增材制造方向
45楼2014-10-04 15:40:26
 已阅   回复此楼   关注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的回帖

jk射手座

新虫 (著名写手)
引用回帖:
46楼: Originally posted by iamikaruk at 2014-10-04 22:11:56
算macroscopic polarization需要用两个POSCAR，一个是参考的顺电相，一个是铁电相。
如果你的模拟体系不是slab的话，不用设置IDIPOL...

哦这样啊，可是如果不设置IDIPOL，LCALCPOL算出来的好像是xyz三个方向极化的矢量和，而我只需要得到c方向的就可以呀，这样也不用设置IDIPOL=3吗？！



另外一个问题，就是POSCAR，还是不太确定，参考的顺电相和铁电想POSCAR坐标要选哪个？不知道我算自洽后的坐标合不合理。。。因为参照英文说明书LCALCPOL-tag里的那个例子算的，结果不是很理想，我觉得我算的盲目了
赞一下 回复此楼
主做电弧增材制造方向
47楼2014-10-05 08:19:42
 已阅   回复此楼   关注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的回帖

jk射手座

新虫 (著名写手)
引用回帖:
48楼: Originally posted by iamikaruk at 2014-10-05 09:25:10
LCALCPOL算出来的应该是个矢量，不是矢量和。
你最好先去看看现代极化理论的review文章，免得盲目计算。...

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轴上的极化呢
赞一下 回复此楼
主做电弧增材制造方向
49楼2014-10-05 10:31:23
 已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 jk射手座 的主题更新
565/6返回列表上一页123456下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭