密度汎関数理論

概説

密度汎関数理論はその概念の根源をトーマス–フェルミ模型に持つものの、DFTは2つのホーエンベルク–コーンの定理（H–K）によって強固な理論的基盤の上に置かれた。最初のH–K定理は、磁場がない場合の非縮退基底状態についてのみ成り立っていたが、以後これらを包含するために一般化されてきた。
H–Kの第1定理は、多電子系の基底状態の性質が3つの空間座標だけに依存する電子密度によって一意に決定されることを論証する。これは、電子密度の汎関数に使用することによって、3つの空間座標について3N個の空間座標を持つN個の電子の多体問題を軽減するための土台を築く。この定理は、時間依存密度汎関数法（TDDFT）を開発するための時間依存定義域へ拡張することができる。TDDFTは励起状態を記述するために使うことができる。
H–Kの第2定理は、系についてのエネルギー汎関数を定義し、正しい基底状態電子密度がこのエネルギー汎関数を最小化することを示す。
コーン–シャムDFT（KS DFT）の枠組みの中では、静的外部ポテンシャル中で相互作用のある電子の扱いにくい多体問題が、有効ポテンシャル中を移動する相互作用のない電子の扱いやすい問題に軽減される。有効ポテンシャルは外部ポテンシャルと電子間のクーロン相互作用（例えば、交換相互作用や相関相互作用）の効果を含む。後者の2つの相互作用のモデル化がKS DFT内での難しさとなる。最も単純な近似は局所密度近似（LDA）であり、これは一様な電子ガスについての厳密な交換エネルギーに基づいている。このエネルギーはトーマス–フェルミ模型や、一様な電子ガスについての相関エネルギーへの当て嵌めから得ることができる。相互作用のない系は解くのが比較的簡単であり、波動関数はオービタルのスレイター行列式として表わすことができ。そのうえ、こういった系の運動エネルギー汎関数は厳密に分かる。全エネルギー汎関数の交換-相関部分は依然として不明であり、近似しなければならない。
KS DFTよりも知られていないが、ほぼ間違いなく最初のH-K定理の精神により密接に関係している別の手法が、オービタルフリー密度汎関数理論（OFDFT）である。OFDFTでは、近似汎関数が相互作用のない系の運動エネルギーについても使われる。

ホーヘンベルク・コーンの定理

電子密度を用いた物理量の計算が原理的に可能であることは1964年にヴァルター・コーンとピエール・ホーエンベルクによって示された。
ある外部ポテンシャルのもとにあるN個の電子系を考える（例えば分子の原子核の配置が決まれば、それらの原子核が電子に及ぼす静電ポテンシャルは決まる）。いま、この系の基底状態の電子密度ρだけがわかっているとする。ホーヘンベルク・コーンの第1定理によれば、ある系の基底状態の電子密度ρが決まると、それを基底状態にもつ外部ポテンシャルがもし存在すればそれはただ1通りに定まる。また電子数Nも電子密度を全空間に渡って積分することで求めることができる。その外部ポテンシャルと電子数から導かれるハミルトニアンHのシュレーディンガー方程式を解けば、その外部ポテンシャルのもとで許される電子系の波動関数Ψがわかるので、あらゆる物理量をそこから求めることができる。つまり、基底状態の電子密度から、系の（励起状態に関わる量も含めて）あらゆる物理量は原理的には計算できることになる。物理量を電子密度から計算する方法を密度汎関数法というが、この定理はそれを正当化するものである。3次元空間内のN電子系の波動関数は各電子について3個、合計3N個の座標変数に依存する関数となる。一方、電子密度は電子が何個になろうとも3個の座標変数に依存するだけであり、取り扱い易さに雲泥の差がある。
また、ホーヘンベルク・コーンの第2定理によれば、外部ポテンシャルをパラメータにもつ電子密度の汎関数（ホーヘンベルク・コーンのエネルギー汎関数）が存在して、この汎関数は与えられた外部ポテンシャルのもとでの基底状態の電子密度で最小値を持ち、基底状態のエネルギーを与える。つまりの定義域のに対して
がなりたつ。よって電子密度関数を変化させて最小のエネルギーを与える電子密度を探索すれば基底状態の電子密度を求めることができる。

スピン密度汎関数理論

ホーヘンベルグ・コーンの定理を拡張して、スピン密度汎関数理論を得ることができる。
いまスピンの量子化軸をz方向にとり、その方向に外部磁場がかけられているとする。ハミルトニアンにゼーマン項を導入すると元来のホーヘンベルグ・コーンの第一定理と同様の議論で、外部ポテンシャルおよび外部磁場は基底状態の電子スピン密度の汎関数であることが示される。また同第二定理で示されているようなホーヘンベルグ・コーンのエネルギースピン密度汎関数も構成することができる。
スピン密度汎関数理論に置けるコーン・シャム理論の構成も容易である。この枠組みで、LDAに対応する交換相関エネルギーに対する近似は特に局所スピン密度近似 (Local Spin Density Approximation, LSDA) と呼ばれることもある。
しばしばスピン密度汎関数理論は密度汎関数理論と特に区別されずに呼ばれ、LSDAも単にLDAと呼ばれることが多い。

適用

実際にはコーン・シャム理論は調べる系に応じていくつかの異なった方法で用いられている。固体の計算では局所密度近似は平面波基底などを用いた手法で未だに使われている。これは電子気体からのアプローチが無限の大きさの固体に広がる非局在電子には適切であるためだと考えられる。しかし分子の計算ではより複雑な手法が必要となり、数多の交換-相関エネルギー汎関数が考えだされてきた。そのうちのいくつかは一様電子気体近似と相反するが、電子密度が一様となる極限ではLDAに帰着しなくてはならない。物理学者のあいだで、おそらくもっとも用いられている汎関数は修正の加えられたPerdew-Burke-Ernzerhofの汎関数であろう。これは自由電子気体のエネルギーを一般化勾配を用いてパラメータ化したもので、自由に決められるパラメーターを持たない。しかし、この方法は気体相の分子では熱量的に正確さを欠く。化学の分野でよく用いられるのはBLYP（Beckeの交換エネルギー表式とLee、Yang、Parrらの相関エネルギー表式を用いていることに由来する）である。B3LYPはさらによく使われるハイブリッド汎関数とよばれる種類の汎関数である。ハイブリッド汎関数では交換エネルギーの汎関数（B3LYPの場合はBeckeの交換汎関数を用いる）はハートリー・フォック理論の交換項と組み合わせられるが、B3LYPの場合3つのパラメーターによって交換相関汎関数が混合される。調整できるパラメーターは一般的にはいくつかの「練習用」の分子にフィッティングすることで決められる。このような汎関数を用いて得られた結果は大抵の場合十分に正確であるのだが、精度を改良するような系統的な手法は存在しない（このことは波動関数を用いた配置間相互作用や連結クラスター法といった伝統的な手法とは好対照である）。したがって、現在の密度汎関数理論のアプローチでは他の手法や実験の結果と比べないと計算の誤差を見積もることができない。

磁場の効果を取り入れるための一般化

これまで述べてきた理論はベクトルポテンシャル（すなわち磁場）が存在する場合にはそのまま用いることができず、状況に応じていくらかの破綻を生じることになる。そのような場合には基底状態の電子密度と波動関数の対応は失われる。磁場の効果を取り入れるための一般化の方法として電流密度汎関数理論 (CDFT) と磁場密度汎関数理論 (BDFT) の2つがあげられる。どちらの理論も交換-相関エネルギー汎関数を一般化して電荷密度以外の効果も取り入れる必要がある。VignaleとRasoltによって確立された電流密度汎関数理論では、汎関数は電荷密度と常磁性電流密度の両方に依存し、Salsbury, Grayce, Harrisらによって確立された磁場密度汎関数理論 (BDFT) では汎関数は電荷密度と磁場に依存し、磁場の形状に依存することもありえる。どちらの理論においてもLDAに相当する近似を超えるような手法が容易に実装できないという問題を抱えている。