这个太难了，还得请各位功底深厚的大伽写个代码才是！

如果硬把轨道做成LCAO的，AO的中心和分子质心不是同一个中心，两个球谐函数可能会牵扯到椭圆坐标，应该会比较麻烦。
这种情况反而是离散成格点用Lebedev的方法做会方便一些，不过均匀格点转换成球形格点应该需要插值。
原文应该是用平面波做的，如果赝波函数比较光滑的话，利用Rayleigh展开把平面波展开成球面波，利用球谐函数和球形Bessel函数的正交性质，应该可以解析或半数值的方式来处理这些积分。

zhangfuqiang 发表于 2020-4-12 13:35
想得到这张图。轨道一个个找，好办。但不知如何得到supterAtom的S,P，D等等轨道的PDOS。公式在SI中Referenc ...

这没法直接由Multiwfn得到，对于整个体系而言的S,P,D等等，得自己写个代码，将分子轨道向不同类型的球谐函数投影来得到

量子化学程序能算什么体系Multiwfn原理上就能分析什么体系
必须说清楚“做”到底指做什么分析