賈維燁, 葛崇員, 孫樹生, 成泰民

(沈陽化工大學(xué) 數(shù)理系, 遼寧 沈陽 110142)

低維磁性系統(tǒng)存在許多奇異的特性,因此這一領(lǐng)域的研究引起了科學(xué)家們廣泛的興趣[1-5].1997年Shiomi[1]等人通過實驗合成出自旋為1/2交替的亞鐵磁性復(fù)合物.2005年Kikuchi[2]等人通過實驗發(fā)現(xiàn)Cu3(CO3)2(OH)2材料的磁化強度在低溫下隨外磁場的變化呈現(xiàn)1/3的磁化平臺,在低溫區(qū)磁化率隨溫度變化出現(xiàn)雙峰等特性.2008年Rule[4]等人利用非彈性中子散射實驗研究了藍(lán)銅礦Cu3(CO3)2(OH)2材料的諸交換積分，并通過實驗也發(fā)現(xiàn)了1/3的磁化平臺的存在.2011年Jeschke[5]利用密度泛函理論從理論上研究了自旋為1/2的棱型鏈磁性質(zhì)及諸交換積分，并利用第一性原理解釋了這一阻挫材料的形成機理.

長程相互作用對系統(tǒng)熱力學(xué)性質(zhì)的影響近年來也被人們廣泛關(guān)注[6-7],外場中具有長程相互作用自旋為1/2的亞鐵磁棱型X-Y鏈的熱力學(xué)性質(zhì)也得到了相應(yīng)研究.本文采用Jordan-Wigner(J-W)變換[8]分析方法對外場中Z方向具有均勻長程相互作用自旋為1/2的亞鐵磁棱型X-Y鏈進行研究.利用不變本征算符(IEO)法[9-14]得到體系的元激發(fā)能譜，對體系的Hamiltonian進行對角化,并利用量子統(tǒng)計理論推導(dǎo)出體系配分函數(shù)和內(nèi)能.

1 系統(tǒng)的Hamiltonian對角化及體系的內(nèi)能

二聚體分子和單基體分子交替排列構(gòu)成的自旋為1/2的亞鐵磁棱型X-Y鏈如圖1所示.其中S2,l和S3,l表示二聚體分子的量子自由基,S1,l表示單基體分子的量子自由基.JF=-ηJ(η,J>0)表示二聚體分子內(nèi)的鐵磁相互作用,JAF=μJ(μ,J>0)表示二聚體分子與單基體分子間的反鐵磁相互作用,并且假設(shè)這樣的循環(huán)結(jié)構(gòu)具有N個(l=1,2,…,N).

圖1 自旋為1/2的亞鐵磁棱型X-Y鏈系統(tǒng)

上述相互作用利用Heisenberg相互作用理論模型描述,那么系統(tǒng)的Hamiltonian可用如下形式表示：

(1)

式(1)中：Si,l是自旋1/2的量子算符,S2，l和S3,l由二聚體分子內(nèi)的鐵磁相互作用JF<0耦合連接,其他作用都由二聚體分子與單基體分子間的反鐵磁相互作用JAF>0耦合連接,uB為玻爾磁子,g為朗德因子,H表示外磁場強度.

對式(1)的自旋升降算符進行Jordan-Wigner(J-W)變換[8,11-14]：

(2)

(3)

(4)

其中：

(5)

其中：h=guBH,h.c.表示前項的厄米共軛,并且γk=exp(ik).

令系統(tǒng)哈密頓量的IEO算符為：

τ3,σa3,k, (σ=1,2,3)

(6)

由式(6)可得:

(7)

根據(jù)式(6)及利用IEO法可得自旋為1/2的亞鐵磁棱型X-Y鏈系統(tǒng)的元激發(fā)能量:

(8)

將式(6)和(7)代入式(5)可得：

(9)

根據(jù)式(4)和(9),利用量子統(tǒng)計理論可得系統(tǒng)的配分函數(shù)：

(10)

其中β=1(kBT) .

由式(10)可得系統(tǒng)的內(nèi)能：

(11)

2 討 論

HC1=ηJ/(2guB)

(12)

(13)

Hpeak=ηJ/(4guB)

(14)

當(dāng)Hpeak≤H≤HC2時，Hpeak是外磁場下體系的磁化強度隨溫度的變化曲線出現(xiàn)峰值的臨界磁場強度.

根據(jù)式(8)、(11)通過數(shù)值計算可得圖2～圖5結(jié)果.

由圖2、圖3可知，體系的內(nèi)能隨外加磁場的增加逐漸減少.由圖2可知:當(dāng)體系處于絕對零度(kBT/JAF=0)時,二聚體分子的鐵磁相互作用JF與單基體與二聚體分子之間的反鐵磁相互作用JAF之比(-JF/JAF)越大，在低外磁場區(qū)體系的內(nèi)能越小.但是在高外磁場區(qū)(即guBH/JAF>2.5)該比值(-JF/JAF)對體系的內(nèi)能幾乎無影響，在此表現(xiàn)為對于不同的JF/JAF體系內(nèi)能隨外磁場變化曲線，在高外磁場區(qū)重分為一條曲線.由圖3可知：在有限溫度下(kBT/JAF=0.5)，在低外磁場區(qū)(-JF/JAF)越大，也出現(xiàn)體系內(nèi)能減少的現(xiàn)象.但是與絕對零度不同的是，對于不同的-JF/JAF，體系內(nèi)能在有限溫度下隨外磁場的變化曲線在更高的外磁場下(即guBH/JAF=3.3)出現(xiàn)重合，且在大外磁場下(即guBH/JAF>4.0)時，曲線重合出現(xiàn)分裂的跡象.這說明在此存在如下4種形式的能量之間的競爭：與溫度相關(guān)的熱無序能(kBT/JAF)、與外磁場強度相關(guān)的自旋磁矩勢能(guBH/JAF)；使二聚體分子內(nèi)電子自旋平行排列的鐵磁交換作用能(JF)；使二聚體與單基體分子間電子自旋反平行排列的反鐵磁交換作用能(JAF).

圖2 在絕對零度下(kBT/JAF=0)對于不同 JF/JAF體系內(nèi)能隨外磁場的變化

圖3 在有限溫度下(kBT/JAF=0.5)對于不同 JF/JAF體系內(nèi)能隨外磁場的變化

圖4 在無外磁場下(guBH/JAF=0)對于不同的 JF/JAF體系內(nèi)能隨溫度的變化

圖5 在有限外磁場下(guBH/JAF=0.2)對于不同的 JF/JAF體系內(nèi)能隨溫度的變化

由圖4、圖5可知：體系的內(nèi)能隨溫度的升高逐漸增大，并且在低溫區(qū)出現(xiàn)U～T曲線向上彎曲和在高溫區(qū)出現(xiàn)向下彎曲的現(xiàn)象.這說明在低溫區(qū)，自旋為1/2的亞鐵磁棱型X-Y鏈體系容易出現(xiàn)比熱異常的峰，從而將出現(xiàn)比熱雙峰現(xiàn)象.但是，比較圖4和圖5可知在有限外加磁場下，體系的內(nèi)能比無外加磁場時更小.

3 結(jié) 論

在有限溫度與有限外磁場下,對于自旋為1/2的亞鐵磁棱型X-Y鏈體系內(nèi)能的研究得出如下主要結(jié)果:

(2) 外磁場強度HHC2時,出現(xiàn)能隙.

(3) 體系的內(nèi)能隨外加磁場的增加逐漸減少；體系的內(nèi)能隨溫度的升高逐漸增大.

(4) 絕對零度下，在高外磁場區(qū)二聚體分子內(nèi)電子自旋平行排列的鐵磁交換作用能JF與二聚體與單基體分子間電子自旋反平行排列的反鐵磁交換作用能JAF之比(-JF/JAF)對體系的內(nèi)能幾乎無影響.但在有限溫度、大外磁場下時，曲線重合出現(xiàn)分裂跡象.說明在該體系存在不同形式能量之間的競爭.

參考文獻：

[1] Shiomi D,Nishizawa M,Sato K,et al.A Prerequisite for Purely Organic Molecule-based Ferrimagnetics:Breakdown of Simple Classical Pictures[J].J.Phys.Chem.B,1997,101(17):3342-3348.

[2] Kikuchi H,Fujii Y,Chiba M,et al.Experimental Observation of the 1/3 Magnetization Plateau in the Diamond-chain Compound Cu3(CO3)2(OH)2[J].Phys.Rev.Lett.,2005,94(22):227201.

[3] Maekawa K,Shiomi D,Ise T,et al.Theoretical Study on Spin Alignments in Ferromagnetic Heterospin Chains with Competing Exchange Interactions:A Generalized Ferrimagnetic System Containing Organic Biradicals in the Singlet Ground State[J].J.Phys.Chem.B,2005,109(19):9299-9304.

[4] Rule K C,Wolter A U B,Sullow S,et al.Nature of the Spin Dynamics and 1/3 Magnetization Plateau in Azurite[J].Phys.Rev.Lett，2008,100(11):117202.

[5] Jeschke H,Opahle I,Kandpal H,et al.Multistep Approach to Microscopic Models for Frustrated Quantum Magnets:The Case of the Natural Mineral Azurite[J].Phys.Rev.Lett,2011,106(21):217201.

[6] Haldane F D M.Nonlinear Field Theory of Large Spin Heisenberg Antiferromagnets:Semiclassically Quantized Solitons of the One-dimensional Easy-axis Néel State[J].Phys.Rev.Lett.,1983,50(15):1153-1156.

[7] Maisinger K,Schollwock U,Brehmer S,et al.Thermodynamics of the (1,1/2) Ferrimagnet in Finite Magnetic Fields[J].Phys.Rev.B,1998,58(10):R5908-R5911.

[8] Batista C D,Ortiz G.Generalized Jordan-Wigner Tr-ansformations[J].Phys.Rev.Lett,2001,86(6):1082-1085.

[9] Fan H Y,Li C.Invariant ‘Eigen-operator′ of the Square of Schr?dinger Operator for Deriving Energy-level Gap[J].Phys.Lett.A,2004,321(2):75-78.

[10] Fan H Y,Wu H.Deriving Vibrating Modes of Some Multiatom Molecules by Virtue of the Invariant Eigenoperator Method[J].Modern Physics Letters B,2005,19(26):1361-1366.

[11] 成泰民,孫樹生,賈維燁,等.自旋為1/2的X-Y模型亞鐵磁棱型鏈的熱力學(xué)性質(zhì)[J].低溫物理學(xué)報,2012,34(6):417-424.

[12] 成泰民，葛崇員，孫樹生，等.自旋為1/2的XY模型亞鐵磁棱型鏈的物性和有序-無序競爭[J].物理學(xué)報,2012,61(18):187502.

[13] 成泰民,王蘊鵬.低溫下自旋為1/2的X-Y亞鐵磁棱型鏈的元激發(fā)譜[J].沈陽化工大學(xué)學(xué)報,2012,26(3):270-273.

[14] 成泰民,王蘊鵬,葛崇員,等.低溫下自旋為1/2的一維亞鐵磁棱型鏈系統(tǒng)的元激發(fā)譜[J].大學(xué)物理,2012,31(1):5-8.

[15] Jacobs I S.Spin-flopping in MnF2by High Magnetic Fields[J].J.Appl.Phys.,1961,32:S61-S62.

[16] Karadamoglou J,Papanicolaou N.Bulk and Surface Spin-flop Transitions in an Antiferromagnetic XYZ Chain[J].Phys.Rev.B,1999,60(13):9477.

[17] Sakai Toru.Spin Flop in One-dimensional Quantum Antiferromagnets[J].Phys.Rev.B,1999,60(9):6238.