求助使用Multiwfn和GaussView分析静电力，空间电荷，空间力

sakura1109

背景更新：在寻找关于C-H键氧化提取过渡态的过程中，发现在小分子同一侧的2个H原子（7、14）能顺利找到过渡态，另一侧（6、15）死活找不到，四个H原子最初认为应该是完全等效的，也不是关于R、S构型的选择性问题，想要解释这个问题，考察了Hirshfeld电荷和Mulliken电荷发现确实有点细微差距，Hirshfeld电荷0.024左右，Mulliken电荷约0.09（7、14）；另一侧Hirshfeld电荷同侧0.038左右，Mulliken电荷约0.11（6、15），（其实有点怀疑这个问题真的有探讨和解答的必要吗，还是只是单纯没有找到？）
如果不适用静电力，空间电荷，空间力来解决，进一步应该考察什么性质呢？
原帖：在Phys. Chem. Chem. Phys., 2017,19, 1496-1503 这篇文献中看到使用了Multiwfn和 GaussView分析静电力，空间电荷，空间力，想问问社长可以再详细展开讲讲如何做出这样的图吗？（在Multiwfn手册中有搜索electrostatic force、steric charge、steric force确实没有找到）想看看能否用这个方法解决疑问


或者有相关的文献、帖子、方法推荐吗？非常感谢！

sobereva

把你的体系结构图贴出来，描述清楚具体反应，别人才可能针对性地回答你

你说的那篇文章里GaussView只不过是基于Multiwfn计算的格点数据作图，GaussView自身没有分析功能。这个分析耗时极高，对于你的问题也未必有什么帮助。

sakura1109

sobereva 发表于 2022-5-10 14:05
把你的体系结构图贴出来，描述清楚具体反应，别人才可能针对性地回答你

你说的那篇文章里GaussView只不 ...

好的社长，已编辑

sobereva

找不到过渡态的问题并不是另一面的H的电子结构相对于7、14有什么差异，应当只是纯粹因为找过渡态的时候上面酰胺部分碍事，和另一个反应物分子有相互作用所致。

当前问题没必要纠结另一侧怎么找不到（虽然恰当用一些优化技巧应该能找到），7、14能找到过渡态就足够了，就已经证明能反应了。本身这个体系两侧的H是化学等价的，而且上面酰胺部分又容易翻转。