24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 733  |  回复: 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射手座 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
[考研] 265求调剂 +11 小木虫085600 2026-04-06 12/600 2026-04-07 08:37 by 蓝云思雨
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +3 努力奋斗112 2026-04-06 3/150 2026-04-07 08:35 by jp9609
[考研] 285求调剂 +8 AZMK 2026-04-05 13/650 2026-04-07 08:14 by xihu1109
[考研] 08600生物与医药-327 +9 18755400796 2026-04-05 9/450 2026-04-06 22:35 by 52305043001
[考研] 372分,材料与化工,一志愿湖南大学,求调剂 +3 蓝笺片 2026-04-01 3/150 2026-04-06 09:04 by 无际的草原
[考研] 一志愿哈尔滨工业大学085600英一数二337分求调剂 +12 lyz0427 2026-04-03 12/600 2026-04-06 06:37 by houyaoxu
[考研] 考研生物学考A区211,初试322,科目生化和生物综合,求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] +5 雾与海 2026-04-02 6/300 2026-04-04 19:53 by 蓝云思雨
[考研] 085602 找调剂 +4 逆时针快乐 2026-04-02 4/200 2026-04-04 19:32 by 蓝云思雨
[考研] 292求调剂 +11 2022080213 2026-04-04 13/650 2026-04-04 18:38 by macy2011
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 266求调剂 +3 08电气工程 2026-04-03 3/150 2026-04-03 14:05 by 1753564080
[基金申请] 请问共同通讯和共同一作的认可度问题 10+4 psa1234 2026-04-01 10/500 2026-04-03 11:08 by Kittylucky
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 293求调剂 +4 珂珂乐 2026-04-02 4/200 2026-04-02 20:10 by 6781022
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-01 11/550 2026-04-02 10:52 by lnilvy
[考研] 348环境工程调剂 +3 吴彦祖24k 2026-04-01 3/150 2026-04-02 09:14 by nanaliuyun
信息提示
请填处理意见