贝里曲率热点怎样从 VASP 和 Wannier 后处理中找到？

引言

贝里曲率（Berry curvature）常被画成红蓝相间的 k 空间图，但真正要处理的是一套后处理判断：热点是不是来自目标能带，Wannier 插值有没有偏离 VASP 能带，积分量是否随 k 网格收敛。贝里曲率热点不是普通态密度峰，而是能带几何在动量空间中的局域贡献。它通常靠近 SOC 打开的微小能隙、反交叉、近简并能带或费米能级附近的强轨道混合区。

本文把重点放在从 VASP 输出到热点坐标的可执行口径：先得到可靠的 SOC 能带，再完成 Wannier 拟合，然后在密 k 网格上扫描 Ω(k)、AHC 或 Chern 数。图 1 中的 MnAlGe 同时给出 SOC 能带、Berry curvature 分布和 AHC 曲线，适合作为完整后处理链条的起点。

图1 MnAlGe 的 SOC 能带、Berry curvature 三维分布、费米面和反常霍尔电导。热点沿近费米能级节点线邻域集中。 DOI: 10.1002/adma.202006301。

热点先从哪类能带附近找？

Berry curvature 的常用 Kubo 形式可以写成式(1)。它把占据态 n 与未占据态 m 之间的速度矩阵元、能量差和占据数放在同一个求和里。能量分母越小，且两个能带之间的矩阵元越强，Ω_n,z(k) 越容易在某些 k 点附近放大。SOC 打开的小能隙和近费米能级反交叉，是寻找热点的第一类候选区域。因此后处理前不要只看总 DOS，还要回到带结构路径和二维/三维 k 面。

Ω_n,z(k)=−2Im∑_m≠n ⟨u_nk|v_x|u_mk⟩⟨u_mk|v_y|u_nk⟩/(ε_mk−ε_nk)²　　式(1)

式(1)中的 u_nk 是 Bloch 态周期部分，v_x 与 v_y 是速度算符，ε_nk 是能带能量。实际计算不需要手工逐项求和，但这个公式能解释为什么热点经常贴着能带反交叉出现。若 VASP 的 SOC 能带在 K、Γ、M 或低对称路径附近出现小能隙，就应在这些区域附近加密 k 面，而不是只沿标准高对称线出图。

图 2 中的 RuO₂Zn₂F₂ 双层说明另一种情况：AB 与 BA 堆垛改变极化方向后，二维布里渊区内正负 Berry curvature 的空间分布同步改变，AHC 曲线也随化学势变化。这里的热点不是由单个结构图直接给出，而是由 SOC、磁序、堆垛对称性和费米能级共同决定。热点定位必须同时记录能带位置、化学势窗口和磁化方向，否则同一张 Berry curvature 图很难和输运积分对应起来。

图2 RuO₂Zn₂F₂ 双层在不同堆垛与磁态下的 AHC 曲线和二维布里渊区 Berry curvature 分布。正负贡献随结构极化切换而改变。 DOI: 10.1002/advs.74923。

对于补偿磁体或反铁磁体系，还要把 Néel 矢量方向纳入比较。图 3 中 Mn₅Si₃ 的 AHC 来自磁性对称性和 SOC 共同作用，净磁矩很小并不意味着霍尔响应必然消失。若只用普通自旋上下通道解释，会漏掉非相对论磁性劈裂、SOC 混合和对称性允许的横向响应。

在催化相关低维材料中，吸附物、空位或应变可能只改变少数近费米能级能带，却足以改变 Berry curvature 的符号分布。以磁性二维载体负载单原子为例，若金属 d 态和载体 p 态在 K 点附近形成反交叉，吸附前后只比较总磁矩会显得变化不大；真正影响 AHC 的部分可能来自该反交叉处的 SOC 混合增强。相反，若吸附导致费米能级远离所有小能隙区域，热点减弱也可以在没有明显几何畸变的情况下出现。

因此第一轮筛查可以把 VASP 的 SOC 能带分成三类区域：穿过 E_F 的金属带、E_F 附近小带隙反交叉、以及远离 E_F 但与目标化学势窗口相邻的能带组。只有前两类通常需要进入密 k 网格热点扫描，第三类更多用于判断掺杂或门电压移动后是否会引入新贡献。

图3 Mn₅Si₃ 的 SOC 费米面、磁性能带和计算 AHC。Néel 矢量与 SOC 打破的简并关系共同决定霍尔响应。 DOI: 10.1038/s41467-024-48493-w。

VASP 到 Wannier 的输入怎样准备？

VASP 阶段先完成结构优化、普通静态计算和 SOC 静态计算。对磁性体系，ISPIN、MAGMOM、LNONCOLLINEAR、LSORBIT、SAXIS、ISYM 与 NBANDS 要保持可追溯；对含 d/f 元素体系，LDAU、LMAXMIX 与 LASPH 的选择也应沿用到最终能带。Wannier 后处理的可靠性首先取决于 SOC 能带本身是否收敛。如果 VASP 能带在费米能级附近仍受 k 点、展宽或磁态影响，后面的热点图只是在放大前端误差。

Wannier 投影要围绕目标能带选择。金属 Kagome、拓扑半金属和磁性薄膜常需要过渡金属 d 轨道、主族 p 轨道或界面 p-d 混合轨道；如果投影太窄，插值能带可能在费米能级附近错开；投影太宽，则 disentanglement 窗口会变复杂。图4说明，投影窗口应先服务于近费米能级带序、反交叉和能量窗口，而不是只追求整张能带图看起来平滑。

图4 Fe₅GeTe₂ 双层不同堆垛下的 SOC 能带、自旋劈裂和电场/应变响应。Wannier 拟合应优先核对近费米能级反交叉和目标窗口。 DOI: 10.1038/s41467-025-67240-3。

win 参考示例包含哪些关键项？

# ===== Wannier90.win 参考示例 =====

num_wann = 44

num_iter = 300

spinors = true

dis_win_min = -6.0

dis_win_max = 4.0

dis_froz_min = -1.5

dis_froz_max = 1.0

begin projections

Mn:d

Ge:p

end projections

这组参数只作为参考示例。num_wann 由目标能带数量决定，spinors = true 对 SOC 自旋子波函数不可缺少；冻结窗口应覆盖费米能级附近准备解释的反交叉，外窗口再容纳完成拟合所需的相邻能带。不同材料的窗口必须依据 VASP 能带图重新选择。

完成 wannier90 之后，把插值能带与 VASP 的 SOC 能带叠在同一路径上。若热点附近的反交叉、带隙和带序没有对齐，Berry curvature 图不进入物理解读。能带拟合误差最大的地方，往往正是 Berry curvature 对结果最敏感的地方。因此要优先检查热点候选区，而不是只比较远离费米能级的整体现象。

如果 VASP 使用非共线磁性与 SOC，WAVECAR、CHGCAR、wannier90.mmn、wannier90.amn 和 wannier90.eig 应来自同一套静态结果。不同磁化方向不要混用同一份 amn 数据；方向改变会改变自旋子态基底，投影矩阵也随之改变。对含真空层的二维 slab，还要保留同一个晶胞定义和倒格矢方向，避免后处理坐标轴与能带路径标注不一致。

投影质量可从两个层面检查。第一层看插值能带是否覆盖目标能量窗口；第二层看每个热点候选区的反交叉是否在同一 k 位置出现。若只在 Γ-M-K-Γ 路径上拟合得好，但二维 k 面上的局部区域出现错位，热点图仍可能偏离真实 VASP 能带。

密 k 网格怎样把热点坐标找出来？

完成 Wannier Hamiltonian 后，在二维或三维 k 网格上输出 Ω_z(k)、Ω_xy(k) 或 Berry curvature dipole。二维材料可从 200×200×1 起步，金属或节点线体系常需要更密采样；三维体系还要在 k_z 切片中定位贡献层。这里的网格是后处理积分网格，不等同于 VASP 自洽 k 点，但 VASP 能带和 Wannier 拟合必须先给出可靠基底。

σ_xy=−e²/ℏ∑_n∫_BZ d³k/(2π)³ f_n,kΩ_n,z(k)　　式(2)

式(2)说明 AHC 来自所有占据态在布里渊区内的积分。热点强不等于总 AHC 一定大，因为正负 Berry curvature 可能互相抵消。图 2 中 AB/BA 堆垛的正负区域分布不同，图 6 中 AB/BA 堆垛的 AHC 符号也相反；这些现象都要求把局部热点和积分曲线一起看。热点图负责定位贡献区域，AHC 或 Chern 数负责检验积分后的净效应。

实际定位时，可以把每个 k 点的 Ω(k)、能带编号、能量 ε_n,k、相对费米能级距离和坐标写入结果。筛选 |Ω| 最大的一小部分点，再把它们映射回能带图和费米面图。若热点集中在同一个反交叉附近，说明它有明确的能带来源；若热点散落在多个远离费米能级区域，则要检查窗口、投影和化学势范围。

二维体系的 Chern 数可由式(3)给出。绝缘体系若费米能级位于全局带隙内，积分结果应接近整数；金属体系通常讨论 AHC 随化学势的曲线，不应把非整数平台直接称为拓扑不变量。Chern 数检验要求全局带隙和积分收敛同时成立。若带隙很小，k 网格、展宽和能带交叉位置都要单独检查。

C=(1/2π)∑_n∈occ∫_BZΩ_n,z(k)d²k　　式(3)

二维图中的颜色范围也要固定。若每张图都自动按最大值缩放，吸附前后的热点强度无法横向比较。更可靠的处理是保留同一个 Ω(k) 色标，并在附加结果中列出 |Ω| 最大 k 点的坐标、能量和能带编号。这样既能看出热点位置是否迁移，也能看出峰值是否只是因为色标改变而显得更强。

三维体系可按 k_z 切片处理。先固定若干 k_z 面输出 Ω_z(k_x,k_y)，再把贡献最大的切片返回到三维费米面。MnAlGe 中节点线附近的二维 Berry curvature 分布就是典型情况：热点不是沿单一路径出现，而是在三维布里渊区中形成片状或线状高贡献区域。

图5 NbIrTe₄ 双层的 Berry curvature dipole 与动量分辨 Berry curvature。局部能带倾斜和化学势窗口决定主要贡献区域。 DOI: 10.1038/s41467-024-47643-4。

热点结果怎样和物理图像互相验证？

热点定位完成后，应回到三个对象：能带、费米面和对称性。图 1 中的 MnAlGe 把节点线附近 SOC 小能隙、Berry curvature 分布和 AHC 曲线放在同一组结果中，说明热点来自特定近费米能级能带，而不是孤立颜色块。若论文只给一张 Ω(k) 热图，本地结果也应补上对应能带路径或局部放大图。

当化学势移动时，AHC 或 Berry curvature dipole 可能显著改变。图 5 把 NbIrTe₄ 的局部能带、Ω_c(k) 与 D_a 曲线联系起来，说明热点贡献处于窄能量窗口内。对催化或表面吸附体系，吸附物、电荷掺杂或外场可能移动 E_F，这时需要比较吸附前后热点是否仍来自同一类能带。

对比结构必须采用同一后处理口径。以 AB/BA 双层堆垛为例，若两个结构的 AHC 符号相反，应检查它们是否使用相同投影、相同能量零点、相同磁化方向定义和相同密 k 网格。图 6 中符号翻转和电场操控相关，因此局部 Berry curvature 正负区域与输运响应要一起报告。

还有一类情况来自磁性方向或 Néel 矢量改变。若改变 SAXIS 后热点从 K 附近转到 Γ-M 路径附近，需要确认最终磁矩方向与输入方向一致，并检查对应能带反交叉是否同步移动。热点位置改变只有和能带、对称性、磁序输出同时吻合时，才支持物理归因。如果只是颜色图变化而能带拟合没有通过，结果应先回到投影窗口和 SOC 静态计算。

实验对照可从 AHC 符号、化学势依赖和温度趋势切入，但 DFT 后处理不应直接替代实验归因。若计算显示 E_F 附近存在正负贡献快速切换，样品中的载流子浓度、缺陷和应变都会影响观测值。此时局部热点图的作用，是指出哪些能带最值得与 ARPES、输运或磁性测量共同比较。

图6 Fe₅GeTe₂ 双层的可翻转异常输运响应。AB/BA 堆垛的 AHC 符号变化可作为 Berry curvature 积分结果的对照。 DOI: 10.1038/s41467-025-67240-3。

哪些结果适合进入最终判断？

最终结果至少包含五类信息：VASP 的 SOC 能带和磁态设置、Wannier 拟合能带、Ω(k) 热点图、AHC 或 Chern 数积分曲线、热点对应的 k 坐标和能带编号。若研究的是二维磁性材料，还应写清磁化方向和 Chern 数；若研究的是金属反常霍尔响应，应给出 σ_xy(E) 随化学势变化，并说明 E_F 所在位置。

图7 Mn₅Si₃ 的磁化和霍尔电阻测量结果。计算得到的 AHC 符号、温度区间和磁序响应，应回到这类可观测量中核对。 DOI: 10.1038/s41467-024-48493-w。

图 7 把 Mn₅Si₃ 的磁化、霍尔电阻和温度响应放在一起，说明：计算可以指出哪一组能带贡献最大，但结论仍要和可观测的 AHC 符号、温度区间和磁序变化相互核对。

k 网格收敛要直接看积分量。将后处理网格从 100×100×1 加密到 200×200×1、300×300×1，比较 σ_xy 或 C 的变化；若热点尖锐，可在热点区域使用自适应加密或更高均匀网格。只要积分量没有收敛，热点颜色再鲜明也不能支撑定量结论。

能量窗口也要随目标问题调整。讨论费米能级附近输运时，冻结窗口覆盖 E_F 附近反交叉即可；讨论掺杂或门电压时，窗口应覆盖化学势移动范围。若引入吸附、缺陷或应变，先确认结构改变没有产生新的磁性亚稳态，再比较热点迁移。

某个热点来自哪两条 SOC 能带的反交叉，位于哪个 k 区域，对 σ_xy 或 C 的贡献在加密网格下是否稳定。可交付的 Berry curvature 结果，应同时说明局部来源和积分后效应。只有这两部分同时清楚，热点图才不是单纯的可视化结果。

本章要点总结

要点：Berry curvature 热点通常来自 SOC 小能隙、反交叉、近简并能带或强轨道混合区域。

要点：VASP 的 SOC 能带、磁化方向和电子收敛是 Wannier 后处理的前置条件。

要点：Wannier 拟合必须优先检查热点候选区，不能只看整条路径的表面重合。

要点：Ω(k) 定位局部贡献，AHC 或 Chern 数检验积分后的净效应。