Tensor Renormalization-Group study of the surface critical behavior of a frustrated two-layer Ising model

基本信息

• 论文ID: 2510.21269
• 标题: Tensor Renormalization-Group study of the surface critical behavior of a frustrated two-layer Ising model
• 作者: C. Chatelain (Université de Lorraine, CNRS, LPCT, France)
• 分类: cond-mat.stat-mech (统计力学)
• 发表时间: 2025年10月24日 (arXiv预印本)
• 论文链接: https://arxiv.org/abs/2510.21269

摘要

本文研究了由受挫自旋-自旋相互作用耦合的两个二维Ising模型副本的表面临界行为。已知该层间耦合是边缘的(marginal)，体相临界行为属于Ashkin-Teller (AT)普适类，与$J_1$-$J_2$ Ising模型相同。作者通过张量重整化群(Tensor Renormalization-Group, TRG)数值计算研究了表面临界行为，将键权重张量重整化群(Bond-Weight Tensor Renormalization Group, BTRG)算法扩展到带边界的系统。研究发现，AT模型中表面磁化标度维数的二重简并在受挫双层Ising模型(F2LIM)中被解除。这种分裂由F2LIM中单个Ising副本自旋反转下$\mathbb{Z}_2$对称性的破缺所解释。F2LIM的两个不同表面磁化标度维数$x_1^s$和$x_2^s$满足简单的对偶关系：$x_1^s = 1/4x_2^s$。

研究背景与动机

要解决的问题

本文研究受挫$J_1$-$J_2$ Ising模型及其等价模型——受挫双层Ising模型(F2LIM)的表面临界行为，这是继体相临界行为研究之后的自然延伸。

问题的重要性

1. 理论重要性：表面临界行为通常比体相临界行为更丰富。对于同一体相普适类，可能存在多个表面普适类（普通、特殊、非常规和表面相变）。
2. 方法学意义：从共形场论(CFT)角度，表面临界行为由单个Virasoro代数的幺正不可约表示描述，而非体相中的全纯和反全纯两个代数。
3. 物理争议：近年来，TRG计算对$J_1$-$J_2$ Ising模型的一些经典Monte Carlo结果提出质疑，需要更精确的数值方法验证。

现有方法的局限性

1. Monte Carlo方法：在受挫系统中可能存在系统误差
2. 小系统精确对角化：转移矩阵方法仅能处理小尺寸条带
3. 现有TRG算法：主要针对周期边界条件，缺乏处理开放和固定边界条件的系统化方法

研究动机

作者在之前工作8中构造了F2LIM模型，与$J_1$-$J_2$ Ising模型共享相同标度极限，但更适合TRG计算。本文将研究扩展到表面临界行为，特别关注：

• F2LIM与AT模型表面临界行为的差异
• 对称性破缺对表面标度维数的影响
• 通过更大系统尺寸获得更精确的数值结果

核心贡献

1. 算法扩展：将BTRG算法扩展到带边界的系统，可处理开放边界条件(OBC)和固定边界条件(FBC)
2. 临界线精确定位：通过研究畴壁自由能的有限尺寸标度，获得F2LIM更精确的临界线（精度$O(10^{-5})$）
3. 对称性破缺发现：揭示F2LIM中$\mathbb{Z}_2$对称性在边界处的破缺导致表面磁化标度维数简并解除
4. 对偶关系：发现两个表面磁化标度维数满足关系$x_1^s = 1/4x_2^s$，类似于Coulomb气体中的电磁对偶
5. 系统性谱分析：提供了AT模型和F2LIM在不同边界条件下完整的低能谱分析

方法详解

任务定义

输入：二维经典自旋模型的哈密顿量及边界条件 输出：

• 配分函数及自由能
• 转移矩阵的本征值谱
• 通过gap-exponent关系提取的标度维数

约束条件：

• 系统在临界点或临界线上
• 共形不变性成立（二阶相变）

BTRG边界算法架构

1. 张量分解基础

系统配分函数表示为张量网络： $e^{-\beta H[s]} = T_{s_i s_j s_k s_l} T_{s_n s_m s_p s_q} \cdots$

其中每个秩-4张量$T$对应方格子上一个格点的玻尔兹曼权重，四个指标对应格点四角的自旋变量。

2. 边界处理

在左右边界引入秩-2张量$J_s$，表示表面两个自旋之间的耦合玻尔兹曼权重。

3. 奇异值分解(SVD)

对张量进行矩阵重构并分解： $T_{s_1 s_2; s_3 s_4} = \sum_n U_{s_1 s_2; n} \sqrt{\Lambda_n} \times \sqrt{\Lambda_n} (V^+)_{n; s_3 s_4}$

保留$\chi$个最大奇异值以控制计算复杂度。

4. BTRG改进

使用修正的SVD分解： $T_{s_1 s_2; s_3 s_4} = \sum_n U_{s_1 s_2; n} (\Lambda_n)^{1-k/2} (\Lambda_n)^k (\Lambda_n)^{1-k/2} (V^+)_{n; s_3 s_4}$

最优参数为$k = -1/2$，显著提高标度维数估计精度。

5. 迭代粗粒化

• 第一步：垂直和水平分解张量，形成秩-3张量
• 第二步：收缩4张量和8张量格点，生成新的秩-4张量
• 边界处理：左右边界张量$T_L$和$T_R$通过收缩3张量格点生成

经过$n$次迭代，每个张量对应$2^{n+1} \times 2^{n+1}$自旋块的玻尔兹曼权重。

标度维数提取

Gap-exponent关系

从转移矩阵最大本征值$\Lambda_0 > \Lambda_1 > \cdots > \Lambda_n$计算能隙： $E_n - E_0 = \ln(\Lambda_n/\Lambda_0)$

周期边界条件(PBC)： $x_n = \frac{L}{2\pi}(E_n - E_0)$

开放/固定边界条件(OBC/FBC)： $x_n^s = \frac{L}{\pi}(E_n - E_0)$

因子差异源于条带映射到复平面上半部分（OBC/FBC）vs 整个复平面（PBC）。

边缘算符的对数修正

对于边缘算符（如4态Potts模型）： $x_n = \frac{L}{2\pi}(E_n - E_0) + \frac{d_n}{\ln L}$

技术创新点

1. 边界张量网络构造：系统化地处理左右边界张量，保持与体相张量迭代的一致性
2. 多类型边界条件：统一框架处理OBC和多种FBC（相同/混合）
3. 畴壁自由能方法：通过相同/混合FBC的自由能差$\Delta F$精确定位临界点，避免中心荷提取的复杂拟合
4. 高精度计算：保留$\chi=64$个态，系统尺寸可达$L=2^{10}$个自旋

实验设置

模型定义

Ashkin-Teller模型

哈密顿量： $-\beta H = J_2 \sum_{i,j} \sigma^A_{i,j}[\sigma^A_{i+1,j} + \sigma^A_{i,j+1}] + J_2 \sum_{i,j} \sigma^B_{i,j}[\sigma^B_{i+1,j} + \sigma^B_{i,j+1}]$$+ J_1 \sum_{i,j} \sigma^A_{i,j}\sigma^B_{i,j}[\sigma^A_{i+1,j}\sigma^B_{i+1,j} + \sigma^A_{i,j+1}\sigma^B_{i,j+1}]$

• 对称性：二面体群$D_4$
• 临界线（自对偶线）：$e^{-2J_1} = \sinh 2J_2$
• 临界指数参数化：$y \in [0, 3/2]$，通过$\cos(\pi y/2) = \frac{1}{2}[e^{4J_1} - 1]$

受挫双层Ising模型(F2LIM)

哈密顿量： $-\beta H = J_2 \sum_{i,j} \sigma^A_{i,j}[\sigma^A_{i+1,j} + \sigma^A_{i,j+1}] + J_2 \sum_{i,j} \sigma^B_{i,j}[\sigma^B_{i+1,j} + \sigma^B_{i,j+1}]$$+ J_1 \sum_{i,j} \sigma^A_{i,j}[\sigma^B_{i+1,j} - \sigma^B_{i,j+1}] + J_1 \sum_{i,j} \sigma^B_{i,j}[\sigma^A_{i+1,j} - \sigma^A_{i,j+1}]$

关键差异：层间耦合在一个方向铁磁、另一个方向反铁磁（受挫）

评价指标

1. 标度维数精度：与精确值（AT模型）或理论预测的偏差
2. 临界点定位精度：$O(10^{-5})$
3. 能级识别：初级算符、后代算符的正确分类

计算参数

• 保留态数：$\chi = 16, 32, 64$
• 迭代次数：6-10次（对应$L = 2^7$到$2^{11}$）
• BTRG参数：$k = -1/2$
• 本征值求解：ARPACK库的Arnoldi算法

边界条件类型

1. OBC：开放边界
2. 相同FBC：左右边界自旋状态相同
   • 第一组：$(\sigma^A, \sigma^B) = (\uparrow, \uparrow)$两侧
   • 第二组：$(\sigma^A, \sigma^B) = (\uparrow, \downarrow)$两侧
3. 混合FBC：左右边界自旋状态相反
   • 第一组：左$(\uparrow, \uparrow)$，右$(\downarrow, \downarrow)$
   • 第二组：左$(\uparrow, \downarrow)$，右$(\downarrow, \uparrow)$

实验结果

AT模型验证

体相标度维数（PBC）

在临界线上测量的标度维数（图4左）：

• $x_\sigma = 1/8$（二重简并）：磁化密度
• $x_{\sigma\tau} = 1/(8-4y)$：极化密度
• $x_\varepsilon = 1/(2-y) = 4x_{\sigma\tau}$：能量密度
• $x_\sigma + 1$（二重简并）：磁化密度的第一后代

精度分析：

• 近Ising点（$J_1 \approx 0$）：与精确值偏差<1%
• 近4态Potts点（$J_1 \approx 0.2$）：出现系统偏差，因对数修正

表面标度维数（OBC）

图4右显示：

• $x_\sigma^s = 1/(2x_\varepsilon)$（二重简并）：表面磁化
• $x_{\sigma\tau}^s = 1$：表面极化
• $x_\sigma^s + 1$（二重简并）：第一后代

与文献20,21,22的精确预测吻合良好（$J_1 \lesssim 0.15$）。

固定边界条件

相同FBC（图5左）：观察到标度维数$2, 4x_\sigma^s, 3, 4x_\sigma^s+1, 4$

混合FBC（图5右）：仅出现整数标度维数$1, 2, 3, \ldots$（恒等算符的后代）

F2LIM临界线确定

方法创新

通过畴壁自由能$\Delta F$（相同FBC与混合FBC的自由能差）确定临界点：

• 顺磁相：$\Delta F$随迭代次数减小（关联长度有限）
• 铁磁相：$\Delta F$随迭代次数发散（畴壁稳定存在）
• 临界点：不同迭代次数曲线的交叉点（图6）

临界线结果

图7显示$0 < J_1 < J_2$区域的临界线：

• $\chi=32$与$\chi=64$结果差异$O(10^{-4})$
• 与之前工作8（中心荷方法）差异$O(10^{-4})$
• 三临界点位置：$J_1 \approx 0.425$（从$x_{\sigma\tau} = 1/8$估计）
• $J_1 > 0.425$：一阶相变区

F2LIM标度维数分析

体相标度维数（PBC）

图8左显示：

• 最低标度维数接近$x_\sigma = 1/8$
• 关键发现：第二低标度维数与第一个出现系统偏离，打破二重简并
• 第三标度维数类似AT的$x_{\sigma\tau}$
• 拟合$x_{\sigma\tau}$得到参数$y$，进而推导其他标度维数

表面标度维数（OBC）——核心发现

图8右展示最重要结果：

简并解除：出现两个新能级，无AT模型对应

• 两者均从Ising点$x_\sigma^s = 1/2$出发
• 能级分裂随$J_1$快速增大

对偶关系： $x_2^s \approx \frac{1}{4x_1^s}$

数值验证极好（图8右虚线）。这类似Coulomb气体中电磁对偶$x_e(n) = n^2/(2g)$和$x_m(m) = m^2 g/2$（$n=m=1$），但耦合常数$g = 1/(2x_1^s) = x_2^s/2$小于AT模型预期值$2/x_\varepsilon = x_\sigma^s$。

其他能级：

• 第三标度维数与体相$x_\varepsilon = 4x_{\sigma\tau}$吻合
• 在一阶相变区（$J_1 > 0.425$）仍可测量，但gap-exponent关系不再适用

固定边界条件

第一组相同FBC（图9左）：

• 两个最低标度维数：$2x_\sigma^s + 1$和$4x_\sigma^s$（在Ising点简并）
• 第三个：$2x_\sigma^s + 2$

第一组混合FBC（图9右）：

• 两个最低标度维数：$2x_\sigma^s$
• 第三个：$2x_\sigma^s + 1$
• 关键差异：不出现整数标度维数（与AT模型不同）

第二组相同FBC（图10左）：

• 最低：$2x_\varepsilon = 8x_{\sigma\tau}$
• 第二（$J_1 \gtrsim 0.25$后变第三）：$4x_\sigma^s$
• 第三（$J_1 \lesssim 0.25$）：接近$3x_\varepsilon$，更可能是$2x_\varepsilon + 1$（第一后代）

第二组混合FBC（图10右）：

• 最低：$x_\varepsilon$
• 第三：$2x_\varepsilon$或$x_\varepsilon + 1$
• 第二：无明确解释（Ising点为$2x_\sigma^s = 1$但增长过快）

关键实验发现总结

1. 对称性破缺机制：AT模型哈密顿量在单个Ising副本自旋反转$(\sigma^A, \sigma^B) \to (-\sigma^A, \sigma^B)$下不变（4自旋耦合为偶）。F2LIM仅在自旋反转+90°旋转组合下不变，但边界不等价（两边界全铁磁层间耦合，另两边界全反铁磁），导致边界对称性破缺。
2. 能谱结构：F2LIM的OBC谱包含AT模型的所有主要特征，但增加了两个分裂的表面磁化模式。
3. 普适性：尽管微观对称性不同，F2LIM的大部分临界行为（除表面磁化简并）仍与AT模型一致。

相关工作

受挫Ising模型研究

1. Monte Carlo方法：
   • Kalz等4,6，Jin等5,7：支持$J_1$-$J_2$模型属于AT普适类
   • 争议：相变阶数、三临界点位置
2. TRG方法：
   • Li和Yang1，Yoshiyama和Hukushima2，Gangat3：质疑MC结果
   • Chatelain8：构造F2LIM模型，确认AT普适类

表面临界现象

1. 理论框架：
   • Binder和Hohenberg13，Cardy14,16-18：表面临界现象一般理论
   • Pleimling15：综述
2. AT模型表面行为：
   • von Gehlen和Rittenberg20,21：小条带转移矩阵对角化
   • Alcaraz等22：XXZ链Bethe ansatz解，给出精确表面指数
3. TRG边界算法：
   • Iino等29,30：HoTRG和TNR的边界扩展
   • 本文：BTRG的边界扩展

共形场论

1. 体相CFT：双Virasoro代数，$c=1$（AT和F2LIM）
2. 边界CFT：单Virasoro代数，不可约表示确定表面谱
3. Coulomb气体：电磁对偶关系$x_e x_m = 1/4$

结论与讨论

主要结论

1. 算法成就：成功将BTRG扩展到带边界系统，可处理多种边界条件，适用于一般2D经典自旋模型
2. AT模型验证：
   • 精确复现已知体相和表面临界指数
   • 验证了不同FBC下的能谱结构
   • 确认了表面磁化标度维数的二重简并
3. F2LIM新发现：
   • 临界线精确定位（精度$10^{-5}$）
   • 发现表面磁化简并解除，产生两个标度维数$x_1^s$和$x_2^s$
   • 建立对偶关系$x_1^s = 1/(4x_2^s)$
   • 解释为边界$\mathbb{Z}_2$对称性破缺的结果
4. 物理机制：F2LIM的受挫层间耦合（一方向铁磁、另一方向反铁磁）导致边界不等价，破坏单副本自旋反转对称性

局限性

1. 理论解释不完整：
   • 两个表面标度维数$x_1^s, x_2^s$对应的Virasoro不可约表示未确定
   • 关联的表面序参量未识别
   • 对偶关系中的耦合常数$g$与预期值不符，缺乏理论解释
2. 数值精度限制：
   • 保留态数$\chi=64$限制了高激发态精度
   • 近4态Potts点对数修正难以处理
   • 一阶相变区（$J_1 > 0.425$）的结果可靠性存疑
3. 能谱识别：
   • 第二组混合FBC的第二标度维数无法解释
   • 部分高激发态是主算符还是后代算符存在歧义
4. 系统尺寸：虽然达到$L=2^{10}$，但有限尺寸效应仍存在，特别是在临界线端点附近

未来方向

1. 理论分析：
   • 通过共形场论确定$x_1^s, x_2^s$的不可约表示
   • 理解对偶关系的深层起源
   • 构造相应的表面序参量
2. 数值改进：
   • 增大保留态数$\chi$
   • 发展处理对数修正的系统方法
   • 扩展到其他边界条件（如自由、周期混合）
3. 模型推广：
   • 研究其他受挫多层模型
   • 探索三维受挫系统的表面行为
   • 应用于量子受挫模型
4. 实验关联：寻找实现F2LIM的实际物理系统（如磁性薄膜、冷原子系统）

深度评价

优点

1. 方法学创新：
   • BTRG边界扩展系统且实用，图1清晰展示算法流程
   • 畴壁自由能方法定位临界点简洁有效，避免复杂拟合
   • 统一框架处理多种边界条件
2. 数值质量：
   • AT模型验证充分，与精确值吻合优秀
   • 临界线精度达$10^{-5}$，优于前期工作
   • 系统尺寸显著大于转移矩阵对角化方法
3. 物理洞察：
   • 对称性破缺解释清晰且令人信服
   • 对偶关系$x_1^s = 1/(4x_2^s)$是重要发现，提示深层结构
   • 系统比较AT和F2LIM揭示微观差异的宏观效应
4. 呈现清晰：
   • 图表丰富（10个主图）且信息量大
   • 表1-2总结标度维数便于参考
   • 数学推导和物理论证平衡得当

不足

1. 理论深度：
   • 对偶关系缺乏第一性原理推导
   • 耦合常数$g$偏差未充分讨论
   • 未探讨与其他可积模型（如XXZ链）的联系
2. 实验设计：
   • 未系统研究$\chi$依赖性（仅16, 32, 64三个点）
   • 缺少对不同BTRG参数$k$的敏感性分析
   • 一阶相变区数据处理不够谨慎
3. 结果完整性：
   • 部分能级无法解释（如第二组混合FBC）
   • 高激发态分析不足
   • 未讨论表面相变与体相相变的关系
4. 可复现性：
   • 缺少具体代码或详细实现细节
   • 某些拟合过程（如$x_{\sigma\tau}$的抛物线拟合）参数未给出
   • 误差估计不够系统

影响力

1. 方法学贡献：
   • BTRG边界算法将成为研究2D表面临界现象的标准工具
   • 畴壁自由能方法可推广到其他系统
   • 为量子系统（通过量子-经典映射）提供技术基础
2. 物理贡献：
   • 解决AT/F2LIM表面行为的长期问题
   • 对偶关系可能启发新的理论框架
   • 为受挫系统表面工程提供理论指导
3. 实用价值：
   • 算法可直接应用于磁性材料表面研究
   • 为量子计算中的边界效应提供参考
   • 临界线数据可作为其他方法的基准
4. 局限：
   • 限于2D经典系统
   • 理论不完整性限制深入应用
   • 需要高性能计算资源

适用场景

1. 直接适用：
   • 2D经典自旋模型的表面临界行为
   • 受挫磁性系统的相变研究
   • 多层耦合系统的临界现象
2. 潜在应用：
   • 量子自旋链（通过转移矩阵形式）
   • 统计场论模型的数值验证
   • 共形场论边界算符的数值识别
3. 不适用：
   • 强耦合量子系统（需要量子TRG）
   • 三维系统（计算复杂度过高）
   • 一阶相变的定量描述

参考文献（关键文献）

8 C. Chatelain, PRE 111, 024109 (2025) - 作者前期工作，构造F2LIM模型

9 D. Adachi et al., PRB 105, L060402 (2022) - BTRG算法原始文献

22 F.C. Alcaraz et al., J. Phys. A 20, 6397 (1987) - AT模型表面指数精确解

29 S. Iino et al., PRB 100, 035449 (2019) - 边界TRG方法（HoTRG版本）

32 J. L. Cardy, J. Phys. A 17, L385 (1984) - Gap-exponent关系经典文献

总体评价：这是一篇高质量的统计力学数值研究论文，在方法学和物理发现上均有重要贡献。BTRG边界算法的扩展具有广泛应用前景，表面磁化简并解除及对偶关系的发现揭示了受挫系统的新物理。尽管理论解释尚不完整，但数值结果可靠且富有启发性。论文适合发表在物理学顶级期刊，将对凝聚态理论和计算物理领域产生持续影响。