杂化泛函（PBE0/HSE06）完全解析：交换项构造、作用范围及对能带/局域电子/磁态计算的影响

说明：本文华算科技主要介绍杂化泛函的交换项构造、PBE0 与 HSE06 的作用范围，以及它们对能带、局域电子和磁态计算的影响。

密度泛函理论把多电子体系的基态能量写成电子密度的泛函。交换关联能负责描述同自旋电子的交换效应和电子运动之间的关联；它的精确形式未知，常用的 PBE 以局域密度及其梯度建立近似。杂化泛函仍属于 DFT，改变的是交换关联能的近似形式，没有把整个计算换成 Hartree–Fock 方法。

Hartree–Fock 交换由占据轨道构造。轨道在空间两处的振幅共同进入交换积分，所以这一项带有非局域的轨道依赖；PBE 交换则依据某点附近的密度和梯度计算。把一定比例的前者与半局域交换相加，再保留相应的半局域关联，就得到常见杂化形式。这里的“精确交换”只指给定单行列式下的交换项，关联能仍需近似。

PBE0 的交换关联能可写为 E_xc^PBE0＝0.25E_x^HF＋0.75E_x^PBE＋E_c^PBE。式中的 25% 是交换项的混合权重，不能读成“总能量有四分之一来自 Hartree–Fock”。这个比例来自绝热连接近似的构造，并非每种材料的普适最优值。占据轨道随自洽迭代改变，Fock 交换也随之更新；杂化计算因而通常比 PBE 消耗更多时间和内存。

单电子体系中，理想交换应抵消电子与自身的经典库仑作用。半局域近似对这种抵消并不严格，常使局域电子分布得过于分散。引入部分 Fock 交换可减轻相关误差，但固定比例并不消除所有自相互作用误差，多电子体系中的电子关联也不会因混入交换项而自动得到精确处理。

周期固体还涉及远距离轨道之间的交换积分。为控制长程作用，屏蔽型杂化泛函把库仑核按电子间距分成短程和长程两部分。图中的蓝线对应随距离增大而衰减的短程屏蔽函数，黑线对应未屏蔽核的权重；橙线和绿线来自该论文提出的长程近似及扩展函数，不能当成 HSE06 自身的长程交换曲线。

短程与长程的划分发生在交换相互作用的距离依赖上，不是按晶胞里的“近邻原子”和“远邻原子”各选一套泛函。波函数即使分布在同一晶胞，轨道对积分仍包含不同的电子间距；屏蔽参数改变的是这些距离上的 Fock 交换权重。

PBE0 对全部距离采用 25% Fock 交换。HSE06 保留短程的 25%，长程交换由 PBE 形式承担，常用屏蔽参数为约 0.11 a₀⁻¹。两者的关联项仍采用 PBE。PBE0 的长程 Fock 项持续存在，HSE06 的 Fock 项则随距离增加衰减；短程权重虽然相同，两个交换算符仍会给出不同的轨道能级。

一维有机–无机碘铅晶体含 [Pb₂I₆]²⁻ 链和 DPNDI 有机组分。晶体中的价带顶主要来自铅、碘组成的无机网络，导带底偏向有机配体；两类态的空间分布和轨道性质不同。计算同一结构的能带时，改变交换算符会移动带边能级，并改变部分轨道态之间的相对位置。

在这套已弛豫结构上，论文给出的间接带隙为：PBEsol 0.60 eV，加入自旋轨道耦合后为 0.44 eV；HSE06 为 1.26 eV，PBE0 为 1.94 eV。后两个杂化结果没有同时加入自旋轨道耦合，四个数字并非只改变泛函这一变量。对含 Pb、I 的晶体，比较实验光学吸收时还需区分自旋轨道作用、激子效应和温度条件。

能带图中蓝色占据带与橙色未占据带之间的最小能量差对应计算带隙，零点设在各计算的价带顶。图中的能带色彩不表示电荷增减；轨道归属来自同一研究的投影态密度与带边分布。PBE0 数值高于 HSE06，反映该结构对长程交换处理敏感，不能据此认定所有半导体都遵循同一差值。

杂化泛函给出的是广义 Kohn–Sham 本征值差。实验吸收边还受电子–空穴相互作用、声子、缺陷和样品温度影响；光电子谱测得的能级也采用不同参照。带隙接近某个实验数值，只支持相应计算条件下的数值比较，无法单凭这一项确认带边轨道与载流子局域状态。

在 PBE0 家族中调节精确交换比例 α，价带顶与导带底常呈现不同的移动幅度。BiVO₄、MgO 和 α-SiO₂ 的计算把带边能级对齐到平均静电势；图中上方横轴为 α，下方横轴为 DFT＋U 的 U 值，不能把两条横轴当作同一物理量。红色竖线对应按分段线性条件确定的参数，图中的两种方法各自改变了不同的能量泛函项。

带边移动之外，局域载流子是否形成极化子更依赖电子密度与晶格畸变的共同自洽。BiVO₄ 自陷电子集中在局部 V–O 环境；MgO 中自陷空穴和杂质俘获空穴具有不同局部结构。增加 Fock 交换通常会改变电子分散与局域之间的能量竞争，却没有固定的“α 越大越局域、结果越正确”规则。

相应图像上方是四种极化子的等密度面，下方沿 z 方向给出平面积分后的密度。经过各自参数校准，PBE0(α_k) 与 DFT＋U_k 在这些例子中的密度曲线接近；这指向特定体系与校准条件下的局域分布，两种方法的交换处理仍不相同。局域电子所在原子、邻近键长和形成能应与密度剖面对应。

含过渡金属的晶体可能存在多个自洽磁态，同一几何结构换用泛函后，低能磁态的相对位置也会改变。Ni₂O₂ 的计算将每个收敛态画成一个点，横轴为局域磁矩之和，纵轴是相对各自基态的每原子能量。PBEsol、PBE0、HSE 与 PBEsol＋U 的点群分布不同，其中 PBE0、HSE 与 PBEsol＋U 都出现相近的低能解。

Ni₂O₂ 在多种初始电子态下自洽，形成图中不同局部极小值。若只保留一次自洽得到的磁矩和能带，可能遗漏同一泛函下的其他低能解；图中的颜色记录占据带与基态的差异，不能解释成“颜色越深磁性越强”。磁态比较还需要在同一结构、原胞和自旋约束下核对总能。

混合比例与屏蔽长度规定了交换算符在不同电子间距上的权重。标准 HSE06 的固定参数适用于不少半导体，但介电屏蔽、窄带 d 态和缺陷电荷态在不同材料中差异很大。将 25% 当作可迁移的材料常数，会掩盖带边、缺陷能级与磁态对参数的变化；调节 α 使带隙吻合实验，还会改变局域磁矩或缺陷电荷分布。

不同程序即使写着同一个 PBE0 或 HSE，带隙也可能存在数值差异。跨程序基准图以全势 LAPW 结果为参照，按 PBE、PBE0、HSE 分行列出 QE、VASP 和 FHI-aims 的偏差分布，虚线标出 0.1 eV。横轴是程序间带隙差值的绝对值，不是相对实验的误差，也不是某个泛函固有的误差。

在这类对比中，晶格常数、内部原子坐标、赝势或全电子表示、屏蔽参数的单位换算、k 点及自旋轨道耦合设置都会影响读数。图中的实心柱限定晶格常数差不超过 0.001 Å 的材料。若两个 HSE 结果差了几十分之一电子伏特，应核对它们是否使用同一结构与相同的屏蔽定义，再讨论交换项的物理影响。

对于缺陷或极化子，带隙校准只能限定带边间距；电荷转变能级需核对不同电荷态的总能、有限尺寸修正和局域结构。对于磁性氧化物，能带与磁矩还受自旋初态和可能的强关联影响。杂化泛函改变交换能的计算方式，其数值结果仍属于指定结构、占据态和计算条件下的预测。

以含 Pb、I 的一维晶体为例，若报告 HSE06 的 1.26 eV 带隙，记录中应注明它对应PBEsol 弛豫结构且未加入自旋轨道耦合；同文的 0.44 eV 则来自 PBEsol 加自旋轨道耦合。两组数值对应不同哈密顿量，不能用它们的差值充当纯粹的精确交换修正量。