应《网络安全法》要求，自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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京化工085600 310分求调剂 +16 0856材料与化工3 2026-04-04 18/900 2026-04-06 01:54 by BruceLiu320
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 22:23 by hemengdong
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 070300化学求调剂 +17 小黄鸭宝 2026-03-30 17/850 2026-04-05 12:03 by 宁馨哈哈
[考研] 271分求调剂学校 +12 zph158488！ 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 材料科学与工程调剂 +19 深V宿舍吧 2026-03-30 20/1000 2026-04-04 22:13 by hemengdong
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-04 20:37 by 蓝云思雨
[考研] 301求调剂 +18 骆驼男人 2026-04-02 18/900 2026-04-04 20:33 by 蓝云思雨
[考研] 一志愿北交大材料工程总分358 +6 cs0106 2026-04-03 6/300 2026-04-04 11:20 by w_xuqing
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
[考研] 274求调剂 +9 顺理成张 2026-04-03 10/500 2026-04-03 15:10 by 啊俊！
[考研] 工科341分调剂 +3 洛多罗 2026-04-03 3/150 2026-04-03 14:20 by 1753564080
[考研] 工科 267求调剂 +5 wanwan00 2026-04-02 7/350 2026-04-03 14:14 by zhangdingwa
[考研] 求调剂！生物与医药专硕 +4 逆转陆先生 2026-04-01 5/250 2026-04-03 08:33 by Jaylen.
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-04-02 3/150 2026-04-02 15:06 by cal0306
[考研] 0703一志愿南师大334求调剂 +4 seven7yu 2026-03-30 4/200 2026-04-01 16:10 by oooqiao
[考研] 0855机械初试280求调剂 +3 kazenotori 2026-03-31 3/150 2026-04-01 10:08 by JourneyLucky
[考研] 301求调剂 +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by 我的船我的海
[考研] 土木304求调剂 +5 顶级擦擦 2026-03-31 5/250 2026-04-01 08:15 by fdcxdystjk￥
信息提示
关闭
请填处理意见
关闭