计算化学公社

标题: 求助Ru2+离子利用Multiwfn单独计算4d，4p，4s轨道的电子径向分布函数 [打印本页]

作者
Author: chenqing 时间: 2026-1-10 17:03
标题: 求助Ru2+离子利用Multiwfn单独计算4d，4p，4s轨道的电子径向分布函数
本帖最后由 chenqing 于 2026-1-10 17:01 编辑

求助Ru2+离子利用Multiwfn单独计算4d，4p，4s轨道的电子径向分布函数出现问题，感觉是设置轨道占据数目时出现问题

Ru2+离子的计算采用Gaussian的泛函和基组是b3lyp\def2tzvp
之后采用Multiwfn计算不同轨道的电子径向分布函数，流程如下
导入计算之后的fch文件6 //修改波函数26 //修改轨道占据数0 //选择所有轨道0 //占据数清零1 //选择1号轨道（对应4s）1 //占据数设1.0q //返回-1 //返回主界面200 //其它功能（part 2）5 //绘制实空间函数的径向分布函数图

0//计算径向分布函数及其积分曲线

1//绘制径向分布函数（RDF）出现图片如下图所示：4s轨道
(, 下载次数 Times of downloads: 0)
同样的方式计算，4p轨道，4d轨道，结果如下：
(, 下载次数 Times of downloads: 0) ， (, 下载次数 Times of downloads: 0)
我是为了复现这篇文献https://doi.org/10.1126/science.aav9104的工作，如下图所示，其中A图显示的是二价铁离子的3d，3p，3s的RDF，我利用上述的过程已经成功复现，但是复现下面的B图二价钌离子的4s，4p，4d轨道的RDF显示却是出现了问题，Multiwfn计算结果如上面三个图所示，一直想不明白原因，所以想请各位大佬帮忙，解决下这个问题.
(, 下载次数 Times of downloads: 0)

作者
Author: KazusaT 时间: 2026-1-10 22:39
这样操作得到的应该是ρ吧，文献说的径向分布应该是4π*r^2*ρ，你可以参考sob老师的这篇博文
使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响

作者
Author: sobereva 时间: 2026-1-10 23:02
先把横坐标和文献设得一致再比
并且弄清楚文献用的是什么基组。def2-TZVP对Ru是赝势基组，绘制这种图应当用全电子基组

作者
Author: chenqing 时间: 2026-1-13 10:12

KazusaT 发表于 2026-1-10 22:39
这样操作得到的应该是ρ吧，文献说的径向分布应该是4π*r^2*ρ，你可以参考sob老师的这篇博文
使用Multiwf ...

这篇文献介绍得很详细，受益匪浅，感谢帮助。

作者
Author: chenqing 时间: 2026-1-13 10:14

sobereva 发表于 2026-1-10 23:02
先把横坐标和文献设得一致再比
并且弄清楚文献用的是什么基组。def2-TZVP对Ru是赝势基组，绘制这种图应当 ...

感谢解惑，这个确实是基组的问题，经过选择全电子基组，问题解决了，感谢！

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