ガウス軌道の概要
計算化学や
分子物理学の分野において、ガウス軌道(Gaussian orbital)またはガウス型軌道(Gaussian type orbital, GTO)は、
原子軌道や
LCAO法に基づく
分子軌道を計算するための重要な関数です。この手法は、電子の位置やエネルギーを正確に予測するために広く利用されています。
ガウス軌道の導入
1950年にフランシス・ボーイズが初めてこの概念を導入しました。従来の
スレーター軌道に代わり、ガウス型軌道を用いることにより、分子の
電子状態の計算が効率化されることが保証されました。ガウス型関数の積法則によって、異なる原子を中心とするGTOを、有限個のガウス関数の和として表現できるため、計算が簡潔化され、特定の条件にもとづいて2中心積分や1中心積分に変換できることが可能です。この結果、計算時間が4〜5桁も短縮される場合があり、ガウス軌道の利用が広がる理由となっています。
数学的構成
ガウス基底関数の形式は、一般的な動径部分と角度部分に分かれています。具体的には、
$$Φ(r) = R_l(r)Y_{lm}(θ, ϕ)$$
という形で、球面調和関数に基づく記述がなされています。ここで、$Y_{lm}(θ, ϕ)$は球面調和関数を示し、動径部分$R_l(r)$はそれぞれ以下のようになります。
$$R_l(r) = A(l, α)r^le^{-αr}$$
($A(l, α)$は規格化定数)
$$R_l(r) = B(l, α)r^le^{-αr^2}$$
($B(l, α)$はガウス軌道用の規格化定数)
原始ガウス軌道は単独では
原子軌道を十分に表現できないため、複数の原始ガウス軌道を重ね合わせることで
スレーター軌道に近づけます。
$$R_l(r) = r^l∑_{p=1}^P c_p A(l, α_p)e^{-α_pr^2}$$
ここで、$c_p$は各原始ガウス軌道に対する縮約係数です。これにより、計算の効率がさらに向上します。
ガウス軌道の変種
原子は球対称性を持つため、ガウス軌道の角度部分は球面調和関数で表現されますが、計算が複雑になることがあります。そこで、デカルト座標系を利用したデカルトガウス軌道がしばしば使われます。
エルミート多項式の利用により計算が簡便になり、エルミートガウス軌道として計算も行われることがあります。
分子積分とその進展
ガウス軌道に関連する分子積分の計算方法は1966年に竹田らによって提示され、その後も多くの研究が行われています。特に1978年には、McMurchieとDavidsonによって漸化関係式が導入され、計算時間が大幅に短縮されました。その後も、PopleとHehreが局所座標法を開発し、計算精度を向上させています。そして、今も様々な研究が続けられ、効率的な計算手法が模索されています。
POLYATOMの革新
POLYATOMは、このようなガウス軌道を活用し、第一原理
分子軌道計算を行うための画期的な計算パッケージです。MITのSolid State and Molecular Theory Groupによって開発され、ガウス軌道の特性を生かした計算が可能となりました。
ガウス軌道は
計算化学において基本的かつ重要な要素であり、その発展と応用は今後も注目されるでしょう。