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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问LCALCPOL和LBERRY计算的极化有什么区别
北京石油化工学院2026年研究生招生接收调剂公告
565/6返回列表上一页123456下一页
查看: 4959  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285求调剂 +7 AZMK 2026-04-04 9/450 2026-04-06 00:06 by 永字号
[考研] 材料调剂 +14 壹贰贰亿 2026-04-04 14/700 2026-04-05 23:31 by 来看流星雨10
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 22:23 by hemengdong
[考研] 一志愿南昌大学，085600，344分求调剂 +6 调剂上岸玘 2026-04-05 6/300 2026-04-05 22:11 by 789风
[考研] 求调剂 +6 chenxrlkx 2026-04-05 6/300 2026-04-05 22:05 by imissbao
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 312求调剂 +3 Say Never 2026-04-04 3/150 2026-04-05 11:19 by guoweigw
[考研] 413求调剂 +4 柯某某 2026-03-31 4/200 2026-04-04 22:18 by 学员6BFVa3
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 085601一志愿北理325分求调剂 +6 找调剂，， 2026-04-02 6/300 2026-04-03 22:20 by 柟�?
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
[考研] 320求调剂 +5 振—TZ 2026-04-02 5/250 2026-04-03 14:42 by fxue1114
[考研] 08工科，295，接受跨专业调剂 +8 lmnlzy 2026-03-30 8/400 2026-04-03 13:08 by nalakaiqi
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 312求调剂 +6 小小墨123 2026-04-02 7/350 2026-04-03 07:32 by jsw79
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 一志愿北交大材料工程总分358 +8 cs0106 2026-04-01 9/450 2026-04-02 10:36 by 不吃魚的貓
信息提示
关闭
请填处理意见
关闭