如何得到spdf轨道上原子电子价密度

最瘦二师兄

请教各位大神，有没有什么方法可以从高斯的输出文件中得到“for each spin the atomic electron valence density (integrated) per L-value”，就像ADF那样，分别输出是s, p, d, f轨道上的原子电子价密度，如下图所示。

sobereva

含有基函数信息的波函数文件（对于Gaussian是fch/fchk文件）载入Multiwfn（http://sobereva.com/multiwfn）后输入
7
5
1

然后会看到Mulliken布居分析结果，包括各个角动量轨道上的电子数


如果只需要得到价电子的，载入文件后输入
6
34  //把内层轨道占据数设0
-1
7
5
1

最瘦二师兄

sobereva 发表于 2021-7-20 21:22
含有基函数信息的波函数文件（对于Gaussian是fch/fchk文件）载入Multiwfn（http://sobereva.com/multiwfn） ...

谢谢sob老师，按照您教的第二种方法（得到价电子的），我得到了不同轨道轨道上的电子数，但所有alpha轨道上的布局都是0 (图population_1)，全部的净电荷为301 (图ner charge_1)，这显然很不合理。我猜这是由于我把内层轨道电子数设置为0了（我这是开壳层体系）。所以我又尝试了第一种方法，这次显示是合理的（图population_2），而且正6价的W(2)原子p轨道上的电子总数只有5.25，接近6，我可把它当成W6+最外层p轨道的电子数吗？但是如果看正4价的Si(1)话的s轨道上的电子数的话，把它当成最外层又不合理。所以我想问该如何计算开壳层价电子在不同角动量轨道上的布局？

sobereva

最瘦二师兄 发表于 2021-7-21 15:12
谢谢sob老师，按照您教的第二种方法（得到价电子的），我得到了不同轨道轨道上的电子数，但所有alpha轨道 ...

把输入文件压缩后上传我看看

W用了赝势，直接给出的布居数就是对应于价电子的情况。Si没用赝势，显然布居数还包含了内层电子

最瘦二师兄

sobereva 发表于 2021-7-22 01:01
把输入文件压缩后上传我看看

W用了赝势，直接给出的布居数就是对应于价电子的情况。Si没用赝势，显然 ...

我试着上传过，但是显示上传失败，还有其他方法可以分享给您吗？

sobereva

最瘦二师兄 发表于 2021-7-23 13:10
我试着上传过，但是显示上传失败，还有其他方法可以分享给您吗？

传网盘，贴链接