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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问LCALCPOL和LBERRY计算的极化有什么区别
565/6返回列表上一页123456下一页
查看: 4932  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +9 不知道叫什么！ 2026-03-15 10/500 2026-03-16 10:04 by houyaoxu
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 304求调剂 +4 素年祭语 2026-03-15 4/200 2026-03-16 09:42 by 闲人终南山
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 331求调剂（0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 一志愿天津大学，英一数二305分求调剂，四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 考研材料与化工，求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 0703，333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-12 5/250 2026-03-13 10:56 by houyaoxu
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
信息提示
关闭
请填处理意见
关闭