李曉杰,劉延春,李新軍

(齊魯理工學(xué)院，濟南 250200)

1 引 言

自1966年，Blume-Capel(BC)模型被創(chuàng)建之后[1,2]，人們已經(jīng)利用不同方法對多種晶格上的BC模型的磁化性質(zhì)、熱力學(xué)性質(zhì)、磁滯回線和相圖進(jìn)行了研究.Y.F.Zhang小組研究了外磁場和交換相互作用都遵循雙模隨機分布時，簡立方晶格中混合自旋系統(tǒng)的臨界行為[3].同年里，該小組又研究了外磁場遵循三模隨機分布時，簡立方晶格中混合自旋系統(tǒng)的相變行為[4].文獻(xiàn)[5]利用有效場理論研究了簡立方晶格中BC模型的補償行為和磁化過程.文獻(xiàn)[6]中研究了外磁場服從雙峰離散分布時蜂巢晶格的相變性質(zhì)，發(fā)現(xiàn)外磁場、晶場和自旋間交換相互作用影響系統(tǒng)的相變并且系統(tǒng)出現(xiàn)重入現(xiàn)象；文獻(xiàn)[7]的研究表明稀釋晶場對蜂巢晶格系統(tǒng)磁學(xué)性質(zhì)和相變的影響，結(jié)果顯示當(dāng)晶場滿足稀釋分布時對系統(tǒng)的相變沒有影響并且系統(tǒng)不會出現(xiàn)三臨界現(xiàn)象.近幾年來，納米管漸漸成為磁熱性質(zhì)研究領(lǐng)域的一個熱點，實驗與理論方面都已經(jīng)取得一定的研究結(jié)果.實驗上，毛瑞等人以植物纖維素為模板，制備出了納米管狀SnO2材料，測試結(jié)果顯示，此SnO2納米管狀材料能夠提高鋰離子的擴散速率，有效解決解充電放電過程中電極材料體積膨脹問題[8].文獻(xiàn)[9]中發(fā)現(xiàn)Fe-Ni磁性納米管具有明顯的各項異性；理論上，Zaim小組研究了外磁場滿足三模分布時納米管上自旋為1Ising模型的磁性質(zhì)和相圖[10]，結(jié)果表明該系統(tǒng)具有一階相變、三相臨界點和二階相變并出現(xiàn)重入現(xiàn)象.Osman等人分別討論了納米管中純自旋系統(tǒng)的磁化強度、磁化率、熱力學(xué)性質(zhì)和相圖[11]，討論了晶場對系統(tǒng)磁熱性質(zhì)的影響，結(jié)果顯示系統(tǒng)存在一階相變和二階相變.文獻(xiàn)[12]利用有效場理論研究混合自旋納米管系統(tǒng)的磁熱特性，結(jié)果顯示系統(tǒng)存在三臨界點.文獻(xiàn)[13]研究了Ising納米管系統(tǒng)中BC模型的磁滯行為.T.Kaneyoshi 討論了納米管系統(tǒng)中外殼層和內(nèi)殼層最近鄰自旋間交換相互作用對磁化率的影響[14].文獻(xiàn)[15]顯示雙模隨機晶場中BC模型的磁化強度和相變性質(zhì)，得到了系統(tǒng)的磁化強度與溫度和雙模隨機晶場的關(guān)系以及相圖，結(jié)果表明系統(tǒng)在雙模隨機晶場中會表現(xiàn)出不同的磁學(xué)性質(zhì)和相變行為.文獻(xiàn)[16]討論了納米管上BEG模型的熱力學(xué)和相變性質(zhì)，研究發(fā)現(xiàn)納米管系統(tǒng)存在三臨界點.文獻(xiàn)[17]研究了稀釋晶場作用下納米管中BC模型的磁化性質(zhì)，結(jié)果表明，稀釋晶場作用下系統(tǒng)的內(nèi)能、比熱和自由能呈現(xiàn)出不同的磁學(xué)性質(zhì).文獻(xiàn)[18]利用基于密度泛函理論的第一性原理方法研究了稀土金屬La吸附摻雜BN納米管的儲氫性能.據(jù)我們所知，迄今人們還沒有研究外磁場、晶場參數(shù)比值和雙模隨機交錯晶場對納米管系統(tǒng)平均磁化強度的影響.為了弄清楚外磁場、晶場參數(shù)比值和雙模隨機交錯晶場對納米管系統(tǒng)磁化性質(zhì)的影響，本文利用有效場理論對納米管上混合自旋BC模型格點的平均磁化強度與雙模隨機交錯晶場取值概率、外磁場、晶場參數(shù)比值和晶場的關(guān)系進(jìn)行了研究，給出了平均磁化強度隨溫度的變化曲線.

2 模型和方法

無限長納米管由內(nèi)殼層和外殼層兩部分組成，見圖1.圖1(a)為納米管的立體示意圖，圖1(b)給出其橫截面示意圖.為了更清晰地顯示不同格點上具有相同配位數(shù)的磁性原子，用圓圈、方塊和三角形分別代表配位數(shù)為5、6和7的磁性原子.外殼層每個磁性原子的自旋為1，內(nèi)殼層每個磁性原子的自旋為1/2.圖中原子間的連線代表最近鄰磁性原子之間存在交換相互作用，它們的大小分別為J1、J2和J.

圖1 納米管示意圖.(a) 立體圖，(b) 截面圖，圓圈和方塊代表外殼層磁性原子，三角代表內(nèi)殼層磁性原子，實線表示最近鄰自旋間的交換相互作用.Fig.1 Theschematic pictures of nanotube:(a) perspective view of the cylindrical nanotube,(b) its cross section.The circles and squares respectively represent magnetic atoms at the surface shell.The triangles are magnetic atoms constituting the core shell.The bonds connecting the magnetic atoms represent the nearest-neighbor exchange interactions.

納米管上BC模型的哈密頓量可表述為

(1)

其中S=-1；σ=±1/2.J1代表外殼層最近鄰自旋間的交換相互作用，J2代表內(nèi)殼層最近鄰自旋間的交換相互作用，J2代表外殼層和最近鄰的內(nèi)殼層自旋間的交換相互作用，h代表外磁場，Di代表作用在格點i上的隨機交錯晶場，它們都滿足雙模晶場分布

P(Di)=pδ(Di-D)+(1-p)δ(Di-αD),

(2)

其中p(0≤p≤1)表示隨機晶場取值為D的概率，1-p代表示隨機晶場取值為αD的概率，α(-1≤α<0)為一無量綱參數(shù)，代表晶場強度的比值.當(dāng)p=1或0時，含雙模交錯晶場的BC模型退化為含恒定晶場的BC模型.

利用有效場理論[19-21]可以得到外殼層格點自旋磁化強度m1和m2，內(nèi)殼層格點自旋磁化強度mc的自洽方程：

(3a)

(3b)

(3c)

其中函數(shù)F(x)定義為

=pf(x,D)+(1-p)f(x,αD),

(4)

(5)

函數(shù)f(x,Di)的表達(dá)式為

(6)

其中β=1/kBT，T是絕對溫度，kB是玻爾茲曼常數(shù).另外，為了整體上描述系統(tǒng)的磁學(xué)性質(zhì)可引入平均磁化強度[22]，它的定義為

(7)

3 結(jié)果與討論

為了便于研究，不特別說明，本文余下內(nèi)容均取J1/J=J2/J=1.為了不失一般性，令外磁場h、晶場強度D和等效溫度kBT以J為單位，通過求解方程(3)，給出了雙模隨機交錯晶場作用下系統(tǒng)平均磁化強度隨溫度的變化曲線.

3.1 類對稱交錯晶場作用下，外磁場為零時，概率對系統(tǒng)平均磁化強度的影響

圖2 外磁場為0，晶場強度比值為-1，晶場參數(shù)D/J分別為(a) 10.0，(b) 2.7，(c) -2.7，(d) -10.0時，系統(tǒng)平均磁化強度隨溫度的變化曲線，曲線上的值為隨機晶場概率p.Fig.2 The temperature dependence of the magnetization is presented with some selected values of crystal field D/J (a) 10.0，(b) 2.7，(c) -2.7，(d) -10.0,when the external magnetic field is 0 and the value of D/J=-1.The real number on each curve denotes the value of p.

圖3 概率p=0.3和0.9時，系統(tǒng)平均磁化強度隨溫度的變化曲線，曲線上的值為晶場參數(shù)D/JFig.3 The temperature dependence of the magnetization is presented with some selected values of crystal field D/J,when the external magnetic field is not zero.The real number on each curve denotes the value of D/J

3.2 類對稱交錯晶場作用下，晶場參數(shù)和外磁場對系統(tǒng)平均磁化強度的影響

為了研究不同晶場強度下系統(tǒng)平均磁化強度隨溫度的變化情況，圖3(a)和(b)分別給出了p=0.3和0.9時，雙模隨機交錯晶場分布下平均磁化強度隨溫度的變化曲線.結(jié)果表明，p較小時，負(fù)晶場對磁化強度阻礙作用不明顯；晶場參數(shù)取正值時，系統(tǒng)發(fā)生一級相變；系統(tǒng)平均磁化強度出現(xiàn)缺失現(xiàn)象(平均磁化強度為零)；當(dāng)D/J=-3.93且溫度較低時，磁化強度隨溫度的變化非常明顯，通過研究，我們認(rèn)為是因為外殼層格點(藍(lán)色圓點代表的磁性原子)的配位數(shù)較小的原因.隨著p增大，負(fù)晶場對系統(tǒng)磁化強度抑制作用明顯，晶場參數(shù)取負(fù)值時，系統(tǒng)發(fā)生一級相變.雖然不同晶場強度作用下，系統(tǒng)平均磁化強度呈現(xiàn)出不同的磁化現(xiàn)象，但總體變化趨勢相似.為了研究外磁場對平均磁化強度缺失現(xiàn)象的影響，我們給出了圖3(c)和(d).從圖中可以看出，外磁場能夠減弱平均磁化強度缺失現(xiàn)象，但是缺失部分平均磁化強度呈現(xiàn)負(fù)值.

3.3 類非對稱交錯晶場作用下，外磁場為零時，概率對系統(tǒng)平均磁化強度的影響

圖4(a)—(f) 給出了晶場滿足類非對稱交錯晶場(例如α=-0.5)分布時，系統(tǒng)平均磁化強度隨溫度的變化曲線.圖4(a)和(b)及圖4(d)和(f)顯示，系統(tǒng)平均磁化強度隨溫度的變化曲線與圖2相似.但是當(dāng)D/J=-5.8時，系統(tǒng)平均磁化強度隨溫度變化曲線局部出現(xiàn)負(fù)值(例如p=0.7,0.9,1.0)，見圖4(e)，說明原子摻雜程度(晶場強度比值α)不同，會使系統(tǒng)呈現(xiàn)出不同的磁化現(xiàn)象.

圖4 外磁場為0，晶場參數(shù)D/J分別為(a) 14，(b) 5.8，(c) 2.7，(d) -2.7，(e) -5.8，(f) -14時，系統(tǒng)平均磁化強度隨溫度的變化曲線，曲線上的值為隨機晶場概率p.Fig.4 The temperature dependence of the magnetization is presented with some selected values of D/J=14,5.8,2.7,-2.7,-5.8,-14,when the external magnetic field is 0.The real number on each curve denotes the value of p.

3.4 類非對稱交錯晶場作用下，晶場參數(shù)和外磁場對系統(tǒng)平均磁化強度的影響

為了研究不同晶場強度比值對系統(tǒng)平均磁化強度的影響，圖5給出了類非對稱交錯晶場作用下平均磁化強度隨溫度的變化曲線.結(jié)果顯示：概率較小和較大時，系統(tǒng)都存在一級相變和二級相變；外磁場較小時，沒有對系統(tǒng)平均磁化強度產(chǎn)生明顯的影響.但是系統(tǒng)平均磁化強度出現(xiàn)缺失和負(fù)值的現(xiàn)象更加明顯.我們認(rèn)為導(dǎo)致這種現(xiàn)象的物理原理是不同原子摻雜，使得系統(tǒng)在一定條件下平均磁化強度出現(xiàn)負(fù)值.

圖5 概率 p分別為0.3，0.9時，系統(tǒng)磁化強度隨溫度的變化曲線，曲線上的值為晶場參數(shù)D/JFig.5 The temperature dependence of the magnetization is presented with some selected values of the probability of p =0.3,0.9,when the external magnetic field is not zero.The real number on each curve denotes the value of D/J.

4 結(jié) 論

本文利用有效場理論研究了雙模隨機交錯晶場作用下spin-1和spin-1/2混合自旋BC模型系統(tǒng)的平均磁化強度隨溫度的變化情況.結(jié)果表明，系統(tǒng)格點的平均磁化強度與隨機交錯晶場的取值概率p、晶場強度比值α、晶場參數(shù)D/J、外磁場、溫度以及交換相互作用密切相關(guān).概率P、外磁場、晶場強度比值和D等諸多因素相互競爭，使系統(tǒng)表現(xiàn)出比恒定晶場作用的BC模型更為豐富的磁化現(xiàn)象：系統(tǒng)格點的平均磁化強度隨溫度的變化曲線呈現(xiàn)缺失和負(fù)值的現(xiàn)象；P、α和D/J會影響系統(tǒng)的一級相變和二級相變.