邱榮科,楊凱麗
(沈陽工業大學 理學院,沈陽 110870)
表面科學是物理學的一個重要分支,也是自然科學中發展最快的一個交叉學科.對表面性質具有重要影響的現象為集體激發效應,例如磁性薄膜中的自旋微擾可以引起表面自旋波的傳播.關于磁性薄膜表面和界面磁性的研究對提高記錄密度和改善許多磁性儀器的性能具有重要作用.沿磁性薄膜表面傳播的表面自旋波具有大量顯著特征,可以應用于自旋波電子學領域.磁性薄膜的動態特性和自旋波特性是自旋電子學的重要組成部分,且與納米磁記憶器件、傳感器和邏輯器件等新型自旋電子學器件的性能密切相關[1-2].磁性薄膜表現出的新特性引起了人們的極大興趣.許多研究者利用自旋波共振、光散射等實驗手段對磁性薄膜的結構、磁性等進行了探討[3-4].相關線性理論[5]表明,對于表面自旋自由的純交換耦合鐵磁體而言,并不存在表面自旋波.人們從動力學觀點出發,利用轉移矩陣方法和相關函數理論等對磁性薄膜的動力學性質、自旋波譜和色散關系進行了研究[6-7].Mamica[5]研究表明,在磁性薄膜中傳播的交換自旋波表現出一些動力學效應,即僅由波傳播引起的效應主要集中在四個現象上:表面(和次表面)局部化、體能帶縮為一個能級、自旋波譜中模序的反轉以及薄膜動態分離為兩個子系統.王煥等[8]研究了由鐵磁與反鐵磁材料構成的異質雙層薄膜結構在界面耦合處的界面自旋波及其相關物理特性,得到界面自旋波的存在條件與色散關系.Flaviano等[9]采用絕熱近似方法分析了W(110)表面的自旋波色散現象,并利用格林函數的第一原理電子結構分析了磁交換相互作用.張光富等[10]利用微磁學模擬方法研究了磁性納米薄膜末端形狀對自旋波模式特性的影響,結果表明,納米膜的末端邊緣形狀對自旋波模式特性影響不大,但隨著磁場的增大,自旋波頻率呈現不同變化.本文試圖從理論角度采用量子格林函數方法研究溫度、外磁場、表面交換耦合和表面各向異性對鐵磁性薄膜結構中表面自旋波頻率的影響.
采用具有單格點各向異性的海森堡模型研究鐵磁性薄膜中表面自旋波的頻率.鐵磁性薄膜中含有12個原子層并沿著z方向堆砌,其中,鐵磁性薄膜的平面結構與表面布里淵區如圖1所示.

圖1 鐵磁性薄膜平面結構與表面布里淵區Fig.1 Planar structure and surface Brillouin zone of ferromagnetic thin film
哈密頓量可以表示為
(1)


(2)
式中:S+和S-分別為自旋產生和湮滅算符;a為常數;ω為自旋波頻率.
采用Tyablikov退耦近似方法處理交換耦合項,采用Anderson-Callen退耦近似方法處理單格點各向異性項.建立格林函數的運動方程,其中格林函數的奇異點對應自旋波譜的解,具體建立過程參見文獻[12].在磁性薄膜表面自旋波頻率的計算中,自旋量子數與薄膜中間磁性離子的最近鄰交換耦合常數均取值為1,為了簡化計算,其他物理量均取與交換耦合常數作比值后獲得的約化數值(無量綱單位).
圖2為溫度對鐵磁性薄膜中表面自旋波頻率的影響.在計算過程中相關參數取值為:B=0.15、Ds=0.05、D=1.在布里淵區表面自旋波存在的Γ-K區域選取k1=4.188 8、k2=3.708 6、k3=3.228 5的三個波矢進行分析.圖2a中Js=0.75,此時薄膜中存在聲學表面自旋波,而圖2b中Js=2,此時薄膜中存在光學表面自旋波.此外,圖2橫坐標為約化溫度且τ=T/Tc,本文中聲學表面自旋波對應的Tc為系統參數S=1、Js=0.75、J=1、B=0.15、Ds=0.05、D=1時薄膜的居里溫度,而光學表面自旋波對應的Tc為系統參數S=1、Js=2、B=0.15、Ds=0.05、D=1時薄膜的居里溫度.
由圖2可見,聲學和光學表面自旋波頻率均隨著溫度的增加而降低.文獻[13]表明,當溫度升高時會發生自旋波頻率軟化現象,這與本文結論一致.由圖2a可見,在低溫區波矢越大,聲學表面自旋波頻率越大;在高溫區隨著波矢的增大,自旋波頻率增大幅度不明顯.由圖2b可見,在低溫區隨著波矢的增大,光學表面自旋波頻率隨之增大;而在高溫區隨著波矢的增大,自旋波頻率增加幅度同樣不明顯.對比圖2a、b可知,當約化溫度為0.3~0.5時,光學表面自旋波頻率隨著溫度的增加而急速降低.

圖2 表面自旋波頻率隨溫度的變化Fig.2 Change of frequency of surface spin waves with temperature
圖3為外磁場磁感應強度對表面自旋波頻率的影響,且Ds=0.05、D=1、τ=0.5.圖3a中Js=0.75,而圖3b中Js=2.由圖3可見,聲學和光學表面自旋波頻率均隨著外磁場磁感應強度的增加而增加.這是由于當磁場增強時,磁振子的能量增大.對比圖3a、b可知,光學表面自旋波隨著磁場加強而增大的幅度大于聲學表面自旋波.當外磁場磁感應強度處于0.15~0.5區域時,隨著波矢的變化,聲學和光學表面自旋波頻率的變化并不明顯.當外磁場磁感應強度較大時,波矢對聲學和光學表面自旋波頻率的影響較大.這是由于當外磁場強度較大時,能譜較高,因而波矢的改變對能譜的影響較大.

圖3 表面自旋波頻率隨外磁場的變化Fig.3 Change of frequency of surface spin waves with external magnetic field
圖4為表面交換耦合參數對表面自旋波頻率的影響,且B=0.15、Ds=0.05、D=1、τ=0.5、J=1.由圖4可見,聲學和光學表面自旋波頻率均隨著表面交換耦合參數的增加而增加.這是由于當表面交換耦合作用增強時,磁振子的能量增大.對比圖4a、b可知,表面交換耦合對聲學表面自旋波的影響較大.

圖4 表面自旋波頻率隨表面交換耦合的變化Fig.4 Change of frequency of surface spin waves with surface exchange coupling
圖5為表面各向異性參數對表面自旋波頻率的影響,且τ=0.5、B=0.15、J=1.圖5a、b中Js分別為0.75和2.由圖5可見,聲學和光學表面自旋波頻率均隨著表面各向異性參數的增加而增加.這是由于當表面各向異性參數增加時,磁振子能量增大.對比圖5a、b可知,表面各向異性對聲學表面自旋波頻率的影響較大.

圖5 表面自旋波頻率隨表面各向異性的變化Fig.5 Change of frequency of surface spin waves with surface anisotropy
圖6為在布里淵區Γ-M區域中薄膜厚度對聲學和光學表面自旋波頻率的影響.其中,B=0.15、J=1、τ=0.5.圖6a、b中Js分別為0.75和2.在布里淵區Γ-M區域中,當波矢等于或大于臨界波矢kc時,薄膜中出現表面自旋波.由圖6a可見,當層數大于5后,聲學表面自旋波頻率隨著薄膜厚度的增加而稍微增大.由圖6b可見,當層數大于5時,光學表面自旋波頻率隨著薄膜厚度的增加而稍微減小.在布里淵區Γ-M區域中,光學表面自旋波頻率大于聲學表面自旋波頻率.觀察圖6可以發現,不論是聲學還是光學表面自旋波,4層與其他層具有明顯區別,且本文結果與文獻[14]一致,即表面自旋波頻率實際上與5個原子層以上的薄膜厚度無關.

圖6 表面自旋波頻率隨薄膜厚度的變化Fig.6 Change of frequency of surface spin waves with thickness of thin film
本文中外磁場的增強使得表面自旋波的頻率增大,這與文獻[15-16]中有關共振頻率的結論相符.文獻[17]研究了磁性薄膜表面自旋波的色散關系,結果表明,表面自旋波頻率隨著表面交換耦合和表面各向異性的增強而增大,這與本文結果一致.文獻[18]研究了薄膜層數對共振頻率變化的影響,其結果亦與本文結論相符.
采用格林函數方法,討論了鐵磁性薄膜中溫度、外磁場、表面交換耦合、表面各向異性和薄膜厚度對表面自旋波頻率的影響.隨著溫度的增加,聲學和光學表面自旋波頻率減小,但在低溫區聲學表面自旋波的變化比光學表面自旋波的變化大,在高溫區則相反.當外磁場磁感應強度增加時,聲學和光學表面自旋波頻率均增加,但外磁場對光學表面自旋波的影響較大.隨著表面交換耦合作用的增加,聲學和光學表面自旋波頻率也隨之明顯升高,但表面交換耦合對聲學表面自旋波的影響較大.隨著表面各向異性的增大,聲學和光學表面自旋波頻率均增大,且表面各向異性對聲學表面自旋波的影響較大.聲學和光學表面自旋波頻率隨波矢的增大而增大.當薄膜原子層數大于5時,光學表面自旋波頻率隨著薄膜厚度的增加而稍微減小,而聲學表面自旋波頻率則隨薄膜厚度的增加而稍微增大.