贝里曲率图怎么看？波函数几何、坐标色标、对称约束与霍尔响应积分全解

说明：本文华算科技主要介绍贝里曲率的波函数几何来源、坐标与色标含义、晶体和磁结构对曲率分布的约束，以及电荷霍尔、谷霍尔和非线性霍尔响应中的积分对象。

贝里曲率的定义怎样联系波函数？

贝里曲率描述布洛赫波函数随波矢变化的几何特征。晶体中的电子态写成 ψ_n,k(r)=e^ik·ru_n,k(r)，周期部分 u_n,k包含轨道混合与相对相位。能带 E_n(k)记录本征能量；即使两组能带外观相近，波函数在相邻 k 点之间的变化仍可能不同。沿闭合动量回路比较波函数相位，可得到贝里相位。

对孤立能带，贝里联络定义为 A_n=i⟨u_n,k|∇_ku_n,k⟩，曲率为 Ω_n=∇_k×A_n。二维小回路的相位近似等于Ω_n,z乘以回路在 k_x–k_y平面包围的面积，符号随回路取向而定。波函数整体相位的重新选择会改变联络，孤立带的曲率保持不变。两分量模型还可把本征态的组成映射为球面方向，借此显示波函数的几何变化。

蜂窝晶格模型以 A、B 两个子晶格作为基底，赝自旋方向反映两分量的权重和相对相位。图1中的箭头表示模型向量 d̂(k)，颜色编码赝自旋的 z 分量。方向在倒易空间怎样转动，决定该两带模型的曲率；a、b、c 三种分布具有不同的球面覆盖方式。箭头所指的是波函数组成，样品内的电子运动方向要由速度与外场另行计算。

第一性原理计算从电子哈密顿量取得波函数，再由相邻 k 点波函数的重叠相位或带间矩阵元求曲率。二维材料常画 Ω_z；三维体系有三个分量，图题中的 Ω_xy或 Ω_yx采用张量记号时，还涉及下标次序和符号约定。严格简并处应处理整个简并子空间；任意挑选其中一条带，可能使单带结果随基底选择变化。

石墨烯/RbMnCl₃的 a 面板用黑线表示能量，红色虚线读取右侧曲率纵轴；b 面板才是二维分布。a 面板的单位写为 10⁵ a.u.，b 面板只标 Max、Min，未给出绝对刻度。红色虚线的线色用于识别变量，峰的正负由零轴判断。c、d 面板分别讨论空穴掺杂和层间距对能带的影响，应与曲率面板分开解释。

普通单带曲率若使用 k 的 Å⁻¹坐标，Ω 的单位通常为 Ų；约化坐标按倒格矢给出分数值，积分时带有坐标变换的面积或体积因子。曲率图还可能画单条带、选定带组或占据态求和。这三种对象采用不同权重，跨论文比较峰高时，带序号、化学势、分量和单位都须一致。高对称折线仅采样一组路径，二维色图覆盖的区域可比该路径宽得多。

正负分布怎样受对称性和带间耦合控制？

普通电荷曲率在时间反演作用下，把 k 与 −k 的对应电子态配成相反符号；反演操作则把两处配成同号。平衡态的时间反演体系有相同的对应态占据，电荷霍尔积分相消。两种对称兼具时，占据子空间的电荷曲率迹为零；自旋分辨量采用另一种加权方式。局部正负区域可以很大，而全区净值保持为零。

蜂窝双层紧束缚模型中，层间势差 2V 与子晶格势差 Δ引起不同的峰形变化。增大 V 时，尖峰从能谷中心向外迁移；增大 Δ 时，峰仍靠近能谷中心，但逐渐展宽、降低。图3的明暗区域显示两谷的反号分布。a 面板积分的是按自旋和能谷划分的区域，峰的尾部交叠后，所得数值可以偏离整数。

两谷曲率相反时，同号电荷在外场下可产生反向的横向速度，形成谷流。谷霍尔响应区分 K、K′ 的贡献，电荷霍尔响应合计它们。若散射把两个能谷强烈混合，谷标签的寿命与可测谷信号会改变。图中的谷分区积分讨论低能谷的贡献，二维绝缘体的总陈数则积分整个布里渊区。

Mn₃Sn 的非共线磁结构显示，磁有序与镜面对称的组合会选择允许的霍尔分量。共面自旋构型沿指定镜面路径的 Ω_yx受对称约束；自旋向面外倾斜 3° 后，节线附近能带分裂，带上的颜色出现明显正负变化。横轴仍是 Γ–K–H 路径，纵轴仍是能量；曲率由附在能带上的色标表示。

该计算在倾斜构型中移除自旋轨道耦合后，Ω_yx和相应霍尔分量仍保留显著响应。三角形上非共面自旋产生标量自旋手性 S₁·(S₂×S₃)，改变电子波函数的相位结构；另外的张量分量对 SOC 有不同依赖。小净磁矩、强曲率与有限霍尔响应可以在反铁磁体系中并存，具体分量由磁结构及磁对称操作决定。

带间表达式中，曲率由耦合矩阵元的虚部决定符号，并以能带间隔的平方作为分母。小间隙且耦合非零的区域容易出现尖峰；矩阵元若受对称性限制为零，小能量差也可对应零曲率。WSe₂/VSe₂异质结构把化学势固定为 −0.75 eV，比较外电场改变前后的曲率和分层投影能带，显示同一占据基准下热点可重新分配。

左列绿线和品红线分别标记两种电场条件，正负值由Ω_z的零轴确定，单位为 Bohr²。右侧投影能带的蓝、红表示自旋成分，黑色水平线标记 −0.75 eV 化学势。尖峰沿动量路径很窄，数值积分需在其附近充分采样，并检查插值能带能否再现原始 DFT 的局部间隙；平滑绘图无法给出积分误差。

电荷霍尔、陈数与曲率偶极怎样使用这些分布？

内禀电荷霍尔电导对占据态曲率作带符号求和：σ_xy^int=−(e²/ℏ)∑_n∫[d^dk/(2π)^d]f_n,kΩ_n,z。f 是费米分布，d 是积分维数。二维得到面电导，三维得到体电导；二维薄膜若用 S cm⁻¹报告，还隐含厚度换算。正峰和负峰按各自面积、占据及采样权重相加，峰高排序与电导排序可以不同。

Co₃Sn₂S₂的计算在固定磁化方向下改变 Co 的磁矩，Weyl 点的位置和霍尔电导随磁矩共同变化。图6的 e 面板是积分后的电导，c 面板显示指定能量处的费米面，f 面板标记 Weyl 点的移动。改变磁矩已经改变哈密顿量，所得响应不能仅按一张固定曲率图中热点的数量估算。

c 面板的 High、Low 色标对应费米面图像强度，白色和绿色点标记节点手性；它没有给出 Ω 的色标。Weyl 节点的拓扑电荷由包围节点的闭合曲面上曲率通量计算，三维色图中的单一截面只显示该曲面的一部分。节点远离常用高对称路径时，折线能带也可能遗漏它。

二维绝缘体的总陈数为 C=(1/2π)∑_n∈occ∫_BZΩ_n,zd²k；在这里的符号约定下，σ_xy=−Ce²/h。全区能隙和完整占据子空间保证这一带拓扑积分的整数特征。金属部分占据的霍尔积分通常不量子化。石墨烯/RbMnCl₃模型得到的非平庸间隙原先低于费米能级，空穴掺杂把费米能级调入间隙，才具备讨论量子化电荷霍尔的占据条件。

时间反演保持而反演破缺的金属可出现二阶非线性霍尔响应。相应曲率偶极写为 D_i=∑_n∫[d²k/(2π)²]f_n,k∂_kᵢΩ_n,z，也可经分部积分写成费米面附近的加权量。曲率反号配对仍使线性电荷霍尔相消，但曲率导数及速度权重允许非零偶极；完整三重旋转对称还会把二维偶极约束为零。由偶极计算二阶电导时，响应系数还包含弛豫时间与驱动频率。

扭转双层 WSe₂模型用0.5% 单轴应变破坏三重旋转对称，得到有限 D_x、D_y。图7上方在空穴填充 n_h与层间势差 V_z平面画 D_x，单位为 Å；下方才在变形布里渊区画最高价带的 Ω，单位为 10⁵ Å²。两种色图的横纵坐标、单位和计算对象均不同。

该图计算温度为 1.5 K，扭角为 1.4°；V_z在模型中以 meV 表示层间势差能量，不能把 4.2 meV 直接读取为电场强度。d、e 面板给出来自同一单层能谷的最高价带在 V_z=3、5 meV 时的曲率，时间反演伙伴贡献要在总响应中另行计入。在这一模型里，n_h接近 2、V_z接近 4.2 meV 时，两条价带接触并交换曲率，偶极出现增强和变号。