24小时热门版块排行榜    

【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 13215  |  回复: 9

焦。。。耳

新虫 (初入文坛)

[交流] VASP计算极化强度已有8人参与

在网上看到的,就复制下来了
Born有效电荷可以通过移动体系原子,然后计算极化强度差,然后除以位移乘以体积得到,或者通过加电场,计算体系原子受力得到;可以用你写的公式计算体系的极化强度,但是前提是你要先找一个参考相,就是所谓的中心对称相,通过比较获得原子位移,然后计算极化,而一般来讲,参考相的选取不是唯一的,不同的参考相会使你计算得到的polarization略有差别,但是差别应该不会很大,假如你选的参考相不太离谱,这是因为born电荷有个求和规则制约,也就是说角标相同的分量加和到一起应该为零。
berry phase的计算同样要选参考相,因为berry phase的计算结果中含有一个任意相,这会使得假如你计算中心对称的体系,都可能计算出电极化,所以如果你用berry phase,你要扣除体系中的量子数,严格的方法也是要找参考相,从参考相变化到铁电相,一点一点变,然后把每个变化的电极化都计算了,扣除量子数,得到铁电相极化!

Berryphase 计算出的单位是e.A,这是电偶极矩的单位,极化强度得用这个值除以体积,注意量纲的转化,经过简化后,你就用Berryphase计算的数值乘以16再除以体积得到单位是C/m2。如果想用uC/cm2,自己可以在转化,这很简单了。


如何用vasp计算Born有效电荷

VASP是采用的Berry phase方法来计算半导体和绝缘体材料的Born有效电荷。这种方法的介绍可以参考VASP手册上提到的文献。

如VASP手册上的介绍,在采用Berry phase方法来计算某个原子某个方向的Born有效电荷时,有两大步:
1) 计算离子未移动时(即平衡状态时)的极化(polarization); 2)手动移动所要计算的原子的位置,计算由原子偏离平衡态时的极化。每一大步中有包括了四个小步骤: i)自恰计算体系以得到收敛的电荷密度(CHG, CHGCAR文件)和波函数(WAVECAR); ii)计算G1方向上的电子极化; iii)计算G2方向上的电子极化; iv)计算G3方向上的电子极化。
最后从每一小步中的OUTCAR文件中找出极化各个部分的贡献: ev (电子项), bp (Berry phase项), ion (离子项),进行简单的数值加减计算,按差分方法计算得到Born有效电荷。 对每类原子的各个方向按同样的步骤进行计算处理,在计算时可根据晶体的对称性,只需要计算Born有效电荷张量的某些分量即可。

下面以计算AlAs的As原子的Born有效电荷为例来介绍具体操作步骤。
1、计算AlAs平衡态时的极化:
i)自恰计算平衡态时的AlAs,主要的输入文件:
INCAR:
SYSTEM = AlAs
ENCUT = 300
ISTART = 0
ICHARG = 2
ISMEAR = 0; SIGMA = 0.002
EDIFF = 1E-5
PREC = Accurate
(注 释: 这里采用ISMEAR=0; SIGMA=0.002,是因为所计算的体系是半导体或绝缘体,电子的占有数实际上固定的,非0就是2或1了,通过将ISMAR=0(即Gaussian 方法来确定电子占有数)且SIGMA设置为一个很小的值来达到一种类似PWSCF程序中的occupation-fixed 方法的效果。如果参看VASP的源代码,可以看到:代码里面会根据INCAR中设置了LBERRY = .TRUE.而自动将ISMEAR和SIGMA进行调整为默认的值:ISMEAR = 0; SIGMA = 0.00)。

KPOINTS:
4x4x4
0
Monma
4 4 4
0 0 0

POSCAR:
AlAs
5.57223
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1 1
Direct
0.000 0.0000 0.0000
0.250 0.2500 0.2500


POTCAR,采用的是LDA的ultrasoft pseudopotential。
自恰计算完之后将得到的CHGCAR, CHG和WAVECAR保存。

ii)计算沿着G1方向上的电子极化:
输入文件如下:
INCAR:
SYSTEM = AlAs
ENCUT = 300
#ISTART = 0
#ICHARG = 2
ISMEAR = 0; SIGMA = 0.002
EDIFF = 1E-5
PREC = Accurate
LBERRY = .TRUE.
IGPAR = 1
NPPSTR = 8
DIPOL = 0.5 0.5 0.5

(注 释:要采用Berry phase方法来计算电子极化时,需在INCAR中设置LBERRY=.TRUE.,另外设置所要计算的方向即设置IGPAR (可赋的值为1,2,3分别表示G1, G2和G3方向), 沿着IGPAR方向上的一串k点的个数 即设置NPPSTR(注意的是,在平衡态时的计算和原子移动之后的计算中,它值应该是一样的),最后设置在计算离子的dipole时的参考点即设置 DIPOL(注意的是,它的设置需要使得原子移动前后的原子都在这个参考点的一侧。比如这个例子中 Al处于(0,0,0),As处于(0.25, 0.25, 0.25)位置,而将DIPOL设置为(0.5, 0.5, 0.5)和(0.125, 0.125, 0.125)都是可以的,但是在考虑移动Al原子时,不要将原子移动原胞之外即偏移量为负数;另外也不要将DIPOL设置在所要移动的原子上,如果是这样的话,则会导致移动该原子后,该原子不在DIPOL的同一侧,使得原子移动之后的极化的Berry-phase项要比平衡态时的大很多。)

POSCAR和KPOINTS同i)中的。
将第i)计算得到的CHGCAR, CHG, WAVECAR拷贝过来,进行计算。
计算完后保存OUTCAR,从OUTCAR文件中找出极化的各个项的贡献,可以采用下面的命令:
grep '' OUTCAR

如下面的数据
Expectation value term: 《R》ev
《R》 = ( -0.00022, -0.00117, 0.00074 ) electrons Angst
Berry-Phase term: 《R》bp
《R》 = ( 0.00000, -0.00455, -0.00455 ) electrons Angst
ionic term: 《R》ion
《R》 = ( 15.32363, 15.32363, 15.32363 ) electrons Angst

(注释:后面的三列数据分别表示的是x, y, z方向的。)

iii)和iv)步同ii)类似,只是分别将IGPAR = 1改为IGPAR = 2和IGPAR = 3。
下面是计算G2方向时(即IGPAR=2)得到的数据:
Expectation value term: 《R》ev
《R》 = ( -0.00117, -0.00022, 0.00074 ) electrons Angst
Berry-Phase term: 《R》bp
《R》 = ( -0.00455, 0.00000, -0.00455 ) electrons Angst
ionic term: 《R》ion
《R》 = ( 15.32363, 15.32363, 15.32363 ) electrons Angst

下面是计算G3方向时(即IGPAR=3)得到的数据:
Expectation value term: 《R》ev
《R》 = ( -0.00117, 0.00074, -0.00022 ) electrons Angst
Berry-Phase term: 《R》bp
《R》 = ( -0.00455, -0.00455, 0.00000 ) electrons Angst
ionic term: 《R》ion
《R》 = ( 15.32363, 15.32363, 15.32363 ) electrons Angst

下面按公式计算计平衡态时《R》ev在G1、G2和G3方向上的平均值:
《R》 ev, undistorted, average =
[( -0.00022, -0.00117, 0.00074 )+
( -0.00117, -0.00022, 0.00074 )+
( -0.00117, 0.00074, -0.00022 ) ]/3.0
= ( -0.00022-0.00117-0.00117, -0.00117-0.00022+0.00074, 0.00074+0.00074-0.00022 ) /3.0
=(-0.000853333, -0.000216667, 0.00042)

把平衡态时各个方向计算时得到《R》bp项求和:
《R》bp, undistorted, all =
( 0.00000, -0.00455, -0.00455 ) +
( -0.00455, 0.00000, -0.00455 ) +
( -0.00455, -0.00455, 0.00000 )
= ( -0.00455*2, -0.00455*2, -0.00455*2 )
= ( -0.0091, -0.0091, -0.0091)

那么平衡态时电子极化项 《R》el, undistorted就是:
《R》 el, undistorted =《R》 ev, undistorted, average + 《R》bp, undistorted, all
=(-0.000853333, -0.000216667, 0.00042) +( -0.0091, -0.0091, -0.0091)
=(-0.000853333 -0.0091, -0.000216667 -0.0091, 0.00042 -0.0091)
=(-0.00995333, -0.00931667, -0.00868)

平衡态时离子极化项《R》ion, undistorted,在任意一个方向计算时得到的数据,即:
《R》ion, undistorted = ( 15.32363, 15.32363, 15.32363 )

2. 计算Al原子偏离平衡位置时的极化:
i)其中的INCAR, KPOINTS同第1 i)步中的,只是将POSCAR中 As原子的位置沿着某个方向略移动一个小量,如下:
POSCAR
AlAs
5.57223
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1 1
Direct
0.000 0.0000 0.0000
0.250 0.2500 0.2300
(注释:这种设置是导致As在x和y方向上都移动了0.05572 angstrom,这个可以自己算出来,或者直接从移动前后两次计算的OUTCAR中的原子cartesian坐标进行简单的加减计算之后得到。)

然后作自恰计算,计算之后保存得到CHGCAR, CHG, WAVECAR。

ii)、计算原子移动之后,体系沿着G1方向的极化:
INCAR和KPOINTS的输入文件同第1 ii)步中的,将上一步保存下来的CHG, CHCAR和WAVECAR拷贝过来,然后计算得到OUTCAR,在保存后,从中可以找到相关的极化计算数据,通过下面的命令来查找:
grep '' OUTCAR

此例子的计算数据:
Expectation value term: 《R》ev
《R》 = ( -0.00091, -0.00100, -0.00030 ) electrons Angst
Berry-Phase term: 《R》bp
《R》 = ( 0.00000, -0.00475, -0.00475 ) electrons Angst
ionic term: 《R》ion
《R》= ( 15.60224, 15.60224, 15.32363 ) electrons Angst

iii)& iv).计算原子移动后体系沿着G2和G3方向的极化:
INCAR和KPOINTS中的分别与第1 iii)&iv)步中的一样。
此例子计算得到的数据:
沿着G2方向的:
Expectation value term: 《R》ev
《R》 = ( -0.00100, -0.00091, -0.00030 ) electrons Angst
Berry-Phase term: 《R》bp
《R》= ( -0.00475, 0.00000, -0.00475 ) electrons Angst
ionic term: 《R》ion
《R》 = ( 15.60224, 15.60224, 15.32363 ) electrons Angst

沿着G3方向的:
Expectation value term: 《R》ev
《R》 = ( -0.00233, 0.00032, -0.00029 ) electrons Angst
Berry-Phase term: 《R》bp
《R》 = ( -0.40026, -0.40026, 0.00000 ) electrons Angst
ionic term: 《R》ion
《R》 = ( 15.60224, 15.60224, 15.32363 ) electrons Angst

同样按公式计算出原子移动之后的极化:
将G1, G2, G3各个方向计算得到《R》ev项进行平均:
《R》ev, distorted, average =
[ ( -0.00091, -0.00100, -0.00030 ) +
( -0.00100, -0.00091, -0.00030 ) +
( -0.00233, 0.00032, -0.00029 )]/3.0
= ( -0.00091-0.00100-0.00233, -0.00100-0.00091+0.00032, -0.00030-0.00030--0.00029)/3.0
=(-0.00141333, -0.00053, -0.000103333)
将G1, G2, G3各个方向计算时得到《R》bp项求和:
《R》bp, distorted, all =
( 0.00000, -0.00475, -0.00475 )+
( -0.00475, 0.00000, -0.00475 )+
( -0.40026, -0.40026, 0.00000 )
=(0.00000-0.00475-0.40026, -0.00475+0.00000-0.40026, -0.00475 -0.00475+0.00000)
=(-0.40501, -0.40501, -0.0095)

由此得到原子移动之后电子极化项 《R》el, distorted:
《R》el, distorted= 《R》ev, distorted, average + 《R》bp, distorted, all
=(-0.00141333, -0.00053, -0.000103333)+(-0.40501, -0.40501, -0.0095)
=(-0.00141333-0.40501, -0.00053-0.40501, -0.000103333-0.0095)
=(-0.406423, -0.40554, -0.00960333)

原子移动后离子极化项《R》ion, distorted取任意一个方向上计算得到的数据《R》ion:
《R》ion, distorted = ( 15.60224, 15.60224, 15.32363 )


由此得到由于As原子的移动引起的极化变化为:
Delta 《R》= 《R》ion, distorted + 《R》el, distorted - 《R》ion, undistorted - 《R》el, undistorted
= ( 15.60224, 15.60224, 15.32363 )+(-0.406423, -0.40554, -0.00960333) - ( 15.32363, 15.32363, 15.32363 )-(-0.00995333, -0.00931667, -0.00868)
=( 15.60224-0.406423-15.32363+0.00995333, 15.60224-0.40554-15.32363+0.00931667, 15.32363-0.00960333-15.32363+0.00868)
=(-0.11786, -0.117613, -0.00092333)

在前面已经提到了这个是在As分别在沿着x, y方向上移动了0.05572 angstrom时的计算结果。那么由此可计算出As原子的沿着x,y方向上的Born有效电荷(实际上由于AlAs是立方晶体,As原子的Born有效电荷沿着xx, yy, zz方向的都相等,xz, yz, xy方向上的都为0)按差分公式可以计算得到:
Z*_As = Delta《R》/ Delta x_As
=-0.11786/0.05572
=-2.11522

由于原胞所有原子的Born有效电荷各个方向的加起来要满足电荷中性的条件,AlAs中只有Al和As两类(也只有2个原子)原子,因此也可以得到Al的Born有效电荷为2.11522,这与ABINIT中的教程中计算的数据很接近。

实际上面两大步的4个小步骤,可采下面的脚本来连续计算完成。

要注意的就是DIPOL 的设置,以及原子移动的选取:不要将DIPOL设置在某个原子位置上,原子移动的选取要保证原子在移动前后都是DIPOL的同一侧。
回复此楼

» 收录本帖的淘贴专辑推荐

计算

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

刘想许婧

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你好,请问你有试过lcalcpol计算极化吗?你第二个POSCAR在移动原子时,如果有多个原子,是不是每个原子都沿同一个方向移动相同的距离?
2楼2018-05-18 09:44:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

持之以恒sc

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你好,DIPOL得设置,有什么要规则呢?

发自小木虫Android客户端
3楼2018-10-30 12:34:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

持之以恒sc

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你好,请问移动前后原子在DIPOL同一侧,这个原子前后的移动,指的是怎样的的移动呢?是优化你过程中移动的前后吗?

发自小木虫Android客户端
4楼2018-10-30 15:53:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

W.Z.Zhou

捐助贵宾 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
您好,最近我想计算体系的极化强度,没有取参考相,计算出来的不符合常识
请问,扣除体系中的量子数是指什么,量子数怎么确定?参考相应该怎么选取呢,比如CO分子,是取C和O不同的距离吗?
希望你能帮助一下我,万分感谢!
幸福也罢,苦难也罢,都来自外界,坚毅则来自内心,来自一个人自身的努力。
5楼2019-08-05 11:05:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

简单才是美

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by W.Z.Zhou at 2019-08-05 11:05:52
您好,最近我想计算体系的极化强度,没有取参考相,计算出来的不符合常识
请问,扣除体系中的量子数是指什么,量子数怎么确定?参考相应该怎么选取呢,比如CO分子,是取C和O不同的距离吗?
希望你能帮助一下我,万 ...

同学,现在会计算极化强度了吗,我最近也要算这个,还是有些疑惑,能不能留个联系方式,我们交流一下

发自小木虫Android客户端
6楼2019-10-24 16:09:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

采采芣苡

木虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by 简单才是美 at 2019-10-24 16:09:53
同学,现在会计算极化强度了吗,我最近也要算这个,还是有些疑惑,能不能留个联系方式,我们交流一下
...

同学,你现在会算极化强度了吗?我也在弄,有些不解,能不能加个qq,交流一下?,2303391842
7楼2021-01-28 10:26:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nuli的菜

银虫 (著名写手)

8楼2021-01-29 21:02:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

warm_sake

新虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
根据那个坐标怎么求的晶格常数啊
9楼2021-02-27 15:03:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

源来汝瓷

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2楼: Originally posted by 刘想许婧 at 2018-05-18 09:44:06
你好,请问你有试过lcalcpol计算极化吗?你第二个POSCAR在移动原子时,如果有多个原子,是不是每个原子都沿同一个方向移动相同的距离?

你好,请问你怎么解决这个问题的,可以回复一下吗
10楼2021-04-19 09:48:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 焦。。。耳 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境学硕287求调剂 5+3 上岸1111144 2024-03-28 6/300 2024-03-29 18:23 by 小damn
[考研] 化学工程与技术总分320求调剂 +8 Wjrjxmt 2024-03-28 9/450 2024-03-29 16:27 by Wk202212
[教师之家] 冷漠的大学 +17 otani 2024-03-27 17/850 2024-03-29 16:21 by sjtu2012
[高分子] 透明聚脲涂层 +3 顺利毕业的小白 2024-03-25 6/300 2024-03-29 13:59 by softwindy
[考研] 化学和化工专业调剂 +3 我爱真心 2024-03-28 5/250 2024-03-29 13:14 by 许昌666
[基金申请] MDPI的期刊 +5 zju2000 2024-03-26 5/250 2024-03-28 23:00 by zzahkj
[考研] 310求调剂 +4 cxr2024 2024-03-28 5/250 2024-03-28 22:52 by ms629
[考研] 考研304化工求调剂 +8 abc123456* 2024-03-27 10/500 2024-03-28 21:59 by qlm5820
[基金申请] 今年的国社科报名通知怎么还不发布? +6 屡战屡败 2024-03-27 6/300 2024-03-28 20:11 by 风声鹤草
[考研] 总分323,一志愿华东理工 +10 努_力66 2024-03-23 10/500 2024-03-28 20:10 by 麦田里的守望者
[考研] 302求调剂 +6 范智乔ε 2024-03-26 6/300 2024-03-28 14:55 by qlm5820
[考研] 一志愿武汉理工材料学,学硕,300分求调剂 +9 想向上pao 2024-03-26 11/550 2024-03-28 14:09 by w7238495
[教师之家] 刚入职的青年教师 +4 ¥笑傲江湖¥ 2024-03-26 6/300 2024-03-28 12:52 by ¥笑傲江湖¥
[考研] 环境科学与工程初试303分求调剂 +4 来瓣儿柚子 2024-03-26 4/200 2024-03-28 11:04 by 小雨哥
[有机交流] 溴素溴化反应生成溴苄 100+4 tanjishuang 2024-03-25 6/300 2024-03-28 09:37 by 高树桩
[考研] 312求调剂 +6 牛肉丸敬告 2024-03-22 6/300 2024-03-27 16:40 by tkdmxl
[考研] 化学工程与技术320求调剂 +7 Wjrjxmt 2024-03-26 8/400 2024-03-27 08:50 by 实验室@
[考研] 09食加 285A区求调剂 211本科 +3 心之漾漾呀 2024-03-26 4/200 2024-03-27 01:47 by 心之漾漾呀
[考研] 335求调剂 +6 云怎么挽留 2024-03-26 7/350 2024-03-26 18:29 by Whisper890
[考研] 一志愿北京工业大学,材料工程,总分301求调剂 +9 hss0418 2024-03-24 14/700 2024-03-25 20:28 by EBSD
信息提示
请填处理意见