计算化学公社

标题: Multiwfn对cp2k电子密度进行盆分析的AIM电荷怎么获得 [打印本页]

作者
Author: Jerry_Fan 时间: 2025-4-2 11:39
标题: Multiwfn对cp2k电子密度进行盆分析的AIM电荷怎么获得
本帖最后由 Jerry_Fan 于 2025-4-2 15:14 编辑

各位老师好，书接上回，一些cp2k计算设置和体系介绍在这一篇帖子（CP2K计算的电荷问题请教 - 第一性原理 (First Principle) - 计算化学公社）有详细的说明，inp文件和之前一样（文件2）。计算的体系简而言之就是硅酸三钙和H2O的界面体系，现在就是在最后使用multiwfn进行盆分析的时候，有些不太明白结构中的结果。
在这篇求助：求助：CP2K如何计算bader电荷 - 第一性原理 (First Principle) - 计算化学公社中，社长贴出的图还是很明显的能够获得Na和Cl离子的电荷，但是在我的图中（图1），原子周围的盆积分似乎都不太一样，请问，我应该怎么获得每个原子的bader电荷？
比如我想要获得Ca1的电荷（图1），那么它周围有240、246、254、292几个盆，在bader.txt（文件1）中他们的值都不一样（0.7188862992、2.8252704528、0.2459707539、1.3708382772），而且我Ca使用的基组是DZVP-MOLOPT-SR-GTH-q10，也就是说用哪个都得不到想要的1.5左右的bader电荷值。这时候应该怎么办？
因为这是想要和文献对比验证一下可行性，文献中只给出了bader电荷分析。

.molden文件太大，我把它和部分输入输出文件放在百度网盘，链接如下：https://pan.baidu.com/s/1oy_ZzNeGNI0DGeFidyNJgw 提取码: w9n2

作者
Author: sobereva 时间: 2025-4-2 14:29
明明是Multiwfn基于CP2K的波函数/电子密度做的盆分析，什么时候成了“CP2K盆分析”了？整个帖子居然连最关键的做盆分析的Multiwfn的名字都不提

自觉把当前是什么体系交代清楚、上传CP2K的输入文件

作者
Author: Jerry_Fan 时间: 2025-4-2 15:15

sobereva 发表于 2025-4-2 14:29
明明是Multiwfn基于CP2K的波函数/电子密度做的盆分析，什么时候成了“CP2K盆分析”了？整个帖子居然连最关 ...

好的社长，我重新编辑了帖子并且上传了文件如上

作者
Author: sobereva 时间: 2025-4-2 17:06
输入文件里不要用Ow、Hw这种原子名，否则Multiwfn根本无法识别那是什么原子。一律用元素名
并且强烈建议用CP2K直接产生的电子密度cube文件当盆分析的输入文件
另外，当前盒子不是严格正交的。没有理由不用正交盒子。目前Multiwfn的盆分析不支持非正交盒子（以后会支持）

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