第二章 理論與計算方法
§2-1 固態材料的電子結構理論
§2-1-1 密度泛函理論
Density functional theory (DFT) 是一種研究多電子體系電子結構 的量子力學方法。密度泛函理論在物理和化學上都有廣泛的應用。古 典電子結構理論的近似方法，特別是 Hartree-Fock 方法(也有人稱為自 我一致場，self-consistent field method (SCF))，必須處理複雜的多電子 波函數。而密度泛函理論最主要的貢獻就是用電子密度取代波函數做 為研究的基本量，這取代了原來的多體(many-body)計算。因為多電 子波函數有 3N 個變數（N 為電子數，每個電子包含三個空間變 數），相較之下電子密度僅是三個變數的函數，無論在概念上還是實 際上都更方便處理。1964 年時，Hohenberg 與 Kohn 兩位學者提出 Hohenberg-Kohn 定理，他們發表系統基態的總能是電子密度的泛函，
系統的描述可由電子密度分布函數來決定。
4
隔年，Kohn 和 Sham 兩位學者從過去的 Hartree-Fock 與 Hohenberg -Kohn 理論中整理並提出了 Kohn-Sham 方程式，可用來估算系統總
5 和 GGA(Generalized gradient approximation)近似。
6
§ 2-1-2 局部密度近似法 (Local Density Approximation, LDA)
雖然在Kohn-Sham的密度泛函數理論架構下，多電子系統的基態
7
雖然在大部分的固態系統中都可得到滿意的結果，但像是半導體表面 上，對於金屬塊材計算或是表面化學反應，例如：分子鍵的生成或解 離，這種電子密度變化較大的系統下，LDA 近似法就失去其準確性，
故 1981 年 Gross、Driezler 及 Perdew 等學者提出廣義梯度近似法 (Generalized gradient approximation, GGA)，以彌補 LDA 不足之處。
§ 2-1-3 廣義梯度近似法 (Generalized Gradient Approximation,
GGA)
8
§ 2-1-4 空間週期性 (periodic boundary condition)
根據能帶理論，固態系統中的電子不再受限只能在原位置附近運 動，而可以在整個系統內部自由移動。由於固態材料具有空間週期 性，因此我們在電腦模擬時可以利用它的特性來簡化計算。簡單來 說，我們可利用其最小重複單元，即單位晶格 (unit cell) 來模擬在空 間中無限延伸的材料，只要可以解出單位晶格中的電子特性，就可了 解整個材料的性質。如圖2-1，假設原子處在 r 時的位能為 U，因空 間週期性的關係，因此在下一個單位晶格的 r + a 位置時，位能同樣 是U。因此我們在處理單位晶格的計算時，也等同於處理無限延伸之 真實系統。
圖 2-1 單位晶格的無限延伸，其中 a 為單位晶格的大小
9
§ 2-1-5 布洛赫定理(Bloch Theorem)
藉由 Bloch Theorem 可用來說明了一個空間週期性系統的物理量
10
11
§ 2-1-6 虛位勢 (pseudopotential)
對於一個 Ab Initio 計算而言，想要得到精確的結果就必須要計算
12 的問題，可以用一個平滑的虛波函數
pseudo( )r (pseudo-wavefunction) 來代替，有了這個虛波函數之後，我們可以再從薛丁格方程式得13
，此時再利用薛丁格方程式的反置(Inverted Schrödinger equation)來解出虛位勢Vpseudo( )r 來。一般在固態材料的計算中常看14
§ 2-1-7 VASP 計算軟體
VASP (Vienna Ab-initio Simulation Package) 是使用虛擬電位勢與平面 波基組進行計算之軟體，其採用週期性邊界條件處理原子、分子以及表面 體系和固體。此軟體應用在我們研究中，可計算材料的結構參數如鍵長、
晶格常數及原子位置等，也可計算材料的電子結構如電子態密度
。
進行 VASP 的計算時，必須在資料夾中準備下列四個檔案: POSCAR、KPOINTS、POTCAR、INCAR。其中 POSCAR 寫入各原子的座標位置，
KPOINTS 則決定了 k 點的取法，POTCAR 包含了各原子所使用的虛位勢，
INCAR 則是計算條件的設定。
15
§2-2 擾動彈簧模型(Nudged Elastic Band; NEB)
在一個多維度的參數系統中，存在很多局部的能量低點 (Local minimum)，想從一個能量低點走到另一個能量低點時，必須爬過一個 能障，不同的反應路徑所通過的能障也不一樣。而我們最感興趣的是如 何找到最低能障路徑。
圖 2-3 二維位勢圖
擾動彈簧模型(Nudged Elastic Band; NEB)即是為了解決此種問題而發 展出來的方法。圖 2-3 中的 R 點與 P 點為局部的最小值(Local minima)，
為了找出從反應物 R 變成產物 P 這個反應路徑的最小能障，首先將 R 與 P 當作初始態與最終態，而 R 與 P 之間的虛線是用來連接初始態和 最終態結構，接著利用等分法將這條虛線切為四等分，其中包含三個插
16
入的映像(image)，再來將這三個映像(image)最佳化後得到其最小的能 量路徑，即 R 與 P 之間的實線路徑，而這條實線路徑會有一個馬鞍點 (Saddle point)，R 點(初始態)和馬鞍點之間的能量差即這個反應的能 障，如圖 2-4 所示。
圖 2-4 位能曲面圖
17
§2-3 態密度(Density of state, DOS)
態密度(Density of state, DOS)為一個相當重要的電子能態分部 分析工具，不管在化學或是物理上都佔有相當的份量。對於物理學家而 言，可以藉由態密度分析塊材的物理性質，而對於化學家而言，更是可 以將分子軌域的概念延伸到態密度中，藉此來幫助處理表面化學的問 題。而一般態密度分為總態密度(Total density of state，TDOS)與部分態 密度(Partial density of state，PDOS)，部分態密度主要是根據選取的原 子來做投影，對於用來分析鍵結變化有相當大的幫助。
Density (states/eV) gaseous ethanol
圖 2-5 氣相乙醇分子的態密度
18
§2-4 電子局域化函數(Electron localization function, ELF)
為了討論 Fe(111)、W(111)，與 Fe-W(111)雙金屬表面上的電荷分布，
我們利用電子局域化函數(ELF)分析表面的變化情況。在 VASP 的計算當中 將參數 LELF = .TRUE.開啟就可以計算出 ELF 值，在文中標記為 Z 值，再 利用軟體尋找所需要的切面，可以得到該切面的 ELF 分布圖。
電子局域化函數為 Becke 和 Edgecombew 為了觀測原子和分子之間的電 子局域化所提出的方法4-6，其計算出的 ELF 值(Z)介於零到一之間。當 Z = 1，表示電子為完全的局域化(perfect localization);反之，若 Z = 0 則表示電 子為非局域化(delocalization)，藉著分析 Z 值的大小可以了解到化學鍵的組 成與電子局域化的現象。
19 純金屬會有不同的結果；而塊材排列方式(ensemble effect) 6的影響則是在於 當不同晶格常數的金屬混合後會得到一個新的晶格，所以會造成金屬原子 彼此間的距離和純金屬時不一樣，而對於分子的吸附造成不同的效果。近
來文獻13-21指出將反應性較差的金屬(如 Cu、Ag、Au)與反應性較好的金屬
(如 Pd、Rh)混和可得到比純金屬還要有更好的催化效果。常見的有 RhCu13,16