24小时热门版块排行榜

>论坛更新日志 (2920)
>虫友互识 (231)
>导师招生 (172)
>文献求助 (171)
>论文投稿 (68)
>休闲灌水 (52)
>硕博家园 (50)
>公派出国 (48)
>考博 (40)
>招聘信息布告栏 (32)
>论文道贺祈福 (27)
>催化 (25)
>教师之家 (24)
>找工作 (24)
>考研 (22)
>博后之家 (21)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jk射手座

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 9882.6
散金: 1955
红花: 2
帖子: 2281
在线: 175.2小时
虫号: 2777009
注册: 2013-11-04
性别: GG
专业: 零件成形制造

[求助] 请问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就可以了吧？

问得比较多，最近的问题导致开始怀疑计算这一系列的方法的正确性，所以希望看到的路过的版主、大神们，当然还有偶尔出没的侯老师帮解下疑！小弟万分感谢啊！

回复此楼

» 猜你喜欢

求助:我三月中下旬出站，青基依托单位怎么办？已经有8人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
磺酰氟产物，毕不了业了！已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

主做电弧增材制造方向

1楼 2014-06-24 15:19:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jk射手座

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 9882.6
散金: 1955
红花: 2
帖子: 2281
在线: 175.2小时
虫号: 2777009
注册: 2013-11-04
性别: GG
专业: 零件成形制造

引用回帖:

15楼: Originally posted by iamikaruk at 2014-07-06 22:12:14
不属于DFPT，所以你还是需要铁电相-参考相来得到宏观极化...

师兄，我今天翻了下我的话题偶然发现上次漏看了“参考相”，我结合我的输入文件，想了一下，是不是说我通过改变晶格常数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

想跟师兄确认一下这样是否正确可行