コーン–シャム方程式

コーン–シャム方程式（コーン–シャムほうていしき、テンプレート:Lang-en）とは、量子化学の、特に密度汎関数理論で用いられる、相互作用のある粒子（典型的には電子）からなる任意の既知の系と同じ密度を生成する、相互作用のない粒子からなる仮想的な系（コーン–シャム系 (テンプレート:Lang)）のシュレーディンガー方程式のことである。コーン–シャム方程式は、相互作用のない粒子がその中を動く局所的な有効（仮想）外部ポテンシャルによって定義される。このポテンシャルはコーン–シャムポテンシャル (テンプレート:Lang) と呼ばれ、典型的には テンプレート:Math や テンプレート:Math と表される。コーン–シャム系中の粒子は相互作用のないフェルミオンであるため、コーン–シャム波動関数 (テンプレート:Lang) は方程式

$(-\frac{{\hslash }^2}{2m}{\mathrm{\nabla }}^2+v_{\mathrm{e}\mathrm{f}\mathrm{f}}(𝒓)){\varphi }_i(𝒓)={\epsilon }_i{\varphi }_i(𝒓)$

の最低エネルギー解である軌道の集合から構築される単一スレイター行列式となる。この固有値方程式はコーン–シャム方程式の典型的な表現である。ここで、テンプレート:Math はコーン–シャム軌道 (Kohn–Sham orbital) テンプレート:Math に対応した軌道エネルギーである。テンプレート:Mvar 粒子系に対する密度は次の式で与えられる：

$\rho (𝒓)=\sum _i^N|{\varphi }_i(𝒓){|}^2$

コーン–シャム方程式という名称は、1965年にカリフォルニア大学サンディエゴ校でこの概念を導入したウォルター・コーンとリュウ・シャム（Lu Jeu Sham, 沈呂九）の2人の名に因む。

コーン–シャムポテンシャル

密度汎関数理論において、系の全エネルギーは電荷密度の汎関数として表される：

$E[\rho ]=T_{\mathrm{s}}[\rho ]+\int \mathrm{d}𝒓v_{\mathrm{e}\mathrm{x}\mathrm{t}}(𝒓)\rho (𝒓)+V_{\mathrm{H}}[\rho ]+E_{\mathrm{x}\mathrm{c}}[\rho ]$

ここで テンプレート:Math はコーン–シャム運動エネルギー (Kohn–Sham kinetic energy) であり、コーン–シャム軌道を用いて表される：

$T_{\mathrm{s}}[\rho ]=\sum _{i=1}^N\int \mathrm{d}𝒓{\varphi }_i^{*}(𝒓)(-\frac{{\hslash }^2}{2m}{\mathrm{\nabla }}^2){\varphi }_i(𝒓)$

テンプレート:Math は相互作用のある系に作用する外部ポテンシャル（少なくとも、分子系においては、電子-原子核相互作用）であり、テンプレート:Math はハートリーエネルギー（クーロンエネルギー）

$V_{\mathrm{H}}=\frac{e^2}{2}\int \mathrm{d}𝒓\int \mathrm{d}{𝒓}^{\prime }\frac{\rho (𝒓)\rho ({𝒓}^{\prime })}{|𝒓-{𝒓}^{\prime }|}$

であり、テンプレート:Math は交換相関エネルギーである。コーン–シャム方程式は、全エネルギー表現を軌道の集合について変化させることによって求められ、コーン–シャムポテンシャルを次の形

$v_{\mathrm{e}\mathrm{f}\mathrm{f}}(𝒓)=v_{\mathrm{e}\mathrm{x}\mathrm{t}}(𝒓)+e^2\int \frac{\rho ({𝒓}^{\prime })}{|𝒓-{𝒓}^{\prime }|}\mathrm{d}{𝒓}^{\prime }+\frac{\delta E_{\mathrm{x}\mathrm{c}}[\rho ]}{\delta \rho (𝒓)}$

として得る。ここで右辺の最後の項

$v_{\mathrm{x}\mathrm{c}}(𝒓)\equiv \frac{\delta E_{\mathrm{x}\mathrm{c}}[\rho ]}{\delta \rho (𝒓)}$

は交換相関ポテンシャルである。この項、そして対応するエネルギー表現は、密度汎関数理論に対するコーン–シャムのアプローチの中で未知となっている唯一のものである。軌道を変化させない近似がハリス汎関数理論である。

コーン–シャム軌道エネルギー テンプレート:Math は、一般に、物理的な意味をあまり持たない（クープマンズの定理参照）。軌道エネルギーの合計は全エネルギーに関係しており、

$E=\sum _i^N{\epsilon }_i-V_{\mathrm{H}}[\rho ]+E_{\mathrm{x}\mathrm{c}}[\rho ]-\int \frac{\delta E_{\mathrm{x}\mathrm{c}}[\rho ]}{\delta \rho (𝒓)}\rho (𝒓)\mathrm{d}𝒓$

軌道エネルギーが、より一般的な制限テンプレート:仮リンクの場合において一意でないため、この方程式は軌道エネルギーの特定の選択についても適用できる（クープマンズの定理参照）。