24小时热门版块排行榜

返回列表

jk射手座

新虫 (著名写手)

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

[求助] 关于算压电极化强度的结果，求教前辈们如何分析之

最近学习vasp算极化强度，但是结果看得不是很懂，下面我把我的各文件以及结果列出来，请路过的大牛们帮忙看下：
INCAR：
SYSTEM = PKT
ENCUT=500
PREC=Accurate
ISMEAR =0
SIGMA=0.05
EDIFF=1E-6
ISTART=0
ICHARG=2
IBRION=2
LCALCEPS=.TRUE.

KPOINTS:
Monkhorst Pack
0
Monkhorst Pack
2 2 1
0 0 0

POSCAR：
data_PKT
1.00000000000000
   3.9300000000000002 0.0000000000000000 0.0000000000000000
   0.0000000000000000 3.9300000000000002 0.0000000000000000
   0.0000000000000000 0.0000000000000000 7.9299999999999997
Pb Ti K Ta O
   1    1    1    1    6
Selective dynamics
Direct
  0.0000000000000000  0.0000000000000000  0.0365263276214520 F F T
  0.5000000000000000  0.5000000000000000  0.2624523516524320 F F T
  0.0000000000000000  0.0000000000000000  0.4962184331355944 F F T
  0.5000000000000000  0.5000000000000000  0.7536897875049664 F F T
  0.5000000000000000  0.0000000000000000  0.2345480460876364 F F T
  0.0000000000000000  0.5000000000000000  0.2345480460876364 F F T
  0.5000000000000000  0.5000000000000000 -0.0041182693364607 F F T
  0.5000000000000000  0.0000000000000000  0.7458680596303440 F F T
  0.0000000000000000  0.5000000000000000  0.7458680596303440 F F T
  0.5000000000000000  0.5000000000000000  0.4943991579860550 F F T

POTCAR是Pb Ti K Ta O 的PAW赝势 _sv的，按顺序叠加得到的文件

下面是我算出的结果，用grep 'moment' OUTCAR 提取的结果：

Ionic dipole moment: p[ion]=(    0.00000    0.00000 -48.72645 ) electrons Angst
Total electronic dipole moment: p[elc]=(    0.00000    0.00000    4.23883 ) electrons Angst
Ionic dipole moment: p[ion]=(    0.00000    0.00000 -48.72645 ) electrons Angst
Total electronic dipole moment: p[elc]=(    0.03863    0.00000    4.23885 ) electrons Angst
Total change in dipole moment: dp[elc]=(    0.03863    0.00000    0.00002 ) electrons Angst

Ionic dipole moment: p[ion]=(    0.00000    0.00000 -48.72645 ) electrons Angst
Total electronic dipole moment: p[elc]=(    0.00000    0.00000    4.23882 ) electrons Angst
Total change in dipole moment: dp[elc]=(    0.03863    0.00000    0.00002 ) electrons Angst

Ionic dipole moment: p[ion]=(    0.00000    0.00000 -48.72645 ) electrons Angst
Total electronic dipole moment: p[elc]=(    0.00000    0.03863    4.23886 ) electrons Angst
Total change in dipole moment: dp[elc]=(    0.00000    0.03863    0.00003 ) electrons Angst

Ionic dipole moment: p[ion]=(    0.00000    0.00000 -48.72645 ) electrons Angst
Total electronic dipole moment: p[elc]=(    0.00000    0.00000    4.23883 ) electrons Angst
Total change in dipole moment: dp[elc]=(    0.00000    0.03863    0.00003 ) electrons Angst

Ionic dipole moment: p[ion]=(    0.00000    0.00000 -48.72645 ) electrons Angst
Total electronic dipole moment: p[elc]=(    0.00000    0.00000    4.27186 ) electrons Angst
Total change in dipole moment: dp[elc]=(    0.00000    0.00000    0.03303 ) electrons Angst
请问，以上中那个Total change in dipole moment的数据怎么用，前两项我知道是离子和电子的偶极矩，但是Total change in dipole moment的数据是做差还是做和得来的，小弟没太看懂。之前在INCAR中设置过LCALPOL 参数，得到的结果不像这个这么多，只有离子和电子的偶极矩两项，简单做差得到uij算极化强度，这个要怎么算呢，我看手册上还有一个LEPSILON，我想都试试，就是不知道三者算极化有什么区别，结果怎么处理，小弟先在此谢过了！

回复此楼