垂直磁各向異性自旋閥結(jié)構(gòu)中磁性相圖

王日興， 吳彥

垂直磁各向異性自旋閥結(jié)構(gòu)中磁性相圖

王日興， 吳彥

（湖南文理學(xué)院 電氣與信息工程學(xué)院， 湖南 常德， 415000）

通過(guò)線(xiàn)性展開(kāi)包含自旋轉(zhuǎn)移矩項(xiàng)的Landau- Lifshitz-Gilbert（LLG）方程， 并且使用穩(wěn)定性分析的方法， 建立了垂直磁各向異性自旋閥結(jié)構(gòu)中的磁性相圖。研究發(fā)現(xiàn)： 通過(guò)調(diào)節(jié)外磁場(chǎng)強(qiáng)度的大小以及直流電流密度，可實(shí)現(xiàn)磁矩從穩(wěn)定態(tài)到進(jìn)動(dòng)態(tài)之間的轉(zhuǎn)化， 以及在不同穩(wěn)定態(tài)之間的翻轉(zhuǎn)。

自旋閥結(jié)構(gòu)； 垂直磁各向異性； 穩(wěn)定性分析； 自旋轉(zhuǎn)移矩

1996年， L Berger和J C Slonczewski在理論上預(yù)言了自旋閥結(jié)構(gòu)中的自旋轉(zhuǎn)移效應(yīng)［1-2］， 由于其在信息的存儲(chǔ)和處理以及微波激發(fā)中具有重要的應(yīng)用前景， 一直以來(lái)是學(xué)術(shù)界研究的熱點(diǎn)領(lǐng)域［3-12］。當(dāng)垂直于具有釘扎層/隔離層/自由層的自旋閥結(jié)構(gòu)施加電流時(shí)， 局域磁矩和自旋極化電流將產(chǎn)生相互作用，被釘扎層自旋極化的電流流過(guò)自由層時(shí)會(huì)給自由層的磁矩施加自旋轉(zhuǎn)移矩， 從而驅(qū)動(dòng)自由層磁矩產(chǎn)生進(jìn)動(dòng)， 甚至發(fā)生翻轉(zhuǎn)［3-5］。磁矩產(chǎn)生進(jìn)動(dòng)在高頻微波振蕩器方面具有重大的應(yīng)用前景［3-5］， 而磁矩的翻轉(zhuǎn)預(yù)計(jì)將會(huì)給磁性隨機(jī)存儲(chǔ)器的數(shù)據(jù)寫(xiě)入帶來(lái)一種全新的方式［6-7］。

在自旋閥結(jié)構(gòu)中的磁動(dòng)力學(xué)研究中， 穩(wěn)定性分析方法［8］是一種非常重要的研究方法。通過(guò)對(duì)自旋閥結(jié)構(gòu)進(jìn)行穩(wěn)定性分析， 可以獲得穩(wěn)定磁性狀態(tài)的失穩(wěn)條件［9-11］， 同時(shí)也可以建立磁性相圖， 得到不同磁性狀態(tài)之間轉(zhuǎn)換的臨界電流密度。與具有平行磁各向異性的自旋閥結(jié)構(gòu)相比較， 具有垂直磁各向異性的自旋閥有更多的優(yōu)點(diǎn)［12］， 例如： 更大的磁電阻率， 更好的熱穩(wěn)定性以及更快速的磁矩翻轉(zhuǎn)。

本文以釘扎層和自由層都具有垂直磁各向異性的自旋閥結(jié)構(gòu)為理論模型， 通過(guò)線(xiàn)性展開(kāi)包含自旋轉(zhuǎn)移矩項(xiàng)的 Landau-Lifshitz-Gilbert（LLG）方程， 同時(shí)使用穩(wěn)定性分析的方法， 建立了以垂直磁各向異性自旋閥結(jié)構(gòu)的磁性狀態(tài)相圖， 相圖以電流密度和外磁場(chǎng)大小為控制參數(shù)。發(fā)現(xiàn)通過(guò)調(diào)節(jié)電流密度和外磁場(chǎng)的大小， 不僅可以得到準(zhǔn)平行和準(zhǔn)反平行穩(wěn)定態(tài)、伸出膜面的進(jìn)動(dòng)態(tài)和雙穩(wěn)態(tài)， 而且還可以實(shí)現(xiàn)不同磁性狀態(tài)之間的轉(zhuǎn)換。此外， 通過(guò)對(duì)微分方程進(jìn)行數(shù)值求解， 給出了磁矩進(jìn)動(dòng)和磁矩翻轉(zhuǎn)隨時(shí)間的演化

1　模型和穩(wěn)定性分析

本文所采用的理論模型如圖1 所示。一個(gè)由3層膜構(gòu)成的自旋閥結(jié)構(gòu)， 最上面的釘扎層和最下面的自由層均為鐵磁層， 中間的隔離層為非磁金屬層。3層膜膜面處于x—z平面，釘扎層和自由層的磁各向異性沿y軸方向， 即垂直于膜面。當(dāng)沿y軸方向施加電流時(shí)， 電流流經(jīng)釘扎層后被極化為和釘扎層具有相同方向的自旋極化電流， 自旋極化電流流過(guò)自由層時(shí)對(duì)自由層磁矩產(chǎn)生自旋轉(zhuǎn)移矩。此時(shí)， 可以用一個(gè)包含自旋轉(zhuǎn)移矩項(xiàng)的無(wú)量綱的 LLG方程來(lái)描述自由層磁矩的磁動(dòng)力學(xué)行為：

圖1　理論模型和坐標(biāo)系

式中， m為自由層磁矩的單位矢量， α為吉爾伯阻尼常量，ye為沿y軸方向的單位向量。， 其中：，kH為各向異性場(chǎng)；

yH 為沿y軸方向的外磁場(chǎng)； J為電流密度， 當(dāng)電流從釘扎層流向自由層時(shí)， 電流密度J取正值， 反之J取負(fù)值。

其中： P為自旋極化率； e為基本電荷； ?為普朗克常量；sM 和d分別為自由層的飽和磁化強(qiáng)度和厚度。為重新定義的時(shí)間變量。

在球坐標(biāo)系中將方程（1）展開(kāi)， 得到關(guān)于θ和φ的微分方程

從方程（1）容易得出m=ey和-ey為系統(tǒng)的平衡點(diǎn)， 對(duì)應(yīng)于球坐標(biāo)系中的（π/2， 3π/2）。在微擾的作用下， 磁矩將在這2個(gè)平衡位置附近振動(dòng)。可以令將其代入方程組（2）并進(jìn)行線(xiàn)性化， 可以得到關(guān)于δθ和δφ的線(xiàn)性微分方程假設(shè)磁矩做周期性振蕩， 令， 其中Δθ和Δφ為振幅， 將其代入式（3）， 可以得到行列式

式（4）中， iμω=。展開(kāi)行列式， 可以得到關(guān)于μ的一元二次方程：

利用式（6）和（7）， 建立垂直磁各向異性自旋閥結(jié)構(gòu)的磁性相圖。相關(guān)的材料參數(shù)為［12］： 吉爾伯特阻尼參數(shù)0.1α=； 旋磁比； 各向異性場(chǎng)； 飽和磁矩Ms=650 Gs； 自旋極化率0.35P=； 自由層的厚度為3 nmd=。對(duì)于此材料，

2　磁性相圖

以式（6）和（7）為基礎(chǔ)， 在圖2中， 以電流密度J和外磁場(chǎng)yH 為控制參數(shù)， 給出了垂直磁各向異性自旋閥結(jié)構(gòu)的磁性狀態(tài)相圖。圖2中的實(shí)線(xiàn)和虛線(xiàn)分別為平衡點(diǎn)和的失穩(wěn)電流。當(dāng)電流密度時(shí)， 自由層磁矩將穩(wěn)定在y軸附近； 當(dāng) J＞J2時(shí)， 自由層磁矩穩(wěn)定在-y軸附近， 分別用準(zhǔn)平行穩(wěn)定態(tài)和準(zhǔn)反平行穩(wěn)定態(tài)來(lái)描述這2種磁性狀態(tài)。當(dāng) J1＜J＜J2時(shí)， 自由層磁矩失穩(wěn)， 此時(shí)磁矩表現(xiàn)為伸出膜面的進(jìn)動(dòng)態(tài)。在其它區(qū)域， 磁矩表現(xiàn)為雙穩(wěn)態(tài)， 即準(zhǔn)平行穩(wěn)定態(tài)和準(zhǔn)反平行穩(wěn)定態(tài)共存， 因此磁矩的磁性狀態(tài)與它的初始磁化方向有關(guān)。如果磁矩的初始磁化方向沿y軸時(shí)， 磁矩為準(zhǔn)平行穩(wěn)定態(tài)， 而當(dāng)自由層磁矩的初始方向沿y- 軸時(shí)， 磁矩為準(zhǔn)反平行穩(wěn)定態(tài)。

磁矩從穩(wěn)定態(tài)到進(jìn)動(dòng)態(tài)的轉(zhuǎn)化以及在2個(gè)穩(wěn)定態(tài)之間的翻轉(zhuǎn)可以通過(guò)調(diào)節(jié)電流密度和外磁場(chǎng)的大小來(lái)實(shí)現(xiàn)。如圖 2所示， 當(dāng)外磁場(chǎng)yH和電流密度J處于“進(jìn)動(dòng)態(tài)”區(qū)域時(shí)， 磁矩表現(xiàn)為伸出膜面的進(jìn)動(dòng)態(tài)， 且與磁矩的初始方向無(wú)關(guān)。當(dāng)外磁場(chǎng)yH和電流密度J位于“準(zhǔn)反平行穩(wěn)定態(tài)”時(shí)， 由于磁矩最終會(huì)穩(wěn)定在y-軸方向， 因此， 當(dāng)磁矩初始狀態(tài)沿y方向， 磁矩將從y軸方向翻轉(zhuǎn)到y(tǒng)-軸方向； 而當(dāng)外磁場(chǎng)Hy和電流密度J位于“準(zhǔn)平行穩(wěn)定態(tài)”時(shí)， 由于磁矩最終會(huì)穩(wěn)定在y軸方向， 因此， 當(dāng)磁矩初始狀態(tài)沿-y方向， 可以實(shí)現(xiàn)磁矩從-y方向到y(tǒng)方向的翻轉(zhuǎn)。

圖2　電流密度J和外磁場(chǎng)　yH 為控制參數(shù)的磁性相圖

圖3　進(jìn)動(dòng)態(tài)的演化軌跡

為了驗(yàn)證磁性相圖是否正確， 在圖2的不同區(qū)域， 取a、b、c 三點(diǎn)， 通過(guò)數(shù)值求解微分方程（2）， 分別給出了磁矩的演化軌跡。圖3為圖2中a點(diǎn)的磁矩隨時(shí)間的演化軌跡。由圖3可知， 不管磁矩的初始磁化方向沿y軸方向還是y- 軸方向，經(jīng)歷一段時(shí)間， 磁矩將繞著y軸或者-y軸方向旋轉(zhuǎn)， 即表現(xiàn)為伸出膜面的進(jìn)動(dòng)態(tài)。圖4（a）為圖2中b點(diǎn)的磁矩隨時(shí)間的演化軌跡， 磁矩的初始磁化方向沿y軸， 即mx=0， my= 1， mz=0。經(jīng)歷一段時(shí)間， 從圖 4可知， mx=0， my=-1， mz=0， 即磁矩最終穩(wěn)定在- y軸方向， 也就是說(shuō)實(shí)現(xiàn)了磁矩從y軸方向向-y軸方向的翻轉(zhuǎn)。圖 4（b）為圖 2中 c點(diǎn)的磁矩隨時(shí)間的演化軌跡， 自由層磁矩的初始磁化方向沿-y軸， 即mx=0， my=-1， mz=0。從圖4（b）可知， 經(jīng)歷一段時(shí)間， mx=0， my=1， mz=0， 即磁矩最后穩(wěn)定在y軸方向， 也就是說(shuō)實(shí)現(xiàn)了磁矩從-y軸方向向y軸方向的翻轉(zhuǎn)。自由層磁矩隨時(shí)間的演化軌跡圖說(shuō)明利用線(xiàn)性穩(wěn)定性分析方法得到的垂直磁各向異性自旋閥結(jié)構(gòu)的磁性相圖是正確的。

圖4　磁矩翻轉(zhuǎn)的演化軌跡

3　結(jié)論

本文以釘扎層和自由層都具有垂直磁各向異性的自旋閥為研究對(duì)象， 通過(guò)對(duì)包含自旋轉(zhuǎn)移矩項(xiàng)的LLG方程線(xiàn)性展開(kāi)， 并且使用穩(wěn)定性分析的方法， 建立了垂直磁各向異性自旋閥結(jié)構(gòu)的磁性相圖。相圖中包括穩(wěn)定態(tài)和進(jìn)動(dòng)態(tài)等各種不同的磁性狀態(tài)。通過(guò)調(diào)節(jié)外磁場(chǎng)和電流密度的大小， 可以實(shí)現(xiàn)磁矩從穩(wěn)定態(tài)到進(jìn)動(dòng)態(tài)的轉(zhuǎn)化以及磁矩在準(zhǔn)平行態(tài)和準(zhǔn)反平行態(tài)之間的翻轉(zhuǎn)。通過(guò)數(shù)值求解微分方程， 給出了不同磁性狀態(tài)的演化軌跡。這些研究將有助于提高基于自旋轉(zhuǎn)移矩的磁矩翻轉(zhuǎn)和微波激發(fā)的效率。

［1］ Berger L. Emission of spin waves by a magnetic multilayer traversed by a current ［J］. Physical Review B， 1996， 54（18）：9 353-9 358.

［2］ Slonczewski J C. Current-driven excitation of magnetic multilayers ［J］. Journal of Magnetism and Magnetic Materials，1996， 159： 1-7.

［3］ Kiselev S I， Sankey J C， Krivorotov I N， et al. Microwave oscillations of a nanomagnet driven by a spin-polarized current ［J］. Nature， 2003， 425（25）： 380-383.

［4］ 方彬， 曾中明. 自旋納米振蕩器［J］. 中國(guó)科學(xué)， 2014， 59（19）： 1 804-1 811.

［5］ 董浩， 陳培毅， 鄧寧. 對(duì)稱(chēng)磁多層結(jié)構(gòu)納米自旋傳輸矩微波振蕩器［J］. 微納電子技術(shù)， 2010， 47（1）： 1-5.

［6］ Myers E B， Ralph D C， Katine J A， et al. Current-Induced Switching of Domains in Magnetic Multilayer Devices ［J］. Science， 1999， 285（6）： 867-870.

［7］ 韓秀峰. 一種新型磁隨機(jī)存取存儲(chǔ)器原理型器件的設(shè)計(jì)與研制［J］. 科學(xué)通報(bào)， 2007， 52（10）： 122.

［8］ 劉秉正， 彭建華. 非線(xiàn)性動(dòng)力學(xué)［M］. 北京： 高等教育出版社， 2007： 31-34.

［9］ Bertotti G， Serpico C， Mayergoyz I D， et al. Magnetization switching and microwave oscillations in nanomagnets driven by Spin-Polarized currents ［J］. Phys Rev Lett， 2005， 94（12）： 127206（1-4）.

［10］ He P B， Wang R X， Li Z D， et al. Current-driven magnetization dynamics in magnetic trilayers with a tilted spin polarizer ［J］. Eur Phys J B， 2010， 73： 417-421.

［11］ Wang R X， He P B， Li Z D， et al. Phase diagram of magnetic multilayers with tilted dual spin torques ［J］. J Appl Phys，2011， 109（3）： 039905（1-7）.

［12］ Mangin S， Ravelosona D， Katine J A， et al. Current-inducedmagnetization reversal in nanopillars with perpendicular anisotropy ［J］. Nature Materials， 2006（5）： 210-215.

（責(zé)任編校： 劉曉霞）

Research on magnetic phase diagram in spin valve structure with perpendicular anisotropy

Wang Rixing， Wu Yan
（College of Electrical and Information Engineering， Hunan University of Arts and Science， Changde 415000， China）

By spreading the Landau-Lifshitz-Gilbert （LLG） equation including the spin transfer-torque term and using the stability analysis method， the magnetic phase diagram is obtained. By adjusting the dc current density and the magnitude of external magnetic field， the switching from stable states to precession states and the magnetization reversal between different stable states can be realized.

spin valve structure； perpendicular anisotropy； stability analysis； spin-transfer torque

O 482.5

1672-6146（2016）03-0026-04

10.3969/j.issn.1672-6146.2016.03.006

王日興， wangrixing1982@sina.com。

2016-06-02

國(guó)家自然科學(xué)基金（11347132）； 湖南省教育廳一般項(xiàng)目（14C0807）； 湖南文理學(xué)院校級(jí)科研項(xiàng)目（14YB02）。軌跡。