胡俊丽,段美玲
(中北大学 理学院,山西 太原030051)
重正化群(RG)方法在统计物理中的应用为临界现象的研究提供了有效途径,成为探讨具有标度不变和自相似性的复杂系统的一个有力工具[1-4],它不仅从微观理论的角度为唯象的标度假设[5]提供依据,同时还为计算临界指数提供了很好的近似方法.其近似程度的好坏,可以通过它所确定的临界温度Tc及临界指数来判定[6].对一维、二维晶格模型的讨论,许多文献提到过[7-15],但从实际材料来考虑,层状晶格临界性质的讨论(包括考虑一些因素对临界性质的影响)应该意义更大一些,本文着重计算了层状正方晶格临界指数,讨论各向异性的影响;同时考虑能级和温度作为权重因子的影响,推导出了更为准确的重正化群解.
在层状正方晶格点阵模型中,每层正方晶格点阵如图1 所示,取5个格点组成一个Kadanoff集团[16],只考虑最近邻相互作用,在相邻层只考虑对应格点间的相互作用,每层有N 个自旋,共有M 层,则伊辛模型的哈密顿量为
图1 五格点Kadanoff集团Fig.1 Kadanoff clusters with five-lattice point
其中
将Heff分为两部分:Heff=Heff0+V,其中Heff0只包括各集团内部格点的相互作用,V 表示集团与集团间的相互作用,以及在外磁场中的塞曼能,即
采用中心点正负定集团自旋正负,以S′I表示第I 个集团的自旋,则S′I的一个取值对应16个内部状态,分别用16个内部自由度来表示.定义:式中:∑{δI}表示集团自旋分布{S′I}确定后对内部自由度求和,则
通过计算可得一阶项:
由式(4),(5)可得:一阶近似下,性化矩阵为:
其中
讨论: 当l=M (M >1)时,所得临界指数与单层情况下(l=1)完全相同,这是由于各向异性削弱了层间相互作用的影响,导致这样的结果,即临界指数值只与原胞中自旋的个数有关,与空间维数无关.这里若考虑能级对概率的影响,会发现所得结果更接近严格解,方便起见,对二维正方格子伊辛模型结果进行修正.
考虑正方格子伊辛模型中集团出现的概率,对任一点自旋Si,其最近邻4个格点上的自旋为S1~S4,Si=±1,最近邻的相互作用能E =-J∑ij
SiSj,其中J 为交换积分,取J >0,自旋平行为低能态[16].由Si及其近邻格点上的自旋指向,计算Si反转前(不妨设Si=1)的自旋作用能E1=-J(S1+S2+S3+S4),与 反 转 后(此 时Si=-1)的自旋作用能E2=J(S1+S2+S3+S4)比较,设ΔE =E2-E1表示Si反转前后自旋相互作用能的增量.
取Si发生变化的概率P 满足:
当E0时,其中,kB为玻尔兹曼常数.
当ΔE <0 时,P =1.
表1 给出二维情况时Si的所有5种近邻情形及相应的反转概率,其中X 表示自旋指向与Si相同的最近邻格点数.
表1 格点自旋的近邻情形及相应的反转概率(取演化前Si=1)Tab.1 Near condition of lattice point spin and corresponding inverse probability(Si=1before inverse)
对S′I=1(取演化前Si=1,演化后Si=-1),可得16个内部自由度所对应集团的出现概率同 理
同时可得出修正后临界指数值,为比较起见,将修正前后临界指数与二维伊辛模型严格解列表,如表2 所示.
表2 临界指数Tab.2 Critical indices
从表2 看出,引入权重因子后的临界指数有一定的改善,尤其临界指数α,γ,ν更接近严格解.
通过计算层状正方晶格伊辛模型的临界指数,得出的结果与二维情况下的值完全相同,具有相同的临界行为;同时,考虑能级与温度对概率的影响,在重正化变换中引入权重因子,重新推导重正化群解,发现修正后的结果更接近严格解.
[1]Hooyberghs J,Vanderzande C.Real-space renormalization for reaction-diffusion system[J].Phys A:Math Gen,2000,33(5):907-919.
[2]Degenhard A.Non-perturbative real-space renormalization group scheme[J].Phys A:Math Gen,2000,33(5):617-630.
[3]Abdusalam H A.Renormalization group method and Julia sets[J].Chaos.Solitons and Fractals,2001,12(2):423-428.
[4]Kubo R.Statistical mechanics[M].North-Holland Publishing Company Amsterdam,1965:303.
[5]West.Bruce J.Comments on the renormalization group,scaling and measure of complexity[J].Chaos.Solitons and Fractals,2004,20(1):33-44.
[6]Ma S.Modern theory of critical phenomena[M].New York:W.A.Benjamin,1976.
[7]Zhao Xin,Zhou Yunson.To deal with tetragonal ising model by using renormalization group method[J].Journal of Capital Normal University,1997,18(2):37-43.
[8]沈抗存.关于二维伊辛模型严格解中的一个问题[J].湖南大学学报,1990,17(3):81-87.Shen Kangcun.One problem of two-dimensional ising model exact solution[J].Journal of Hunan University,1990,17(3):81-87.(in Chinese)
[9]柴玉珍,李娟娟,李庆士.一维Ising模型的重正化[J].太原理工大学学报,1999,30(1):93-95.Chai Yuzhen,Li Juanjuan,Li Qingshi.Renormalization of one-dimensional Ising model[J].Journal of Taiyuan University of Technology,1999,30(1):93-95(in Chinese)
[10]胡俊丽,许丽萍.自旋玻璃重正化群解的修正[J].中北大学学报(自然科学版),2009,30(6):495-498.Hu Junli,Xu Liping.Revision of renormalization group solution for spin glass[J].Journal of North University of China(Natural Science Edition),2009,30(6):495-498.(in Chinese)
[11]布和.二维六角形晶格伊辛模型的重正化群解[J].大学物理,2001,20(11):12-15.Bu He.The solution of the renormalization group of a planar hexagonal crystal lattice Ising model[J].Col-lege Physics,2001,20(11):12-15.(in Chinese)
[12]章国顺.对《二维六角形晶格伊辛模型的重正化群解》一文的进一步计算[J].大学物理,2006,25(8):24-25.Zhang Guoshun.Deliberation of“The Solution of the Renormalization Group of Planar Hexagonal Crystal Lattice Ising Model”[J].College Physics,2006,25(8):24-25.(in Chinese)
[13]李晓寒.二维Ising模型重正化群方法及数值模拟分析[D].重庆:重庆大学,2005.
[14]王宗笠,李晓寒.Ising模型集团大小对重正化结果影响的分析[J].安徽大学学报,2005,29(4):50-53.Wang Zongli,Li Xiaohan.The analysis of the influence of the group size to the result of renormalization group in Ising model[J].Journal of Anhui University,2005,29(4):50-53.(in Chinese)
[15]龚平,陈耀锋,周详,等.二维晶格上Kadanoff分块及重正化群方法[J].武汉大学学报,2003,49(1):71-74.Gong Ping,Chen Yaofeng,Zhou Xiang,et al.Kadanoff blocks and renormalization group approach to two-dimen sional lattice[J]. Wuhan University,2003,49(1):71-74.(in Chinese)
[16]北京大学物理系《量子统计物理学》编写组.量子统计物理学[M].北京:北京大学出版社,1987.