密度泛函理論

**密度泛函理論** （英語：density functional theory，簡稱 DFT）是一種研究多電子體系電子結構的量子力學方法。密度泛函理論在物理和化學上都有廣泛的應用，特別是用來研究分子和凝聚體的性質，是凝聚體物理和計算化學領域最常用的方法之一。

## 理論概述

電子結構理論的古典方法，特別是 Hartree-Fock 方法和後 Hartree-Fock 方法，是基於複雜的多電子波函數的。密度泛函理論的主要目標就是用電子密度取代波函數做為研究的基本量。因為多電子波函數有 $3N$ 個變量（ $N$ 為電子數，每個電子包含三個空間變量），而電子密度僅是三個變量的函數，無論在概念上還是實際上都更方便處理。

雖然密度泛函理論的概念起源於 Thomas-Fermi 模型，但直到 Hohenberg-Kohn 定理提出之後才有了堅實的理論依據。Hohenberg-Kohn 第一定理指出體系的基態能量僅僅是電子密度的泛函。

Hohenberg-Kohn 第二定理證明了以基態密度為變量，將體系能量通過變分得到最小值之後就得到了基態能量。

HK 理論最初只適用於沒有磁場存在的基態，現在已經被推廣。最初的 Hohenberg-Kohn 定理僅僅指出了一一對應關係的存在，但是沒有提供任何這種精確的對應關係。正是在這些精確的對應關係中存在著近似（這個理論可以被推廣到時間相關領域，從而用來計算激發態的性質 [6]）。

密度泛函理論最普遍的應用是通過 Kohn-Sham 方法實現的。 在 Kohn-Sham DFT 的框架中，複雜的多體問題（由於處在一個外部靜電位中的電子交互作用而產生的）被簡化成一個沒有交互作用的電子在有效勢場中運動的問題。這個有效勢場包括了外部勢場以及電子間庫侖交互作用的影響，例如交換和關聯作用。處理交換關聯作用是 KS DFT 的難點，目前尚沒有精確求解交換相關能 $E_{XC}$ 的方法。最簡單的近似求解方法是局域密度近似 (LDA)。LDA 近似用均勻電子氣來計算體系的交換能（均勻電子氣的交換能是可以精確求解的），而採用對自由電子氣進行擬合的方法來處理關聯能。

自 1970 年以來，密度泛函理論在固態物理學計算中得到廣泛的應用。多數情況下，與其它解決量子力學多體問題的方法相比，採用局域密度近似的密度泛函理論給出了非常令人滿意的結果，同時固態計算相比實驗的費用要少。儘管如此，人們普遍認為量子化學計算不能給出足夠精確的結果，直到二十世紀九十年代，理論中所採用的近似被重新提煉成更好的交換關聯作用模型。密度泛函理論是目前多種領域中電子結構計算的領先方法。 密度泛函理論儘管得到改進，但是描述分子間作用力 ，特別是范德華力，或者計算半導體的能隙還有一定困難。

## 早期模型: Thomas-Fermi 模型

密度泛函理論可以上溯到由 Thomas 和 Fermi 在 1920 年代發展的 _Thomas-Fermi_ 模型。他們將一個原子的動能表示成電子密度的泛函，並加上原子核 - 電子和電子 - 電子交互作用（兩種作用都可以通過電子密度來表達）的古典表達來計算原子的能量。

_Thomas-Fermi_ 模型是很重要的第一步，但是由於沒有考慮 Hartree-Fock 理論指出的原子交換能，Thomas-Fermi 方程式的精度受到限制。1928 年保羅・狄拉克在該模型基礎上增加了一個交換能泛函項。

然而，在大多數應用中 Thomas-Fermi-Dirac 理論表現得非常不夠準確。其中最大的誤差來自動能的表示，然後是交換能中的誤差，以及對電子相關作用的完全忽略。

## 導出過程和表達式

在通常的多體問題電子結構的計算中，原子核可以看作靜止不動的（波恩 - 歐本海默近似），這樣電子可看作在原子核產生的靜電位 $\,\!V$ 中運動。電子的定態可由滿足多體薛丁格方程式的波函數 $\Psi ({\vec {r}}_{1},\dots ,{\vec {r}}_{N})$ 描述：

　
: $H\Psi =\left[{T}+{V}+{U}\right]\Psi =\left[\sum _{i}^{N}-{\frac {\hbar ^{2}}{2m}}\nabla _{i}^{2}+\sum _{i}^{N}V({\vec {r}}_{i})+\sum _{i<j}U({\vec {r}}_{i},{\vec {r}}_{j})\right]\Psi =E\Psi $

其中 $\,\!N$ 為電子數目， $\,\!U$ 為電子間的交互作用勢。算符 $\,\!T$ 和 $\,\!U$ 稱為普適算符，它們在所有系統中都相同，而算符 $\,\!V$ 則依賴於系統，為非普適的。可以看出，單粒子問題和比較複雜的多粒子問題的區別在於交換作用項 $\,\!U$ 。目前有很多成熟的方法來解多體薛丁格方程式，例如：物理學裡使用的圖形微擾理論和量子化學里使用的基於斯萊特行列式中波函數系統展開的組態交互作用（CI）方法。然而，這些方法的問題在於較大的計算量，很難用於大規模複雜系統的計算。

相比之下，密度函理論將含 $\,\!U$ 的多體問題轉化為不含 $\,\!U$ 的單體問題上，成為解決此類問題的一個有效方法。在密度泛函理論中，最關鍵的變量為粒子密度 $n({\vec {r}})$ ，它由下式給出

　
: $n({\vec {r}})=N\int {\rm {d}}^{3}r_{2}\int {\rm {d}}^{3}r_{3}\cdots \int {\rm {d}}^{3}r_{N}\Psi ^{*}({\vec {r}},{\vec {r}}_{2},\dots ,{\vec {r}}_{N})\Psi ({\vec {r}},{\vec {r}}_{2},\dots ,{\vec {r}}_{N}).$

皮埃爾・奧昂貝格和沃爾特・科恩在 1964 年提出，上面的關係可以反過來，即給出基態電子密度 $n_{0}({\vec {r}})$ ，原則上可以計算出對應的基態波函數 $\Psi _{0}({\vec {r}}_{1},\dots ,{\vec {r}}_{N})$ 。也就是說， $\,\!\Psi _{0}$ 是 $\,\!n_{0}$ 的唯一泛函，即

　
: $\,\!\Psi _{0}=\Psi _{0}[n_{0}]$

對應地，所有其它基態可觀測量 $\,\!O$ 均為 $\,\!n_{0}$ 的泛函

　
: $\left\langle O\right\rangle [n_{0}]=\left\langle \Psi _{0}[n_{0}]\left|O\right|\Psi _{0}[n_{0}]\right\rangle .$

進而可以得出，基態能量也是 $\,\!n_{0}$ 的泛函

　
: $E_{0}=E[n_{0}]=\left\langle \Psi _{0}[n_{0}]\left|T+V+U\right|\Psi _{0}[n_{0}]\right\rangle $ ,

其中外勢場的貢獻 $\left\langle \Psi _{0}[n_{0}]\left|V\right|\Psi _{0}[n_{0}]\right\rangle $ 可以用密度表示成

　
: $V[n(r)]=\int V({\vec {r}})n({\vec {r}}){\rm {d}}^{3}r.$

泛函 $\,\!T[n(r)]$ 和 $\,\!U[n]$ 稱為普適泛函，而 $\,\!V[n]$ 顯然不是普適的，它取決於所考慮的系統。對於確定的系統，即 $\,\!V$ 已知，需要將泛函

　
: $E[n(r)]=T[n(r)]+U[n(r)]+\int V({\vec {r}})n({\vec {r}}){\rm {d}}^{3}r$

對於 $n({\vec {r}})$ 求極小值。這裡假定能夠得出 $\,\!T[n(r)]$ 和 $\,\!U[n(r)]$ 的表達式。對能量泛函求極值可以得到基態電子密度 $\,\!n_{0}$ ，進而求得所有基態可觀測量。

對能量泛函 $\,\!E[n(r)]$ 求變分極值可以用不定算子的拉格朗日方法，這由科恩和沈呂九在 1965 年完成。這裡我們使用如下結論：上面方程式中的泛函可以寫成一個無交互作用的體系的密度泛函

　
: $E_{s}[n(r)]=\left\langle \Psi _{s}[n]\left|T_{s}+V_{s}\right|\Psi _{s}[n(r)]\right\rangle ,$

其中 $\,\!T_{s}$ 為無交互作用的動能， $\,\!V_{s}$ 為粒子運動感受到的外勢場。顯然， $n_{s}({\vec {r}})\equiv n({\vec {r}})$ ，若 $\,\!V_{s}$ 取為

　
: $V_{s}=V+U+\left(T-T_{s}\right).$

這樣，可以解這個輔助的無交互作用體系的科恩 - 沈呂九方程式

　
: $\left[-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}+V_{s}({\vec {r}})\right]\phi _{i}({\vec {r}})=\epsilon _{i}\phi ({\vec {r}}),$

可以得到一系列的電子軌域 $\,\!\phi _{i}$ ，並由此求得原來的多體體系的電子密度 $n({\vec {r}})$

　
: $n({\vec {r}})\equiv n_{s}({\vec {r}})=\sum _{i}^{N}\left|\phi _{i}({\vec {r}})\right|^{2}.$

等效的單粒子勢 $\,\!V_{s}$ 可以表示成

　
: $V_{s}=V+\int {\frac {e^{2}n_{s}({\vec {r}}\,')}{|{\vec {r}}-{\vec {r}}\,'|}}{\rm {d}}^{3}r'+V_{\rm {XC}}[n_{s}({\vec {r}})],$

其中第二項為描述電子間庫侖斥力的哈特里項，最後一項 $\,\!V_{\rm {XC}}$ 叫做交換關聯勢，包含所有多粒子的交互作用。由於哈特里項和交換關聯項 $\,\!V_{\rm {XC}}$ 都依賴於 $n({\vec {r}})$ , $n({\vec {r}})$ 又依賴於 $\,\!\phi _{i}$ , 而 $\,\!\phi _{i}$ 又依賴於 $\,\!V_{s}$ , 科恩 - 沈呂九方程式的求解需要用自洽方法。通常首先假設一個初始的 $n({\vec {r}})$ , 然後計算對應的 $\,\!V_{s}$ 並求解科恩 - 沈呂九方程式中的 $\,\!\phi _{i}$ 。進而可以計算出新的密度分布，並開始新一輪計算。此過程不斷重複，直到計算結果收斂。

## 參考資料

[1] P. Hohenberg and W. Kohn, Phys. Rev. _136_ (1964) B864 

[2] W. Kohn and L. J. Sham, Phys. Rev. _140_ (1965) A1133 

[3] A. D. Becke, J. Chem. Phys. _98_ (1993) 5648 

[4] C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B _37_ (1988) 785 

[5] P. J. Stephens, F. J. Devlin, C. F. Chabalowski, and M. J. Frisch, J. Phys. Chem. _98_ (1994) 11623

## 相關閱讀

* Klaus Capelle, A bird's-eye view of density-functional theory（頁面存檔備份，存於網際網路檔案館）