计算化学公社

标题: 如何使用Multiwfn用原子电荷和波函数计算同一分子在相同表面上静电势 [打印本页]

作者
Author: ZDcccA 时间: 2026-1-29 15:54
标题: 如何使用Multiwfn用原子电荷和波函数计算同一分子在相同表面上静电势
各位老师：
  我希望比较RESP电荷和波函数电子密度方法计算同一个短肽分子表面静电势的差别，在两种方法中需要使用同一个表面，但不太清楚如何在Multiwfn（Version 3.8(dev)2025-Oct-2）中设定同一表面。
  我尝试都使用Promolecular electron density划定面但发现若输入文件不是.chg而是波函数文件，Promolecular electron density的选项似乎没有。
  还请各位老师赐教！

作者
Author: sobereva 时间: 2026-1-29 17:34
仔细看Multiwfn手册4.7.8节，专门说了怎么考察原子电荷对静电势的重现性

作者
Author: ZDcccA 时间: 2026-1-30 11:51

sobereva 发表于 2026-1-29 17:34
仔细看Multiwfn手册4.7.8节，专门说了怎么考察原子电荷对静电势的重现性

谢谢社长！仔细看了手册4.7.8节，获益匪浅。
另外想礼貌追问一下，我的目标是在PB模型框架下探究原子电荷的有效性，溶质-溶剂界面选择的是solvent-excluded surface (SES)，是不是比较原子电荷和波函数电子密度计算的在SES表面上的静电势会更合理一些？但是使用主功能7中的功能13似乎只能是范德华表面。是否有办法指定SES呢，还是需要使用外部程序算好SES表面再导入？

作者
Author: sobereva 时间: 2026-1-31 04:23

ZDcccA 发表于 2026-1-30 11:51
谢谢社长！仔细看了手册4.7.8节，获益匪浅。
另外想礼貌追问一下，我的目标是在PB模型框架下探究原子电 ...

不一定更合理。如果你想尝试的话，可以让Multiwfn基于原子电荷或基于波函数计算外部程序产生的SES表面顶点上的静电势。Multiwfn主功能5可以载入文本文件里记录的一批外部坐标，一次性对它们计算，参考Multiwfn手册3.6节中Special case: Calculate data for a set of arbitrarily distributed points部分。

作者
Author: ZDcccA 时间: 2026-2-1 18:45

sobereva 发表于 2026-1-31 04:23
不一定更合理。如果你想尝试的话，可以让Multiwfn基于原子电荷或基于波函数计算外部程序产生的SES表面顶 ...

明白了，多谢社长的耐心指导！

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