计算化学公社

标题: 求助使用Multiwfn计算亲电描述符一直为NaN [打印本页]

作者
Author: zohn_ecnu 时间: 2025-10-20 14:45
标题: 求助使用Multiwfn计算亲电描述符一直为NaN
本帖最后由 zohn_ecnu 于 2025-10-31 22:22 编辑

各位老师好，
我想请教关于使用Multiwfn计算亲电描述符为NaN的问题

1、我使用windows的Gaussian09，在B3LYP/6-31g(d,p) scrf=(smd,solvent=water) em=gd3bj这些关键词下自行计算苯酚的单点能，分别计算N+1，N，N-1，N-2四个电子状态的单点能，产生.chk文件，然后经Gaussian09转化成.fch格式
2、其次启动Multiwfn进行各种量的计算
3、程序开始计算，但生成的文件中亲电描述符为NaN： Electrophilic descriptor (epsilon): NaN Hartree, NaN eV

作者
Author: sobereva 时间: 2025-10-20 23:17
根本没算对
N+1是体系多一个电子的状态，当然应该是-1电荷的阴离子，你却当成了1 2来算
认真把下文末尾看了
使用Multiwfn超级方便地计算出概念密度泛函理论中定义的各种量
http://sobereva.com/484（http://bbs.keinsci.com/thread-13328-1-1.html）

作者
Author: zohn_ecnu 时间: 2025-10-21 17:16

sobereva 发表于 2025-10-20 23:17
根本没算对
N+1是体系多一个电子的状态，当然应该是-1电荷的阴离子，你却当成了1 2来算
认真把下文末尾看 ...

好的，非常感谢社长！

作者
Author: zohn_ecnu 时间: 2025-10-31 22:21

sobereva 发表于 2025-10-20 23:17
根本没算对
N+1是体系多一个电子的状态，当然应该是-1电荷的阴离子，你却当成了1 2来算
认真把下文末尾看 ...

社长，我又算了一下亚硫酸根阴离子的CDFT定义的各种量，电荷和自选多重度已按照教程（http://sobereva.com/484）设置好，但这次又出现了 Electrophilic descriptor (epsilon): NaN Hartree, NaN eV，请问这又是什么原因呢？

作者
Author: sobereva 时间: 2025-11-1 05:07

zohn_ecnu 发表于 2025-10-31 22:21
社长，我又算了一下亚硫酸根阴离子的CDFT定义的各种量，电荷和自选多重度已按照教程（http://sobereva.co ...

这是Electrophilic descriptor (epsilon)指数定义的原因，可以理解为其局限性。如Multiwfn手册所示，这个指数有个根号项，当里面的值为负值的时候，就没法计算了。
(, 下载次数 Times of downloads: 0)

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