计算化学公社

标题: 计算空间中任一点的电子密度求助 [打印本页]

作者
Author: jtravers70 时间: 2020-6-1 04:07
标题: 计算空间中任一点的电子密度求助
各位大侠，

我想计算空间中任一点的电子密度。我用了Sobereva 老师的程序（http://sobereva.com/182），但发现一个奇怪的现象。我计算了一基态小分子（甲烷），无对称性。我发现在电子密度小的地方， Sobereva 老师的程序给出的结果和高斯cubegen 产生的结果非常一致，但在密度大的地方，Sobereva 老师的程序给出的结果远大于高斯cubegen 产生的结果（>5倍）。我用的是高斯16A03 产生的wfu 文件。各位大侠，能给些建议吗？谢谢！

作者
Author: sobereva 时间: 2020-6-1 07:45
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反映出帖子具体内容，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“电子密度空间分布”改了，以后务必注意

不要用那个程序，本来就没有任何实用性，纯粹是教学目的。直接用Multiwfn算。进入主功能1，输入坐标就能算出来密度
不懂Multiwfn的话看
Multiwfn FAQ
http://sobereva.com/452
Multiwfn入门tips
http://sobereva.com/167

不存在叫做wfu的文件，要么wfn要么wfx，仔细看
详谈Multiwfn支持的输入文件类型、产生方法以及相互转换
http://sobereva.com/379（http://bbs.keinsci.com/thread-6020-1-1.html）

对于wfn文件、没使用赝势的情况，Multiwfn算出来的密度肯定和cubegen一样，不一样一定是操作错了，诸如坐标输错了等等，没细节没法说。根本没必要用cubegen，因为Multiwfn算电子密度速度远快于cubegen，功能也远比cubegen强大得多、使用上灵活得多。

作者
Author: jtravers70 时间: 2020-6-1 20:53

sobereva 发表于 2020-6-1 07:45
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反 ...

谢谢sobereva老师！我用Multiwfn！同时能不能请教一下原来的哪个程序哪里有问题？

作者
Author: jtravers70 时间: 2020-6-2 04:06

sobereva 发表于 2020-6-1 07:45
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反 ...

Sobereva老师，我用了Multiwfu 3.7 计算电子密度，我发现了同样现象。在电子密度小的的地方，Multiwfu 和cubegen的结果非常一致。电子密度大的的地方，Multiwfu 的结果远大于cubegen的结果. 我使用了高斯fchk，在Multiwfu中选了1，然后输入了坐标（bohr)。麻烦Sobereva老师给一些建议。谢谢！

作者
Author: sobereva 时间: 2020-6-2 08:53

jtravers70 发表于 2020-6-2 04:06
Sobereva老师，我用了Multiwfu 3.7 计算电子密度，我发现了同样现象。在电子密度小的的地方，Multiwfu ...

世界上没有叫Multiwfu的程序
n和u辨别不清楚的话说明该重新配眼镜了

作者
Author: jtravers70 时间: 2020-6-2 09:37
不好意思！写错了

作者
Author: sobereva 时间: 2020-6-3 21:18
我不知道你是具体怎么用cubegen计算某个点的坐标的，有fchk文件，有具体的操作步骤才能具体回复

作者
Author: jtravers70 时间: 2020-6-3 22:26
谢谢Sob老师回复！我已经找到了问题的根源。高斯cubegen必须用fdensity,否则它不考虑内核的电子！用fdensity后cubegen和multiwfn一致！

作者
Author: xxzj 时间: 2021-4-16 11:11

sobereva 发表于 2020-6-1 07:45
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反 ...

老师，您好，当卟啉中心一个金属和N2形成配位键后，金属会将一部分d电子转移到N2的Π*轨道上，如何考察这个转移的量呢？有什么方法可以计算得到吗？

作者
Author: sobereva 时间: 2021-4-21 15:38

xxzj 发表于 2021-4-16 11:11
老师，您好，当卟啉中心一个金属和N2形成配位键后，金属会将一部分d电子转移到N2的Π*轨道上，如何考察这 ...

此问题和此帖主题无关，应当单独发帖，或者发到http://bbs.keinsci.com/thread-10181-1-1.html的回复里

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