分子のエネルギー計算

多中心積分ができるようになったので、後は結構楽でした。分子の計算もできます。

↑入力ファイルはこんな感じで、ＨとＯを配置。

水分子の電子密度。左はエネルギー順位で、↑↓はスピンの向きを示す。

↑これはナトリウム原子の電子密度。電子軌道も書いてます。
ハチの字は2p軌道です。

ナトリウムのエネルギー順位。ちなみに、このときの軌道エネルギーの総和は-161.84(hartree)で
本に載っている値-161.85とほとんど一緒でした。Gaussian Lobesでも結構いけることがわかります。