24小时热门版块排行榜

返回列表

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

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就可以了吧？

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

回复此楼

» 猜你喜欢

2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
孩子确诊有中度注意力缺陷已经有6人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

主做电弧增材制造方向

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

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