非球対称ポテンシャルでは、波動関数を球面調和函数で展開すると角運動量を添字とする行列になることを示した。
koideforest.hatenadiary.com
角度積分するのに掛かる時間を、必要な最大の軌道角運動量と動径メッシュの数をそれぞれ、とすると、全ての行列要素を求める時間はと表せる。
は、例えば、まで必要とすると、で4つ必要であり、それを表せている。
次に、直接求めるのではなく、ポテンシャルを角運動量(球面調和関数)で展開することを考える。
この展開に掛かる時間はと表せる。
これを使って、波動関数を角運動量展開した時の行列要素は、以下のように表せる。
はGaunt係数である。は、最大での2倍が必要であるため、Gaunt係数を求める時間をとすると、行列要素全てを求める時間はとなる。
ポテンシャルを展開するところと合わせれば、と表せる。
したがって、
の時、Gaunt積分の方法の方が速くなる。
ただ、直感的には、どっちもあまり変わらなそうな気もする。