小木虫

当前位置：首页 > 第一原理 >请教关于电极化的数据处理，从波恩有效电荷怎么得到宏观极化值

请教关于电极化的数据处理，从波恩有效电荷怎么得到宏观极化值

作者 wd12309812: 来源: 小木虫 2200 44 举报帖子

通过DFPT计算得到波恩有效电荷，继续求电极化，我计算公式[latex]\sum_{j\beta} Z_{\alpha \beta} \Delta u_{j\beta}[/latex]得到的数值是1.90，体积是115A^3，想请问极化值应该是多少。我看网上有人说应该是1.90*16/115=0.26，这时单位是多少呢？文献中给出的数值是87.3uC/cm2我不知道我是不是计算有误

求助大家返回小木虫查看更多

今日热帖

精华评论

wd12309812

引用回帖:
7楼: Originally posted by tsglss at 2018-01-21 18:32:44
BiFeO3优化后的晶格常数多大？

请问这张图里，每种原子的BECs是怎么得到的呢？
您看这样对不对
Z_xx+Z_yx+Z_zx=Zx
Z_xy+Z_yy+Z_zy=Zy
Z_xz+Z_yz+Z_zz=Zz
sqrt(Zx^2+Zy^2+Zz^2)=Z*

屏幕快照 2018-01-21 下午2.29.30.png
tsglss
tsglss
wd12309812

引用回帖:
10楼: Originally posted by tsglss at 2018-01-21 21:56:04
好像不对，应对矩阵对角化。请问为什么平方和再开方呢？又不是矢量。
...

谢谢回复。能不能说得详细一点呢
矩阵对角化之后，怎么得到这个原子的BECs?
wd12309812

引用回帖:
11楼: Originally posted by tsglss at 2018-01-21 22:05:22
计算波恩有效电荷有专门的公式的（虽然是张量）。在vasp 输出文件里有没有对每个原子对应的有效电荷呢？你在找一找。
...

有的，VASP输出的波恩有效电荷是：
ion 1
1    4.33126 -0.28082 -0.23113
2 -0.29921    4.63435 -0.08390
3 -0.18631 -0.14509    4.77567
ion 2
1    4.33126 -0.31291 -0.18538
2 -0.25887    4.63604 -0.15157
3 -0.22042 -0.08472    4.77567
ion 3
1    3.89460 -0.23805 -0.36669
2 -0.45738    4.26339    0.09968
3 -0.08856 -0.37214    4.37951
ion 4
1    3.89460 -0.42580 -0.09905
2 -0.21371    4.36077 -0.35700
3 -0.36448    0.11617    4.37951
ion 5
1 -2.68144 -0.04598    0.60679
2 -0.03188 -3.39386 -0.04688
3    0.68150 -0.27430 -2.70427
ion 6
1 -3.24044    0.37394    0.19789
2    0.21091 -2.30640    0.54407
3    0.23327    0.35124 -3.23048
ion 7
1 -2.30398    0.41202 -0.52204
2    0.33333 -3.24990    0.01468
3 -0.33265 -0.10346 -3.22043
ion 8
1 -2.30398 -0.35064    0.56511
2 -0.20305 -3.24591 -0.10404
3    0.42124    0.01857 -3.23048
ion 9
1 -2.68143    0.55490 -0.24976
2    0.63481 -2.87793 -0.46257
3 -0.26027 -0.23156 -3.22043
ion 10
1 -3.24044    0.31334    0.28427
2    0.28506 -2.82056    0.54754
3    0.11669    0.72530 -2.70426
，
tsglss

再问，你算出来的偶极距（1.9）是软件自动给算出来的？还是你自己手动计算的？另外，VASP输出的波恩有效电荷分别对应的是哪个离子，你清楚不？
wd12309812

每个原子位移：
ion1 0.00616278 0.006125918 0.006198892
ion2 0.00616084 0.00612622 0.006198639
ion3 0.028559795 0.028526006 0.028595215
ion4 0.028559976 0.028524817 0.02859425
ion5 -0.004131282 0.03419816 -0.139840591
ion6 0.034232197 -0.139908205 -0.004103957
ion7 -0.139872621 -0.00416877 0.034264675
ion8 -0.139875187 0.03419972 -0.004101392
ion9 -0.004136411 -0.139908122 0.034266009
ion10 0.034230577 -0.00416658 -0.13984008

麻烦您帮我看一下，我是在铁电相优化后的晶胞中，直接放入中心对称相的原子坐标，测量每个原子在两种相中笛卡尔坐标之差得到的。这样计算没有问题吧……
现在一个问题是宏观极化值不一致，另一个问题是每个原子的BECs不会计算
麻烦您了！如果方便的话能留一个QQ号吗，谢谢！！
wd12309812

引用回帖:
14楼: Originally posted by tsglss at 2018-01-21 23:32:11
再问，你算出来的偶极距（1.9）是软件自动给算出来的？还是你自己手动计算的？另外，VASP输出的波恩有效电荷分别对应的是哪个离子，你清楚不？

1.9我自己手动算的，ion1 ion2是Bi，ion3 ion4是Fe，余下6个是O