24小时热门版块排行榜

返回列表

拍拍熊出隐刀

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 4502.1
散金: 1
红花: 3
帖子: 176
在线: 256.5小时
虫号: 1238512
注册: 2011-03-19
性别: GG
专业: 凝聚态物性I:结构、力学和

[求助] 如何计算材料的极化强度

在网上查了很多资料，大多数都是介绍计算Berry有效电荷，很少关于材料本身的自发极化强度的。
自己感觉如果是只计算材料自身自发极化强度的话，就只需要计算平衡态时的极化，然后对G1、G2、G3三个方向的电子项求平均，在加上三个方向上的离子项和berry项，最后除以体积就可以了。
但是尝试算了一下NaNO2的自发极化强度，跟文献报道的计算结果差距很大啊，求大神帮忙看看。并能否告知计算bulk材料自身极化强度的准确方法（非berry有效电荷）。

关于NaNO2的极化强度，文献报道是：采用实验晶格常数体积固定为104 A时，极化0.116 C/m2 ，  J. Phys. Soc. Jpn., Vol. 71, No. 1, January, 2002
我自己的计算结果，自己做了结构优化，体积为111.4 A。
G1：
                              e<r>_ev=(    0.00000 -0.00553 -0.00001 ) e*Angst
                              e<r>_bp=(    0.00000    0.00000    0.00000 ) e*Angst
Total electronic dipole moment: p[elc]=(    0.00000 -0.00553 -0.00001 ) e*Angst
         ionic dipole moment: p[ion]=( 43.70475 -0.10591 32.79778 ) e*Angst

G2：
                              e<r>_ev=(    0.00000 -0.00553    0.00000 ) e*Angst
                              e<r>_bp=(    0.00000 -0.62608    0.00000 ) e*Angst
Total electronic dipole moment: p[elc]=(    0.00000 -0.63161    0.00000 ) e*Angst
         ionic dipole moment: p[ion]=( 43.70475 -0.10591 32.79778 ) e*Angst

G3：
                              e<r>_ev=( -0.00002 -0.00553    0.00000 ) e*Angst
                              e<r>_bp=(    0.00000    0.00000    0.00000 ) e*Angst
Total electronic dipole moment: p[elc]=( -0.00002 -0.00553    0.00000 ) e*Angst
         ionic dipole moment: p[ion]=( 43.70475 -0.10591 32.79778 ) e*Angst

回复此楼