24小时热门版块排行榜    

Znn3bq.jpeg
查看: 734  |  回复: 2

jk射手座

新虫 (著名写手)

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

最近学习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,我想都试试,就是不知道三者算极化有什么区别,结果怎么处理,小弟先在此谢过了!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

主做电弧增材制造方向
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)

自己先顶一下,在线等!
主做电弧增材制造方向
2楼2014-02-26 10:10:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jk射手座

新虫 (著名写手)

晕,都科研去了吗,没人呢
主做电弧增材制造方向
3楼2014-02-26 10:56:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jk射手座 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中农微生物,288分,三年实验经历 +9 代fish 2026-04-09 9/450 2026-04-09 21:46 by bljnqdcc
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 化工求调剂! +34 RichLi_ 2026-04-06 34/1700 2026-04-09 20:39 by zhouxiaoyu
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[考研] 材料299专硕求调剂 +10 +21 2026-04-09 10/500 2026-04-09 17:34 by 1753564080
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +3 m永o不v言o弃m 2026-04-09 3/150 2026-04-09 17:17 by 紫曦紫棋
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 085404 293求调剂 +7 勇远库爱314 2026-04-08 7/350 2026-04-09 16:02 by 猪会飞
[考研] 一志愿上海大学生物学346 +4 上海大学346调剂 2026-04-03 4/200 2026-04-09 10:52 by yiminglu
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 320分人工智能调剂 +9 振—TZ 2026-04-03 10/500 2026-04-08 19:56 by 振—TZ
[考研] 生物医药调剂|SCI中科院三区一作+多项科研成果 +8 likangxing 2026-04-07 11/550 2026-04-08 00:02 by lys0704
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
[考研] 293求调剂 +3 勇远库爱314 2026-04-06 3/150 2026-04-07 11:15 by hugr
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 309求调剂 +4 快乐的小白鸽 2026-04-04 5/250 2026-04-04 15:55 by cql1109
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 学硕288调剂!!! +3 小王xw123 2026-04-03 3/150 2026-04-03 21:20 by 啵啵啵0119
信息提示
请填处理意见