第十章重整化群理论 10.1 引言 相变理论的核心问题:求临界指数和标度律,阐明相变普适类的根源。 最直接的方法? ---- 求配分函数!因为配分函数包含了统计平衡系统的几乎全部热力学信息。遗憾的是,除了理想气体和少数几个有相互作用的气体,严格求解配分函数十分困难! 朗道平均场理论?---- 可用来求临界指数和标度律,但结果多数情况下与实验结果不符!它只对 的系统适用。而且有两个缺陷(1)它假设自由能是序参量的解析函数,于是可用序参量的幂级数来展开,但这与临界点附近热力学量有奇异性看似是矛盾的;(2)它忽略了涨落,但涨落在发生相变时是很重要的。 标度理论?---- 认为自由能可写为广义齐次函数,是一个形式理论,只可以求出标度律,但不能求出临界指数的值! 数学上的一些严格/近似方法 --- 只针对一些具体的系统求解,不能一般的处理求临界指数的问题。 出路? 考虑系统在相变点附近的对称性!不去求配分函数,而是去寻找保持系统不变的对称变换!(由于系统在这时有标度不变性)普适的临界指数应该对对称性的性质给出描述(普适类的根源)。 这些对称变换的集合形成了一个半群,即重整化群。
10.2 卡丹诺夫变换,块自旋 由于系统在相变点附近有标度不变性,为了考察对称变换的性质,我们将改变观看原系统的(尺度)大小。 唯一相关的尺度是关联长度ξ。系统的临界性质于是不依赖于系统在短距离内的详细细节,只依赖于长程涨落。因此我们可以采用粗粒化(coarse graining)技术把短距离内的细节平均掉。由于系统的最小尺寸为晶格常数 ,我们知道当 为常数时,做这样的尺度变化不会改变系统的性质! 由此我们引入块自旋变换(卡丹诺夫变换):在临界点附近,重要的是大块内的平均自旋而不是格点上的单个自旋,因此可以用只含块自旋的有效哈密顿量来描述系统! 以Ising模型为例,(格点)哈密顿量为 把晶格分为 大小的单元块,定义块自旋 ,它可写为原来格点的平均: 变换后新的哈密顿量为: 这就保持了对称性,只需 显然还有
这么做的一些问题:一般而言系统的哈密顿量可能会发生改变(如果不是仅有最近邻相互作用的话),但我们以后会看到这些多余的耦合不会改变临界指数。 和标度理论的关系: (1)系统自由能:其奇异部分在变换前后有: 为单个自由度的自由能。 令 即得 ,这即是我们前面讨论过的广义齐次函数。 (2)关联函数: 块自旋关联函数可定义为(两点距离为 ): 由 易知 带入到关联函数定义式有: 这即是我们前面讨论过的广义齐次函数。一般地,有
有限尺度标度(finite-size scaling)理论 这里我们考虑一种特别情形,即系统尺度 有限时 (Lꞌ 是系统实际线度大小, 是晶格常数)在临界点附近的性质。这种情形常见于我们对系统进行Monte Carlo模拟时。 这时系统自由能是含L的解析函数,但当 时其某阶导数可能出现奇异性。当L有限时,我们仍然可以把可能出现奇异性的、有自相似性的部分称为自由能的奇异部分。保持系统实际线度Lꞌ 不变,我们进行块自旋变换。和前面讨论类似,单个自由度的自由能的奇异部分变换为: 这里我们也可以认为 热力学极限对应于 的情形, 是相变的临界点。为简单起见,我们下面仅讨论外场h=0的情形。由标度理论,取 我们有磁化率: 其中 是无限大系统在临界点附近的关联长度。类似对比热我们也有: 有限尺度下磁化率χ和比热C随温度的变化: 当 时,关联长度感觉不到系统边界,有限系统的行为和无限系统类似: 当 时,有限系统的行为开始偏离无限系统,有限尺度效应出现; 对更小的 t,有限系统行为不由临界点确定。由于自由能是解析函数, χ和C会出现圆滑的峰。即在 t 的有限变化范围内可近似认为是常数。由此有 当 这时χ的标度形式可能用下式更方便(C的类似):
有限尺度下磁化率χ示意图: 对有限系统,我们可以把磁化率χ或比热C出现极大值的温度认为是有限系统的“临界温度” 。一般而言,这和无限系统的临界温度 并不相等。 假设 的极大值出现在 处,我们有(其中 是常数): 对有限系统我们假设周期性或固定边界条件时,一般有 ,自由边界条件时一般有 (考虑边界对涨落的影响,涨落足够大时破坏有序态产生相变)
10.3 重整化群的定义 仍以Ising模型为例,一般的哈密顿量可写为: 其中α表示某类相互作用(如最近邻、次近邻、与外界作用等)。由上节可知,重整化群变换包括两部分(1)块自旋变换以缩小分辨率;(2)重标变换以恢复到和原有模型一致。 块自旋 和第I块中原有格点的关系为一映射: ,因此可定义函数 对整个系统则可定义: 易知有 用块自旋表示的哈密顿量为: 做变换时我们应保持配分函数不变,因此有 于是重整化变换满足(*):
10.4 重整化群变换的不动点 这时每个自由度的自由能(奇异部分)为: 重整化群变换(*)说明耦合常数 和 间有某种确定的函数关系: 重整化群变换(*)说明耦合常数 和 间有某种确定的函数关系: 这个函数关系易知满足 ,且单位元为 。这个群变换不含逆元(由于粗粒化,细节已经失去了,逆过程是不可能的),因此重整化群是个半群。 重整化群变换实际可由一个生成元 来构造。群的元素为: 这时下列变换实际上与n无关(变换的递推关系): 10.4 重整化群变换的不动点 我们考虑耦合常数K组成的参数空间。如果这个参数空间的一个点在重整化群变换下不变,我们则称这个点是变换的不动点,即这个点满足: 这个点在物理上是非常重要的。注意在这个点关联长度在重整化群变换下有: 要满足这个条件必有 这说明物理上的临界点与非平庸不动点有关。
现在我们来看一下系统在不动点附近的行为。 设 在 附近,我们可以写为: 是一个小量。把变换 在 附近展开到线性项,我们有: 这里 是变换矩阵。由此我们发现 变换矩阵一般不是对称阵,因此其左右本征矢可以不一样。考虑变换矩阵的任一左本征矢 , 而 为其本征值。由定义我们有: 现在我们定义一个标度场 这个标度场经重整化群变换后变为(**): 这个变换对任意的s都成立。另外对任意的s和n有: ,要满足这一条件,必须有: 下面利用(**)式我们讨论本征值的大小: (1)本征值大于一,我们称为关涉本征值,相应本征矢为关涉本征矢,在粗粒化过程中标度场将被不断放大,点也不断远离不动点; (2)本征值小于一,我们称为非关涉本征值,相应本征矢为非关涉本征矢,在粗粒化过程中标度场将被不断缩小,点也不断靠近不动点; (3)本征值等于一,我们称为边缘本征值,相应本征矢为边缘本征矢,它依赖于系统细节。
(1) (2) 非关涉本征矢组成了一个“面”,在这个面上的任一点,经过粗粒化后都将趋于不动点!我们称这个面为临界面。下面是两个例子: 沿临界面上的轴 和 , 在变换下不断接近 , 在这两个方向上是稳定不动点;在 方向上则是不稳定不动点。 (2) 一维情形。 是不稳定不动点。 临界点和不动点的关系: 临界点必与不稳定不动点有关!进行一次重整化群变换,易知关联长度 即系统远离不动点.。 (1) 从物理上看,一个临界系统只能存在两个关涉本征矢,分别为温度和磁场,因此它们称为关涉场。在它们的轴上,粗粒化的方向远离临界面。由此在不动点附近有 (1)当 时: 于是当 时, (2)当 时:当s增加时,开始第四项以后起主要作用, ;但s足够大之后,第二项起主要作用, 又迅速背离 (如右图)。
标度指数: 考虑自由能的奇异部分,由前面的结果我们有: 借助标量场u,我们可把上式改写为: 略去上式里的非关涉本征矢部分即有 这即是我们在标度理论里用到的广义齐次函数。因此我们可以通过求本征值 来获得标度指数! 普适性: 若几个系统的临界哈密顿量都在同一个临界面上且经过重整 化群变换到达同一个不动点,则它们的临界行为属于同一个 普适类。如右图,考虑最近邻和次近邻 的Ising模型: 在临界面(线)上, 和 分别对应于仅含最近邻和仅含 次近邻相互作用的临界汉密顿量,它们的临界行为都由不动 点 决定。
10.5 一维Ising模型(实空间重整化,精确结果) 对一维Ising模型我们可以得到严格结果。粗粒变换可取为:spin 2-spin block,即 由上章的结果,我们有:若令 这里 则配分函数可写为: 对块自旋,我们要求它的矩阵与原先的有相同的形式,即: 这里我们加上了一个参数C,因为解上述矩阵的话有三个方程,两个未知数一般是不够的。 解之得 而粗粒变换可视为参数 的变换。不动点为: (1) 相互作用无穷大/零温,零外场。这个点是不稳定不动点; (2) 相互作用为零/无穷高温度,任意大小外场。这个点是稳定不动点。 临界指数:可通过求解矩阵 的本征值而获得。
一维Ising模型耦合常数和外场在粗粒变换下的性质
10.6 三角晶格上Ising模型的重整化群解(实空间重整化,累积展开) 在很多情况下我们很难获得精确解。在这节里我们用实空间上的累积展开法来获得重整化群近似解。三角晶格上的铁磁Ising系统的哈密顿量和配分函数分别为(N为晶格总数): 首先对晶格进行粗粒变换,我们采用右图所示的变换(实虚) ,这样晶格常数由 变为 ,即重标因子 块自旋 的 取值仍然是 。 的取值由多数法则确定: 这里 表示第I块里的三个格点自旋。 (1) 这四种位形可记为 (2) 可记为 由此,配分函数可写为: 粗粒变换后,新的晶格哈密顿量可通过对 的部分求和获得: 求和只对有确定值 的 进行。 由于块自旋间有相互作用,要实现上面的变换并不简单。为此我们把哈密顿量分为块内自旋间的相互作用 和不同(相邻)块的相互作用及外场的贡献V两部分:
于是有 这里 表示 对 的统计平均。而 为一个自旋块对求和的贡献,M是自旋块的个数,易知 直接计算可得,对 和对 ,均有 然后来求 令V为小量,做展开: 由此可得 最后的展开即为累积展开。如取其中的有限项,可得重整化的哈密顿量的近似表达式。最简单的只保留到 项。此时相邻块间相互作用为右下图所示,于是 因此有 和(***):
下面来求 ,当 时, 当 时, 于是有 带入到(***)式中可得 其中(aa): 于是有 上面第一项为常数,它贡献了系统自由能的解析部分,因此可以略去。 不动点:由(aa)式可发现有三个不动点: 前两个为平庸的稳定不动点,分别对应于温度无穷大和零; 第三个是不稳定不动点,代表非零有限的临界温度。 右图列出了K分量在重整化群变换里的变化趋势。
10.7 动量空间重整化群 临界指数:一阶近似下重整化群变换的变换矩阵为: 这个矩阵的本征值为: 由此可解得 由 和 与六个临界指数的关系即可求得它们的值(见杨展如书304页)。 10.7 动量空间重整化群 我们除了可以在实空间里做粗粒变换,也可以在动量空间里做,而且后者在很多情形下可能更方便。我们先和实空间的情形对比,找到动量空间重整化群的具体步骤。 实空间最基本的动力学变量是自旋 ,为转到动量空间,我们可以做傅里叶变换: 动量空间里的自旋 满足(逆傅里叶变换): 若 为实数,则 一般地,我们有:
动量空间和实空间的粗粒变换 在实空间里,动力学量其实是粗粒化的(coarse grained),它受限于晶格的最小单位 。因此波长不能小于 ,相应地,波矢的最大值为 右图对比了动量空间和实空间的粗粒变换。 可以看出,在实空间里粗粒化相当于在动量空间里改变截断(cutoff)的位置。这可能更为方便。 因此,动量空间重整化群变换分为以下三步: (1)把哈密顿量对应于波矢 的部分积分掉,只剩下 的部分(粗粒化): 配分函数: 这里 (2)把截断从 恢复为Λ (重标晶格长度),即 (3)把 的形式写为与 相同的形式,从而找到耦合常数的递推关系,即相当于找到: 后面的步骤与以前相同。
10.8 高斯模型(动量空间重整化群解) 高斯模型(精确结果): 高斯模型是Ising模型的一个推广,它的自旋取值是连续的且可为任意实数。为使计算结果收敛,我们需要引入权重函数W,它可取为高斯型: 而高斯哈密顿量为: 配分函数为: 由上我们可定义高斯模型的“有效”哈密顿量: 现在把有效哈密顿量换到动量空间。我们利用: 这里 为晶格常数, 是d维超立方体的元胞体积。波矢的截断值为 对动量的求和遍及k空间的第一布里渊区,具体有:
记 我们发现 类似有 于是配分函数可写为: 现在来求 由上,系统的精确的自由能为: 在上式里,第二项求和后与动量无关。第一项显然对所有的k均有 若自由能有奇异性,则奇异性必发生在 处。而由上易知 于是临界温度由下式确定:
重整化群解: 由于 对自由能的奇异性起决定作用,我们可以把前页获得的 的表达式在 附近做展开,并保留到二阶项: 一般的有效哈密顿量可以近似为: 上式最后一项来自于外磁场: 把哈密顿量中的求和换为积分,积分限为 ,并引入记号: 我们有: 这里最后的近似由于奇异性在 我们用一个内切超球代替了第一布里渊区。而配分函数为: 上式对 的积分是泛函积分。我们将以上面两个式子为基础进行重整化群变换。
第一步:粗粒化,积分掉短波部分。 我们把k分为长波区和短波区两部分。长波区为: 短波区为: 而自旋也可分为长波部分和短波部分: 比如令 相应地,哈密顿量和配分函数都可分为两部分: 这里 于是系统的自由能可写为: 由于自由能奇异性只在 处出现,因此 应与自由能的解析部分有关,而对临界行为不产生原则性影响,我们可以将其略去。于是粗粒变换后的哈密顿量变为: 第二步:重标变换。 粗粒变换后自旋和积分区域变了,因此我们做重标: 这样 的取值变回 自旋可以做重标: 稍后确定。
第三步:找到耦合常数的递推关系。 把重标后的变量代到哈密顿量中得: 令 并重新记 为 我们得 与重整化群变换前的哈密顿量对比,我们发现只要令 二者就一致。这就确定了 于是我们发现重整化群的递推关系为: 不动点: 由上面的递推关系易知系统有一个不稳定不动点: 对应的临界温度满足: 临界指数: (对角的)变换矩阵的本征值易知为: 因此 于是可得临界指数为: