Theory of local orbital magnetization: local Berry curvature

技术摘要：局部轨道磁化与局部贝里曲率理论

问题陈述
尽管 2005 年已建立了周期性晶体轨道磁化的连贯量子理论（通过轨道矩和贝里曲率定义磁化强度），但对于局部轨道磁化的微观理论在不同几何结构下的适用性仍处于空白。现有的框架（如针对晶体的现代理论）依赖于胞内平均值，从而掩盖了子晶格和单点位的磁性纹理。此外，虽然 Středa 公式将周期性系统中的轨道磁化强度与霍尔电导率联系起来，但目前尚不清楚拓扑信息是如何在有限系统、带状结构或非晶体结构中进行局部编码的。以往的实空间方法（如 Bianco-Resta 局部 Chern 标记）是从全局总磁化强度出发来定义局部磁化强度的，这引入了与表面电流和铰链效应相关的规范歧义。本研究旨在建立一种统一的、规范不变的热力学理论，用以解决周期性晶体、带状结构及有限开边界系统中原子尺度的轨道磁化问题。

方法论
作者开发了一种基于磁场依赖局部态密度 (LDOS) $\rho(r, \varepsilon, B)$ 微扰展开的热力学理论。

• 形式体系： 从局部热势 $\Xi(r, \mu, B) = -\int d\varepsilon F(\varepsilon)\rho(r, \varepsilon, B)$ 出发，局部轨道磁化被定义为零场极限 $M(r, \mu) = -\partial_B \Xi(r)|_{B=0}$。
• 模型： 该理论针对非相互作用无自旋电子紧束缚模型构建，通过 Peierls 替换引入磁场。LDOS 使用规范不变微扰理论进行评估。
• 推导： 核心结果是一个涉及留数格林函数 $G(\varepsilon)$ 和速度算符 $\hat{v}$ 的表达式：
  $$M_{orb}(r, \mu) = \frac{e}{2\pi\hbar} \int d\varepsilon F(\varepsilon) \text{Re}\langle r | G [(\hat{v}G) \times (\hat{v}G)] | r \rangle$$
  该表达式针对三种不同的几何结构进行了评估：具有开边界的有限系统、无限周期晶体（解析至子晶格）以及带状结构（ribbon）。

主要贡献与结果

1. 局部贝里曲率的识别：
   该理论揭示了一个此前未被识别的量：局部贝里曲率 $\Omega_n(r)$（针对有限系统）或 $\Omega_{nk}(r_\alpha)$（针对子晶格）。

• 在有限系统中，局部磁化强度可分解为两项：投影轨道矩 $m_n(r) = |\langle r|n\rangle|^2 m_n$ 和局部贝里曲率项 $F(\varepsilon_n)\Omega_n(r)$。
• 从物理上看，轨道矩描述的是磁场诱导的能级移动（与位置无关），而局部贝里曲率则控制着电子谱权重在实空间的磁场诱导重分布。它满足求和规则 $\sum_r \Omega_n(r) = 0$，表明它仅重新分配权重而不改变状态的总归一化。

2. 跨几何结构的统一描述：
   该形式体系为以下情况提供了连贯的描述：

• 有限系统： 局部磁化强度在对所有位点求和时，可以恢复传统的总轨道磁化强度，此时局部贝里曲率的贡献相互抵消。
• 周期性晶体： 该理论可在子晶格尺度上解析磁化强度。局部贝里曲率分解为传统布洛赫贝里曲率的投影以及一个纯几何贡献 $\Omega_{nk}^{geom}(r_\alpha)$。该几何项在单位晶胞内的子晶格间重新分配谱权重，并满足 $\sum_\alpha \Omega_{nk}^{geom}(r_\alpha) = 0$。
• 带状结构： 这种分为轨道矩和局部贝里曲率的机制同样适用，提供了一个统一的框架。

3. 拓扑编码与斜率量子化：
   作者证明了绝缘隙内的局部轨道磁化强度斜率是量子化的，并且在带状结构和有限几何结构中局部满足 Středa 关系 $\partial_\mu M = \frac{e}{h}C$。

• 该斜率受局部贝里曲率支配。
• 在两能带系统中，斜率包含一个与子晶格无关的拓扑贡献（Chern 数的一半）和一个源于 $\Omega_{nk}^{geom}$ 的与子晶格相关的几何贡献。
• 几何贡献在粒子-空穴对称系统中消失，但在平凡能隙中非常显著，从而能够区分不同类型的轨道磁序（铁磁、反铁磁和亚铁磁）在子晶格层面的表现。

4. 与 Bianco-Resta 方法的比较：
   论文将结果与 Bianco-Resta 局部 Chern 标记进行了对比。

   • 计算效率： 新方法的时间复杂度为 $N^2$，而 Bianco-Resta 方法为 $N^3$。
   • 概念差异： Bianco-Resta 方法从全局总量构建局部磁化强度，这会导致与表面磁化相关的规范歧义。而本研究通过依赖于磁场依赖的 LDOS 直接导出局部响应，提供了一个唯一的、规范不变的局部磁化定义。
   • 物理区别： 两种方法得到的子晶格纹理存在显著差异。关键在于，Bianco-Resta 方法无法解释与几何部分相关的能隙内磁化强度斜率，而这一部分在平凡相中非常显著。

5. 体相拓扑的解剖：
   该理论将拓扑响应分解为能带态和能隙态的贡献。

   • 局部（在体相内部）： 拓扑响应完全源于能带态的局部贝里曲率；能隙态对局部体相拓扑响应没有贡献。
   • 宏观层面： 量子化斜率通过不同的微观机制在不同几何结构下得以恢复。在带状结构中，它由能带态的贝里曲率承载。在有限开边界系统中，贝里曲率贡献在求和后消失，量子化斜率通过能隙态的轨道矩得以恢复。

意义
本文声称确立了局部轨道磁化作为一个真正的微观可观测物理量，将轨道磁化的现代理论从单位晶胞平均量扩展到了原子尺度。通过识别局部贝里曲率，这项工作提供了有限系统中体相拓扑的描述，并揭示了拓扑响应可以独立于边缘态而在体相内部进行局部编码。该框架解析了单点位和子晶格尺度的轨道磁性纹理，提供了一个统一的热力学视角，阐明了不同的微观机制（轨道矩 vs. 贝里曲率重分布）是如何协同产生宏观量子化响应的。这为研究缺乏平移对称性的系统（如分子、非晶态材料、准晶体和莫尔结构）中的轨道磁性纹理开辟了道路。