| 查看: 10072 | 回复: 52 | |||||||||
[交流]
Berry phase 计算铁电极化的步骤 confirmed by 大牛Ederer 已有18人参与
|
|||||||||
|
四年之前曾经发帖求助Berry phase 如何计算铁电极化的问题，那个帖子之后有一些同行回帖或者短消息我， 抱歉我都没有回复， 主要是本人对小木虫使用很有限，而且功能不熟悉，每次登陆自己账户都是为了下载资源，之后匆匆下线。 今天为了回馈多年来小木虫上对自己的帮助，特此分享下自己在使用Berry phase方法上的一些经验。 可能大家都知道根据定义，Berry phase 算极化是用铁电相的极化（包含离子项+电子项）减去对应中心对称相的极化， 但是，这个方法过于粗糙，而且会遇到至少两个问题： 1.中心对称相可能不绝缘，导致无法使用Berry phase， 这时会报错：Error in subroutine BERRY: did not find all determinants（可能不是这句话，记不清了）。有时+U可以让体系绝缘，但某些情况再怎么+U都不能绝缘，比如BiFeO3 四方中心对称相P4/mmm； 2.铁电相和中心对称相的离子项差异可能会非常大，导致二者差值也远超过正常的极化值。 这两个问题曾经困扰了我一段时间，后来我发邮件请教了一位大牛Ederer，他传授给我一套正确的计算步奏： 1.首先我们需要构建一系列的中间过渡相。假设铁电相中各个原子，相对于中心对称相的对应原子，位移为100%， 那么需要构建比如原子位移为90%，80%。。。。。10%这一系列的中间过渡相。注意这里中间过渡相建立好之后，原子位置应当固定，不需要再弛豫。 2. 然后计算100%，90%，80%。。。这一系列相的极化，保证DOPOL设置一样。 但在计算中我们可能发现，比如，从50%到0%体系不绝缘了，那么这些相没法计算极化了，也不需要计算。 3.以P(100%)，P(90%), P(80%)......为y轴，百分比为x轴 （参照PHYSICAL REVIEW B 71, 014113）作图，即使中间缺失了从50%到0%的点（假设它们不绝缘无法计算），只要其余的点是有一个很好的trend变化，然后那条曲线一直延伸到原点就好了。 注意，原点处P（0%）不一定是0！而是polarization quantum (eR/V) 的整数(或整数+1/2)倍，所取的P（0%）应当是最为符合那条曲线trend的值。大功告成，极化值就是P（100%）-P（0%）。 4.假如那些点没有一个好trend，中间有的点有跳跃，那么一定是由于polarization quantum 造成的，这时候需要适当加减polarization quantum的倍数，让那些跳跃的点回到那条有trend的曲线上来。 利用此方法，可以顺利求出R3C BiFeO3的极化值，而且从100%的R3c 到 0%的R-3c， 中间所有过渡相都是绝缘。 也可以顺利求出四方相P4mm的极化，但是从100%的P4mm到0%的P/4mmm中间有些相是不绝缘的，那些相的极化值缺失，但是没有关系，需用剩余的点足以得到一条曲线。BaTiO3的极化计算更为简单，我印象中都是所有相都是绝缘的，而且没有polarization quantum 造成的跳跃。 另外附上我的一些论文，都是用此方法计算的，一篇纯计算的发表在Sci. Rep. (本来投的PRB，但跟一个reviewer fight 很久失败，其中原因说多都是泪  ）。另外一篇是计算+实验发表在ACS. Appl. Mater. Interfaces, 还有一篇帮我师兄计算的文章发表在PRB，详情见https://scholar.google.com.sg/ci ... h-CN&oi=sra， 也可以去我的researchgate下载http://www.researchgate.net/profile/Zhen_Fan3。最后祝大家新年发好文！  |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理计算-tg | 第一性原理 | 技术贴 | soft study |
| 计算 | vasp |
» 本帖已获得的红花（最新10朵）
zyl雨田» 猜你喜欢
急需调剂
已经有3人回复
-
22408 312求调剂
已经有21人回复
-
恳请有学校收留
已经有4人回复
-
297，工科调剂?
已经有9人回复
-
300求调剂
已经有4人回复
-
291求调剂
已经有5人回复
-
085404 22408 309分求调剂
已经有10人回复
-
296求调剂
已经有13人回复
-
通信工程求调剂！！！
已经有6人回复
-
294求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 多铁材料铁电极化（Berryphase方法）计算问题
已经有6人回复
|
你点那个MINSUP，然后choose 铁电相，比如R3C, 然后就会列出一堆supergroup：http://www.cryst.ehu.es/cgi-bin/ ... y=&gnum=161， R-3C就在其中。至于哪些相是中心对称的，参看：http://wenku.baidu.com/view/ae42f76fb84ae45c3b358c50 |
» 本帖已获得的红花（最新10朵）
25楼2015-03-17 14:40:47
|
我们这里是已知了BFO,BTO等材料是铁电材料，只是在计算时发现某些中心对称相是不绝缘导致没法算，但文中的方法可以克服这一困难。至于Ga2O3，我觉得你还是先要调研下文献，如果实验上Ga2O3是很导电的，那肯定不是铁电体，因为铁电要求材料绝缘，否则加电压都变成漏电流跑掉了。 如果实验上它是铁电的，在计算时候需要找它对应的中心对称相，可以把Ga2O3的space group输入到http://www.cryst.ehu.es/，根据sub/super-group 的关系去寻找对应中心对称相。 |
» 本帖已获得的红花（最新10朵）
14楼2015-03-03 14:56:52
|
在Berry phase计算方法里面，中心对称结构的极化可以是polarization quantum的整数或整数+1/2倍，所以它可以是零，也可以是一个值，要符合那条有trend的曲线（具体参考那篇PRB）。 |
» 本帖已获得的红花（最新10朵）
34楼2015-03-18 10:47:01
41楼2015-03-20 12:39:43
|
像你这种情况，你就把I41/AMD里面的原子往高对称的地方移动，比如0.5,0.5,0.52,你就移到0.5,0.5,0.5.然后用material studio找对称性。如果还是不行，我也没办法了，你自己再查查文献看看 |
» 本帖已获得的红花（最新10朵）
44楼2015-03-29 08:39:46
47楼2015-05-14 02:20:29
wuli8
荣誉版主 (知名作家)
…………
- 1ST强帖: 2
- 应助: 35 (小学生)
- 贵宾: 12.924
- 金币: 20190.4
- 散金: 15888
- 红花: 88
- 沙发: 4
- 帖子: 7840
- 在线: 1114.6小时
- 虫号: 465889
- 注册: 2007-11-23
- 专业: 物理学I
- 管辖: 计算模拟
2楼2015-02-28 14:42:44
gongjijun
木虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 4161.8
- 散金: 48
- 红花: 5
- 帖子: 507
- 在线: 846小时
- 虫号: 935467
- 注册: 2009-12-30
- 专业: 凝聚态物性 II ：电子结构
3楼2015-03-01 17:15:18
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3524.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
4楼2015-03-02 00:37:27
5楼2015-03-02 09:13:27
6楼2015-03-02 09:17:08
wbcui
至尊木虫 (知名作家)
- 应助: 221 (大学生)
- 贵宾: 0.349
- 金币: 16296.6
- 散金: 359
- 红花: 73
- 帖子: 8321
- 在线: 919小时
- 虫号: 543804
- 注册: 2008-04-11
- 专业: 金属功能材料
7楼2015-03-02 10:00:53
stonezhoujun
金虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 4911.8
- 散金: 613
- 红花: 2
- 帖子: 2726
- 在线: 499.8小时
- 虫号: 256246
- 注册: 2006-06-03
- 性别: GG
- 专业: 无机非金属类电介质与电解
8楼2015-03-02 11:21:06
wangqj1
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 金币: 7724.5
- 帖子: 2164
- 在线: 485.7小时
- 虫号: 485098
- 注册: 2007-12-27
- 专业: 凝聚态物性 II ：电子结构
9楼2015-03-02 14:55:03
10楼2015-03-02 19:23:32
5