求助Ru2+离子利用Multiwfn单独计算4d，4p，4s轨道的电子径向分布函数

chenqing

求助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轨道

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

KazusaT

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

sobereva

先把横坐标和文献设得一致再比
并且弄清楚文献用的是什么基组。def2-TZVP对Ru是赝势基组，绘制这种图应当用全电子基组

chenqing

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

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

chenqing

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

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