2. 量子材料与电子结构基础

要理解量子材料的电子结构计算,不能只把软件当作“黑箱”:输入晶体结构,输出一张能带图。真正决定计算结果物理含义的,是五个层层递进的核心对象——波函数、哈密顿量、能带、态密度、费米面。它们不是并列的名词,而是一条从“微观量子态”到“宏观可测量”的逻辑链。本节沿着这条链,说明每个对象是什么、在计算中如何出现、以及它对应真实材料中的什么物理意义。

一、波函数:电子的量子态描述

在量子力学中,电子的状态由波函数 ψ(r) 描述,|ψ(r)|² 给出在位置 r 找到电子的概率密度。晶体是周期性结构,因此电子波函数必须满足布洛赫定理:

ψ_nk(r) = e^{ik·r} · u_nk(r)

其中:
  k   —— 布洛赫波矢,属于第一布里渊区
  n   —— 能带指标
  u_nk(r) —— 与晶格同周期的周期函数,u_nk(r+R)=u_nk(r)

这个形式的关键含义是:电子不是局域在某个原子上,而是以“调幅平面波”的形式在整个晶体中延展。k 是标记电子量子态的量子数,n 是同一 k 下不同能级的编号。计算中我们并不直接存储连续的 ψ(r),而是把它展开成基组(平面波、原子轨道或数值轨道)的线性组合,系数就是实际求解的未知量。因此,波函数是电子结构计算的“解”本身,能带、态密度、费米面都是从它派生出来的。

二、哈密顿量:决定一切的本征方程

波函数由薛定谔方程确定,其核心是哈密顿量 H:

H ψ_nk(r) = E_nk ψ_nk(r)

H = T_e + V_ion + V_H + V_xc

  T_e   —— 电子动能
  V_ion —— 电子与原子核(离子实)的相互作用
  V_H   —— 电子间库仑排斥(Hartree 项)
  V_xc  —— 交换关联势(密度泛函理论中近似处理)

在密度泛函理论(DFT)框架下,多体问题被转化为一组单电子方程,V_xc 承担了所有多体效应的近似。哈密顿量的具体形式决定了计算方法的精度与代价:平面波赝势方法适合金属与简单半导体,全电子方法适合含 f 电子的强关联体系,而杂化泛函或 DFT+U 则用于修正带隙和局域电子。对量子材料而言,哈密顿量的选择往往比后续数据处理更影响结论——例如拓扑绝缘体的能带反转、莫特绝缘体的能隙,都高度依赖交换关联与自旋轨道耦合是否被正确纳入。

三、能带:E–k 色散关系

把本征值 E_nk 沿布里渊区高对称路径画出,就得到能带结构。它回答的是:电子能量如何随动量 k 变化。几个必须掌握的判读要点:

  • 带隙:价带顶与导带底之间的能量差。直接带隙指两者在同一 k 点,间接带隙指 k 不同。硅是间接带隙(约 1.1 eV),砷化镓是直接带隙(约 1.4 eV),这直接决定其发光效率。
  • 能带交叉与反转:能带在特定 k 点交叉,或 s、p 轨道能级次序反转,是拓扑材料(如 Bi₂Se₃)的标志性特征。
  • 有效质量:由能带曲率 m* = ħ²/(d²E/dk²) 给出。曲率越大,有效质量越小,载流子迁移率越高。
  • 费米能级位置:决定材料是金属、半导体还是绝缘体。

能带图是电子结构最直观的“指纹”,但它只沿高对称路径取样,可能遗漏布里渊区其他位置的能谷或鞍点,因此需要与态密度配合判读。

四、态密度:能量维度上的电子分布

态密度 DOS(E) 定义为单位能量区间内的电子态数目。它把整个布里渊区的信息压缩到能量轴上:

DOS(E) = Σ_n ∫_{BZ} δ(E - E_nk) dk / (2π)³

费米能级处的态密度 DOS(E_F) 决定:
  · 金属的电子比热  C_e ∝ DOS(E_F)·T
  · 泡利顺磁磁化率  χ ∝ DOS(E_F)
  · 超导转变温度(BCS 框架下)

DOS 曲线上的峰对应能带的平坦区域(范霍夫奇点),谷对应色散陡峭的区域。对量子材料,DOS(E_F) 的高低直接关联电子关联强度与磁性倾向:DOS(E_F) 高,体系易发生斯通纳失稳而出现铁磁或超导。投影态密度(PDOS)进一步把每个能量区间分解到各原子、各轨道,用于判断某个能带究竟来自谁的贡献——例如判断过渡金属氧化物中 O-2p 与金属-3d 的杂化程度。

五、费米面:金属输运与拓扑的舞台

在绝对零度下,电子填满到费米能级。在 k 空间中,满足 E_nk = E_F 的所有 k 点构成费米面。它只对金属有意义;半导体和绝缘体的费米能级落在带隙中,没有费米面。

费米面的形状决定了金属的输运性质:

  • 费米面为球形 → 近自由电子行为,如碱金属。
  • 费米面为柱状或管状 → 准二维输运,如铜氧化物超导体。
  • 费米面存在嵌套(nesting)→ 易发生自旋密度波或电荷密度波,如铬。
  • 费米面拓扑(是否包围时间反演不变点)→ 决定材料是普通金属还是拓扑半金属。

在计算中,费米面通常通过密集 k 网格上的本征值插值得到,需要足够细的网格才能准确刻画其形状与连通性。对拓扑材料,还需计算贝里曲率与陈数等拓扑不变量,判断费米面附近的能带是否具有非平庸拓扑。

从波函数到费米面:电子结构计算的逻辑链 哈密顿量 H T_e + V_ion + V_H + V_xc 决定计算精度与方法选择 波函数 ψ_nk(r) 布洛赫形式 e^{ik·r}·u_nk(r) 基组展开系数为求解未知量 本征方程 Hψ = Eψ 自洽迭代求解 能带结构 E_n(k) 沿布里渊区高对称路径的本征值色散 Γ X M Γ R E_g 导带 价带 态密度 DOS(E) 能量维度上的态数目 E DOS 范霍夫奇点 费米面 E_n(k) = E_F 仅金属存在;形状决定输运与拓扑性质 球形 近自由电子 柱状 准二维输运 嵌套 密度波失稳 量子态求解 色散关系 宏观物性

六、五个概念如何协同工作

把上述对象串起来看:哈密顿量定义问题,波函数是本征解,能带是解在 k 空间的色散,态密度是解在能量空间的统计,费米面是解在占据边界的几何。一次完整的电子结构分析通常按以下顺序推进:

  1. 确定晶体结构与磁性构型,选择交换关联泛函和赝势。
  2. 设置 k 点网格与截断能,进行自洽计算,得到收敛的电荷密度与势。
  3. 沿高对称路径做非自洽计算,输出能带;同时用更密网格计算 DOS 与 PDOS。
  4. 定位费米能级,判断金属/半导体/绝缘体,检查带隙与实验值偏差。
  5. 若为金属,提取费米面并分析其拓扑;若涉及拓扑材料,进一步计算贝里曲率与拓扑不变量。
判读提示:能带图与 DOS 必须交叉验证。若能带显示某处有交叉但 DOS 在该能量处没有明显态密度,通常是 k 路径取样不足或展宽参数设置不当。费米面计算对 k 网格密度极为敏感,金属体系建议网格不低于 24×24×24,二维材料在垂直方向至少取 1 个 k 点但面内需加密。
常见误区:不要把 DFT 计算的带隙直接等同于实验光学带隙。DFT 通常低估带隙 30%–50%,需用杂化泛函、GW 方法或剪刀修正。此外,费米能级在计算中常被设为 0 或价带顶,实际材料中它由掺杂和温度决定,解读时须明确参照约定。