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 网格上的本征值插值得到,需要足够细的网格才能准确刻画其形状与连通性。对拓扑材料,还需计算贝里曲率与陈数等拓扑不变量,判断费米面附近的能带是否具有非平庸拓扑。
六、五个概念如何协同工作
把上述对象串起来看:哈密顿量定义问题,波函数是本征解,能带是解在 k 空间的色散,态密度是解在能量空间的统计,费米面是解在占据边界的几何。一次完整的电子结构分析通常按以下顺序推进:
- 确定晶体结构与磁性构型,选择交换关联泛函和赝势。
- 设置 k 点网格与截断能,进行自洽计算,得到收敛的电荷密度与势。
- 沿高对称路径做非自洽计算,输出能带;同时用更密网格计算 DOS 与 PDOS。
- 定位费米能级,判断金属/半导体/绝缘体,检查带隙与实验值偏差。
- 若为金属,提取费米面并分析其拓扑;若涉及拓扑材料,进一步计算贝里曲率与拓扑不变量。