以第一原理方法研究量子尺度效应.docVIP

物理雙月刊（廿七卷四期）2005年8月 PAGE 559 以第一原理方法研究量子尺度效應 對金屬表面場激發電流的影響 文/黃聖峰、李筱文、梁贊全、陳子亭 摘要 1928年R. H. Fowler與L. W. Nordheim提出其場激發電流的計算方法，此方法對金屬表面場激發電流提供很不錯的描述，但由於其場激發電流只與表面之功函數有關，故對磁性表面的場激發電流自旋偏極化行為則束手無策。我們發展了一套以局部密度泛涵近似法為基礎的場激發電流第一原理計算方法，並以之研究鐵(001)與(110)表面之場激發電流，本文將介紹此計算方法及相關的計算結果，並探討量子尺度效應對場激發電流的影響。 一、前言 關於場激發電流的理論計算，最早可追朔至1928年的Fowler-Nordheim理論[1]，利用簡單的量子井模擬金屬表面之位勢行為，並以一維自由電子氣來模擬其場激發電流的特性，此方法對金屬表面場激發電流提供不錯的描述，其著名的Fowler-Nordheim 方程式 為 (1) 其中與存在一個斜率約為的簡單關係，如圖一所示[2]。 圖一：鎢 圖一：鎢(116)、(111)、(211)及(110)表面之Fowler-Nordheim圖形及其功函數。 一般認為場激發電流主要是由費米能級附近的電子提供，若在費米能級附近自旋向上的表層局部能態密度與自旋向下的表層局部能態密度有差異時，即表示其穿隧電流具偏極性。自旋偏極化P定義為不同自旋方向電流之相差與總和之比值 (2) 1979年Campagna等人由實驗中得到鐵和鎳在不同表面自旋偏極化之數值如表1所示[3]。 表1 實驗測得鐵與鎳在不同表面之自旋偏極化行為 Spin Polzrization of the Field Emission Current Ni (hkl) P (%) Fe (hkl) P (%) 100 -3 ± 1 110 -5 ± 10 110 +26 ± 1 100 +26 ± 5 111 -6 ± 1 111 +20 ± 5 321 +7.3 ± 0.5 掃描式電子穿隧顯微鏡(STM)[4]讓我們可以得到原子尺度的表面結構。而自旋極化的掃描式電子穿隧顯微鏡(SPSTM)[5]的出現，使產生原子尺度解析之表面磁結構的夢想露出曙光。由於鎢的高熔點特質，大部分的STM探針皆以鎢製成，但是鎢沒有磁性，所以不能用以探測磁性材料表面之磁結構；鐵、鈷、鎳等3d過渡金屬具有磁性，但它們的熔點太低，不能做成很好的探針。若以少數層的鐵、鈷或鎳覆蓋在鎢的針尖上，則有可能產生既具磁性，熔點又高的探針，產生自旋偏極化的場激發電流，並用以探測磁性表面的磁結構。如此一來，對磁性表面及磁性薄膜場激發電流的理論計算便變得非常重要。 因為Fowler-Nordheim理論是基於一維自由電子氣模型，其場激發電流只與其功函數有關，故不能描述磁性表面場激發電流之自旋偏極化行為，也無法理解基板對磁性薄膜場激發電流自旋偏極化的影響。局部密度泛函近似法(Local Density Functional Approximation, LDA)能對非強相關性材料提供非常好的描述[6]，我們以LDA為基礎發展了一套場激發電流計算方法[7]，採用週期重複的厚層超晶格來模擬金屬表面行為，而模擬外加一垂直表面方向電場時，則以一平面電荷分布置於真空層中央調整其位勢，如圖二所示。 圖二：厚層超晶格之平均位勢示意圖， 圖二：厚層超晶格之平均位勢示意圖，(a)為一般金屬表面模擬計算，(b)為在真空層中央加上平面電荷模擬外加電場。 (a) (b) 透過自洽計算，我們可得到此金屬表面在外場下之表面原子結構、電子能態分佈與其對應的波函數、功函數以及表面電子所感受到的庫倫位勢等之精確計算結果。本計算方法是利用以上的資訊計算各能態對應的場激發電流。其原理很簡單，計算量不大且易於分析，可以探討不同k點及不同能階對電流的貢獻。我們利用此方法討論在弱電場下，金屬表面的場激發電流隨厚層層數變化的影響，並就不同k點觀察各個k點對電流的貢獻。而若將此方法應用在自旋偏極化的LDA計算中，即可分別處理不同自旋方向的場激發電流性質，進而研究其場激發電流的自旋偏極化行為。 當我們以N層的厚層來描述金屬表面時，垂直方向的週期性被破壞，所以其對應的能態被量化，這就是所謂的量子尺度效應[8]，此時表面的原子結構、功函數以及電子能態皆會隨著N的改變而改變。當然這些物理量皆會影響該表面的場激發電流，也就是說，場激發電流也會隨著N的改變而改變。由於量子尺度效應對自旋向上與向下電子的作用不一樣，所以場激發電流的自旋偏極化也會隨著厚層數而改變。本文