24小时热门版块排行榜    

查看: 728  |  回复: 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射手座 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 268求调剂 +7 好运连绵不绝 2026-03-12 8/400 2026-03-17 20:28 by xilongliang
[考研] 070300化学319求调剂 +4 锦鲤0909 2026-03-17 4/200 2026-03-17 18:21 by 重科小霸王
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 271求调剂 +12 生如夏花… 2026-03-11 14/700 2026-03-17 10:56 by lovewei0727
[考研] 267一志愿南京工业大学0817化工求调剂 +6 SUICHILD 2026-03-12 6/300 2026-03-17 09:24 by 雾散后相遇lc
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
信息提示
请填处理意见