24小时热门版块排行榜    

查看: 829  |  回复: 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射手座 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 15:18 by KM5EcsNQRBPn
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:45 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 14:30 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:21 by KM5EcsNQRBPn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见