一簇金刚石晶格上S 4模型的相变*

尹训昌 刘万芳 马业万 孔祥木 闻军 章礼华

1) (安庆师范大学物理与电气工程学院, 安庆 246011)

2) (曲阜师范大学物理工程学院, 曲阜 273165)

(2018 年7 月7日收到; 2018 年11 月29日收到修改稿)

采用重整化群和累积展开的方法, 研究了一簇金刚石晶格上 S4 模型的相变, 求得了系统的临界点. 结果表明: 当分支数 m=2 和 m＞12 时, 该系统只存在一个Gauss不动点 K∗=b2/2 , u∗2=0 ; 当分支数3≤m≤12时, 该系统不仅有Gauss不动点, 还存在一个Wilson-Fisher不动点, 并且后一个不动点对系统的临界特性产生决定性的影响.

1 引 言

相变是凝聚态物理学中的一个重要研究领域,它是自然界中的一种普遍现象, 如人们熟知的水的气态和液态转变就属于相变. 根据自旋模型的自旋取值不同, 可分为离散型和连续型两种不同的模型. 离散模型只能取一些分离的有限数值, Ising模型和Potts模型是两种最有代表性的离散模型. 而连续模型则允许在整个实数空间取连续值, 如S4模型和Gauss模型. 起初人们主要研究具有平移对称性的平移对称晶格上自旋模型的相变, 随着研究的逐渐深入, 人们对具有高度自相似性的分形晶格产生了浓厚的兴趣. 20世纪80年代, Gefen等[1-4]率先研究了几种简单有规分形上离散模型的相变,求得了一系列有意义的结论. 此后, 分形晶格上的相变问题得到了越来越多的重视, 出现了许多有趣的成果[5-21]. 1999年, 孔祥木等[13]研究了金刚石晶格(2DH)上Gauss模型的相变, 得到了该系统的临界点. 前人的研究发现: 分形晶格上的Ising模型只有当分形维数大于或等于2时才存在有限温度的相变; 分形晶格上的Gauss模型存在有限温度的相变, 临界点数值的变化依赖于分形维数; 分形晶格上S4模型总是存在一个Gauss不动点,Wilson-Fisher不动点是否出现主要取决于分形维数的数值. 本文主要研究具有m个分支数目的金刚石晶格 (mDH)上S4模型的相变, 和 Gauss模型相比较,S4模型具有一个很大的优点, 即增加了一个表示四次方项的参数, 它更加接近自然界中的真实铁磁系统. 因此, 研究S4模型的相变问题可以更好地解释自然界中的相变. 采用实数空间重整化群变换并结合累积展开的方法, 我们研究了一簇金刚石晶格上S4模型的相变.

2 mDH晶格和S4模型

为了便于理解, 我们以两个分支的金刚石晶格(2DH)为例阐述mDH的生成过程. 它的构造过程见图1, 先定义由两个格点和一个键组成的图形作为一个基元(即n=0 级), 再把6个基元依次排列形成一个生成元(即n=1 级), 然后生成元的每一个键再用生成元本身替换, 这样的步骤一直重复下去, 最后得到的晶格为金刚石晶格. 用同样的方法, 可以得到mDH晶格. 该晶格属于非均匀分形, 即格点上的配位数因格点的位置不同而不同.它的分形维数和分岔度分别为df=1+lnm/ln3 和R=∞. 本文研究的金刚石晶格本质上是由迭代产生的具有自相似性的分形晶格, 只不过从外形上看与金刚石的形状相类似, 我们形象地把它简称为金刚石晶格, 它不同于固体物理中的具有周期性的金刚石晶格.

图1 金刚石晶格的构造过程Fig.1. Constructional procedure of diamond lattice.

mDH晶格上S4模型的约化的有效哈密顿量为

其中参数K表示格点自旋si和sj之间的约化相互作用,bi表示格点i位置上的Gauss分布常数, 参数ui代表格点i位置上的四次方项相互作用. 为了求解此种晶格上的相变, 假设Gauss分布常数bi和四次方项相互作用ui取决于格点i的配位数qi, 即满足等式

3 重整化群变换和累积展开

为了方便描述, 选取mDH晶格的生成元来说明重整化群过程(图2). 为了避免混淆, 重整化群变换前后的格点位置用不同的字母和数字来表示,且每个格点上的配位数与格点的位置相关. 通过分析可知2,···,m). 根据(1)式, 生成元的有效哈密顿量写为

其中

图2 mDH晶格的重整化群过程Fig.2. The renormalization group procedure of mDH lattice.

其中大写字母P表示一个与格点自旋无关的重整化群变换常数, 而H′则代表变换后的有效哈密顿量.

为了方便计算, 把有效哈密顿量分解成H0和V两个表达式, 并把表达式V看作H0的微扰进行计算. 根据正则系统配分函数的定义, 容易得到

下面定义部分迹 (PT) 为

重新改写(8)式为

其中

(10)式称为累积平均.

因为V是一个十分小的微扰项, 做级数展开为

把(11)式代入(9)式, 部分迹改写为

消 约 掉 内部 格 点 1,2,···,m和的自旋后, 该系统的有效的哈密顿量重新写为

根据 (12)式得到

其中

利用(5)式和(10)式, 计算得到

其中

应用(5), (10)和(15)式, 可以得到

式中

通过 (13)—(16)式, 经过重整化群变换后, 得到系统的有效的哈密顿量为

为了满足重整化群变换前后系统的有效哈密顿量形势保持一致的特点, 对自旋变量进行重新标度, 定义自旋变量重新标度后系统的有效哈密顿量为

其中

(20)式和(21)式为该系统重整化群变换的递推关系.

4 计算结果

根据上面得到的递推关系, 通过计算得到下面的结果: 1)在mDH晶格的分支数表m=2 和两种情况下, 该晶格上的S4模型只存在一个不动点, 即Gauss不动点2)当mDH晶格的分支数满足 3 ≤m≤12 时, 该晶格上的S4模型不但有Gauss不动点, 而且还发现一个Wilson-Fisher不动点, 不同的分支数m所对应的Wilson-Fisher不动点的数值如表1所列. 根据重整化群的标度变换理论, 可以计算出当分支数满足 3 ≤m≤12 时, Wilson-Fisher不动点附近关联长度的临界指数ν的数值(表1).

表1 mDH晶格上 S4 模型的Wilson-Fisher不动点和关联长度 ν的数值Table 1. Wilson-Fisher fixed point of S4 model on mDH lattices and the value of correlation length ν.

5 结 论

应用实空间重整化群和累积展开的方法, 研究了mDH晶格上S4模型的相变, 得到了如下结论:1)当分支数即分形维数df=1.63和时, 该系统只存在 Gauss不动点其临界点与该晶格上Gauss模型的临界点相同, 表明这两个系统属于同一个普适类; 2)当分支数符合 3 ≤m≤12 , 即分形维数2≤df≤3.26时, 该晶格上的S4模型同时存在Gauss不动点和Wilson-Fisher不动点, 并且后一个不动点对系统的临界特性的影响起主导作用. 进一步研究发现,mDH晶格上S4模型的临界点的个数取决于分形维数, 这与平移对称晶格上的S4模型的相变依赖于空间维数十分类似.