強(qiáng)進(jìn) 何開宙 劉東妮 盧啟海 韓根亮 宋玉哲 王向謙

(甘肅省傳感器與傳感技術(shù)重點(diǎn)實(shí)驗(yàn)室,甘肅省科學(xué)院傳感技術(shù)研究所,蘭州 730000)

具有手型和極性雙重屬性的磁渦旋結(jié)構(gòu),自發(fā)現(xiàn)以來(lái)就被視為有望成為下一代自旋電子學(xué)器件的可能載體之一.尤其近年來(lái)隨著Dzyaloshinskii—Moriya 相互作用的發(fā)現(xiàn),在具有中心反演對(duì)稱性破缺或強(qiáng)自旋軌道耦合的體系中,使得一般的面內(nèi)磁渦旋具有更加多樣的動(dòng)力學(xué)行為.本文通過(guò)微磁學(xué)模擬的方法,系統(tǒng)研究了在等邊三角形結(jié)構(gòu)中磁渦旋能夠穩(wěn)定存在的條件,并在此基礎(chǔ)上分別通過(guò)施加面內(nèi)和面外兩個(gè)方向的微波磁場(chǎng)來(lái)激勵(lì)其振蕩,其中除常見的面內(nèi)旋轉(zhuǎn)模式和面外呼吸模式外,還存在高頻微波磁場(chǎng)下的分裂模式,以及呼吸旋轉(zhuǎn)同步的自旋波模式.最后,通過(guò)改變體系中Dzyaloshinskii—Moriya 有效場(chǎng)的強(qiáng)度來(lái)改變整個(gè)三角形中的磁結(jié)構(gòu),進(jìn)而調(diào)控不同自旋波模式的本征頻率.本文結(jié)果對(duì)研究自旋波模式的多樣性具有一定的借鑒意義,并且為多類型的自旋波模式能夠在相關(guān)自旋電子學(xué)器件的研發(fā)提供更多選擇.

1 引言

磁渦旋作為納米級(jí)磁結(jié)構(gòu)的一種,具有體積小、易于集成、易調(diào)控、驅(qū)動(dòng)電流密度低、熱損耗低等優(yōu)勢(shì),由于其潛在的應(yīng)用價(jià)值和研究意義自發(fā)現(xiàn)起就受到越來(lái)越多的關(guān)注和青睞[1-5].磁渦旋的存在是由于體系內(nèi)交換能和靜磁能相互競(jìng)爭(zhēng)的結(jié)果[6,7].結(jié)構(gòu)上通常包含面內(nèi)磁矩和垂直于面內(nèi)的中心區(qū)域磁矩,因此通常具有手性與極性的雙重屬性,其中手性與面內(nèi)磁矩排布方式有關(guān),可分為順時(shí)針與逆時(shí)針兩種,而極性則是用+1,—1 分別代表磁渦旋核中心區(qū)域磁矩沿z軸向上和向下[8,9].通過(guò)驅(qū)動(dòng)磁渦旋極性翻轉(zhuǎn),可作為信息載體應(yīng)用于存儲(chǔ)器件中,具有讀寫速度快、易擦除和改寫等優(yōu)勢(shì)[10-12].而且有望應(yīng)用于下一代自旋電子學(xué)器件中,如基于磁渦旋所設(shè)計(jì)的自旋納米振蕩器,可實(shí)現(xiàn)持續(xù)輸出高頻微波信號(hào)[13-15].為進(jìn)一步增強(qiáng)磁渦旋的可應(yīng)用性,在具有對(duì)稱性破缺或者強(qiáng)自旋軌道耦合的體系中引入Dzyaloshinskii-Moriya 相互作用(DMI)[16-19],可通過(guò)改變磁渦旋結(jié)構(gòu)來(lái)調(diào)控其動(dòng)力學(xué)過(guò)程.在缺乏界面反演對(duì)稱性的異質(zhì)結(jié)構(gòu)體系中,DM 有效場(chǎng)起著迫使相鄰磁矩垂直排列的作用,所以DM 相互作用的存在可使得磁渦旋中面內(nèi)磁矩被拉到垂直方向上來(lái).本文通過(guò)改變體系中DM 有效場(chǎng)強(qiáng)度來(lái)改變?nèi)切沃写艤u旋結(jié)構(gòu),三角形納米磁結(jié)構(gòu)的尖角位置能產(chǎn)生強(qiáng)大的靜磁場(chǎng).由于形狀對(duì)稱性而產(chǎn)生的平面內(nèi)簡(jiǎn)并的旋轉(zhuǎn)對(duì)稱性消失,這能幫助提高磁力顯微鏡的對(duì)比度,從而有利于觀察到更加清晰的圖形,這有助于在基于三角形納米磁結(jié)構(gòu)所設(shè)計(jì)的磁存儲(chǔ)設(shè)備中實(shí)現(xiàn)更好的信號(hào)編碼[20,21],另外,三角形納米磁結(jié)構(gòu)在自旋波動(dòng)力學(xué)方面的研究也表現(xiàn)出獨(dú)特的磁特性[22].在得到穩(wěn)定的磁渦旋后,在不同方向上施加微波磁場(chǎng)來(lái)激勵(lì)其振蕩,通過(guò)快速傅里葉變換得到本征頻率后,進(jìn)一步探究其在不同本征頻率下的自旋波模式以及不同模式的變化過(guò)程,最后,改變體系中DM 有效場(chǎng)的強(qiáng)度,實(shí)現(xiàn)了調(diào)控不同自旋波模式下的本征頻率.這些結(jié)果為磁渦旋在自旋電子學(xué)相關(guān)領(lǐng)域的發(fā)展及應(yīng)用提供了一定的借鑒意義.

2 微磁學(xué)模擬

通過(guò)微磁學(xué)模擬的方法進(jìn)行計(jì)算,研究了等邊三角形結(jié)構(gòu)中磁渦旋的自旋波模式,其相關(guān)磁化動(dòng)力學(xué)研究則是通過(guò)Landau-Lifshitz-Gilbert 方程來(lái)表達(dá)[23]:

對(duì)于一般的磁渦旋來(lái)說(shuō),由于體系中交換場(chǎng)的作用,相鄰磁矩之間趨于平行排列,使得磁矩平躺于面內(nèi),在垂直于面內(nèi)方向是沒(méi)有分量的,然而在考慮DM 相互作用后磁矩會(huì)向著z方向發(fā)生偏轉(zhuǎn),這是因?yàn)镈M 有效場(chǎng)的作用是迫使相鄰磁矩垂直排列,所以在二者的競(jìng)爭(zhēng)之下,相鄰磁矩之間會(huì)處于面外向面內(nèi)過(guò)渡的一種狀態(tài),當(dāng)改變DMI 強(qiáng)度時(shí),就會(huì)出現(xiàn)如圖1 所示的不同磁結(jié)構(gòu).當(dāng)D=1.5 mJ/m2時(shí),磁矩出現(xiàn)沿z方向從向上到向下過(guò)渡,當(dāng)D=2.0 mJ/m2時(shí),逐漸出現(xiàn)成核跡象,只有當(dāng)DMI 常數(shù)大于 2.5 mJ/m2時(shí),才會(huì)在三角形中心附近形成渦旋結(jié)構(gòu).

圖1 不同DMI 有效場(chǎng)強(qiáng)度下的三角形結(jié)構(gòu)中磁渦旋基態(tài)Fig.1.Ground state of magnetic vortex in triangular structure under different DMI constant.

3 模擬結(jié)果與分析

首先,研究了磁渦旋分別在面內(nèi)(x)和面外(z)方向的本征頻率,即在兩個(gè)方向上分別施加一個(gè)sinc 形式的磁場(chǎng),如圖2(a)所示,其具體表達(dá)形式為hi(t)=h0sinc[2πf(t ?t0)]=h0sin[2πf(t?t0)]/[2πf(t ?t0)][25],h0為振幅,固定為10 mT,f為截止頻率,設(shè)為100 GHz.t0=1ns,表示該激勵(lì)磁場(chǎng)作用在1 ns 處,i代表x或者z方向.該sinc 函數(shù)形式磁場(chǎng)形式與激發(fā)三角形中磁渦旋的自旋波模式?jīng)]有關(guān)系,施加該函數(shù)形式微波磁場(chǎng)的目的是激勵(lì)磁矩振動(dòng),待磁矩弛豫穩(wěn)定后通過(guò)快速傅里葉變換得到動(dòng)態(tài)磁化率虛部隨頻率的變化情況.在探究x和z兩個(gè)方向的本征頻率時(shí),該動(dòng)力學(xué)過(guò)程總時(shí)長(zhǎng)皆為100 ns,在整個(gè)過(guò)程中,每10 ps采集一次數(shù)據(jù)(磁化分布mi和激勵(lì)場(chǎng)hi),并通過(guò)快速傅里葉變換從時(shí)域轉(zhuǎn)換到頻域,磁化率虛部Imχi峰值所在處所對(duì)應(yīng)頻率即為本征頻率,在該頻率下磁渦旋對(duì)微波具有很強(qiáng)的響應(yīng).如圖2 所示,選擇D=2.5 mJ/m2時(shí)的穩(wěn)定磁結(jié)構(gòu)為初始態(tài),圖2(b)和(c)分別為面內(nèi)和面外本征頻率分布情況.因?yàn)檠豿方向和y方向施加微波激勵(lì)磁場(chǎng)后能夠得到相同的本征頻率,所以以x方向?yàn)槔紤]面內(nèi)自旋波的激發(fā).當(dāng)所加激勵(lì)磁場(chǎng)沿面內(nèi)方向時(shí),為面內(nèi)自旋波模式,其本征頻率分別為0.37 GHz,0.74 GHz,4.03 GHz 和5.57 GHz,可以看出當(dāng)頻率為0.37 GHz 時(shí)對(duì)微波的響應(yīng)明顯大于其他頻率.當(dāng)激勵(lì)磁場(chǎng)沿面外方向時(shí),對(duì)應(yīng)面外自旋波模式的本征頻率分別為1.98 GHz,5.68 GHz 和9.67 GHz,可以看出峰值最強(qiáng)處頻率為1.98 GHz.

圖2 (a) sinc 形式的磁場(chǎng)隨時(shí)間變化情況;(b)和(c)分別為激勵(lì)磁場(chǎng)在x 和z 方向上時(shí),磁渦旋的動(dòng)態(tài)磁化率虛部Imχx 和Imχz Fig.2.(a) The relationship between time and sinc-type microwave field;(b) and (c) are the imaginary part of dynamic susceptibility Imχx and Imχz when the exciting magnetic field is applied along in-plane direction and out-of-plane direction,respectively.

在得到磁渦旋本征頻率的基礎(chǔ)上,進(jìn)一步研究了其在各個(gè)本征頻率下不同自旋波模式的類型.如圖3 所示,分別展示了面內(nèi)和面外兩個(gè)方向上在不同本征頻率下傅里葉變換后的振幅和相位,在各自本征頻率下,第1 行表示x,y,z方向上的振幅強(qiáng)度,第2 行表示3 個(gè)方向上的相位分布.對(duì)不同自旋波模式的提取主要在于觀察磁結(jié)構(gòu)在各自本征頻率的微波磁場(chǎng)下振幅及相位的變化情況,對(duì)于常見的面內(nèi)旋轉(zhuǎn)模式來(lái)說(shuō),主要表現(xiàn)為在磁平面上其相位能夠連續(xù)從π 變化到—π,而對(duì)于面外呼吸模式來(lái)說(shuō),主要表現(xiàn)為磁疇壁區(qū)域周期性地?cái)U(kuò)張和收縮,并且其振幅和相位是關(guān)于徑向?qū)ΨQ[26,27].本文對(duì)于面內(nèi)自旋波模式來(lái)說(shuō),當(dāng)f1=0.37 GHz 時(shí),由垂直于面內(nèi)方向觀察到共振區(qū)域主要集中在邊界處,同時(shí)結(jié)合x及y方向傅里葉空間振幅分布可知三角形磁結(jié)構(gòu)中大部分磁矩都對(duì)微波磁場(chǎng)產(chǎn)生了響應(yīng),所以在該頻率下對(duì)微波磁場(chǎng)的吸收功率最大,另外,相位連續(xù)從π 變化到了—π,為面內(nèi)逆時(shí)針旋轉(zhuǎn)模式.當(dāng)f2=0.74 GHz 時(shí),共振區(qū)域依然主要集中在邊界處,但是結(jié)合x及y方向振幅分布情況可知對(duì)微波磁場(chǎng)產(chǎn)生響應(yīng)的磁矩區(qū)域相對(duì)有所減少,所以在該頻率下的吸收峰值有所降低,而且此時(shí)出現(xiàn)了分疇旋轉(zhuǎn)的現(xiàn)象,空間相位出現(xiàn)了兩次從π 到—π的變化,為面內(nèi)逆時(shí)針旋轉(zhuǎn)模式.當(dāng)f3=4.03 GHz 和f4=5.57 GHz 時(shí),結(jié)合x,y,z方向上傅里葉空間振幅分布情況可知,共振區(qū)域明顯減小,所以在動(dòng)態(tài)磁化率虛部隨頻率的變化情況中吸收功率最低.而且分疇現(xiàn)象更加明顯,也就是說(shuō)此時(shí)高頻率微波磁場(chǎng)改變了三角形磁結(jié)構(gòu),使得中心處的磁渦旋出現(xiàn)了分裂,但依然可以從空間相位的變化看出是從π 變化到了—π,可定義為分裂的面內(nèi)逆時(shí)針旋轉(zhuǎn)模式.對(duì)于面外自旋波模式來(lái)說(shuō),當(dāng)f5=1.98 GHz 和f6=5.68 GHz 時(shí),共振區(qū)域的分布是關(guān)于三角形中心對(duì)稱的,邊界處振幅較強(qiáng),而中心區(qū)域?qū)ξ⒉ǖ捻憫?yīng)較弱,從空間相位的變化情況來(lái)看,可以發(fā)現(xiàn)三角形中心的磁渦旋分裂為3 個(gè)區(qū)域,而且同步逆時(shí)針旋轉(zhuǎn),空間相位都是從π 變化到了—π,同時(shí)可以發(fā)現(xiàn)空間相位相同處距離三角形中心的距離相等,表現(xiàn)出呼吸模式,綜上可以說(shuō)明在f5=1.98 GHz 和f6=5.68 GHz 時(shí),為面外呼吸逆時(shí)針旋轉(zhuǎn)模式.當(dāng)f7=9.67 GHz 時(shí),此時(shí)在三角形結(jié)構(gòu)在具有高本征頻率的微波磁場(chǎng)作用下發(fā)生了變化,磁疇分裂為多個(gè)區(qū)域,但中心區(qū)域的空間相位依然是從π 變化到了—π,可定義為分裂的面外逆時(shí)針旋轉(zhuǎn)模式.對(duì)比3 種面外自旋波模式,同時(shí)結(jié)合在此三個(gè)本征頻率下x,y,z方向上傅里葉空間振幅的分布情況,亮色的區(qū)域面積相差不大,表明對(duì)微波磁場(chǎng)產(chǎn)生共振的磁矩區(qū)域面積相接近,所以在面外自旋波激發(fā)過(guò)程中,這3 個(gè)本征頻率的吸收峰值相差不大,這與圖2 中的結(jié)果相符合[28,29].

圖3 磁渦旋在x, y, z 方向上不同自旋波模式的空間振幅及相位分布Fig.3.The spatial amplitude and phase distributions of the different spin wave mode of magnetic vortex in the x,y,and z directions.

為了進(jìn)一步驗(yàn)證上述7 種模式,分別在面內(nèi)和面外兩個(gè)方向上施加了不同本征頻率的正弦微波磁場(chǎng)來(lái)激勵(lì)整個(gè)三角形結(jié)構(gòu),該微波磁場(chǎng)的形式為:Bi(t)=B0sin(2πft),i代表x或者z方向,f為本征頻率.圖4 記錄了不同時(shí)刻下相對(duì)于初始時(shí)刻凈余磁矩的變化情況.采用δmz(t)=mz(t)?mz(0)描述凈余磁矩,其中mz(0)為初始時(shí)刻靜態(tài)三角形結(jié)構(gòu)的空間磁矩.對(duì)于面內(nèi)自旋波模式來(lái)說(shuō),當(dāng)f1=0.37 GHz 時(shí),磁矩隨時(shí)間演化呈現(xiàn)出明顯的逆時(shí)針旋轉(zhuǎn),為面內(nèi)逆時(shí)針旋轉(zhuǎn)模式(Mode 1).當(dāng)f2=0.74 GHz 時(shí),渦旋在逆時(shí)針旋轉(zhuǎn)的同時(shí)出現(xiàn)了分裂,整體依然是面內(nèi)逆時(shí)針旋轉(zhuǎn)模式(Mode 2).當(dāng)f3=4.03 GHz 和f4=5.57 GHz 時(shí),可以看到整個(gè)三角形結(jié)構(gòu)中的磁矩分布比較混亂,分裂為多個(gè)區(qū)域,磁矩整體在做逆時(shí)針旋轉(zhuǎn),為分裂的面內(nèi)逆時(shí)針旋轉(zhuǎn)模式(Mode 3 和Mode 4).對(duì)于面外自旋波模式來(lái)說(shuō),當(dāng)f5=1.98 GHz 和f6=5.68 GHz時(shí),三角形中心的渦旋區(qū)域分裂為三部分,呈現(xiàn)出明顯的周期性變化,該區(qū)域在不斷地收縮和擴(kuò)張的同時(shí)也在做逆時(shí)針旋轉(zhuǎn),為面外呼吸逆時(shí)針旋轉(zhuǎn)模式(Mode 5 和Mode 6).當(dāng)f7=9.67 GHz 時(shí),三角形中心的渦旋分裂為多個(gè)區(qū)域,但是磁矩整體還在做逆時(shí)針旋轉(zhuǎn),為分裂的面外逆時(shí)針旋轉(zhuǎn)模式(Mode 7).

圖4 磁渦旋在對(duì)應(yīng)本征頻率的微波磁場(chǎng)下不同時(shí)刻的凈余磁矩(δmz)Fig.4.The net magnetization (δmz) of the magnetic vortex at different time under the microwave magnetic field,which corresponding to the eigenfrequencies.

最后,在保證中心處渦旋穩(wěn)定存在的前提下,改變了體系中DMI的強(qiáng)度,進(jìn)而調(diào)控三角形磁渦旋結(jié)構(gòu)中自旋波模式的本征頻率.DM 有效場(chǎng)起著改變?nèi)切沃行拇艤u旋結(jié)構(gòu)的作用,如圖5(a)所示,當(dāng)DMI 數(shù)值增大時(shí),整個(gè)三角形結(jié)構(gòu)中磁渦旋所占區(qū)域的比例有所降低,這主要是由于DMI有效場(chǎng)的作用是將磁矩拉到垂直于面內(nèi)的方向上來(lái),所以隨著DMI 數(shù)值的增大面內(nèi)方向上的磁矩逐漸向z方向過(guò)渡,同時(shí)受到邊界的限制作用,所以渦旋所在區(qū)域的面積逐漸減小.當(dāng)三角形磁結(jié)構(gòu)隨著DMI 數(shù)值的變化而改變時(shí),不同磁結(jié)構(gòu)對(duì)微波的響應(yīng)情況也是不同的,從而影響不同自旋波模式的本振頻率,如圖5(b)和(c)所示,分別為面內(nèi)和面外各自旋波模式隨DMI 有效場(chǎng)強(qiáng)度變化時(shí)所呈現(xiàn)的功率譜密度相圖.由圖3 結(jié)果可知,三角形結(jié)構(gòu)中的磁渦旋具有4 種面內(nèi)自旋波模式和3 種面外自旋波模式,對(duì)于Mode 1,只有三角形中心處小部分磁矩對(duì)微波磁場(chǎng)沒(méi)有產(chǎn)生響應(yīng),其余大部分磁矩都產(chǎn)生了不同程度的振蕩,隨著DMI 數(shù)值的增大,三角形中心磁渦旋區(qū)域減小,因此對(duì)微波產(chǎn)生振蕩的磁矩區(qū)域增大,因此該模式下的本征頻率隨DMI 數(shù)值增大而增大.對(duì)于Mode 2,三角形中心磁渦旋在微波磁場(chǎng)的作用下出現(xiàn)了分疇旋轉(zhuǎn)的現(xiàn)象,而疇壁的產(chǎn)生會(huì)使得該部分區(qū)域磁矩整體對(duì)微波的響應(yīng)減弱,并且隨著DMI的增大,疇壁區(qū)域在z方向的分量越大,從而導(dǎo)致磁矩整體對(duì)微波磁場(chǎng)產(chǎn)生響應(yīng)的減弱效果越明顯,致使本征頻率降低.對(duì)于Mode 3 和Mode 4,對(duì)微波磁場(chǎng)產(chǎn)生共振的磁矩主要集中在三角形頂點(diǎn)處,而磁渦旋結(jié)構(gòu)在隨著DMI 變化的過(guò)程中,對(duì)頂點(diǎn)處磁結(jié)構(gòu)影響不大,所以這兩種模式下的本征頻率能夠相對(duì)維持穩(wěn)定.對(duì)于Mode 5 和Mode 6,出現(xiàn)了更加明顯的分疇現(xiàn)象,理同Mode 2,本征頻率出現(xiàn)下降.對(duì)于Mode 7,共振區(qū)域出現(xiàn)在磁渦旋附近,所以在隨著DMI 增大致使磁渦旋面積減小時(shí),共振區(qū)域面積增大,本征頻率隨之提高.

圖5 (a) 磁渦旋所在區(qū)域占三角形結(jié)構(gòu)整體面積之比隨DMI 系數(shù)變化情況;(b)和(c)分別為激勵(lì)磁場(chǎng)施加在面內(nèi)和面外方向上時(shí),磁渦旋的DMI 系數(shù)與共振頻率的關(guān)系相圖Fig.5.(a) The relationship between DMI constant and the ratio of the area where the magnetic vortex is located to the overall area of the triangular structure.Phase diagram of magnetic vortex as a function of DMI constant and resonance frequency,where (b)and (c) represent the excitation field applied in in-plane and out-of-plane directions,respectively.

4 結(jié)論

本文通過(guò)微磁學(xué)模擬的方法,對(duì)于穩(wěn)定的等邊三角形磁結(jié)構(gòu),分別在面內(nèi)和面外方向上施加微波磁場(chǎng)后,可激發(fā)出多種不同的自旋波模式.在面內(nèi)微波磁場(chǎng)的作用下,除常見的面內(nèi)旋轉(zhuǎn)模式外,在高本征頻率的微波磁場(chǎng)作用下,該體系中還存在分裂的面內(nèi)旋轉(zhuǎn)模式.在面外微波磁場(chǎng)的作用下,存在著面外自旋波呼吸旋轉(zhuǎn)模式,即發(fā)生周期性呼吸的同時(shí)也在做旋轉(zhuǎn),此外,面外高頻微波磁場(chǎng)依然會(huì)使得整個(gè)磁結(jié)構(gòu)發(fā)生分裂,出現(xiàn)分裂的面外旋轉(zhuǎn)模式.最后,通過(guò)改變體系中DMI 強(qiáng)度來(lái)改變整個(gè)三角形的磁結(jié)構(gòu),從而起到調(diào)控不同自旋波模式的本征頻率的作用.本文結(jié)果為利用自旋波共振來(lái)操作磁渦旋提供了可能,并且也為磁渦旋在自旋電子學(xué)中的應(yīng)用奠定了一定的基礎(chǔ).