本章目标
(1) 理解卡达诺夫的块自旋思想:临界点附近,系统在尺度变换下"看起来一样";(2) 对一维伊辛模型严格地做一次重整化变换,看清不动点与重整化流;(3) 对二维伊辛模型做一个近似的重整化变换,得到有限的临界温度与临界指数 ν;(4) 理解不动点的线性化如何给出临界指数,以及普适性从何而来;(5) 了解威尔逊–费希尔不动点与 ε 展开。
§20.1尺度变换的思想
在临界点 ξ=∞,涨落出现在一切尺度上:大的有序区域中套着小的有序区域,小的里面又套着更小的。系统没有特征长度,在统计意义上是自相似的。
卡达诺夫(1966)的想法:把自旋分成边长为 b(以晶格常数为单位)的小块,每块用一个"块自旋"代替(例如取块内多数自旋的方向),块自旋之间有某种有效相互作用 K′。新的晶格常数是原来的 b 倍,所以以新的晶格常数为单位,关联长度变为 ξ′=ξ/b。这样的变换 K→K′=R(K) 称为重整化变换;它保持配分函数(长波长的物理)不变,只是去掉了短尺度的细节。
反复变换时:若 ξ 有限,它不断缩小,系统流向完全无序(K→0)或完全有序(K→∞)的平凡状态;若 ξ=∞,变换后仍然是 ∞。所以临界点对应于重整化变换的不动点 K∗=R(K∗)。
§20.2一维伊辛模型的抽取变换
零场一维链,K=βJ。对每隔一个的自旋求和("抽取"),例如对 s2:
s2=±1∑eKs1s2+Ks2s3=2coshK(s1+s3)≡AeK′s1s3
右边是剩下的自旋之间的新的近邻耦合。分别取 s1=s3 与 s1=−s3:2cosh2K=AeK′,2=Ae−K′。两式相除得 e2K′=cosh2K,即
tanhK′=tanh2K,A=2cosh2K(20.1)
(由 tanhK′=e2K′+1e2K′−1=cosh2K+1cosh2K−1=2cosh2K2sinh2K。)这次变换是严格的:新的链仍然只有近邻相互作用,b=2。
重整化流。因为 tanhK<1,tanhK′<tanhK:每做一次变换,耦合都变弱,最终流向 K=0(无穷高温)。不动点只有两个:K∗=0(稳定,吸引一切有限的 K)与 K∗=∞(T=0,不稳定)。没有有限温度的不动点,所以没有有限温度的相变——与 §18.2 一致。
K′(K),(20.1)K′=K迭代的路径 K (变换以后)—
ξ0/ξ=2n—
图 20.1一维伊辛模型的抽取变换 tanhK′=tanh2K。曲线始终在对角线下方:每变换一次耦合就变弱,关联长度减半,最终流向 K=0。两条线只在 K=0 和 K=∞ 相交,没有有限温度的不动点,所以一维没有相变。
检验:由 (18.9),ξ(K)=−1/lntanhK,所以 ξ(K′)=−1/lntanh2K=ξ(K)/2:恰好是 ξ′=ξ/b。
§20.3二维伊辛模型的近似重整化
在正方晶格上,把格点按国际象棋棋盘那样分成两套,对其中一套("黑格")的自旋求和。剩下的"白格"自旋组成一个转过 45° 的正方晶格,晶格常数为原来的 2 倍(b=2)。被求和的自旋 s0 有 4 个近邻 s1,s2,s3,s4(按顺序绕一圈):
s0∑eKs0(s1+s2+s3+s4)=2coshK(s1+s2+s3+s4)=Aexp[a′(s1s2+s2s3+s3s4+s4s1)+c′(s1s3+s2s4)+d′s1s2s3s4]
右边是能够写出的最一般的对称形式。比较四种情形:(i) 四个自旋相同:ln(2cosh4K)=lnA+4a′+2c′+d′;(ii) 三同一异:ln(2cosh2K)=lnA−d′;(iii) 两两相邻同号(如 ++−−):ln2=lnA−2c′+d′;(iv) 交错(+−+−):ln2=lnA−4a′+2c′+d′。由 (iii)、(iv) 得 a′=c′;(i) 减 (iv) 得 8a′=lncosh4K;(ii)、(iii) 给出 d′=c′−21lncosh2K:
a′=c′=81lncosh4K,d′=81lncosh4K−21lncosh2K(20.2)
在新的晶格中,每一对近邻(如 s1、s2)同时与两个被求和的自旋相邻,所以新的近邻耦合为 K1=2a′=41lncosh4K;次近邻耦合(如 s1、s3)为 K2=c′=81lncosh4K;此外还出现了四自旋耦合 K3=d′。变换后的哈密顿量比原来的复杂——这是重整化变换的普遍特点:耦合常数的空间不断扩大。
近似。略去四自旋耦合,并把次近邻耦合并入近邻耦合(在有序态中两者都使自旋平行,每个自旋有 4 个近邻与 4 个次近邻,贡献相近):
K′≈K1+K2=83lncosh4K(20.3)
不动点。数值求解 K∗=83lncosh4K∗,得 K∗≈0.507,即 kBTc/J≈1.97(严格值 2.269)。K<K∗ 时流向 K=0(无序相),K>K∗ 时流向 K=∞(有序相):近似的变换正确地给出了有限温度的相变。
临界指数。在不动点附近线性化:δK′=λδK,λ=dKdK′K∗=23tanh4K∗≈1.449>1(偏离不动点的扰动被放大)。设 ξ=C∣K−K∗∣−ν。由 ξ(K)=bξ(K′):C∣δK∣−ν=bC∣λδK∣−ν,即 λν=b:
ν=lnλlnb=ln1.449ln2≈0.935(20.4)
严格值是 ν=1。一个如此粗糙的近似给出了相当接近的结果,而且它的思路可以系统地改进(保留更多耦合)。
§20.4不动点、相关变量与普适性
一般地,把所有耦合常数写成一个矢量 K=(K1,K2,…),重整化变换为 K′=R(K)。在不动点 K∗ 附近线性化:δKa′=∑bMabδKb。把 δK 按矩阵 M 的本征矢量展开,各分量(称为标度场 ui)在每次变换中独立地乘以本征值 λi=byi:ui′=byiui。
- yi>0:相关变量,变换中被放大。要到达临界点,必须把它们精确地调为零。伊辛模型有两个:温度 t 与磁场 h。
- yi<0:无关变量,变换中逐渐消失。次近邻耦合、晶格结构、流体中缺少的粒子–空穴对称性等微观细节都属于这一类。
- yi=0:边缘变量,通常导致对数修正。
普适性的来源:微观上不同的系统,只要它们的重整化流趋向同一个不动点,长距离的行为就由这个不动点的相关本征值决定,与那些逐渐消失的无关变量无关。于是它们有相同的临界指数。
标度律的推导。变换后自旋数减少为 1/bd,而总自由能不变,所以每个自旋的自由能奇异部分满足
fs(t,h)=b−dfs(bytt, byhh)(20.5)
b 是任意的,取 b=∣t∣−1/yt,得 fs=∣t∣d/ytfs(±1, h/∣t∣yh/yt)——这正是维多姆形式 (19.21),其中 2−α=d/yt,Δ=yh/yt。又由 ξ′=ξ/b 与 t′=bytt 得 ν=1/yt,于是 2−α=dν:超标度关系 (19.24) 被推导出来了。进一步,m=−∂fs/∂h∝∣t∣(d−yh)/yt,所以 β=(d−yh)/yt;同理 γ=(2yh−d)/yt,δ=yh/(d−yh)。所有临界指数都由两个数 yt、yh 决定。
§20.5威尔逊–费希尔不动点与 ε 展开
威尔逊(1971)在动量空间中对金兹堡–朗道泛函 (19.10) 实施重整化:把波矢在 Λ/b<k<Λ 范围内的涨落积分掉,再把长度与场重新标度。
为什么 4 维是特殊的(量纲分析)。为避免与尺度因子 b 混淆,这里把 (19.10) 中四次项的系数记为 u。在临界点 a=0。把长度放大 b 倍 r=br′,要求梯度项 ∫ddrg∣∇m∣2 的形式不变:它按 bd−2m2 变化,所以场必须按 m=b(2−d)/2m′ 重新标度。于是四次项 ∫ddrum4 按 bd⋅b2(2−d)=b4−d 变化:
yu=4−d(20.6)
d>4 时四次项是无关的:长距离行为由 u=0 的"高斯不动点"决定,平均场指数是严格的。d<4 时四次项是相关的,流向一个新的不动点(威尔逊–费希尔不动点),它的位置正比于 ε=4−d。这与金兹堡判据 (19.19) 得到的上临界维数完全一致。
结果引用
对 n 分量的序参量,ν=21+4(n+8)n+2ε+O(ε2),η=2(n+8)2n+2ε2+O(ε3)(威尔逊与费希尔,1972)。
对三维伊辛模型(n=1,ε=1),一阶结果 ν≈0.583,η≈0.019;计算到高阶并作适当的重求和后得到 ν=0.630。2010 年代发展起来的共形自举方法给出 ν=0.629971(4),是目前最精确的结果之一。
§20.6本章小结
- 临界点是重整化变换的不动点;偏离不动点的相关方向决定临界指数,ν=lnb/lnλ。
- 一维伊辛模型的抽取变换是严格的,它只有 T=0 与 T=∞ 两个不动点;二维的近似变换给出有限的 Tc 与 ν≈0.94。
- 无关变量在变换中消失,这是普适性的来源;自由能的标度关系 (20.5) 推出了维多姆标度与超标度。
- 四次项的标度维数为 4−d:d>4 时平均场严格,d<4 时流向威尔逊–费希尔不动点。
自测题
- 验证 (20.1) 中的 A=2cosh2K,并由此写出每个自旋的自由能满足的递推关系。[答:−βf(K)=21lnA+21[−βf(K′)]]
- 数值求解 (20.3) 的不动点并计算 ν。
- 二维伊辛模型的严格值为 yt=1、yh=15/8。用 §20.4 的公式求 α、β、γ、δ,并与表 §19.6 比较。[答:0、1/8、7/4、15]
这一篇已记为读完。