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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 请问LCALCPOL和LBERRY计算的极化有什么区别
201/2返回列表12下一页
查看: 4696  |  回复: 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

木虫 (著名写手)
★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-07-08 09:25:30
看了下tag，似乎lcalcpol是自动把lberry的所有工作给做了。
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
11楼2014-07-05 22:40:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
12楼: Originally posted by jk射手座 at 2014-07-06 08:31:42
恩恩，其实，我开始用的就是LCALC，尽管BERRY很早就听过，但一直没实践过，我的vasp是5.2的。最近试过BERRY，但是结果没提取出来不知道怎么搞得，LCALC的结果 grep  'moment' 一下就出来了...

有那个应该不用berry了吧。
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
13楼2014-07-06 09:09:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-07-08 09:25:42
引用回帖:
14楼: Originally posted by jk射手座 at 2014-07-06 10:44:34
嗯是的，我还想确认下几个事情，这个tag 是否属于密度泛函 微扰理论的；；如果算应变产生的极化怎么加应变上去，我直接按应变量改了晶格常数c，a没变，不知道是否正确。。其余的计算极化值、ezz和BEC就可以一路算下 ...

不属于DFPT，所以你还是需要铁电相-参考相来得到宏观极化
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
15楼2014-07-06 22:12:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-07-08 09:24:54
引用回帖:
16楼: Originally posted by jk射手座 at 2014-07-07 09:23:45
了解了师兄；那DFPT只有那个LEPSILON了吧，请问你用过 LEPSILON 或者LCALCEPS 计算BEC 吗，这两种方法是否是现在很常用的计算BEC的方法，我在论坛里没搜到太多有用的信息。。...

用过LEPSILON，用这个算BEC就行了。
赞一下(1人) 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
17楼2014-07-07 11:34:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
19楼: Originally posted by jk射手座 at 2014-07-08 10:06:48
师兄你回复的不是应助贴，我没法评分呀！那我再问个问题吧，LEPSILON算完之后是怎么提取结果的呀，这个还没用过！...

不用应助了。
输出结果里面有BEC张量
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

jk射手座
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
20楼2014-07-08 11:41:03
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
25楼: Originally posted by jk射手座 at 2014-08-19 10:57:46
师兄，我今天翻了下我的话题偶然发现上次漏看了“参考相”， 我结合我的输入文件，想了一下，是不是说我通过改变晶格常数c来体现应变后的相有一个POSCAR， 原来的相有一个POSCAR， 也就是 POSCAR1 POSCAR2 ？ 我之 ...

参考相和铁电相最好能够保证晶格参数一致
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
27楼2014-08-19 16:13:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
28楼: Originally posted by jk射手座 at 2014-08-20 09:08:18
额有一点没太懂，师兄，晶格参数一致是说a 和c 都要一致吗？ 但是因为我算的是应变情况下的极化行为，所以只有 应变ε=0的时候，POSCAR1 POSCAR2的 里的ac才会是一样的，如果要施加应变，我是保持a不变，只改变c， ...

a和c要一致，应变的情况下就生造一个有应变的顺电相
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
29楼2014-08-20 12:00:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
★ ★ ★
jk射手座(fzx2008代发): 金币+3, 谢谢指导 2014-09-21 19:33:17
引用回帖:
30楼: Originally posted by jk射手座 at 2014-08-20 23:05:57
额。。还是没太懂，能以这个举个例子吗师兄
data_PKT
1.0
3.95 0.00 0.00
0.00 3.95 0.00
0.00 0.00 8.12
1 1 1 1 6
Selective dynamics
………………
比方说这个是没应变的，有应变的顺电相怎么 ...

铁电相：
1.0
3.95 0.0 0.0
0.0 3.95 0.0
0.0 0.0 4.0

Ba 0.0 0.0 0.0
Ti 0.5 0.5 0.51
……
参考相：
1.0
3.95 0.0 0.0
0.0 3.95 0.0
0.0 0.0 4.0

Ba 0.0 0.0 0.0
Ti 0.5 0.5 0.5
……
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
31楼2014-08-23 15:27:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
32楼: Originally posted by jk射手座 at 2014-08-24 09:18:11
哦师兄我明白你的意思了，和英文手册里的那个例子一样，移动相应原子的坐标是吗？但是这样可以提现出应变吗...

可以
回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
33楼2014-08-24 23:24:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamikaruk

木虫 (著名写手)
引用回帖:
34楼: Originally posted by jk射手座 at 2014-08-26 09:00:01
额。。还是没太想懂，是怎么体现的呢。。晶格常数ac不变不是没有应力吗？...

参考相和铁电相同样应变，但是就是原子位置不一样。
赞一下 回复此楼
[img]http://osu.ppy.sh/stat2/iamikaruk-0.png[/img]
35楼2014-08-26 17:10:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jk射手座 的主题更新
201/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭