24小时热门版块排行榜    

Znn3bq.jpeg
查看: 746  |  回复: 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 xq83 2026-04-22 3/150 2026-04-22 03:56 by 啊哒哒哒叨
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 一志愿A区211,22408 321求调剂 +7 随心所欲☆ 2026-04-15 8/400 2026-04-21 08:22 by Equinoxhua
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
信息提示
请填处理意见