蒙特卡羅模擬人工自旋冰的退磁過(guò)程

虞麗菊 孟祥雨 李俊琴 曹杰峰 王 勇 敬 超 吳衍青 邰仁忠

?

蒙特卡羅模擬人工自旋冰的退磁過(guò)程

虞麗菊1孟祥雨2李俊琴2曹杰峰2王 勇2敬 超1吳衍青2邰仁忠2

1（上海大學(xué) 物理系 上海 200444）2（中國(guó)科學(xué)院上海應(yīng)用物理研究所 張江園區(qū) 上海 201204）

近年來(lái)，人工自旋冰(Artificial spin ice, ASI)因具有潛在的應(yīng)用價(jià)值以及理論研究?jī)r(jià)值而倍受關(guān)注，成為磁學(xué)領(lǐng)域的研究熱點(diǎn)之一。以往的研究集中在Square和Kagome兩種晶格結(jié)構(gòu)，由于能量勢(shì)壘和阻挫的大量存在，實(shí)驗(yàn)上難以通過(guò)退磁等方法達(dá)到其內(nèi)部相互作用決定的能量基態(tài)。本文基于Kagome晶格結(jié)構(gòu)提出了一種新的人工自旋冰結(jié)構(gòu)，并基于蒙特卡羅模擬了三種不同的退磁方法對(duì)體系能量的影響。結(jié)果表明，退磁過(guò)程中旋轉(zhuǎn)磁場(chǎng)方向可更有效地降低體系能量，并誘導(dǎo)多種磁結(jié)構(gòu)的出現(xiàn)。

人工自旋冰，蒙特卡羅模擬，退磁

自旋冰是Harris等[1]在研究鈦酸鹽Ho2Ti2O7時(shí)，發(fā)現(xiàn)的一種自旋阻挫新象。該體系中晶格結(jié)構(gòu)由四面體基本單元構(gòu)成，每個(gè)Ho原子周圍存在4個(gè)Ti原子，其中兩個(gè)自旋指向四面體內(nèi)部，另外兩個(gè)指向外部，遵循“冰法則”排列[2]。這種鈦酸鹽存在絕對(duì)零度剩余熵和阻挫現(xiàn)象，并且其剩熵幅值與Pauling預(yù)言的冰氫剩熵幅值大小一致[3]，因此稱之為自旋冰。

阻挫指體系中等價(jià)單元之間存在相互競(jìng)爭(zhēng)，但每一方都無(wú)法取得支配地位的現(xiàn)象，廣泛存在于自然界及凝聚態(tài)體系中。研究阻挫現(xiàn)象對(duì)于計(jì)算機(jī)錯(cuò)誤校正碼[4]、蛋白質(zhì)動(dòng)力學(xué)[5]以及生物神經(jīng)網(wǎng)絡(luò)[6]的模型研究具有重要幫助。磁學(xué)中主要研究幾何阻挫[7]，如圖1所示，幾何晶格結(jié)構(gòu)和近鄰交換作用不相容，位于格點(diǎn)的磁矩?zé)o論如何排列，都不能滿足所有相鄰磁矩間相互作用能極小，形成阻挫。

圖1 幾何阻挫結(jié)構(gòu)示意圖Fig.1 Structure charts of geometrical frustration.

研究表明，阻挫的存在導(dǎo)致了自旋冰體系基態(tài)高度簡(jiǎn)并，是凝聚態(tài)領(lǐng)域中十分重要的研究體系。然而，鈦酸鹽等天然自旋冰存在樣品制備困難、表征平臺(tái)少以及分辨能力有限等缺點(diǎn)，導(dǎo)致相關(guān)的研究受阻。2006年，Wang等[8]采用電子束刻蝕的方法，制備了世界上首個(gè)人工自旋冰體系。該課題組在硅片基底上沉積單個(gè)尺寸為220 nm×80 nm× 25nm的坡莫合金，形成準(zhǔn)二維Square晶格陣列，并使用原子力顯微鏡(Atomic Force Microscope, AFM)和磁力顯微鏡(Magnetic Force Microscope, MFM)觀測(cè)了陣列的磁結(jié)構(gòu)成像。磁體為沿矩形長(zhǎng)軸方向排列的單疇?wèi)B(tài)，每4個(gè)磁體形成一個(gè)阻挫作用點(diǎn)。圖2為Square和Kagome兩種人工自旋冰體系的晶格示意圖。

圖2 兩種人工自旋冰結(jié)構(gòu)示意圖(a) Square晶格，(b) Kagome晶格Fig.2 Two Lattice structures of ASI.(a) Square lattice, (b) Kagome lattice

該項(xiàng)工作實(shí)現(xiàn)了自旋冰體系的人工制備和成像觀測(cè)，很大程度上降低了阻挫現(xiàn)象的研究難度。此外，這種有序的周期結(jié)構(gòu)可作為多體相互作用系統(tǒng)的統(tǒng)計(jì)力學(xué)研究平臺(tái)[9]，其內(nèi)部單元之間存在很強(qiáng)的作用力，磁耦合能和阻挫的相關(guān)研究還將有望給圖案化磁記錄技術(shù)提供參考。

實(shí)現(xiàn)以上研究和應(yīng)用首要解決的問(wèn)題就是尋找人工自旋冰體系的基態(tài)。理論認(rèn)為，體系總是在能量最小時(shí)處于平衡穩(wěn)定的狀態(tài)，而總能量是由內(nèi)部單元耦合相互作用決定的。由于阻挫和勢(shì)壘的大量存在，制備完成的樣品往往處于無(wú)序的能量激發(fā)態(tài)。使用傳統(tǒng)的退火處理效果不甚理想，原因在于自旋冰體系近鄰磁矩間的耦合作用能高達(dá)105K[10]，要在室溫附近對(duì)其進(jìn)行退火處理難度很大。Mrogan等[11]通過(guò)在樣品生長(zhǎng)過(guò)程中退火誘導(dǎo)Square自旋冰體系達(dá)到了有序基態(tài)。不少工作對(duì)自旋冰體系的退磁處理進(jìn)行了研究，Wang等[12]認(rèn)為如果在變換磁場(chǎng)的同時(shí)旋轉(zhuǎn)樣品，可使Square體系較大范圍接近能量基態(tài)。Mengotti等[13]使用類似的方法研究了Kagome自旋冰的退磁，結(jié)果表明，對(duì)單個(gè)六邊形進(jìn)行退磁處理，較易出現(xiàn)自旋首尾相接的能量基態(tài)，但隨著六邊形的數(shù)目增大，近鄰磁體形成阻挫，很難再出現(xiàn)這種基態(tài)。

圖3 渦旋人工自旋冰結(jié)構(gòu)示意圖(a) 單元結(jié)構(gòu)，中間形成六邊形，另外三個(gè)磁體分別沿[10]，[?11]，[?1?1]三個(gè)方向，(b) 每三個(gè)渦旋單元交匯處形成一個(gè)阻挫作用點(diǎn)，(c) 晶格示意圖Fig.3 Structure charts of Vortex ASI.(a) Unit structure, there is a hexagon in the middle, the other three along [10], [?11], [?1?1] respectively,(b) Every three units forming a frustrated Vortex,(c) Lattice structure

為降低阻挫對(duì)體系的干擾，我們基于Kagome晶格設(shè)計(jì)了一種渦旋自旋冰結(jié)構(gòu)，并采用蒙特卡羅模擬方法來(lái)研究三種不同的退磁手段對(duì)體系的退磁效果。如圖3所示，每個(gè)渦旋結(jié)構(gòu)單元由9個(gè)矩形坡莫合金磁體組成，磁體尺寸為220 nm×80nm× 25nm。單元中間6個(gè)磁體形成一個(gè)六邊形結(jié)構(gòu)，而在單元交接處，每3個(gè)磁體形成一個(gè)阻挫作用點(diǎn)。我們模擬了包括21×21個(gè)與41×41個(gè)渦旋單元結(jié)構(gòu)的體系，結(jié)果基本相同，因此以下以21×21為例進(jìn)行分析。

由于矩形磁體沿長(zhǎng)軸方向形成單疇?wèi)B(tài)，有兩種取向，可近似為伊辛自旋。對(duì)于渦旋單元結(jié)構(gòu)中心的六邊形，可以形成26共64種自旋組態(tài)。然而近鄰磁體之間存在耦合作用，當(dāng)渦旋結(jié)構(gòu)處于能量基態(tài)時(shí)，將形成自旋首尾相接的兩種磁組態(tài)：自旋順時(shí)針(Type-CW)和自旋逆時(shí)針(Type-CCW)，是兩種簡(jiǎn)并的能量基態(tài)。對(duì)于阻挫作用點(diǎn)，其基態(tài)自旋結(jié)構(gòu)為“兩進(jìn)一出”或“兩出一進(jìn)”，共有6種簡(jiǎn)并態(tài)。激發(fā)態(tài)呈現(xiàn)為自旋三個(gè)朝內(nèi)(Type-In)或者三個(gè)朝外(Type-Out)，是兩種能量簡(jiǎn)并態(tài)(圖4)。體系處于理想基態(tài)時(shí)(圖6(a))，六邊形形成Type-CW與Type-CCW交替排列，在阻挫作用點(diǎn)，則是“兩進(jìn)一出”和“兩出一進(jìn)”交替排列，而不會(huì)出現(xiàn)Type-In或Type-Out的激發(fā)態(tài)結(jié)構(gòu)。

圖4 渦旋人工自旋冰的4種典型結(jié)構(gòu)Fig.4 Four typical structures of Vortex ASI.

1 蒙特卡羅模擬

本文基于蒙特卡羅模擬三種不同的退磁方法對(duì)體系的退磁效果。其主要思想是：假設(shè)整個(gè)體系共有個(gè)自旋構(gòu)成一個(gè)組態(tài)，隨機(jī)分配給體系一個(gè)初始態(tài)，假設(shè)體系中的部分自旋發(fā)生翻轉(zhuǎn)而產(chǎn)生新的組態(tài)，通過(guò)近鄰磁耦合作用計(jì)算能量變化Δ及玻爾茲曼權(quán)重指數(shù)：0=exp(?Δ/B)，決定是否翻轉(zhuǎn)，其中：B為玻爾茲曼常數(shù)，為溫度，模擬中取為300 K；不斷重復(fù)上述過(guò)程，直到體系能量趨于穩(wěn)定。對(duì)于近鄰磁體之間的相互作用，本文選用磁偶極子模型[10]。該模型將每個(gè)磁體看成理想的點(diǎn)磁矩，磁矩之間的作用能大小為：

式中：0為玻爾茲曼常數(shù)；和指相互作用的兩個(gè)磁矩，其值由坡莫合金的磁化強(qiáng)度=800×103A?m?1計(jì)算得到；代表兩個(gè)磁矩中心的距離矢量。然而，當(dāng)近鄰磁矩間距較小時(shí)，磁體的體積不可忽略。本文改進(jìn)了磁偶極子模型，將每個(gè)220 nm× 80nm×25 nm的坡莫合金磁體細(xì)分成許多個(gè)大小為22 nm×10 nm×5 nm的單元，總的近鄰磁體之間的耦合能采用求和的方式得到。施加外磁場(chǎng)時(shí)，還需考慮體系的塞曼能：

由于兩個(gè)磁體之間的耦合能會(huì)隨著它們之間的距離增大而快速減小。對(duì)于一個(gè)給定的磁體，當(dāng)另一個(gè)磁體與它之間的耦合能小于來(lái)自最近鄰磁體耦合能的3%時(shí)，忽略不計(jì)。對(duì)于單個(gè)矩形磁體，由于形狀各向異性能的影響，退磁過(guò)程中需要克服一定的能量勢(shì)壘才能發(fā)生翻轉(zhuǎn)。根據(jù)文獻(xiàn)[14]，對(duì)于尺寸為220 nm×80 nm×25 nm的坡莫合金，形狀各向異性能大小為：0=3.17′10?18J（約為2.3′105K）?？紤]漲落的存在，計(jì)算中對(duì)每個(gè)磁體的勢(shì)壘進(jìn)行隨機(jī)分配：0j0，其中：為無(wú)量綱隨機(jī)分布，期望值=1，標(biāo)準(zhǔn)差選取了0.1?0.3之間的一系列參數(shù)進(jìn)行試驗(yàn)，對(duì)結(jié)果的影響很小，最終選取=0.25。

本工作中，基于能量最小原理的蒙特卡羅模擬計(jì)算包含以下具體步驟：1) 給體系隨機(jī)分配一個(gè)初始狀態(tài)；2) 隨機(jī)假設(shè)每個(gè)單元中的個(gè)磁體進(jìn)行翻轉(zhuǎn)，其中服從概率分布：P=(1/2)，(=1,2,…,9)；3) 計(jì)算翻轉(zhuǎn)所需的能量變化：Δ=E?E+Σ0，其中：E和E分別為翻轉(zhuǎn)后和翻轉(zhuǎn)前的體系總能量。如果Δ≤0，表明翻轉(zhuǎn)后體系能量降低，發(fā)生翻轉(zhuǎn)，回到第2步；4) 如果Δ>0，表明翻轉(zhuǎn)后體系能量將增加，計(jì)算玻爾茲曼權(quán)重指數(shù)：0=exp(?Δ/B)，產(chǎn)生一個(gè)隨機(jī)數(shù)，∈(0,1)，如果0>，發(fā)生翻轉(zhuǎn)，回到第2步；5) 如果0≤，不發(fā)生翻轉(zhuǎn)，回到第2步；6) 多次重復(fù)以上過(guò)程，直到系統(tǒng)達(dá)到穩(wěn)定的能量狀態(tài)。

本工作中我們一共采用三種退磁方法。方法1(圖5(a))：首先沿[10]方向(圖3)施加0.15 T的初始飽和磁場(chǎng)，等待一定的時(shí)間使系統(tǒng)達(dá)到平衡，然后快速變化到反方向，每次反向磁場(chǎng)的大小都減少一定幅度Δ，其值等于初始飽和磁場(chǎng)除以循環(huán)總數(shù)，待退磁完成所有的循環(huán)，磁場(chǎng)剛好降到零。方法2(圖5(b))：在每個(gè)循環(huán)中，磁場(chǎng)大小按一定的步長(zhǎng)變化，達(dá)到最大值時(shí)無(wú)需等待系統(tǒng)平衡，直接進(jìn)行下一個(gè)循環(huán)。方法3(圖5(c))：磁場(chǎng)的角度沿順時(shí)針（或逆時(shí)針）方向勻速變化，同時(shí)磁場(chǎng)的大小按一定的步長(zhǎng)均勻減小，直至磁場(chǎng)為零。

圖5 三種退磁方法示意圖(a) 方法1，(b) 方法2，(c) 方法3Fig.5 Demagnetization scheme.(a) Method 1, (b) Method 2, (c) Method 3

為尋找最優(yōu)的參數(shù)，我們對(duì)每種方法采取了一系列的循環(huán)次數(shù)和步數(shù)進(jìn)行模擬，并以退磁后體系的總能量為指標(biāo)，比較而得較優(yōu)參數(shù)。對(duì)于方法1，選取了50000、100000、150000、200000這4種循環(huán)次數(shù)，體系的退磁總能量先降低，然后趨于穩(wěn)定，磁場(chǎng)的大小快速變化，步數(shù)為1；方法2的循環(huán)次數(shù)取1000，步數(shù)分別取100?400，間隔為100，隨著步數(shù)增大，退磁總能量并沒(méi)有降低，取100即可；而對(duì)于方法3，循環(huán)次數(shù)取500和1000，步數(shù)分別取180?720，間隔為180，結(jié)果表明當(dāng)循環(huán)次數(shù)取1000，步數(shù)取540時(shí)，退磁總能量已降到最低，更精細(xì)化的參數(shù)對(duì)結(jié)果幾乎沒(méi)有影響。

2 結(jié)果與討論

首先，選取晶格常數(shù)=500 nm的體系進(jìn)行模擬計(jì)算，得到三種方法退磁后體系的自旋組態(tài)(圖6)。圖6(a)為理想的體系基態(tài)，中間六邊形形成Type-CW與Type-CCW交替排列的磁結(jié)構(gòu)，兩種結(jié)構(gòu)出現(xiàn)的頻率都為50%，在阻挫作用點(diǎn)，則是“兩進(jìn)一出”和“兩出一進(jìn)”交替排列。與有序基態(tài)相比，方法3更容易形成Type-CW和Type-CCW兩種基態(tài)結(jié)構(gòu)，方法2次之，由于近鄰單元之間存在復(fù)雜的相互作用，Type-CCW和Type-CW并不是交替排列的，而有可能同一類型連續(xù)出現(xiàn)，這是由于近鄰磁矩之間存在較強(qiáng)的耦合作用能，引起磁體連續(xù)翻轉(zhuǎn)造成的。圖7對(duì)三種方法退磁后體系中4種結(jié)構(gòu)的出現(xiàn)頻率進(jìn)行了統(tǒng)計(jì)，對(duì)應(yīng)圖6(b)?(d)。在理想情況下的最低能量狀態(tài)(圖6(a))，每一個(gè)六角結(jié)構(gòu)均為CW或者CCW，CW和CCW交替出現(xiàn)，各自占六角結(jié)構(gòu)總數(shù)的50%。在實(shí)際的退磁過(guò)程中，體系中磁體的翻轉(zhuǎn)是一個(gè)漲落過(guò)程，不可能達(dá)到最低能量狀態(tài)，有些六角結(jié)構(gòu)既不是CW也不是CCW，而是一種無(wú)序狀態(tài)。這種情況在方法1中尤其明顯，只有約10%的六角結(jié)構(gòu)是CW和CCW；方法2比方法1好一點(diǎn)，CW和CCW在所有六角結(jié)構(gòu)中的占比達(dá)到了20%左右；方法3最理想，達(dá)到了40%左右，體系的總能量也最接近基態(tài)，同時(shí)由于大量Type-CW和Type-CCW結(jié)構(gòu)的出現(xiàn)，誘導(dǎo)近鄰的阻挫作用點(diǎn)出現(xiàn)Type-In以及Type-Out的能量激發(fā)態(tài)，在所有三角結(jié)構(gòu)中的占比接近10%。比較三種方法退磁后的體系總能量(圖8)，方法3仍最大程度接近體系的基態(tài)能量1.947×10?15J。

圖6 蒙特卡羅計(jì)算體系的基態(tài)自旋排布，以及三種方法退磁后的自旋排布，對(duì)應(yīng)圖4中的4種結(jié)構(gòu)類型(a) 理想基態(tài)，(b) 方法1，(c) 方法2，(d) 方法3Fig.6 Structure chart for ground state and demagnetization results of three methods corresponding to Fig.4.(a) Ground state, (b) Method 1, (c) Method 2, (d) Method 3

另外在圖7中可以看出，在三種退磁方法中，Type-CW和Type-CCW、Type-In和Type-Out的數(shù)目并不嚴(yán)格相等。這是由于磁體的翻轉(zhuǎn)具有隨機(jī)性，進(jìn)行特定的退磁過(guò)程后，會(huì)產(chǎn)生對(duì)稱性破缺，各種結(jié)構(gòu)的數(shù)目不相等，體系不是處于最低能量狀態(tài)。并且由于每個(gè)磁體的翻轉(zhuǎn)存在勢(shì)壘，當(dāng)熱漲落的能量不足以克服勢(shì)壘時(shí)，整個(gè)體系會(huì)達(dá)到平衡狀態(tài)，這種平衡狀態(tài)不具有嚴(yán)格的對(duì)稱性。

圖7 三種方法退磁后4種結(jié)構(gòu)類型的統(tǒng)計(jì)頻率Fig.7 Statistical frequency of four types by three methods.

隨后調(diào)整體系的晶格常數(shù)為=400 nm和=600 nm，同樣模擬了三種方法的退磁效果。對(duì)于這兩種晶格參數(shù)的體系，方法3仍然是最理想的退磁方法，能將體系總能量降得最低(表1)。統(tǒng)一采用方法3退磁后，對(duì)等于400 nm、500 nm、600 nm的體系進(jìn)行比較，4種結(jié)構(gòu)類型的分布如圖9所示。當(dāng)=400 nm時(shí)，近鄰磁體之間距離更近，相互作用更強(qiáng)，退磁后容易形成能量較低的兩種結(jié)構(gòu)類型Type-CW和Type-CCW，頻率在40%左右，幾乎不會(huì)出現(xiàn)能量激發(fā)態(tài)結(jié)構(gòu)。隨著晶格常數(shù)增大，Type-CW和Type-CCW的出現(xiàn)頻率減小，Type-In和Type-Out兩種激發(fā)態(tài)結(jié)構(gòu)的出現(xiàn)頻率則在增大，=600 nm時(shí)，分別達(dá)到12.5%和12.8%。上述結(jié)果表明，晶格常數(shù)可有效地調(diào)控內(nèi)部相互作用，后者是驅(qū)動(dòng)體系抵達(dá)有序能量基態(tài)的主要因素，而退磁可以部分消除能量勢(shì)壘的影響，某種程度上誘導(dǎo)自旋實(shí)現(xiàn)“自由翻轉(zhuǎn)”，使體系更容易接近能量基態(tài)。

圖8 三種方法退磁后體系的總能量Fig.8 System energy after demagnetization by three methods.

表1 不同晶格常數(shù)的體系退磁后的總能量
Table 1 System energy of different lattices after demagnetization.

晶格常數(shù)

/ nm

方法1

Method 1 / 10J

方法2

Method 2 / 10J

方法3

Method 3 / 10J

最低能量

Min energy / 10J

400

?4.040

?4.143

?4.159

?4.227

500

?1.799

?1.811

?1.818

?1.947

600

?9.842

?9.761

?9.890

?1.094

圖9 方法3退磁后4種結(jié)構(gòu)的頻率分布Fig.9 Statistical frequency by method 3.

3 結(jié)語(yǔ)

本文基于蒙特卡羅方法模擬了三種方法對(duì)人工自旋冰體系的退磁效果，通過(guò)調(diào)整外磁場(chǎng)的循環(huán)次數(shù)、步數(shù)以及場(chǎng)方向，確定旋轉(zhuǎn)磁場(chǎng)退磁為最佳退磁方案。選取循環(huán)次數(shù)為1000，步數(shù)為540，磁場(chǎng)沿面內(nèi)旋轉(zhuǎn)，當(dāng)一個(gè)循環(huán)完成時(shí)回到原點(diǎn)。作為消除磁性材料勢(shì)壘的常用方法，退磁可使人工自旋冰體系中的自旋磁矩克服能量障礙，實(shí)現(xiàn)一定程度的“自由翻轉(zhuǎn)”，最終通過(guò)近鄰磁矩間的相互作用決定體系的自旋組態(tài)。無(wú)論是作為多體相互作用系統(tǒng)的統(tǒng)計(jì)力學(xué)平臺(tái)，還是圖案化磁記錄技術(shù)以及阻挫現(xiàn)象的研究體系，退磁處理都是人工自旋冰研究過(guò)程中首要的、關(guān)鍵的一步。我們將在隨后的工作中采用電子束曝光等手段[15]制備本文提出的人工自旋冰結(jié)構(gòu)，對(duì)其進(jìn)行磁性表征，研究其豐富的物理性質(zhì)并與本文的理論結(jié)果進(jìn)行比較。

1 Harris M J, Bramwell S T, Mcmorrow D F,. Geometrical frustration in the ferromagnetic pyrochlore Ho2Ti2O7[J]. Physical Review Letters, 1997, 79(13): 2554?2557. DOI: 10.1103/PhysRevLett.79.2554

2 Wills A S, Ballou R, Lacroix C. Model of localized highly frustrated ferromagnetism: the kagomspin ice[J].2002, 66(14): 144407. DOI: 10.1103/ PhysRevB.66.144407

3 Pauling L. The structure and entropy of ice and of other crystals with some randomness of atomic arrangement[J]. Journal of the American Chemical Society, 1935, 57(12): 2680?2684. DOI: 10.1021/ja01315a102

4 Soulas N. Spin-glass models as error-correcting codes[J]. Nature, 1989, 333(6227): 693?695. DOI: 10.1038/ 339693a0

5 Bryngelson J D, Wolynes P G. Spin glasses and the statistical mechanics of protein folding[J]. Proceedings of the National Academy of Sciences of the United States of America, 1987, 84(21): 7524?7528. DOI: 10.1073/pnas. 84.21.7524

6 Toulouse G, Dehaene S, Changeux J P. Spin-glass model of learning by selection[J]. Proceedings of the National Academy of Sciences of the United States of America, 1986, 83(6): 1695?1698. DOI: 10.1073/pnas.83.6.1695

7 范誠(chéng). 磁阻挫材料RMnO3和Dy2Ti2O7的單晶生長(zhǎng)與低溫物性[D]. 合肥: 中國(guó)科學(xué)技術(shù)大學(xué), 2013
FAN Cheng. Crystal growth and low-temperature physical properties of magnetically frustrated RMnO3and Dy2Ti2O7[D]. Hefei: University of Science and Technology of China, 2013

8 Wang R F, Nisoli C, Freitas R S,. Artificial ‘spin ice’ in a geometrically frustrated lattice of nanoscale ferromagnetic islands[J]. Nature, 2006, 439(7074): 303?306. DOI: 10. 1038/nature04447

9 Lammert P E, Ke X L, Li J,. Direct entropy determination and application to artificial spin ice[J]. Nature Physics, 2010, 6(10): 786?789. DOI: 10.1038/ NPHYS1728

10 Qi Y, Brintlinger T, Cumings J. Direct observation of the ice rule in an artificial Kagome spin ice[J]. Physical Review B, 2008, 77(9): 094418. DOI: 10.1103/PhysRevB. 77.094418

11 Mrogan J P, Stein A, Langrige S,. Thermal ground-state ordering and elementary excitations in artificial magnetic square ice[J]. Nature Physics, 2010, 7: 75?79. DOI: 10.1038/NPHYS1853

12 Wang R F, Li J, McConville W,. Demagnetization protocols for frustrated interacting nanomagnet arrays[J]. Journal of Applied Physics, 2007, 101(9): 1?15. DOI: 10.1063/1.2712528

13 Mengotti E, Heyderman L J, Rodríguez A F,. Building blocks of an artificial kagome spin ice: photoemission electron microscopy of arrays of ferromagnetic islands[J]. Physical Review B, 2008, 78: 144402. DOI: 10.1103/PhysRevB.78.144402

14 Moller G, Moessner R. Artificial square ice and related dipolar nanoarrays[J]. Physical Review Letters, 2006, 96(23): 237202. DOI: 10.1103/PhysRevLett.96.237202

15 陳媚, 于懷娜, 趙俊, 等. 光子晶體制備方法對(duì)同步輻射閃爍體探測(cè)器成像分辨率的影響[J]. 核技術(shù), 2015, 38(7): 070101. DOI: 10.11889/j.0253-3219.2015.hjs.38. 070101
CHEN Mei, YU Huaina, ZHAO Jun,. Effect of photonic crystal prepared by different methods on the spatial resolution of scintillator-based synchrotron radiation imaging detector[J]. Nuclear Techniques, 2015, 38(7): 070101. DOI: 10.11889/j.0253-3219.2015.hjs.38. 070101

國(guó)家自然科學(xué)基金(No.11475251、No.11275255、No.11225527、No.5137111)、上海市科委基礎(chǔ)研究重點(diǎn)項(xiàng)目(No.13JC402400)資助

Supported by National Natural Science Foundation of China (No.11475251, No.11275255, No.11225527, No.5137111), and the Key Project of Shanghai Municipal Science and Technology Commission (No.13JC402400)

Monte Carlo simulation on demagnetization of artificial spin ice

YU Liju1MENG Xiangyu2LI Junqin2CAO Jiefeng2WANG Yong2JING Chao1WU Yanqing2TAI Renzhong2

1(Department of Physics, Shanghai University, Shanghai 200444, China)2(Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Zhangjiang Campus, Shanghai 201204, China)

Background: Recently, the artificial spin ice (ASI) draws a wide attention for its exotic phenomenon and potential applications. Previous works mainly studied the Square and Kagome lattice systems, and extensive researches concentrated upon how to achieve the ground state. It turned out to be experimentally difficult to get the absolute ground state of Kagome ASI because of the high frustration. Purpose: Here we propose a new structure based on Kagome lattice, which is just locally frustrated. In order to get good approximation of its ground state, we employed a Monte Carlo simulation on demagnetization process with three methods. Methods: The Monte Carlo simulation was carried on a system including 21×21 cells, and an extra field had been employed to drive the spins overturning. Results & Conclusion: It turned out that the rotating magnetic field was more effective to get to the ground state, and induced the emergence of various magnetic structures.

ASI, Monte Carlo simulation, Demagnetization

YU Liju, female, born in 1989, graduated from Guangxi University in 2013, master student, focusing on magnetic nanostructure

WANG Yong, E-mail: wangyong@sinap.ac.cn; WU Yanqing, E-mail: wuyanqing@sinap.ac.cn; TAI Renzhong, E-mail: tairenzhong@sinap.ac.cn

TL65，O469

10.11889/j.0253-3219.2016.hjs.39.060502

虞麗菊，女，1989年出生，2013年畢業(yè)于廣西大學(xué)，現(xiàn)為碩士研究生，研究領(lǐng)域?yàn)榇判约{米結(jié)構(gòu)

王勇，E-mail: wangyong@sinap.ac.cn；吳衍青，E-mail: wuyanqing@sinap.ac.cn；邰仁忠，E-mail: tairenzhong@sinap.ac.cn

2016-01-25，

2016-04-12