计算化学公社

标题: 求助slab模型的偶极矩计算 [打印本页]

作者
Author: diffcwmol 时间: 2020-5-11 20:17
标题: 求助slab模型的偶极矩计算
请问各位同仁，我在使用VASP计算一个电中性的表面时，遇到如下问题：
使用tag，如LDIPOL=.TRUE.，IDIPOL=3，从OUTCAR里面读取dipolmoment数值。

但是计算结果和我用高斯计算得到的值符号相反的，虽然functional和basis是不同的，但我想偶极矩的正负性至少应该相同吧。

所以我觉得OUTCAR里面输出的dipolmoment应该不是表面的偶极矩（可能是偶极矫正加的补偿偶极矩？），请问有没有同仁碰到过类似的问题？

作者
Author: Xian 时间: 2022-7-6 19:11
我也想slab模型的偶极矩，请问楼主的问题有答案了吗？谢谢。

作者
Author: diffcwmol 时间: 2022-7-14 16:25

Xian 发表于 2022-7-6 19:11
我也想slab模型的偶极矩，请问楼主的问题有答案了吗？谢谢。

就是偶极矫正加的补偿偶极矩

作者
Author: Xian 时间: 2022-7-21 16:44

diffcwmol 发表于 2022-7-14 16:25
就是偶极矫正加的补偿偶极矩

好的呢，谢谢啦

作者
Author: Jack 时间: 2022-11-22 19:17

diffcwmol 发表于 2022-7-14 16:25
就是偶极矫正加的补偿偶极矩

这表示实际偶极矩与OUTCAR中的补偿偶极矩互为相反数么？

作者
Author: pupu1616 时间: 2023-8-7 09:57

Jack 发表于 2022-11-22 19:17
这表示实际偶极矩与OUTCAR中的补偿偶极矩互为相反数么？

请问这个问题解决了吗

作者
Author: pupu1616 时间: 2023-8-7 10:00

diffcwmol 发表于 2022-7-14 16:25
就是偶极矫正加的补偿偶极矩

请问就是OUTCAT里
DIPCOR: dipole corrections for dipol
direction  3 min pos 326,
dipolmoment          0.000000    0.000000    -0.204589 electrons x Angstroem
Tr[quadrupol]    -1182.295529
其中-0.204589就是偶极矩吗。还是要相反数

作者
Author: Jack 时间: 2023-8-7 14:36

pupu1616 发表于 2023-8-7 10:00
请问就是OUTCAT里
DIPCOR: dipole corrections for dipol
direction 3 min pos 326,

是的

作者
Author: pupu1616 时间: 2023-8-7 16:35

Jack 发表于 2023-8-7 14:36
是的

前辈。那这意思是它偶极矩大小是0.204589，负号意味着方向是Z轴负方向，那正电中心在靠下面，负电中心在靠上面，粗略是这样分析吗

欢迎光临计算化学公社 (http://bbs.keinsci.com/)