Tensor network method for solving the Ising model with a magnetic field

基本信息

• 论文ID: 2501.01098
• 标题: Tensor network method for solving the Ising model with a magnetic field
• 作者: Myung-Hoon Chung (College of Science and Technology, Hongik University)
• 分类: cond-mat.stat-mech, quant-ph
• 发表时间: January 3, 2025
• 论文链接: https://arxiv.org/abs/2501.01098

摘要

本文使用张量网络方法研究二维方格晶格上带磁场的伊辛铁磁体和反铁磁体。重点关注规范固定的作用，将配分函数表示为张量网络形式。铁磁体和反铁磁体的张量具有不同的对称性质。配分函数的张量网络被解释为一维量子哈密顿量的多重乘积。使用无限密度矩阵重整化群(iDMRG)来收缩二维张量网络。文章并排展示了伊辛铁磁体和反铁磁体的磁化强度和纠缠熵的数值结果。为了确定温度和磁场参数空间中的临界线，使用一维量子态的半链纠缠熵。纠缠熵精确地指示出反铁磁情况下形成抛物线形状的临界线，而铁磁情况下显示临界点。

研究背景与动机

问题背景

1. 计算复杂性挑战：量子多体系统的计算复杂度随系统大小和维度呈指数增长，需要有效的数值方法处理
2. 张量网络的发展：自1990年代以来，张量网络技术已成为处理多体系统的基础工具，独立地在统计物理、量子多体系统等领域被发现
3. 伊辛模型的重要性：二维伊辛模型是理论研究量子关联和经典统计力学的原型系统

研究动机

1. 方法论创新：探索张量网络方法在处理带磁场的二维伊辛模型中的应用
2. 对称性分析：研究规范固定在张量网络构造中的作用，特别是铁磁和反铁磁情况下的不同对称性
3. 相变检测：利用纠缠熵作为标记来检测量子相变，这是量子信息概念在统计物理中的重要应用

核心贡献

1. 统一框架：同时研究铁磁和反铁磁伊辛模型的张量网络表示，揭示了两者在对称性上的本质差异
2. 规范固定分析：深入分析了规范固定在张量网络构造中的关键作用，展示了不同规范选择的等价性
3. 纠缠熵应用：成功使用半链纠缠熵确定临界线，特别是在反铁磁情况下发现抛物线形状的临界区域
4. 数值方法优化：采用iDMRG方法有效收缩二维张量网络，避免了粗粒化方法在处理缺陷张量时的困难

方法详解

任务定义

研究二维方格晶格上的伊辛模型，包括：

• 输入：温度参数β、磁场强度B、耦合常数J
• 输出：磁化强度、纠缠熵、临界线位置
• 约束：处理铁磁(J>0)和反铁磁(J>0但相互作用为排斥性)两种情况

张量网络构造

配分函数表示

伊辛模型的配分函数：

Z = Σ{σ} e^(-βH({σ}))

其中哈密顿量为：

• 铁磁情况：$H_f = -J\sum_{\langle ij\rangle} \sigma_i\sigma_j - B\sum_i \sigma_i$
• 反铁磁情况：$H_a = J\sum_{\langle ij\rangle} \sigma_i\sigma_j - B\sum_i \sigma_i$

奇异值分解

通过奇异值分解将相互作用项分解：

e^(βJσ_iσ_j) = Σ_m W_σim W_σjm
e^(-βJσ_iσ_j) = Σ_m V_σim W_σjm

其中：

W = (√cosh(βJ)    √sinh(βJ))
    (√cosh(βJ)   -√sinh(βJ))

V = (√cosh(βJ)   -√sinh(βJ))
    (√cosh(βJ)    √sinh(βJ))

四腿张量构造

• 铁磁张量：$T^f_{ruld} = \sum_σ e^{βBσ} W_{σr}W_{σu}W_{σl}W_{σd}$
• 反铁磁张量：$T^a_{ruld} = \sum_σ e^{βBσ} W_{σr}V_{σu}V_{σl}W_{σd}$

关键观察：$T^f_{ruld}$具有完全的指标置换对称性，而$T^a_{ruld}$不具有此对称性。

iDMRG算法实现

传递矩阵方法

将二维统计模型映射为一维量子模型：

... T|A⟩ T|B⟩ T|A⟩ T|B⟩ ... = λ_max × ... |A⟩|B⟩|A⟩|B⟩ ...

其中T为行到行的传递矩阵算符。

矩阵乘积态表示

边界态假设为矩阵乘积态(MPS)，通过iDMRG确定三腿张量A、B和Schmidt系数λ_i。

规范变换

对于反铁磁情况，通过规范变换：

W_σr V_σu V_σl W_σd → V_σr W_σu W_σl V_σd

确保上下收缩过程的一致性。

实验设置

参数设置

• 耦合常数：J = 1（固定）
• 环境键维数：χ = 60
• 温度范围：T ∈ 2.10, 2.40（铁磁），T ∈ 2.18, 2.34（反铁磁）
• 磁场范围：B ∈ -0.10, 0.10（铁磁），B ∈ -1.5, 1.5（反铁磁）

评价指标

1. 磁化强度：通过缺陷张量网络与配分函数的比值计算
2. 半链纠缠熵：$S_h = -\sum_i \lambda_i^2 \log \lambda_i^2$
3. 临界温度对比：与Onsager精确解$T_c = 2.269185314...$比较

基准对比

• Onsager的二维伊辛模型精确解（零磁场情况）
• Yang的铁磁情况精确磁化强度
• Müller-Hartmann和Zittartz的反铁磁临界线结果

实验结果

主要结果

磁化强度

1. 铁磁情况：在零磁场下与Yang的精确结果吻合良好
2. 反铁磁情况：零磁场时总磁化强度在T=2.27以上消失，以下呈现棋盘图案但总磁化为零

纠缠熵分析

1. 铁磁情况：在T≈2.27和零磁场附近出现尖锐峰值，对应临界点
2. 反铁磁情况：显示抛物线区域，所有大熵值几乎相同，形成临界线

临界行为对比

• 铁磁体：展现点状临界行为（临界点）
• 反铁磁体：展现线状临界行为（抛物线形临界线）

数值精度验证

零磁场下的临界温度T_c = 2.269185314与理论值高度一致，验证了数值方法的准确性。

相关工作

张量网络发展

1. 历史发展：从Baxter模型、AKLT态到White的DMRG算法
2. 理论统一：Östlund和Rommer证明DMRG得到的基态就是矩阵乘积态
3. 量子信息视角：从纠缠角度重新理解DMRG，发展出MPO和MPS框架

二维张量网络收缩方法

1. 粗粒化方法：张量重整化群、高阶张量重整化群、张量网络重整化
2. 边界处理方法：变分均匀矩阵乘积态、角传递矩阵重整化群、iDMRG

伊辛模型研究

• Onsager的精确解（1944年）
• 纠缠熵在相变检测中的应用
• 反铁磁伊辛模型的临界行为研究

结论与讨论

主要结论

1. 对称性差异：铁磁和反铁磁伊辛模型的张量网络表示具有本质不同的对称性
2. 纠缠熵标记：纠缠熵成功作为相变的标记，比磁化强度更清晰地指示临界行为
3. 临界行为差异：铁磁体显示点状临界点，反铁磁体显示抛物线形临界线

局限性

1. 键维数限制：χ = 60的键维数可能限制了高磁场区域的数值稳定性
2. 计算复杂度：iDMRG方法在处理大系统时仍面临计算资源限制
3. 适用范围：方法主要适用于可表示为张量网络的配分函数

未来方向

1. 更高键维数：增加键维数以提高数值精度和稳定性
2. 其他模型扩展：将方法推广到其他统计物理模型
3. 量子相变研究：深入探索纠缠熵在检测各类量子相变中的应用

深度评价

优点

1. 方法创新性：巧妙地利用规范固定揭示铁磁和反铁磁模型的本质差异
2. 数值精度高：与已知精确解的对比验证了方法的可靠性
3. 物理洞察深刻：纠缠熵作为相变标记的应用展现了量子信息与统计物理的深度结合
4. 表述清晰：论文逻辑清晰，技术细节充分

不足

1. 参数范围有限：反铁磁情况下大磁场区域的数据不稳定，提示需要更长的模拟时间
2. 理论分析不足：缺乏对为什么纠缠熵能够精确指示临界线的理论解释
3. 计算效率：与其他二维张量网络收缩方法的效率对比分析不充分

影响力

1. 学术贡献：为张量网络方法在统计物理中的应用提供了重要案例
2. 方法价值：所提出的纠缠熵检测相变的方法具有普适性
3. 可复现性：技术细节充分，便于其他研究者复现和扩展

适用场景

1. 二维经典统计模型：适用于可表示为张量网络的各类二维统计物理模型
2. 量子相变研究：为研究量子相变提供了有效的数值工具
3. 材料科学应用：可用于研究磁性材料的相变行为

参考文献

论文引用了26篇重要文献，涵盖了张量网络理论、DMRG方法、伊辛模型精确解等关键领域的经典工作，为研究提供了坚实的理论基础。

总体评价：这是一篇高质量的理论物理论文，成功地将张量网络方法应用于带磁场的二维伊辛模型研究，在方法论和物理洞察方面都有重要贡献。特别是利用纠缠熵检测相变的思路具有重要的学术价值和应用前景。