挑戰激發態：多體微擾理論之GW近似
文/薛宏中
當第一原理(first-principles)方法，在計算材料的基態性質上大放異彩之後，我們開始期待征服另一座高山—激發態。事實上，許多重要的實驗，都與材料電子結構的激發態關係密切。然而，材料中多體效應與激發態之關係為何？現階段之基態理論有何不足之處？何謂激發態之GW計算？新的激發態計算方法真的是救世主嗎？GW計算的未來發展又為何？本文在儘量不涉及複雜數學語言的情況下，將一一探討這些問題。
 |
一、材料中如影隨形的多體效應
現代人總是強調人際關係的重要，思考如何加強人與人之間互動與交流；在微觀的原子世界裡，也存在著相似難題。古希臘文明即揭櫫物質構成之基本元素─原子；然而，單就材料中如天文數字般（每莫耳約6X1023個）的原子數目，就已經讓物理學家對材料的微觀研究捉襟見肘了，更遑論原子間還存在著複雜的交互作用；但是，這些障礙並非完全無法克服：在1950年以前，理論學家利用正則轉換(canonical transformation)，嘗試將具有交互作用之多體系統，轉換成多個幾乎無關之獨立系統，以降低複雜程度。另一方面，1956年之後發展的量子場論(quantum field theory)中，以近乎完全獨立之準粒子(quasi- particle)，來近似在多體效應下之真實粒子(或稱為裸粒子(bare particle)，再利用費曼圖(Feynman diagram) ［註：一種非常清晰直覺的圖解分析工具，按照費曼法則（Feynman Rule），理論物理學者可以準確而容易地進行一系列複雜的量子場論計算！］計算這些準粒子的特性。
那到底什麼是準粒子呢？其實準粒子包含兩部份：裸粒子及交互作用在此粒子周圍所形成之效應；就好像是多體效應為裸粒子披上一件外衣一樣，遮蔽了原來裸粒子之間較強的交互作用。如此一來，這些”穿上衣服”的裸粒子即便在多體環境之下，也可視為獨立粒子。［註：Richard
D. Mattuck[1]對準粒子有
圖一、(a)獨立(裸)粒子與”獨立馬”；(b)準粒子與”準馬”？！
一個更生動的比擬：準粒子就如同在黃沙中絕塵而去的駿馬，你只能看到被塵土＂遮蔽＂下，約略的”馬形”，而無法細究此名駒之毛色與肌理！(如圖一)］舉個例子：將一個帶負電的電子，放進一電中性環境中(包括均勻分佈的眾多負電荷所構成之電子氣，與相同數量正電荷所形成之背景)。此外加電子將藉由庫侖斥力，將其周圍的電子氣排開；而在此處分佈之背景正電荷，即猶如一帶正電的電子雲，圍繞著外加電子形成準電子(quasi-electron)。此一動態的產生機制，預言了準電子是有生命期限的(life time)；而在電子之間原有的庫侖斥力，也因遮蔽效應可以幾乎忘了它的存在！因此準電子與獨立電子(裸電子)的特性(例如特徵能量)，存在著相當的差異。帶著多體效應基因的準粒子圖像，不僅可應用於原子核物理、固態物理、鐵磁性物理等多體系統之基態(ground-state)研究，更能適當地描述激發態(excited-state)特性(例如各種激發子：exciton、plasmon、magnon等等)。
另一方面，奠基於密度泛函理論(Density Functional Theory─DFT)的第一原理電子結構計算，近年來蓬勃發展(請參考本刊文章)，成功地拓展了在原子尺寸下研究材料科學精確度與視野，也提供人們預測與設計新穎未來材料的能力。過去，各種重要材料的發現與發明(例如：半導體、超導體、奈米材料等等)，都是靠著理論計算與精密實驗的通力合作。與實驗的配合上，計算材料科學中最重要的課題之一，就是如何更準確描述多電子系統其電子結構之基態與激發態。因此，結合第一原理計算與準粒子圖像(多體效應)，將是材料科學研究中的一項利器。以下，筆者將藉本篇短文，介紹第一原理材料計算中，發展最成功的激發態計算方法：GW近似法。
二、基態密度泛函理論的困境
在介紹第一原理激發態計算之前，我們先談談密度泛函理論[2]。以密度為參數，來描述量子力學系統，大概可追溯至Thomas-Fermi模型(1927提出)；然而，時下流行的第一原理計算，所依據之密度泛函理論，卻是濫觴於Hohenberg及Kohn于1964所提出之理論：多電子系統中的所有性質，都可視為系統基態電子密度函數r0(r)的函數(或稱為泛函)。至於這個抽象的函數確切形式，則是由Kohn與Sham在1965年才決定。Kohn-Sham方法是利用等效之獨立電子系統(具備等效位勢effective potential)，來取代原多體問題；而其中攸關電子間之多體效應，則以交換相干(exchange-correlation)泛函來描述。在實際計算上，通常以局部之密度泛函(例如Local Density
Approximation—LDA)，近似交換相干項。這一頗具巧思的安排，使得Kohn-Sham方程式能透過自洽(self consistent)的方法，有效地找到多電子系統”真正”(exact)的基態電子密度，與對應之基態總能(total energy)；進而讓此一方法，在材料基態特性(例如：鍵長、結合能、聲子譜、等等)計算上，獲得空前的成功。事實上，這樣的成果的確有點出人意表(嚴格來說，以均勻電子氣模型來描述交換相干項，是過於粗糙的近似)。當我們為了與實驗結果比較，而企圖將Kohn-Sham方法推廣到激發態計算時，卻遇到挫敗；最著名的例子，就是半導體與絕緣體之電子結構能帶(band structure)計算：雖然，自洽Kohn-Sham方法(LDA近似)所得之能帶色散(energy dispersion)行為，與實驗結果接近；然而，其計算之主能隙(fundamental band-gap)值(如圖二(a))，卻都嚴重低估。此外，所得之能帶寬度(band width)，也與實驗結果相差甚遠。另一方面，密度泛函理論計算之光學頻譜(optical spectrum)，更是與實驗測量結果大相逕庭；此乃肇因於密度泛函理論之自洽Kohn-Sham方法，
圖二、(a)各種半導體之主能隙值：密度泛函理論LDA計算(橘色)、GW近似計算(藍色)，與實驗值(粉紅色直線)之比較；(b) 鍺(Ge)之能帶結構(摘錄自[4])：實驗量測(實心與空心符號)，與GW近似計算(黑實線)。
僅對於系統基態密度有嚴格的定義(事實上，Kohn-Sham本徵能量只是對應於虛擬之獨立電子，而非擁有多體效應內涵之準粒子)。因此，密度泛函理論也被稱為基態理論。顧名思義，對上述激發態性質，Kohn-Sham方法完全是束手無策。這對實驗物理學家而言，無疑是澆了一盆冷水，因為許多新奇且複雜的物理機制，都與系統激發態性質密切相關。正因如此，引發了理論物理學家發展第一原理激發態計算的強烈動機。有鑑於過去第一原理的成功經驗，激發態計算大多根據基態密度泛函理論做修正，特別是針對交換相干項中之相干泛函(相較於此一具備非局域性(non-local)及多體效應的函數而言，我們對於代表電子自旋之間交互作用的交換能量(exchange energy)形式，則是清楚多了)。近年來，第一原理激發態計算的發展態勢，有如風起雲湧一般。根據理論背景之嚴謹程度，及計算代價之高低，激發態計算包括以下形式：依據哈伯(Hubbard)模型而外加高度局域性作用U的LDA+U計算法、自我作用修正法(Self-Interaction
Correction)、計算準粒子之GW近似法[3]、與時間有關之密度泛函理論(Time-Dependent DFT)[3]、考慮非局域性泛函的平均密度近似(Average Density
Approximation)與權重密度近似法(Weighted Density Approximation)、等等。各種方法皆有其利弊，限於篇幅，無法一一詳述。接下來，筆者將探討其中之GW近似計算。
三、平均場近似與多體微擾理論GW計算
在多電子系統中，準粒子(即為準電子)與獨立粒子之間最關鍵的差異，就是多體作用之”遮蔽效應”。因此，多體計算就是要處理那件披在裸電子身上的”外衣”。然而有些問題，在深入細節之前，必須先釐清：何種函數可以適切的表達準粒子呢？如何計算遮蔽效應的強弱？遮蔽效應對於系統能量的改變又是如何呢？
格林函數(Green’s function)[2]
在量子場論的數學語言當中，通常使用格林函數(以G代表)計算：當粒子由初始狀態，穿過多體系統後到達新的狀態時，所有可能路徑的機率期望值，進而可計算其頻譜函數(spectral function)與生命期限。而其中單粒子(single-particle)格林函數，可用來描述N個粒子系統，經由粒子的產生(或湮滅)，演進(propagate)到N+1個粒子系統(或N-1個粒子系統)時，多粒子系統基態期望值。此一圖像，恰巧與前述之準電子圖像(外加電子，進入N個電子所形成之電子氣中，造成準電子)相似；就像是專為其量身打造的數學模型一樣，格林函數當然是描述準粒子的首選。［註：若考慮系統總粒子數守恆，也就是同時產生及湮滅粒子(例如成對的電子─電洞(electron-hole))時，可採用雙粒子(two-particle)格林函數］
介電函數(Dielectric function)
在基礎電磁學中，我們知道當電磁波進入不同材料時，會因材料特性(可視為材料中，所有電子對入射電磁波的集體響應(response))不同而改變其強度，並以材料介電係數(dielectric constant)的大小，來此估計此效應。我們也可以同樣的想法，度量準粒子間遮蔽效應的強弱。例如：在準電子之間，”遮蔽”庫侖位勢(screened Coulomb potential，通常以W代表)可視為原獨立電子(裸電子)間庫侖位勢之修正結果(原庫侖位勢除以介電函數)。因此，如何準確且有效的決定系統之介電函數，將成為進行準粒子計算的重要工作之一。
自我作用能(Self-energy)
在準電子圖像中，我們曾經提及其與獨立電子之間的差異。事實上，當外加電子與原均勻分佈電子氣發生交互作用，而在其自身周圍產生一帶正電之電子雲時，此電子雲也會反過來影響外加電子的狀態，而產生能量改變；就好像是外加電子，透過遮蔽在周遭的電子雲，作用在自己身上。此一能量變化，就稱為自我作用能(通常以S代表)。因此，自我作用能S就是多體效應在系統總能上的貢獻。依據密度泛函的經驗，我們將多電子系統的總能量，切割成獨立電子系統的能量，加上電子的自我作用能(對應於密度泛函中的交換相干能)。只要能準確計算自我作用能，此多電子系統便可迎刃而解。因為牽涉到遮蔽效應，自我作用能S當然與W有關。然而，想要解決具有非局欲性、與能量相關且複雜的S泛函，的確相當棘手。
No comments:
Post a Comment