Intertwined Orders, Quantum Criticality and Skyrmions in Tunable Topological Bands

基本信息

• 论文ID: 2510.21913
• 标题: Intertwined Orders, Quantum Criticality and Skyrmions in Tunable Topological Bands
• 作者: Xuepeng Wang (Cornell University), Johannes S. Hofmann (Max-Planck-Institut für Physik Komplexer Systeme), Debanjan Chowdhury (Cornell University)
• 分类: cond-mat.str-el (Strongly Correlated Electrons), cond-mat.supr-con (Superconductivity)
• 提交时间: 2025年10月24日
• 论文链接: https://arxiv.org/abs/2510.21913

摘要

本文研究了可调拓扑能带中的交织序、量子临界性和skyrmion激发。作为姊妹篇（arXiv:2507.22971）的延续，作者使用量子蒙特卡罗方法和场论分析研究了Chern反铁磁绝缘体与超导体之间的量子相变。数值结果表明，在临界点附近出现了紧急的SO(5)对称性，这强烈暗示了退禁闭量子（伪）临界性的存在。文章还详细比较了中性集体模式、单电子和复合自旋极化子激发的动量分辨谱函数。

研究背景与动机

核心问题

本研究旨在解决强关联拓扑平带系统中的多个基本物理问题：

1. 交织序的量子相变：Chern反铁磁绝缘体（CAF）与超导体（SC）之间的相变性质
2. Skyrmion激发的稳定性：在非Landau能级极限下，skyrmion激发如何在拓扑平带中存在
3. 退禁闭量子临界性：是否可以在具有拓扑性质的电子模型中实现退禁闭量子临界点（DQCP）

研究重要性

1. 实验相关性：莫尔材料（如扭转双层石墨烯）为研究拓扑、强关联和超导的相互作用提供了高度可调的实验平台
2. 理论挑战：在非投影极限的Chern能带中实现DQCP是一个开放性问题，现有研究主要集中在自旋模型或Dirac点系统
3. 新物理机制：理解skyrmion如何介导超导性，以及拓扑如何影响量子临界行为

现有方法的局限性

1. 投影极限的限制：大多数DQCP研究局限于完全投影的最低Landau能级或活跃Chern能带
2. 符号问题：包含Coulomb排斥的完整模型在一般填充下存在符号问题，限制了QMC方法的应用
3. 平均场近似的不足：Hartree-Fock方法无法准确捕获skyrmion激发的能量学

核心贡献

1. 发现直接连续相变：通过PQMC（投影量子蒙特卡罗）计算，在相互作用各向异性和带宽调控下，观察到CAF与SC之间的单一连续量子相变
2. 紧急SO(5)对称性的数值证据：
   • 在临界点附近，连接CAF和SC序参量的算符Δ_T近似守恒
   • 临界指数（ν ≈ 0.75-0.79, η ≈ 0.20-0.28）偏离传统3D-XY和3D-O(3)普适类
3. 系统的谱函数分析：
   • 比较了DQMC和Bethe-Salpeter方程（基于Hartree-Fock）对中性magnon激发和带电激发的描述
   • 揭示了量子涨落对稳定S≥3/2 skyrmion激发的重要性
4. 低能有效场论：
   • 推导了SO(5)非线性σ模型（NLσM）配合level-1 Wess-Zumino-Witten（WZW）项
   • WZW项的量子化系数反映了底层拓扑，赋予CAF序参量的拓扑skyrmion结构2e电荷
5. 理论与数值的一致性：在JA = 6|JH|时，投影有效理论显示紧急SO(5)对称性，与数值观察的临界行为一致

方法详解

任务定义

研究时间反演对称的自旋全Chern能带模型在半填充（ν=2）附近的量子相变：

• 输入参数：带宽（通过平坦度比F = W/Egap调控）、相互作用各向异性（|JH|/JA）
• 输出：相图、临界指数、序参量关联函数、激发谱
• 约束条件：保持sign-problem-free（JA ≥ 2|JH|），时间反演对称性

模型架构

1. 微观哈密顿量

模型包含动能项和相互作用项：

动能项（每个valley τ = ±）：

H_kin^(τ) = Σ_k ψ†_{k,τ} [B⁰_{k,τ}η₀ + B_{k,τ}·η] σ₀ ψ_{k,τ}

其中：

• ψ†_{k,τ}包含sublattice (η = A,B)和spin (σ =↑,↓)自由度
• 在π-flux方格子上实现，具有最近邻跳跃t、交错次近邻跳跃t₂ = t/√2、第五近邻跳跃t₅（调控平坦度）
• 两个valley通过时间反演关联：T = iτ_yK，实现Chern数C = τ

相互作用项：

H_intra-valley = J_H Σ_{r,τ} S^τ_r · S^τ_r  (J_H < 0, 铁磁)
H_inter-valley = J_A Σ_r S^+_r · S^-_r     (J_A > 0, 反铁磁)

其中S^τ_r是valley-τ的自旋密度算符。

2. 序参量定义

• CAF序参量：O^α_CAF = ψ†_r σ^α τ^z ψ_r (α = 1,2,3)
• SC序参量：O_SC = c†{r+↑}c†{r-↓} - c†{r+↓}c†{r-↑} (自旋单态、valley三重态)
• SO(5)超自旋：N = (ReO_SC, O^α_CAF, ImO_SC)

3. 数值方法：投影量子蒙特卡罗（PQMC）

• 基态制备：|Ψ^Θ_G⟩ = e^{-ΘH}|Ψ_trial⟩
• 虚时演化：Trotter分解，Δτ = Θ/N_Trotter
• 辅助场采样：离散Hubbard-Stratonovich变换，单自旋翻转更新
• 有限尺寸标度：系统尺寸L = 6, 8, 10, 12；投影长度Θ ∝ L（假设动力学指数z = 1）

4. 互补方法：Hartree-Fock + Bethe-Salpeter

• 平均场分解：在CAF对称破缺态上构建
• 准粒子谱：求解HF重整化的单粒子哈密顿量
• 集体激发：通过Bethe-Salpeter方程求解两粒子激发谱
  • Magnon：b†{k,q;+} = c†{k+q,+↑}c_{k,+↓}
  • Spin-polaron：s†+ = c†{ke1;+↑}b†_{ke2,q;+}（电子与magnon束缚态）

技术创新点

1. Sign-problem-free模型设计：
   • 巧妙选择JA ≥ 2|JH|确保QMC可行性
   • 在保持拓扑性质的同时实现大规模数值模拟
2. 双参数相变研究：
   • 相互作用各向异性调控：固定F = 0.01，变化|JH|/JA
   • 带宽调控：固定|JH|/JA = 0.5，变化F
   • 两种路径均观察到类似的临界行为
3. 紧急对称性探测：
   • 通过虚时关联函数χ_ΔT(τ)直接测量紧急守恒荷
   • χ_ΔT(τ) ≈ const在临界点附近，表明H, Δ_T ≈ 0
4. 有效场论推导：
   • 从微观哈密顿量出发，投影到活跃Chern能带
   • 积分掉gapped费米子自由度，推导WZW项
   • 系数量子化反映拓扑不变量π₄(SO(5)/SO(4)) = Z

实验设置

数值参数

• 系统尺寸：L × L晶格，L = 6, 8, 10, 12（部分计算）
• 填充因子：ν = 2（半填充活跃能带）
• 平坦度比：F = 0.01（相互作用调控）；F = 0.01-0.2（带宽调控）
• 相互作用参数：
  • 各向异性调控：JA固定，|JH|/JA = 0.1-0.5
  • 带宽调控：|JH|/JA = 0.5固定
• 温度：T = 0（通过投影实现，Θ = 20L）
• Trotter步数：典型N_Trotter = 200-400

评价指标

1. 关联长度：

ξ_O = 1/(2sin(π/L)) √[S_O(q=0)/S_O(q=(2π/L,0)) - 1]

2. 静态磁化率：

S_O(q) = 1/L² Σ_{r,r'} e^{-i(r-r')·q} ⟨O†(r')O(r)⟩

3. 临界指数：
   • 关联长度指数ν：ξ_O/L = f_ξ(J_O) + L^{-ω}g_ξ(J_O)，J_O = (J-J_c)/J_cL^{1/ν}
   • 反常维度η：S_O(q=0) = L^{1-η}h_{S,ξ}(ξ_O/L)
4. 紧急守恒荷：

χ_ΔT(τ) = Σ_α ⟨Ψ^Θ_G| e^{τH}Δ†_{Tα}e^{-τH}Δ_{Tα} |Ψ^Θ_G⟩ / ⟨Ψ^Θ_G|Ψ^Θ_G⟩

对比方法

1. Hartree-Fock平均场理论：提供准粒子和集体激发的基准
2. Bethe-Salpeter方程：在HF基态上计算复合激发
3. 有效场论预测：SO(5) NLσM + WZW项的理论期望

实验结果

主要结果

1. 相互作用各向异性调控的相变

• 临界点：
  • SC出现：|J^{c1}_H|/J_A = 0.32(1)
  • CAF消失：|J^{c2}_H|/J_A = 0.33(1)
  • 两者在误差范围内重合，支持单一连续相变
• 临界指数（图6）：
  • SC：ν_SC = 0.75(4)，η_SC = 0.20(5)
  • CAF：ν_CAF = 0.79(6)，η_CAF = 0.28(3)
  • 明显偏离3D-XY（η ≈ 0.038）和3D-O(3)（η ≈ 0.037）普适类
• 标度塌缩：ξ_O/L和S_O(q=0)数据在不同系统尺寸下展现良好的标度塌缩（图6b-e）

2. 带宽调控的相变

• 临界点（图7）：
  • SC出现：F^{c1} = 0.166(5)
  • CAF消失：F^{c2} = 0.160(5)
  • 同样观察到直接相变
• 临界指数：
  • SC：ν_SC = 0.65(3)，η_SC = 0.25(2)
  • CAF：ν_CAF = 0.78(6)，η_CAF = 0.27(2)
  • ν_SC显著偏离各向异性调控的结果，可能与电荷涨落有关

3. 紧急SO(5)对称性的证据

• 虚时关联函数（图6f）：
  • 远离临界点（|J_H|/J_A = 0.1, 0.5）：χ_ΔT(τ)呈指数衰减
  • 临界点附近（|J_H|/J_A = 0.34）：χ_ΔT(τ) ≈ const，表明Δ_T近似守恒
  • 这直接证明了SO(5)对称性的紧急性

Hartree-Fock与QMC的比较

1. 准粒子激发（图2a）

• HF和DQMC的单电子谱函数定性一致
• 相互作用诱导的能隙大致相符
• HF略微高估了准粒子速度

2. Magnon激发（图2b）

• Goldstone模式在Γ点无能隙（符合预期）
• HF计算的magnon速度略高于QMC结果
• 整体色散关系定性一致

3. Spin-polaron激发（图2c-d）

关键发现：HF方法失败

• JA = 0情况（图2c）：
  • S = 3/2 spin-polaron（蓝点）与S = 1/2电子准粒子（黑叉）能量简并
  • 最低能带电激发仍是电子准粒子（红线）
• |J_H|/J_A = 0.5情况（图2d）：
  • S = 3/2 spin-polaron明显高于电子准粒子能量
  • 与QMC结果矛盾：QMC观察到S ≥ 3/2 skyrmion是最低能带电激发

物理解释：

• HF + Bethe-Salpeter只包含"baby-skyrmion"（一个电子+一个自旋翻转）
• 量子涨落对电子-magnon有效耦合的重整化至关重要
• 需要超越单圈近似才能正确描述skyrmion能量学

消融实验

虽然论文未明确标注为"消融实验"，但以下分析具有类似作用：

1. 相互作用通道的贡献（表I，方程18-19）：
   • 在J_A = 6|J_H|时，投影哈密顿量显示完美SO(5)对称性
   • 不同相互作用项（U₀, U_αβ, V_μ）的系数在此点统一
2. 形状因子的影响（附录B）：
   • 时间反演对称性导致形状因子分解：F_{k,p,q} = u+(k+q)u+(k)u+(p-q)u+(p)
   • 这确保了紧急对称性从原子极限继承到平带模型
3. WZW项的拓扑起源（方程20）：
   • 系数量子化为3C/(8π²)，其中C是Chern数
   • 直接连接了底层拓扑与有效场论

相关工作

莫尔材料中的强关联物理

1. 扭转双层石墨烯（TBG）：
   • 魔角TBG中的关联绝缘体和超导9-14
   • Chern绝缘体和反常霍尔效应17, 21, 27
   • Kekulé螺旋序和skyrmion晶体43, 46, 47
2. 理论模型：
   • 投影模型：Bultinck等22，Christos等24
   • 重费米子图像：Song和Bernevig50
   • 非局域矩和Mott半金属：Ledwith等51

Skyrmion物理

1. 量子霍尔体系：
   • 经典工作：Sondhi等62，Barrett等63
   • Skyrmion作为拓扑激发64, 65
2. 拓扑平带中的skyrmion：
   • Khalaf和Vishwanath：baby skyrmion57
   • Chatterjee等：DMRG证据59
   • 本工作：非Landau能级极限的稳定性

退禁闭量子临界性（DQCP）

1. 自旋模型：
   • Néel-VBS转变：Senthil等67, 68
   • J-Q模型：Sandvik等75-77, 103, 104
2. 电子模型：
   • QSH-SC转变：Liu等69-71
   • 双层石墨烯：Lee和Sachdev95, 96
3. 数值方法：
   • Fuzzy sphere正则化：Zhou等72, 73，Chen等80, 81
   • 本工作贡献：在非投影Chern能带中实现DQCP

本文相比相关工作的优势

1. 微观现实性：从电子哈密顿量出发，而非纯自旋模型
2. 拓扑非平凡性：活跃能带具有Chern数C = ±1
3. 非投影极限：F = 0.01-0.2远离完全平带极限
4. Sign-problem-free：使QMC大规模计算成为可能
5. 系统比较：数值（PQMC）、半解析（HF+BS）、场论三位一体

结论与讨论

主要结论

1. 直接连续相变：在CAF和SC之间观察到单一连续量子相变，无中间相
   • 相互作用各向异性调控：|J^c_H|/J_A ≈ 0.32-0.33
   • 带宽调控：F^c ≈ 0.16
2. 紧急SO(5)对称性：
   • 临界指数（ν ≈ 0.75-0.79，η ≈ 0.20-0.28）偏离传统普适类
   • 虚时关联函数χ_ΔT(τ) ≈ const在临界点
   • 与SO(5) NLσM + level-1 WZW理论一致
3. 退禁闭（伪）临界性：
   • 数值结果强烈暗示DQCP
   • 可能是"伪临界"行为（由复数CFT导致的弱一阶相变）
   • 需要更大系统尺寸最终确定
4. 量子涨落的关键作用：
   • HF方法无法捕获S ≥ 3/2 skyrmion的能量学
   • DQMC正确描述了skyrmion介导的超导

局限性

1. 系统尺寸：
   • 最大L = 12，可能仍存在显著有限尺寸效应
   • 临界指数ν的估计显示较大漂移
   • 无法完全排除弱一阶相变的可能性
2. 模型简化：
   • 忽略了长程Coulomb相互作用
   • 后者会引入符号问题，限制了研究范围
   • 实际材料中电荷涨落可能重要
3. 相互作用参数限制：
   • Sign-problem-free条件（J_A ≥ 2|J_H|）限制了参数空间
   • J_A = 0的单Chern扇区无法用QMC研究
4. 动力学性质：
   • 主要关注基态和虚时关联函数
   • 实时动力学（如光学电导）未涉及
5. 与实验的联系：
   • 模型参数与实际莫尔材料的对应关系需进一步建立
   • 实验观测临界行为的可行性未讨论

未来方向

1. 更大系统尺寸的计算：
   • 确定相变的真实性质（连续 vs. 弱一阶）
   • 减少临界指数的不确定性
2. 包含长程Coulomb作用：
   • 发展新的数值方法克服符号问题
   • 研究电荷涨落对相变的影响
3. 实验实现路径：
   • 莫尔材料（扭转双层石墨烯、TMD异质结）
   • 冷原子光晶格系统115-119
   • 需要精确控制相互作用和带宽
4. 其他紧急对称性：
   • 探索不同相互作用导致的其他DQCP
   • 如SU(4)、Sp(n)对称性73
5. 动力学响应：
   • 光学电导和低能求和规则52
   • 实时skyrmion动力学
6. 拓扑激发的传输性质：
   • Skyrmion霍尔效应
   • 热电响应

深度评价

优点

1. 方法学创新：
   • 三管齐下：PQMC（非微扰）+ HF+BS（半解析）+ 场论（解析），相互验证
   • Sign-problem-free设计：巧妙的相互作用选择使大规模QMC成为可能
   • 紧急对称性的直接探测：χ_ΔT(τ)方法简洁有效
2. 物理洞察深刻：
   • 揭示HF局限性：清晰展示量子涨落对skyrmion稳定性的重要性（图2c-d）
   • 拓扑-关联-超导链条：从Chern数到WZW项到skyrmion凝聚，逻辑完整
   • 双参数相图：带宽和相互作用调控提供互补视角
3. 数值工作扎实：
   • 系统尺寸L = 6-12，多个参数点
   • 标度塌缩分析规范（图6-7）
   • 独立的crossing-point分析（附录C）验证结果
4. 理论推导严谨：
   • WZW项系数的推导详细（附录B）
   • 形状因子分解的对称性分析清晰
   • SO(5)旋子表示完整（方程16-17）
5. 写作清晰：
   • 结构合理：方法-结果-理论-讨论
   • 图表信息丰富（图1-7，表I）
   • 技术细节放入附录，主文流畅

不足

1. 临界指数的不确定性：
   • ν_SC在两种调控下差异显著（0.65 vs. 0.75）
   • 误差棒较大（如ν_CAF = 0.79(6)）
   • 有限尺寸漂移明显，需要L > 12的数据
2. 伪临界 vs. 真临界的模糊性：
   • 文中多次使用"pseudo-criticality"但未明确定义
   • 与fuzzy sphere结果72的对比不够深入
   • 弱一阶相变的可能性未充分排除
3. 实验连接薄弱：
   • 模型参数（J_H, J_A, t_5）与实际材料的对应关系缺失
   • 未讨论实验观测临界行为的可行性（温度、能量尺度）
   • 与TBG实验现象的定量比较不足
4. 动力学性质缺失：
   • 仅有虚时关联函数，无实时动力学
   • 解析延拓的不确定性未充分讨论（图2a-b的"diffuse curves"）
   • 光学响应、输运性质未涉及
5. 理论预测的可检验性：
   • SO(5)对称性如何在实验中测量？
   • Skyrmion的直接观测证据？
   • 临界指数的实验测量方案？

影响力

1. 对领域的贡献：
   • 概念突破：在非投影拓扑平带中实现DQCP
   • 方法示范：sign-problem-free模型设计范式
   • 理论框架：连接拓扑、强关联、DQCP的统一图景
2. 实用价值：
   • 材料设计：为莫尔材料的相图理解提供指导
   • 数值工具：PQMC + 标度分析的标准流程
   • 场论工具箱：SO(5) + WZW可推广到其他系统
3. 可复现性：
   • 代码开源：使用ALF package121（公开可用）
   • 参数详细：系统尺寸、Trotter步数、投影长度均给出
   • 数据公开性：未明确但图表数据点清晰
4. 后续研究方向：
   • 已有姊妹篇66研究skyrmion介导的超导
   • 可能催生实验验证工作
   • 为其他拓扑平带系统（如twisted TMD）提供模板

适用场景

1. 直接适用：
   • 扭转双层石墨烯（魔角附近）
   • 扭转双层WSe₂（已观测到超导54, 55）
   • ABC堆叠三层石墨烯/hBN49
2. 概念推广：
   • 其他Chern绝缘体平台（拓扑绝缘体表面态）
   • 冷原子实现的合成拓扑能带115-119
   • 重费米子系统中的拓扑Kondo绝缘体
3. 方法论借鉴：
   • 需要研究DQCP的其他量子材料
   • 交织序的量子相变（如铜氧化物高温超导体）
   • 拓扑序与对称性破缺序的竞争
4. 不适用场景：
   • 强无序系统（模型假设平移不变性）
   • 开放系统（方法基于基态投影）
   • 经典相变（量子涨落至关重要）

参考文献（精选关键文献）

方法学基础

• 67, 68 T. Senthil et al., "Deconfined quantum critical points", Science/PRB 2004（DQCP原始提出）
• 98 F. Parisen Toldin et al., "Finite-size scaling of RG-invariant observables", PRB 2015（标度分析方法）
• 121 F. F. Assaad et al., "The ALF project", SciPost Phys. Codebases 2022（QMC代码）

相关理论工作

• 22 N. Bultinck et al., "Ground state and hidden symmetry of magic-angle graphene", PRX 2020
• 56 E. Khalaf et al., "Charged skyrmions and topological origin of superconductivity", Sci. Adv. 2021
• 72 Z. Zhou et al., "SO(5) deconfined phase transition under the fuzzy-sphere microscope", PRX 2024

姊妹篇

• 66 X. Wang et al., "Spin-Polaron Mediated Superconductivity in Doped Chern Antiferromagnets", arXiv:2507.22971

总体评价：这是一篇高质量的凝聚态理论论文，在拓扑平带中实现退禁闭量子临界性方面取得了重要进展。方法学上结合了数值精确性和理论洞察力，物理图景清晰。主要局限在于有限系统尺寸和与实验的距离。对于强关联拓扑系统、莫尔材料和量子临界性领域具有重要参考价值。建议后续工作聚焦于更大尺寸计算和实验可观测量的预测。