24小时热门版块排行榜

返回列表

jk射手座

新虫 (著名写手)

应助: 0 (幼儿园)
金币: 9882.6
散金: 1955
红花: 2
帖子: 2281
在线: 175.2小时
虫号: 2777009
注册: 2013-11-04
性别: GG
专业: 零件成形制造

[求助] 求教LCALCPOL计算极化的结果怎么处理的，多出来一行

我先贴下我的INCAR：
SYSTEM=PTOLAO
ENCUT=450
ISMEAR=-5
SIGMA=0.02
EDIFF=1E-3
ISTART=0
ICHARG=2
IBRION=2
PREC=Accurate
LCALCPOL=.TRUE.
DIPOL=0.3 0.3 0.3
IDIPOL=3

LDAU=.TRUE.
LDAUTYPE=2
LDAUL=-1 -1 3 -1 -1
LDAUU=0 0 6 0 0
LDAUJ=0 0 0 0 0
LDAUPRINT=2

K网格是自动  6 6 4 的
POTCAR用的是PAW的
POSCAR有两个，
其一：
data_PTOLAO
1.00000000000000
   3.8100000000000001 0.0000000000000000 0.0000000000000000
   0.0000000000000000 3.8100000000000001 0.0000000000000000
   0.0000000000000000 0.0000000000000000 8.1400000000000000
Pb Ti La Al O
   1    1    1    1    6
Selective dynamics
Direct
0.00  0.00  0.00 F F T
0.50  0.50  0.25 F F T
0.00  0.00  0.50 F F T
0.50  0.50  0.75 F F T
0.50  0.00  0.25 F F T
0.00  0.50  0.25 F F T
0.50  0.50  0.00 F F T
0.50  0.00  0.75 F F T
0.00  0.50  0.75 F F T
0.50  0.50  0.50 F F T

其二：
data_PTOLAO
1.00000000000000
   3.8100000000000001 0.0000000000000000 0.0000000000000000
   0.0000000000000000 3.8100000000000001 0.0000000000000000
   0.0000000000000000 0.0000000000000000 7.9772000000000000
Pb Ti La Al O
   1    1    1    1    6
Selective dynamics
Direct
  0.0000000000000000  0.0000000000000000 -0.0444396041984977 F F T
  0.5000000000000000  0.5000000000000000  0.2327315837951796 F F T
  0.0000000000000000  0.0000000000000000  0.4971976926301170 F F T
  0.5000000000000000  0.5000000000000000  0.7423875204210301 F F T
  0.5000000000000000  0.0000000000000000  0.2778793310366013 F F T
  0.0000000000000000  0.5000000000000000  0.2778793310366013 F F T
  0.5000000000000000  0.5000000000000000  0.0092994307667890 F F T
  0.5000000000000000  0.0000000000000000  0.7482258333968913 F F T
  0.0000000000000000  0.5000000000000000  0.7482258333968913 F F T
  0.5000000000000000  0.5000000000000000  0.5106130477183970 F F T

然后计算出的结果我提取的 'moment'.提取出来是酱婶儿的：

angular momentum for each species LDAUL = -1 -1 3 -1 -1
dipolmoment          0.000000    0.000000 -21.359368 electrons x Angstroem
         Ionic dipole moment: p[ion]=( 12.57300 12.57300 -20.04233 ) electrons Angst
Total electronic dipole moment: p[elc]=(    0.00000    0.00000    0.28231 ) electrons Angst
[root@node16 -1.5calc]# grep 'moment' OUTCAR1
angular momentum for each species LDAUL = -1 -1 3 -1 -1
dipolmoment          0.000000    0.000000 -24.451270 electrons x Angstroem
         Ionic dipole moment: p[ion]=( 12.57300 12.57300 -24.01300 ) electrons Angst
Total electronic dipole moment: p[elc]=(    0.00000    0.00000    0.15750 ) electrons Angst
[root@node16 -1.5calc]# grep 'volume' OUTCAR2
  volume/ion in A,a.u.             =    11.64       78.54
  volume of cell :    116.39
  volume of cell :    116.39
请问 electrons x Angstroem那行的dipolemoment 和下边两行电子与离子的偶极矩是啥关系啊？多设置了一个IDIPOLE=3之后就多了 electrons x Angstroem的那一行；
搜了些相关的说明与讨论，如果我理解为： electrons x Angstroem那行dipolemoment是总P（也就是p[ion]+p[elc]）在Z方向上的分量，这样理解对吗？我计算了下，分别是3.1027和4.20832，那么其与z轴的夹角为arccos(3.1027/4.20832),大概比45°小一点；我算的是两个四方相组合的超晶格，总P的方向应该离体对角线不远，不知这样理解对不？

请各路大神帮忙指点下！！！！
（另外还想问下，LBERRY大家都知道是BERRY PHASE方法，那么设置LCALCPOL算极化，这个方法有名字吗？？）

[ Last edited by jk射手座 on 2014-7-2 at 10:12 ]

回复此楼