努力学习英语
用berry pha方法计算BiFeO3极化强度的方法（总的思路是：算出顺电相时（即铁电畸变前的berry pha 和铁电相时的berry pha，然后二者相减）
第一步：对铁电相时的R3c结构的BiFeO3进行结构优化，在结构优化的基础上进行静态计算，以便得到CHGCAR 和 WAVECAR 文件。
第二步：计算R3c结构时的berry pha。（此步要用到第一步静态时得到的CHGCAR 和 WAVECAR文件， 如要计算x方向时的berry pha时INCAR中的部分关键设置为 ISTART=1; ICHARG= 11; ISMEAR= -5; LBERRY = TRUE;LGPAR=1;NPPSTR=8;DIPOL=0.125 0.125 0.125 ; ISPIN=2），同样的方法计算Y 和Z方向的berry pha（设置基本与x是相同，只是将LGPAR分别设为2 和3）。
第三步：将第一步优化后得到的晶格常数a b c α等值作为BiFeO3顺电相R-3C时的初始晶格常数，进行结构优化，以便得到顺电相R-3c时的CHGCAR 和WAVECAR 。
第四步：计算顺电相R-3C时的berry pha（方法同第二步）
第五步：以上步骤一共可以得到六个berry pha 的outcar文件（铁电相R3C时x y z三分方向
outcar及顺电相x y z三个方向的outcar）。分别为：
R3C--X
Total electronic dipole moment: p[elc]=( 2.38612 -0.03401 -0.02221 ) e*Angst清凉的近义词
ionic dipole moment: p[ion]=( 39.24057 22.26457 15.64253 ) e*Angst
R3C—Y
Total electronic dipole moment: p[elc]=( 1.19475 2.06575 -0.02221 ) e*Angst局域网ip地址
ionic dipole moment: p[ion]=( 39.24057 22.26457 15.64253彩虹谷 ) e*Angst
R3C—Z
Total electronic dipole moment: p[elc]=( 1.19475 0.67947 1.95093 ) e*Angst
ionic dipole moment: p[ion]=( 39.24057女生网名唯美小清新 22.26457 15.64253 ) e*Angst
R-3C—x
Total electronic dipole moment: p[elc]=( 2.83940 0.00000 0.00000 ) e*Angst
ionic dipole moment: p[ion]=( 43.08737 24.48986 17.21690 ) e*Angst小儿消积止咳口服液
R-3C—y
Total electronic dipole moment: p[elc]=( 1.45279 2.43958 0.00000 ) e*Angst
游戏切水果 ionic dipole moment: p[ion]=( 43.08737 24.48986 17.21690 ) e*Angst
R-3C—Z
Total electronic dipole moment: p[elc]=( -1.45279 -0.82573 -2.29559 ) e*Angst
肠胃感冒症状 ionic dipole moment: p[ion]=( 43.08737 24.48986 17.21690 ) e*Angst
第六步： ΔP=R3C – R-3C = (122.49-132.1 69.50-75.08 48.83-49.35)
=(-0.96 -6.3 -0.52)
最后总的极化强度= 将上面三个方向的值平方相加在开方，即-0.96 -6.3 -0.52三个数的平方和再开方即为最后的极化强度
