Berry curvature 读图指南：k空间分布、局部热点如何积分成霍尔电导

说明：本文华算科技主要介绍贝里曲率的计算来源、动量空间曲线与色图的读法，以及局部热点如何进入霍尔响应。

周期晶体的电子态由能带序号 n 和波矢 k 标记。沿倒易空间中的封闭小回路移动时，布洛赫波函数会积累一段几何相位。回路缩小后，单位动量面积对应的相位变化就是贝里曲率。二维材料常画垂直于晶面的分量 Ω_n,z(k)；横坐标标记的 k 属于布里渊区，而非样品里的实空间位置。

对一条孤立能带，周期部分波函数 u_n,k 给出贝里联络 A_n(k)=i⟨u_n,k|∇_ku_n,k⟩，二维曲率写成 Ω_n,z=∂A_n,y/∂k_x−∂A_n,x/∂k_y。波函数相位的任选方式会改变联络数值，曲率仍保持规范不变；严格简并附近应处理整个简并子空间，单带相位在该处无法连续指定。

第一性原理计算给出各 k 点的本征能量和波函数；相邻态的重叠相位或速度算符的带间矩阵元可用于构造曲率。能带图画的是 E_n(k)，贝里曲率图画的是 Ω_n,z(k)；同一个高对称点可以在能量图中普通，却在曲率图里形成尖峰，原因还要查该处相邻能带的耦合。

二维 Rashba 模型的能带在外磁场打开 Γ 点简并后，曲率集中到 Γ 附近；晶体场改变轨道多重态时，另一组模型在 Γ 周围形成正负相间的斑块。曲率的红蓝区域表示符号和幅值，图旁的能带曲面表示能量。两类图共用动量坐标，却记录不同物理量。

一维曲线通常沿 K–Γ–K′ 等高对称路径取样：横轴是路径上的波矢，纵轴是某条带或占据态求和后的 Ω_z。二维色图则以 k_x、k_y 为坐标，颜色数值要由该图的色标零点、单位和正负方向读取。单能带曲率常用 Ų；经过电导系数加权的占据态求和图可能采用另一套单位，跨图比较峰值前须核对图注定义。

曲率在三维倒易空间有三个分量。二维薄膜常画 Ω_z，它描述 k_x–k_y 小回路的相位变化；三维论文可能固定一个 k_z 截面，或另画 Ω_x、Ω_y。截面位置和分量下标与颜色同样重要。坐标若以倒格矢为单位给出 0 到 1 的相对值，图上的 0.5 并非 0.5 Å⁻¹。

单层 SnTe 的色图在 X 能谷附近有相邻的正负曲率区。图中的蓝色虚线确定截线位置，旁边的曲线才是沿该线的 Ω 值。铁电极化 P 翻转后，同一 k 点的曲率符号随之反转；曲线峰高随极化幅度变化，不能把色块面积直接当成霍尔电导。

时间反演对称给出 Ω_z(−k)=−Ω_z(k)；反演对称给出 Ω_z(−k)=Ω_z(k)。兼具这两种对称性的非简并带，其普通贝里曲率在全区为零。SnTe 的铁电位移破坏反演对称，正负斑块仍受时间反演配对约束，局部峰值与全区净积分因而是两种不同结果。

MoS₂ 的层数对照把能带、曲率曲线和布里渊区色图排在同一图版。单层与三层具有非零能谷曲率，双层恢复反演对称后曲线接近零；导带谷底在单层位于 K 附近，在双层、三层转向 Q 附近。谷底位置与费米能级穿过哪条子带会改变载流子实际采样的曲率区域。

蓝色曲线合计费米能级以下全部占据态，橙色曲线只合计价带；两线在导带底附近的差值对应新增占据的导带贡献。曲线与零轴的相对位置给出曲率符号，线条的蓝色仅用于区分计算对象。高对称折线只记录布里渊区内的一组点，二维积分使用的是整片倒易空间。

带间表达式中，Ω_n,z 含有速度矩阵元 v_x^nmv_y^mn 的虚部，并由两带能量差的平方加权。相邻能带隔得很近且矩阵元非零时，局部曲率容易变大；峰的符号由矩阵元的相位决定，轨道成分、对称性和自旋轨道耦合均会影响它。

Ni 的 (110) 镜面内，原先受镜面对称保护的节线靠近费米能级。自旋轨道耦合打开其中一段交叉后，占据态曲率在相应位置形成狭窄热点。图版同时给出 U₂PN₂ 的 (100) 镜面与能带截线；节线位置和尖峰位置的空间重合比孤立的一块亮色更能约束热点来源。

Ni 图中的两个坐标是倒格矢方向的相对坐标，只截取镜面里含主要热点的一块区域。该论文把占据态曲率与电导系数合并，色标单位写成 S cm⁻¹ Å³。该数值与 SnTe 的 Ų 峰值采用不同口径，积分维数和归一化方式均已改变。

窄热点对动量网格尤其敏感：粗网格可能从峰两侧穿过，得到较小积分，亦可能把少数大值赋予过宽的权重。论文对 DFT 网格、Wannier 插值、展宽和响应网格分别做了收敛统计。不同网格下的霍尔电导差值量化积分误差；色图外观的平滑程度没有提供收敛数值。

展宽 Γ 增大时，能量差很小的带间响应会分散到更宽的能量范围，曲率尖峰的高度和积分值可能改变；Wannier 能带若在交叉附近偏离原 DFT 能带，热点位置也会偏移。原始能带与插值能带的局部差异，以及积分值随网格和展宽的变化，可用于区分物理热点与数值采样造成的斑点。

内禀反常霍尔电导来自占据态曲率的带符号积分，常写作 σ_xy^int=−e²/ℏ ∑_n∫_BZ[d^dk/(2π)^d] f_n,kΩ_n,z(k)。其中 f_n,k 选出费米能级以下的态。正负区域以符号相加，一处很亮的热点与反号热点相抵后，宏观电导仍可能很小；σ_xy 的整体符号还随张量分量和曲率约定而变。

V 掺杂 Sb₂Te₃ 薄膜的计算把能带、曲率色图和 σ_xy 随化学势的曲线并列。外加电场改变带间隙附近的曲率分布，化学势从 −30 meV 移到 0 和 +42 meV 时，积分纳入的电子态随之改变。同一材料的局部热点形状与测得的霍尔符号之间还隔着能级占据这一步。

若二维绝缘体在整个布里渊区都有能隙，全部占据带的积分 C=(1/2π)∑_n∈occ∫_BZΩ_n,zd²k 给出整数陈数；按上式的符号约定，σ_xy=−Ce²/h。金属的费米能级切过能带，部分占据使这一积分通常不再量子化。陈数对应完整布里渊区和完整占据带，曲率图上一块孤立色斑没有相同含义。

图题若写成自旋贝里曲率，响应算符已经涉及自旋流；若写成贝里曲率偶极，计算对象则含 Ω 对波矢的导数及费米面附近的占据权重。普通电荷曲率、带自旋权重的曲率和曲率偶极分别进入不同输运系数。时间反演体系的普通电荷霍尔积分可以为零，非对称费米面附近的曲率偶极却仍可能有限。

时间反演保持而反演破缺的单层 MoS₂ 中，K 与 K′ 能谷的曲率相反，平衡态电荷霍尔积分相消；谷选择激发或非平衡载流子分布可留下谷霍尔响应。磁性薄膜失去时间反演约束时，正负区域的占据权重可以失衡。霍尔电导曲线的过零点反映费米能级移动或能带重排后的曲率权重变化。

反常霍尔示意图用红、蓝标注的是样品两侧的电势，并非动量空间色图中的 Ω 正负；图内正负曲率积分对应不同的霍尔电导符号。电阻率 ρ_yx 与电导率 σ_xy 的符号约定要分别核对，实验电阻率还含普通霍尔项和散射相关贡献。

在 V 掺杂 Sb_2−xV_xTe₃ 中，ρ_yx(0)–ρ_xx(0) 关系记录不同掺杂薄膜的输运标度；将计算曲率与实验霍尔数据对应时，磁化方向、载流子浓度和电阻率张量的换算口径都需与测量条件一致。