丁直，晏世雷

（蘇州大學物理科學與技術學院，江蘇蘇州 215006）

二模磁場中隨機晶場作用混合自旋橫向Ising模型的臨界行為

丁直，晏世雷

（蘇州大學物理科學與技術學院，江蘇蘇州 215006）

在有效場理論和切斷近似的框架內，研究了二模磁場中隨機晶場作用的混合自旋橫向Ising模型的臨界行為.在T-h相圖中，晶場濃度對有序相的影響需要將二模磁場分段考慮；大的晶場和適當?shù)木鰸舛扔欣诙壷厝胂嘧儸F(xiàn)象；但橫場對三臨界點的作用值得更多研究.在T-D相圖中，負晶場比例較大降低有序相，正晶場比例較大增大有序相；強二模磁場和橫場都抑制有序相，但強二模磁場有利于三臨界點的出現(xiàn)，而橫場抑制三臨界點.在T-Ω相圖中，晶場和晶場濃度對有序相范圍有顯著影響.

混合自旋橫向Ising模型；臨界行為；隨機晶場；二模磁場

Jiang等人研究了晶場作用自旋為1的橫向Ising模型，給出適用一般配位數(shù)的二級相變和三臨界點表達式，并考察了晶場對相變的影響[1].晶場作用的橫向Ising模型是橫向Ising模型（TIM）和Blume-Capel模型（BCM）的發(fā)展.該模型中包含了近鄰格點的交換相互作用、晶格排列引起的晶場作用，以及描述量子隧穿效應的橫場作用.隨后該模型受到了其他研究者的重視[2-5].Yan等人將晶場作用單自旋橫向Ising模型發(fā)展為混合自旋系統(tǒng)[6]，一個重要結論是：補償溫度不僅受晶場控制，還受橫場影響.其他研究者的工作已涉及到系統(tǒng)中鍵、晶場和橫場隨機性的影響，獲得很多有意義的結果[7-9].最近，Li等人在不同橫場作用的系統(tǒng)中發(fā)現(xiàn)過飽和磁化行為[10].總的來說，隨機性的引入給系統(tǒng)性質帶來豐富的變化.本文在有效場理論和切斷近似的框架內，研究簡立方晶格情形下二模磁場中隨機晶場作用的混合自旋橫向Ising模型的臨界行為.

1 理論推導

二模磁場中隨機晶場作用的混合自旋橫向Ising模型的哈密頓量為：

這里的格子由子格A和子格B相互嵌套而成，子格A是自旋為1/2的磁性原子，和是z方向和x方向的自旋矩陣；子格B是自旋為1的磁性原子，和是子格B磁性原子z方向和x方向的自旋矩陣.J為最近鄰格點間的交換相互作用，并設J＜0，表明所考慮系統(tǒng)是亞鐵磁耦合，兩個不同子格的磁化方向相反.若J＞0，則系統(tǒng)轉為鐵磁耦合，這里僅考慮亞鐵磁耦合情況.Ω為系統(tǒng)橫場，代表子格A和子格B的量子隧穿效應.Dj是作用在子格B上的晶場，滿足隨機分布

p是隨機晶場濃度.當p=0時，表現(xiàn)為純的負晶場；當p=1時，表現(xiàn)為純的正晶場；0＜p＜1時是混合隨機晶場.hα()

α=i,j是作用于系統(tǒng)的磁場，平行于z軸方向且滿足二模分布：

在有效場理論的框架下，子格A和子格B的子格磁化σ、m具有下列形式：

子格B的四極矩q具有如下形式：

表達式中的＜…＞表示熱力學平均，?=?/?x是微分算符.函數(shù)F(x)、G(x)和U(x)的定義如下：

其中f(x,hi)、g(x,hj,Dj)和u(x,hj,Dj)的表達式是：

這里β=1/kBT，其它涉及量的定義為：

其中系數(shù)為：

因為涉及多自旋相關，要精確求解方程（4）-（6）相當困難，通常做法是引入切斷近似，即：

有i≠j≠…≠k.在切斷近似下，方程（4）-（6）可以表示為：

選取簡立方格子（z=6）作為三維研究對象.聯(lián)立方程（22）-（24），得到子格A磁化σ滿足的自洽方程：

其中

包含的系數(shù)為：

根據(jù)朗道相變理論，二級相變線附近σ很小，因此只需保留線性項，得到相變方程a=1.因為相變方程是（22）-（24）聯(lián)立解的結果，其中包含子格A和B的共同貢獻，因此a=1方程所得到的臨界溫度代表整個亞鐵磁系統(tǒng)的臨界溫度.若b＜0，系統(tǒng)的相變?yōu)槎壪嘧?；若b＞0，系統(tǒng)的相變?yōu)橐患壪嘧?當一條相變線上同時存在二級相變和一級相變時，存在三臨界點.而三臨界點的位置由a=1和b=0所決定.

2 數(shù)值結果和分析

圖1是臨界溫度隨磁場的變化關系，分別對應于（a）D/J=1、Ω/J=0，（b）D/J=2、Ω/J=0，（c）D/J=2、Ω/J=1，不同的曲線對應不同的晶場濃度p.圖1（a）表明，任意晶場濃度下，二級相變線隨二模磁場的增大而單調下降，強二模磁場有效降低臨界溫度.相變線上均存在三臨界點，但大的晶場濃度使三臨界點溫度升高.注意到存在一個臨界二模磁場值h/J=1.84，在h/J＜1.84時，有序相的范圍隨二模磁場的增大而增大，這意味著此時隨機晶場中正晶場的比例越大，越有利于有序相；然而h/J＞1.84時，有序相的范圍隨晶場濃度的增大而減小.即二模磁場分布條件下，隨機晶場濃度的大小對有序相的影響需要分段考慮.從圖1（b）中看到，晶場增大對三臨界點有顯著影響.在0≤p≤0.38范圍內，三臨界點被完全抑制.在0.39≤p≤0.49范圍內，在二級相變線中出現(xiàn)部分一級相變過程，即在一條相變線上存在雙三臨界點，在此范圍內隨著晶場濃度增大，兩個三臨界點分別向高溫和低溫方向發(fā)展.在p≥0.5時，低溫的三臨界點消失，高溫的三臨界點隨晶場濃度的增加繼續(xù)升高.因此在晶場值增大后，較大的負晶場比例（p較?。┠苡行б种迫R界點；適當大小的隨機晶場濃度能導致雙三臨界點的出現(xiàn).與圖1（a）相比，零場時的居里溫度分別向高溫與低溫方向擴展，由于p=0為純負晶場，大的負晶場使得零場時的有序范圍減小，p=1時則反之.在0≤p≤0.05內有序相的范圍擴展到更大的磁場處.另一方面在0.06≤p≤0.46內，觀察到二級重入相變現(xiàn)象，且隨著晶場濃度的增加重入現(xiàn)象逐漸抑制直至消失.因此二級重入相變現(xiàn)象與大的晶場和適當?shù)碾S機濃度有關.注意到基態(tài)時二級相變線在0≤p≤0.05和0.06≤p≤0.38分別聚于橫軸的兩點，說明基態(tài)時在一定的晶場濃度范圍內出現(xiàn)簡并態(tài)，在相圖上即表現(xiàn)為相變線交于同一點.圖1（c）中，在0≤p≤0.32范圍內觀察不到三臨界點，與圖1（b）比較晶場濃度范圍有所減小.以往的工作普遍認為橫場能有效抑制三臨界點，而這里結果表明，在T-h空間中橫場對三臨界點的影響不同，因此對以往的結論需重新認識.零場的居里溫度降低，二級重入相變現(xiàn)象消失，另外零溫時二級相變線的磁場簡并被打開，說明橫場對這些結果有直接影響.綜合圖1，晶場、晶場濃度、二模磁場和橫場這四個因素對居里溫度、三臨界點溫度、二級重入相變、磁場簡并等有重要作用.

圖1 臨界溫度隨磁場的變化關系（曲線上的數(shù)值是隨機晶場濃度p）

圖2 臨界溫度隨晶場的變化關系（曲線上的數(shù)值是隨機晶場濃度p）

圖2是系統(tǒng)在T-D空間的相圖，分別對應于（a）h/J=0.8、Ω/J=0；（b）h/J=1.6、Ω/J=0；（c）h/J=0.8、Ω/J=1，不同的曲線對應不同的隨機晶場濃度p.圖2（a）中，在0≤p≤0.08時存在一個三臨界點，并且隨著晶場濃度的增大而降低；在0.09≤p≤0.12時存在兩個三臨界點，它們隨著晶場濃度的增大而逐漸合并；在0.13≤p≤1時，不存在三臨界點.此時，增大晶場濃度抑制三臨界點.當p≥pC=0.12時，相變線與橫軸沒有交點，這意味著無論晶場多大，低溫時系統(tǒng)始終為有序.p=0.5時居里溫度單調減小，p=0.7時居里溫度先增大后減小，p=0.9時居里溫度單調增大.這是因為負晶場比例較大降低有序相，正晶場比例較大則促進有序相.在圖2（b）中，二模磁場的增大抑制了有序相范圍.在0≤p≤0.73范圍內可觀察到三臨界點，p=0時三臨界點溫度降低，然后三臨界點溫度隨著晶場濃度的增大而升高，晶場濃度對三臨界點的影響與圖2（a）的結果不同.顯然二模磁場的大小對三臨界點有很大影響：弱二模磁場中隨機晶場能抑制三臨界點，強二模磁場中隨機晶場則有利于三臨界點.圖2（c）給出了h/J=0.8、Ω/J=1時的相圖，由于受到橫場的影響，三臨界點完全抑制，有序相范圍減小.圖2清晰地表明，強二模磁場和橫場都抑制有序相，但強二模磁場有利于三臨界點出現(xiàn).

圖3顯示的是系統(tǒng)在T-Ω空間中h/J=0.8時的臨界性質，分別對應于（a）D/J=1和（b）D/J=2，不同的曲線對應不同的晶場濃度p.圖3（a）中，任意晶場濃度下，居里溫度隨著橫場單調降低.隨著晶場濃度增大，系統(tǒng)有序相的范圍增加.圖3（b）與圖3（a）比較，居里溫度分別向低溫和高溫方向發(fā)展，即：隨機濃度小時系統(tǒng)有序相范圍減小，隨機濃度大時有序相范圍增大，這也很清楚地表明負晶場抑制有序相而正晶場促進有序相.

圖3 臨界溫度隨橫場的變化關系（h/J=0.8，曲線上的數(shù)值是隨機晶場濃度p）

3 結論

本文基于有效場理論和切斷近似討論了二模磁場中晶場作用的混合自旋橫向Ising模型的臨界行為.在T-h相圖中，晶場濃度對有序相范圍的影響需要將二模磁場分段考慮；大的晶場和適當?shù)木鰸舛扔欣诙壷厝胂嘧儸F(xiàn)象；橫場可打開基態(tài)的磁場簡并；特別是引入橫場后完全抑制三臨界點的晶場濃度范圍減小，表明T-h空間中橫場對三臨界點的影響是不同的，因此以往的結論需重新認識.在T-D相圖中，負晶場比例較大降低有序相，正晶場比例較大促進有序相；強二模磁場和橫場都抑制有序相，但強二模磁場有利于三臨界點出現(xiàn)，而橫場抑制三臨界點.在T-Ω相圖中，隨著晶場濃度增大，有序相范圍增加；增大晶場后晶場濃度的不同使居里溫度分別向低溫和高溫方向發(fā)展.

[1]Jiang X F,Li J L,Zhong J L,et al.Effect of a crystal field on phase transitions in a spin-1 transverse Ising model[J].Phys Rev B, 1993,827:47-50.

[2]Jiang X F,Zhong J L,Li J L.The random-mixed-bond spin-1 Ising model with a single-ion anisotropy in a transverse field[J].J Phys:Condens Matter,1994,6:523-532.

[3]Deng L L,Yan S L.Phase transition features of the bond-dilution transverse ferromagnetic Ising spin system with random crystal field[J].J Magn Magn Mater,2002,251:138-147.

[4]Htoutou K,Ainane A,Saber M,et al.The site diluted transverse spin-1 Ising model with a longitudinal crystal-field[J].Physica A, 2005,358:184-196.

[5]Akinci U,Yuksel Y,Polat H.Effects of the bond dilution on the phase diagrams of a spin-1 transverse Ising model with crystal field interaction on a honeycomb lattice[J].Physica A,2011,390:541-552.

[6]Yan S L,Yang C Z.The compensation temperature study of the mixed Ising ferrimagnetic spin system with crystal field in a transverse field[J].Z Phys B,1997,103:93-96.

[7]Yan S L,Yang C Z.Effect of the transverse field on the bond-diluted mixed Ising spin system with single-ion anisotropy[J].Phys Rev B,1998,57:3512-3517.

[8]Benayad N,Zerhouni R,Klümper A.Mixed spin transverse Ising model with longitudinal random crystal field interactions[J].EurPhys J B,1998,5:687-695.

[9]Zhang Y N,Yan S L.Thermodynamic properties of random transverse field mixed spin system in the presence of single-ion anisotropy[J].Commun Theor Phys,2003,40:375-380.

[10]Li D R,Yan S L,Zhang Y F.Compensation behaviors and magnetization processes of different transverse fields mixed Blume–Capel model in a magnetic field[J].Sol Stat Commun,2010,150:2186-2189.

The Critical Behaviors of Mixed Spin Transverse Ising Model with the Random Crystal Field in a Bimodal Magnetic Field

DING Zhi,YAN Shi-lei
(School of Physical Science and Technology,Soochow University,Suzhou 215006,China)

Within the framework of the effective field theory(EFT)and cutting approximation,the critical behav?iors of mixed spin transverse Ising model with the random crystal field in a bimodal magnetic field are investigat?ed.InT-hspace,the effect from crystal field concentration to ordered phase range needs to be considered at different areas of bimodal magnetic field.A big crystal field and appropriate random concentration benefit the second phase reentrance phenomenon.The result shows it needs more discussion about the influence of trans?verse field to the tricritical point.InT-Dspace,a high proportion of negative crystal field leads to small or?dered phase range while high proportion of positive crystal field brings large ordered phase range.Both a big bi?modal magnetic field and transverse field can depress ordered phase.But a big bimodal magnetic field can in?crease tricritical point temperature while the transverse field suppresses tricritical point.InT-Dspace,crystal field and crystal field concentration have a marked impact on ordered phase.

critical behaviors;mixed spin transverse Ising model;random crystal field;bimodal magnetic field

O469

A

1008-2794（2013）04-0029-06

2013-05-15

江蘇省高校優(yōu)勢學科建設工程資助項目（YK40000112）

晏世雷，教授，研究方向：凝聚態(tài)理論，E-mail：sly@szda.edu.cn.