散射理论方程_非弹性中子磁散射方法简介之自旋波激发
1. 前言
凝聚態(tài)物質(zhì)是由大量粒子組成,并且粒子間存在較強(qiáng)相互作用的系統(tǒng)。物質(zhì)的靜態(tài)結(jié)構(gòu)由系統(tǒng)的基態(tài)決定,而物質(zhì)的動力學(xué)性質(zhì)與其激發(fā)態(tài)相關(guān)聯(lián)。前文中我們討論了彈性中子散射方法研究凝聚態(tài)物質(zhì)的晶體結(jié)構(gòu)和磁結(jié)構(gòu),在此我們將討論非彈性中子散射方法研究凝聚態(tài)物質(zhì)的動力學(xué)特性。凝聚態(tài)物質(zhì)的動力學(xué)特性通常表現(xiàn)為集體激發(fā)現(xiàn)象,并且與其熱學(xué)性能和磁性能等物理性質(zhì)直接相關(guān),是深入理解結(jié)構(gòu)相變和磁相變以及溫度依賴現(xiàn)象的重要基礎(chǔ)。物質(zhì)的動力學(xué)特性可以通過對系統(tǒng)元激發(fā)的實(shí)驗(yàn)研究直接得到。非彈性中子散射方法是研究凝聚態(tài)物質(zhì)中元激發(fā),包括聲子激發(fā)、自旋波激發(fā)和晶場激發(fā)等的重要實(shí)驗(yàn)方法。
? 在自旋晶格系統(tǒng)中,磁性離子自旋的有序排列狀態(tài)就是系統(tǒng)的基態(tài),比如常見的鐵磁有序、反鐵磁有序和亞鐵磁有序基態(tài)。如果系統(tǒng)的基態(tài)受到微擾,例如某一個(gè)初始格點(diǎn)上的自旋偏離量子化軸的方向,那么系統(tǒng)將處于一個(gè)低能激發(fā)態(tài)。由于相鄰自旋之間存在相互作用就會使臨近初始格點(diǎn)的自旋取向也發(fā)生改變,并且繼續(xù)影響次近鄰格點(diǎn)自旋,最后以集體偏離軸向的運(yùn)行形式擴(kuò)散和傳播。同時(shí)相鄰自旋也會通過自旋間的相互作用試圖使偏離量子化軸方向的初始格點(diǎn)上的自旋方向恢復(fù)軸向狀態(tài)。最終在物質(zhì)的晶格中會形成自旋的震蕩,表現(xiàn)為自旋在格點(diǎn)上類似波動行為的進(jìn)動。這種自旋的波動就是自旋波,自旋波的量子就是磁振子。磁振子和晶格振動對應(yīng)的聲子一樣,都是凝聚態(tài)物質(zhì)中的集體激發(fā),并服從波色統(tǒng)計(jì)分布規(guī)律,因此也是一種玻色子。通過研究自旋波的色散關(guān)系,基于哈密頓量模型獲得自旋相互作用系統(tǒng)的本征態(tài),我們就能得到自旋之間的相互作用特性。
圖1凝聚態(tài)物質(zhì)中的重要元激發(fā):聲子-格波的量子和磁振子-自旋波的量子。
研究自旋波激發(fā)的關(guān)鍵在于得到準(zhǔn)確的色散關(guān)系,即,也就是哈密頓量的本征值和激發(fā)的動量之間的關(guān)系。因此我們研究自旋波色散關(guān)系就必須應(yīng)用能夠與研究體系發(fā)生能量和動量交換的探針。雖然不同的磁性體系具有不同的磁轉(zhuǎn)變溫度,對應(yīng)著不同的磁相互作用強(qiáng)度和磁激發(fā)能量,但是絕大多數(shù)磁性體系的磁激發(fā)能量一般在幾個(gè)毫電子伏特(meV)到幾百個(gè)毫電子伏特之間。而典型的冷中子和熱中子的能量恰好能夠覆蓋幾個(gè)meV到幾百個(gè)meV的能量范圍。 此外,凝聚態(tài)體系中動量與其原子間距相關(guān),一般都落在幾個(gè)?-1的范圍內(nèi),這也正是中子的可探測范圍。加之中子磁偶極矩的存在,使中子與體系中原子磁矩之間產(chǎn)生相互作用,因此中子是研究磁激發(fā)的最有效的微探針。當(dāng)中子射向凝聚態(tài)體系,并被體系散射之后后,通過建立動量轉(zhuǎn)移和能量轉(zhuǎn)移之間的關(guān)系,我們便可以從實(shí)驗(yàn)上得到體系的自旋波激發(fā)色散關(guān)系,并通過建立理論模型對其色散關(guān)系進(jìn)行分析研究。
2. 非彈性中子散射方法
? 前文中我們已經(jīng)討論過,當(dāng)中子在被體系散射的過程中不發(fā)生能量變化的過程為彈性散射,而非彈性散射過程必然伴隨發(fā)生能量的轉(zhuǎn)移。彈性中子散射過程中的入射矢量ki和出射矢量kf方向不同大小相等(圖2)。而在非彈性中子散射過程中,能量在中子和體系之間發(fā)生轉(zhuǎn)移,散射矢量關(guān)系也相應(yīng)地發(fā)生變化(圖3)。入射矢量和出射矢量與散射矢量以及能量之間的關(guān)系為:
圖2對應(yīng)于倒空間中鍺(400)衍射峰的彈性散射過程的矢量關(guān)系(ki = kf)。
圖3對應(yīng)于倒空間中鍺Q= (4.6,0.6,0)位置,能量為10meV的非彈性散射過程的矢量關(guān)系(ki≠ kf)。
通過非彈性中子散射方法可測得的動量轉(zhuǎn)移和能量轉(zhuǎn)移范圍與入射中子的能量密切相關(guān),因此在儀器硬件不受限的條件下首先要保證入射矢量、出射矢量和散射矢量之間能夠形成三角閉環(huán),其次要保證在特定的布里淵區(qū)能夠?qū)崿F(xiàn)預(yù)期的能量轉(zhuǎn)移。非彈性中子散射方法的應(yīng)用,特別是在自旋波激發(fā)方面的應(yīng)用主要可以通過三軸中子譜儀和飛行時(shí)間中子譜儀這兩類譜儀裝置得以實(shí)現(xiàn)。
2.1三軸中子譜儀
? 三軸中子譜儀是最常用的研究自旋波激發(fā)的實(shí)驗(yàn)裝置。三軸中子譜儀的發(fā)明人Bertram N. Brockhouse 因其在中子譜學(xué)發(fā)展中的重要貢獻(xiàn),和發(fā)展了中子衍射技術(shù)的CliffordG. Shull共同獲得了1994年的諾貝爾物理學(xué)獎(jiǎng)。顧名思義,中子三軸譜儀具有三個(gè)旋轉(zhuǎn)軸,在三個(gè)旋轉(zhuǎn)軸上分別放置三個(gè)重要部件,依次是單色器、樣品和分析晶體。旋轉(zhuǎn)單色器可以得到不同波長/能量的入射中子,旋轉(zhuǎn)樣品可以得到不同的散射角和入射中子波矢,旋轉(zhuǎn)分析晶體可以得到不同波長/能量的散射中子。通過三個(gè)軸的聯(lián)動旋轉(zhuǎn)我們就能夠得到在倒空間平面中各個(gè)位置的能量變化。在三軸譜儀的實(shí)際應(yīng)用中,我們一般通過分析晶體固定出射中子的能量,通過改變動量轉(zhuǎn)移或者能量轉(zhuǎn)移來分別獲得固定能量掃描和固定動量掃描實(shí)驗(yàn)結(jié)果。一般來說三軸中子譜儀都是依托能夠提供穩(wěn)定中子束流的反應(yīng)堆中子源搭建,通過中子光學(xué)器件不但可以對中子的能量進(jìn)行選擇并可以有效地聚焦中子增加束流強(qiáng)度。由于三軸譜儀的散射幾何清晰明確,因此一般相比飛行時(shí)間中子譜儀具有更好的空間和能量分辨率。
圖4位于法國格勒諾布爾的勞厄-朗之萬研究所(ILL)的冷中子三軸譜儀IN12的示意圖和實(shí)物圖。
采用三軸中子譜儀在倒空間中的掃描一般都是線掃描(固定Q掃描或固定E掃描),如果想得到整個(gè)倒空間的激發(fā)信號就必須進(jìn)行逐個(gè)線掃描然后整合在一起,這個(gè)測試過程非常耗時(shí),因此人們也發(fā)展了多路器分析技術(shù),通過增加分析晶體建立更多的分析通道,覆蓋更大的測試區(qū)間,提高測試效率。在利用非彈性中子散射技術(shù)研究體系元激發(fā)的過程中,我們需要首先建立實(shí)空間和倒空間的轉(zhuǎn)換關(guān)系并通過矩陣參數(shù)切換,通常非彈實(shí)驗(yàn)的測試范圍還受到三軸譜儀的硬件和空間限制,通過vTAS模擬軟件可以直觀地展示三軸譜儀布局和相應(yīng)測試區(qū)間的矢量關(guān)系,對于理解非彈性中子散射實(shí)驗(yàn)和進(jìn)行實(shí)驗(yàn)設(shè)計(jì)非常有幫助。
根據(jù)入射中子的能量以及相應(yīng)的能量探測范圍,三軸中子譜儀一般可以分為冷三軸和熱三軸。目前我國也已經(jīng)依托研究反應(yīng)堆建設(shè)了數(shù)臺三軸中子譜儀,包括四川綿陽的中國綿陽研究堆(CMRR)的冷中子三軸譜儀“鯤鵬”;中國先進(jìn)研究堆(CARR)的中德熱中子三軸譜儀SV30(固定ki模式),中科院物理所依托中國先進(jìn)研究堆(CARR)建設(shè)的熱中子三軸譜儀“翠竹”,中國人民大學(xué)依托中國先進(jìn)研究堆(CARR)建設(shè)的冷中子三軸極化譜儀“行知”以及中國人民大學(xué)依托中國先進(jìn)研究堆(CARR)建設(shè)的冷中子廣譜譜儀“博雅”。這些譜儀在凝聚態(tài)物質(zhì)的動力學(xué)研究方面已經(jīng)取得了許多優(yōu)秀的科研成果。
三軸中子譜儀目前在技術(shù)方面已經(jīng)發(fā)展地非常成熟,Gen Shirane等在Neutron Scattering witha Triple-Axis Spectrometer Basic Techniques這本書里對三軸譜儀的原理和應(yīng)用做了非常詳盡的描述,特別是對分辨率方程(resolution function)和偽峰(spurious peaks)方面的講解是同時(shí)基于散射原理和實(shí)踐經(jīng)驗(yàn)的全面總結(jié),因此是依托三軸譜儀完成非彈性中子散射實(shí)驗(yàn)的重要參考。具體內(nèi)容可參考該書,本文中就不展開討論了。
2.2飛行時(shí)間中子譜儀
? 由于中子的能量和它的速度平方成正比,因此不同飛行速度的中子具有不同的能量。中子在被材料體系散射后,我們通過測量被樣品散射的中子在樣品和探測器之間已知的固定距離的飛行時(shí)間,就能夠得到被散射中子的能量,從而判斷中子和樣品之間的能量交換,實(shí)現(xiàn)對非彈性中子散射的能量分析。然而采用飛行時(shí)間方法用到的中子脈沖都是微秒級別,在第一束中子到達(dá)探測器之后下一束中子脈沖才能開始飛行,因此這一方法更適用于脈沖式散裂中子源。雖然在反應(yīng)堆式的穩(wěn)態(tài)中子源也可以通過加裝費(fèi)米斬波器應(yīng)用飛行時(shí)間中子譜儀,但是中子的利用效率較低,一般只有2%左右。飛行時(shí)間中子譜儀可以有兩種設(shè)置,一種是直接幾何,另一種是非直接幾何,或者稱為逆幾何(圖5)。在直接幾何非彈譜儀中,入射中子首先通過斬波器被單色化,在和樣品發(fā)生散射后又被不同散射角度的探測器探知。探測器記錄散射中子的角度位置以及到達(dá)探測器的時(shí)間,并據(jù)此得到散射中子的波矢和能量。在非直接幾何譜儀中,寬波段的中子射向樣品,隨后通過一系列分析晶體對散射中子的能量進(jìn)行分析。雖然我國還沒有可用的非彈性中子散射飛行時(shí)間譜儀,但是目前中山大學(xué)和散裂中子源科學(xué)中心正在依托中國散裂中子源聯(lián)合建設(shè)高能直接幾何非彈性中子散射飛行時(shí)間譜儀,該譜儀預(yù)計(jì)2022年底完成譜儀建設(shè)并面向用戶開放。此外,依托中國散裂中子源的非直接幾何譜儀非彈性中子譜儀也在規(guī)劃建設(shè)中。
圖5直接幾何和非直接幾何非彈性中子散射飛行時(shí)間譜儀的中子飛行路徑圖。
3. 自旋波激發(fā)理論模型
? 具有自發(fā)磁矩的鐵磁體中所有的自旋都平行排列,因此其磁晶胞和化學(xué)晶胞大小一致,并且交換常數(shù)大于零。通過鐵磁自旋波理論可以對低能海森堡哈密頓量近似對角化得到自旋波解,其自旋波具有二次方色散特征(圖6)。具有自旋反平行排列特征的反鐵磁體的磁晶胞一般大于其化學(xué)晶胞,并且交換常數(shù)小于零。通過反鐵磁自旋波理論結(jié)合雙子晶格模型可以對低能海森堡哈密頓量近似對角化得到自旋波解,其自旋波具有線性色散特征(圖6)。自旋波激發(fā)理論的詳細(xì)推導(dǎo)可以很容易在固體物理或磁性物理教科書中找到,這里不再詳述。
圖6鐵磁體和反鐵磁體的低溫自旋波色散關(guān)系,D0為各向異性常數(shù)。
4. 自旋波激發(fā)研究實(shí)例
?? 同時(shí)具有鐵磁性和鐵電性的多鐵材料一直以來都是材料科學(xué)和凝聚態(tài)物理研究領(lǐng)域關(guān)注的熱點(diǎn),特別是具有較強(qiáng)磁電耦合的多鐵材料更是蘊(yùn)含了豐富的物理現(xiàn)象和物理內(nèi)涵。MnWO4化合物是一種典型的自旋驅(qū)動鐵電材料,因此深入研究該多鐵材料的自旋動力學(xué)不但有助于理解其多鐵特性并且有助于探索本征的磁電耦合機(jī)理。MnWO4為單斜P2/c結(jié)構(gòu),高自旋的Mn2+離子是晶胞中唯一的磁性離子,具有反鐵磁的磁基態(tài),傳播矢量為k= (±1/4, 1/2, 1/2),磁結(jié)構(gòu)如下圖所示。
圖7MnWO4多鐵化合物的低溫晶體結(jié)構(gòu)和反鐵磁磁結(jié)構(gòu)示意圖。
? 通過MnWO4化合物13.5K的磁轉(zhuǎn)變溫度可以判斷其自旋波激發(fā)主要集中在低能量范圍,因此可以選取冷中子非彈性中子散射譜儀對其自旋波激發(fā)進(jìn)行研究。位于德國慕尼黑海因茨·邁爾-萊布尼茨中心(MLZ)反應(yīng)堆,由德國于利希中子科學(xué)中心(JCNS)運(yùn)行的冷中子三軸譜儀PANDA和位于美國橡樹嶺國家實(shí)驗(yàn)室(ORNL)散裂中子源(SNS)的冷中子飛行時(shí)間譜儀CNCS都是理想的選擇。首先采用浮區(qū)單晶生長方法獲得具有5克左右質(zhì)量的大塊單晶樣品,并用勞厄衍射儀對單晶樣品進(jìn)行定向,使晶體的[102]和[010]方向保持在散射平面內(nèi),這樣就能夠通過非彈性中子實(shí)驗(yàn)獲得磁布里淵區(qū)內(nèi)沿著高對稱方向的各點(diǎn),如圖8所示。在采用冷中子三軸譜儀PANDA展開實(shí)驗(yàn)的過程中我們采用PG(002)單色器,同時(shí)在分析晶體前放置Be過濾器用于過濾高階諧波中子。實(shí)驗(yàn)采用固定kf = 1.2 ?-1的模式,并有散射矢量 Q =(qx,qy,qz)?-1, (HKL) = (qxa/2π,qyb/2π,qzc/2π)。圖9是分別沿[102]和[010]方向測得的磁激發(fā)信號。
圖8MnWO4單晶及其定向示意圖。
圖9沿[102]和[010]方向的磁激發(fā)信號, (a)(c)實(shí)驗(yàn)數(shù)據(jù)(b)(d)理論計(jì)算。
? 通過線性自旋波近似可以對非彈性中子散射數(shù)據(jù)進(jìn)行分析。首先,建立包含自旋交換部分和各向異性部分的海森堡哈密頓量:
其中Jij代表圖7中所示從J1到J12的自旋交換常數(shù)。Ds代表單軸單離子各向異性常數(shù)。MnWO4化合物共線反鐵磁基態(tài)的磁單胞中總共有8個(gè)自旋,其中4個(gè)向上,4個(gè)向下。引入升降算符后可以對量子自旋進(jìn)行線性的霍斯坦因-普里馬科夫變換,有:
式中和以及和都滿足對易關(guān)系。對玻色算符進(jìn)行傅里葉變換有:
獲得動量空間中的玻色哈密頓量為:
其中由4個(gè)自共軛的厄米矩陣組成,有:
為了獲得矩陣元素值需要對包含每個(gè)交換常數(shù)的項(xiàng)進(jìn)行歸納整理,下面是對最近鄰項(xiàng)J1以及其它幾個(gè)交換常數(shù)的推導(dǎo)示例。
最后得到的矩陣元素為:
其中部分矩陣元素存在共軛關(guān)系,歸納整理后得到:
矩陣中包含各交換常數(shù)的矩陣元素為:
在取得玻色哈密頓量矩陣之后,不能直接用對角化厄米的方法對角化哈密頓量矩陣。需要借助不改變玻色對易關(guān)系的轉(zhuǎn)換矩陣實(shí)現(xiàn)對角化,并得到實(shí)數(shù)本征值。這樣就能夠取得在倒空間中分別與自旋波量子數(shù)和自旋波態(tài)向量對應(yīng)的本征值和本征矢量。這里我們用的是線性自旋波近似方法,而沒有考慮自旋波之間的相互作用,如果把相互作用部分也加以考慮的話就需要解算符的非線性方程,在一維情況下該非線性方程有孤子解,代表非線性元激發(fā)。
? ?結(jié)合非彈性中子散射截面就可以在線性自旋波理論模型的基礎(chǔ)上對實(shí)驗(yàn)中子譜進(jìn)行擬合。非彈性中子散射截面為:
其中kf和ki分別為入射矢量和出射矢量,= 0.539×10-12 cm=5.39 fm為單電子磁散射強(qiáng)度,g為朗德因子,為Mn2+的磁形狀因子。為德拜-沃勒因子,是自旋關(guān)聯(lián)方程。由于MnWO4的磁單胞中有8個(gè)Mn自旋,因此在零場下自旋波應(yīng)該具有4支,每一支都是二重簡并的。通過擬合交換常數(shù)我們可以得到和實(shí)驗(yàn)數(shù)據(jù)非常吻合的計(jì)算自旋波激發(fā)譜,如圖9所示,相應(yīng)地我們可以得到一系列代表自旋相互作用大小和性質(zhì)的自旋交換常數(shù)。
圖10在飛行時(shí)間中子譜儀上測得的沿[102]方向通過(1/4,1/2, 1/2)磁峰(a)和(1/4, 1, 1/2)位置(b)的自旋波譜。沿[10-2]方向通過(1/4,1/2, -1/2)磁峰(c)和(1/4, 1, -1/2)位置(d)的自旋波譜。
圖11與圖10中實(shí)驗(yàn)自旋波譜對應(yīng)的計(jì)算自旋波譜。
? 這些自旋波譜也可以通過配備費(fèi)米斬波器的冷中子飛行時(shí)間譜儀獲得(圖10),入射中子的波長為λ = 4.4 ? 。飛行時(shí)間中子譜儀的優(yōu)勢在于能夠快速完成數(shù)據(jù)采集工作,通過單次測量就可以覆蓋大面積Q-E空間,而三軸中子譜儀一次只能測得Q-E空間中的一點(diǎn)。由于具有出色的矩陣運(yùn)算和可視化作圖能力,MATLAB是分析自旋波激發(fā)中子數(shù)據(jù)的優(yōu)秀軟件之一。附件附上了MATLAB自旋波計(jì)算代碼供參考。采用擬合得到的交換常數(shù)我們可以計(jì)算得到布里淵區(qū)中各個(gè)方向的自旋波色散關(guān)系或者各個(gè)能量的自旋波分布情況,例如與實(shí)驗(yàn)自旋波譜對應(yīng)的沿著 [1 0 2]和[10 ?2]方向的自旋波譜(圖11),或者任意幾個(gè)切面的色散關(guān)系(圖12)。
圖12MnWO4在布里淵區(qū)任意面上的自旋波激發(fā)譜。
? ?自旋波是磁有序系統(tǒng)中常見的激發(fā)模式,可以用經(jīng)典的海森堡模型進(jìn)行描述。然而一些復(fù)雜磁系統(tǒng)中的相互作用或耦合會導(dǎo)致不能被海森堡模型所描述的復(fù)雜的磁激發(fā)行為出現(xiàn),比如動態(tài)磁電耦合會導(dǎo)致新的集體激發(fā)現(xiàn)象出現(xiàn),磁性三角晶格中幾何阻搓引起的多重簡并態(tài)磁激發(fā)等,這些都是引起廣泛關(guān)注并值得深入研究的課題。
參考文獻(xiàn)
[1]T. Chatterji (Editor), Neutron Scattering from MagneticMaterials(2006) Elsevier.
[2] Th. Brückel et al., Neutron Scattering (2008), Forschungszentrum Jülich, Jülich.
[3] G. Shirane, S. M.Shapiro, and J. M. Tranquada, Neutron Scattering with a Triple-AxisSpectrometer Basic Techniques (2004) Cambridge University Press.
[4]M. Boehm, et al.,Nuclear Inst. and Methods in Physics Research, A(2013)?697, 40-44.
[5]F. Ye,et al., Phys. Rev. B 83,140401(R) (2011).
[6]Y.Xiao, et al., Phys. Rev. B 93,214428 (2016).
[7]B.-Q.Liu, et al., J. Phys.: Condens. Matter 30, 295401 (2018).
[8] J. Wang, et al., Phys. Rev. B 98, 214425 (2018).
總結(jié)
以上是生活随笔為你收集整理的散射理论方程_非弹性中子磁散射方法简介之自旋波激发的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: boot spring 怎么执行hql_
- 下一篇: dubbo服务接口如何mock_2019