关于cp2k输出的mulliken文件的分析问题

Flights · 发表于 Post on 2020-4-5 10:30:37

下图是我把用cp2k计算MD中的一个片段截取出来做mulliken分析的结果输出，就想请教一下各位大佬，该怎么分析这个输出数据？

sobereva · 发表于 Post on 2020-4-5 15:04:42

仔细看《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）了解原子电荷的基本常识
并且阅读此文了解什么叫自旋布居
谈谈自旋密度、自旋布居以及在Multiwfn中的绘制和计算
http://sobereva.com/353（http://bbs.keinsci.com/thread-4422-1-1.html）

jiejie · 发表于 Post on 2022-7-28 16:36:06

社长大大，《原子电荷计算方法的对比》文章中mulliken电荷部分的原子电荷就是CP2K导出的Net charge对吧。如果net charge是负的代表该原子得到电子对吧？

jiejie · 发表于 Post on 2022-7-28 16:38:24

我按照上面提出的理解，发现mulliken结果和我其他计算结果是完全相反的

jiejie · 发表于 Post on 2022-7-28 17:55:41

sobereva 发表于 2020-4-5 15:04
仔细看《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）了解原 ...

求助社长大大，问题在上面

Eudaimonia · 发表于 Post on 2022-7-28 18:13:25

jiejie 发表于 2022-7-28 16:38
我按照上面提出的理解，发现mulliken结果和我其他计算结果是完全相反的

Mulliken 电荷存在一些严重问题: (1) 交叉项平分的做法物理意义不严格, 没有考虑到原子间的差异性; (2) 基组依赖性非常大. 尤其是使用大基组、含弥散函数基组的情况下问题更为明显

社长文章里有说

jiejie · 发表于 Post on 2022-7-28 18:22:37

Eudaimonia 发表于 2022-7-28 18:13
Mulliken 电荷存在一些严重问题: (1) 交叉项平分的做法物理意义不严格, 没有考虑到原子间的差异性; (2) ...

嗯嗯，谢谢回复，我觉得应该有点问题。那如果直接将态密度积分到费米面，这种方法比较可靠吗？另外，CP2K的input文件中可以写入&PRINT &E_DENSITY_CUBE 然后得到电子密度图，那个图靠谱吗？我通过该电子图看出石墨烯衬底上电子密集的地方，mulliken net charge反而是正的。

sobereva · 发表于 Post on 2022-7-28 19:15:18

jiejie 发表于 2022-7-28 17:55
求助社长大大，问题在上面

提问时候一次性说完整，别分好几个回帖零零散散的描述

sobereva · 发表于 Post on 2022-7-28 19:16:14

jiejie 发表于 2022-7-28 18:22
嗯嗯，谢谢回复，我觉得应该有点问题。那如果直接将态密度积分到费米面，这种方法比较可靠吗？另外，CP2K ...

用Multiwfn基于CP2K的波函数算CM5电荷，远比Mulliken电荷强得多

按下文产生molden文件
使用Multiwfn结合CP2K通过NCI和IGM方法图形化考察固体和表面的弱相互作用
http://sobereva.com/588（http://bbs.keinsci.com/thread-21742-1-1.html）

然后用Multiwfn主功能7的相应子功能计算CM5电荷

jiejie · 发表于 Post on 2022-7-28 21:56:34

sobereva 发表于 2022-7-28 19:16
用Multiwfn基于CP2K的波函数算CM5电荷，远比Mulliken电荷强得多

按下文产生molden文件

嗯嗯！感谢社长！

jiejie · 发表于 Post on 2022-7-30 11:46:38

sobereva 发表于 2022-7-28 19:16
用Multiwfn基于CP2K的波函数算CM5电荷，远比Mulliken电荷强得多

按下文产生molden文件

社长大大，我用7，16计算了cm5电荷得到了chg文件，应该如何做进一步做占据数分析呢，文件也看不懂。第一次用multiwfn，找了些帖子看但没搞懂

sobereva · 发表于 Post on 2022-7-31 03:45:13

jiejie 发表于 2022-7-30 11:46
社长大大，我用7，16计算了cm5电荷得到了chg文件，应该如何做进一步做占据数分析呢，文件也看不懂。第一 ...

原子电荷 = 原子核电荷数 - 原子所带电子数
如果你想得到各个原子轨道的占据数，光是计算诸如CM5这样的原子电荷是没用的，原子电荷是对整个原子而言的。
得到各个原子轨道的占据数可以用Multiwfn做Mulliken分析，参考
利用布居分析判断基函数与原子轨道的对应关系
http://sobereva.com/418（http://bbs.keinsci.com/thread-10094-1-1.html）

lihao2690 · 发表于 Post on 2023-3-20 22:24:05

您好，请问输入什么关键词能够给出mulliken电荷呢

sobereva · 发表于 Post on 2023-3-24 05:27:32

lihao2690 发表于 2023-3-20 22:24
您好，请问输入什么关键词能够给出mulliken电荷呢

Multiwfn创建CP2K输入文件界面里有个选项是计算原子电荷，里面选上Mulliken。另外，PRINT_LEVEL MEDIUM的时候也默认会输出。顺带一提，用Multiwfn载入CP2K的molden文件做Mulliken布居分析，能得到远比Mulliken原子电荷更丰富的信息。

lihao2690 · 发表于 Post on 2023-3-30 10:29:04

sobereva 发表于 2022-7-28 19:16
用Multiwfn基于CP2K的波函数算CM5电荷，远比Mulliken电荷强得多

按下文产生molden文件

谢谢sob老师，发现跑MD的话产生molden文件会使计算速率降低一倍左右，另外每一步会产生一个molden文件，这些如何在multiwfn中进行批处理呢？谢谢

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 关于cp2k输出的mulliken文件的分析问题