| 查看: 5593 | 回复: 55 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
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就可以了吧？ 问得比较多，最近的问题导致开始怀疑计算这一系列的方法的正确性，所以希望看到的路过的版主、大神们，当然还有偶尔出没的侯老师帮解下疑！小弟万分感谢啊！ |
回复此楼» 猜你喜欢
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有12人回复
-
极限拉扯，如何抉择
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
jk射手座
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 9882.6
- 散金: 1955
- 红花: 2
- 帖子: 2281
- 在线: 175.2小时
- 虫号: 2777009
- 注册: 2013-11-04
- 性别: GG
- 专业: 零件成形制造
|
师兄，我今天翻了下我的话题偶然发现上次漏看了“参考相”， 我结合我的输入文件，想了一下，是不是说我通过改变晶格常数c来体现应变后的相有一个POSCAR， 原来的相有一个POSCAR， 也就是 POSCAR1 POSCAR2 ？ 我之前用LCALCPOL 一直是拿脚本算这两个POSCAR的，我给你看下，举个例子： POSCAR1：（这个是PKT超晶格没受应力时的晶格参数和坐标POSCAR） 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 Direct 0.00000 0.00000 0.00000 F F T 0.50000 0.50000 0.25000 F F T 0.00000 0.00000 0.50000 F F T 0.50000 0.50000 0.75000 F F T 0.50000 0.00000 0.25000 F F T 0.00000 0.50000 0.25000 F F T 0.50000 0.50000 0.00000 F F T 0.50000 0.00000 0.75000 F F T 0.00000 0.50000 0.75000 F F T 0.50000 0.50000 0.50000 F F T POSCAR2：（这个是PKT发生 -1.5%的应变时的晶格参数，也就是按照(c-c0)/c0=ε得到的c=7.9982原子坐标我用的是relax之后的原子坐标）： data_PKT 1.00000000000000 3.9500000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 3.9500000000000001 0.0000000000000000 0.0000000000000000 0.0000000000000000 7.9982000000000000 Pb Ti K Ta O 1 1 1 1 6 Selective dynamics Direct 0.0000000000000000 0.0000000000000000 -0.0444396041984977 F F T 0.5000000000000000 0.5000000000000000 0.2327315837951796 F F T 0.0000000000000000 0.0000000000000000 0.4971976926301170 F F T 0.5000000000000000 0.5000000000000000 0.7423875204210301 F F T 0.5000000000000000 0.0000000000000000 0.2778793310366013 F F T 0.0000000000000000 0.5000000000000000 0.2778793310366013 F F T 0.5000000000000000 0.5000000000000000 0.0092994307667890 F F T 0.5000000000000000 0.0000000000000000 0.7482258333968913 F F T 0.0000000000000000 0.5000000000000000 0.7482258333968913 F F T 0.5000000000000000 0.5000000000000000 0.5106130477183970 F F T 想跟师兄确认一下这样是否正确可行  |
25楼2014-08-19 10:57:46
jk射手座
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 9882.6
- 散金: 1955
- 红花: 2
- 帖子: 2281
- 在线: 175.2小时
- 虫号: 2777009
- 注册: 2013-11-04
- 性别: GG
- 专业: 零件成形制造
怎么没人做这方面的计算吗？小弟实在物理底子差，看过很多帖子但讲的不是很直白，这块不搞懂做起模拟来很别扭，希望有前辈点拨，金币可以再加，拜托拜托！！！ 2楼2014-06-26 19:11:24
elliot.cc
铁虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 53.5
- 帖子: 49
- 在线: 18.1小时
- 虫号: 2784792
- 注册: 2013-11-07
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2014-06-28 14:02:33
jk射手座
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 9882.6
- 散金: 1955
- 红花: 2
- 帖子: 2281
- 在线: 175.2小时
- 虫号: 2777009
- 注册: 2013-11-04
- 性别: GG
- 专业: 零件成形制造
4楼2014-06-28 16:15:45
