计算化学公社

标题: 能否使用Multiwfn输出空间中一大批点的Total ESP [打印本页]

作者
Author: Kalinite 时间: 2026-1-1 20:24
标题: 能否使用Multiwfn输出空间中一大批点的Total ESP
Multiwfn的主功能1可以输出空间中任意一点的各种性质，其中也包含该点的ESP。
我想通过Multiwfn计算一大批点（通常数万个）的Total ESP，如果反复使用主功能1，因为主功能1还计算和输出了许多其他的信息，实际操作起来数万个点需要计算几十分钟。有没有办法调用Multiwfn只计算ESP的值，从而加速完成计算大量点的ESP？

（ORCA在带点电荷计算时是可以输出这部分信息的，但我不想使用ORCA）

作者
Author: sobereva 时间: 2026-1-1 22:27
Multiwfn主功能5可以一次性计算一批自定义位置的特定属性，手册里有介绍：
(, 下载次数 Times of downloads: 29)

作者
Author: Kalinite 时间: 2026-1-4 14:10

sobereva 发表于 2026-1-1 22:27
Multiwfn主功能5可以一次性计算一批自定义位置的特定属性，手册里有介绍：

解决了，感谢卢老师，几分钟就能算完了！

作者
Author: Kalinite 时间: 7 day ago

sobereva 发表于 2026-1-1 22:27
Multiwfn主功能5可以一次性计算一批自定义位置的特定属性，手册里有介绍：

在做QM/MM的时候，静电嵌入方法把MM原子作为背景电荷处理，计算梯度的时候需要背景电荷的梯度。ORCA在写engrad关键词的时候，可以直接输出.pcgrad，其中包含了点电荷的梯度。Gaussian并不能这么做，但理论上可以用Multiwfn做差分，计算6N个点处的ESP，差分得到每个点受到QM区的力、再加上点电荷之间的相互作用，从而从波函数出发就能尽量重现.pcgrad的结果。
实际做的时候，因为点电荷数量往往有上万个，6N达到了近10万个点，计算需要花费约半小时，已经超过了QM区计算的耗时。请问卢老师，是否还有进一步缩短时间的做法呢。

作者
Author: sobereva 时间: 7 day ago

Kalinite 发表于 2026-4-26 06:24
在做QM/MM的时候，静电嵌入方法把MM原子作为背景电荷处理，计算梯度的时候需要背景电荷的梯度。ORCA在写e ...

一种做法是对QM区计算拟合静电势电荷，靠点电荷之间相互作用公式来近似计算MM与QM区原子之间的静电作用力

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