Multiwfn计算Bader电荷输出异常

星殇~浮梦 · 发表于 Post on 2025-10-4 17:10:02

本帖最后由星殇~浮梦于 2025-10-4 19:20 编辑

在使用Multiwfn计算Bader电荷时，输出“ Warning: Unable to determine the attractor    2 belongs to which atom!
If this is a non-nuclear attractor, simply press ENTER button to continue. If you used pseudopotential and this attractor corresponds to the cluster of all maxima of its valence electron, then input the index of this atom (e.g. 9). Else you should input q to return and regenerate basins with smaller grid spacing”
输入分别是“17,1,1,2,7,2,1”
计算的体系是直链，额外带一个电子的烷烃，C<13。

Multiwfn输入的文件是orca计算得出的.gbw文件。命令行内容是“!  wB97X-D3BJ def2-TZVPD def2/J RIJCOSX TightSCF defgrid3 SP ”

在使用B3LYP，PBE泛函的时候，计算Bader电荷正常，同时，在wB97X下，计算不额外带电子的价态时0的体系时正常输出，在价态-1的体系中，计算C数目在1-4,6,8,10时也正常输出，没有报错。计算C数目是5,7,9,11,12时报错。

尝试了增大隔点精度，在“ Please select a method for setting up grid
-10 Set grid extension distance for mode 1~6, current: Detect rho isosurface
1 Low quality grid, spacing=0.20 Bohr, cost: 1x
2 Medium quality grid, spacing=0.10 Bohr, cost: 8x
3 High quality grid, spacing=0.06 Bohr, cost: 36x
4 Lunatic quality grid, spacing=0.04 Bohr, cost: 120x
5 Only input grid spacing, automatically set other parameters
6 Only input the number of points in X,Y,Z, automatically set other parameters
7 Input original point, translation vector and the number of points
8 Set center position, grid spacing and box length
9 Use grid setting of another cube file
10 Set box of grid data visually using a GUI window”

之中选择3,4。最后也是没有改善。仍然输出此部分内容

sobereva · 发表于 Post on 2025-10-4 19:18:00

上传压缩后的波函数文件，以及对应的ORCA输入文件

没有叫Mulitiwfn的程序

星殇~浮梦 · 发表于 Post on 2025-10-4 19:26:44

sobereva 发表于 2025-10-4 19:18
上传压缩后的波函数文件，以及对应的ORCA输入文件

没有叫Mulitiwfn的程序

不好意思sob老师，拼写错误已更改

星殇~浮梦 · 发表于 Post on 2025-10-4 19:28:42

压缩后的波函数文件和对应的orca输入文件如下

sp_w_opt_5.zip (1.28 MB, 下载次数 Times of downloads: 2)

sp_w_opt_5.inp (1.05 KB, 下载次数 Times of downloads: 4)

sobereva · 发表于 Post on 2025-10-5 05:00:15

上传molden文件不要传gbw，否则就像给别人提供Gaussian的chk文件一样，别人的版本和你的不同时根本没法转换，而且你又没说你用的是ORCA什么版本

并且你提供的输入文件的关键词和! wB97X-D3BJ def2-TZVPD def2/J RIJCOSX TightSCF defgrid3 SP不同

星殇~浮梦 · 发表于 Post on 2025-10-5 10:21:05

老师我使用的是orca6.0,这个molden文件是Multiwfn调用orca_2mkl生成的。

星殇~浮梦 · 发表于 Post on 2025-10-5 10:22:29

下面是您要求的文件

sobereva · 发表于 Post on 2025-10-5 13:02:29

你的体系的里德堡轨道的电子分布得很弥散，导致出现了额外的吸引子，对应的盆的电子密度>0.001 a.u.部分如下所示。因此没法算AIM电荷

这体现出AIM电荷在定义上的不足。还有不少其它情况也没法算AIM电荷，诸如KrH+，我写的原子电荷综述文章里专门说了
一篇深入浅出、完整全面介绍原子电荷的综述文章已发表！
http://sobereva.com/714（http://bbs.keinsci.com/thread-46067-1-1.html）

以后回复我的帖子的时候点击帖子下面的回复按钮，否则>99%的可能性我注意不到

星殇~浮梦 · 发表于 Post on 2025-10-5 15:14:05

sobereva 发表于 2025-10-5 13:02
你的体系的里德堡轨道的电子分布得很弥散，导致出现了额外的吸引子，对应的盆的电子密度>0.001 a.u.部分如 ...

好的，谢谢sob老师解答

DoubleRice · 发表于 Post on 2025-12-22 21:21:46

你好，我想问一下multiwfn生成bader电荷的步骤是什么呀，没有搜到老师的教程

sobereva · 发表于 Post on 2025-12-23 05:22:23

DoubleRice 发表于 2025-12-22 21:21
你好，我想问一下multiwfn生成bader电荷的步骤是什么呀，没有搜到老师的教程

Multiwfn手册4.17.1节明确给了对分子算AIM电荷的例子，Multiwfn quick start.pdf里明确写了：

用于周期性体系，北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里的“电子结构的分析部分”也给了例子

PS：bader电荷是对AIM电荷的谬称

DoubleRice · 发表于 Post on 2025-12-23 08:32:04

sobereva 发表于 2025-12-23 05:22
Multiwfn手册4.17.1节明确给了对分子算AIM电荷的例子，Multiwfn quick start.pdf里明确写了：

收到！非常感谢sober老师！

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

[Multiwfn使用咨询] Multiwfn计算Bader电荷输出异常

浏览过的版块